Effective Hamiltonian for liquid-vapor interfaces

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.

Surface tension of bulky colloids, capillarity under gravity, and the microscopic origin of the Kardar-Parisi-Zhang equation

Experimental measurements of the surface tension of colloidal interfaces have long been in conflict with computer simulations. In this Letter we show that the surface tension of colloids as measured by surface fluctuations picks up a gravity-dependent contribution which removes the discrepancy. The presence of this term puts a strong constraint on the structure of the interface which allows one to identify corrections to the fundamental equation of equilibrium capillarity and deduce bottom up the microscopic origin of a growth model with close relation to the Kardar-Parisi-Zhang equation.

Local structure of liquid/vapour interfaces approaching the critical point

Investigating the structure of fluid interfaces at high temperatures is a particularly delicate task that requires effective ways of discriminating liquid from vapour and identifying the location of the liquid phase boundary, thereby allowing to distinguish intrinsic from capillary fluctuations. Several numerical approaches require introducing a coarse-graining length scale, often heuristically chosen to be the molecular size, to determine the location of the liquid phase boundary. Here, we propose an alternative rationale for choosing this coarse-graining length scale; we require the average position of the local liquid phase dividing surface to match its flat, macroscopic counterpart. We show that this approach provides additional insight into the structure of the liquid/vapour interface, suggesting the presence of another length scale beyond the bulk correlation one that plays an important role in determining the interface structure.

Surface viscosity in simple liquids

The response of Newtonian liquids to small perturbations is usually considered to be fully described by homogeneous transport coefficients like shear and dilatational viscosity. However, the presence of strong density gradients at the liquid/vapor boundary of fluids hints at the possible existence of an inhomogeneous viscosity. Here, we show that a surface viscosity emerges from the collective dynamics of interfacial layers in molecular simulations of simple liquids. We estimate the surface viscosity to be 8-16 times smaller than that of the bulk fluid at the thermodynamic point considered. This result can have important implications for reactions at liquid surfaces in atmospheric chemistry and catalysis.

Comparative study of force-based classical density functional theory

We reexamine results obtained with the recently proposed density functional theory framework based on forces (force-DFT) [S. M. Tschopp et al., Phys. Rev. E 106, 014115 (2022)2470-004510.1103/PhysRevE.106.014115]. We compare inhomogeneous density profiles for hard sphere fluids to results from both standard density functional theory and from computer simulations. Test situations include the equilibrium hard sphere fluid adsorbed against a planar hard wall and the dynamical relaxation of hard spheres in a switched harmonic potential. The comparison to grand canonical Monte Carlo simulation profiles shows that equilibrium force-DFT alone does not improve upon results obtained with the standard Rosenfeld functional. Similar behavior holds for the relaxation dynamics, where we use our event-driven Brownian dynamics data as benchmark. Based on an appropriate linear combination of standard and force-DFT results, we investigate a simple hybrid scheme which rectifies these deficiencies in both the equilibrium and the dynamical case. We explicitly demonstrate that although the hybrid method is based on the original Rosenfeld fundamental measure functional, its performance is comparable to that of the more advanced White Bear theory.

Interfacial stiffness of nematic-smectic B interface in Gay-Berne liquid crystals using capillary wave theory

The interfacial stiffness for nematic-smectic B (nm-smB) interface in a liquid crystalline (LC) material is calculated using Capillary Wave Theory (CWT) and molecular dynamics simulations. The Gay-Berne (GB) pair potential with parameters κ, κ', μ, and ν equal to 3, 5, 2, and 1 is used to model the LC material. Using a smart three-step recipe, we have obtained an nm-smB phase coexistence in our simulations where the nm and smB directors are nearly parallel to each other and perpendicular to the interface normal. The density profiles are used to compute the nm-smB coexisting density range, the interfacial width, and its position. The smectic phase is differentiated from the nematic phase by using the local bond order parameter (q₆q₆), which has helped us to demonstrate that the interface is indeed rough. Finally, the interfacial stiffness of the nm-smB interface is computed by following the CWT analysis and is found to be γ̃_nm-smB=0.39861k_BT/σ_ee ²=0.04429/σ_ss ², where σ_ee and σ_ss are the length and diameter of the GB LC particles.

