Kinetic theory of nucleation based on a first passage time analysis: improvement by the density-functional theory

A recent kinetic theory of nucleation [see, e.g., E. Ruckenstein and B. Nowakowski, J. Colloid Interface Sci. 137, 583 (1990)] is based on molecular interactions and avoids the traditional thermodynamics. The rate of emission of molecules from a cluster is found via a first passage time analysis. This time is calculated by solving the single-molecule master equation for the probability distribution function of a surface molecule located in the potential field created by the cluster. The liquid cluster was assumed to have sharp boundaries and uniform density. In the present paper, this assumption is removed by using the density-functional theory to find the density profiles. Thus, more accurate calculations of the potential field created by the cluster, its emission rate, and nucleation rate are obtained. The modified theory is illustrated by numerical calculations for a molecular pair interaction potential combining the dispersive attraction with the hard-sphere repulsion.     

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.     

A kinetic approach to the theory of heterogeneous nucleation on soluble particles during the deliquescence stage

Deliquescence is the dissolution of a solid nucleus in a liquid film formed on the nucleus due to vapor condensation. Previously, the kinetics of deliquescence was examined in the framework of the capillarity approximation which involves the thermodynamic interfacial tensions for a thin film and the approximation of uniform density therein. In the present paper we propose a kinetic approach to the theory of deliquescence which avoids the use of the above macroscopic quantities for thin films. The rates of emission of molecules from the liquid film into the vapor and from the solid core into the liquid film are determined through a first passage time analysis whereas the respective rates of absorption are calculated through the gas kinetic theory. The first passage time is obtained by solving the single-molecule master equation for the probability distribution of a "surface" molecule moving in a potential field created by the cluster. Furthermore, the time evolution of the liquid film around the solid core is described by means of two mass balance equations which involve the rates of absorption and emission of molecules by the film at its two interfaces. When the deliquescence of an ensemble of solid particles occurs by means of large fluctuations, the time evolution of the distribution of composite droplets (liquid film+solid core) with respect to the independent variables of state is governed by a Fokker-Planck kinetic equation. When both the vapor and the solid soluble particles are single component, this equation has the form of the kinetic equation of binary nucleation. A steady-state solution for this equation is obtained by the method of separation of variables. The theory is illustrated with numerical calculation regarding the deliquescence of spherical particles in a water vapor with intermolecular interactions of the Lennard-Jones kind. The new approach allows one to qualitatively explain an important feature of experimental data on deliquescence, namely the occurrence of nonsharp deliquescence, a feature that the previous deliquescence theory based on classical thermodynamics could not account for.     

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.     

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.     

Molecular dynamics of homogeneous nucleation in the vapor phase of Lennard-Jones. III. Effect of carrier gas pressure

A molecular dynamics simulation of vapor phase nucleation has been performed with 40,000 Lennard-Jones particles for the target gas and 0-160,000 particles for the carrier gas. Three carrier gas models are adopted, including a soft-core model, a Lennard-Jones model, and a modified Lennard-Jones model in which the attractive interaction can be adjusted. The effect of the carrier-gas pressure is assessed through computing and comparing the rate of nucleation and cluster size distribution. It is found that the effect of the carrier-gas pressure can be strongly dependent on the carrier-gas model. A positive effect (enhancement of the nucleation rate) is found with the soft-core potential model, whereas negligible effect is found with the Lennard-Jones potential model. For the modified Lennard-Jones potential with a weak attractive interaction, the carrier-gas effect is positive. However, the effect is negligible with a stronger attractive interaction between the target and carrier-gas particles. A reason for the negligible effect is that the carrier-gas particles are adsorbed on the cluster surface when the density of target and carrier-gas particles are comparable. When the density of carrier-gas particles are four times that of the target particles, the carrier-gas particles tend to mix with the target particles in the clusters.     

Kinetic theory of heterogeneous nucleation; effect of nonuniform density in the nuclei

The heterogeneous nucleation of a liquid from a vapor in contact with a planar solid surface or a solid surface with cavities is examined on the basis of the kinetic theory of nucleation developed by Nowakowski and Ruckenstein [J. Phys. Chem. 96 (1992) 2313] which is extended to nonuniform fluid density distribution (FDD) in the nucleus. The latter is determined under the assumption that at each moment the FDD in the nucleus is provided by the density functional theory (DFT) for a nanodrop. As a result of this assumption, the theory does not require to consider that the contact angle which the nucleus makes with the solid surface and the density of the nucleus are independent parameters since they are provided by the DFT. For all considered cases, the nucleation rate is higher in the cavities than on a planar surface and increases with increasing strength of the fluid-solid interactions and decreasing cavity radius. The difference is small at high supersaturations (small critical nuclei), but becomes larger at low supersaturations when the critical nucleus has a size comparable with the size of the cavity. The nonuniformity of the FDD in the nucleus decreases the nucleation rate when compared to the uniform FDD.     

Photochemically-induced crystallization of protein

