Comprehensive Study of Crystal Growth from Solution

Crystal growth is an essential phase in crystallization kinetics. The rate of crystal growth provides significant information for the design and control of crystallization processes; nevertheless, obtaining accurate growth rate data is still challenging due to a number of factors that prevail in crystal growth. In industrial crystallization, crystals are generally grown from multi-componentand multi-particle solutions under complicated hydrodynamic conditions; thus, it is crucial to increase the general understanding of the growth kinetics in these systems. The aim of this work is to develop a model of the crystal growth rate from solution. An extensive literature review of crystal growth focuses on themodelling of growth kinetics and thermodynamics, and new measuring techniques that have been introduced in the field of crystallization. The growth of a singlecrystal is investigated in binary and ternary systems. The binary system consists of potassium dihydrogen phosphate (KDP, crystallizing solute) and water (solvent), and the ternary system includes KDP, water and an organic admixture. The studied admixtures, urea, ethanol and 1-propanol, are employed at relatively highconcentrations (of up to 5.0 molal). The influence of the admixtures on the solution thermodynamics is studied using the Pitzer activity coefficient model. Theprediction method of the ternary solubility in the studied systems is introduced and verified. The growth rate of the KDP (101) face in the studied systems aremeasured in the growth cell as a function of supersaturation, the admixture concentration, the solution velocity over a crystal and temperature. In addition, the surface morphology of the KDP (101) face is studied using ex situ atomic force microscopy (AFM). The crystal growth rate in the ternary systems is modelled on the basis of the two-step growth model that contains the Maxwell-Stefan (MS) equations and a surface-reaction model. This model is used together with measuredcrystal growth rate data to develop a new method for the evaluation of the model parameters. The validation of the model is justified with experiments. The crystal growth rate in an imperfectly mixed suspension crystallizer is investigatedusing computational fluid dynamics (CFD). A solid-liquid suspension flow that includes multi-sized particles is described by the multi-fluid model as well as by a standard k-epsilon turbulence model and an interface momentum transfer model. The local crystal growth rate is determined from calculated flow information in a diffusion-controlled crystal growth regime. The calculated results are evaluated experimentally.

Epäpuhtauksien vaikutus monokaliumfosfaattikiteen kasvuun

Työssä on tutkittu epäpuhtauksien vaikutusta kastelulannoitesuolojen monokaliumfosfaatin, kaliumnitraatin ja ureafosfaatin kiteytyksessä. Kirjallisuusosassa on käsitelty kastelulannoiteprosessit ja epäpuhtauksien vaikutus kastelulannoitteiden valmistuksessa. Kiteytys ja kiteenkasvu on esitetty perusyhtälöin, joissa on otettu epäpuhtauksien vaikutus huomioon. Tarkemmin on perehdytty monokaliumfosfaatin kiteytykseen ja kolmenarvoisten kationeiden, Al3+, Fe3+ ja Cr3+, vaikutukseen kiteiden kasvuun. Kolmenarvoiset metalli-ionit adsorboituvat kiteen pintaan haitaten kiteenkasvua, mikä vaikuttaa erityisesti kiteen prismapinnan kasvuun. Lisäksi on esitelty muita kiteenkasvuun vaikuttavia olosuhteita. Lopuksi on käsitelty kompleksinmuodostajia metalli-ionien haitallisten vaikutusten ehkäisijöinä. Kokeellisessa osassa suoritettiin liukoisuuskokeita monokaliumfosfaatin liukoisuuden selvittämiseksi eri pH-olosuhteissa. Suoritetuissa yksikidekokeissa tutkittiin pH:n ja kolmenarvoisten kationeiden; Al3+, Fe3+ ja Cr3+, vaikutus monokaliumfosfaattikiteen pituus- ja leveyskasvuun ja kidemuotoon eri ylikylläisyyksillä. Lisäksi tutkittiin voidaanko lämpötilaa ja pH muuttamalla tai pyrofosfaattia lisäämällä poistaa raudan kasvua inhiboima vaikutus. Kiteytyslämpötilaa nostamalla voidaan poistaa raudan haitallinen vaikutus kiteen kasvuun.

Morfologian hallinta nikkelin sulfidisaostuksessa