How surface stress transforms surface profiles and adhesion of rough elastic bodies

The surface of soft solids carries a surface stress that tends to flatten surface profiles. For example, surface features on a soft solid, fabricated by moulding against a stiff-patterned substrate, tend to flatten upon removal from the mould. In this work, we derive a transfer function in an explicit form that, given any initial surface profile, shows how to compute the shape of the corresponding flattened profile. We provide analytical results for several applications including flattening of one-dimensional and two-dimensional periodic structures, qualitative changes to the surface roughness spectrum, and how that strongly influences adhesion.

Self-Sustaining 3D Thin Liquid Films in Ambient Environments

Thin liquid films (TLF) have fundamental and technological importance ranging from the thermodynamics of cell membranes to the safety of light-water cooled nuclear reactors. The creation of stable water TLFs, however, is very difficult. In this paper, the realization of thin liquid films of water with custom 3D geometries that persist indefinitely in ambient environments is reported. The wetting films are generated using microscale "mounts" fed by microfluidic channels with small feature sizes and large aspect ratios. These devices are fabricated with a custom 3D printer and resin, which were developed to print high resolution microfluidic geometries as detailed in Reference 26. By modifying the 3D-printed polymer to be hydrophilic and taking advantage of well-known wetting principles and capillary effects, self-sustaining microscale "water fountains" are constructed that continuously replenish water lost to evaporation while relying on surface tension to stabilize their shape. To the authors' knowledge, this is the first demonstration of stable sub-micron thin liquid films (TLFs) of pure water on curved 3D geometries.

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.

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.

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.

Effect of an external field on capillary waves in a dipolar fluid

The role of an external field on capillary waves at the liquid-vapor interface of a dipolar fluid is investigated using molecular dynamics simulations. For fields parallel to the interface, the interfacial width squared increases linearly with respect to the logarithm of the size of the interface across all field strengths tested. The value of the slope decreases with increasing field strength, indicating that the field dampens the capillary waves. With the inclusion of the parallel field, the surface stiffness increases with increasing field strength faster than the surface tension. For fields perpendicular to the interface, the interfacial width squared is linear with respect to the logarithm of the size of the interface for small field strengths, and the surface stiffness is less than the surface tension. Above a critical field strength that decreases as the size of the interface increases, the interface becomes unstable due to the increased amplitude of the capillary waves.

Predictions of homogeneous nucleation rates for n-alkanes accounting for the diffuse phase interface and capillary waves

Homogeneous droplet nucleation has been studied for almost a century but has not yet been fully understood. In this work, we used the density gradient theory (DGT) and considered the influence of capillary waves (CWs) on the predicted size-dependent surface tensions and nucleation rates for selected n-alkanes. The DGT model was completed by an equation of state (EoS) based on the perturbed-chain statistical associating fluid theory and compared to the classical nucleation theory and the Peng-Robinson EoS. It was found that the critical clusters are practically free of CWs because they are so small that even the smallest wavelengths of CWs do not fit into their finite dimensions. The CWs contribute to the entropy of the system and thus decrease the surface tension. A correction for the effect of CWs on the surface tension is presented. The effect of the different EoSs is relatively small because by a fortuitous coincidence their predictions are similar in the relevant range of critical cluster sizes. The difference of the DGT predictions to the classical nucleation theory computations is important but not decisive. Of the effects investigated, the most pronounced is the suppression of CWs which causes a sizable decrease of the predicted nucleation rates. The major difference between experimental nucleation rate data and theoretical predictions remains in the temperature dependence. For normal alkanes, this discrepancy is much stronger than observed, e.g., for water. Theoretical corrections developed here have a minor influence on the temperature dependency. We provide empirical equations correcting the predicted nucleation rates to values comparable with experiments.

Capillary wave theory of adsorbed liquid films and the structure of the liquid-vapor interface