I demonstrate photochemically induced crystallization of metastable hen egg-white lysozyme solution by weak UV irradiation for several tens seconds. The most effective irradiation time range is 10–60 s, and in this range the enzyme activity is maintained. Intermediates, neutral radicals at tryptophan residual produced by one-photon absorption, enhance nucleation. When the intermediate is selectively excited by visible light, the intermediate is denatured. At that time the light-induced nucleation is inhibited. This result indicates the intermediate induces nucleation. The radical forms lysozyme dimer that is detected by an SDS-PAGE electrophoresis experiment. An addition of polyethylene glycol (PEG) greatly enhances light-induced nucleation. PEG affects to shorten the intermediate radical lifetime, which suggests that PEG assists to form dimer. We consider that the photochemical dimer behaves as smallest cluster to grow critical nucleus. The smallest cluster formation is the rate determining step in classical nucleation theory due to surface energy disadvantage. The photochemical dimer is formed by a covalent bond, and the nucleation is initiated from stable dimer. The nucleation enhancement is reasonably explained. The present researches results point out the development of a new method for controlling nucleation and growth that could be applied for structural genomics and pharmaceutical industry for instance.     

Steady-state nucleation rate and flux of composite nucleus at saddle point

The steady-state nucleation rate and flux of composite nucleus at the saddle point is studied by extending the theory of binary nucleation. The Fokker-Planck equation that describes the nucleation flux is derived using the Master equation for the growth of the composite nucleus, which consists of the core of the final stable phase surrounded by a wetting layer of the intermediate metastable phase nucleated from a metastable parent phase recently evaluated by Iwamatsu [J. Chem. Phys. 134, 164508 (2011)]. The Fokker-Planck equation is similar to that used in the theory of binary nucleation, but the non-diagonal elements exist in the reaction rate matrix. First, the general solution for the steady-state nucleation rate and the direction of nucleation flux is derived. Next, this information is then used to study the nucleation of composite nucleus at the saddle point. The dependence of steady-state nucleation rate as well as the direction of nucleation flux on the reaction rate in addition to the free-energy surface is studied using a model free-energy surface. The direction of nucleation current deviates from the steepest-descent direction of the free-energy surface. The results show the importance of two reaction rate constants: one from the metastable environment to the intermediate metastable phase and the other from the metastable intermediate phase to the stable new phase. On the other hand, the gradient of the potential Φ or the Kramers crossover function (the commitment or splitting probability) is relatively insensitive to reaction rates or free-energy surface.     

A thermodynamically consistent kinetic framework for binary nucleation

The traditional theory for binary homogeneous nucleation follows the classical derivation of the nucleation rate in the supposition of a hypothetical constrained-equilibrium distribution in the calculation of the cluster evaporation rate. This model enables calculation of the nucleation rate, but requires evaluation of the cluster distribution and cluster properties for an unstable equilibrium with supersaturated vapor. An alternate derivation of the classical homomolecular nucleation rate eliminated the need for this nonphysical approximation by calculating the evaporative flux at full thermodynamic equilibrium. The present paper develops that approach for binary nucleation; the framework is readily extended to ternary nucleation. In this analysis, the evaporative flux is evaluated by applying mass balance at full thermodynamic equilibrium of the system under study. This approach eliminates both the need for evaluating cluster properties in an unstable constrained-equilibrium state and ambiguity in the normalization constant required in the nucleation-rate expression. Moreover, it naturally spans the entire composition range between the two pure monomers. The cluster fluxes derived using this new model are similar in form to those of classical derivations, so previously developed methods for evaluation of the net nucleation rate can be applied directly to the new formulation.     

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.     

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.     

To the theory of homogeneous nucleation: Cluster energy

An attempt is made to critically analyze the modern state of the theory of homogeneous nucleation as concerns its ability to describe experiments with high accuracy. An analysis of the experimental data led us to conclude that the dependence of the nucleation rate on supersaturation and temperature T was not described by the theory, which underestimates the critical cluster size compared with the Gibbs-Thomson equation. The possibility of applying density functional theory (one of the latest achievements in the theory of homogeneous nucleation) was questioned. Within this theory, the Gibbs-Thomson equation remains valid even outside the classic capillary approximation. It is suggested that, to bring theory in consistency with experiment, certain fundamental propositions of the theory of nucleation should be revised. The inclusion of an additional contribution to the Gibbs energy of a cluster caused by the size dependence of the specific heat capacity of the cluster decreases the critical cluster size compared with the value calculated by the Gibbs-Thomson equation. The calculated dependence of nucleation rate on supersaturation was in agreement with the experimental results.     

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.     

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.     

Inverse density-functional theory as an interpretive tool for measuring colloid-surface interactions in dense systems

