Toward a Quantitative Theory of Ultrasmall Liquid Droplets and Vapor—Liquid Nucleation

Microscopic derivation of the thin film equation using the Mori-Zwanzig formalism

The hydrodynamics of thin films is typically described using macroscopic models whose connection to the microscopic particle dynamics is a subject of ongoing research. Existing methods based on density functional theory provide a good description of static thin films but are not sufficient for understanding nonequilibrium dynamics. In this work, we present a microscopic derivation of the thin film equation using the Mori-Zwanzig projection operator formalism. This method allows to directly obtain the correct gradient dynamics structure along with microscopic expressions for mobility and free energy. Our results are verified against molecular dynamics simulations for both simple fluids and polymers.

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.

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

It is argued that to arrive at a quantitative description of the surface tension of a liquid drop as a function of its inverse radius, it is necessary to include the bending rigidity k and Gaussian rigidity k in its description. New formulae for k and k in the context of density functional theory with a non-local, integral expression for the interaction between molecules are presented. These expressions are used to investigate the influence of the choice of Gibbs dividing surface, and it is shown that for a one-component system, the equimolar surface has a special status in the sense that both k and k are then the least sensitive to a change in the location of the dividing surface. Furthermore, the equimolar value for k corresponds to its maximum value and the equimolar value for k corresponds to its minimum value. An explicit evaluation using a short-ranged interaction potential between molecules shows that k is negative with a value around minus 0.5-1.0 kBT and that k is positive with a value that is a bit more than half the magnitude of k. Finally, for dispersion forces between molecules, we show that a term proportional to log(R)/R(2) replaces the rigidity constants and we determine the (universal) proportionality constants.

A perspective on the interfacial properties of nanoscopic liquid drops

The structural and interfacial properties of nanoscopic liquid drops are assessed by means of mechanical, thermodynamical, and statistical mechanical approaches that are discussed in detail, including original developments at both the macroscopic level and the microscopic level of density functional theory (DFT). With a novel analysis we show that a purely macroscopic (static) mechanical treatment can lead to a qualitatively reasonable description of the surface tension and the Tolman length of a liquid drop; the latter parameter, which characterizes the curvature dependence of the tension, is found to be negative and has a magnitude of about a half of the molecular dimension. A mechanical slant cannot, however, be considered satisfactory for small finite-size systems where fluctuation effects are significant. From the opposite perspective, a curvature expansion of the macroscopic thermodynamic properties (density and chemical potential) is then used to demonstrate that a purely thermodynamic approach of this type cannot in itself correctly account for the curvature correction of the surface tension of liquid drops. We emphasize that any approach, e.g., classical nucleation theory, which is based on a purely macroscopic viewpoint, does not lead to a reliable representation when the radius of the drop becomes microscopic. The description of the enhanced inhomogeneity exhibited by small drops (particularly in the dense interior) necessitates a treatment at the molecular level to account for finite-size and surface effects correctly. The so-called mechanical route, which corresponds to a molecular-level extension of the macroscopic theory of elasticity and is particularly popular in molecular dynamics simulation, also appears to be unreliable due to the inherent ambiguity in the definition of the microscopic pressure tensor, an observation which has been known for decades but is frequently ignored. The union of the theory of capillarity (developed in the nineteenth century by Gibbs and then promoted by Tolman) with a microscopic DFT treatment allows for a direct and unambiguous description of the interfacial properties of drops of arbitrary size; DFT provides all of the bulk and surface characteristics of the system that are required to uniquely define its thermodynamic properties. In this vein, we propose a non-local mean-field DFT for Lennard-Jones (LJ) fluids to examine drops of varying size. A comparison of the predictions of our DFT with recent simulation data based on a second-order fluctuation analysis (Sampayo et al 2010 J. Chem. Phys. 132 141101) reveals the consistency of the two treatments. This observation highlights the significance of fluctuation effects in small drops, which give rise to additional entropic (thermal non-mechanical) contributions, in contrast to what one observes in the case of planar interfaces which are governed by the laws of mechanical equilibrium. A small negative Tolman length (which is found to be about a tenth of the molecular diameter) and a non-monotonic behaviour of the surface tension with the drop radius are predicted for the LJ fluid. Finally, the limits of the validity of the Tolman approach, the effect of the range of the intermolecular potential, and the behaviour of bubbles are briefly discussed.

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.

Inkjet metrology: high-accuracy mass measurements of microdroplets produced by a drop-on-demand dispenser

We describe gravimetric methods for measuring the mass of droplets generated by a drop-on-demand (DOD) microdispenser. Droplets are deposited, either continuously at a known frequency or as a burst of known number, into a cylinder positioned on a submicrogram balance. Mass measurements are acquired precisely by computer, and results are corrected for evaporation. Capabilities are demonstrated using isobutyl alcohol droplets. For ejection rates greater than 100 Hz, the repeatability of droplet mass measurements was 0.2%, while the combined relative standard uncertainty (u(c)) was 0.9%. When bursts of droplets were dispensed, the limit of quantitation was 72 microg (1490 droplets) with u(c) = 1.0%. Individual droplet size in a burst was evaluated by high-speed videography. Diameters were consistent from the tenth droplet onward, and the mass of an individual droplet was best estimated by the average droplet mass with a combined uncertainty of about 1%. Diameters of the first several droplets were anomalous, but their contribution was accounted for when dispensing bursts. Above the limits of quantitation, the gravimetric methods provided statistically equivalent results and permit detailed study of operational factors that influence droplet mass during dispensing, including the development of reliable microassays and standard materials using DOD technologies.

Cluster formation and bulk phase behavior of colloidal dispersions

We investigated cluster formation in model colloids where the interparticle potential beyond the collision diameter includes a short-ranged attraction and a longer-ranged repulsion. The structure and thermodynamic properties of the model system can be accurately described by the first-order mean-spherical approximation (FMSA) and a perturbation theory for bulk phases and by a nonlocal density-functional theory (DFT) for stable/metastable clusters. In corroboration with recent experiments and molecular simulations, we predicted that the bulk to cluster transition may occur at both one-phase and two-phase regions of the colloidal phase diagram. While formation of colloidal clusters appears uncorrelated with the bulk phase transitions, the local phenomena may nevertheless play an important role in the dynamics of stable colloids and the kinetics of colloidal phase transitions.

Nucleation and hysteresis of vapor-liquid phase transitions in confined spaces: effects of fluid-wall interaction

In this work, we propose a method to stabilize a nucleus in the framework of lattice density-functional theory (LDFT) by imposing a suitable constraint. Using this method, the shape of critical nucleus and height of the nucleation barrier can be determined without using a predefined nucleus as input. As an application of this method, we study the nucleation behavior of vapor-liquid transition in nanosquare pores with infinite length and relate the observed hysteresis loop on an adsorption isotherm to the nucleation mechanism. According to the dependence of hysteresis and the nucleation mechanism on the fluid-wall interaction, w , in this work, we have classified w into three regions ( w>0.9 , 0.1< or =w< or =0.9 , and w<0.1 ), which are denoted as strongly, moderately, and weakly attractive fluid-wall interaction, respectively. The dependence of hysteresis on the fluid-wall interaction is interpreted by the different nucleation mechanisms. Our constrained LDFT calculations also show that the different transition paths may induce different nucleation behaviors. The transition path dependence should be considered if morphological transition of nuclei exists during a nucleation process.