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.     

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.     

Test of classical nucleation theory on deeply supercooled high-pressure simulated silica

We test classical nucleation theory (CNT) in the case of simulations of deeply supercooled, high density liquid silica, as modeled by the van Beest-Kramer-van Santen potential. We find that at density rho=4.38 gcm(3), spontaneous nucleation of crystalline stishovite occurs in conventional molecular dynamics simulations at temperature T=3000 K, and we evaluate the nucleation rate J directly at this T via "brute force" sampling of nucleation events in numerous independent runs. We then use parallel, constrained Monte Carlo simulations to evaluate DeltaG(n), the free energy to form a crystalline embryo containing n silicon atoms, at T=3000, 3100, 3200, and 3300 K. By comparing the form of DeltaG(n) to CNT, we test the ability of CNT to reproduce the observed behavior as we approach the regime where spontaneous nucleation occurs on simulation time scales. We find that the prediction of CNT for the n dependence of DeltaG(n) fits reasonably well to the data at all T studied. Deltamu, the chemical potential difference between bulk liquid and stishovite, is evaluated as a fit parameter in our analysis of the form of DeltaG(n). Compared to directly determined values of Deltamu extracted from previous work, the fitted values agree only at T=3300 K; at lower T the fitted values increasingly overestimate Deltamu as T decreases. We find that n(*), the size of the critical nucleus, is approximately ten silicon atoms at T=3300 K. At 3000 K, n(*) decreases to approximately 3, and at such small sizes methodological challenges arise in the evaluation of DeltaG(n) when using standard techniques; indeed even the thermodynamic stability of the supercooled liquid comes into question under these conditions. We therefore present a modified approach that permits an estimation of DeltaG(n) at 3000 K. Finally, we directly evaluate at T=3000 K the kinetic prefactors in the CNT expression for J, and find physically reasonable values; e.g., the diffusion length that Si atoms must travel in order to move from the liquid to the crystal embryo is approximately 0.2 nm. We are thereby able to compare the results for J at 3000 K obtained both directly and based on CNT, and find that they agree within an order of magnitude. In sum, our work quantifies how certain predictions of CNT (e.g., for Deltamu) break down in this deeply supercooled limit, while others [the n dependence of DeltaG(n)] are not as adversely affected.     

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.     

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.     

Experimental study of homogeneous nucleation from the bismuth supersaturated vapor: evaluation of the surface tension of critical nucleus

The homogeneous nucleation of bismuth supersaturated vapor is studied in a laminar flow quartz tube nucleation chamber. The concentration, size, and morphology of outcoming aerosol particles are analyzed by a transmission electron microscope (TEM) and an automatic diffusion battery (ADB). The wall deposit morphology is studied by scanning electron microscopy. The rate of wall deposition is measured by the light absorption technique and direct weighting of the wall deposits. The confines of the nucleation region are determined in the "supersaturation cut-off" measurements inserting a metal grid into the nucleation zone and monitoring the outlet aerosol concentration response. Using the above experimental techniques, the nucleation rate, supersaturation, and nucleation temperature are measured. The surface tension of the critical nucleus and the radius of the surface of tension are determined from the measured nucleation parameters. To this aim an analytical formula for the nucleation rate is used, derived from author's previous papers based on the Gibbs formula for the work of formation of critical nucleus and the translation-rotation correction. A more accurate approach is also applied to determine the surface tension of critical drop from the experimentally measured bismuth mass flow, temperature profiles, ADB, and TEM data solving an inverse problem by numerical simulation. The simulation of the vapor to particles conversion is carried out in the framework of the explicit finite difference scheme accounting the nucleation, vapor to particles and vapor to wall deposition, and particle to wall deposition, coagulation. The nucleation rate is determined from simulations to be in the range of 10(9)-10(11) cm(-3) s(-1) for the supersaturation of Bi(2) dimers being 10(17)-10(7) and the nucleation temperature 330-570 K, respectively. The surface tension σ(S) of the bismuth critical nucleus is found to be in the range of 455-487 mN/m for the radius of the surface of tension from 0.36 to 0.48 nm. The function σ(S) changes weakly with the radius of critical nucleus. The value of σ(S) is from 14% to 24% higher than the surface tension of a flat surface.     

