Nucleation near the spinodal: Limitations of mean field density functional theory

Free energy of critical droplets-from the binodal to the spinodal

Arguably, the main challenge of nucleation theory is to accurately evaluate the work of formation of a critical embryo in the new phase, which governs the nucleation rate. In Classical Nucleation Theory (CNT), this work of formation is estimated using the capillarity approximation, which relies on the value of the planar surface tension. This approximation has been blamed for the large discrepancies between predictions from CNT and experiments. In this work, we present a study of the free energy of formation of critical clusters of the Lennard-Jones fluid truncated and shifted at 2.5σ using Monte Carlo simulations, density gradient theory, and density functional theory. We find that density gradient theory and density functional theory accurately reproduce molecular simulation results for critical droplet sizes and their free energies. The capillarity approximation grossly overestimates the free energy of small droplets. The incorporation of curvature corrections up to the second order with the Helfrich expansion greatly remedies this and performs very well for most of the experimentally accessible regions. However, it is imprecise for the smallest droplets and largest metastabilities since it does not account for a vanishing nucleation barrier at the spinodal. To remedy this, we propose a scaling function that uses all relevant ingredients without adding fitting parameters. The scaling function reproduces accurately the free energy of the formation of critical droplets for the entire metastability range and all temperatures examined and deviates from density gradient theory by less than one k_BT.

Formation free energy of an i-mer at spinodal

In statistical mechanics, the formation free energy of an i-mer can be understood as the Gibbs free energy change in a system consisting of pure monomers after and prior to the formation of the i-mer. For molecules interacting via Lennard-Jones potential, we have computed the formation free energy of a Stillinger i-mer [F. H. Stillinger, J. Chem. Phys. 38, 1486 (1963)] and a ten Wolde-Frenkel (tWF) [P. R. ten Wolde and D. Frenkel, J. Chem. Phys. 109, 9901 (1998)] i-mer at spinodal at reduced temperatures from 0.7 to 1.2. It turns out that the size of a critical Stillinger i-mer remains finite and its formation free energy is on the order of k_BT, and the size of a critical tWF i-mer remains finite and its formation free energy is even higher. This can be explained by Binder's theory [K. Binder, Phys. Rev. A 29, 341 (1984)] that for a system, when approaching spinodal, if the Ginzburg criterion is not satisfied, a gradual transition will take place from nucleation to spinodal decomposition, where the free-energy barrier height is on the order of k_BT.

Predictions of homogeneous nucleation rates for n-alkanes accounting for the diffuse phase interface and capillary waves

Homogeneous droplet nucleation has been studied for almost a century but has not yet been fully understood. In this work, we used the density gradient theory (DGT) and considered the influence of capillary waves (CWs) on the predicted size-dependent surface tensions and nucleation rates for selected n-alkanes. The DGT model was completed by an equation of state (EoS) based on the perturbed-chain statistical associating fluid theory and compared to the classical nucleation theory and the Peng-Robinson EoS. It was found that the critical clusters are practically free of CWs because they are so small that even the smallest wavelengths of CWs do not fit into their finite dimensions. The CWs contribute to the entropy of the system and thus decrease the surface tension. A correction for the effect of CWs on the surface tension is presented. The effect of the different EoSs is relatively small because by a fortuitous coincidence their predictions are similar in the relevant range of critical cluster sizes. The difference of the DGT predictions to the classical nucleation theory computations is important but not decisive. Of the effects investigated, the most pronounced is the suppression of CWs which causes a sizable decrease of the predicted nucleation rates. The major difference between experimental nucleation rate data and theoretical predictions remains in the temperature dependence. For normal alkanes, this discrepancy is much stronger than observed, e.g., for water. Theoretical corrections developed here have a minor influence on the temperature dependency. We provide empirical equations correcting the predicted nucleation rates to values comparable with experiments.

Lifetime of metastable states in a Ginzburg-Landau system: Numerical simulations at large driving forces

We developed a "brute-force" simulation method and conducted numerical "experiments" on homogeneous nucleation in an isotropic system at large driving forces (not small supersaturations) using the stochastic Ginzburg-Landau approach. Interactions in the system are described by the asymmetric (no external field), athermal (temperature-independent driving force), tangential (simple phase diagram) Hamiltonian, which has two independent "drivers" of the phase transition: supersaturation and thermal noise. We obtained the probability distribution function of the lifetime of the metastable state and analyzed its mean value as a function of the supersaturation, noise strength, and volume. We also proved the nucleation theorem in the mean-field approximation. The results allowed us to find the thermodynamic properties of the barrier state and conclude that at large driving forces the fluctuating volumes are not independent.

On the interfacial thermodynamics of nanoscale droplets and bubbles

