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.     

Two interacting particles in a spherical pore

In this work we analytically evaluate, for the first time, the exact canonical partition function for two interacting spherical particles into a spherical pore. The interaction with the spherical substrate and between particles is described by an attractive square-well and a square-shoulder potential. In addition, we obtain exact expressions for both the one particle and an averaged two particle density distribution. We develop a thermodynamic approach to few-body systems by introducing a method based on thermodynamic measures [I. Urrutia, J. Chem. Phys. 134, 104503 (2010)] for nonhard interaction potentials. This analysis enables us to obtain expressions for the pressure, the surface tension, and the equivalent magnitudes for the total and Gaussian curvatures. As a by-product, we solve systems composed of two particles outside a fixed spherical obstacle. We study the low density limit for a many-body system confined to a spherical cavity and a many-body system surrounding a spherical obstacle. From this analysis we derive the exact first order dependence of the surface tension and Tolman length. Our findings show that the Tolman length goes to zero in the case of a purely hard wall spherical substrate, but contains a zero order term in density for square-well and square-shoulder wall-fluid potentials. This suggests that any nonhard wall-fluid potential should produce a non-null zero order term in the Tolman length.     

On the determination of the structure and tension of the interface between a fluid and a curved hard wall

The structure and tension of the interface between a fluid and a spherically shaped hard wall are studied theoretically. The authors show the equivalence of different expressions for the surface tension and Tolman length using the squared-gradient model and density functional theory with a nonlocal, integral expression for the interaction between molecules. Even though both these models yield equilibrium density profiles that do not satisfy the wall theorem, they still obey the basic requirement of mechanical equilibrium. The authors trace back the origin of the difference between these two observations to the (lack of) continuity of the cavity function at the hard wall.     

An exact formalism to study the thermodynamic properties of hard-sphere systems under spherical confinement

This paper presents a modified grand canonical ensemble which provides a new simple and efficient scheme to study few-body fluid-like inhomogeneous systems under confinement. The new formalism is implemented to investigate the exact thermodynamic properties of a hard sphere (HS) fluid-like system with up to three particles confined in a spherical cavity. In addition, the partition function of this system was used to analyze the surface thermodynamic properties of the many-HS system and to derive the exact curvature dependence of both the surface tension and adsorption in powers of the density. The expressions for the surface tension and the adsorption were also obtained for the many-HS system outside of a fixed hard spherical object. We used these results to derive the dependence of the fluid-substrate Tolman length up to first order in density.     

Scaled particle theory for hard sphere pairs. II. Numerical analysis

We use the extension of scaled particle theory presented in the accompanying paper [F. H. Stillinger et al., J. Chem. Phys. 125, 204504 (2006)] to calculate numerically the pair correlation function of the hard sphere fluid over the density range 0< or =rhosigma(3)< or =0.96. Comparison with computer simulation results reveals that the new theory is able to capture accurately the fluid's structure across the entire density range examined. The pressure predicted via the virial route is systematically lower than simulation results, while that obtained using the compressibility route is lower than simulation predictions for rhosigma(3)< or =0.67 and higher than simulation predictions for rhosigma(3)> or =0.67. Numerical predictions are also presented for the surface tension and Tolman length of the hard sphere fluid.     

Small solid particles and liquid lenses at fluid/fluid interfaces

The behaviour of small solid particles and liquid droplets at fluid interfaces is of wide interest, in part because of the roles they play in the stability of foams and emulsions. Here we focus on solid particles at liquid interfaces, both singly and in highly structured monolayers. We briefly mention small oil lenses on water in connection with the determination of line tension, τ. Particles are surface-active in the sense that they often adhere quite strongly to liquid surfaces, although of course they are not usually amphiphilic. The three-phase contact line around a particle at an interface is associated with an excess free energy resulting in a tendency of the line to contract (positive τ, which is a 1D analogue of surface tension) or to expand (negative τ). Positive line tension acts so as to push the contact angle of a particle with the fluid interface further away from 90°, i.e. to force the particle towards the more “wetting” of the two bulk phases. It also leads to activation barriers to entry and departure of particles from an interface. The behaviour of particle monolayers at octane/water interfaces is also discussed . It is found that, for monodisperse spherical polystyrene particles containing ionisable sulphate groups at the surface, highly ordered monolayers are formed. This appears to result from very long range electrostatic repulsion mediated through the oil phase. Surface pressure–surface area isotherms are discussed for particle monolayers and it is shown, using light microscopy, that at monolayer “collapse” particles are not expelled from the monolayers but rather the monolayer folds, remaining intact. This has an important bearing on methods, involving the use of the Langmuir trough, for the experimental determination of contact angles and line tensions in particulate systems. Received: 18 July 1999/Accepted: 30 August 1999     

