Hit and miss of classical nucleation theory as revealed by a molecular simulation study of crystal nucleation in supercooled sulfur hexafluoride

Classical nucleation theory pictures the homogeneous nucleation of a crystal as the formation of a spherical crystalline embryo, possessing the properties of the macroscopic crystal, inside a parent supercooled liquid. In this work we study crystal nucleation in moderately supercooled sulfur hexafluoride by umbrella sampling simulations. The nucleation free energy evolves from 5.2kBT at T=170 K to 39.1kBT at T=195 K. The corresponding critical nucleus size ranges from 40 molecules at T=170 K to 266 molecules at T=195 K. Both nucleation free energy and critical nucleus size are shown to evolve with temperature according to the equations derived from the classical nucleation theory. Inspecting the obtained nuclei we show, however, that they present quite anisotropic shapes in opposition to the spherical assumption of the theory. Moreover, even though the critical nuclei possess the structure of the stable bcc plastic phase, the only mechanically stable crystal phase for SF6 in the temperature range investigated, they are shown to be less ordered than the corresponding macroscopic crystal. Their crystalline order is nevertheless shown to increase regularly with their size. This is confirmed by a study of a nucleus growth from a critical size to a size of the order of 10(4) molecules. Similarly to the fact that it does not affect the temperature dependence of the nucleation free energy and of the critical nucleus size, the ordering of the nucleus with size does not affect the growth rate of the nucleus.     

Nucleation of water vapor in microcracks on the surface of β-AgI aerosol particles: 2. Thermodynamics of nucleation

Free energy, entropy, and the work of formation of condensation nuclei at 260 K in microcracks of β-AgI crystal structure at the initial stage of nucleation preceding crystallization are calculated by the Monte Carlo method. Unlike ideal crystal surface, nuclei in microcracks are thermodynamically stable and the barrier of free energy of nucleation is absent. Conditions of microcrack are favorable for the crystallization that qualitatively changes the regime and rate of nucleation. Stable size of nuclei at the humidity corresponding to natural atmosphere is sufficient for the filling of nanoscopic microcracks and the attainment of substrate surface. The probability of nucleus formation in microcracks by the fluctuation mechanism is incomparably higher than the probability of their formation on the defect-free surface. High crystallization ability of the particles of βAgI aerosol is ensured by multiple surface microcracks acting as active sites in combination with its complementary crystal structure. The efficiency of aerosols as stimulants of the nucleation of water vapor at negative Celsius temperatures is determined by the surface density and geometry of nanoscopic cracks and fissures on the particle surface.     

磁场影响聚合物本体结晶动力学的机理分析

在经典的热力学理论基础上,探讨了磁场对聚合物本体结晶过程的成核与生长的影响,建立了相关结晶动力学理论方程.初步认为,磁场产生的＂磁结晶效应＂可能是由于晶相与非晶相之间磁化率差异导致了两相之间磁化能的差异,也可能由于聚合物体系在结晶前会形成一种有序相,减小了体系的熵值,进而改变了结晶过程中的体系自由能,影响其成核与晶体生长,乃至整个结晶动力学方程.利用Matlab软件结合PLLA的各结晶参数值,绘制了结晶自由能与各成核临界参数之间的函数图像.结果表明,在低过冷度下,较小的自由能扰动可能导致较大的晶核临界参数变化.     

Atomistic simulation of the homogeneous nucleation and of the growth of N2 crystallites

We report on a computer simulation study of the early stages of the crystallization of molecular nitrogen. First, we study how homogeneous nucleation takes place in supercooled liquid N(2) for a moderate degree of supercooling. Using the umbrella sampling technique, we determine the free energy barrier of formation for a critical nucleus of N(2). We show that, in accord with Ostwald's rule of stages, the structure of the critical nucleus is predominantly that of a metastable polymorph (alpha-N(2) for the state point investigated). We then monitor the evolution of several critical nuclei through a series of unbiased molecular dynamics trajectories. The growth of N(2) crystallites is accompanied by a structural evolution toward the stable polymorph beta-N(2). The microscopic mechanism underlying this evolution qualitatively differs from that reported previously. We do not observe any dissolution or reorganization of the alpha-like core of the nucleus. On the contrary, we show that alpha-like and beta-like blocks coexist in postcritical nuclei. We relate the structural evolution to a greater adsorption rate of beta-like molecules on the surface and show that this transition actually starts well within the precritical regime. We also carefully investigate the effect of the system size on the height of the free energy barrier of nucleation and on the structure and size of the critical nucleus.     

Crystal nucleation of colloidal hard dumbbells