Diplomityön tarkoituksena oli tutkia nikkelin sulfidisaostuksessa syntyvien kiteiden morfologiaa ja siihen vaikuttavia parametreja. Syntyvien kiteiden kasvua ja morfologiaa tutkittiin kiteen muodostumisen ja kasvun teorioiden avulla. Saostuksen olosuhteet, kuten lämpötila, paine ja pH vaikuttavat muodostuvien kiteiden morfologiaan. Muilla parametreilla, kuten liuoksen ylikylläisyydellä, epäpuhtauksilla, lisäaineilla, sekoituksella ja reaktioajalla on myös suuri merkitys. Kokeiden avulla haluttiin liuoskoostumuksen, saostusolosuhteiden ja muiden komponenttien vaikutusta nikkelisulfidikiteiden morfologiaan. Kokeissa käytettiin kahta eri sulfidilähdettä: natriumvetysulfidia ja rikkivetyä. Puolipanoskokeissa nikkelipitoisuus oli 1,5 g/l, paine 101,3 kPa ja sekoitusnopeus 650 rpm. Saostuskokeet tehtiin natriumsulfaatti- 5 g/l ja ammoniumsulfaattiliuoksissa 300 g/l. Saostuskokeissa muuttujia olivat saostimen konsentraatio ja määrä, rauta- ja magne-siumepäpuhtaudet, lämpötila ja lisäaineet. Diplomityön kokeellisessa osassa morfologiaa tutkittiin suoraan valomikroskoopin ja pyyhkäisyelektronimikroskoopin (SEM) avulla. Morfologiaa tutkittiin myös epäsuorasti laskeutumisnopeuden, keskimääräisen partikkelikoon, ja ominaispinta-alamittausten avulla. Saostimen pitoisuuden vaikutukset partikkelimuotoon olivat pieniä, mutta vaikutukset ominaispinta-alaan ja partikkelikokoon olivat suuria. Natriumlauryylisul-faatti ja EDTA ohjasivat partikkelien rakennetta levymäisemmäksi, joka johti hitaaseen laskeutumisnopeuteen. Polyakryylihappo lisäaineena muuttaa partikkelien morfologiaa kuutiomaisemmaksi. Flokkulanttien ja raudan morfologiset vaikutukset olivat pieniä. Partikkelikoko ja omaispinta-ala pienenivät selvästi magnesiumpitoisuuden kasvaessa. Lämpötilan kasvattaminen lisäsi epäsäännöllisten kiteiden määrää ja muodostuneet kiteet olivat enemmän neulamaisia.

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

Towards desired crystalline product properties: In-situ monitoring of batch crystallization

The objective of industrial crystallization is to obtain a crystalline product which has the desired crystal size distribution, mean crystal size, crystal shape, purity, polymorphic and pseudopolymorphic form. Effective control of the product quality requires an understanding of the thermodynamics of the crystallizing system and the effects of operation parameters on the crystalline product properties. Therefore, obtaining reliable in-line information about crystal properties and supersaturation, which is the driving force of crystallization, would be very advantageous. Advanced techniques, such asRaman spectroscopy, attenuated total reflection Fourier transform infrared (ATR FTIR) spectroscopy, and in-line imaging techniques, offer great potential for obtaining reliable information during crystallization, and thus giving a better understanding of the fundamental mechanisms (nucleation and crystal growth) involved. In the present work, the relative stability of anhydrate and dihydrate carbamazepine in mixed solvents containing water and ethanol were investigated. The kinetics of the solvent mediated phase transformation of the anhydrate to hydrate in the mixed solvents was studied using an in-line Raman immersion probe. The effects of the operation parameters in terms of solvent composition, temperature and the use of certain additives on the phase transformation kineticswere explored. Comparison of the off-line measured solute concentration and the solid-phase composition measured by in-line Raman spectroscopy allowedthe identification of the fundamental processes during the phase transformation. The effects of thermodynamic and kinetic factors on the anhydrate/hydrate phase of carbamazepine crystals during cooling crystallization were also investigated. The effect of certain additives on the batch cooling crystallization of potassium dihydrogen phosphate (KDP) wasinvestigated. The crystal growth rate of a certain crystal face was determined from images taken with an in-line video microscope. An in-line image processing method was developed to characterize the size and shape of thecrystals. An ATR FTIR and a laser reflection particle size analyzer were used to study the effects of cooling modes and seeding parameters onthe final crystal size distribution of an organic compound C15. Based on the obtained results, an operation condition was proposed which gives improved product property in terms of increased mean crystal size and narrowersize distribution.

Sulfatiatsolin leijupetikiteytys