Molecular Dynamics Studies of the Kinetics of Phase Changes in Clusters IV： Crystal Nucleation from Molten （NaCl）256 and （NaCl）500 Clusters

Molecular dynamics computer simulation based on the Born-Mayer-Huggins potential function has been carried out to study the effects of duster size and temperature on the nucleation rate of sodium chloride dusters in the temperature range of 580 K to 630 K. Clusters with 256 and 500 NaCl molecules have been studied and the results have been compared with those obtained from 108 molecule dusters. The melting point (MP) of the clusters were observed to increase with the size of the clusters and can be well described by a linear equation MP =1107(37)-1229(23)N^-1/3(N is the number of molecules in the duster).The nucleation rate was found to decrease with increasing the duster size or temperature. Various nucleation theories have been used to interpret the nucleation rates obtained from this molecular dynamics simulation. It is possible to use a constant diffuse interface thickness to interpret the nucleation rate from the diffuse interface theory in the temperature range of this study. However, the interfacinl free energy estimated from classical nucleation theory and diffuse interface theory increases too fast with increasing the temperature while that from Gran-Gunton theory does not change with changing temperatures.The sizes of critical nuclei estimated from all the theories are smaller than those estimated from our simulations.     

Nucleation rate isotherms of argon from molecular dynamics simulations

We report six nucleation rate isotherms of vapor-liquid nucleation of Lennard-Jones argon from molecular dynamics simulations. The isotherms span three orders of magnitude in nucleation rates, 10(23)相似文献   

Water vapor nucleation on an infinite substrate with complementary structure: 1. The mechanism of nucleation

The nucleation of water vapor on the infinite surface of a silver iodide crystal at 260 K is simulated. Long-range electrostatic and polarization interactions are taken into account by the Ewald method. The free energy and work of equilibrium formation of nuclei are calculated at the molecular level by the method of bicanonical statistical ensemble. It is shown that, at the initial stage, the substrate is completely covered with a water monolayer. The substrate tends to decrease by two orders of magnitude the vapor pressure required to form the critical nucleus of a monomolecular film with a size of 10² molecules, the nucleation rate being increased by tens of orders of magnitude as compared to homogeneous nucleation. The saturation pressure above the adsorbed monomolecular film is 12 times lower than that above the flat ice surface. The free energy at the edges of “spots” per unit length is 1.4 × 10^?11 J/m. The critical size of the spot increases with a decrease in vapor pressure as the inverse second power of the logarithm of pressure.     

Enhanced protein crystallization around the metastable critical point

We report on a computer-simulation study of homogeneous crystal nucleation in a model for globular proteins. We find that the presence of a metastable vapour-liquid critical point drastically changes the pathway for the formation of a critical nucleus. But what is more important, the large density fluctuations near the critical point also lowers the free-energy barrier to nucleation and hence increases the nucleation rate. As␣the location of the vapour-liquid critical point can be controlled by changing the solvent conditions, our simulation results suggest a guided approach to protein crystallization. Received: 4 June 1998/Accepted: 3 September 1998 / Published online: 10 December 1998     

Nucleation of NaCl nanoparticles in supercritical water: molecular dynamics simulations

Formation of NaCl nanoparticles in supercritical water is studied using molecular dynamics simulation method. We have simulated particle nucleation and growth in NaCl-H2O fluids, with salt concentration of 5.1 wt %, in the temperature and density range of 673-1073 K and 0.17-0.34 g/cm(3), respectively. The cluster size distributions, the size of critical nuclei and cluster lifetimes are reported. The size distribution of emerging clusters shows a very strong dependence on the system's density, with larger clusters forming at lower densities. Clusters consisting of approximately 14-24 ions appear critical for the thermodynamic states examined. The local structures of critical clusters are found to be amorphous. The lifetime values for clusters containing more than 20 ions are in the range of 10-50 ps. We have calculated the NaCl nucleation rates, which appear to be on the order of 10(28) cm(-3) s(-1).     

Comparison between the classical theory predictions and molecular simulation results for heterogeneous nucleation of argon

