Molecular simulation of the homogeneous crystal nucleation of carbon dioxide

We report on a molecular simulation study of the homogeneous nucleation of CO2 in the supercooled liquid at low pressure (P = 5 MPa) and for degrees of supercooling ranging from 32% to 60%. In all cases, regardless of the degree of supercooling, the structure of the crystal nuclei is that of the Pa3 phase, the thermodynamically stable phase. For the more moderate degree of supercooling of 32%, the nucleation is an activated process and requires a method to sample states of high free energy. In this work, we apply a series of bias potentials, which promote the ordering of the centers of mass of the molecules and allow us to gradually grow crystal nuclei. The reliability of the results so obtained is assessed by studying the evolution of the nuclei in the absence of any bias potential, and by determining their probability of growth. We estimate that the size of the critical nucleus, for which the probability of growth is 0.5, is approximately 240 molecules. Throughout the nucleation process, the crystal nuclei clearly exhibit a Pa3 structure, in apparent contradiction with Ostwald's rule of stages. The other polymorphs have a much larger free energy. This makes their formation highly unlikely and accounts for the fact that the nucleation of CO2 proceeds directly in the stable Pa3 structure.     

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.     

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.     

Polymorph selection during the crystallization of Yukawa systems

Using molecular-dynamics simulations, we study the crystallization of supercooled liquids of charge-stabilized colloidal suspensions, modeled by the Yukawa (screened-Coulomb) potential. By modifying the value of the screening parameter lambda, we are able to invert the stability of the body-centered cubic (bcc) and face-centered cubic (fcc) polymorphs and study the crystal nucleation and growth in the domain of stability of each polymorph. We show that the crystallization mechanism strongly depends on the value of lambda. When bcc is the stable polymorph (lambda=3), the crystallization mechanism is straightforward. Both kinetics and thermodynamics favor the formation of the bcc particles and polymorph selection takes place early during the nucleation step. When fcc is the stable polymorph (lambda=10), the molecular mechanism is much more complex. First, kinetics favor the formation of bcc particles during the nucleation step. The growth of the post-critical nucleus proceeds through the successive cross-nucleation of the stable fcc polymorph on the metastable hcp polymorph as well as of the hcp polymorph on the fcc polymorph. As a result, polymorph selection occurs much later, i.e., during the growth step, than for lambda=3. We then extend our findings established in the case of homogeneous crystal nucleation to a situation of practical interest, i.e., when a seed of the stable polymorph is used. We demonstrate that the growth from the (111) face of a perfect fcc crystal into the melt proceeds through the same mechanisms.     

Physical basis for constrained lattice density functional theory

To study nucleation phenomena in an open system, a constrained lattice density functional theory (LDFT) method has been developed before to identify the unstable directions of grand potential functional and to stabilize nuclei by imposing a suitable constraint. In this work, we answer several questions about the method on a fundamental level, and give a firmer basis for the constrained LDFT method. First, we demonstrate that the nucleus structure and free energy barrier from a volume constraint method are equivalent to those from a surface constraint method. Then, we show that for the critical nucleus, the constrained LDFT method in fact produces a bias-free solution for both the nucleus structure and nucleation barrier. Finally, we give a physical interpretation of the Lagrange multiplier in the constraint method, which provides the generalized force to stabilize a nucleus in an open system. The Lagrange multiplier is found to consist of two parts: part I of the constraint produces an effective pressure, and part II imposes a constraint to counteract the supersaturation.     

Neutron scattering and monte carlo determination of the variation of the critical nucleus size with quench depth

We have used a combination of neutron scattering experiments and Monte Carlo simulations to study the initial stages of first-order phase transitions. We focus on quenches wherein the nascent phase is formed by homogeneous nucleation, and we approach the spinodal, i.e., the quench depth at which the original phase becomes unstable. In this regime, we show how critical nuclei sizes are determined from neutron scattering structure factors. Prevailing thought is that the size of the critical nucleus should increase with increasing quench depth and diverge at the spinodal. To the contrary, our experiments and simulations indicate that the critical nucleus size decreases monotonically as quench depth is increased and is finite at the spinodal.     

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.     

Role of liquid polymorphism during the crystallization of silicon

Using molecular simulation, we establish the pivotal role played by liquid polymorphs during the crystallization of silicon. When undercooled at a temperature 20% below the melting point, a silicon melt is under the form of the highly coordinated, high-density liquid (HDL) polymorph. We find that crystallization starts with the formation, within the HDL liquid, of a nanosized droplet of the least stable liquid polymorph, known as the almost tetracoordinated low-density liquid (LDL) polymorph. We then show that the crystalline embryo forms within the LDL droplet, close to the interface with the surrounding HDL liquid, thereby following a pathway associated with a much lower free energy barrier than the direct formation of the crystalline embryo from the HDL liquid would have required. This implies that, for substances exhibiting liquid polymorphs, theories, like the classical nucleation theory, and empirical rules, like Ostwald's rule, should be modified to account for the role of liquid polymorphs in the nucleation process.     

