Pair correlation functions and the wavevector-dependent surface tension in a simple density functional treatment of the liquid-vapour interface

Structure of liquid-vapor interfaces: Perspectives from liquid state theory, large-scale simulations, and potential grazing-incidence x-ray diffraction

Grazing-incidence x-ray diffraction (GIXRD) is a scattering technique that allows one to characterize the structure of fluid interfaces down to the molecular scale, including the measurement of surface tension and interface roughness. However, the corresponding standard data analysis at nonzero wave numbers has been criticized as to be inconclusive because the scattering intensity is polluted by the unavoidable scattering from the bulk. Here, we overcome this ambiguity by proposing a physically consistent model of the bulk contribution based on a minimal set of assumptions of experimental relevance. To this end, we derive an explicit integral expression for the background scattering, which can be determined numerically from the static structure factors of the coexisting bulk phases as independent input. Concerning the interpretation of GIXRD data inferred from computer simulations, we extend the model to account also for the finite sizes of the bulk phases, which are unavoidable in simulations. The corresponding leading-order correction beyond the dominant contribution to the scattered intensity is revealed by asymptotic analysis, which is characterized by the competition between the linear system size and the x-ray penetration depth in the case of simulations. Specifically, we have calculated the expected GIXRD intensity for scattering at the planar liquid-vapor interface of Lennard-Jones fluids with truncated pair interactions via extensive, high-precision computer simulations. The reported data cover interfacial and bulk properties of fluid states along the whole liquid-vapor coexistence line. A sensitivity analysis shows that our findings are robust with respect to the detailed definition of the mean interface position. We conclude that previous claims of an enhanced surface tension at mesoscopic scales are amenable to unambiguous tests via scattering experiments.

Finite-size corrections for the static structure factor of a liquid slab with open boundaries

The presence of a confining boundary can modify the local structure of a liquid markedly. In addition, small samples of finite size are known to exhibit systematic deviations of thermodynamic quantities relative to their bulk values. Here, we consider the static structure factor of a liquid sample in slab geometry with open boundaries at the surfaces, which can be thought of as virtually cutting out the sample from a macroscopically large, homogeneous fluid. This situation is a relevant limit for the interpretation of grazing-incidence diffraction experiments at liquid interfaces and films. We derive an exact, closed expression for the slab structure factor, with the bulk structure factor as the only input. This shows that such free boundary conditions cause significant differences between the two structure factors, in particular, at small wavenumbers. An asymptotic analysis of this result yields the scaling exponent and an accurate, useful approximation of these finite-size corrections. Furthermore, the open boundaries permit the interpretation of the slab as an open system, supporting particle exchange with a reservoir. We relate the slab structure factor to the particle number fluctuations and discuss conditions under which the subvolume of the slab represents a grand canonical ensemble with chemical potential μ and temperature T. Thus, the open slab serves as a test-bed for the small-system thermodynamics in a μT reservoir. We provide a microscopically justified and exact result for the size dependence of the isothermal compressibility. Our findings are corroborated by simulation data for Lennard-Jones liquids at two representative temperatures.

Microscopic determination of correlations in the fluid interfacial region in the presence of liquid-gas asymmetry

In a recent article, we showed how the properties of the density-density correlation function and its integral, the local structure factor, in the fluid interfacial region, in systems with short-ranged forces, can be understood microscopically by considering the resonances of the local structure factor [A. O. Parry and C. Rascón, Nat. Phys. 15, 287 (2019)NPAHAX1745-247310.1038/s41567-018-0361-z]. Here, we illustrate, using mean-field square-gradient theory and the more microscopic Sullivan density functional model, and how this approach generalizes when there is liquid-gas asymmetry, i.e., when the bulk correlation lengths of the coexisting liquid and gas phases are different. In particular, we are able to express the correlation function exactly as a simple average of contributions arising from two effective Ising-symmetric systems referred to as the symmetric gas and symmetric liquid. When combined with our earlier results, this generates analytical approximations for the correlation function and the local structure factor, which are near indistinguishable from the numerical solution to the Ornstein-Zernike equations over the whole range of wave vectors. Our results highlight how asymmetry affects the correlation function structure and describes the crossover from a long-ranged Goldstone mode to short-ranged properties determined by the local density as the wave vector increases.

