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

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.

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.

Surface reconstruction amendment to the intrinsic sampling method

The intrinsic sampling method (ISM) is a powerful tool that allows the exploration of interfacial properties from molecular simulations by fitting a function that represents the local boundary between two phases. However, owing to the non-physical nature of an "intrinsic" surface, there remains an ambiguity surrounding the comparison of theoretical properties with the physical world. It is therefore important that the ISM remains internally consistent when reproducing simulated properties which match experiments, such as the surface tension or interfacial density distribution. We show that the current ISM procedure causes an over-fitting of the surface to molecules in the interface region, leading to a biased distribution of curvature at these molecular coordinates. We assert that this biased distribution is a cause of the disparity between predicted interfacial densities upon convolution to a laboratory frame, an artefact which has been known to exist since the development of the ISM. We present an improvement to the fitting procedure of the ISM in an attempt to alleviate the ambiguity surrounding the true nature of an intrinsic surface. Our "surface reconstruction" method is able to amend the shape of the interface so as to reproduce the global curvature distribution at all sampled molecular coordinates. We present the effects that this method has on the ISM predicted structure of a simulated Lennard-Jones fluid air-liquid interface. Additionally, we report an unexpected relationship between surface thermodynamic predictions of our reconstructed ISM surfaces and those of extended capillary wave theory, which is of current interest.

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.

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.

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.