Työssä tutkittiin leijupetikiteyttimen toimivuutta sulfatiatsolikiteiden kasvunopeuksien mittaamisessa. Sulfatiatsoli on lääkeaine, jota käytetään antibioottina. Kirjallisuusosassa on käsitelty kiteytyksen perusteita sekä olosuhteiden vaikutuksia kidemorfologiaan. Koska kyseessä on lääkeaine, on työssä myös selvitetty tekijöitä, jotka vaikuttavat syntyvän kiteen polymorfimuotoon. Näitä ovat mm. siemenkiteen polymorfimuoto sekä käytetty liuotin. Kirjallisuusosassa on myös esitelty teollisia ja laboratoriomittakaavan leijupetikiteyttimiä. Kokeellisessa osassa on keskitytty testaamaan leijupetikiteyttimen toimintaa sekä etsimään sopivia olosuhteita kasvukokeiden suorittamiselle. Kokeissa tutkittiin ylikylläisyyden vaikutusta sulfatiatsolikiteen kasvuun. Kiteen kasvun seurantaan pyrittiin etsimään sopivia hiukkaskokoanalysaattoreita. Mittauksissa käytetyt laitteistot olivat PIA 4000 mod LUT on-line-videomikroskooppi ja laserdiffraktioon perustuva Coulter LS 130 off-line-partikkelikokoanalysaattori.

NMR in Crystals and Powders of Topazes with Different Colours

The present work is a part of the large project with purpose to investigate microstructure and electronic structure of natural topazes using NMR method. To reach this task we determined the relative contents of fluorine and hydrogen in crystals blue, colorless, wine and wine irradiated topazes. Then we determined the electric field gradients in site of aluminium atoms by NMR method, calculated EFG using ab initio method, and measured relaxation time dependence on heating temperature for blue, colorless, Swiss blue and sky blue topazes. Nuclear magnetic resonance (NMR) is an effective method to investigate the local structure in the crystal. The NMR study of the single crystal gives detailed information especially about the local crystal structure. As a result of this work we have received practical data, which is possible to use in future for making personal dosimetry and for preparation of mullite, which is widely used in traditional and advanced ceramic materials.

Galvanomagnetic properties of CdSb

Main aim of this work was preparation of a computer program for investigation of galvanomagnetic effects in solid state materials. These effects were investigated in magnetic field up to 6 T at temperatures 4.6 and 80.5 K. Two CdSb samples with Ni shallow impurities (concentration of impurity was 5% by mass) and one undoped CdSb single crystal were studied. Obtained results were compared with previous experimental results for these samples, and showed their identity.

Effect of hydrodynamics on modelling, monitoring and control of crystallization

Crystallization is a purification method used to obtain crystalline product of a certain crystal size. It is one of the oldest industrial unit processes and commonly used in modern industry due to its good purification capability from rather impure solutions with reasonably low energy consumption. However, the process is extremely challenging to model and control because it involves inhomogeneous mixing and many simultaneous phenomena such as nucleation, crystal growth and agglomeration. All these phenomena are dependent on supersaturation, i.e. the difference between actual liquid phase concentration and solubility. Homogeneous mass and heat transfer in the crystallizer would greatly simplify modelling and control of crystallization processes, such conditions are, however, not the reality, especially in industrial scale processes. Consequently, the hydrodynamics of crystallizers, i.e. the combination of mixing, feed and product removal flows, and recycling of the suspension, needs to be thoroughly investigated. Understanding of hydrodynamics is important in crystallization, especially inlargerscale equipment where uniform flow conditions are difficult to attain. It is also important to understand different size scales of mixing; micro-, meso- and macromixing. Fast processes, like nucleation and chemical reactions, are typically highly dependent on micro- and mesomixing but macromixing, which equalizes the concentrations of all the species within the entire crystallizer, cannot be disregarded. This study investigates the influence of hydrodynamics on crystallization processes. Modelling of crystallizers with the mixed suspension mixed product removal (MSMPR) theory (ideal mixing), computational fluid dynamics (CFD), and a compartmental multiblock model is compared. The importance of proper verification of CFD and multiblock models is demonstrated. In addition, the influence of different hydrodynamic conditions on reactive crystallization process control is studied. Finally, the effect of extreme local supersaturation is studied using power ultrasound to initiate nucleation. The present work shows that mixing and chemical feeding conditions clearly affect induction time and cluster formation, nucleation, growth kinetics, and agglomeration. Consequently, the properties of crystalline end products, e.g. crystal size and crystal habit, can be influenced by management of mixing and feeding conditions. Impurities may have varying impacts on crystallization processes. As an example, manganese ions were shown to replace magnesium ions in the crystal lattice of magnesium sulphate heptahydrate, increasing the crystal growth rate significantly, whereas sodium ions showed no interaction at all. Modelling of continuous crystallization based on MSMPR theory showed that the model is feasible in a small laboratoryscale crystallizer, whereas in larger pilot- and industrial-scale crystallizers hydrodynamic effects should be taken into account. For that reason, CFD and multiblock modelling are shown to be effective tools for modelling crystallization with inhomogeneous mixing. The present work shows also that selection of the measurement point, or points in the case of multiprobe systems, is crucial when process analytical technology (PAT) is used to control larger scale crystallization. The thesis concludes by describing how control of local supersaturation by highly localized ultrasound was successfully applied to induce nucleation and to control polymorphism in reactive crystallization of L-glutamic acid.

