Subscriber access provided by READING UNIV
Review
Determination methods for crystal nucleation kinetics in solutions Yan Xiao, Jingkang Wang, Xin Huang, Huanhuan Shi, Yanan Zhou, Shuyi Zong, Hongxun Hao, Ying Bao, and Qiuxiang Yin Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.7b01223 • Publication Date (Web): 30 Nov 2017 Downloaded from http://pubs.acs.org on November 30, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Crystal Growth & Design is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
Determination methods for crystal nucleation kinetics in solutions Yan Xiao a,b, Jingkang Wang a,b, Xin Huang a,b, Huanhuan Shi a,b, Yanan Zhou a,b, Shuyi Zong a,b, Hongxun Hao a,b *, Ying Bao a,b, Qiuxiang Yin a,b
a
State Key Laboratory of Chemical Engineering, School of Chemical Engineering and
Technology, Tianjin University, Tianjin 300072, P. R. China b
Collaborative Innovation Center of Chemical Science and Chemical Engineering
(Tianjin), Tianjin 300072, P.R. China
Corresponding Author *Tel: +86-022-27405754. Fax: +86-022-27374971. E-mail:
[email protected] 1
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Abstract: Crystal nucleation process in solutions plays an important role in both the engineering processes and the physical science. To seek the truth of nucleation, researchers have worked almost a century from either the kinetic viewpoint or the solution chemistry viewpoint. In this work, the nucleation rate measurement methods, including deterministic method, droplet based method, double-pulse method, microfluidic method and stirred small volume solution method, were reviewed and their advantages and disadvantages were also discussed. Furthermore, problems in the nucleation kinetic investigation were also listed in the end.
Key words: Crystal nucleation, Nucleation rate measurement, Deterministic method, Small volume solution methods
Introduction Nucleation phenomenon exists everywhere. It plays a significant role not only in the engineering processes, including condensation, boiling, crystallization, sublimation, etc., but also in the biology science, involved in the development of teeth, bone, kidney stones, gallstones, as well as the protein crystallization related diseases. Besides, it is also responsible for the formation of crystalline materials on earth and certain aerosols which cause air pollution. Despite the importance of nucleation process, it has not been fully understood.
1–6
According to the origin of the nuclei,
nucleation can be divided into primary nucleation, which means that nuclei are formed from absolute pure supersaturated multi-component system (or supercooled melt system for single component), known as homogeneous nucleation, or are catalyzed by the presence of existing foreign surfaces (such as dust, container walls, impellers), known as heterogenous nucleation, and secondary nucleation, where nucleation is induced by the prior presence of seeds from the same material. 7,8 Nucleation is the first step of crystallization process, which may influence the crystal size and its distribution, polymorph, purity of the crystalline products in many circumstances. Since many factors will affect the nucleation process, it is still quite 2
ACS Paragon Plus Environment
Page 2 of 36
Page 3 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
difficult to design an ideal nucleation process to obtain products with desired physicochemical properties. Recently, although the development of spectroscopic techniques9–11 as well as the molecular simulation methodology12,13 enables us to investigate the molecules assembly state in solution, which may be related to the mature macroscopic crystal structures, there is still a long way to directly observe and characterize the nucleation process.14,15 Besides, many researchers have also worked a lot on the process modelling16 and the non-steady state investigation17,18 through solving mathematical models based on the theory of intermediate stage of crystal growth. Modelling results showed good agreements with the supersaturation decay curve and particle size distributions for several proteins19,20 whose data were easily found from literatures. However, more experimental data is needed to help verify and further modify the theoretical results. In classical nucleation theory, nucleation process was simplified as a reversible chain reaction, molecules successively attach to and detach from the clusters of different sizes, and the frequency of clusters past the critical size is denoted as the nucleation rate.21 Therefore, molecular understanding of nucleation process could be inferred from the accurate nucleation rate measurements. For decades, investigations on nucleation kinetics remains one of the most important ways to gain information of nucleation.22,23 Roger et al and Peter et al
2,24–26
have
reviewed the classical nucleation theory and two-step nucleation theory to interpret nucleation rate data for different supersaturations and solvents, or even different compounds, to distract molecular scale information of nucleation process. The prevalent nucleation rate from solutions is defined as the number of crystals nucleated per unit volume (or area) during per unit time.27 Starting from the definition, nucleation rate can be estimated directly by measuring the number or number density of crystals versus time profile. Alternatively, the time needed for the appearance of “first crystal” can also be used to calculate the nucleation rate. Two broad classes of experimental methods are mainly used to characterize the elapsed time period, polythermal method, which measures the metastable zone width (MSZW), temperature (time) difference, at different constant cooling rates for the saturated 3
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
solution to nucleate, and isothermal method, which determines the induction time needed for the commencing of nucleation at constant temperature for supersaturated solutions with different supersaturation ratios. And the “first crystal” is usually treated as the point at which the number density or the volume fraction of crystals reaches some critical values. Therefore, both the MSZW and induction time are related to the nucleation and growth rate of crystallization process, which could provide the kinetics information indirectly.28–31 Many methods have been proposed to acquire massive nucleation rate data of solutions. As for the detection of nucleation process, because of the tiny size of critical nucleus and the limitation of analytical technique threshold, the detection could not be done directly and instantly. Therefore, indirect methods, such as monitoring the system property changes, seem to be better choices. In the publications, microscopes32 can be used to observe or count the number of crystals which already grow into big size, while laser33, conductivity meter34, thermometer35, and even bare eyes can also be used to monitor the system property changes (refractivity, conductivity, temperature etc.) which are caused by nucleation and the appearance of massive crystals. As to the devices used for measuring nucleation kinetics, jacketed vessels36, quasi 2-dimensional cells37, well plates38, electrodynamic levitation trap apparatus39 are widely employed, and the medium or high throughput devices, such as Crystal 1640, Crystalline parallel reactors, microfluidic devices41, are getting more popular for investigating nucleation process because of easy operating and capability of indicating the stochastic nature of nucleation. The measurement of nucleation kinetics in solutions can also be divided into deterministic (also called traditional) method and probability method, based on the multi nuclei model and single nucleus model, respectively.42 Given the great importance of crystal nucleation kinetics in solutions, this work will summarize the methods used to obtain the crystal nucleation rate. Their corresponding pros and cons will be discussed. The future development of the crystal nucleation kinetics estimation will also be discussed in the end. 4
ACS Paragon Plus Environment
Page 4 of 36
Page 5 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
Deterministic method For large or industrial scale crystallizers, crystal nucleation process tends to be a deterministic (repeatable) process, as a combination of primary nucleation, secondary nucleation, crystal growth, crystal break and agglomeration, and so on. Also, agitation is an essential process, which enhances nucleation in most of the crystallization engineering process. Deterministic method is the most widely used one to estimate nucleation kinetics for the scale up or controlling of crystallization process for many years as a result of its similar complexity with the real industrial process. According to the properties measured to characterize the nucleation process, it can be divided into three subcategories: particle number based method, induction time or MSZW based method and population balance equation (PBE), mass balance equation combined with empirical kinetic equations based method. In 1978, Dugua et al34 measured the nucleation rates of sodium perborate from aqueous solutions with and without the addition of sodium salt of butyl ester oleic acid as the surfactant. Nucleation rates were estimated from the total precipitated crystal mass profile by keeping the supersaturation constant artificially (solution concentration stayed the same and was monitored by using conductivity method), which could be related to the number of crystals, crystal morphology and crystal growth rate. In pure solvents, the nucleation rate was a staircase-form function of time, while the rate was reduced and became more irregular with the presence of surfactants. Bhamidi et al43,44 measured the crystal nucleation rates of hen egg-white lysozyme at different ionic strength using a particle counter, which used the principle of near angle light scatter to detect particles and could directly export the number density of crystals in the crystallizing solution. The nucleation rates, determined from the slope of the crystal number density versus time plots, were correlated by the classical nucleation theory expression. The systematic deviation of model predictions toward the lower protein concentration region was illustrated by the authors as that two nucleation mechanisms may exist: heterogeneous nucleation at lower protein concentration and 5
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 36
homogeneous nucleation at higher protein concentration. Therefore, the data was split into two regions and the corresponding kinetic parameters were estimated separately. Benefiting from the convenient colloid particle length scale (set by particle size, ~μm) and time scale (set by particle diffusion, ~ms to days), studies of colloids crystallization could be investigated using light scattering techniques45. Cheng et al33,46 used time solved Bragg light scattering to characterize crystal nucleation and crystal growth of colloidal hard spheres in a microgravity environment on a space shuttle. Also, they used UV-visible transmission spectroscopy to investigate the nucleation behavior of poly-N-isopropylacrylamide (PNIPAM) particles, which could be described by empirical power law relations. In-situ focused beam reflectance measurement
(FBRM)
was
used
to
evaluate
the
nucleation
rate
of
hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) from γ-butyrolactone in unseeded batch cooling crystallization by Kim et al36. The nucleation rates were estimated from the constant slope of total counts/s versus time plots recorded by FBRM, shown in Fig. 1, before which a conversion from the counts/s to the actual number density of crystals was performed. By introducing a temperature term concerning the dependence of kinetics on saturation temperature, a modified classical nucleation rate model was proposed to correlate the nucleation rate data and showed a better prediction for the MSZW in the cooling crystallization.
Fig.1 Measurements of tonset, Tonset, and JFBRM from the operating profile in the unseeded batch cooling crystallization of RDX (Tsat,i = 343.15 K and q = 0.1 K/min). (Reproduced with permission from ref 36. Copyright 2013 John Wiley and Sons.) 6
ACS Paragon Plus Environment
Page 7 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
Liu et al47 measured the nucleation kinetics of L-arginine trifluoroacetate (LATF) in water through induction time measurements, whose reciprocal was correlated to the nucleation rates. Results showed that nucleation is more probable for higher supersaturation. Similarly, nucleation rates of nesquehonite from water and water/dimethylformamide (DMF) solvent mixtures were determined by the reciprocal of induction time values to examine solvent effect on the nucleation kinetics by Zhao et al48 (Fig. 2), while the induction time was calculated from the in-situ conductivity detection of the solution. Results showed a positive correlation between nucleation rate and DMF concentration. Analysis based on the classical nucleation theory revealed an unexpected increase of surface energy with the increase of DMF content, implying that kinetic parameter, rather than the thermodynamic parameter resulted in the enhancement of nucleation. Further analysis suggested that kinetic acceleration was possible with high surface energy if the nucleation process underwent a cluster aggregation mechanism instead of critical nucleus formation. Shiau et al49 investigated the solvent influence on the pre-exponential factor and interfacial energy by measuring the MSZW data of salicylamide in different solvents. An integral model based on the classical nucleation theory and Stokes-Einstein equation was developed to determine the intrinsic nucleation rate constant and interfacial energy. Results indicated that the interfacial energy went up with the increase of solvent boiling point and the nucleation constant showed a good linear relationship with the interfacial energy.
7
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 36
Fig.2 Measured induction time of nesquehonite nucleation in pure H2O and H2O–DMF binary solvents at different saturation levels (A) before and (B) after solubility correction. The trend lines are drawn only for the purpose of guiding the eye. (Reproduced with permission from ref 48. Copyright 2013 Elsevier.)
PBE model and the mass balance equation, together with all kinds of in-situ process analytical technologies (PAT), are widely used in the estimation of nucleation and growth kinetics for industrial crystallization process. Trifkovic et al50 estimated the nucleation and growth parameters of paracetamol in isopropanol-water antisolvent batch crystallization process by using the nonlinear regression method for the moments of the crystal population density. The moments were calculated from the measured chord length distribution (CLD) generated by FBRM, while the concentration
profile was calculated through
the
ATR-FTIR
spectroscopy
measurements, as shown in the schematic experimental setup, Fig. 3. Both the 8
ACS Paragon Plus Environment
Page 9 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
nucleation rate and growth rate expressions adopted the power law correlations51, derived from the classical nucleation theory (CNT). Open-loop experimental results validated the crystallization model, using the estimated kinetic parameters. And the single- and multi-objective optimization on the antisolvent flow rate profiles by using the obtained model improved both the particle size and its distribution, as well as the yield. In the same year, Nagy et al52 used the same method established the kinetic model for the cooling crystallization process of paracetamol from water by MSZW measurements. The combination of data obtained from RAMAN which was used to confirm the polymorph composition in the precipitated crystals, ATR-FTIR which was used to monitor the concentration, and FBRM which was used to gather the crystal size distribution information, was used by Trifkovic et al53 to determine the solvent mediated polymorph transformation kinetics for Buspirone Hydrochloride (BUS-HCl) by separating the transformation experiment into three processes: the dissolution of metastable Form II, the nucleation of stable Form I and the growth of stable Form I.
Fig.3 Experimental setup. (Reproduced with permission from ref 50. Copyright 2008 American Chemical Society.)
9
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
As we can see, experiment results show no big disperse and tend to be constant in the deterministic method. Therefore, several measurements will be enough to estimate the kinetic parameters and establish a process model. In terms of the setup, it is relatively simple and the macroscopic size and stirring provides similar fluid dynamic conditions with the scaled-up crystallizers. Also, this method is crystallization method independent and it can be used in cooling crystallization, antisolvent crystallization, and evaporation crystallization, either seeded or un-seeded crystallization, batch or continuous crystallization, and so on. However, despite the advantages of time-saving and easy to scale up, the deterministic method is hard to operate at the desired temperature or supersaturation under exact uniform fluid dynamic conditions as a result of the large dimensions of the setup and transport phenomenon, rendering the system difficult to operate, control and analyze. Also, many assumptions are applied in this kind of method,50 such as, the crystallizer is well mixed, crystal agglomeration and breakage phenomenon are neglected or simplified in the PBE model, no growth dispersion happens, the crystals are born at size zero, etc. Usually, for simplicity, the crystal morphology or size is described by only one variable. Furthermore, the poor conversion from the chord length distribution (CLD) and total counts/s measured by FBRM to the crystal size distribution (CSD) and crystal number density restricted the wider use of FBRM and the PBE model method.54–56
Droplet based method Crystal nucleation process usually happens quite fast so that it is not very easy and suitable to testify the prediction of nucleation theory, as a result of the unknown and inevitable catalytic surfaces or centers existing in the bulk system. Also, the inevitable existence of secondary nucleation affects the investigation of primary nucleation. Theoretically, if the system is subdivided into sufficiently small droplets, a large portion of the droplets should be free or at least having substantially reduced amount of heterogeneous nucleant, which provide a method for investigating the homogenous 10
ACS Paragon Plus Environment
Page 10 of 36
Page 11 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
nucleation process in supersaturated systems. The small droplet based technique has been widely used to gather nucleation rate data for different systems. Vonnegut57 was the first to recognize that miniaturization of a macroscopic sample into non-communicating micro samples will suppress the effects of heterogeneous nucleation. The isothermal nucleation rates of supercooled liquid tin were determined using dilatometric measurements. The crystallization systems (Fig. 4) was composed of many small droplets (1-10 micron diameter) dispersed by manual shaking and isolated by an oxide film. Results from different temperatures showed that the nucleation rates had a very large negative temperature coefficient corresponding to an activation energy of the order of ~ 2 × 105 calories. In 1952, Turnbull58 investigated the solidification kinetics of supercooled liquid mercury using the same method with droplets diameter of 2-8 microns. (a)
11
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
(b)
Fig.4 (a) Dilatometer. (b) Fraction of tin drops remaining unfrozen as function of time from dilatometer data. (Reproduced with permission from ref 57. Copyright 1948 Elsevier.)
In addition to the volume change during the solidification of supercooled liquid, the heat released from crystallization process can also be measured to characterize nucleation rate. Rasmussen et al35 applied differential scanning calorimetry (DSC) to monitor the heat effect during the solidification of an alloy (90% tin-10% bismuth) emulsion droplet, stabilized in organic carried fluid. By rapid cooling to isothermal holding temperature or using a constant cooling rate, the temperature dependence of nucleation rate could be derived. In 1959, White et al59 extended this technique to nucleation process from solutions. By using a droplet generating device, small droplets of potassium nitrate aqueous solution were produced continuously and isolated by mineral oil. The nucleation kinetics was recorded by observing the fraction of uncrystallized droplets with time by using microscope. Results showed that with ordinary solutions, nucleation was catalyzed by impurities, which was heterogeneous nucleation mechanism, while solutions without catalyzing centers which were removed by using a partial crystallization technique, turned to follow homogeneous nucleation process. In 1978, Keller et al60 applied this method to acquire the reaction crystallization 12
ACS Paragon Plus Environment
Page 12 of 36
Page 13 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
nucleation rate of CaSO4•2H2O aqueous solution by using an in-situ generation of insoluble substance from very soluble precursors protocol. Droplet based technique does diminish the effect of unknown catalytical centers in the system and testify the rationality of classical nucleation theory. It can also save materials used to measure nucleation kinetics, which is suitable for expensive products. The easy operation of dispersing droplets places less limits for the studied system, which is specifically suitable for high melting point materials.35,61 However, almost all the works ignored the crystal growth time which is defined as the elapsed time from the creation of nuclei to growing to detectable crystals. Besides, the generated droplets usually had a relatively large size variation which would influence the accuracy of the experimental results, as illustrated by Vonnegut and Filippone et al.57,61
Double-pulse method One of the main reasons that nucleation is difficult to study is the nano-sized critical nuclei, which makes the observation and locating of them almost impossible with the techniques to date.14,15 According to the idea proposed by Tammann38 in 1922, this question could be solved by separating the nucleation and growth processes of crystallization, which is called double-pulse method and is schematically shown in Fig. 5(a)62. Its essence is that, during the first stage, crystals are only nucleated by creating high supersaturations for a relatively short time (minutes, sometimes hours) and crystal growth is not fast enough to exhaust the materials and decrease the supersaturation. Then, the supersaturation is rapidly lowered to the metastable zone, where the system is unable to produce further nuclei and only growth of the existed nuclei could happen. For slow growing molecules, the growth time could be set as several hours or even days. By counting the number of macroscopic crystals formed during a certain nucleation time, nucleation rate could be estimated, as shown in Fig. 5(b)63. This method eliminates the error resulted from ignoring growth times and has been successfully applied to many systems, especially proteins. 13
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 36
(a)
(b)
Fig.5 (a) Schematic representation of the classical principle of separation in time of the nucleation and growth stages. (Reproduced with permission from ref 62. Copyright 2010 American Chemical Society.) (b) n [crystallite number] vs. t [minutes] dependencies for surfaces rendered hydrophobic (by means of hexamethyl-disilazane). The corresponding concentration ratios σ = C/Csat are shown in the graph’s; C is the actual concentration, Csat the equilibrium one. (Reproduced with permission from ref 63. Copyright 1999 Elsevier.)
In 1954, Dunning and Shipman64 used a counting cell which was a flat glass vessel, 13 cm diameter and 1 cm deep, to obtain the nucleation rate of sucrose solution at different supersaturations and temperatures. The nucleation rate was obtained from the slope of crystal number vs. time, while the intercept on the time axis was the induction time. Based on Volmer’s theory, the rate controlling step was 14
ACS Paragon Plus Environment
Page 15 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
considered to be the frequency of molecular clusters growing beyond the critical size and the corresponding activation energy was estimated. Even though they did not mention the growth time for the critical nuclei to be observed macroscopically, the time interval was considered to be included in the induction time they denoted. Then in 1957, Dunning and Notley32 used a double-pulse method to investigate the nucleation rate of crystalline cyclonite from aqueous acetone in a jacketed vessel. The number of crystals were counted by direct microscopic examination or indirectly integration of the size distribution curves of the crystals. Results from different temperatures and binary solvent mixtures indicated the nature change of the critical nuclei. Galkin and Vekilov38,65 applied the double-pulse method to the direct measurement of nucleation rates of lysozyme. Teflon well plate containing 400 droplets of 0.7 μL was used to conduct the kinetics determination experiments. Similar with the droplet based method, nucleation events showed stochasticity in small volume solutions. Mean crystal number for a certain nucleation time was used to calculate the nucleation rates for different concentrations of protein and precipitant. Results showed that protein crystal nucleation may occur at or even beyond the boundary of applicability of classical, continuum nucleation models. Nanev et al used a quasi-two-dimensional cell, shown in Fig. 666, to study the effect of different surface templates on lysozyme nucleation kinetics based on the double-pulse methodology. Steady state nucleation rates were calculated from the number of nuclei vs. time plots. It was found that hydrophobic surfaces (coated with hexamethyl-disilazane)63 could promote nucleation rate while templates of ploy-L-lysine67 could suppress the nucleation rate, when compared with bare glass cells. Besides, they modified the experiment method by adding an ultrasonic field to the glass cell and found that the external ultrasonic field doubled the nucleation rate of lysozyme, when compared to conditions without external field68. In 2004, Nanev et al69 modified the classical double-pulse method by using a supersaturation gradient (generated by a temperature gradient) along an insulin solution contained in a glass capillary tube, and 15
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
successfully obtained nucleation rates at different supersaturations simultaneously. This method turned out to be time, labor and material saving. Furthermore, using the same method, they also measured the crystal nucleation rate at different places of the cell respectively: bulk solution, glass support, air/solution interface, solution/glass/air boundary for insulin and lysozyme. Interestingly, it was found that no true homogeneous nucleation but heterogenous (surface) crystal nucleation happened even in the bulk solution, where biomaterial presented in the protein solutions could serve as foreign particles62,66.
Fig.6 Quasi two-dimensional cell for crystallization investigations. It consists of two optically parallel glass plates (1) and two inlets (2). The protein solution is injected in the gap (δ) between the plates (1). This gap was changed stepwise. (Reproduced with permission from ref 66. Copyright 2010 John Wiley and Sons.)
Double-pulse method detaches the crystal nucleation and growth processes ingeniously, avoiding the observation of invisible nuclei and the systematic error caused by ignoring the growth time. However, the direct measurement of nucleation rate by counting crystals is a labor and time consuming work and it also needs the growth rate to be slow enough to apply the separation strategy38. Furthermore, the determination of experiment conditions is not as easy as the theory62,66. The growth temperature needs to make sure that no new nuclei will form and the original nuclei can only grow, even though the Ostwald ripen may happen inevitably to dissolve some of the crystals, which makes the nucleation rate underestimated. Different solution concentrations may need different growth conditions. To count the crystal number clearly and accurately, the crystals should be separated from each other, which limits the supersaturation range we could study. All of these problems need to 16
ACS Paragon Plus Environment
Page 16 of 36
Page 17 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
be solved through preliminary experiments, which add more work to researchers.
Microfluidic method The emerging of high throughput microfluidic technique provides a new method for crystallization process investigation, such as polymorph or cocrystal screening, single crystal preparation, nanocrystal preparation, nucleation or crystallization kinetics, etc.70–75 It has been used to crystalize a variety of inorganic and organic compounds, including water, salts, amino acids, proteins and active pharmaceutical ingredients. Apart from the good control of transport phenomena (enhanced mass and heat transfer), little gravity effect, confinement individual space and few impurities, the large quantities of nearly monodisperse droplets make microfluidic a convenient and high-efficient technique to conduct a statistical analysis of stochastic nucleation process, with higher accuracy. 41,76,77 Laval et al78 developed a microfluidic device (shown in Fig. 7) which allows us to produce monodisperse droplets (~ 100 nL) and control their temperatures as they flow. With specific channel geometries, continuous droplets of potassium nitrate (KNO3) dispersed in silicone oil were formed and flowed along the channel. By using a stereo microscope and a CCD camera, movies of observation areas under crossed polarizers were made to estimate the probability of droplets with crystals at different nucleation times (corresponding to different chip area for the flowing droplets). Nucleation rates at different supersaturations and temperatures were derived from the probability temporal evolution curve. The results showed good agreement with the classical nucleation theory.
17
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Fig.7 Design of the microdevice. The oil and aqueous phases are injected, respectively, in inlets 1 and 2. Oil can also be injected in inlet 3 to increase the velocity of the droplets. The top-right image shows the formation of the droplets. The two dashed areas are temperature controlled at T1 and T2. The dotted area in the middle of the serpentine corresponds to the observation zone through a stereo microscope. (Reproduced with permission from ref 78. Copyright 2007 Elsevier.)
Gong et al77 applied microfluidic device to measure the nucleation rate of a thermoresponsive colloidal poly-N-isopropylacrylamide (PNIPAM) system. The colloidal PNIPAM suspension was injected into a microfluidic flow-focusing device to generate monodispersed droplets in oil. Temperature control can fine tune the volume fraction of PNIPAM particles, while the flow rate can be changed to control the volume of droplets. An optical microscope with a pair of crossed polarizers was employed to detect the nucleation and the percentage of crystallized droplets (shown in Fig. 8) within the investigated area over time. Therefore, the nucleation rate could be calculated from the linear line of crystallized droplets percentage versus time. Results showed that, in bigger volume droplets, multiple crystallites existed and interacted with each other, which affected the accuracy of the measurement of 18
ACS Paragon Plus Environment
Page 18 of 36
Page 19 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
nucleation rate.
Fig.8 Emulsion crystallization observed between a pair of crossed polarizers on an optical microscope. (a) PNIPAM droplets (500 µm size) at 23.6 °C. (b) PNIPAM droplets (100 µm size) at 23.6 °C. The multiple crystallites in the larger droplet system indicate interactions among nuclei at high nucleation volume. Arrows indicate the crystallized droplets in the smaller droplet system. (Reproduced with permission from ref 77. Copyright 2007 American Chemical Society.)
Microfluidics could also be used to study the nucleation kinetics of evaporation crystallization and precipitation process. Braatz et al74,79 determined the nucleation kinetics of paracetamol and glycine crystals under conditions of time-varying supersaturation in a droplet-based microfluidics, named as Evaporation-based 19
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
microwell crystallization platform. Vitry et al80 developed a tailored microfluidic device which could synchronously generate drops of two different reagents and could ensure the total mixing process to happen within few milliseconds. By carefully design experimental conditions and apparatus parameters, high nucleation rates of neodymium oxalate crystals formed from reaction between neodymium nitrate and oxalic acid were measured successfully. Of course, the double-pulse method can be used together with the microfluidics device. Ildefonso et al76 used a simply constructed microfluidics setup to measure the nucleation kinetics of lysozyme, shown in Fig. 9. The nucleation time dependence of mean number of crystals in 250 nL droplets was used to estimate the steady state nucleation rate. Compared to the work of Galkin and Vekilov38, the diminished volume of solution widened the range of measurable supersaturation or nucleation rate. Rossi et al41 employed the double-pulse method to investigate the nucleation rate of adipic acid solutions by using a microfluidic device under both stagnant and flowing conditions. Instead of the number of crystals, the number of nucleated droplets was counted for different length of nucleation time to calculate the nucleation rate. It was found that nucleation was enhanced in moving droplets. It was considered to be resulted from the internal circulation in the flowing droplets, which promoted the attachment frequencies of molecules to the nuclei.
20
ACS Paragon Plus Environment
Page 20 of 36
Page 21 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
Fig.9 (a) Image of the device, (b) zoom of the inlets of the plug factory showing droplet formation, (c) image of the stored droplets, and (d−f) examples of lysozyme droplets observed after 20 h, at Tnucleation = 20 °C and Tgrowth = 40 °C (channel width 500 µm). (Reproduced with permission from ref 76. Copyright 2011 American Chemical Society.)
Microfluidic technique perfectly solves the size inconformity problem of the droplet method. It also has the advantages of large data quantity and material saving, which is helpful to accurate measurement of the nucleation rates. Therefore, it has been widely used to successfully depict the classical nucleation theory expression. However, this method also has some drawbacks81, mainly associated to the design and manufacture of microfluidic device82–84, such as the device material need to be transparent and compatible to the measured system, the temperature control and measurement strategy is hard for temperature-varying device, the small diameter of the pipe containing the droplets makes clog happen easily, the cleaning of the system need to be careful. Also, the massive data processing is time and labor consuming. Without combining with the double-pulse method, microfluidic method alone fails to include the interference of crystal growth time.
21
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Stirred small volume solutions Droplet based method, double pulse method and microfluidic methods mentioned above have been widely used to gather crystal nucleation kinetics, and could be related to the classical nucleation theory or two-step nucleation theory expressions to dig the nucleation mechanism from molecule viewpoint. However, the kinetic parameters cannot be used in the industrial manufacturing process, due to the lack of stirring process. As we all know, stirring is necessary and inevitable in crystallization process to achieve good mass and heat transfer and can cause secondary nucleation which will enhance nucleation rate.40 Recently, two new methods are developed to measure the nucleation kinetics in stirred small volume solutions, mainly for cooling crystallization. These two methods make use of the stochastic nature of crystal nucleation, which are reflected by the variations in induction time (Fig. 10a) and metastable zone width (Fig. 10b), respectively.85 With multiple repeats under the same experiment conditions, nucleation rate can be characterized from the distribution of induction time or metastable zone width values.
Fig.10 (a) Two subsequent induction time measurements of the same sample of INA in ethanol at a supersaturation ratio of 1.44 (102 mg/mL) at 25°C. The measured induction times are 2191 and 2826 s. (b) Two subsequent MSZW measurements of the same sample of INA in ethanol at a concentration of 102 mg/mL with a cooling rate of 0.4 °C/min. The measured MSZW are 12 and 10°C. (Reproduced with permission from ref 85. Copyright 2013 American Chemical Society.)
22
ACS Paragon Plus Environment
Page 22 of 36
Page 23 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
Fig.11 (a) Induction time of 80 experimental data points for m-ABA solutions at the supersaturation ratio S = 1.96. (b) Probability distribution of the induction times of these 80 points. (Reproduced with permission from ref 40. Copyright 2010 American Chemical Society.)
In 2011, Jiang et al40 developed the induction time probability distribution method using medium throughput Crystal16 device (can hold 16 HPLC vials, Avantium) for the measurement of nucleation rate in stirred small volume solutions. The experiment procedure is dispersing equal amount (usually 1-1.6 mL) of bulk solutions into each HPLC vial, then setting the hold temperature above the saturation point to completely dissolve the solute for an hour or so, then crash cooling to a low temperature value to create the desired supersaturation and maintaining for hours. In-line light transmission was used to detect the onset of crystallization for each vial. Induction time was taken as the difference between the time at which the transmission started to decrease and the time at which the low temperature was reached, which means the designed supersaturation is created. By repeating the thermal cycle several times, 80 or more induction times from the same conditions could be used to calculate its probability distribution (Fig. 1140). Based on the single nucleus mechanism14,42,86, which has been observed in many stirred small volume solutions, and the assumption that the appearance of nucleus (often referred as a critical nucleus) is a random process, independent of each other and follows Poisson Distribution40, the probability curve could be fitted with a derived exponential expression, correlated to the 23
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
nucleation rate and solution volume. It is worth noting that the time elapse from the occurrence of nucleation and the detection of crystals is defined as the growth time and estimated as the shortest induction time detected in the fitting equation. Using this method, the nucleation rates of m-aminobenzoic acid (m-ABA) and L-histidine (L-His) at different supersaturations were successfully determined and it followed the trend expected from classical nucleation theory. Results indicated that heterogeneous nucleation happened and the concentration of active foreign particles was relatively low than prediction. Sullivan et al23 used the same method to determine the nucleation rates of p-aminobenzoic acid and benzoic acid in polar and nonpolar solvents and linked the results to the solution chemistry investigated by FTIR and computational chemistry. For this system, solute dimerization and desolvation were found to be rate-determining processes in the overall nucleation pathway. The nucleation rates of racemic diprophylline (DPL) stable Form RI from dimethylformamide and metastable form RII from isopropanol were also determined from the induction time distributions by Brandel et al14. The difference in nucleation behavior in the two solvents originates from the energy barrier for nucleation and the longer induction time corresponds to higher energy barrier. Mealey et al87 applied the induction time method to acquire nucleation kinetics of risperidone in seven different solvents. The FTIR spectroscopy and Density Functional Theory (DFT) were employed to investigate the solvent-solute interaction, and it was found that the stronger the solvent bound to the risperidone molecule in solution, the slower the nucleation became. In addition to the induction time, metastable zone width (MSZW) is another widely used property to characterize nucleation rate, which is also called polythermal method. Usually, MSZW is defined as the temperature difference between the supersaturation point and the nucleation point during a cooling crystallization process with a constant cooling rate. As a result of the stochastic nature of nucleation, MSZW is also not a determined parameter in small volume solutions and different MSZW corresponds to different temperature. Therefore, a MSZW distribution inherently contains the nucleation rate – temperature dependency85. Kadam et al42 performed the 24
ACS Paragon Plus Environment
Page 24 of 36
Page 25 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
MSZW measurements for paracetamol in water at 1 mL and 1 L scales to investigate the scale up of the stirred small volume nucleation kinetics. As we can see from Fig. 1242, large MSZW variations were displayed in 1 mL experiments, while almost no variations were found for 1 L scale ones. Interestingly, the measured MSZW values at 1 L scale corresponded to the onset of the MSZW distribution at 1 mL scale, implying a simple scale up rule for laboratory MSZW data to larger scales. In their another work88, MSZW data were measured for paracetamol - water and isonicotinamide ethanol systems at four different scales from 500 mL to 1 L. Together with the results from their previous work mentioned above, it showed that the spread of MSZW increased roughly inversely proportional to the volume. The nonreproducible MSZW for smaller scale was explained by single nucleus mechanism and the reproducible result for larger scale was explained by the deterministic population balance model.
Fig.12 MSZW (∆T) for paracetamol concentration of 0.0150 g/mL (a) and 0.0470 g/mL (b) at a cooling rate of 0.5 °C/min for different volumes. The solid line indicates the smallest MSZW detected with light obscuration at 1 mL while the dashed line indicates the smallest MSZW detected by visual observation at volumes between 500 25
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
mL and 1 L. (Reproduced with permission from ref 42. Copyright 2012 Elsevier.)
By using isonicotinamide - ethanol as the model system, both the induction time distribution (at constant temperature and different concentrations) and MSZW distribution (at different cooling rates and different concentrations) were measured and compared by Kulkarni et al85. Mealey et al15 also used both methods to investigate the crystal nucleation of salicylic acid in different solvents. Although the MSZW method is less labor intensive and could give the nucleation kinetic parameters from only one distribution, it failed to consider the temperature dependence of the parameters and the changing temperature could cause a temperature difference between the vials and the heating block, which is difficult to account for. In comparison, the induction time method could provide more accurate data and the data is much easier to analysis, although it is much more labor intensive. The induction time or MSZW distribution method to obtain crystal nucleation kinetics in small volume solutions saves material and reduces the experimental cost, especially for expensive pharmaceutical molecules or proteins. More importantly, it provides the potential to relate the lab scale data to the industrial scale crystallizer by introducing stirring in the measurements. However, to extract nucleation rate from such a large amount of data, a distribution curve, is kind of tricky and difficult. ter Horst et al40,89 applied a Poisson Distribution to describe the appearance of independent nuclei by assuming the supersaturation is constant during the nucleation process and used a single nucleus mechanism to illustrate the process from the onset of nucleation to the detection of crystals. Rasmuson et al15 adopted lognormal cumulative probability function and a mean value of the induction time or MSZW value to qualitatively compare or calculate the nucleation rate. Mazzotti et al90 presented two different models to describe a batch cooling crystallization process, where the interplay between the stochastic nucleation and the deterministic crystal growth was described differently. All of these data treatment methods need much work and different assumptions. Besides, in the stirred small volume solutions, 26
ACS Paragon Plus Environment
Page 26 of 36
Page 27 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
sometimes, one can only get one polymorph of the polymorphic compounds, which makes it impossible to evaluate the nucleation kinetics of other polymorphs.23 The “crowning” problem23,89, solids form on the glass vial - solution boundary, may make the detection of crystals postpone or even fail and the solution concentration could not maintain constant. Also, the crystals can arouse secondary nucleation. To guarantee the effective detection of nucleation, the amount of crystals should fall in an appropriate range.
Conclusions and outlook Given the significance of crystal nucleation, its kinetics study has remained a hot subject for decades. Five measurement methods, including deterministic method, droplet based method, double-pulse method, microfluidic method and stirred small volume solution method were reviewed in this work. For the scaling up of lab scale crystallization experiments, the deterministic method, which covers the stochastic nucleation phenomenon using its large experimental setup, could be employed to estimate the nucleation and growth kinetics. The droplet based method diminishes the concentration of catalytic centers by subdividing the system into small droplets and works well for high melting point metals. But it fails to generate droplets with uniform size. The microfluidic technique helps to solve the droplet size nonuniform problem although the design of microfluidic device need to be improved. The double-pulse method, which are widely used in the protein crystallization research, applies the idea that separating the crystal nucleation and growth by controlling the supersaturation using the metastable zone concept to make the invisible critical nuclei visible and nucleation rates can be obtained directly by counting the number of crystals in a certain time interval. As for the stirred small volume solution method, both the stochastic nature of nucleation, shown in the dispersion of measured induction time or MSZW values, and the effect of stirring on crystallization process are embodied. The latter four methods are all performed in relatively small systems or static conditions and the crystallization process is simplified to primary nucleation, 27
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
(secondary nucleation, agglomeration, breakage processes are ignored) or single nucleus mechanism. The kinetics results can help to verify the classical nucleation theory or develop new nucleation theories, helping us to better understand nucleation process mechanisms.23,65,91 Therefore, to choose an appropriate method for estimating nucleation kinetics, the system properties, the device requirements and even the purposes to get the information should be taken into account. To further understand and uncover the mystery of the nucleation process, there are still a lot of work to do in the future. Firstly, to improve the accuracy of nucleation kinetics, the analytical techniques with higher accuracy and higher resolution should be developed.92,93 If the critical nuclei could be observed or characterized more directly, nucleation rate could be obtained directly, instead of monitoring various system properties change which could cause extra uncertainty. Secondly, large amounts of experiment data need to be analyzed to estimate the nucleation kinetics, especially for the data with large dispersion. Middle value, average value, Poisson distribution assumption, lognormal cumulative probability function fitting, etc. have been used by different researchers. Qualitatively, they all can be used to compare nucleation rate. But quantitatively, which method could provide the most reliable nucleation rate remains unknown. Thirdly, for most of the nucleation rate data correlated by classical nucleation theory, the fitting result for high supersaturation points is always not as good as lower ones (model values are lower than experiment values).57 This might be caused by the uneven distribution of experimental points (more points for low supersaturation or low nucleation rate)89, or might be resulted from the intrinsic feature of nucleation, which means a tiny supersaturation change could cause huge difference in the nucleation rate for large supersaturation.64 This could also imply that the classical nucleation theory is so simplified that it cannot cover such a big nucleation rate or supersaturation range. Modified or new forms of nucleation kinetic expressions or nucleation mechanisms need to be developed to describe the nucleation process.94,95 Fourthly, in most of the crystallization processes, heterogeneous nucleation rather than homogeneous nucleation is proved to be 28
ACS Paragon Plus Environment
Page 28 of 36
Page 29 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
dominant, as a result of lower energy barrier. However, currently, we have little information about the catalytical centers in the system and this clouds further investigation of nucleation process. Nowadays, researchers have put a lot of effort on functional surface nucleation, which is an important step to understand the formation of critical nuclei.96–98 Finally, the kinetic data obtained from small volume or static methods is much easier analyzed by nucleation theories and may provide molecular mechanism of nucleation when combined with the solution chemistry or molecular simulation, while the deterministic method tends to be a good method to establish a kinetic crystallization model, which could be used for scaling up or optimizing the industrial manufacture processes. However, further connection between the small volume methods and the deterministic method is missing. With a better understanding of nucleation process, the data from different methods and different scales should be analyzed together and give reliable kinetic results. In return, the reliable nucleation rates can also help to further understand the nucleation puzzle.
Acknowledgements This research is financially supported by National Natural Science Foundation of China (No. 21376165) and National Key Research and Development Program of China (No. 2016YFB0600504.)
References
(1)
Miller, R. C.; Anderson, R. J.; Kassner, J. L.; Hagen, D. E. J. Chem. Phys. 1983, 78, 3204–3211.
(2)
Davey, R. J.; Schroeder, S. L. M.; Horst, J. H. ter. Angew. Chem. Int. Ed. 2013, 52, 2166–2179.
(3)
Davey, R. J.; Back, K. R.; Sullivan, R. A. Faraday Discuss 2015, 179, 9–26.
(4)
Vekilov, P. G. Cryst. Growth Des. 2010, 10, 5007–5019.
(5)
Erdemir, D.; Lee, A. Y.; Myerson, A. S. Acc. Chem. Res. 2009, 42, 621–629. 29
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
(6)
Dubrovskii, V. G. Nucleation theory and growth of nanostructures; Springer, Berlin, 2014.
(7)
Mullin, J. W. Crystallization (Fourth Edition); Butterworth-Heinemann: Oxford, 2001.
(8)
Davey, R. J.; Garside, J. From Molecules to Crystallizers: An Introduction to Crystallization; Oxford University Press: Oxford, UK, 2000.
(9)
Spitaleri, A.; Hunter, C. A.; McCabe, J. F.; Packer, M. J.; Cockroft, S. L. CrystEngComm 2004, 6, 489–493.
(10) Kulkarni, S. A.; McGarrity, E. S.; Meekes, H.; Horst, J. H. ter. Chem. Commun. 2012, 48, 4983–4985. (11) Parveen, S.; Davey, R. J.; Dent, G.; Pritchard, R. G. Chem. Commun. 2005, 1531–1533. (12) Chen, J.; Trout, B. L. J. Phys. Chem. B 2008, 112, 7794–7802. (13) Belenguer, A. M.; Lampronti, G. I.; Cruz-Cabeza, A. J.; Hunter, C. A.; Sanders, J. K. M. Chem. Sci. 2016, 7, 6617–6627. (14) Brandel, C.; Horst, J. H. ter. Faraday Discuss 2015, 179, 199–214. (15) Mealey, D.; Croker, D. M.; Rasmuson, Å. C. CrystEngComm 2015, 17, 3961–3973. (16) Alexandrov, D. V.; Malygin, A. P. J. Phys. Math. Theor. 2013, 46, 455101. (17) Shneidman, V. A. Phys. Rev. E 2010, 82, 031603. (18) Alexandrov, D. V.; Nizovtseva, I. G. Proc. R. Soc. A 2013, 470, 20130647. (19) Barlow, D. A. J. Cryst. Growth 2009, 311, 2480–2483. (20) Barlow, D. A. J. Cryst. Growth 2017, 470, 8–14. (21) Kashchiev, D.; van Rosmalen, G. M. Cryst. Res. Technol. 2003, 38, 555–574. (22) Schmitt, J. L.; Adams, G. W.; Zalabsky, R. A. J. Chem. Phys. 1982, 77, 2089–2097. (23) Sullivan, R. A.; Davey, R. J.; Sadiq, G.; Dent, G.; Back, K. R.; Horst, J. H. ter; Toroz, D.; Hammond, R. B. Cryst. Growth Des. 2014, 14, 2689–2696. (24) Vekilov, P. G. Cryst. Growth Des. 2004, 4, 671–685. 30
ACS Paragon Plus Environment
Page 30 of 36
Page 31 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
(25) Karthika, S.; Radhakrishnan, T. K.; Kalaichelvi, P. Cryst. Growth Des. 2016, 16, 6663–6681. (26) Cruz-Cabeza, A. J.; Davey, R. J.; Sachithananthan, S. S.; Smith, R.; Tang, S. K.; Vetter, T.; Xiao, Y. Chem. Commun. 2017, 53, 7905–7908. (27) Jiang, S. Crystallization Kinetics in Polymorphic Organic Compounds; Enschede, The Netherlands, 2009. (28) Kubota, N. J. Cryst. Growth 2010, 312, 548–554. (29) Camacho Corzo, D. M.; Borissova, A.; Hammond, R. B.; Kashchiev, D.; Roberts, K. J.; Lewtas, K.; More, I. CrystEngComm 2014, 16, 974–991. (30) Kashchiev, D.; Borissova, A.; Hammond, R. B.; Roberts, K. J. J. Cryst. Growth 2010, 312, 698–704. (31) Shneidman, V. A. J. Chem. Phys. 2014, 141, 051101. (32) Dunning, W. J.; Notley, N. T. Z. Für Elektrochem. Berichte Bunsenges. Für Phys. Chem. 1957, 61, 55–59. (33) Cheng, Z.; Chaikin, P. M.; Zhu, J.; Russel, W. B.; Meyer, W. V. Phys. Rev. Lett. 2002, 88, 015501. (34) Dugua, J.; Simon, B. J. Cryst. Growth 1978, 44, 265–279. (35) Rasmussen, D. H.; Loper. Jr., C. R. Acta Metall. 1976, 24, 117–123. (36) Kim, J.-W.; Kim, J.; Lee, K.-D.; Koo, K.-K. Cryst. Res. Technol. 2013, 48, 1097–1105. (37) Tsekova, D. S. Cryst. Growth Des. 2009, 9, 1312–1317. (38) Galkin, O.; Vekilov, P. G. J. Phys. Chem. B 1999, 103, 10965–10971. (39) Knezic, D.; Zaccaro, J.; Myerson, A. S. J. Phys. Chem. B 2004, 108, 10672–10677. (40) Jiang, S.; Horst, J. H. ter. Cryst. Growth Des. 2011, 11, 256–261. (41) Rossi, D.; Gavriilidis, A.; Kuhn, S.; Candel, M. A.; Jones, A. G.; Price, C.; Mazzei, L. Cryst. Growth Des. 2015, 15, 1784–1791. (42) Kadam, S. S.; Kulkarni, S. A.; Coloma Ribera, R.; Stankiewicz, A. I.; Horst, J. H. ter; Kramer, H. J. M. Chem. Eng. Sci. 2012, 72, 10–19. 31
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
(43) Bhamidi, V.; Skrzypczak-Jankun, E.; Schall, C. A. J. Cryst. Growth 2001, 232, 77–85. (44) Bhamidi, V.; Varanasi, S.; Schall, C. A. Cryst. Growth Des. 2002, 2, 395–400. (45) Wunderlich, B. Crystal Nucleation, Growth, Annealing; Academic Press: London, 1979; Vol. 2. (46) Tang, S.; Hu, Z.; Cheng, Z.; Wu, J. Langmuir 2004, 20, 8858–8864. (47) Liu, X.; Wang, Z.; Duan, A.; Zhang, G.; Wang, X.; Sun, Z.; Zhu, L.; Yu, G.; Sun, G.; Xu, D. J. Cryst. Growth 2008, 310, 2590–2592. (48) Zhao, L.; Zhu, C.; Ji, J.; Chen, J.; Teng, H. H. Geochim. Cosmochim. Acta 2013, 106, 192–202. (49) Shiau, L.-D. CrystEngComm 2016, 18, 6358–6364. (50) Trifkovic, M.; Sheikhzadeh, M.; Rohani, S. Ind. Eng. Chem. Res. 2008, 47, 1586–1595. (51) Nývlt, J. J. Cryst. Growth 1968, 3-4, 377–383. (52) Nagy, Z. K.; Fujiwara, M.; Woo, X. Y.; Braatz, R. D. Ind. Eng. Chem. Res. 2008, 47, 1245–1252. (53) Trifkovic, M.; Rohani, S.; Sheikhzadeh, M. J. Cryst. Process Technol. 2012, 02, 31–43. (54) Li, H.; Grover, M. A.; Kawajiri, Y.; Rousseau, R. W. Chem. Eng. Sci. 2013, 89, 142–151. (55) Agimelen, O. S.; Jawor-Baczynska, A.; McGinty, J.; Dziewierz, J.; Tachtatzis, C.; Cleary, A.; Haley, I.; Michie, C.; Andonovic, I.; Sefcik, J.; Mulholland, A. J. Chem. Eng. Sci. 2016, 144, 87–100. (56) Unno, J.; Hirasawa, I. Adv. Chem. Eng. Sci. 2017, 07, 91–107. (57) Vonnegut, B. J. Colloid Sci. 1948, 3, 563–569. (58) Turnbull, D. J. Chem. Phys. 1952, 20, 411–424. (59) White, M. L.; Frost, A. A. J. Colloid Sci. 1959, 14, 247–251. (60) Keller, D. M.; Massey, R. E.; Hileman, JR., O. E. Can. J. Chem. 1978, 56, 831–838. 32
ACS Paragon Plus Environment
Page 32 of 36
Page 33 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
(61) Filipponi, A.; Malvestuto, M. Meas. Sci. Technol. 2003, 14, 875–882. (62) Nanev, C. N.; Hodzhaoglu, F. V.; Dimitrov, I. L. Cryst. Growth Des. 2011, 11, 196–202. (63) Tsekova, D.; Dimitrova, S.; Nanev, C. N. J. Cryst. Growth 1999, 196, 226–233. (64) Dunning, W. J.; Shipman, A. J. Proc Agric Ind. 10th Int. Conf. 1954, 1448–1456. (65) Galkin, O.; Vekilov, P. G. J. Am. Chem. Soc. 2000, 122, 156–163. (66) Hodzhaoglu, F. V.; Nanev, C. N. Cryst. Res. Technol. 2010, 45, 281–291. (67) Nanev, C. N.; Tsekova, D. Cryst. Res. Technol. 2000, 35, 189–195. (68) Nanev, C. N.; Penkova, A. J. Cryst. Growth 2001, 232, 285–293. (69) Penkova, A.; Dimitrov, I.; Nanev, C. Ann. N. Y. Acad. Sci. 2004, 1027, 56–63. (70) Frenz, L.; Harrak, A. El; Pauly, M.; Bégin-Colin, S.; Griffiths, A. D.; Baret, J.-C. Angew. Chem. Int. Ed. 2008, 47, 6817–6820. (71) Song, Y.; Modrow, H.; Henry, L. L.; Saw, C. K.; Doomes, E. E.; Palshin, V.; Hormes, J.; Kumar, C. S. S. R. Chem. Mater. 2006, 18, 2817–2827. (72) Zheng, B.; Roach, L. S.; Ismagilov, R. F. J. Am. Chem. Soc. 2003, 125, 11170–11171. (73) Sanjoh, A.; Tsukihara, T. J. Cryst. Growth 1999, 196, 691–702. (74) Goh, L.; Chen, K.; Bhamidi, V.; He, G.; Kee, N. C. S.; Kenis, P. J. A.; Zukoski, C. F.; Braatz, R. D. Cryst. Growth Des. 2010, 10, 2515–2521. (75) Shi, H.; Xiao, Y.; Ferguson, S.; Huang, X.; Wang, N.; Hao, H. Lab Chip 2017, 17, 2167–2185. (76) Ildefonso, M.; Candoni, N.; Veesler, S. Cryst. Growth Des. 2011, 11, 1527–1530. (77) Gong, T.; Shen, J.; Hu, Z.; Marquez, M.; Cheng, Z. Langmuir 2007, 23, 2919–2923. (78) Laval, P.; Salmon, J.-B.; Joanicot, M. J. Cryst. Growth 2007, 303, 622–628. (79) Chen, K.; Goh, L.; He, G.; Kenis, P. J. A.; Zukoski, C. F.; Braatz, R. D. Chem. Eng. Sci. 2012, 77, 235–241. 33
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
(80) Vitry, Y.; Teychené, S.; Charton, S.; Lamadie, F.; Biscans, B. Chem. Eng. Sci. 2015, 133, 54–61. (81) Stan, C. A.; Schneider, G. F.; Shevkoplyas, S. S.; Hashimoto, M.; Ibanescu, M.; Wiley, B. J.; Whitesides, G. M. Lab. Chip 2009, 9, 2293–2305. (82) Yu, Y.; Wang, X.; Oberthür, D.; Meyer, A.; Perbandt, M.; Duan, L.; Kang, Q. J. Appl. Crystallogr. 2012, 45, 53–60. (83) Song, K.-Y.; Zhang, W.-J.; Gupta, M. M. J. Manuf. Sci. Eng. 2012, 134, 044504. (84) Tanyeri, M.; Ranka, M.; Sittipolkul, N.; Schroeder, C. M. Lab. Chip 2011, 11, 1786–1794. (85) Kulkarni, S. A.; Kadam, S. S.; Meekes, H.; Stankiewicz, A. I.; Horst, J. H. ter. Cryst. Growth Des. 2013, 13, 2435–2440. (86) Kulkarni, S. A.; Meekes, H.; Horst, J. H. ter. Cryst. Growth Des. 2014, 14, 1493–1499. (87) Mealey, D.; Zeglinski, J.; Khamar, D.; Rasmuson, Å. C. Faraday Discuss 2015, 179, 309–328. (88) Kadam, S. S.; Kramer, H. J. M.; Horst, J. H. ter. Cryst. Growth Des. 2011, 11, 1271–1277. (89) Xiao, Y.; Tang, S. K.; Hao, H.; Davey, R. J.; Vetter, T. Cryst. Growth Des. 2017, 17, 2852–2863. (90) Maggioni, G. M.; Mazzotti, M. Faraday Discuss 2015, 179, 359–382. (91) Khamar, D.; Zeglinski, J.; Mealey, D.; Rasmuson, Å. C. J. Am. Chem. Soc. 2014, 136, 11664–11673. (92) Tsoutsouva, M. G.; Duffar, T.; Garnier, C.; Fournier, G. Cryst. Res. Technol. 2015, 50, 55–61. (93) Nagy, Z. K.; Fevotte, G.; Kramer, H.; Simon, L. L. Chem. Eng. Res. Des. 2013, 91, 1903–1922. (94) Slezov, V. V. Kinetics of First-Order Phase Transitions; Wiley-VCH: Weinheim, 2009. 34
ACS Paragon Plus Environment
Page 34 of 36
Page 35 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Crystal Growth & Design
(95) Kelton, K. F.; Greer, A. L. Nucleation in condensed matter: applications in materials and biology; Elsevier: Amsterdam, 2010. (96) Diao, Y.; Helgeson, M. E.; Myerson, A. S.; Hatton, T. A.; Doyle, P. S.; Trout, B. L. J. Am. Chem. Soc. 2011, 133, 3756–3759. (97) Cui, Y.; Myerson, A. S. Cryst. Growth Des. 2014, 14, 5152–5157. (98) Diao, Y.; Myerson, A. S.; Hatton, T. A.; Trout, B. L. Langmuir 2011, 27, 5324–5334.
35
ACS Paragon Plus Environment
Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
For Table of Contents Use Only
Determination methods for crystal nucleation kinetics in solutions Yan Xiao a,b, Jingkang Wang a,b, Xin Huang a,b, Huanhuan Shi a,b, Yanan Zhou a,b, Shuyi Zong a,b, Hongxun Hao a,b *, Ying Bao a,b, Qiuxiang Yin a,b
Because of the tiny size of critical nucleus and the limitation of analytical technique threshold, detection of crystal nucleation process could not be done directly and instantly. Therefore, indirect methods including deterministic method, droplet based method, double-pulse method, microfluidic method and stirred small volume solution method, have been developed to estimate the nucleation kinetics for different systems and purposes.
36
ACS Paragon Plus Environment
Page 36 of 36