Using computer simulations, we investigate the homogeneous crystal nucleation in suspensions of colloidal hard dumbbells. The free energy barriers are determined by Monte Carlo simulations using the umbrella sampling technique. We calculate the nucleation rates for the plastic crystal and the aperiodic crystal phase using the kinetic prefactor as determined from event driven molecular dynamics simulations. We find good agreement with the nucleation rates determined from spontaneous nucleation events observed in event driven molecular dynamics simulations within error bars of one order of magnitude. We study the effect of aspect ratio of the dumbbells on the nucleation of plastic and aperiodic crystal phases, and we also determine the structure of the critical nuclei. Moreover, we find that the nucleation of the aligned close-packed crystal structure is strongly suppressed by a high free energy barrier at low supersaturations and slow dynamics at high supersaturations.     

Supercooled states in the benzene-naphthalene system

Precrystallization supercooled states in the benzene-naphthalene system were studied by thermal analysis methods. A sharp change in the character of the crystallization of benzene and naphthalene, from quasi-equilibrium without supercooling to nonequilibrium-explosive with supercooling depending on liquid phase superheating with respect to the temperature of fusion was observed. Such a transition occurred monotonically for alloys. The phase diagram of the benzene-naphthalene system including metastable regions with supercooling under normal crystallization conditions was constructed. The rules governing Gibbs energy changes and mean supercoolings in crystallization depending on the concentrations of the components were determined. The supercooling values were also used to calculate such nucleation parameters as the size of critical nuclei, work of their formation, and number of unit cells and molecules in them. The molecular structure of benzene and naphthalene crystals and the structural changes of these substances during fusion that influenced the character of their subsequent crystallization were described.     

CO32-对羟基乙叉二膦酸体系铜电结晶的影响

采用循环伏安和计时安培法研究了羟基乙叉二膦酸(HEDPA)镀铜液中铜在玻碳电极上电结晶的初期行为。结果表明：羟基乙叉二膦酸(HEDPA)镀铜体系中,铜的电沉积过程经历了晶核形成过程;当溶液中不含CO32-时,其电结晶按连续三维成核方式进行,而CO32-的加入,使得铜电结晶按瞬时三维成核方式进行;成核数密度都随着电位的提高而增加。这可能是CO32- -以第二配体形式进入HEDPA和Cu2+构成的络合结构,从而形成更稳定的络合物吸附在电极表面所致。     

Surface area controlled heterogeneous nucleation

Heterogeneous nucleation of liquid from a gas phase on nanoparticles has been studied under various saturation ratios and nuclei size. The probability of liquid droplet nucleation, especially at a low degree of deviation from equilibrium, was measured for both atmospheric aerosol particles and engineered nanoparticles Cr(2)O(3). The concept of a critical saturation ratio and the validity of the one-to-one relationship between the nuclei number and the number of droplets were examined. A transient zone between no nucleation and established nucleation termed the surface area controlled nucleation was observed. In this zone, the probability of stable phase formation is determined by the surface area of nuclei. There are two distinctive features of the surface area controlled nucleation: the nucleation probability is much less than 1 and is proportional to the surface area of nuclei. For condensation particle counters (CPCs) counting nanoparticles, these features mean that counts measured are proportional to the surface area of nanoparticles and, therefore, the CPCs counts can be calibrated to measure the surface area.     

A simulation test of the optical Kerr mechanism for laser-induced nucleation

Recent experiments have demonstrated that intense, nanosecond laser pulses can induce crystal nucleation from supersaturated solutions that are transparent at the incident wavelengths, a phenomenon termed nonphotochemical laser-induced nucleation (NPLIN). Previous work has proposed that this effect is due to the alignment of solute molecules in solution due to the electric field of the applied laser light, promoting crystalline order. We have used simulations of NPLIN to examine how an orientational bias in solution affects nucleation with Monte Carlo simulations of a Potts lattice gas model. We examine this effect within both a classical, one-step nucleation framework as well as in the context of two-step nucleation. Our results indicate that an orientational bias can reduce the free energy barrier to nucleation within the one-step picture as well as promote the crystallization of amorphous precritical nuclei (the rate-determining step in the two-step picture). However, these effects are only present with field strengths that are much greater than those used in experiments.     

Formation of crystal nuclei near critical supersaturation in small volumes

This work deals with the nucleation of crystals in confined systems in response to the recent high interest in research on crystallization in emulsion and microemulsion droplets. In these confined systems, crystallization often occurs at high supercooling; thus, nucleation determines the overall crystallization process. A decrease in the volume of the confined mother phase leads to the higher supercooling needed for the phase transition. We have numerically solved kinetic equations in order to determine the conditions under which the first crystal nuclei are formed by homogeneous and heterogeneous nucleation from supercooled melt and supersaturated solution, depending on the volume of the mother phase. Supersaturation (or supercooling) increases with decreasing volume of the mother phase. The nucleation barrier depends linearly on the logarithm of volume of the mother phase in all cases under consideration, as follows from the numerical solution of kinetic equations.     

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: II. Nucleation in the metastable liquid immiscibility region