In this paper we try to work out in detail the implications of a microscopic theory for capillary waves under the assumption that the density is given along lines normal to the interface. Within this approximation, which may be justified in terms of symmetry arguments, the Fisk-Widom scaling of the density profile holds for frozen realizations of the interface profile. Upon thermal averaging of capillary wave fluctuations, the resulting density profile yields results consistent with renormalization group calculations in the one-loop approximation. The thermal average over capillary waves may be expressed in terms of a modified convolution approximation where normals to the interface are Gaussian distributed. In the absence of an external field we show that the phenomenological density profile applied to the square-gradient free energy functional recovers the capillary wave Hamiltonian exactly. We extend the theory to the case of liquid films adsorbed on a substrate. For systems with short-range forces, we recover an effective interface Hamiltonian with a film height dependent surface tension that stems from the distortion of the liquid-vapor interface by the substrate, in agreement with the Fisher-Jin theory of short-range wetting. In the presence of long-range interactions, the surface tension picks up an explicit dependence on the external field and recovers the wave vector dependent logarithmic contribution observed by Napiorkowski and Dietrich. Using an error function for the intrinsic density profile, we obtain closed expressions for the surface tension and the interface width. We show the external field contribution to the surface tension may be given in terms of the film's disjoining pressure. From literature values of the Hamaker constant, it is found that the fluid-substrate forces may be able to double the surface tension for films in the nanometer range. The film height dependence of the surface tension described here is in full agreement with results of the capillary wave spectrum obtained recently in computer simulations, and the predicted translation mode of surface fluctuations reproduces to linear order in field strength an exact solution of the density correlation function for the Landau-Ginzburg-Wilson Hamiltonian in an external field.

Thermal fluctuations of an interface near a contact line

The effect of thermal fluctuations near a contact line of a liquid interface partially wetting an impenetrable substrate is studied analytically and numerically. Promoting both the interface profile and the contact line position to random variables, we explore the equilibrium properties of the corresponding fluctuating contact line problem based on an interfacial Hamiltonian involving a "contact" binding potential. To facilitate an analytical treatment, we consider the case of a one-dimensional interface. The effective boundary condition at the contact line is determined by a dimensionless parameter that encodes the relative importance of thermal energy and substrate energy at the microscopic scale. We find that this parameter controls the transition from a partial wetting to a pseudopartial wetting state, the latter being characterized by a thin prewetting film of fixed thickness. In the partial wetting regime, instead, the profile typically approaches the substrate via an exponentially thinning prewetting film. We show that, independently of the physics at the microscopic scale, Young's angle is recovered sufficiently far from the substrate. The fluctuations of the interface and of the contact line give rise to an effective disjoining pressure, exponentially decreasing with height. Fluctuations therefore provide a regularization of the singular contact forces occurring in the corresponding deterministic problem.

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.

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

Order Parameters and Algorithmic Approaches for Detection and Demarcation of Interfaces in Hydrate-Fluid and Ice-Fluid Systems

Some aspects of the use of order parameter fields in molecular dynamics simulations to delimit solid phases containing water, namely ice and hydrate, in both hydrophilic and hydrophobic fluids are examined; this includes the influences of rectangular meshes and of filtering on the quality of these parameters. Three order parameters are studied: the mass density, ρ; an angular tetrahedrality measure, Sg (Chau and Hardwick, Mol. Phys. 1998, 93, 511); and the water-dimer dihedral angle, F4 (Rodger et al. Fluid Phase Equilib. 1996, 116, 326). The parameters are studied to find their ability to distinguish between bulk phases, their consistency in different environments, their noise susceptibility, and their ability to demarcate the interface region. Spatial sampling and filtering are covered in detail, and some temporal features are illustrated by using autocorrelation maps. The parameters are employed to determine the position of interfaces as functions of time and, with the capillary wave fluctuation method (Hoyt et al. Phys. Rev. Lett. 2001, 86, 5530; Math. Comput. Simul. 2010, 80, 1382), to estimate solid-fluid interfacial stiffnesses, with partial success for the hydrophilic/hydrophobic-type interfaces.

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