Recent advances in optical microscopy, such as total internal reflection and confocal scanning laser techniques, now permit the direct three-dimensional tracking of large numbers of colloidal particles both near and far from interfaces. A novel application of this technology, currently being developed by one of the authors under the name of diffusing colloidal probe microscopy (DCPM), is to use colloidal particles as probes of the energetic characteristics of a surface. A major theoretical challenge in implementing DCPM is to obtain the potential energy of a single particle in the external field created by the surface, from the measured particle trajectories in a dense colloidal system. In this paper we develop an approach based on an inversion of density-functional theory (DFT), where we calculate the single-particle-surface potential from the experimentally measured equilibrium density profile in a nondilute colloidal fluid. The underlying DFT formulation is based on the recent work of Zhou and Ruckenstein [Zhou and Ruckenstein, J. Chem. Phys. 112, 8079 (2000)]. For model hard-sphere and Lennard-Jones systems, using Monte Carlo simulation to provide the "experimental" density profiles, we found that the inversion procedure reproduces the true particle-surface-potential energy to an accuracy within typical DCPM experimental limitations (approximately 0.1 kT) at low to moderate colloidal densities. The choice of DFT closures also significantly affects the accuracy.     

Binary homogeneous nucleation in water-succinic acid and water-glutaric acid systems

Binary homogeneous nucleation of water-succinic acid and water-glutaric acid systems have been investigated. The numerical approach was based on the classical nucleation theory. Usually, nucleation is discussed in terms of kinetics, but the thermodynamics involved is undoubtedly equally important. In this paper we studied the above mentioned binary systems giving a quantitative insight into the nucleation process and a detailed consideration of the thermodynamics involved. Both diacids in study are in solid state at room temperature. They behave in environment according to their liquid state properties because of the absence of crystalline lattice energies, and therefore their subcooled liquid state thermodynamics have to be considered. The lack of consistent thermodynamic data for pure organic components and their aqueous solutions represent a high source of uncertainty. However, the present simulations indicate that in atmospheric conditions these binary systems will not form new particles.     

Curvature dependence of surface free energy of liquid drops and bubbles: A simulation study

We study the excess free energy due to phase coexistence of fluids by Monte Carlo simulations using successive umbrella sampling in finite L×L×L boxes with periodic boundary conditions. Both the vapor-liquid phase coexistence of a simple Lennard-Jones fluid and the coexistence between A-rich and B-rich phases of a symmetric binary (AB) Lennard-Jones mixture are studied, varying the density ρ in the simple fluid or the relative concentration x(A) of A in the binary mixture, respectively. The character of phase coexistence changes from a spherical droplet (or bubble) of the minority phase (near the coexistence curve) to a cylindrical droplet (or bubble) and finally (in the center of the miscibility gap) to a slablike configuration of two parallel flat interfaces. Extending the analysis of Schrader et al., [Phys. Rev. E 79, 061104 (2009)], we extract the surface free energy γ(R) of both spherical and cylindrical droplets and bubbles in the vapor-liquid case and present evidence that for R→∞ the leading order (Tolman) correction for droplets has sign opposite to the case of bubbles, consistent with the Tolman length being independent on the sign of curvature. For the symmetric binary mixture, the expected nonexistence of the Tolman length is confirmed. In all cases and for a range of radii R relevant for nucleation theory, γ(R) deviates strongly from γ(∞) which can be accounted for by a term of order γ(∞)/γ(R)-1∝R(-2). Our results for the simple Lennard-Jones fluid are also compared to results from density functional theory, and we find qualitative agreement in the behavior of γ(R) as well as in the sign and magnitude of the Tolman length.     

On the closure conjectures for the Gibbsian approximation model of a binary droplet

Within the framework of Gibbsian thermodynamics, a binary droplet is regarded to consist of a uniform interior and dividing surface. The properties of the droplet interior are those of the bulk liquid solution, but the dividing surface is a fictitious phase whose chemical potentials cannot be rigorously determined. The state of the nucleus interior and free energy of nucleus formation can be found without knowing the surface chemical potentials, but the latter are still needed to determine the state of the whole nucleus (including the dividing surface) and develop the kinetics of nucleation. Thus it is necessary to recur to additional conjectures in order to build a complete, thermodynamic, and kinetic theory of nucleation within the framework of the Gibbsian approximation. Here we consider and analyze the problem of closing the Gibbsian approximation droplet model. We identify micro- and Gamma-closure conjectures concerning the surface chemical potentials and excess surface coverages, respectively, for the droplet surface of tension. With these two closure conjectures, the Gibbsian approximation model of a binary droplet becomes complete so that one can determine both the surface and internal characteristics of the whole nucleus and develop the kinetic theory, based on this model. Theoretical results are illustrated by numerical evaluations for binary nucleation in a water-methanol vapor mixture at T=298.15 K. Numerical results show a striking increase in the droplet surface tension with decreasing droplet size at constant overall droplet composition. A comparison of the Gibbsian approximation with density functional calculations for a model surfactant system indicate that the excess surface coverages from the Gibbsian approximation are accurate enough for large droplets and droplets that are not too concentrated with respect to the solute.