Understanding the in vitro dissolution rate of glasses with respect to future clinical applications

Glass is a unique material with a long history. Several glass products are used daily in our everyday life, often unnoticed. Glass can be found not only in obvious applications such as tableware, windows, and light bulbs, but also in tennis rackets, windmill turbine blades, optical devices, and medical implants. The glasses used at present as implants are inorganic silica-based melt-derived compositions mainly for hard-tissue repair as bone graft substitute in dentistry and orthopedics. The degree of glass reactivity desired varies according to implantation situation and it is vital that the ion release from any glasses used in medical applications is controlled. Understanding the in vitro dissolution rate of glasses provides a first approximation of their behavior in vivo. Specific studies concerning dissolution properties of bioactive glasses have been relatively scarce and mostly concentrated to static condition studies. The motivation behind this work was to develop a simple and accurate method for quantifying the in vitro dissolution rate of highly different types of glass compositions with interest for future clinical applications. By combining information from various experimental conditions, a better knowledge of glass dissolution and the suitability of different glasses for different medical applications can be obtained. Thus, two traditional and one novel approach were utilized in this thesis to study glass dissolution. The chemical durability of silicate glasses was tested in water and TRIS-buffered solution at static and dynamic conditions. The traditional in vitro testing with a TRISbuffered solution under static conditions works well with bioactive or with readily dissolving glasses, and it is easy to follow the ion dissolution reactions. However, in the buffered solution no marked differences between the more durable glasses were observed. The hydrolytic resistance of the glasses was studied using the standard procedure ISO 719. The relative scale given by the standard failed to provide any relevant information when bioactive glasses were studied. However, the clear differences in the hydrolytic resistance values imply that the method could be used as a rapid test to get an overall idea of the biodegradability of glasses. The standard method combined with the ion concentration and pH measurements gives a better estimate of the hydrolytic resistance because of the high silicon amount released from a glass. A sensitive on-line analysis method utilizing inductively coupled plasma optical emission spectrometer and a flow-through micro-volume pH electrode was developed to study the initial dissolution of biocompatible glasses. This approach was found suitable for compositions within a large range of chemical durability. With this approach, the initial dissolution of all ions could be measured simultaneously and quantitatively, which gave a good overall idea of the initial dissolution rates for the individual ions and the dissolution mechanism. These types of results with glass dissolution were presented for the first time during the course of writing this thesis. Based on the initial dissolution patterns obtained with the novel approach using TRIS, the experimental glasses could be divided into four distinct categories. The initial dissolution patterns of glasses correlated well with the anticipated bioactivity. Moreover, the normalized surface-specific mass loss rates and the different in vivo models and the actual in vivo data correlated well. The results suggest that this type of approach can be used for prescreening the suitability of novel glass compositions for future clinical applications. Furthermore, the results shed light on the possible bioactivity of glasses. An additional goal in this thesis was to gain insight into the phase changes occurring during various heat treatments of glasses with three selected compositions. Engineering-type T-T-T curves for glasses 1-98 and 13-93 were stablished. The information gained is essential in manufacturing amorphous porous implants or for drawing of continuous fibers of the glasses. Although both glasses can be hot worked to amorphous products at carefully controlled conditions, 1-98 showed one magnitude greater nucleation and crystal growth rate than 13-93. Thus, 13-93 is better suited than 1-98 for working processes which require long residence times at high temperatures. It was also shown that amorphous and partially crystalline porous implants can be sintered from bioactive glass S53P4. Surface crystallization of S53P4, forming Na2O∙CaO∙2SiO2, was observed to start at 650°C. The secondary crystals of Na2Ca4(PO4)2SiO4, reported for the first time in this thesis, were detected at higher temperatures, from 850°C to 1000°C. The crystal phases formed affected the dissolution behavior of the implants in simulated body fluid. This study opens up new possibilities for using S53P4 to manufacture various structures, while tailoring their bioactivity by controlling the proportions of the different phases. The results obtained in this thesis give valuable additional information and tools to the state of the art for designing glasses with respect to future clinical applications. With the knowledge gained we can identify different dissolution patters and use this information to improve the tuning of glass compositions. In addition, the novel online analysis approach provides an excellent opportunity to further enhance our knowledge of glass behavior in simulated body conditions.

