Molecular dynamics simulations of nucleation from vapor to solid composed of Lennard-Jones molecules

We performed molecular dynamics (MD) simulations of nucleation from vapor at temperatures below the triple point for systems consisting of 10(4)-10(5) Lennard-Jones (L-J) type molecules in order to test nucleation theories at relatively low temperatures. Simulations are performed for a wide range of initial supersaturation ratio (S(0) ? 10-10(8)) and temperature (kT = 0.2-0.6ε), where ε and k are the depth of the L-J potential and the Boltzmann constant, respectively. Clusters are nucleated as supercooled liquid droplets because of their small size. Crystallization of the supercooled liquid nuclei is observed after their growth slows. The classical nucleation theory (CNT) significantly underestimates the nucleation rates (or the number density of critical clusters) in the low-T region. The semi-phenomenological (SP) model, which corrects the CNT prediction of the formation energy of clusters using the second virial coefficient of a vapor, reproduces the nucleation rate and the cluster size distributions with good accuracy in the low-T region, as well as in the higher-T cases considered in our previous study. The sticking probability of vapor molecules onto the clusters is also obtained in the present MD simulations. Using the obtained values of sticking probability in the SP model, we can further refine the accuracy of the SP model.     

Molecular Dynamics Studies of the Kinetics of Phase Changes in Clusters II: Crystal Nucleation from Molten (RbCl)256 and (RbCl)500 Clusters

Molecular dynamics computer simulations have been carried out to study the effects of cluster size and temperature on the nucleation rate of rubidium chloride clusters in the temperature range of 500-650 K. Clusters with 256 and 500 RbCl molecules have been studied and the results are compared with those obtained from 108 molecule clusters. The melting point (MP) of the clusters was observed to increase with the size of the clusters and can be described by a linear equation MP=997-405 N^−1/3, where N is the number of molecules in the cluster. The nucleation rate is found to decrease with increasing cluster size or increasing nucleation temperature. Both classical nucleation theory and diffuse interface theory are used to interpret our observed results.     

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.     

Homogeneous nucleation rates of higher n-alcohols measured in a laminar flow diffusion chamber

Nucleation rate isotherms of n-butanol, n-pentanol, n-hexanol, n-heptanol, and n-octanol were measured in a laminar flow diffusion chamber using helium as carrier gas. The measurements were made at 250-310 K, corresponding to reduced temperatures of 0.43-0.50, and at atmospheric pressure. Experimental nucleation rate range was from 10(3) to 10(7) cm(-3) s(-1). The expression and accuracy of thermodynamic parameters, in particular equilibrium vapor pressure, were found to have a significant effect on calculated nucleation rates. The results were compared to the classical nucleation theory (CNT), the self-consistency corrected classical theory (SCC) and the Hale's scaled model of the CNT. The average ratio between the experimental and theoretical nucleation rates for all alcohols used was 1.5x10(3) when the CNT was used, and 0.2x10(-1) when the SCC was used and 0.7x10(-1) when the Hale's scaled theory was used. The average values represent all the alcohols used at the same reduced temperatures. The average ratio was about the same throughout the temperature range, although J(exp)/J(the) calculated with the Hale's scaled theory increased slightly with increasing temperature. The saturation ratio dependency was predicted closest to experiment with the classical nucleation theory. The nucleation rates were compared to those found in the literature. The measurements were in reasonable agreement with each other. The molecular content of critical alcohol clusters was between 35 and 80 molecules. At a fixed reduced temperature, the number of molecules in a critical cluster decreased as a function of alcohol carbon chain length. The number of molecules in critical clusters was compared to those predicted by the Kelvin equation. The theory predicted the critical cluster sizes well.     

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.     

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq)

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.     

Mean-field kinetic nucleation theory

A new semiphenomenological model of homogeneous vapor-liquid nucleation is proposed in which the cluster kinetics follows the "kinetic approach to nucleation" and the thermodynamic part is based on the revised Fisher droplet model with the mean-field argument for the cluster configuration integral. The theory is nonperturbative in a cluster size and as such is valid for all clusters down to monomers. It contains two surface tensions: macroscopic (planar) and microscopic. The latter is a temperature dependent quantity related to the vapor compressibility factor at saturation. For Lennard-Jones fluids the microscopic surface tension possesses a universal behavior with the parameters found from the mean-field density functional calculations. The theory is verified against nucleation experiments for argon, nitrogen, water, and mercury, demonstrating very good agreement with experimental data. Classical nucleation theory fails to predict experimental results when a critical cluster becomes small.     

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.     

