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.     

Nucleation rates in a new phenomenological model

In this paper we develop a new theory to evaluate the nucleation rate in the framework of the EMLD-DNT model. Beyond the model, our theory deals with cluster translation and exclusion, effects that have been virtually ignored in classical nucleation theory. We apply the model to the case of 1-pentanol, and compare the predictions with experimental results. We find an excellent agreement between the nucleation rate predicted by our theory and experimental data. The distinguishing feature of the model is its ability to predict successfully the rate of formation of the critical nucleus without the use of an intermolecular potential, employing only macroscopic thermodynamic properties.     

Accounting for fluctuations in cluster parameters upon the description of nucleation in supersaturated vapor

A theory is proposed for stationary homogeneous nucleation in supersaturated vapor in which a modified expression for the rate of cluster evaporation was used to calculate the equilibrium distribution over the nucleus sizes and the rates of their formation. This rate was determined by the extrapolation to the region of small sizes of the corresponding expression for the macroscopic droplet derived according to thermodynamic notions that take fluctuations into account. Modified dependences of the size of critical nucleus and the rate of nucleation on the supersaturation and the temperature are determined and compared with the data of the classical theory of nucleation and experimental results.     

Crossover model for the work of critical cluster formation in nucleation theory

We propose a relation for the work of critical cluster formation in nucleation theory W for the systems with long-range interparticle interactions. The method of bridge functions is used to combine the system behavior at sufficiently small quenches, adequately predicted by the classical nucleation theory, with nonclassical effects at deep quenches in the vicinity of the thermodynamic spinodal, described within the framework of the field theoretical approach with an appropriate Ginzburg-Landau functional. The crossover between the two types of nucleation behavior takes place in the vicinity of the kinetic spinodal where the lifetime of a metastable state is of the order of the relaxation time to local equilibrium. We argue that the kinetic spinodal corresponds to the minimum of the excess number of molecules in the critical cluster. This conjecture leads to the form of W containing no adjustable parameters. The barrier scaling function Gamma = W/W(cl), where W(cl) is the classical nucleation barrier, depends parametrically on temperature through the dimensionless combination of material properties. The results for argon nucleation are presented.     

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.     

Hydrogen Bonding and Vaporization Thermodynamics in Hexafluoroisopropanol-Acetone and -Methanol Mixtures. A Joined Cluster Analysis and Molecular Dynamic Study

Binary mixtures of hexafluoroisopropanol with either methanol or acetone are analyzed via classical molecular dynamics simulations and quantum cluster equilibrium calculations. In particular, their populations and thermodynamic properties are investigated with the binary quantum cluster equilibrium method, using our in-house code Peacemaker 2.8, upgraded with temperature-dependent parameters. A novel approach, where the final density from classical molecular dynamics, has been used to generate the necessary reference isobars. The hydrogen bond network in both type of mixtures at molar fraction of hexafluoroisopropanol of 0.2, 0.5, and 0.8 respectively is investigated via the molecular dynamics trajectories and the cluster results. In particular, the populations show that mixed clusters are preferred in both systems even at 0.2 molar fractions of hexafluoroisopropanol. Enthalpies and entropies of vaporization are calculated for the neat and mixed systems and found to be in good agreement with experimental values.     

Quasi-unary homogeneous nucleation of H2SO4-H2O

We show that the binary homogeneous nucleation (BHN) of H2SO4-H2O can be treated as quasi-unary nucleation of H2SO4 in equilibrium with H2O vapor. A scheme to calculate the evaporation coefficient of H2SO4 molecules from H2SO4-H2O clusters is presented and a kinetic model to simulate the quasi-unary nucleation of H2SO4-H2O is developed. In the kinetic model, the growth and evaporation of sulfuric acid clusters of various sizes are explicitly simulated. The kinetic quasi-unary nucleation model does not have two well-recognized problems associated with the classical BHN theory (violation of the mass action law and mismatch of the cluster distribution for monomers) and is appropriate for the situations where the assumption of equilibrium cluster distribution is invalid. The nucleation rates predicted with our quasi-unary kinetic model are consistent with recent experimental nucleation experiments in all the cases studied, while the most recent version of the classical BHN model systematically overpredicts the nucleation rates. The hydration of sulfuric acid clusters, which is not considered in the classical model but is accounted for implicitly in our kinetic quasi-unary model, is likely to be one of physical mechanisms that lead to lower nucleation rates. Further investigation is needed to understand exactly what cause the difference between the kinetic quasi-unary model and the classical BHN model.     

The molecular approach to heterogeneous nucleation