Microcracks and Vortices in Superconducting Thin Films

In this work, superconducting YBa₂ Cu₃O_6+x (YBCO) thin films have been studied with the experimental focus on the anisotropy of BaZrO₃ (BZO) doped YBCOthin films and the theoretical focus on modelling flux pinning by numerically solving Ginzburg- Landau equations. Also, the structural properties of undoped YBCO thin films grown on NdGaO₃ (NGO) and MgO substrates were investigated. The thin film samples were made by pulsed laser ablation on single crystal substrates. The structural properties of the thin films were characterized by X-ray diffraction and atomic force microscope measurements. The superconducting properties were investigated with a magnetometer and also with transport measurements in pulsed magnetic field up to 30 T. Flux pinning was modelled by restricting the value of the order parameter inside the columnar pinning sites and then solving the Ginzburg-Landau equations numerically with the restrictions in place. The computations were done with a parallel code on a supercomputer. The YBCO thin films were seen to develop microcracks when grown on NGO or MgO substrates. The microcrack formation was connected to the structure of the YBCO thin films in both cases. Additionally, the microcracks can be avoided by careful optimization of the deposition parameters and the film thickness. The BZO doping of the YBCO thin films was seen to decrease the effective electron mass anisotropy, which was seen by fitting the Blatter scaling to the angle dependence of the upper critical field. The Ginzburg-Landau simulations were able to reproduce the measured magnetic field dependence of the critical current density for BZO doped and undoped YBCO. The simulations showed that in addition to the large density also the large size of the BZO nanorods is a key factor behind the change in the power law behaviour between BZO doped and undoped YBCO. Additionally, the Ginzburg-Landau equations were solved for type I thin films where giant vortices were seen to appear depending on the film thickness. The simulations predicted that singly quantized vortices are stable in type I films up to quite large thicknesses and that the size of the vortices increases with decreasing film thickness, in a way that is similar to the behaviour of the interaction length of Pearl vortices.

The Role of Structural Imperfections in Sr2FeMoO6 Thin Films

In this work, Sr2FeMoO6 (SFMO) thin films were studied with the main focus on their magnetic and magneto-transport properties. The fabrication process of pulsed laser deposited SFMO films was first optimized. Then the effects of strain, film thickness and substrate were thoroughly investigated. In addition to these external factors, the effect of intrinsic defects on the magnetic properties of SFMO were also clarified. Secondly, the magnetoresistivity mechanims of SFMO films were studied and a semiempirical model of the temperature dependence of resistivity was introduced. The films were grown on single crystal substrates using a ceramic target made with sol-gel method. The structural characterization of the films were carried out with X-ray diffraction, atomic force microscopy, transmission electron microscopy and high kinetic energy photoelectron spectroscopy. The magnetic properties were measured with SQUID magnetometer and the magneto-transport properties by magnetometer with a resistivity option. SFMO films with the best combination of structural and magnetic properties were grown in Ar atmosphere at 1050 °C . Their magnetic properties could not be improved by the ex situ post-annealing treatments aside from the treatments in ultra-high vacuum conditions. The optimal film thickness was found to be around 150 nm and only small improvement in the magnetic properties with decreasing strain was observed. Instead, the magnetic properties were observed to be highly dependent on the choice of the substrate due to the lattice mismatch induced defects, which are best avoided by using the SrTiO3 substrate. The large difference in the Curie temperature and the saturation magnetization between the SFMO thin film and polycrystalline bulk samples was connected to the antisite disorder and oxygen vacancies. Thus, the Curie temperature of SFMO thin films could be improved by increasing the amount of oxygen vacancies for example with ultra-high vacuum treatments or improving the B-site ordering by further optimization of the deposition parameters. The magneto-transport properties of SFMO thin films do not follow any conventional models, but the temperature dependence of resistivity was succesfully described with a model of two spin channel system. Also, evidences that the resistivity-temperature behaviour of SFMO thin films is dominated by the structural defects, which reduce the band gap in the majority spin band were found. Moreover, the magnetic field response of the resistivity in SFMO thin films were found to be superposition of different mechanisms that seems to be related to the structural changes in the film.