Pore formation in fluctuating membranes

We study the nucleation of a single pore in a fluctuating lipid membrane, specifically taking into account the membrane fluctuations, as well as the shape fluctuations of the pore. For large enough pores, the nucleation free energy is well-described by shifts in the effective membrane surface tension and the pore line tension. Using our framework, we derive the stability criteria for the various pore formation regimes. In addition to the well-known large-tension regime from the classical nucleation theory of pores, we also find a low-tension regime in which the effective line and surface tensions can change sign from their bare values. The latter scenario takes place at sufficiently high temperatures, where the opening of a stable pore of finite size is entropically favorable.     

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.     

Sensitivity of nucleation phenomena on range of interaction potential

Theoretical and computational investigations of nucleation have been plagued by the sensitivity of the phase diagram to the range of the interaction potential. As the surface tension depends strongly on the range of interaction potential and as the classical nucleation theory (CNT) predicts the free energy barrier to be directly proportional to the cube of the surface tension, one expects a strong sensitivity of nucleation barrier to the range of the potential; however, CNT leaves many aspects unexplored. We find for gas-liquid nucleation in Lennard-Jones system that on increasing the range of interaction the kinetic spinodal (KS) (where the mechanism of nucleation changes from activated to barrierless) shifts deeper into the metastable region. Therefore the system remains metastable for larger value of supersaturation and this allows one to explore the high metastable region without encountering the KS. On increasing the range of interaction, both the critical cluster size and pre-critical minima in the free energy surface of kth largest cluster, at respective kinetic spinodals, shift towards smaller cluster size. In order to separate surface tension contribution to the increase in the barrier from other non-trivial factors, we introduce a new scaling form for surface tension and use it to capture both the temperature and the interaction range dependence of surface tension. Surprisingly, we find only a weak non-trivial contribution from other factors to the free energy barrier of nucleation.     

Gradient theory computation of the radius-dependent surface tension and nucleation rate for n-nonane clusters

The Van der Waals-Cahn-Hilliard gradient theory (GT) is applied to determine the structure and the work of formation of clusters in supersaturated n-nonane vapor. The results are analyzed as functions of the difference of pressures of the liquid phase and vapor phase in chemical equilibrium, which is a measure for the supersaturation. The surface tension as a function of pressure difference shows first a weak maximum and then decreases monotonically. The computed Tolman length is in agreement with earlier results [L. Granasy, J. Chem. Phys. 109, 9660 (1998)] obtained with a different equation of state. A method based on the Gibbs adsorption equation is developed to check the consistency of GT results (or other simulation techniques providing the work of formation and excess number of molecules), and to enable an efficient interpolation. A cluster model is devised based on the density profile of the planar phase interface. Using this model we analyze the dependency of the surface tension on the pressure difference. We find three major contributions: (i) the effect of asymmetry of the density profile resulting into a linear increase of the surface tension, (ii) the effect of finite thickness of the phase interface resulting into a negative quadratic term, and (iii) the effect of buildup of a low-density tail of the density profile, also contributing as a negative quadratic term. Contributions (i)-(iii) fully explain the dependency of the surface tension on the pressure difference, including the range relevant to nucleation experiments. Contributions (i) and (ii) can be predicted from the planar density profile. The work of formation of noncritical clusters is derived and the nucleation rate is computed. The computed nucleation rates are closer to the experimental nucleation rate results than the classical Becker-D?ring theory, and also the dependence on supersaturation is better predicted.     

Numerical approaches to determine the interface tension of curved interfaces from free energy calculations

A recently proposed method to obtain the surface free energy σ(R) of spherical droplets and bubbles of fluids, using a thermodynamic analysis of two-phase coexistence in finite boxes at fixed total density, is reconsidered and extended. Building on a comprehensive review of the basic thermodynamic theory, it is shown that from this analysis one can extract both the equimolar radius R(e) as well as the radius R(s) of the surface of tension. Hence the free energy barrier that needs to be overcome in nucleation events where critical droplets and bubbles are formed can be reliably estimated for the range of radii that is of physical interest. It is found that the conventional theory of nucleation, where the interface tension of planar liquid-vapor interfaces is used to predict nucleation barriers, leads to a significant overestimation, and this failure is particularly large for bubbles. Furthermore, different routes to estimate the effective radius-dependent Tolman length δ(R(s)) from simulations in the canonical ensemble are discussed. Thus we obtain an instructive exemplification of the basic quantities and relations of the thermodynamic theory of metastable droplets/bubbles using simulations. However, the simulation results for δ(R(s)) employing a truncated Lennard-Jones system suffer to some extent from unexplained finite size effects, while no such finite size effects are found in corresponding density functional calculations. The numerical results are compatible with the expectation that δ(R(s) → ∞) is slightly negative and of the order of one tenth of a Lennard-Jones diameter, but much larger systems need to be simulated to allow more precise estimates of δ(R(s) → ∞).     