We present a new self-consistent thermodynamic formalism for the interfacial properties of nanoscale embryos whose interiors do not exhibit bulklike behavior and are in complete equilibrium with the surrounding mother phase. In contrast to the standard Gibbsian analysis, whereby a bulk reference pressure based on the same temperature and chemical potentials of the mother phase is introduced, our approach naturally incorporates the normal pressure at the center of the embryo as an appropriate reference pressure. While the interfacial properties of small embryos that follow from the use of these two reference pressures are different, both methods yield by construction the same reversible work of embryo formation as well as consistency between their respective thermodynamic and mechanical routes to the surface tension. Hence, there is no a priori reason to select one method over another. Nevertheless, we argue, and demonstrate via a density-functional theory (with the local density approximation) analysis of embryo formation in the pure component Lennard-Jones fluid, that our new method generates more physically appealing trends. For example, within the new approach the surface tension at all locations of the dividing surface vanishes at the spinodal where the density profile spanning the embryo and mother phase becomes completely uniform (only the surface tension at the Gibbs surface of tension vanishes in the Gibbsian method at this same limit). Also, for bubbles, the location of the surface of tension now diverges at the spinodal, similar to the divergent behavior exhibited by the equimolar dividing surface (in the Gibbsian method, the location of the surface of tension vanishes instead). For droplets, the new method allows for the appearance of negative surface tensions (the Gibbsian method always yields positive tensions) when the normal pressures within the interior of the embryo become less than the bulk pressure of the surrounding vapor phase. Such a prediction, which is allowed by thermodynamics, is consistent with the interpretation that the mother phase's attempted compression of the droplet is counterbalanced by the negative surface tension, or free energy cost to decrease the interfacial area. Furthermore, for these same droplets, the surface of tension can no longer be meaningfully defined (the surface of tension always remains well defined in the Gibbsian method). Within the new method, the dividing surface at which the surface tension equals zero emerges as a new lengthscale, which has various thermodynamic analogs to and similar behavior as the surface of tension.

Stability of critical bubble in stretched fluid of square-gradient density-functional model with triple-parabolic free energy

The square-gradient density-functional model with triple-parabolic free energy, which was used previously to study the homogeneous bubble nucleation [M. Iwamatsu, J. Chem. Phys. 129, 104508 (2008)], is used to study the stability of the critical bubble nucleated within the bulk undersaturated stretched fluid. The stability of the bubble is studied by solving the Schrodinger equation for the fluctuation. The negative eigenvalue corresponds to the unstable growing mode of the fluctuation. Our results show that there is only one negative eigenvalue whose eigenfunction represents the fluctuation that corresponds to the isotropically growing or shrinking nucleus. In particular, this negative eigenvalue survives up to the spinodal point. Therefore, the critical bubble is not fractal or ramified near the spinodal.

Low-density/high-density liquid phase transition for model globular proteins

The effect of molecule size (excluded volume) and the range of interaction on the surface tension, phase diagram, and nucleation properties of a model globular protein is investigated using a combination of Monte Carlo simulations and finite temperature classical density functional theory calculations. We use a parametrized potential that can vary smoothly from the standard Lennard-Jones interaction characteristic of simple fluids to the ten Wolde-Frenkel model for the effective interaction of globular proteins in solution. We find that the large excluded volume characteristic of large macromolecules such as proteins is the dominant effect in determining the liquid-vapor surface tension and nucleation properties. The variation of the range of the potential is important only in the case of small excluded volumes such as for simple fluids. The DFT calculations are then used to study the homogeneous nucleation of the high-density phase from the low-density phase including the nucleation barriers, nucleation pathways, and rate. It is found that the nucleation barriers are typically only a few k(B)T and that the nucleation rates are substantially higher than would be predicted by classical nucleation theory.

Crossover from nucleation to spinodal decomposition in a condensing vapor

The mechanism controlling the initial step of a phase transition has a tremendous influence on the emerging phase. We study the crossover from a purely nucleation-controlled transition toward spinodal decomposition in a condensing Lennard-Jones vapor using molecular dynamics simulations. We analyze both the kinetics and at the same time the thermodynamics by directly reconstructing the free energy of cluster formation. We estimate the location of the spinodal, which lies at much deeper supersaturations than expected. Moreover, the nucleation barriers we find differ only by a constant from the classical nucleation theory predictions and are in very good agreement with semiempirical scaling relations. In the regime from very small barriers to the spinodal, growth controls the rate of the transition but not its nature because the activation barrier has not yet vanished. Finally, we discuss in detail the influence of the chosen reaction coordinate on the interpretation of such simulation results.

A limit of stability in supercooled liquid clusters