氟化钾团簇的相变、成核和重结晶的分子动力学模拟

我们利用Born-Mayer-Huggins相互作用势函数对(KF)_N(N=108，256，500和864)团簇进行了分子动力学(MD)模拟。为了避免周期性边界条件对相变、成核和重结晶的干扰作用，对体系采用了自由边界。基于MD模拟结果，对团簇的熔化温度、熔化焓、扩散系数、成核速率、固液界面自由能、临界核大小等进行了计算和讨论。在对(KF)₈₆₄双晶团簇的热退火模拟中，观察到了固态的重结晶和晶粒的生长。经典的成核理论成功地解释了(KF)₈₆₄双晶团簇的重结晶MD模拟结果。     

Influence of thermostats and carrier gas on simulations of nucleation

We investigate the influence of carrier gas and thermostat on molecular dynamics (MD) simulations of nucleation. The task of keeping the temperature constant in MD simulations is not trivial and an inefficient thermalization may have a strong influence on the results. Different thermostating mechanisms have been proposed and used in the past. In particular, we analyze the efficiency of velocity rescaling, Nose-Hoover, and a carrier gas (mimicking the experimental situation) by extensive MD simulations. Since nucleation is highly sensitive to temperature, one would expect that small variations in temperature might lead to differences in nucleation rates of up to several orders of magnitude. Surprisingly, the results indicate that the choice of the thermostating method in a simulation does not have--at least in the case of Lennard-Jones argon--a very significant influence on the nucleation rate. These findings are interpreted in the context of the classical theory of Feder et al. [Adv. Phys. 15, 111 (1966)] by analyzing the temperature distribution of the nucleating clusters. We find that the distribution of cluster temperatures is non-Gaussian and that subcritically sized clusters are colder while postcritically sized clusters are warmer than the bath temperature. However, the average temperature of all clusters is found to be always higher than the bath temperature.     

Density functional theory of inhomogeneous liquids. IV. Squared-gradient approximation and classical nucleation theory

The squared-gradient approximation to the modified-core Van der Waals density functional theory model is developed. A simple, explicit expression for the SGA coefficient involving only the bulk equation of state and the interaction potential is given. The model is solved for planar interfaces and spherical clusters and is shown to be quantitatively accurate in comparison to computer simulations. An approximate technique for solving the SGA based on piecewise-linear density profiles is introduced and is shown to give reasonable zeroth-order approximations to the numerical solution of the model. The piecewise-linear models of spherical clusters are shown to be a natural extension of classical nucleation theory and serve to clarify some of the nonclassical effects previously observed in liquid-vapor nucleation. Nucleation pathways are investigated using both constrained energy-minimization and steepest-descent techniques.     

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.     

Heterogeneous multicomponent nucleation theorems for the analysis of nanoclusters

In this paper we present a new form of the nucleation theorems applicable to heterogeneous nucleation. These heterogeneous nucleation theorems allow, for the first time, direct determination of properties of nanoclusters formed on pre-existing particles from measured heterogeneous nucleation probabilities. The theorems can be used to analyze the size (first theorem) and the energetics (second theorem) of heterogeneous clusters independent of any specific nucleation model. We apply the first theorem to the study of small water and n-propanol clusters formed at the surface of 8 nm silver particles. According to the experiments the size of the two-component critical clusters is found to be below 90 molecules, and only less than 20 molecules for pure water, less than 300 molecules for pure n-propanol. These values are drastically smaller than the ones predicted by the classical nucleation theory, which clearly indicates that the nucleating clusters are too small to be quantitatively described using a macroscopic theory.     

Microscopic simulations of molecular cluster decay: Does the carrier gas affect evaporation?