In the second part of our paper, we address crystal nucleation in the metastable liquid miscibility region of eutectic systems that is always present, though experimentally often inaccessible. While this situation resembles the one seen in single component crystal nucleation in the presence of a metastable vapor-liquid critical point addressed in previous works, it is more complex because of the fact that here two crystal phases of significantly different compositions may nucleate. Accordingly, at a fixed temperature below the critical point, six different types of nuclei may form: two liquid-liquid nuclei: two solid-liquid nuclei; and two types of composite nuclei, in which the crystalline core has a liquid "skirt," whose composition falls in between the compositions of the solid and the initial liquid phases, in addition to nuclei with concentric alternating composition shells of prohibitively high free energy. We discuss crystalline phase selection via exploring/identifying the possible pathways for crystal nucleation.     

Nucleation kinetics of folded chain crystals of polyethylene on active centers

Kinetic equations describing formation of nuclei of a new phase on active centers from supercooled melt taking into account exhaustion of active centers are solved numerically. Basic characteristics of nucleation process (total number of nuclei and nucleation rate) of folded chain crystals of polyethylene at low supercooling are determined and compared with measured data. Our model gives good coincidence with experimental measurements of the total number density of nuclei. Nucleation rate reaches some quasistationary limit at sufficiently long time, which is approximately 20% of the stationary nucleation rate determined by the standard way.     

Thermodynamics of crystalline nanonucleus formation from liquid phase

The energy of crystal nucleation from liquid phase was considered, with the following two stages taken into account: (1) the formation of metastable supercooled melt (solution), containing pre-nuclei with intermediate amorphous (quasicrystalline) structure, and (2) the transformation of amorphous clusters into solid crystalline nuclei having different structures. With growth of a nucleus the nucleation energy profile manifests 2–3 maxima corresponding to these stages, and the kinetics of the non-stationary nucleation has five characteristic variations.     

Rate of homogeneous crystal nucleation in molten NaCl

We report a numerical simulation of the rate of crystal nucleation of sodium chloride from its melt at moderate supercooling. In this regime nucleation is too slow to be studied with "brute force" molecular-dynamics simulations. The melting temperature of ("Tosi Fumi") NaCl is approximately 1060 K. We studied crystal nucleation at T = 800 and 825 K. We observe that the critical nucleus formed during the nucleation process has the crystal structure of bulk NaCl. Interestingly, the critical nucleus is clearly faceted, the nuclei have a cubical shape. We have computed the crystal-nucleation rate using two completely different approaches, one based on an estimate of the rate of diffusive crossing of the nucleation barrier, the other based on the forward flux sampling and transition interface sampling methods. We find that the two methods yield the same result within an order of magnitude. However, when we compare the extrapolated simulation data with the only available experimental results for NaCl nucleation, we observe a discrepancy of nearly five orders of magnitude. We discuss the possible causes for this discrepancy.     

Structure prediction, disorder and dynamics in a DMSO solvate of carbamazepine

We have applied crystal structure prediction methods to understand and predict the formation of a DMSO solvate of the anti-convulsant drug carbamazepine (CBZ), in which the DMSO molecules are disordered. Crystal structure prediction calculations on the 1:1 CBZ:DMSO solvate revealed the generation of two similar low energy structures which differ only in the orientation of the DMSO molecules. Analysis of crystal energy landscapes generated at 0 K suggests the possibility of solvent disorder. A combined computational and experimental study of the changes in the orientation of the DMSO within the crystal structure revealed that the nature of the disorder changes with temperature. At low temperature, the DMSO disorder is static whilst at high temperature the DMSO configurations can interconvert by a 180° rotation of the DMSO molecules within the lattice. This 180° rotation of the DMSO molecules drives a phase change from a high temperature dynamically disordered phase to a low temperature phase with static disorder. Crystallisation of a DMSO solvate of the related molecule epoxycarbamazepine resulted in a different degree of DMSO disorder in the crystal structure, despite the similarity of the carbamazepine and epoxycarbamazepine molecules. We believe consideration of disorder and its contribution to entropy and crystal free energies at temperature other than 0 K is fundamental for the accuracy of future energy rankings in crystal structure prediction calculations of similar solvated structures.     

Order-parameter-based Monte Carlo simulation of crystallization