We have performed Monte Carlo simulations of homogeneous and heterogeneous nucleations of Lennard-Jones argon clusters. The simulation results were interpreted using the major concept posing a difference between the homogeneous and heterogeneous classical nucleation theories-the contact parameter. Our results show that the multiplication concept of the classical heterogeneous nucleation theory describes the cluster-substrate interaction surprisingly well even for small molecular clusters. However, in the case of argon nucleating on a rigid monolayer of fcc(111) substrate at T=60 K, the argon-substrate atom interaction being approximately one-third as strong as the argon-argon interaction, the use of the classical theory concept results in an underestimation of the heterogeneous nucleation rate by two to three orders of magnitude even for large clusters. The main contribution to this discrepancy is induced by the failure of the classical theory of homogeneous nucleation to predict the energy involved in bringing one molecule from the vapor to the cluster for clusters containing less than approximately 15 molecules.     

Ion-induced nucleation in solution: promotion of solute nucleation in charged levitated droplets

We have investigated the nucleation and growth of sodium chloride in both single quiescent charged droplets and charged droplet populations that were levitated in an electrodynamic levitation trap (EDLT). In both cases, the magnitude of a droplet's net excess charge (ions(DNEC)) influenced NaCl nucleation and growth, albeit in different capacities. We have termed the phenomenon ion-induced nucleation in solution. For single quiescent levitated droplets, an increase in ions(DNEC) resulted in a significant promotion of NaCl nucleation, as determined by the number of crystals observed. For levitated droplet populations, a change in NaCl crystal habit, from regular cubic shapes to dome-shaped dendrites, was observed once a surface charge density threshold of -9 x 10(-4) e.nm(-2) was surpassed. Although promotion of NaCl nucleation was observed for droplet population experiments, this can be attributed in part to the increased rate of solvent evaporation observed for levitated droplet populations having a high net charge. Promotion of nucleation was also observed for two organic acids, 2,4,6-trihydroxyacetophenone monohydrate (THAP) and alpha-cyano-4-hydroxycinnamic acid (CHCA). These results are of direct relevance to processes that occur in both soft-ionization techniques for mass spectrometry and to a variety of industrial processes. To this end, we have demonstrated the use of ion-induced nucleation in solution to form ammonium nitrate particles from levitated droplets to be used in in vitro toxicology studies of ambient particle types.     

Homogeneous ice nucleation from supercooled water

Homogeneous ice nucleation from supercooled water was studied in the temperature range of 220-240 K through combining the forward flux sampling method (Allen et al., J. Chem. Phys., 2006, 124, 024102) with molecular dynamics simulations (FFS/MD), based on a recently developed coarse-grained water model (mW) (Molinero et al., J. Phys. Chem. B, 2009, 113, 4008). The calculated ice nucleation rates display a strong temperature dependence, ranging from 2.148 ± 0.635 × 10(25) m(-3) s(-1) at 220 K to 1.672 ± 0.970 × 10(-7) m(-3) s(-1) at 240 K. These rates can be fitted according to the classical nucleation theory, yielding an estimate of the effective ice-water interface energy γ(ls) of 31.01 ± 0.21 mJ m(-2) for the mW water model. Compared to experiments, our calculation underestimates the homogeneous ice nucleation rate by a few orders of magnitude. Possible reasons for the discrepancy are discussed. The nucleating ice embryo contains both cubic ice Ic and hexagonal ice Ih, with the fraction of each structure being roughly 50% when the critical size is reached. In particular, a novel defect structure containing nearly five-fold twin boundaries is identified in the ice clusters formed during nucleation. The way such defect structure is formed is found to be different from mechanisms proposed for the formation of the same defect in metallic nanoparticles and thin film. The quasi five-fold twin boundary structure found here is expected to occur in the crystallization of a wide range of materials with the diamond cubic structure, including ice.     

Stresses inside critical nuclei