We develop a kinetic theory of cluster decay by considering the stochastic motion of molecules within an effective potential of mean force (PMF) due to the cluster. We perform molecular dynamics simulations on a 50-atom argon cluster to determine the mean radial force on a component atom and hence the confining potential of mean force. Comparisons between isolated clusters and clusters thermostatted through the presence of a 100-atom helium carrier gas show that the heat bath has only a slight effect upon the PMF. This confirms the validity of calculations of cluster properties using isolated cluster simulations. The PMF is used to calculate the atomic evaporation rate from these clusters, and results are compared with the predictions of the capillarity approximation together with detailed balance, both components of the classical theory of aerosol nucleation.     

Implementation of dynamical nucleation theory effective fragment potentials method for modeling aerosol chemistry

In this work, the dynamical nucleation theory (DNT) model using the ab initio based effective fragment potential (EFP) is implemented for evaluating the evaporation rate constant and molecular properties of molecular clusters. Predicting the nucleation rates of aerosol particles in different chemical environments is a key step toward understanding the dynamics of complex aerosol chemistry. Therefore, molecular scale models of nanoclusters are required to understand the macroscopic nucleation process. On the basis of variational transition state theory, DNT provides an efficient approach to predict nucleation kinetics. While most DNT Monte Carlo simulations use analytic potentials to model critical sized clusters, or use ab initio potentials to model very small clusters, the DNTEFP Monte Carlo method presented here can treat up to critical sized clusters using the effective fragment potential (EFP), a rigorous nonempirical intermolecular model potential based on ab initio electronic structure theory calculations, improvable in a systematic manner. The DNTEFP method is applied to study the evaporation rates, energetics, and structure factors of multicomponent clusters containing water and isoprene. The most probable topology of the transition state characterizing the evaporation of one water molecule from a water hexamer at 243 K is predicted to be a conformer that contains six hydrogen bonds, with a topology that corresponds to two water molecules stacked on top of a quadrangular (H(2)O)(4) cluster. For the water hexamer in the presence of isoprene, an increase in the cluster size and a 3-fold increase in the evaporation rate are predicted relative to the reaction in which one water molecule evaporates from a water hexamer cluster.     

Argon nucleation: bringing together theory, simulations, and experiment

We present an overview of the current status of experimental, theoretical, molecular dynamics (MD), and density functional theory (DFT) studies of argon vapor-to-liquid nucleation. Since the experimental temperature-supersaturation domain does not overlap with the corresponding MD and DFT domains, separate comparisons have been made: theory versus experiment and theory versus MD and DFT. Three general theoretical models are discussed: Classical nucleation theory (CNT), mean-field kinetic nucleation theory (MKNT), and extended modified liquid drop model-dynamical nucleation theory (EMLD-DNT). The comparisons are carried out for the area below the MKNT pseudospinodal line. The agreement for the nucleation rate between the nonclassical models and the MD simulations is very good--within 1-2 orders of magnitude--while the CNT deviates from simulations by about 3-5 orders of magnitude. Perfect agreement is demonstrated between DFT results and predictions of MKNT (within one order of magnitude), whereas CNT and EMLD-DNT show approximately the same deviation of about 3-5 orders of magnitude. At the same time the agreement between all theoretical models and experiment remains poor--4-8 orders of magnitude for MKNT, 12-14 orders for EMLD-DNT, and up to 26 orders for CNT. We discuss possible reasons for this discrepancy and the ways to carry out experiment and simulations within the common temperature-supersaturation domain in order to produce a unified picture of argon nucleation.     

Kinetic theory of binary nucleation based on a first passage time analysis

The binary classical nucleation theory (BCNT) is based on the Gibbsian thermodynamics and applies the macroscopic concept of surface tension to nanosize clusters. This leads to severe inconsistencies and large discrepancies between theoretical predictions and experimental results regarding the nucleation rate. We present an alternative approach to the kinetics of binary nucleation which avoids the use of classical thermodynamics for clusters. The new approach is an extension to binary mixtures of the kinetic theory previously developed by Narsimhan and Ruckenstein and Ruckenstein and Nowakowski [J. Colloid Interface Sci. 128, 549 (1989); 137, 583 (1990)] for unary nucleation which is based on molecular interactions and in which the rate of emission of molecules from a cluster is determined via a mean first passage time analysis. This time is calculated by solving the single-molecule master equation for the probability distribution of a "surface" molecule moving in a potential field created by the cluster. The starting master equation is a Fokker-Planck equation for the probability distribution of a surface molecule with respect to its phase coordinates. Owing to the hierarchy of characteristic time scales in the evolution of the molecule, this equation can be reduced to the Smoluchowski equation for the distribution function involving only the spatial coordinates. The new theory is combined with density functional theory methods to determine the density profiles. This is essential for nucleation in binary systems particularly when one of the components is surface active. Knowing these profiles, one can determine the potential fields created by the cluster, its rate of emission of molecules, and the nucleation rate more accurately than by using the uniform density approximation. The new theory is illustrated by numerical calculations for a model binary mixture of Lennard-Jones monomers and rigidly bonded dimers of Lennard-Jones atoms. The amphiphilic character of the dimer component (i.e., its surface activity) is induced by the asymmetry in the interaction between a monomer and the two different sites of a dimer. The inconsistencies of the BCNT are avoided in the new theory.     