We study the density-density correlation function G(r, r') in the interfacial region of a fluid (or Ising-like magnet) with short-ranged interactions using square gradient density functional theory. Adopting a simple double parabola approximation for the bulk free-energy density, we first show that the parallel Fourier transform G(z, z'; q) and local structure factor S(z; q) separate into bulk and excess contributions. We attempt to account for both contributions by deriving an interfacial Hamiltonian, characterised by a wavevector dependent surface tension σ(q), and then reconstructing density correlations from correlations in the interface position. We show that the standard crossing criterion identification of the interface, as a surface of fixed density (or magnetization), does not explain the separation of G(z, z'; q) and the form of the excess contribution. We propose an alternative definition of the interface position based on the properties of correlations between points that 'float' with the surface and show that this describes the full q and z dependence of the excess contributions to both G and S. However, neither the 'crossing-criterion' nor the new 'floating interface' definition of σ(q) are quantities directly measurable from the total structure factor S(tot)(q) which contains additional q dependence arising from the non-local relation between fluctuations in the interfacial position and local density. Since it is the total structure factor that is measured experimentally or in simulations, our results have repercussions for earlier attempts to extract and interpret σ(q).

Logarithmic finite-size effects on interfacial free energies: phenomenological theory and Monte Carlo studies

The computation of interfacial free energies between coexisting phases (e.g., saturated vapor and liquid) by computer simulation methods is still a challenging problem due to the difficulty of an atomistic identification of an interface and interfacial fluctuations on all length scales. The approach to estimate the interfacial tension from the free-energy excess of a system with interfaces relative to corresponding single-phase systems does not suffer from the first problem but still suffers from the latter. Considering d-dimensional systems with interfacial area L(d-1) and linear dimension L(z) in the direction perpendicular to the interface, it is argued that the interfacial fluctuations cause logarithmic finite-size effects of order ln(L)/L(d-1) and order ln(L(z))/L(d-1), in addition to regular corrections (with leading-order const/L(d-1)). A phenomenological theory predicts that the prefactors of the logarithmic terms are universal (but depend on the applied boundary conditions and the considered statistical ensemble). The physical origin of these corrections are the translational entropy of the interface as a whole, "domain breathing" (coupling of interfacial fluctuations to the bulk order parameter fluctuations of the coexisting domains), and capillary waves. Using a new variant of the ensemble switch method, interfacial tensions are found from Monte Carlo simulations of d = 2 and d = 3 Ising models and a Lennard-Jones fluid. The simulation results are fully consistent with the theoretical predictions.

Disjoining pressure and the film-height-dependent surface tension of thin liquid films: new insight from capillary wave fluctuations

In this paper we review simulation and experimental studies of thermal capillary wave fluctuations as an ideal means for probing the underlying disjoining pressure and surface tensions, and more generally, fine details of the Interfacial Hamiltonian Model. We discuss recent simulation results that reveal a film-height-dependent surface tension not accounted for in the classical Interfacial Hamiltonian Model. We show how this observation may be explained bottom-up from sound principles of statistical thermodynamics and discuss some of its implications.

Effect of dispersion forces on the capillary-wave fluctuations of liquid surfaces

We present molecular dynamics evidence for the nonanalytic effects of the long-range dispersion forces on the capillary waves fluctuations of a Lennard-Jones liquid surface. The results of the intrinsic sampling method, for the analysis of the instantaneous interfacial shape, are obtained in large systems for several cut-off distances of the potential tail, and they show good agreement with the theoretical prediction by Napiórkowski and Dietrich, based on a density functional analysis. The enhancement of the capillary waves is quantified to be within 1% for a simple liquid near its triple point.

Putting water on a lattice: the importance of long wavelength density fluctuations in theories of hydrophobic and interfacial phenomena

The physics of air-water interfaces plays a central role in modern theories of the hydrophobic effect. Implementing these theories, however, has been hampered by the difficulty of addressing fluctuations in the shape of such soft interfaces. We show that this challenge is a fundamental consequence of mapping long wavelength density variations onto discrete degrees of freedom. Drawing from studies of surface roughness in lattice models, we account for the resulting nonlinearities simply but accurately. Simulations show that this approach captures complex solvation behaviors quantitatively.

Intrinsic fluid interfaces and nonlocality

We present results of an extensive molecular dynamics simulation of the structure and fluctuations of a liquid-gas interface, close to its triple point, in a system with cutoff Lennard-Jones interactions. The equilibrium density profile, averaged and (shape dependent) constrained intrinsic density profiles together with the fluctuations of the interfacial shape are extracted using an intrinsic sampling method. The correlation between fluctuations in the interfacial shape and in the intrinsic density show that the latter is not due to rigid translations of some underlying profile, as is most commonly assumed. Instead, over the whole range of wavelengths from the system size down to the molecular diameter, we see wave-vector dependent behavior in good agreement with a nonlocal interfacial Hamiltonian theory specifying the shape dependence of the intrinsic profiles.

Capillary fluctuations and film-height-dependent surface tension of an adsorbed liquid film

Our understanding of both structure and dynamics of adsorbed liquids heavily relies on the capillary wave Hamiltonian, but a thorough test of this model is still lacking. Here we study the capillary wave fluctuations of a liquid film with short-range forces adsorbed on a solid exhibiting van der Waals interactions. We show for the first time that the measured capillary wave spectrum right above the first order wetting transition provides an interface potential consistent with independent calculations from thermodynamic integration. However, the surface tension exhibits an oscillatory film thick dependence which reveals a hitherto unnoticed capillary wave broadening mechanism beyond mere interfacial displacements.

Mode expansion for the density profiles of crystal-fluid interfaces: hard spheres as a test case

We present a technique for analyzing the full three-dimensional density profiles of planar crystal-fluid interfaces in terms of density modes. These density modes can also be related to crystallinity order parameter profiles which are used in coarse-grained, phase field type models of the statics and dynamics of crystal-fluid interfaces and are an alternative to crystallinity order parameters extracted from simulations using local crystallinity criteria. We illustrate our results for the hard sphere system using finely resolved, three-dimensional density profiles from a density functional theory of fundamental measure type.

Intrinsic profiles and the structure of liquid surfaces

We present a brief review of the advances made in the characterization of liquid surfaces over the last decade. We focus particularly on the links between the capillary wave theory, the density functional formalism and the direct evaluation of the intrinsic density profiles from computer simulations. A new perspective of the liquid surfaces is appearing, with a sharper view of their molecular structure, which opens new challenges for theoretical and experimental studies. Novel results on the intrinsic interfacial structure of molten salt liquid-vapor interfaces are presented.

Spectral properties of thermal fluctuations on simple liquid surfaces below shot-noise levels

We study the spectral properties of thermal fluctuations on simple liquid surfaces, sometimes called ripplons. Analytical properties of the spectral function are investigated and are shown to be composed of regions with simple analytic behavior with respect to the frequency or the wave number. The derived expressions are compared to spectral measurements performed orders of magnitude below shot-noise levels, which is achieved using a novel noise reduction method. The agreement between the theory of thermal surface fluctuations and the experiment is found to be excellent, elucidating the spectral properties of the surface fluctuations. The measurement method requires relatively only a small sample both spatially (few μm) and temporally (~20 s). The method also requires relatively weak light power (~0.5 mW) so that it has a broad range of applicability, including local measurements, investigations of time-dependent phenomena, and noninvasive measurements.

The order of filling transitions in acute wedges

Using a square-gradient density functional model we test the prediction that the filling transition for a fluid in a wedge geometry changes from continuous to first-order as the wedge becomes more acute. Our numerical findings confirm such a change of order, but the value of the tilt angle at which it occurs, α* ≈ 45°, is considerably smaller than the original theoretical prediction. We critically reassess this work, which was based on allowing for the self-interaction of the fluid interface, and argue that the interfacial curvature and effective wavevector dependent surface tension can further lower the predicted value of α*, in keeping with our numerical findings. Interfacial fluctuation effects, occurring beyond mean-field level, are also discussed using effective Hamiltonian theory and are shown to substantially increase the value of α*.

Effective Hamiltonian of liquid-vapor curved interfaces in mean field

We analyze a one-component simple fluid in a liquid-vapor coexistence state, which forms an arbitrarily curved interface. By using an approach based on density functional theory, we obtain an exact and simple expression for the grand potential at the level of mean field approximation that depends on the density profile and the short-range interaction potential. By introducing the step-function approximation for the density profile, and using general geometric arguments, we expand the grand potential in powers of the principal curvatures of the surface and find consistency with the Helfrich phenomenological model in the second order approximation.