Density functional theory of a curved liquid-vapour interface: evaluation of the rigidity constants

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.

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.

Mesoscale perspective on the Tolman length

We demonstrate that the multiphase Shan-Chen lattice Boltzmann method (LBM) yields a curvature dependent surface tension σ as computed from three-dimensional hydrostatic droplets and bubbles simulations. Such curvature dependence is routinely characterized, at first order, by the so-called Tolman length δ. LBM allows one to precisely compute σ at the surface of tension R_{s} and determine the Tolman length from the coefficient of the first order correction. The corresponding values of δ display universality for different equations of state, following a power-law scaling near the critical temperature. The Tolman length has been studied so far mainly via computationally demanding Molecular Dynamics simulations or by means of Density Functional Theory approaches playing a pivotal role in extending Classical Nucleation Theory. The present results open a hydrodynamic-compliant mesoscale arena, in which the fundamental role of the Tolman length, alongside real-world applications to cavitation phenomena, can be effectively tackled. All the results can be independently reproduced through the "idea.deploy" framework.

Curvature Corrections Remove the Inconsistencies of Binary Classical Nucleation Theory

The study of nucleation in fluid mixtures exposes challenges beyond those of pure systems. A striking example is homogeneous condensation in highly surface-active water-alcohol mixtures, where classical nucleation theory yields an unphysical, negative number of water molecules in the critical embryo. This flaw has rendered multicomponent nucleation theory useless for many industrial and scientific applications. Here, we show that this inconsistency is removed by properly incorporating the curvature dependence of the surface tension of the mixture into classical nucleation theory for multicomponent systems. The Gibbs adsorption equation is used to explain the origin of the inconsistency by linking the molecules adsorbed at the interface to the curvature corrections of the surface tension. The Tolman length and rigidity constant are determined for several water-alcohol mixtures and used to show that the corrected theory is free of physical inconsistencies and provides accurate predictions of the nucleation rates. In particular, for the ethanol-water and propanol-water mixtures, the average error in the predicted nucleation rates is reduced from 11-15 orders of magnitude to below 1.5. The curvature-corrected nucleation theory opens the door to reliable predictions of nucleation rates in multicomponent systems, which are crucial for applications ranging from atmospheric science to research on volcanos.

Studying vapor-liquid transition using a generalized ensemble

Homogeneous vapor-liquid nucleation is studied using the generalized Replica Exchange Method (gREM). The generalized ensemble allows the study of unstable states that cannot directly be studied in the canonical ensemble. Along with replica exchange, this allows for efficient sampling of the multiple states in a single simulation. Statistical Temperature Weighted Histogram Analysis Method is used for postprocessing to get a continuous free energy curve from bulk vapor to bulk liquid. gREM allows the study of planar, cylindrical, and spherical interfaces in a single simulation. The excess Gibbs free energy for the formation of a spherical liquid droplet in vapor for a Lennard-Jones system is calculated from the free energy curve and compared against the umbrella sampling results. The nucleation free energy barrier obtained from gREM is then used to calculate the nucleation rate without relying on any classification scheme for separating the vapor and liquid.

Measurement of the Curvature-Dependent Surface Tension in Nucleating Colloidal Liquids

The curvature dependence of the surface tension is central to the nucleation of liquids, but remains difficult to access experimentally and predict theoretically. This curvature dependence arises from the curvature-dependent molecular structure, which, for small nuclei, can deviate significantly from that of the planar liquid interface. Simulations and density functional theory have been used to predict this curvature dependence, however with contradicting results. Here, we provide the first direct measurement of the curvature-dependent surface tension in nucleating colloidal liquids. We employ critical Casimir forces to finely adjust colloidal particle interactions and induce liquid nucleation, and image individual nuclei at the particle scale to measure their curvature-dependent surface tension directly from thermally excited surface distortions. Using continuum models, we elucidate the interplay between nucleus structure, particle pair potential, and surface tension. Our results reveal a 20% lower surface tension for nuclei of critical size compared to bulk liquids, leading to 3 orders of magnitude higher nucleation rates, thus highlighting the importance of surface tension curvature corrections for accurate prediction of nucleation rates.

Tolman lengths and rigidity constants of multicomponent fluids: Fundamental theory and numerical examples

The curvature dependence of the surface tension can be described by the Tolman length (first-order correction) and the rigidity constants (second-order corrections) through the Helfrich expansion. We present and explain the general theory for this dependence for multicomponent fluids and calculate the Tolman length and rigidity constants for a hexane-heptane mixture by use of square gradient theory. We show that the Tolman length of multicomponent fluids is independent of the choice of dividing surface and present simple formulae that capture the change in the rigidity constants for different choices of dividing surface. For multicomponent fluids, the Tolman length, the rigidity constants, and the accuracy of the Helfrich expansion depend on the choice of path in composition and pressure space along which droplets and bubbles are considered. For the hexane-heptane mixture, we find that the most accurate choice of path is the direction of constant liquid-phase composition. For this path, the Tolman length and rigidity constants are nearly linear in the mole fraction of the liquid phase, and the Helfrich expansion represents the surface tension of hexane-heptane droplets and bubbles within 0.1% down to radii of 3 nm. The presented framework is applicable to a wide range of fluid mixtures and can be used to accurately represent the surface tension of nanoscopic bubbles and droplets.

Temperature anisotropy at equilibrium reveals nonlocal entropic contributions to interfacial properties