Homogeneous nucleation at high supersaturation and heterogeneous nucleation on microscopic wettable particles: A hybrid thermodynamic/density-functional theory

Homogeneous nucleation at high supersaturation of vapor and heterogeneous nucleation on microscopic wettable particles are studied on the basis of Lennard-Jones model system. A hybrid classical thermodynamics and density-functional theory (DFT) approach is undertaken to treat the nucleation problems. Local-density approximation and weighted-density approximation are employed within the framework of DFT. Special attention is given to the disjoining pressure of small liquid droplets, which is dependent on the thickness of wetting film and radius of the wettable particle. Different contributions to the disjoining pressure are examined using both analytical estimations and numerical DFT calculation. It is shown that van der Waals interaction results in negative contribution to the disjoining pressure. The presence of wettable particles results in positive contribution to the disjoining pressure, which plays the key role in the heterogeneous nucleation. Several definitions of the surface tension of liquid droplets are discussed. Curvature dependence of the surface tension of small liquid droplets is computed. The important characteristics of nucleation, including the formation free energy of the droplet and nucleation barrier height, are obtained.     

Thermodynamics of heterogeneous crystal nucleation in contact and immersion modes

One of the most intriguing problems of heterogeneous crystal nucleation in droplets is its strong enhancement in the contact mode (when the foreign particle is presumably in some kind of contact with the droplet surface) compared to the immersion mode (particle immersed in the droplet). Heterogeneous centers can have different nucleation thresholds when they act in contact or immersion modes. The underlying physical reasons for this enhancement have remained largely unclear. In this paper we present a model for the thermodynamic enhancement of heterogeneous crystal nucleation in the contact mode compared to the immersion one. To determine if and how the surface of a liquid droplet can thermodynamically stimulate its heterogeneous crystallization, we examine crystal nucleation in the immersion and contact modes by deriving and comparing with each other the reversible works of formation of crystal nuclei in these cases. The line tension of a three-phase contact gives rise to additional terms in the formation free energy of a crystal cluster and affects its Wulff (equilibrium) shape. As an illustration, the proposed model is applied to the heterogeneous nucleation of hexagonal ice crystals on generic macroscopic foreign particles in water droplets at T = 253 K. Our results show that the droplet surface does thermodynamically favor the contact mode over the immersion one. Surprisingly, the numerical evaluations suggest that the line tension contribution (from the contact of three water phases (vapor-liquid-crystal)) to this enhancement may be of the same order of magnitude as or even larger than the surface tension contribution.     

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.     

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.     

New approach to the kinetics of heterogeneous unary nucleation on liquid aerosols of a binary solution

The formation of a droplet on a hygroscopic center may occur either in a barrierless way via Kohler activation or via nucleation by overcoming a free energy barrier. Unlike the former, the latter mechanism of this process has been studied very little and only in the framework of the classical nucleation theory based on the capillarity approximation whereby a nucleating droplet behaves like a bulk liquid. In this paper the authors apply another approach to the kinetics of heterogeneous nucleation on liquid binary aerosols, based on a first passage time analysis which avoids the concept of surface tension for tiny droplets involved in nucleation. Liquid aerosols of a binary solution containing a nonvolatile solute are considered. In addition to modeling aerosols formed through the deliquescence of solid soluble particles, the considered aerosols constitute a rough model of "processed" marine aerosols. The theoretical results are illustrated by numerical calculations for the condensation of water vapor on binary aqueous aerosols with nonvolatile nondissociating solute molecules using Lennard-Jones potentials for the molecular interactions.     

Molecular dynamics study of structural and thermodynamic characteristics of nanodroplets of simple fluid

Equilibrium configurations of Lennard-Jones nanodroplets composed of 10–15000 spherically symmetric molecules placed in the center of a spherical container are studied at constant temperature by the molecular dynamics method. The distribution of local density is found and size dependences of density in the center of droplet, first coordination number, and energy surface tension coinciding for equimolecular dividing surface with specific excess free energy of droplet are studied. Radial distribution function is also determined. It is established that the passage of structural characteristics to their macroscopic values is observed for droplets containing as little as about 300 molecules, while, for energy surface tension, analogous passage for energy surface tension occurs for droplets containing 700–6000 molecules.     

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.     

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.     

Line tension effects for liquid droplets on circular surface domains

We study the morphologies of single liquid droplets wetting a substrate in the presence of the line tension of the three-phase contact line. On a homogeneous substrate, the line tension leads to a discontinuous unbinding of the droplet if its volume is decreased below a critical value. For a droplet wetting a structured surface with a circular domain, a line tension contrast gives rise to discontinuous depinning transitions of the contact line from the domain boundary as the droplet volume is varied. We calculate the corresponding free energy bifurcation diagram analytically for axisymmetric droplet shapes. Numerical minimization of the droplet free energy shows that line tension contrasts can stabilize nonaxisymmetric droplet shapes, thus modifying the bifurcation diagram. These latter shapes should be accessible to experiments and can be used to reveal the presence of a line tension contrast.     

