Connection between the virial equation of state and physical clusters in a low density vapor

We carry out Monte Carlo simulations of physical Lennard-Jones and water clusters and show that the number of physical clusters in vapor is directly related to the virial equation of state. This relation holds at temperatures clearly below the critical temperatures, in other words, as long as the cluster-cluster interactions can be neglected--a typical assumption used in theories of nucleation. Above a certain threshold cluster size depending on temperature and interaction potential, the change in cluster work of formation can be calculated analytically with the recently proposed scaling law. The breakdown of the scaling law below the threshold sizes is accurately modeled with the low order virial coefficients. Our results indicate that high order virial coefficients can be analytically calculated from the lower order coefficients when the scaling law for cluster work of formation is valid. The scaling law also allows the calculation of the surface tension and equilibrium vapor density with computationally efficient simulations of physical clusters. Our calculated values are in good agreement with those obtained with other methods. We also present our results for the curvature dependent surface tension of water clusters.     

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.     

Surface tension and scaling of critical nuclei in diatomic and triatomic fluids

Density functional theory has been used to investigate surface tension and scaling of critical clusters in fluids consisting of diatomic and rigid triatomic molecules. The atomic sites are hard spheres with attractive interactions obtained from the tail part of the Lennard-Jones potential. Asymmetry in attractive interactions between the atomic sites has been introduced to cause molecular orientation and oscillatory density profiles at liquid-vapor interfaces. The radial dependence of cluster surface tension in fluids showing modest orientation in unimolecular layer at the interface or no orientation at all resembles the surface tension behavior of clusters in simple monoatomic fluids, although the surface tension maximum becomes more pronounced with increasing chain length of the molecule. Surface tension of clusters having multiple oscillatory layers at the interface shows a prominent maximum at small cluster sizes; however, the surface tension of large clusters is lower than the planar value. The scaling relation for the number of molecules in the critical cluster and the nucleation barrier height developed by McGraw and Laaksonen [Phys. Rev. Lett. 76, 2754 (1996)] are well obeyed for fluids with little structure at liquid-vapor interface. However, fluids having enhanced interfacial structure show some deviation from the particle number scaling, and the barrier height scaling breaks up seriously.     

Computer simulations of nucleation of nanoparticle superclusters from solution

This paper presents simulation studies of nanoparticle supercluster (NPSC) nucleation from a temperature quenched system. The nanoparticles are represented as 5 nm, spherical gold nanoparticles ligated with alkane thiols. The pair potential accounts for the van der Waals interaction between the metallic cores and ligand-ligand and ligand-solvent interactions. Phenomena well-known for molecular systems are observed including a prenucleation induction period, fluctuating prenucleation clusters that predominately add monomers one at a time, a critical nucleus size, and growth of NPSCs from solution in the presence of an equilibrium supernatant, all consistent with classical nucleation theory. However, only the largest prenucleating clusters are dense, and the cluster size can occasionally range greater than the critical size in the prenucleation regime until a cluster with low enough energy occurs, then nucleation ensues. Late in the nucleation process, the clusters display a crystalline structure that is a random mix of face-centered cubic (fcc) and hexagonal close-packed (hcp) lattices and indistinguishable from a randomized icosahedra structure.     

Rethinking the application of the first nucleation theorem to particle formation

The critical cluster is the threshold size above which a cluster will be more likely to grow than to evaporate. In field and laboratory measurements of new particle formation, the number of molecules of a given species in the critical cluster is commonly taken to be the slope of the log-log plot of the formation rate versus the concentration of the species. This analysis is based on an approximate form of the first nucleation theorem, which is derived with the assumption that there are no minima in the free energy surface prior to the maximum at the critical size. However, many atmospherically relevant systems are likely to exhibit such minima, for example, ions surrounded by condensable vapour molecules or certain combinations of acids and bases. We have solved numerically the birth-death equations for both an electrically neutral one-component model system with a local minimum at pre-critical sizes and an ion-induced case. For the ion-induced case, it is verified that the log-log slope of the nucleation rate versus particle concentration plot gives accurately the difference between the cluster sizes at the free energy maximum and minimum, as is expected from the classical form of the ion-induced nucleation rate. However, the results show that applying the nucleation theorem to neutral systems with stable pre-nucleation clusters may lead to erroneous interpretations about the nature of the critical cluster.     