Density gradient theory for fluids has played a key role in the study of interfacial phenomena for a century. In this work, we revisit its fundamentals by examining the vapor-liquid interface of argon, represented by the cut and shifted Lennard-Jones fluid. The starting point has traditionally been a Helmholtz energy functional using mass densities as arguments. By using rather the internal energy as starting point and including the entropy density as an additional argument, following thereby the phenomenological approach from classical thermodynamics, the extended theory suggests that the configurational part of the temperature has different contributions from the parallel and perpendicular directions at the interface, even at equilibrium. We find a similar anisotropy by examining the configurational temperature in molecular dynamics simulations and obtain a qualitative agreement between theory and simulations. The extended theory shows that the temperature anisotropy originates in nonlocal entropic contributions, which are currently missing from the classical theory. The nonlocal entropic contributions discussed in this work are likely to play a role in the description of both equilibrium and nonequilibrium properties of interfaces. At equilibrium, they influence the temperature- and curvature-dependence of the surface tension. Across the vapor-liquid interface of the Lennard Jones fluid, we find that the maximum in the temperature anisotropy coincides precisely with the maximum in the thermal resistivity relative to the equimolar surface, where the integral of the thermal resistivity gives the Kapitza resistance. This links the temperature anisotropy at equilibrium to the Kapitza resistance of the vapor-liquid interface at nonequilibrium.

Bending and Gaussian rigidities of confined soft spheres from second-order virial series

We use virial series to study the equilibrium properties of confined soft-spheres fluids interacting through the inverse-power potentials. The confinement is induced by hard walls with planar, spherical, and cylindrical shapes. We evaluate analytically the coefficients of order two in density of the wall-fluid surface tension γ and analyze the curvature contributions to the free energy. Emphasis is in bending and Gaussian rigidities, which are found analytically at order two in density. Their contribution to γ(R) and the accuracy of different truncation procedures to the low curvature expansion are discussed. Finally, several universal relations that apply to low-density fluids are analyzed.

Reply to "Comment on 'Simple improvements to classical bubble nucleation models' "

We reply to the Comment by Schmelzer and Baidakov [Phys. Rev. E 94, 026801 (2016)].10.1103/PhysRevE.94.026801 They suggest that a more modern approach than the classic description by Tolman is necessary to model the surface tension of curved interfaces. Therefore we now consider the higher-order Helfrich correction, rather than the simpler first-order Tolman correction. Using a recent parametrization of the Helfrich correction provided by Wilhelmsen et al. [J. Chem. Phys. 142, 064706 (2015)]JCPSA60021-960610.1063/1.4907588, we test this description against measurements from our simulations, and find an agreement stronger than what the pure Tolman description offers. Our analyses suggest a necessary correction of order higher than the second for small bubbles with radius ≲1 nm. In addition, we respond to other minor criticism about our results.

Coherent description of transport across the water interface: From nanodroplets to climate models

Transport of mass and energy across the vapor-liquid interface of water is of central importance in a variety of contexts such as climate models, weather forecasts, and power plants. We provide a complete description of the transport properties of the vapor-liquid interface of water with the framework of nonequilibrium thermodynamics. Transport across the planar interface is then described by 3 interface transfer coefficients where 9 more coefficients extend the description to curved interfaces. We obtain all coefficients in the range 260-560 K by taking advantage of water evaporation experiments at low temperatures, nonequilibrium molecular dynamics with the TIP4P/2005 rigid-water-molecule model at high temperatures, and square gradient theory to represent the whole range. Square gradient theory is used to link the region where experiments are possible (low vapor pressures) to the region where nonequilibrium molecular dynamics can be done (high vapor pressures). This enables a description of transport across the planar water interface, interfaces of bubbles, and droplets, as well as interfaces of water structures with complex geometries. The results are likely to improve the description of evaporation and condensation of water at widely different scales; they open a route to improve the understanding of nanodroplets on a small scale and the precision of climate models on a large scale.

Curvature Dependence of the Liquid-Vapor Surface Tension beyond the Tolman Approximation

Surface tension is a macroscopic manifestation of the cohesion of matter, and its value σ_{∞} is readily measured for a flat liquid-vapor interface. For interfaces with a small radius of curvature R, the surface tension might differ from σ_{∞}. The Tolman equation, σ(R)=σ_{∞}/(1+2δ/R), with δ a constant length, is commonly used to describe nanoscale phenomena such as nucleation. Here we report experiments on nucleation of bubbles in ethanol and n-heptane, and their analysis in combination with their counterparts for the nucleation of droplets in supersaturated vapors, and with water data. We show that neither a constant surface tension nor the Tolman equation can consistently describe the data. We also investigate a model including 1/R and 1/R^{2} terms in σ(R). We describe a general procedure to obtain the coefficients of these terms from detailed nucleation experiments. This work explains the conflicting values obtained for the Tolman length in previous analyses, and suggests directions for future work.

Simple improvements to classical bubble nucleation models

We revisit classical nucleation theory (CNT) for the homogeneous bubble nucleation rate and improve the classical formula using a correct prefactor in the nucleation rate. Most of the previous theoretical studies have used the constant prefactor determined by the bubble growth due to the evaporation process from the bubble surface. However, the growth of bubbles is also regulated by the thermal conduction, the viscosity, and the inertia of liquid motion. These effects can decrease the prefactor significantly, especially when the liquid pressure is much smaller than the equilibrium one. The deviation in the nucleation rate between the improved formula and the CNT can be as large as several orders of magnitude. Our improved, accurate prefactor and recent advances in molecular dynamics simulations and laboratory experiments for argon bubble nucleation enable us to precisely constrain the free energy barrier for bubble nucleation. Assuming the correction to the CNT free energy is of the functional form suggested by Tolman, the precise evaluations of the free energy barriers suggest the Tolman length is ≃0.3σ independently of the temperature for argon bubble nucleation, where σ is the unit length of the Lennard-Jones potential. With this Tolman correction and our prefactor one gets accurate bubble nucleation rate predictions in the parameter range probed by current experiments and molecular dynamics simulations.