Correlation-function structure in square-gradient models of the liquid-gas interface: Exact results and reliable approximations

In a recent article, we described how the microscopic structure of density-density correlations in the fluid interfacial region, for systems with short-ranged forces, can be understood by considering the resonances of the local structure factor occurring at specific parallel wave vectors q [Nat. Phys. 15, 287 (2019)1745-247310.1038/s41567-018-0361-z]. Here we investigate this further by comparing approximations for the local structure factor and pair correlation function against three new examples of analytically solvable models within square-gradient theory. Our analysis further demonstrates that these approximations describe the pair correlation function and structure factor across the whole spectrum of wave vectors, encapsulating the crossover from the Goldstone mode divergence (at small q) to bulklike behavior (at larger q). As shown, these approximations are exact for some square-gradient model potentials and never more than a few percent inaccurate for the others. Additionally, we show that they describe very accurately the correlation function structure for a model describing an interface near a tricritical point. In this case, there are no analytical solutions for the correlation functions, but the approximations are nearly indistinguishable from the numerical solutions of the Ornstein-Zernike equation.

Random phase approximation for the non-uniform Yukawa fluid

Mean-field density functional theory can be used to estimate the free energy of non-uniform fluids. The second functional derivative with respect to density of the free energy is related to the direct correlation function of the fluid and, in principle, this can be inverted to find an improved approximation for the pair correlation function and hence the free energy, the so-called 'random phase approximation'. If the repulsive molecular interaction is approximated by the local density approximation and the attractive interaction is assumed to be of the Yukawa form, the problem reduces to that of finding the eigenvalues of Schrödinger-like equations, which, for certain models (such as the 'Φ⁴ model'), can be done analytically in the planar case. The relationship between this approach and field theoretical treatment of the vapour-liquid interface is discussed. The ultraviolet divergence of the expression can be eliminated by separating the first term in the expansion, although quantitative results still depend on the behaviour of the attractive potential in the repulsive core. In the case of a spherical droplet of radius R, correction terms to the free energy involving lnR appear due to (i) cluster translational invariance, (ii) the unstable mode corresponding to droplet growth, and (iii) capillary waves. The net effect of these terms is to modify the classical expression for the nucleation rate by a factor proportional to R ^4/3.

Density correlation in liquid surfaces: Bedeaux-Weeks high order terms and non capillary wave background

We present Molecular Dynamics (MD) simulations of liquid-vapor surfaces, and their Intrinsic Sampling Method analysis, to get a quantitative test for the theoretical prediction of the capillary wave (CW) effects on density correlation done by Bedeaux and Weeks (BW) in 1985. The results are contrasted with Wertheim's proposal which is the first term in BW series and are complemented with a (formally defined and computational accessible) proposal for the background of non-CW fluctuations. Our conclusion is that BW theory is both accurate and needed since it may differ significantly from Wertheim's proposal. We discuss the implications for the analysis of experimental X-ray surface diffraction data and MD simulations.

Curvature corrections to the nonlocal interfacial model for short-ranged forces

In this paper we revisit the derivation of a nonlocal interfacial Hamiltonian model for systems with short-ranged intermolecular forces. Starting from a microscopic Landau-Ginzburg-Wilson Hamiltonian with a double-parabola potential, we reformulate the derivation of the interfacial model using a rigorous boundary integral approach. This is done for three scenarios: a single fluid phase in contact with a nonplanar substrate (i.e., wall); a free interface separating coexisting fluid phases (say, liquid and gas); and finally a liquid-gas interface in contact with a nonplanar confining wall, as is applicable to wetting phenomena. For the first two cases our approaches identifies the correct form of the curvature corrections to the free energy and, for the case of a free interface, it allows us to recast these as an interfacial self-interaction as conjectured previously in the literature. When the interface is in contact with a substrate our approach similarly identifies curvature corrections to the nonlocal binding potential, describing the interaction of the interface and wall, for which we propose a generalized and improved diagrammatic formulation.

Capillary waves as eigenmodes of the density correlation at liquid surfaces