Some estimates of the surface tension of curved surfaces using density functional theory

Density functional theory is used to calculate the surface tension of planar and slightly curved surfaces, which can be written as gamma(R)=gamma(infinity)(1-2delta(infinity)R), where R is the radius of curvature of the surface. Calculations are performed for a Lennard-Jones fluid, split into a hard-sphere repulsive potential and an attractive part. The repulsive part is treated using the local density approximation. The attractive part is treated using a high temperature approximation (HTA) in which the pair correlation function is approximated by the Percus-Yevick pair correlation function of a uniform hard-sphere fluid evaluated at a position-dependent average density. An expression relating the Tolman length delta(infinity) to the density profile of the planar surface is derived. Numerical results are presented for the planar surface tension gamma(infinity) and for delta(infinity) and are compared with those using mean field theory (MFT) and with those using the square-gradient approximation. Values for gamma(infinity) using the HTA are 30%-40% higher than those using MFT. Values for delta(infinity) using the HTA are around -0.1 (in units of the Lennard-Jones parameter sigma) and only weakly dependent on temperature. These values are less negative than the values from MFT. The square-gradient approximation gives reasonable estimates of the more accurate nonlocal results for both the MFT and the HTA.     

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) → ∞).     

Viscoelastic theory for nematic polymer interfaces

A macroscopic theory for the dynamics of compressible nematic polymer‐viscous fluid interfaces is developed from first principles. The theory is used to define and characterize the basic interfacial viscoelastic material properties of the ordered interfaces. The theory is based on a decomposition of the kinematic fields and nematic tensor order parameter that takes into account the symmetry breaking of the interface. The interfacial rate of entropy production used to identify the interfacial viscoelastic modes is given in terms of surface rate of deformation tensor and the surface Jaumann derivative of the tangential component nematic tensor order parameter. The derived surface viscous stress tensor is asymmetric and thus describes surface flow‐induced changes in the tensor order parameter. Consistency with the Boussinesq surface fluid appropriate for Newtonian interfaces is established. The interfacial material functions are identified as the dynamic surface tension, the interfacial dilational viscosities, and the interfacial shear viscosities. The interfacial material functions depend on the surface tensor order parameter and as a consequence anisotropy is their characteristic feature. Two characteristic interfacial tensions and two dilational viscosities are predicted depending on the director orientation. In addition six interfacial shear viscosities arise as the directors sample the velocity, velocity gradient, and vorticity directions. Finally the theory provides for the necessary theoretical tools needed to describe the interfacial dynamics of nematic polymer interfaces, such as capillary instabilities, Marangoni flows, and wetting phenomena.     

Surface activity of lysozyme and dipalmitoyl phosphatidylcholine vesicles at compressed and supercritical fluid interfaces

The surface activities of lysozyme and dipalmitoyl phosphatidylcholine (DPPC) vesicles at aqueous/compressed fluid interfaces are examined via high-pressure interfacial tension measurements using the pendant drop technique. The density and interfacial tension in compressible fluid systems vary significantly with pressure, providing a versatile medium for elucidating interactions between biomolecules and fluid interfaces and a method to elicit pressure-dependent interfacial morphological responses. The effects of lysozyme concentration (0.0008, 0.01, and 1 mg/mL) and pressure (> or = 7 MPa) on the dynamic surface response in the presence of ethane, propane, N2, and CO2 at 298 K were examined. Interfacial lysozyme adsorption reduced the induction phase and quickly led to interfacial tensions consistent with protein conformational changes and monolayer saturation at the compressed fluid interfaces. Protein adsorption, as indicated by surface pressure, correlated with calculated Hamaker constants for the compressed gases, denoting the importance of dispersion interactions. For DPPC at aqueous/compressed or aqueous/supercritical CO2 interfaces (1.8-20.7 MPa, 308 K), 2-3-fold reductions in interfacial tension were observed relative to the pure binary fluid system. The resulting surface pressures infer pressure-dependent morphological changes within the DPPC monolayer.     