A Monte Carlo simulation method is presented for simulation of phase transitions, with emphasis on the study of crystallization. The method relies on a random walk in order parameter Phi(q(N)) space to calculate a free energy profile between the two coexisting phases. The energy and volume data generated over the course of the simulation are subsequently reweighed to identify the precise conditions for phase coexistence. The usefulness of the method is demonstrated in the context of crystallization of a purely repulsive Lennard-Jones system. A systematic analysis of precritical and critical nuclei as a function of supercooling reveals a gradual change from a bcc to a fcc structure inside the crystalline nucleus as it grows at large degrees of supercooling. The method is generally applicable and is expected to find applications in systems for which two or more coexisting phases can be distinguished through one or more order parameters.     

Polymorph selection during the crystallization of softly repulsive spheres: the inverse power law potential

Using hybrid Monte Carlo molecular simulations, we study crystallization from the melt of softly repulsive spheres interacting through an inverse power law potential. We work at fixed supercooling (i.e., at a temperature 25% below the melting temperature) and consider three systems, defined by different values for the inverse power exponent n: n = 5, n = 6.67, and n = 10. Modifying the value of n allows us to study the onset of crystallization in the domain of stability of the body-centered cubic (bcc) phase (n = 5 and n = 6.67) and in the domain of stability of the face-centered cubic (fcc) phase (n = 10). We show that, for the three systems, polymorph selection does not take place during crystal nucleation since the structure of the critical nuclei obtained for the three systems is not well defined. However, our results demonstrate that polymorph selection takes place during the growth step since growth proceeds either into the stable bcc phase for the two smaller values of n (n = 5 and n = 6.67) or into the stable fcc phase for the larger value of n (n = 10). We also show that we did not achieve complete control of polymorphism for n = 10. The growth step gives rise to either slowly growing crystallites composed of two blocks of different structures (the stable fcc form and the metastable bcc form) or rapidly growing crystallites of the metastable bcc form.     

Reorganization and growth of metastable alpha-N2 critical nuclei into stable beta-N2 crystals

We report on results on the crystal nucleation and growth of nitrogen. Using molecular dynamics simulations, we show that while nucleation proceeds into the metastable alpha-phase (i.e., the crystalline phase associated with the lowest free energy barrier of formation), growth of the crystallite proceeds through a reorganization of the nucleus into the thermodynamically stable beta-phase.     

Can amorphous nuclei grow crystalline clathrates? The size and crystallinity of critical clathrate nuclei

Recent studies reveal that amorphous intermediates are involved in the formation of clathrate hydrates under conditions of high driving force, raising two questions: first, how could amorphous nuclei grow into crystalline clathrates and, second, whether amorphous nuclei are intermediates in the formation of clathrate crystals for temperatures close to equilibrium. In this work, we address these two questions through large-scale molecular simulations. We investigate the stability and growth of amorphous and crystalline clathrate nuclei and assess the thermodynamics and kinetic factors that affect the crystallization pathway of clathrates. Our calculations show that the dissociation temperature of amorphous clathrates is just 10% lower than for the crystals, facilitating the formation of metastable amorphous intermediates. We find that, at any temperatures, the critical crystalline nuclei are smaller than critical amorphous nuclei. The temperature dependence of the critical nucleus size is well described by the Gibbs-Thomson relation, from which we extract a liquid-crystal surface tension in excellent agreement with experiments. Our analysis suggests that at high driving force the amorphous nuclei may be kinetically favored over crystalline nuclei because of lower free energy barriers of formation. We investigated the role of the initial structure and size of the nucleus on the subsequent growth of the clathrates, and found that both amorphous and sI crystalline nuclei yield crystalline clathrates. Interestingly, growth of the metastable sII crystal polymorph is always favored over the most stable sI crystal, revealing kinetic control of the growth and indicating that a further step of ripening from sII to sI is needed to reach the stable crystal phase. The latter results are in agreement with the observed metastable formation of sII CO(2) and CH(4) clathrate hydrates and their slow conversion to sI under experimental conditions.     

Kinetics of supersaturated vapor nucleation on molecular condensation nuclei

The kinetics of nucleation is calculated for a supersaturated vapor containing molecular condensation nuclei, that is, foreign molecules able to induce the formation of viable nuclei of a condensed phase by themselves. In contrast to the previous calculation, the possibility of the escape of molecular condensation nuclei from very small clusters containing a few condensed vapor molecules is taken into account. More exact equations are derived for the rate of steady-state nucleation and the concentration of aerosol particles in a quasisteady-state regime of nucleation. The calculation demonstrates that, at a high probability of the escape of a molecular condensation nucleus, the predominating mechanism of cluster formation is the attachment of a molecular condensation nucleus to a cluster formed from vapor molecules rather than their condensation on the nucleus. At the same time, allowances for the possible escape of molecular condensation nuclei from clusters slightly affect the rate of nucleation and the concentration of aerosol particles being formed.