Monte Carlo simulations of critical cluster sizes and nucleation rates of water

We have calculated the critical cluster sizes and homogeneous nucleation rates of water at temperatures and vapor densities corresponding to experiments by Wolk and Strey [J. Phys. Chem B 105, 11683 (2001)]. The calculations have been done with an expanded version of a Monte Carlo method originally developed by Vehkamaki and Ford [J. Chem. Phys. 112, 4193 (2000)]. Their method calculates the statistical growth and decay probabilities of molecular clusters. We have derived a connection between these probabilities and kinetic condensation and evaporation rates, and introduce a new way for the calculation of the work of formation of clusters. Three different interaction potential models of water have been used in the simulations. These include the unpolarizable SPC/E [J. Phys. Chem. 91, 6269 (1987)] and TIP4P [J. Chem. Phys. 79, 926 (1983)] models and a polarizable model by Guillot and Guissani [J. Chem. Phys. 114, 6720 (2001)]. We show that TIP4P produces critical cluster sizes and a temperature and vapor density dependence for the nucleation rate that agree well with the experimental data, although the magnitude of nucleation rate is constantly overestimated by a factor of 2 x 10(4). Guissani and Guillot's model is somewhat less successful, but both the TIP4P and Guillot and Guissani models are able to reproduce a much better experimental temperature dependency of the nucleation rate than the classical nucleation theory. Using SPC/E results in dramatically too small critical clusters and high nucleation rates. The water models give different average binding energies for clusters. We show that stronger binding between cluster molecules suppresses the decay probability of a cluster, while the growth probability is not affected. This explains the differences in results from different water models.     

Impurity-stimulated heterogeneous nucleation of supercooled H2 clusters

The sizes and mass spectra of large (N=1900-13,700 molecules) cold (approximately 3.1 K) H2 clusters have been measured after scattering from CO molecules. Cluster-size measurements after between 2 and 8 collisions indicate that 7% of the H2 molecules are evaporated. This loss agrees with calculations for the number of H2 molecules evaporated by the heat released in the transition from an initial liquid state to a final solid state. Even though heterogeneous nucleation is initiated after only a few collisions with CO molecules, the mass spectra show that additional captured CO molecules coagulate to form large CO clusters with up to n=11 molecules, suggesting that the outer layer is sufficiently liquidlike to facilitate mobility of the CO molecules. Since the calculated H2 cluster temperature (approximately 3.1 K) is below the superfluid transition temperature predicted for pH2 with density between 40% and 80% of the triple-point density, a shell-like region of low density near the cluster surface can be expected to be superfluid.     

Molecular dynamics study of homogeneous nucleation in supercooled clusters of sodium bromide

Molecular dynamics simulations are carried out on clusters comprised of 108, 256, 500, and 864 Na⁺Br⁻ ion pairs represented by a Born–Mayer–Hugins interaction potential. Clusters with free boundaries are chosen in order to avoid the interference with nucleation caused by periodic boundary conditions. The melting point increases with increasing size of cluster as expected. The nucleation has been found to start near the surface of a NaBr cluster. The rates of nucleation in melted clusters are estimated based on classical nucleation theory. The interfacial free energies of 73.3–76.4 mJ/m² in temperature range of 400–550 K derived from kinetics of freezing are in the same order of magnitude as those predicted by Turnbull's relation and Buckle and Ubbelohde's observation. Sizes of critical nuclei are inferred from classical expressions and Voronoi polyhedra analyses.