Hydrophobic solvation of Gay-Berne particles in modified water models

The solvation of large hydrophobic solutes, modeled as repulsive and attractive Gay-Berne oblate ellipsoids, is characterized in several modified water liquids using the SPC/E model as the reference water fluid. We find that small amounts of attraction between the Gay-Berne particle and any model fluid result in wetting of the hydrophobic surface. However, significant differences are found among the modified and SPC/E water models and the critical distances in which they dewet the hydrophobic surfaces of pairs of repulsive Gay-Berne particles. We find that the dewetting trends for repulsive Gay-Berne particles in the various model liquids correlate directly with their surface tensions, the widths of the interfaces they form, and the openness of their network structure. The largest critical separations are found in liquids with the smallest surface tensions and the broadest interfaces as measured by the Egelstaff-Widom length.     

The liquid–gas interface of oxygen–nitrogen solutions 2: Description in the Framework of the van der Waals gradient theory

The van der Waals gradient theory (vdW GT) is used to calculate surface tension, density profiles, adsorption, the Tolman length and to determine the position of dividing surfaces in the liquid–gas interface of an oxygen–nitrogen solution. The Helmholtz energy density (HED) is determined via an equation of state (EOS), unified for a liquid and gas, which describes stable, metastable and two-phase states of solutions. The influence parameters are calculated from data on the surface tension of pure components with the use of the mixing rule. At temperatures T > 100 K the vdW GT describes experimental data on the surface tension of oxygen–nitrogen solutions [V.G. Baidakov, A.M. Kaverin, V.N. Andbaeva, The liquid–gas interface of oxygen–nitrogen solutions: 1. Surface tension, Fluid Phase Equilib. 270 (2008) 116–120] within the experimental error. It is shown that the Tolman length, which determines the dependence of surface tension on the curvature of the dividing surface, depends considerably on the solution concentration.     

Mechanical model for anisotropic curved interfaces with applications to surfactant-laden liquid-liquid crystal interfaces

A mechanical model for anisotropic curved interfaces, applicable to thermodynamically closed surfactant-laden liquid-liquid crystal interfaces is developed. The model takes into account the mechanical effects due to surface bending and surface tilting (anchoring) and incorporates liquid crystal anisotropy into classical fluid membrane mechanics. In the absence of the aligned liquid crystal, the model converges to the fluid membrane mechanical model, and in the absence of surfactant, it converges to the nematic interface mechanical model. Use of the well-known Helfrich-Rapini-Papoular surface energies leads to the Laplace equation for anisotropic curved interfaces, whose material limits are the vesicle shape equation and the liquid crystal Herring equation. Applications of the model to shape selection in liquid drops embedded in aligned nematic liquid crystals illustrates the competition between anchoring and bending and shows how anisotropic surface tension distorts the droplet and how bending tends to restore the spherical shape. This theory presented in this article shows that the interaction of interfacial anchoring and bending creates new regimes in classical fluid membrane mechanics.     

On the curvature dependence of the interfacial tension in a symmetric three-component interface

We consider a symmetric interface between two polymers A(N) and B(N) in a common monomeric solvent S using the mean-field Scheutjens-Fleer self-consistent field theory and focus on the curvature dependence of the interfacial tension. In multi-component systems there is not one unique scenario to curve such an interface. We elaborate on this by keeping either the chemical potential of the solvent or the bulk concentration of the solvent fixed, that is we focus on the semi-grand canonical ensemble case. Following Helfrich, we expand the surface tension as a Taylor series in the curvature parameters and find that there is a non-zero linear dependence of the interfacial tension on the mean curvature in both cases. This implies a finite Tolman length. In a thermodynamic analysis we prove that the non-zero Tolman length is related to the adsorption of solvent at the interface. Similar, but not the same, correlations between the solvent adsorption and the Tolman length are found in the two scenarios. This result indicates that one should be careful with symmetry arguments in a Helfrich analysis, in particular for systems that have a finite interfacial tension: one not only should consider the structural symmetry of the interface, but also consider the constraints that are enforced upon imposing the curvature. The volume fraction of solvent, the chain length N as well as the interaction parameter chi(AB) in the system can be used to take the system in the direction of the critical point. The usual critical behavior is found. Both the width of the interface and the Tolman length diverge, whereas the density difference between the two phases, adsorbed amount of solvent at the interface, interfacial tension, spontaneous curvature, mean bending modulus as well as the Gaussian bending modulus vanish upon approach of the critical point.     