Calculation of solid-liquid interfacial free energy: a classical nucleation theory based approach

We present a simple approach to calculate the solid-liquid interfacial free energy. This new method is based on the classical nucleation theory. Using the molecular dynamics simulation, we employ spherical crystal nuclei embedded in the supercooled liquids to create an ideal model of a homogeneous nucleation. The interfacial free energy is extracted by fitting the relation between the critical nucleus size and the reciprocal of the critical undercooling temperature. The orientationally averaged interfacial free energy is found to be 0.302+/-0.002 (in standard LJ unit). The temperature dependence of the interfacial free energy is also obtained in this work. We find that the interfacial free energy increases slightly with increasing temperature. The positive temperature coefficient of the interfacial free energy is in qualitative agreement with Spaepen's analysis [Solid State Phys. 47, FS181 (1994)] and Turnbull's empirical estimation [J. Appl. Phys. 21, 1022 (1950)].     

Nucleation near the spinodal: limitations of mean field density functional theory

We investigate the diverging size of the critical nucleus near the spinodal using the gradient theory (GT) of van der Waals and Cahn and Hilliard and mean field density functional theory (MFDFT). As is well known, GT predicts that at the spinodal the free energy barrier to nucleation vanishes while the radius of the critical fluctuation diverges. We show numerically that the scaling behavior found by Cahn and Hilliard for these quantities holds quantitatively for both GT and MFDFT. We also show that the excess number of molecules Deltag satisfies Cahn-Hilliard scaling near the spinodal and is consistent with the nucleation theorem. From the latter result, it is clear that the divergence of Deltag is due to the divergence of the mean field isothermal compressibility of the fluid at the spinodal. Finally, we develop a Ginzburg criterion for the validity of the mean field scaling relations. For real fluids with short-range attractive interactions, the near-spinodal scaling behavior occurs in a fluctuation dominated regime for which the mean field theory is invalid. Based on the nucleation theorem and on Wang's treatment of fluctuations near the spinodal in polymer blends, we infer a finite size for the critical nucleus at the pseudospinodal identified by Wang.     

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.     

Crystal nucleation and the solid-liquid interfacial free energy

We present the results of molecular dynamics simulation of crystal nucleation in a supercooled Lennard-Jones liquid. Temperature and baric dependences of the nucleation rate, the Zeldovich factor, nucleus size diffusion coefficient, the radius, and the pressure in a critical crystal nucleus are defined in computer simulation. The data obtained have been used in the framework of classical nucleation theory to calculate the effective surface energy of crystal nuclei γ(e). It is shown that the value of γ(e) at T = const exceeds the value of the interfacial free energy at a flat crystal-liquid interface γ(∞) and γ(e) < γ(∞) at p = const.     

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.     

Free-energy landscape of nucleation with an intermediate metastable phase studied using capillarity approximation

Capillarity approximation is used to study the free-energy landscape of nucleation when an intermediate metastable phase exists. The critical nucleus that corresponds to the saddle point of the free-energy landscape as well as the whole free-energy landscape can be studied using this capillarity approximation, and various scenarios of nucleation and growth can be elucidated. In this study, we consider a model in which a stable solid phase nucleates within a metastable vapor phase when an intermediate metastable liquid phase exists. We predict that a composite critical nucleus that consists of a solid core and a liquid wetting layer as well as pure liquid and pure solid critical nuclei can exist depending not only on the supersaturation of the liquid phase relative to that of the vapor phase but also on the wetting behavior of the liquid surrounding the solid. The existence of liquid critical nucleus indicates that the phase transformation from metastable vapor to stable solid occurs via the intermediate metastable liquid phase, which is quite similar to the scenario of nucleation observed in proteins and colloidal systems. By studying the minimum-free-energy path on the free-energy landscape, we can study the evolution of the composition of solid and liquid within nuclei which is not limited to the critical nucleus.     

Free energy landscapes for homogeneous nucleation of ice for a monatomic water model

We simulate the homogeneous nucleation of ice from supercooled liquid water at 220 K in the isobaric-isothermal ensemble using the MW monatomic water potential. Monte Carlo simulations using umbrella sampling are performed in order to determine the nucleation free energy barrier. We find the Gibbs energy profile to be relatively consistent with that predicted by classical nucleation theory; the free energy barrier to nucleation was determined to be ~18 k(B)T and the critical nucleus comprised ~85 ice particles. Growth from the supercooled liquid gives clusters that are predominantly cubic, whilst starting with a pre-formed subcritical nucleus of cubic or hexagonal ice results in the growth of predominantly that phase of ice only.     

Quasihomogeneous nucleation of amyloid beta yields numerical bounds for the critical radius, the surface tension, and the free energy barrier for nucleus formation