The usual derivation of classical nucleation theory is inappropriate for crystal nucleation. In particular, it leads to a seriously flawed estimate of the pressure inside a critical nucleus. This has consequences for the prediction of possible metastable phases during the nucleation process. In this paper, we reanalyze the theory for crystal nucleation based on the thermodynamics of small crystals suspended in a liquid, due to Mullins (J. Chem. Phys. 1984, 81, 1436). As an illustration of the difference between the classical picture and the present approach, we consider a numerical study of crystal nucleation in binary mixtures of hard spherical colloids with a size ratio of 1:10. The stable crystal phase of this system can be either dense or expanded. We find that, in the vicinity of the solid-solid critical point where the crystallites are highly compressible, small crystal nuclei are less dense than large nuclei. This phenomenon cannot be accounted for by either classical nucleation theory or by the Gibbsian droplet model.     

Nucleation in electrochemical growth of the Ag(100) crystal face: determining the nucleus size via the nucleation theorem

We employ the nucleation theorem for a model-independent determination of the size of the two-dimensional (2D) Ag nucleus with the aid of experimental data for the nucleation-mediated electrochemical growth of the Ag(100) crystal face in aqueous solution of AgNO(3) at 318 K. These data are for the stationary rate of 2D nucleation, for the initial portion of the potentiostatic current transient pertaining to atomically smooth face, and for the galvanostatic current corresponding to stationary growth of the face. It turns out that the 2D nucleus is constituted of 17-64 Ag atoms when the overpotential is in the range of 12-22.4 mV. Upon expressing the overpotential in terms of supersaturation, it is found that the experimental data for the size of the 2D Ag nucleus are in conformity with existing simulation data for the size of the 2D nucleus on the (100) face of Kossel crystal (the simulation nucleus contains 1-30 atoms). It is found as well that the Gibbs-Thomson equation of the classical theory of 2D nucleation describes very well the supersaturation dependence of the size of both the Ag and the simulation nucleus.     

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)].     

Extended study of molecular dynamics simulation of homogeneous vapor-liquid nucleation of water

Using the simple point charge/extended water model, we performed molecular dynamics simulations of homogeneous vapor-liquid nucleation at various values of temperature T and supersaturation S, from which the nucleation rate J, critical nucleus size n(*), and the cluster formation free energy DeltaG were derived. As well as providing lots of simulation data, the results were compared with theories on homogeneous nucleation, including the classical, semi-phenomenological, and scaled models, but none of these gave a satisfactory explanation for our results. It was found that two main factors made the theories fail: (1) The average cluster structure including the nonspherical shape and the core structure that is not like the bulk liquid and (2) the forward rate which is larger than assumed by the theories by about one order of magnitude. The quantitative evaluation of these factors is left for future investigations.     

Formation of FeCl(2)/NaCl-nanoparticles in supercritical water investigated by molecular dynamics simulations: nucleation rates

Nucleation and growth of FeCl(2) in supercritical water containing NaCl at different state points between temperatures of 798 and 873 K and system densities between 0.24 and 0.14 g cm(-3) have been studied by molecular dynamics simulations. The number of NaCl ion pairs was chosen to simulate particle formation in seawater and brine of higher salinity. Rigid SPC/E water was used to model the water molecules while a combination of Coulomb and Lennard-Jones potentials was used for the ions. Two different methods for determination of nucleation rates are applied and their results compared. We find decreasing nucleation rates with both increasing temperature and decreasing system density. Our results are also compared to those we recently obtained in an investigation of pure FeCl(2) from supercritical water. We find both increasing nucleation rates and a decreasing size of the critical cluster with increasing amount of NaCl.     

Determining the nucleation rate from the dimer growth probability

A new method is proposed for the determination of the stationary one-component nucleation rate J with the help of data for the growth probability P2 of a dimer which is the smallest cluster of the nucleating phase. The method is based on an exact formula relating J and P2, and is readily applicable to computer simulations of nucleation. Using the method, the dependence of J on the supersaturation s is determined by kinetic Monte Carlo simulations of two-dimensional (2D) nucleation of monolayers on the (100) face of Kossel crystal. The change of J over nearly 11 orders of magnitude is followed and it is found that the classical nucleation theory overestimates the simulation J values by an s-dependent factor. The 2D nucleus size evaluated via the nucleation theorem is described satisfactorily by the classical Gibbs-Thomson equation and its corrected version accounting for the spinodal limit of 2D nucleation.