We analyze the density correlations in a liquid-vapor surface to establish a quantitative connection between the Density Functional (DF) formalism, Molecular Dynamic (MD) simulations, and the Capillary Wave (CW) theory. Instead of the integrated structure factor, we identify the CW fluctuations as eigenmodes of the correlation function. The square-gradient DF approximation appears as fully consistent with the use of the thermodynamic surface tension to describe the surface fluctuations for any wavevector because it misses the upper cutoff in the surface Hamiltonian from the merging of the CW mode with the non-CW band. This mesoscopic cutoff may be accurately predicted from the main peak in the structure factor of the bulk liquid. We explore the difference between the full density-density correlation mode and the bare CW that represents the correlation between the corrugation of the intrinsic surface and the density at the interfacial region. The non-local decay of the CW effects, predicted from DF analysis and observed in MD simulations with the intrinsic sampling method, is found to characterize the bare CW fluctuations, which also require a wavevector-dependent surface tension.

Structure of electric double layers in capacitive systems and to what extent (classical) density functional theory describes it

Ongoing scientific interest is aimed at the properties and structure of electric double layers (EDLs), which are crucial for capacitive energy storage, water treatment, and energy harvesting technologies like supercapacitors, desalination devices, blue engines, and thermocapacitive heat-to-current converters. A promising tool to describe their physics on a microscopic level is (classical) density functional theory (DFT), which can be applied in order to analyze pair correlations and charge ordering in the primitive model of charged hard spheres. This simple model captures the main properties of ionic liquids and solutions and it predicts many of the phenomena that occur in EDLs. The latter often lead to anomalous response in the differential capacitance of EDLs. This work constructively reviews the powerful theoretical framework of DFT and its recent developments regarding the description of EDLs. It explains to what extent current approaches in DFT describe structural ordering and in-plane transitions in EDLs, which occur when the corresponding electrodes are charged. Further, the review briefly summarizes the history of modeling EDLs, presents applications, and points out limitations and strengths in present theoretical approaches. It concludes that DFT as a sophisticated microscopic theory for ionic systems is expecting a challenging but promising future in both fundamental research and applications in supercapacitive technologies.

Capillary waves and the decay of density correlations at liquid surfaces

Wertheim predicted strong density-density correlations at free liquid surfaces, produced by capillary wave fluctuations of the interface [M. S. Wertheim, J. Chem. Phys. 65, 2377 (1976)JCPSA60021-960610.1063/1.433352]. That prediction has been used to search for a link between capillary wave (CW) theory and density functional (DF) formalism for classical fluids. In particular, Parry et al. have recently analyzed the decaying tails of these CW effects moving away from the interface as a clue for the extended CW theory [A. O. Parry et al., J. Phys.: Condens. Matter 28, 244013 (2016)JCOMEL0953-898410.1088/0953-8984/28/24/244013], beyond the strict long-wavelength limit studied by Wertheim. Some apparently fundamental inconsistencies between the CW and the DF theoretical views of the fluid interfaces arose from the asymptotic analysis of the CW signal. In this paper we revisit the problem of the CW asymptotic decay with a separation of local non-CW surface correlation effects from those that are a truly nonlocal propagation of the CW fluctuations from the surface towards the liquid bulk.

Capillary wave Hamiltonian for the Landau-Ginzburg-Wilson density functional

We study the link between the density functional (DF) formalism and the capillary wave theory (CWT) for liquid surfaces, focused on the Landau-Ginzburg-Wilson (LGW) model, or square gradient DF expansion, with a symmetric double parabola free energy, which has been extensively used in theoretical studies of this problem. We show the equivalence between the non-local DF results of Parry and coworkers and the direct evaluation of the mean square fluctuations of the intrinsic surface, as is done in the intrinsic sampling method for computer simulations. The definition of effective wave-vector dependent surface tensions is reviewed and we obtain new proposals for the LGW model. The surface weight proposed by Blokhuis and the surface mode analysis proposed by Stecki provide consistent and optimal effective definitions for the extended CWT Hamiltonian associated to the DF model. A non-local, or coarse-grained, definition of the intrinsic surface provides the missing element to get the mesoscopic surface Hamiltonian from the molecular DF description, as had been proposed a long time ago by Dietrich and coworkers.

The local structure factor near an interface; beyond extended capillary-wave models