Amyloid aggregates are believed to grow through a nucleation mediated pathway, but important aggregation parameters, such as the nucleation radius, the surface tension of the aggregate, and the free energy barrier toward aggregation, have remained difficult to measure. Homogeneous nucleation theory, if applicable, can directly relate these parameters to measurable quantities. We employ fluorescence correlation spectroscopy to measure the particle size distribution in an aggregating solution of Alzheimer's amyloid beta molecule (Abeta(1-40)) and analyze the data from a homogeneous nucleation theory perspective. We observe a reproducible saturation concentration and a critical dependence of various aspects of the aggregation process on this saturation concentration, which supports the applicability of the nucleation theory to Abeta aggregation. The measured size distributions show a valley between two peaks ranging from 5 to 50 nm, which defines a boundary for the value of the nucleation radius. By carefully controlling the conditions to inhibit heterogeneous nucleation, we can hold off nucleation in a 25 times supersaturated solution for at least up to 3 h at room temperature. This quasi-homogeneous kinetics implies that at room temperature, the surface energy of the Abeta/water interface is > or =4.8 mJ/m(2), the free energy barrier to nucleation (at 25 times supersaturation) is > or =1.93x10(-19) J, and the number of monomers in the nucleus is > or =29.     

Nucleation and size distribution of nucleus during induction period of polyethylene crystallization

The crystallization process from supercooled melt results in the formation of nanosize nuclei in the earlier stage (induction period) through subsequent attachment or detachment of repeating unit to nuclei. The size distribution of nucleus f(N(j),t) in the induction period of nucleation process from the melts has not been experimentally confirmed yet by direct observation. The reason is that the number density of nuclei nu is too small to be detected experimentally. In our previous work, we showed the direct evidence of nucleation experimentally by means of small angle x-ray scattering (SAXS) technique. Further we have succeeded to observe the nucleation and f(N(j),t) of polymer crystallization from the melts by SAXS using synchrotron radiation. We increased nu by adding a nucleating agent to a polymer (polyethylene). The time evolution of f(N(j),t) was observed for the first time.     

Spinodal for the solution-to-crystal phase transformation

The formation of crystalline nuclei from solution has been shown for many systems to occur in two steps: the formation of quasidroplets of a disordered intermediate, followed by the nucleation of ordered crystalline embryos within these droplets. The rate of each step depends on a respective free-energy barrier and on the growth rate of its near-critical clusters. We address experimentally the relative significance of the free-energy barriers and the kinetic factors for the nucleation of crystals from solution using a model protein system. We show that crystal nucleation is 8-10 orders of magnitude slower than the nucleation of dense liquid droplets, i.e., the second step is rate determining. We show that at supersaturations of three or four k(B)T units, crystal nuclei of five, four, or three molecules transform into single-molecule nuclei, i.e., the significant nucleation barrier vanishes below the thermal energy of the molecules. We show that the main factor, which determines the rate of crystal nucleation, is the slow growth of the near-critical ordered clusters within the quasidroplets of the disordered intermediate. Analogous to the spinodal in supersaturated fluids, we define a solution-to-crystal spinodal from the transition to single-molecule crystalline nuclei. We show that heterogeneous nucleation centers accelerate nucleation not only because of the wettinglike effects that lower the nucleation barrier, as envisioned by classical theory, but by helping the kinetics of growth of the ordered crystalline embryos.     

Metadynamics simulations of ice nucleation and growth

The metadynamics method for accelerating rate events in molecular simulations is applied to the problem of ice freezing. We demonstrate homogeneous nucleation and growth of ice at 180 K in the isothermal-isobaric ensemble without the presence of external fields or surfaces. This result represents the first report of continuous and dynamic ice nucleation in a system of freely evolving density. Simulations are conducted using a variety of periodic simulation domains. In all cases the cubic polymorph ice I(c) is grown. The influence of boundary effects on estimates of the nucleation free energy barrier are discussed in relation to differences between this and earlier work.     

Monte Carlo simulation study of droplet nucleation

A new rigorous Monte Carlo simulation approach is employed to study nucleation barriers for droplets in Lennard-Jones fluid. Using the gauge cell method we generate the excess isotherm of critical clusters in the size range from two to six molecular diameters. The ghost field method is employed to compute the cluster free energy and the nucleation barrier with desired precision of (1-2)kT. Based on quantitative results obtained by Monte Carlo simulations, we access the limits of applicability of the capillarity approximation of the classical nucleation theory and the Tolman equation. We show that the capillarity approximation corrected for vapor nonideality and liquid compressibility provides a reasonable assessment for the size of critical clusters in Lennard-Jones fluid; however, its accuracy is not sufficient to predict the nucleation barriers for making practical estimates of the rate of nucleation. The established dependence of the droplet surface tension on the droplet size cannot be approximated by the Tolman equation for small droplets of radius less than four molecular diameters. We confirm the conclusion of ten Wolde and Frenkel [J. Chem. Phys. 109, 9901 (1998)] that integration of the normal component of the Irving-Kirkwood pressure tensor severely underestimates the nucleation barriers for small clusters.