Surface tensions in NaCl-water-air systems from MD simulations

Surface tensions for liquid-vapor (lv), solid-liquid (sl), and solid-vapor (sv) interfaces are calculated from molecular dynamics simulations of the NaCl-water-air system. Three distinct calculation techniques based on thermodynamic properties are used to describe the multicomponent mixtures. Simulations of each bulk phase (including a liquid saturated solution) and various interfaces are carried out at both NPT and NVT conditions. The thermodynamic relation for energy difference between interface and bulk phases provides an upper bound to the surface tension, while the energy-integral and test area methods provide direct estimates. At 1 atm and 300 K, the best predictions for surface tensions are sigmasv (NaCl-air) of 114 mN m(-1), sigmasl (NaCl- soln) of 63 mN m(-1), sigmalv (soln-air) of 82 mN m(-1), and sigmalv (water-air) of 66 mN m(-1). The calculated surface tensions from simulations have uncertainties between 5 and 10%, which are higher than measurements for the liquid interfaces and lower than the measurement uncertainty for the solid interfaces. The calculated upper bounds for surface tensions of liquid interfaces compare well with experimental results but provide no improvement over existing measurements. However, the bounding values for solid interfaces lower uncertainty by as much as a factor of 10 as compared to the indirect experimental measurements currently available. The energy-integral and test area methods appear to underestimate the surface tension of water by 10%, which is consistent with previous studies using similar model potentials. The calculated upper bounds of surface tension show a weakly positive correlation with pressure in the 0.1-100 atm range for liquid-solid, liquid-vapor, and solid-vapor interfaces.     

Orientation of a nanocylinder at a fluid interface

A nanocylinder placed on a fluid interface can assume an end-on or side-on orientation, or it can immerse itself in the surrounding bulk phases. Any of these orientations can satisfy a mechanical force balance when the particle is small enough that gravitational effects are negligible. The orientation is determined by the surface energies of the fluid-solid, fluid-vapor, and vapor-solid surfaces. A comparison of the energy of each state allows phase diagrams to be defined in terms of the scaled aspect ratio x=2L/pir and the contact angle thetao, where L and r denote the nanocylinder length and radius, respectively. Line tension can also influence the orientations by changing the equilibrium contact angle theta and by increasing the energetic cost of the contact line. Phase diagrams accounting for positive line tensions Sigma are also constructed. These phase diagrams can be divided into two classes. In the first, over some range of x and Sigma, nanocylinders can be driven from side-on to end-on orientations with increasing Sigma. This transition terminates at a triple point where the side-on, end-on, and immersed energies are the same. In the second class, there is no triple point and, for a range of Sigma values, nanocylinders of all aspect ratios x prefer an end-on orientation. In all cases, for high enough Sigma, line tension drives a wetting transition similar to that already noted in the literature for spherical particles. The zero line tension predictions are compared favorably to experiment, in which functionalized gold nanowires made by template synthesis are spread at aqueous-gas interfaces, immobilized using a gel-fixation technique, and observed by scanning electron microscopy. The small aspect ratio particles (disks) were in an end-on configuration, while the longer nanowires were in a side-on orientation, in agreement with the theory.     

Homogeneous nucleation of water between 200 and 240 K: new wave tube data and estimation of the Tolman length

