Fundamentals of the Diffuse Interface Theory of Nucleation

Thermodynamic Theory of Curvature-Dependent Surface Tension: Tolman's Theory Revisited

An exact equation for determining the Tolman length (TL) as a function of radius is obtained within the framework of classical thermodynamics and a computational procedure for solving it is proposed. As a result of implementing this procedure, the dependences of the TL and surface tension on radius are obtained for the drop and bubble cases and various equations of state. As one of the results of the thermodynamic study, a new equation for the dependence of surface tension on radius (curvature effect) alternative to the corresponding Tolman equation and associated with the spinodal point is obtained. The fundamental impossibility to determine the curvature effect analytically from the binodal point, i.e., using the Tolman equation, is established; it is calculated only from the spinodal point and is determined by the characteristics of the system at this point. The sign of the TL asymptotic value debated in the literature in recent decades is uniquely determined in the theory: it is negative for drops and positive for bubbles.

Nucleation of melt: From fundamentals to dispersed systems

The most evident aspects of a first order transition of a system from an old to a new phase, are the presence of a discontinuity at the interface between both phases and the thermal effects related to the latent heat exchanged with the surrounding environment. These effects are the result of a sequence of events promoted by thermodynamic conditions persisting over the equilibrium in a metastable state. The breakdown of metastability is promoted by infinitesimal energy fluctuations resulting in the germination of clusters of the new phase that can grow to a critical size (nucleus) and then develop or vanish. Examples of these sequences are common in various technological fields such as combustion, food processing, pharmaceutical manufacturing, condensation, and phase change heat transfer, etc. This work aims to highlight a logical path that leads the readers from the fundamental phenomenology to the most intricated aspects of the nucleation within dispersed systems such as oil-in-water emulsions. Differences between the homogeneous and heterogeneous mechanisms are, under the light of the Classical Nucleation Theory (CNT), presented in bulk and confined systems until defining a minimum confinement size. By collecting insights coming from a rich scientific literature mostly focused on the stability of emulsified systems, the discussion is then on the aspects related to the surface related mechanisms. Two main aspects are then considered: a) the wettability of the nucleating cluster by the surrounding melt; b) the affinity between the adsorbed layer, where a surfactant is located, and the oil melt phase (mainly n-alkanes and triacylglycerols with different moieties). In cases where nucleation is dominating over the dewetting of the nucleus, the contact angle can be considered as a constant value. The affinity in terms of molecular features between the surfactant and the oil phase can promote the template effect. Several factors seem to play a role in this interaction such as the thermal characteristics of the surfactant and comparable dimensions between the molecule (or fractions) of the dispersed compound and the tail of the surfactant.

A self-consistent phase-field approach to implicit solvation of charged molecules with Poisson-Boltzmann electrostatics

Dielectric boundary based implicit-solvent models provide efficient descriptions of coarse-grained effects, particularly the electrostatic effect, of aqueous solvent. Recent years have seen the initial success of a new such model, variational implicit-solvent model (VISM) [Dzubiella, Swanson, and McCammon Phys. Rev. Lett. 96, 087802 (2006) and J. Chem. Phys. 124, 084905 (2006)], in capturing multiple dry and wet hydration states, describing the subtle electrostatic effect in hydrophobic interactions, and providing qualitatively good estimates of solvation free energies. Here, we develop a phase-field VISM to the solvation of charged molecules in aqueous solvent to include more flexibility. In this approach, a stable equilibrium molecular system is described by a phase field that takes one constant value in the solute region and a different constant value in the solvent region, and smoothly changes its value on a thin transition layer representing a smeared solute-solvent interface or dielectric boundary. Such a phase field minimizes an effective solvation free-energy functional that consists of the solute-solvent interfacial energy, solute-solvent van der Waals interaction energy, and electrostatic free energy described by the Poisson-Boltzmann theory. We apply our model and methods to the solvation of single ions, two parallel plates, and protein complexes BphC and p53/MDM2 to demonstrate the capability and efficiency of our approach at different levels. With a diffuse dielectric boundary, our new approach can describe the dielectric asymmetry in the solute-solvent interfacial region. Our theory is developed based on rigorous mathematical studies and is also connected to the Lum-Chandler-Weeks theory (1999). We discuss these connections and possible extensions of our theory and methods.