Characterization of watermarks formed in nano-carpet effect

In the spreading of a water droplet on an aligned silicon nanorod array surface, a precursor rim was detected moving ahead of the contact line. In this process, nanorods were bundled by the capillary force to form clusters, and a watermark developed on the surface after water evaporated. The size of the watermark, R(p)max, corresponding to the maximum radius of the precursor rim, followed a simple power law relationship with the volume of water droplet omega, R(p)max proportional to omega(beta). The scaling exponent beta increased when the nanorod height decreased, but all in the vicinity of 1/3. This behavior was attributed to the competition of evaporation and spreading of a water droplet during the spreading process. The size of the bundled nanorod cluster formed by the capillary force not only depended on the nanorod height but also on the location in the watermark. The cluster size almost remained as a constant near the center, and then it decreased with the distance from the center. This phenomenon can be qualitatively interpreted through the change of the total free energy during the precursor invading the nanorod array, by considering the contribution from the mechanical energy change due to the bending and clustering of nanorods.     

Vapor-to-droplet transition in a Lennard-Jones fluid: simulation study of nucleation barriers using the ghost field method

We report a comprehensive Monte Carlo (MC) simulation study of the vapor-to-droplet transition in Lennard-Jones fluid confined to a spherical container with repulsive walls, which is a case study system to investigate homogeneous nucleation. The focus is made on the application of a modified version of the ghost field method (Vishnyakov, A.; Neimark, A. V. J. Chem. Phys. 2003, 119, 9755) to calculate the nucleation barrier. This method allows one to build up a continuous trajectory of equilibrium states stabilized by the ghost field potential, which connects a reference droplet with a reference vapor state. Two computation schemes are employed for free energy calculations, direct thermodynamic integration along the constructed trajectory and umbrella sampling. The nucleation barriers and the size dependence of the surface tension are reported for droplets containing from 260 to 2000 molecules. The MC simulation study is complemented by a review of the simulation methods applied to computing the nucleation barriers and a detailed analysis of the vapor-to-droplet transition by means of the classical nucleation theory.     

Multicomponent nucleation: thermodynamically consistent description of the nucleation work

A thermodynamically consistent formula is derived for the nucleation work in multicomponent homogeneous nucleation. The derivation relies on the conservative dividing surface which defines the nucleus as having specific surface energy equal to the specific surface energy sigma0 of the interface between the macroscopically large new and old phases at coexistence. Expressions are given for the radius of the nucleus defined by the conservative dividing surface and by the surface of tension. As a side result, the curvature dependence of the surface tension sigmaT of the nucleus defined by the surface of tension is also determined. The analysis is valid for nuclei of any size, i.e., for nucleation in the whole range of conditions between the binodal and the spinodal of the metastable old phase provided the inequality sigmaT < or = sigma0 is satisfied. It is found that under the conditions of validity of the analysis the nucleation rate is higher than the nucleation rate given by the classical nucleation theory. The general results are applied to nucleation of unary liquids or solids in binary gaseous, liquid or solid mixtures.     

Surface area controlled heterogeneous nucleation

Heterogeneous nucleation of liquid from a gas phase on nanoparticles has been studied under various saturation ratios and nuclei size. The probability of liquid droplet nucleation, especially at a low degree of deviation from equilibrium, was measured for both atmospheric aerosol particles and engineered nanoparticles Cr(2)O(3). The concept of a critical saturation ratio and the validity of the one-to-one relationship between the nuclei number and the number of droplets were examined. A transient zone between no nucleation and established nucleation termed the surface area controlled nucleation was observed. In this zone, the probability of stable phase formation is determined by the surface area of nuclei. There are two distinctive features of the surface area controlled nucleation: the nucleation probability is much less than 1 and is proportional to the surface area of nuclei. For condensation particle counters (CPCs) counting nanoparticles, these features mean that counts measured are proportional to the surface area of nanoparticles and, therefore, the CPCs counts can be calibrated to measure the surface area.     

Using the Surface Evolver to model droplet formation processes in membrane emulsification

A model was developed to describe the droplet formation mechanism in membrane emulsification from the point of view of Gibbs free energy with the help of the Surface Evolver, which is an interactive finite element program for the study of interfaces shaped by surface tension. A program to test the model was written and run which allows the user to track the droplet shape as it grows, to identify the point of instability due to free energy, and thus predict droplet size. The inputs of the program are pore geometry, oil-aqueous phase interfacial tension, and contact angle. The model reasonably predicted droplet sizes for oblong-shaped pores under quiescent conditions where the force balance approach is not applicable. The model was validated against experimental conditions from the literature where the average error of the predictions compared to the mean droplet sizes was 8%.