We have measured homogeneous nucleation rates of water at 200-240 K in the carrier gas helium, in the range of 10(13) - 10(17) m(-3) s(-1) using an expansion wave tube. The rates agree well with the results of Wolk and Strey [J. Phys. Chem. B 105, 11683 (2001)] in the range of overlap (220-240 K), and are summarized by the empirical fit J = S exp[4.6 + 0.244T-(906.8 - 2.914T)(ln S)(2)], with J the nucleation rate in m(-3) s(-1), S the supersaturation, and T the temperature in K. We find that the supersaturation dependence of both our rates and those of Wolk and Strey is lower than classical theory predicts, and that the critical cluster is smaller than the classical critical size. These deviations are explained in the framework of the Tolman theory for surface tension, and the "Tolman length" is estimated from our experimental results. We find a positive Tolman length that increases with decreasing temperature, from about 0.1 Angstrom at 260 K to (0.6 +/- 0.4) Angstroms at 200 K. We present a nucleation rate expression that takes the Tolman length into account and show that both the supersaturation and temperature dependence are improved, compared to the classical theory.     

Evaluation of surface tension and Tolman length as a function of droplet radius from experimental nucleation rate and supersaturation ratio: metal vapor homogeneous nucleation

Zinc and silver vapor homogeneous nucleations are studied experimentally at the temperature from 600 to 725 and 870 K, respectively, in a laminar flow diffusion chamber with Ar as a carrier gas at atmospheric pressure. The size, shape, and concentration of aerosol particles outcoming the diffusion chamber are analyzed by a transmission electron microscope and an automatic diffusion battery. The wall deposit is studied by a scanning electron microscope (SEM). Using SEM data the nucleation rate for both Zn and Ag is estimated as 10(10) cm(-3) s(-1). The dependence of critical supersaturation on temperature for Zn and Ag measured in this paper as well as Li, Na, Cs, Ag, Mg, and Hg measured elsewhere is analyzed. To this aim the classical nucleation theory is extended by the dependence of surface tension on the nucleus radius. The preexponent in the formula for the vapor nucleation rate is derived using the formula for the work of formation of noncritical embryo [obtained by Nishioka and Kusaka [J. Chem. Phys. 96, 5370 (1992)] and later by Debenedetti and Reiss [J. Chem. Phys. 108, 5498 (1998)]] and Reiss replacement factor. Using this preexponent and the Gibbs formula for the work of formation of critical nucleus the dependence of surface tension on the radius R(S) of the surface of tension is evaluated from the nucleation data for above-mentioned metals. For the alkali metals and Ag the surface tension was determined to be a strong function of R(S). For the bivalent metals (Zn, Hg, and Mg) the surface tension was independent of radius in the experimental range. A new formula for the Tolman length delta as a function of surface tension and radius R(S) is derived by integration of Gibbs-Tolman-Koenig equation assuming that delta is a monotonic function of radius. The formula derived is more correct than the Tolman formula and convenient for the elaboration of experimental data. Using this formula the values of delta are determined as a function of R(S) from the experimental nucleation data. It is determined that all the metals considered are characterized by strong dependence of delta on radius; for the bivalent metals delta changes sign.     

A molecular dynamics study to determine the solid-liquid interfacial tension using test area simulation method (TASM)

Molecular dynamics (MD) studies on heat transfer from a heated nanoparticle into the surrounding fluid have indicated that the fluid next to a spherical nanoparticle can get heated well above its boiling point without observing a phase change, while a contradicting behavior was observed for a flat surface-fluid interface. Another interesting observation is that the critical heat flux was found to increase with increase in the wetting characteristics of solid. Thus, the interfacial tension or free energy of solid-liquid interface could play a pivotal role in the mechanism of heat transfer. A recent study by Gloor et al. [J. Chem. Phys. 123, 134703 (2005)] has proposed test area simulation method (TASM) for the determination of interfacial tension. The present study involves the determination and the comparison of solid-liquid interfacial tension for planar and spherical interfaces using MD based on TASM and analyze the results. A higher interfacial tension value is observed for spherical nanoparticle fluid interface compared to flat surface fluid interface. The results also indicate that the solid-liquid interfacial tension is a size and temperature dependent property. The results from this study are also expected to give better insights into the possible reasons for the observed differences in the thermal transport for spherical nanoparticle-liquid interface compared to planar-liquid interface.