Pattern formation by phase-field relaxation of bending energy with fixed surface area and volume

We explore a wide variety of patterns of closed surfaces that minimize the elastic bending energy with fixed surface area and volume. To avoid complicated discretization and numerical instabilities for sharp surfaces, we reformulate the underlying constrained minimization problem by constructing phase-field functionals of bending energy with penalty terms for the constraints and develop stable numerical methods to relax these functionals. We report our extensive computational results with different initial surfaces. These results are discussed in terms of the reduced volume and are compared with the known results obtained using the sharp-interface approach. Finally, we discuss the implications of our numerical findings.

Phase-field approach to implicit solvation of biomolecules with Coulomb-field approximation

A phase-field variational implicit-solvent approach is developed for the solvation of charged molecules. The starting point of such an approach is the representation of a solute-solvent interface by a phase field that takes one value in the solute region and another in the solvent region, with a smooth transition from one to the other on a small transition layer. The minimization of an effective free-energy functional of all possible phase fields determines the equilibrium conformations and free energies of an underlying molecular system. All the surface energy, the solute-solvent van der Waals interaction, and the electrostatic interaction are coupled together self-consistently through a phase field. The surface energy results from the minimization of a double-well potential and the gradient of a field. The electrostatic interaction is described by the Coulomb-field approximation. Accurate and efficient methods are designed and implemented to numerically relax an underlying charged molecular system. Applications to single ions, a two-plate system, and a two-domain protein reveal that the new theory and methods can capture capillary evaporation in hydrophobic confinement and corresponding multiple equilibrium states as found in molecular dynamics simulations. Comparisons of the phase-field and the original sharp-interface variational approaches are discussed.

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: I. Transitions in the one-phase liquid region

The phase field theory (PFT) has been applied to predict equilibrium interfacial properties and nucleation barrier in the binary eutectic system Ag-Cu using double well and interpolation functions deduced from a Ginzburg-Landau expansion that considers fcc (face centered cubic) crystal symmetries. The temperature and composition dependent free energies of the liquid and solid phases are taken from CALculation of PHAse Diagrams-type calculations. The model parameters of PFT are fixed so as to recover an interface thickness of approximately 1 nm from molecular dynamics simulations and the interfacial free energies from the experimental dihedral angles available for the pure components. A nontrivial temperature and composition dependence for the equilibrium interfacial free energy is observed. Mapping the possible nucleation pathways, we find that the Ag and Cu rich critical fluctuations compete against each other in the neighborhood of the eutectic composition. The Tolman length is positive and shows a maximum as a function of undercooling. The PFT predictions for the critical undercooling are found to be consistent with experimental results. These results support the view that heterogeneous nucleation took place in the undercooling experiments available at present. We also present calculations using the classical droplet model [classical nucleation theory (CNT)] and a phenomenological diffuse interface theory (DIT). While the predictions of the CNT with a purely entropic interfacial free energy underestimate the critical undercooling, the DIT results appear to be in a reasonable agreement with the PFT predictions.

Analytical density functional theory of homogeneous vapor condensation

Starting from an exact gradient transcription of the perturbative density functional theory of homogeneous vapor condensation, we propose an analytical approximation that reproduces the density profile and free energy of critical fluctuations to high accuracy. For a broad variety of substances, including nonpolar, weakly polar, and metallic liquids, the method predicts nucleation rates that are orders of magnitude closer to experiment than those from the classical approach. The present treatment incorporates detailed molecular theory into macroscopic modeling.

STRUCTURE AND TRANSFORMATION: Large Molecular Clusters as Models of Condensed Matter

This paper reviews investigations of homogeneous nucleation in phase transitions in large molecular clusters. The principal techniques brought to bear are electron diffraction analyses of transformations in clusters formed by condensation of vapor in supersonic expansions and computer simulations of spontaneous phase changes in clusters. Results obtained to date are contrasted with those of larger systems and interpreted in terms of nucleation theory. The review also refers to some unresolved aspects of nucleation theory.