We examine the metastable liquid phase of a supercooled gold nanocluster by studying the free energy landscape using the largest solidlike embryo as an order parameter. Just below freezing, the free energy exhibits a local minimum at small embryo sizes and a maximum at a larger critical embryo size. At T=660 K the free energy becomes a monotonically decreasing function of the order parameter as the liquid phase becomes unstable, indicating that we have reached a limit of stability. In contrast to the mean-field theory predictions for a spinodal, the size of the critical embryo remains finite as the limit of stability is approached. We also calculate the rate of nucleation, independently from our free energy calculations, and observe a rapid increase in its temperature dependence when the free energy barrier is on the order of kT. We suggest that this supports the idea that freezing becomes a barrierless process at low temperatures.

Elucidating the mechanism of nucleation near the gas-liquid spinodal

We have constructed a multidimensional free energy surface of nucleation of the liquid phase from the parent supercooled and supersaturated vapor phase near the gas-liquid spinodal. In particular, we remove the Becker-Doring constraint of having only one growing cluster in the system. Close to the spinodal, the free energy, as a function of the size of the largest cluster, develops surprisingly a minimum at a subcritical cluster size. It is this minimum at intermediate size that is found to be responsible for the barrier towards further growth of the nucleus at large supersaturation. An alternative free energy pathway involving the participation of many subcritical clusters is found near the spinodal where the growth of the nucleus is promoted by a coalescence mechanism. The growth of the stable phase becomes collective and spatially diffuse, and the significance of a "critical nucleus" is lost for deeper quenches.

Bubble formation in lattice Boltzmann immiscible shear flow

Bubble formation and the effect of shear on bubble formation in a van der Waals fluid is investigated by means of lattice Boltzmann mesoscale simulations. In the absence of shear, the maximum number of formed bubbles increases with undercooling but the incubation time before bubble formation decreases dramatically. The results are in agreement with classical phase transition theory. In shear flow, the maximum number of bubbles is not affected by shear but the bubble growth rate is accelerated. The effect of shear on bubble growth rate weakens at large undercoolings. The reasons are twofold. On the one hand the highly undercooled system takes less time to complete phase transition due to the large driving force so that there is less time to accumulate the flow effect. On the other hand the mechanism for bubble growth changes from coarsening to coalescence at large undercoolings.

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.

How to extend hard sphere density functional approximation to nonuniform nonhard sphere fluids: Applicable to both subcritical and supercritical temperature regions

A methodology for the formulation of density functional approximation (DFA) for nonuniform nonhard sphere fluids is proposed by following the spirit of a partitioned density functional approximation [Zhou, Phys. Rev. E 68, 061201 (2003)] and mapping the hard core part onto an effective hard sphere whose high order part of the functional perturbation expansion is treated by existing hard sphere DFAs. The resultant density functional theory (DFT) formalism only needs a second order direct correlation function and pressure of the corresponding coexistence bulk fluid as inputs and therefore can be applicable to both supercritical and subcritical temperature cases. As an example, an adjustable parameter-free version of a recently proposed Lagrangian theorem-based DFA is imported into the present methodology; the resultant DFA is applied to Lennard-Jones fluid under the influence of external fields due to a single hard wall, two hard walls separated by a small distance, a large hard sphere, and a spherical cavity with a hard wall. By comparing theoretical predictions with previous simulation data and those recently supplied for coexistence bulk fluid situated at "dangerous" regions, it was found that the present DFA can predict subtle structure change of the density profile and therefore is the most accurate among all existing DFT approaches. A detailed discussion is given as to why so excellent DFA for nonhard sphere fluids can be drawn forth from the present methodology and how the present methodology differs from previous ones. The methodology can be universal, i.e., it can be combined with any other hard sphere DFAs to construct DFA for other nonhard sphere fluids with a repulsive core.

Structures and adsorption of binary hard-core Yukawa mixtures in a slitlike pore: Grand canonical Monte Carlo simulation and density-functional study

The grand canonical ensemble Monte Carlo simulation and density-functional theory are applied to calculate the structures, local mole fractions, and adsorption isotherms of binary hard-core Yukawa mixtures in a slitlike pore as well as the radial distribution functions of bulk mixtures. The excess Helmholtz energy functional is a combination of the modified fundamental measure theory of Yu and Wu [J. Chem. Phys. 117, 10156 (2002)] for the hard-core contribution and a corrected mean-field theory for the attractive contribution. A comparison of the theoretical results with the results from the Monte Carlo simulations shows that the corrected theory improves the density profiles of binary hard-core Yukawa mixtures in the vicinity of contact over the original mean-field theory. Both the present corrected theory and the simulations suggest that depletion and desorption occur at low temperature, and the local segregation can be observed in most cases. For binary mixtures in the hard slitlike pore, the present corrected theory predicts more accurate surface excesses than the original one does, while in the case of the attractive pore, no improvement is found in the prediction of a surface excess of the smaller molecule.