Gradient theory computation of the radius-dependent surface tension and nucleation rate for n-nonane clusters

The Van der Waals-Cahn-Hilliard gradient theory (GT) is applied to determine the structure and the work of formation of clusters in supersaturated n-nonane vapor. The results are analyzed as functions of the difference of pressures of the liquid phase and vapor phase in chemical equilibrium, which is a measure for the supersaturation. The surface tension as a function of pressure difference shows first a weak maximum and then decreases monotonically. The computed Tolman length is in agreement with earlier results [L. Granasy, J. Chem. Phys. 109, 9660 (1998)] obtained with a different equation of state. A method based on the Gibbs adsorption equation is developed to check the consistency of GT results (or other simulation techniques providing the work of formation and excess number of molecules), and to enable an efficient interpolation. A cluster model is devised based on the density profile of the planar phase interface. Using this model we analyze the dependency of the surface tension on the pressure difference. We find three major contributions: (i) the effect of asymmetry of the density profile resulting into a linear increase of the surface tension, (ii) the effect of finite thickness of the phase interface resulting into a negative quadratic term, and (iii) the effect of buildup of a low-density tail of the density profile, also contributing as a negative quadratic term. Contributions (i)-(iii) fully explain the dependency of the surface tension on the pressure difference, including the range relevant to nucleation experiments. Contributions (i) and (ii) can be predicted from the planar density profile. The work of formation of noncritical clusters is derived and the nucleation rate is computed. The computed nucleation rates are closer to the experimental nucleation rate results than the classical Becker-D?ring theory, and also the dependence on supersaturation is better predicted.     

Heterogeneous condensation of the Lennard-Jones vapor onto a nanoscale seed particle

The heterogeneous condensation of a Lennard-Jones vapor onto a nanoscale seed particle is studied using molecular dynamics simulations. Measuring the nucleation rate and the height of the free energy barrier using the mean first passage time method shows that the presence of a weakly interacting seed has little effect on the work of forming very small cluster embryos but accelerates the rate by lowering the barrier for larger clusters. We suggest that this results from a competition between the energetic and entropic features of cluster formation in the bulk and at the heterogeneity. As the interaction is increased, the free energy of formation is reduced for all cluster sizes. We also develop a simple phenomenological model of film formation on a small seed that captures the general features of the nucleation process for small heterogeneities. A comparison of our simulation results with the model shows that heterogeneous classical nucleation theory provides a good estimate of the critical size of the film but significantly overestimates the size of the barrier.     

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.     

Molecular dynamics investigation of the transient regime in the freezing of salt clusters

Molecular dynamics (MD) investigations of the freezing of supercooled liquids can identify nuclei far smaller than can be detected in laboratory experiments, to date, and consequently can provide information about nucleation so far inaccessible to experiment. In a recent MD study of the freezing of clusters of SeF6, a new method of recording nucleation events was introduced. It involved the observation of times of first appearance of nuclei of the size of n. An advantage of the new approach is that it provides information about the size of the critical nucleus. For nuclei smaller than the critical size, it also avoids the overshoots of nucleation rates that precluded the application of the Miloshev-Wu method in the subcritical region. Kinetic information in the transient regime can be characterized by three parameters, the time lag, the reduced moment, and the steady state nucleation rate. To get some idea of how general the new approach is, a very different system was investigated, that of clusters of NaCl. Two different fitting functions were used to analyze the results. The first one adopted the log-normal probability distribution of Wu. The second function was a modification of Shneidman's analytical solution appropriate for large nuclei. The second function gave a rather good account of MD data for all nuclear sizes and temperatures and gave more stable results in the subcritical region. Several inferences of the sizes of critical nuclei were made. Applying the criterion for n* based on the Zeldovich solution of the Becker-D?ring equations, we estimated the critical nucleus sizes to be 14, 18, and 24 ions for quench temperatures of 640, 690, and 740 K, respectively. Even though the interionic interactions initiating nucleation in salt are very different from the van der Waals interactions in clusters of SeF6, the characteristic aspects of the transient regimes of the two systems were quite similar.     