A molecular approach to heterogeneous nucleation has been developed. The expressions for the equilibrium cluster distribution, the reversible work of the cluster formation, and the nucleation rate have been derived. Two separate statements for the work of formation were formulated. If the equilibrium cluster distribution is normalized on the monomer concentration near the substrate surface, the reversible work of formation is expressed by DeltaG(het) (I) = (F(n) (het)-F(n) (hom))-(F(1) (het)-F(1) (hom)) + DeltaG(hom) where F(n) (het) and F(n) (hom) are the Helmholtz free energies of a cluster interacting with a substrate and a cluster not interacting with the substrate, respectively. If the equilibrium cluster distribution is normalized on the monomer concentration far from the substrate surface, the work of cluster formation is given by DeltaG(het) (II) = (F(n) (het)-F(n) (hom)) + DeltaG(hom). The former expression corresponds to the approach of the classical heterogeneous nucleation theory. The cluster partition function appears to be dependent on the location of a virtual plane, which separates the volume, where the interaction of the clusters with the substrate is effective from the one where interaction is negligible. Our Monte Carlo simulations have shown that the dependence is rather weak and thus the location of the plane is not very important. According to the simulations the variation of the plane position in the range from 20 to 50 Angstroms does not lead to a considerable change of the heterogeneous nucleation rate.     

Extension of the associated perturbed anisotropic chain theory to mixtures with more than one associating component

The Associated Perturbed Anisotropic Chain Theory (APACT), a theory useful for calculation of equilibrium thermodynamic properties for fluid mixtures presented recently by Ikonomou and Donohue has been extended to treat multicomponent mixtures in which two or more components associate. Mixtures containing non-associating components (diluents) also can be treated. Association equilibria for each component hydrogen bonding to itself as well as for cross-association between associating components are considered. Following an approach similar to that used previously in the derivation of APACT, a closed-form equation of state has been derived. Preliminary results show that APACT does well in fitting phase equilibrium data for binary mixtures involving alcohols and water with small values of k_ij. The fits are not as good, however, for systems containing small molecules (water and methanol) and the reasons for this shortcoming are discussed.     

由成核率数据拟合饱和蒸汽压力

Saturated vapor pressure was calculated from the nucleation experimental data using the thermodynamically consistent nucleation theory in which the effect of real gas is considered. The cubic polynomial fit equations of saturation pressure for several substances were obtained based on the calculation. The results of the calculations were compared to those of thermodynamic equilibrium equation and the empirical equation and applied to the predictions of the classical nucleation theory. The results show that the saturation pressures estimated from the nucleation data agree fairly well with those of empirical equations for the substances investigated, and this indicates that the predictions from the classical nucleation theory are close to the experimental data.     

Application of scaling and kinetic equations to helium cluster size distributions: Homogeneous nucleation of a nearly ideal gas

A previously published model of homogeneous nucleation [Villarica et al., J. Chem. Phys. 98, 4610 (1993)] based on the Smoluchowski [Phys. Z. 17, 557 (1916)] equations is used to simulate the experimentally measured size distributions of 4He clusters produced in free jet expansions. The model includes only binary collisions and does not consider evaporative effects, so that binary reactive collisions are rate limiting for formation of all cluster sizes despite the need for stabilization of nascent clusters. The model represents these data very well, accounting in some cases for nearly four orders of magnitude in variation in abundance over cluster sizes ranging up to nearly 100 atoms. The success of the model may be due to particularities of 4He clusters, i.e., their very low coalescence exothermicity, and to the low temperature of 6.7 K at which the data were collected.     

Statistical modeling of equilibrium adsorption of non-ideal mixtures on zeolites

A new model for equilibrium adsorption of a binary mixture on zeolites that takes into account the energy nonuniformity of the adsorption field in the zeolite cavities was developed on the basis of statistical thermodynamics. The nonuniformity of the adsorption field produces rearrangement of molecules in the cavity volume, decreasing the entropy, internal energy, and Helmholz free energy. A procedure for calculation of the thermodynamic functions from the data on the adsorption of pure components was proposed. The limiting cases of maximum ordering of the molecules in the cavity and their random distribution were considered. The approach proposed was exemplified by the substantially non-ideal system nitrogen--argon--zeolite NaX at 160 K. The proposed model describes the behavior of this mixture much better than that of the ideal adsorbed solution theory.     

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.     

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.     

A multigrid method for N-component nucleation

A multigrid algorithm has been developed enabling more efficient solution of the cluster size distribution for N-component nucleation from the Becker-D?ring equations. The theoretical derivation is valid for an arbitrary number of condensing components, making the simulation of many-component nucleating systems feasible. A steady state ternary nucleation problem is defined to demonstrate its efficiency. The results are used as a validation for existing nucleation theories. The non-steady state ternary problem provides useful insight into the initial stages of the nucleation process. We observe that for the ideal mixture the main nucleation flux bypasses the saddle point.     

Prediction of thermodynamic behavior of binary mixtures by applying the crossover approach to Soave–Redlich–Kwong equation of state

Thermodynamic analysis of binary mixtures near the critical region is essential for many chemical process designs. In this research, based on isomorphism principle and incorporating general crossover approach the original Soave–Redlich–Kwong (SRK) equation of state (EOS) was developed for the binary mixtures. We have introduced an additional term in the crossover function in order to take into account the difference between the classical critical parameters and the real critical parameters. The applicability of this crossover EOS was verified against methane–ethane mixture to describe their thermodynamic properties over a wide range of thermodynamic states, because of their wide applications. It is shown that based on this approach we can received too much more accuracy for predicting thermodynamic properties in comparison with classical form of SRK EOS.     

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.