We investigate the local structure factor S (z;q) at a free liquid-gas interface in systems with short-ranged intermolecular forces and determine the corrections to the leading-order, capillary-wave-like, Goldstone mode divergence of S (z;q) known to occur for parallel (i.e. measured along the interface) wavevectors [Formula: see text]. We show from explicit solution of the inhomogeneous Ornstein-Zernike equation that for distances z far from the interface, where the profile decays exponentially, S (z;q) splits unambiguously into bulk and interfacial contributions. On each side of the interface, the interfacial contributions can be characterised by distinct liquid and gas wavevector dependent surface tensions, [Formula: see text] and [Formula: see text], which are determined solely by the bulk two-body and three-body direct correlation functions. At high temperatures, the wavevector dependence simplifies and is determined almost entirely by the appropriate bulk structure factor, leading to positive rigidity coefficients. Our predictions are confirmed by explicit calculation of S (z;q) within square-gradient theory and the Sullivan model. The results for the latter predict a striking temperature dependence for [Formula: see text] and [Formula: see text], and have implications for fluctuation effects. Our results account quantitatively for the findings of a recent very extensive simulation study by Höfling and Dietrich of the total structure factor in the interfacial region, in a system with a cut-off Lennard-Jones potential, in sharp contrast to extended capillary-wave models which failed completely to describe the simulation results.

Anisotropic pair correlations in binary and multicomponent hard-sphere mixtures in the vicinity of a hard wall: A combined density functional theory and simulation study

The fundamental measure approach to classical density functional theory has been shown to be a powerful tool to predict various thermodynamic properties of hard-sphere systems. We employ this approach to determine not only one-particle densities but also two-particle correlations in binary and six-component mixtures of hard spheres in the vicinity of a hard wall. The broken isotropy enables us to carefully test a large variety of theoretically predicted two-particle features by quantitatively comparing them to the results of Brownian dynamics simulations. Specifically, we determine and compare the one-particle density, the total correlation functions, their contact values, and the force distributions acting on a particle. For this purpose, we follow the compressibility route and theoretically calculate the direct correlation functions by taking functional derivatives. We usually observe an excellent agreement between theory and simulations, except for small deviations in cases where local crystal-like order sets in. Our results set the course for further investigations on the consistency of functionals as well as for structural analysis on, e.g., the primitive model. In addition, we demonstrate that due to the suppression of local crystallization, the predictions of six-component mixtures are better than those in bidisperse or monodisperse systems. Finally, we are confident that our results of the structural modulations induced by the wall lead to a deeper understanding of ordering in anisotropic systems in general, the onset of heterogeneous crystallization, caging effects, and glassy dynamics close to a wall, as well as structural properties in systems with confinement.

Mesoscopic Hamiltonian for the fluctuations of adsorbed Lennard-Jones liquid films

We use Monte Carlo simulations of a Lennard-Jones fluid adsorbed on a short-range planar wall substrate to study the fluctuations in the thickness of the wetting layer, and we get a quantitative and consistent characterization of their mesoscopic Hamiltonian, H[ξ]. We have observed important finite-size effects, which were hampering the analysis of previous results obtained with smaller systems. The results presented here support an appealing simple functional form for H[ξ], close but not exactly equal to the theoretical nonlocal proposal made on the basis a generic density-functional analysis by Parry and coworkers. We have analyzed systems under different wetting conditions, as a proof of principle for a method that provides a quantitative bridge between the molecular interactions and the phenomenology of wetting films at mesoscopic scales.

Liquid-gas asymmetry and the wave-vector-dependent surface tension

Attempts to extend the capillary-wave theory of fluid interfacial fluctuations to microscopic wavelengths, by introducing an effective wave-vector (q)-dependent surface tension σeff(q), have encountered difficulties. There is no consensus as to even the shape of σeff(q). By analyzing a simple density functional model of the liquid-gas interface, we identify different schemes for separating microscopic observables into background and interfacial contributions. In order for the backgrounds of the density-density correlation function and local structure factor to have a consistent and physically meaningful interpretation in terms of weighted bulk gas and liquid contributions, the background of the total structure factor must be characterized by a microscopic q-dependent length ζ(q) not identified previously. The necessity of including the q dependence of ζ(q) is illustrated explicitly in our model and has wider implications; i.e., in typical experimental and simulation studies, an indeterminacy in ζ(q) will always be present, reminiscent of the cutoff used in capillary-wave theory. This leads inevitably to a large uncertainty in the q dependence of σeff(q).