Nucleation in cylindrical capillaries

We use a local density functional theory in the square gradient approximation to explore the properties of critical nuclei for the liquid-vapor transition of van der Waals fluids in cylindrical capillaries. The proposed model allows us to investigate the effect of pore size, surface field, and supersaturation on the behavior of the system. Our calculations predict the existence of at least three different pathways for the nucleation of droplets and bubbles in these confined fluids: axisymmetric annular bumps and lenses, and asymmetric droplets. The morphological transition between these different structures is driven by the existence of states of zero compressibility in the capillary. We show that the classical capillarity theory provides surprisingly accurate predictions for the work of formation of critical nuclei in cylindrical pores when line tension contributions to the free energy are taken into account.     

Surface Tension of the Vapor–Liquid Interface with Finite Curvature

Based on the division of particles into internal and surface particles, the expression is derived closing the system of equations of classical thermodynamics for curvature-dependent surface tension, equimolar radius, and radius of tension surface. A solution to this system allows one to find the surface tension of new phase nucleus of any size (including minimal) and any sign of surface curvature. The obtained results indicate the weak size dependence of thermodynamic parameters that are the functions of surface tension; it is shown that Tolman's length cannot be determined using experimental determination of these parameters. It is shown that the work of nucleus formation strongly depends on its size and is the function of effective rather than true surface tension. Numerical simulation of clusters by the molecular dynamics method indicates that the pressure inside a fairly small cluster is described by Laplace's formula with the coefficient of surface tension for the plane surface of a liquid that agrees with the proposed theory.     

Analysis of nucleation ability of cluster configurations with Monte Carlo simulations of argon

We determine the nucleation ability of argon clusters from Monte Carlo simulations. The nucleation rate appears to be defined by a sole characteristic of the clusters, namely, the stability. The stability is calculated as the ratio of grand canonical growth and decay rates and can be assigned to individual cluster configurations. We study the connection between the stability of the cluster configurations and their volume and total potential energy. Neither the potential energy nor the volume of a cluster configuration has a clear relation to its stability, and thus to the nucleation ability. On the other hand, we show that it is possible to use a specific volume for each cluster size to calculate the work of the cluster formation. These clusters with a unique volume have the same average stability as the full set of clusters. Our simulation method allows us to study the effect of possible deviations from equilibrium in the cluster configuration distributions. We argue that the nucleation process itself can produce a source for such a deviation. We show that even a small deviation from equilibrium in the cluster configuration distribution can lead to a dramatic deceleration of the nucleation rate. Although our simulations may overestimate the magnitude of the effect, they give qualitative estimates for its importance.     

Recent developments in the kinetic theory of nucleation

A review of recent progress in the kinetics of nucleation is presented. In the conventional approach to the kinetic theory of nucleation, it is necessary to know the free energy of formation of a new-phase particle as a function of its independent variables at least for near-critical particles. Thus the conventional kinetic theory of nucleation is based on the thermodynamics of the process. The thermodynamics of nucleation can be examined by using various approaches, such as the capillarity approximation, density functional theory, and molecular simulation, each of which has its own advantages and drawbacks. Relatively recently a new approach to the kinetics of nucleation was proposed [Ruckenstein E, Nowakowski B. J Colloid Interface Sci 1990;137:583; Nowakowski B, Ruckenstein E. J Chem Phys 1991;94:8487], which is based on molecular interactions and does not employ the traditional thermodynamics, thus avoiding such a controversial notion as the surface tension of tiny clusters involved in nucleation. In the new kinetic theory the rate of emission of molecules by a new-phase particle is determined with the help of a mean first passage time analysis. This time is calculated by solving the single-molecule master equation for the probability distribution function of a surface layer molecule moving in a potential field created by the rest of the cluster. The new theory was developed for both liquid-to-solid and vapor-to-liquid phase transitions. In the former case the single-molecule master equation is the Fokker-Planck equation in the phase space which can be reduced to the Smoluchowski equation owing to the hierarchy of characteristic time scales. In the latter case, the starting master equation is a Fokker-Planck equation for the probability distribution function of a surface layer molecule with respect to both its energy and phase coordinates. Unlike the case of liquid-to-solid nucleation, this Fokker-Planck equation cannot be reduced to the Smoluchowski equation, but the hierarchy of time scales does allow one to reduce it to the Fokker-Plank equation in the energy space. The new theory provides an equation for the critical radius of a new-phase particle which in the limit of large clusters (low supersaturations) yields the Kelvin equation and hence an expression for the macroscopic surface tension. The theory was illustrated with numerical calculations for a molecular pair interaction potential combining the dispersive attraction with the hard-sphere repulsion. The results for the liquid-to-solid nucleation clearly show that at given supersaturation the nucleation rate depends on the cluster structure (for three cluster structures considered-amorphous, fcc, and icosahedral). For both the liquid-to-solid and vapor-to-liquid nucleation, the predictions of the theory are consistent with the results of classical nucleation theory (CNT) in the limit of large critical clusters (low supersaturations). For small critical clusters the new theory provides higher nucleation rates than CNT. This can be accounted for by the fact that CNT uses the macroscopic interfacial tension which presumably overpredicts the surface tension of small clusters, and hence underpredicts nucleation rates.     

