Multicomponent nucleation: Thermodynamically consistent description of the nucleation work

Hierarchical approximations to the nucleation work in the entire range of metastability

The work W to form a nucleus (also known as the critical nucleus) is a key quantity in the description of nucleation phenomena because of its exponentially strong effect on the nucleation rate. The present study provides a general approximate expression for W, which comprises a hierarchy of approximations to the dependence of W on the experimentally controlled overpressure Δp of a nucleating multicomponent phase. This general expression is used to derive explicit formulas for the lowest-order members of the W(Δp) hierarchy as well as for the respective lowest-order approximations to the Δp dependences of the nucleus surface tension, the nucleus radius, the Gibbs-Tolman length, and the stationary nucleation rate. The second-order and the third-order approximations to the W(Δp) dependence are confronted with available W(Δp) data, and the latter is found to agree very well with the data. The results obtained are applicable to homogeneous single-component or multicomponent nucleation from the binodal to the spinodal of the old phase, i.e., in the entire range of the old-phase metastability.

Nucleation work, surface tension, and Gibbs-Tolman length for nucleus of any size

In the framework of the Gibbs approach to nucleation thermodynamics, expressions are derived for the nucleation work, nucleus size, surface tension, and Gibbs-Tolman length in homogeneous single-component nucleation at a fixed temperature. These expressions are in terms of the experimentally controlled overpressure of the nucleating phase and are valid for the entire overpressure range, i.e., for nucleus of any size. Analysis of available data for bubble and droplet nucleation in Lennard-Jones fluid shows that the theory describes well the data by means of a single free parameter, the Gibbs-Tolman length of the planar liquid/vapor interface. It is found that this length is about one-tenth of the Lennard-Jones molecular-diameter parameter and that it is positive for the bubble nucleus and negative for the droplet nucleus. In a sufficiently narrow temperature range, the nucleation work, nucleus radius, scaled surface tension, and Gibbs-Tolman length are apparently universal functions of scaled overpressure.

Theoretical study of vapor-liquid homogeneous nucleation using stability analysis of a macroscopic phase

Stability analysis is generally used to verify that the solution to phase equilibrium calculations corresponds to a stable state (minimum of the free energy). In this work, tangent plane distance analysis for stability of macroscopic mixtures is also used for analyzing the nucleation process, reconciling thus this analysis with classical nucleation theories. In the context of the revised nucleation theory, the driving force and the nucleation work are expressed as a function of the Lagrange multiplier corresponding to the mole fraction constraint from the minimization problem of stability analysis. Using a van der Waals fluid applied to a ternary mixture, Lagrange multiplier properties are illustrated. In particular, it is shown how the Lagrange multiplier value is equal to one on the binodal and spinodal curves at the same time as the driving force of nucleation vanishes on these curves. Finally, it is shown that, on the spinodal curve, the nucleation work from the revised and generalized nucleation theories are characterized by two different local minima from stability analysis, irrespective of any interfacial tension models.

CaCO3 nucleation by cyanobacteria: laboratory evidence for a passive, surface-induced mechanism

Calcite nucleation on the surface of cyanobacteria of the Synechococcus leopoliensis strain PCC 7942 was investigated to assess the influence of photosynthetic uptake of inorganic carbon and active ion exchange processes across the cell membrane on the nucleation and precipitation mechanisms. We performed long-term precipitation experiments at a constant CO(2) level in ambient air by adding suspensions of previously washed cyanobacteria to solutions of NaHCO(3)/CaCl(2) which were supersaturated with respect to calcite. Induction times between 4 and 110 h were measured over a range of saturation states, Omega, between 8 and 4. The kinetics of CaCO(3) nucleation was compared between experiments: (i) with ongoing photosynthesis, (ii) with cells metabolizing but not undergoing photosynthetic uptake of inorganic carbon and (iii) in darkness without photosynthesis. No significant differences were observed between the three treatments. The results reveal that under low nutrient concentrations and permanent CO(2) supply, photosynthetic uptake of inorganic carbon predominantly uses CO(2) and consequently does not directly influence the nucleation process of CaCO(3) at the surface of S. leopoliensis. Furthermore, ion exchange processes did not affect the kinetics, indicating a passive nucleation process wherein the cell surface or extracellular polymers provided preferential sites for mineral nucleation. The catalyzing effect of the cyanobacteria on calcite nucleation was equivalent to a approximately 18% reduction in the specific interfacial free energy of the calcite nuclei. This result and the ubiquitous abundance of cyanobacteria suggest that this process may have an impact on local and global carbon cycling.

Analysis of experimental data for the nucleation rate of water droplets

A formula for the stationary nucleation rate J is proposed and used for analysis of experimental data for the dependence of J on the supersaturation ratio S in isothermal homogeneous nucleation of water droplets in vapors. It is found that the experimental data are described quite successfully by the proposed formula which is based on (i) the Gibbs presentation of the nucleation work in terms of overpressure, (ii) the Girshick-Chiu [J. Chem. Phys. 93, 1273 (1990); 94, 826 (1991)] self-consistency correction to the equilibrium cluster size distribution, and (iii) the Reguera-Rubi [J. Chem. Phys. 115, 7100 (2001)] kinetic accounting of the nucleus translational-rotational motion. The formula, like that of Wolk and Strey [J. Phys. Chem. B 105, 11683 (2001)], could be used as a semiempirical relation describing the J(S) dependence for nucleation in vapors of single-component droplets or crystals of substances with insufficiently well known equations of state.

Forms and applications of the nucleation theorem

The nucleation theorem is a general relation between the nucleation work, the nucleus size, and the supersaturation or other thermodynamic parameters of the old phase. The theorem appears in different forms, depending not only on the chosen set of independent variables describing the nucleation work but also on which of these variables is changed while the others are held fixed. This paper gives a rigorous, systematic, and comprehensive presentation of various forms of the nucleation theorem and shows how some of them can be applied to concrete cases of nucleation. Both theoretical and experimental applications of the theorem to nucleation in unary, binary, or ternary phases are considered.

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.

Comment on “Multicomponent nucleation: Thermodynamically consistent description of the nucleation work” [J. Chem. Phys. , 3749 (2004)]

In the above cited paper, Kashchiev extends the "thermodynamically consistent" description of nucleation work to phase formation processes in multicomponent systems. It is claimed that this approach is generally applicable widely independent on the kind of phase transformation considered. In contrast to this statement, it is shown in our Comment that Kashchiev's approach is theoretically not well founded. It contains, in addition, a number of assumptions, which are not generally fulfilled. Because of this his method cannot give, in general, a satisfactory determination of the work of critical cluster formation in nucleation theory. An experimental example is given verifying our conclusions.

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.