Homogeneous nucleation at high supersaturation and heterogeneous nucleation on microscopic wettable particles: A hybrid thermodynamic/density-functional theory

Homogeneous nucleation at high supersaturation of vapor and heterogeneous nucleation on microscopic wettable particles are studied on the basis of Lennard-Jones model system. A hybrid classical thermodynamics and density-functional theory (DFT) approach is undertaken to treat the nucleation problems. Local-density approximation and weighted-density approximation are employed within the framework of DFT. Special attention is given to the disjoining pressure of small liquid droplets, which is dependent on the thickness of wetting film and radius of the wettable particle. Different contributions to the disjoining pressure are examined using both analytical estimations and numerical DFT calculation. It is shown that van der Waals interaction results in negative contribution to the disjoining pressure. The presence of wettable particles results in positive contribution to the disjoining pressure, which plays the key role in the heterogeneous nucleation. Several definitions of the surface tension of liquid droplets are discussed. Curvature dependence of the surface tension of small liquid droplets is computed. The important characteristics of nucleation, including the formation free energy of the droplet and nucleation barrier height, are obtained.     

Thermodynamics of epitaxial calcite nucleation on self-assembled monolayers

Classical heterogeneous nucleation theory is used to describe the epitaxial nucleation of calcite on self-assembled monolayers (SAMs). Both spherical and faceted clusters are considered. The use of faceted clusters reveals a useful relation between the shape of very small crystals and the ratio of the heterogeneous and homogeneous nucleation barriers. The experimental approach of this paper concerns the measurement of the threshold driving forces for both homogeneous and heterogeneous nucleation of calcite. This is accomplished by preparing solutions with well-defined driving forces and by measuring the resulting types of nucleation that are observed after a fixed experimental time. The results of the experiments and the theoretical shape analysis are compared, and it is shown that in the experiments no homogeneous nucleation of calcite occurs for driving forces up to at least Deltamu/k(B)T approximately equal to 6.0. A calculation of the critical cluster size for heterogeneous nucleation results in a range of 2-28 growth units and faceted critical clusters from 3-28 growth units, depending on the value of the surface free energy of calcite. These sizes are 50-100 times smaller than the crystalline domain sizes of SAMs and therefore small enough to explain the promoting effect of the substrate.     

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.     

A kinetic model for nonisothermic homogeneous nucleation

A model for isothermal homogeneous nucleation is proposed that improves the classical model. A quasiequilibrium distribution of clusters was calculated on a basis of the Frenkel’-Lothe-Pound theory. The dependence of the free energy of clusters on their size was represented by an interpolation formula relating the free energy of dimers and large clusters to which a notion of macroscopic surface tension is applicable. The nucleation rate and the dependence of the cluster temperature on their size were calculated by balance equations describing the heating of from a cluster due to the condensation of monomers and its cooling due to collisions with an ambient gas. It is shown that the nucleation rate in excess buffer gas is higher than for the pure condensing gas by approximately two orders of magnitude. The model adequately describes the experimental data for the nucleation of methanol supersaturated vapor.