Quantifying Density Fluctuations in Water at a Hydrophobic Surface: Evidence for Critical Drying

What is the best simulation approach for measuring local density fluctuations near solvo-/hydrophobes?

Measurements of local density fluctuations are crucial to characterizing the interfacial properties of equilibrium fluids. A specific case that has been well-explored involves the heightened compressibility of water near hydrophobic entities. Commonly, a spatial profile of local fluctuation strength is constructed from the measurements of the mean and variance of solvent particle number fluctuations in a set of contiguous subvolumes of the system adjacent to the solvo-/hydrophobe. An alternative measure proposed by Evans and Stewart [J. Phys.: Condens. Matter 27, 194111 (2015)] defines a local compressibility profile in terms of the chemical potential derivative of the spatial number density profile. Using Grand canonical Monte Carlo simulation, we compare and contrast the efficacy of these two approaches for a Lennard-Jones solvent at spherical and planar solvophobic interfaces and SPC/E water at a hydrophobic spherical solute. Our principal findings are as follows: (i) the local compressibility profile χ(r) of Evans and Stewart is considerably more sensitive to variations in the strength of local density fluctuations than the spatial fluctuation profile F(r) and can resolve much more detailed structure; and (ii) while the local compressibility profile is essentially independent of the choice of spatial discretization used to construct the profile, the spatial fluctuation profile exhibits a strong systematic dependence on the size of the subvolumes on which the profile is defined. We clarify the origin and nature of this finite-size effect.

Partial and complete wetting of droplets of active Brownian particles

We study wetting droplets formed of active Brownian particles in contact with a repulsive potential barrier, in a wedge geometry. Our numerical results demonstrate a transition between partially wet and completely wet states, as a function of the barrier height, analogous to the corresponding surface phase transition in passive fluids. We analyse partially wet configurations characterised by a nonzero contact angle θ between the droplet surface and the barrier including the average density profile and its fluctuations. These findings are compared with two equilibrium systems: a Lennard-Jones fluid and a simple contour model for a liquid-vapour interface. We locate the wetting transition where cos(θ) = 1, and the neutral state where cos(θ) = 0. We discuss the implications of these results for possible definitions of surface tensions in active fluids.

Local measures of fluctuations in inhomogeneous liquids: statistical mechanics and illustrative applications

We show in detail how three one-body fluctuation profiles, namely the local compressibility, the local thermal susceptibility, and the reduced density, can be obtained from a statistical mechanical many-body description of classical particle-based systems. We present several different and equivalent routes to the definition of each fluctuation profile, facilitating their explicit numerical calculation in inhomogeneous equilibrium systems. This underlying framework is used for the derivation of further properties such as hard wall contact theorems and novel types of inhomogeneous one-body Ornstein-Zernike equations. The practical accessibility of all three fluctuation profiles is exemplified by grand canonical Monte Carlo simulations that we present for hard sphere, Gaussian core and Lennard-Jones fluids in confinement.

Understanding the physics of hydrophobic solvation

Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations. The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jones solvent, and Grand Canonical Monte Carlo (GCMC) simulations of a monatomic water (mw) model. We argue that the density depletion and enhanced fluctuations are near-critical phenomena. Specifically, we show that they can be viewed as remnants of the critical drying surface phase transition that occurs at bulk liquid-vapor coexistence in the macroscopic planar limit, i.e., as the solute radius R_s → ∞. Focusing on the radial density profile ρ(r) and a sensitive spatial measure of fluctuations, the local compressibility profile χ(r), our binding potential analysis provides explicit predictions for the manner in which the key features of ρ(r) and χ(r) scale with R_s, the strength of solute-water attraction ɛ_sf, and the deviation from liquid-vapor coexistence of the chemical potential, δμ. These scaling predictions are confirmed by our DFT calculations and GCMC simulations. As such, our theory provides a firm basis for understanding the physics of hydrophobic solvation.

Force density functional theory in- and out-of-equilibrium

When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force, which plays an essential role in determining both the equilibrium and dynamic properties of interacting fluids. We present a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions. By consideration of local Noether-invariance against spatial distortion of the system we demonstrate the fundamental status of the Yvon-Born-Green (YBG) equation as a local force-balance within the fluid. Using the inhomogeneous Ornstein-Zernike equation we show that the two-body correlations are density functionals and, thus, that the average interparticle force entering the YBG equation is also a functional of the one-body density. The force-based theory we develop provides an alternative to standard density functional theory for the study of inhomogeneous systems both in- and out-of-equilibrium. We compare force-based density profiles to the results of the standard potential-based (dynamical) density functional theory. In-equilibrium, we confirm both analytically and numerically that the standard approach yields profiles that are consistent with the compressibility pressure, whereas the force-density functional gives profiles consistent with the virial pressure. For both approaches we explicitly prove the hard-wall contact theorem that connects the value of the density profile at the hard-wall with the bulk pressure. The structure of the theory offers deep insights into the nature of correlation in dense and inhomogeneous systems.

Why Noether's theorem applies to statistical mechanics

Noether's theorem is familiar to most physicists due its fundamental role in linking the existence of conservation laws to the underlying symmetries of a physical system. Typically the systems are described in the particle-based context of classical mechanics or on the basis of field theory. We have recently shown (2021Commun. Phys.4176) that Noether's reasoning also applies to thermal systems, where fluctuations are paramount and one aims for a statistical mechanical description. Here we give a pedagogical introduction based on the canonical ensemble and apply it explicitly to ideal sedimentation. The relevant mathematical objects, such as the free energy, are viewed as functionals. This vantage point allows for systematic functional differentiation and the resulting identities express properties of both macroscopic average forces and molecularly resolved correlations in many-body systems, both in and out-of-equilibrium, and for active Brownian particles. To provide further background, we briefly describe the variational principles of classical density functional theory, of power functional theory, and of classical mechanics.

The coexistence curve and surface tension of a monatomic water model

We study the monatomic water model of Molinero and Moore the grand canonical ensemble Monte Carlo simulation. Measurements of the probability distribution of the number density obtained via multicanonical sampling and histogram reweighting provide accurate estimates of the temperature dependence of both the liquid-vapor coexistence densities and the surface tension. Using finite-size scaling methods, we locate the liquid-vapor critical point at T_c = 917.6 K, ρ_c = 0.311 g cm^-3. When plotted in scaled variables, in order to test the law of corresponding states, the coexistence curve of monatomic water is close to that of real water. In this respect, it performs better than extended simple point charge (SPC/E), TIP4P, and TIP4P/2005 water.

Nanoconfined Fluids: Uniqueness of Water Compared to Other Liquids

Nanoconfinement can drastically change the behavior of liquids, puzzling us with counterintuitive properties. It is relevant in applications, including decontamination and crystallization control. However, it still lacks a systematic analysis for fluids with different bulk properties. Here we address this gap. We compare, by molecular dynamics simulations, three different liquids in a graphene slit pore: (1) A simple fluid, such as argon, described by a Lennard-Jones potential; (2) an anomalous fluid, such as a liquid metal, modeled with an isotropic core-softened potential; and (3) water, the prototypical anomalous liquid, with directional HBs. We study how the slit-pore width affects the structure, thermodynamics, and dynamics of the fluids. All the fluids show similar oscillating properties by changing the pore size. However, their free-energy minima are quite different in nature: (i) are energy-driven for the simple liquid; (ii) are entropy-driven for the isotropic core-softened potential; and (iii) have a changing nature for water. Indeed, for water, the monolayer minimum is entropy driven, at variance with the simple liquid, while the bilayer minimum is energy driven, at variance with the other anomalous liquid. Also, water has a large increase in diffusion for subnm slit pores, becoming faster than bulk. Instead, the other two fluids have diffusion oscillations much smaller than water, slowing down for decreasing slit-pore width. Our results, clarifying that water confined at the subnm scale behaves differently from other (simple or anomalous) fluids under similar confinement, are possibly relevant in nanopores applications, for example, in water purification from contaminants.

Chasing the Critical Wetting Transition. An Effective Interface Potential Method

Wettablity is one of the important characteristics defining a given surface. Here we show that the effective interface potential method of determining the wetting temperature, originally proposed by MacDowell and Müller for the surfaces exhibiting the first order wetting transition, can also be used to estimate the wetting temperature of the second order (continuous) wetting transition. Some selected other methods of determination of the wetting temperature are also discussed.

Fluctuation Profiles in Inhomogeneous Fluids

Three one-body profiles that correspond to local fluctuations in energy, in entropy, and in particle number are used to describe the equilibrium properties of inhomogeneous classical many-body systems. Local fluctuations are obtained from thermodynamic differentiation of the density profile or equivalently from average microscopic covariances. The fluctuation profiles follow from functional generators and they satisfy Ornstein-Zernike relations. Computer simulations reveal markedly different fluctuations in confined fluids with Lennard-Jones, hard sphere, and Gaussian core interactions.

Mesoscopic method to study water flow in nanochannels with different wettability

Molecular dynamics (MD) simulations is currently the most popular and credible tool to model water flow in nanoscale where the conventional continuum equations break down due to the dominance of fluid-surface interactions. However, current MD simulations are computationally challenging for the water flow in complex tube geometries or a network of nanopores, e.g., membrane, shale matrix, and aquaporins. We present a novel mesoscopic lattice Boltzmann method (LBM) for capturing fluctuated density distribution and a nonparabolic velocity profile of water flow through nanochannels. We incorporated molecular interactions between water and the solid inner wall into LBM formulations. Details of the molecular interactions were translated into true and apparent slippage, which were both correlated to the surface wettability, e.g., contact angle. Our proposed LBM was tested against 47 published cases of water flow through infinite-length nanochannels made of different materials and dimensions-flow rates as high as seven orders of magnitude when compared with predictions of the classical no-slip Hagen-Poiseuille (HP) flow. Using the developed LBM model, we also studied water flow through finite-length nanochannels with tube entrance and exit effects. Results were found to be in good agreement with 44 published finite-length cases in the literature. The proposed LBM model is nearly as accurate as MD simulations for a nanochannel, while being computationally efficient enough to allow implications for much larger and more complex geometrical nanostructures.

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

Enhanced fluctuation for pinned surface nanobubbles

By employing molecular dynamics simulations we investigate the fluctuation of surface nanobubbles immersed in liquid phase. Our simulation results indicate that in comparison with the surrounding liquid and nanobubble interior, the vapor-liquid or gas-liquid interface of nanobubbles always exhibits the largest compressibility, demonstrating the enhanced fluctuation for nanobubble interfaces. We also find that vapor surface nanobubbles and gas surface nanobubbles exhibit different fluctuation behaviors. For vapor nanobubbles that appear in overheated pure liquid, both density fluctuation and interface fluctuation are independent on the external pressure since the internal pressure remains constant at a given temperature. For gas nanobubbles that appear in gas supersaturated solution, the density fluctuation monotonously decreases with the increase of gas concentration, while the interface fluctuation shows a nonmonotonic variation. Departure from the intermediate gas concentration with the minimal interface fluctuation would enhance the fluctuation, which may finally lead to nanobubble destabilization. Finally, our simulation results indicate that the complicated interface fluctuation of surface nanobubbles comprises two different modes: interface deformation and interface oscillation, both of which display similar trends as that of the combined interface fluctuation.

A unified description of hydrophilic and superhydrophobic surfaces in terms of the wetting and drying transitions of liquids

Clarifying the factors that control the contact angle of a liquid on a solid substrate is a long-standing scientific problem pertinent across physics, chemistry, and materials science. Progress has been hampered by the lack of a comprehensive and unified understanding of the physics of wetting and drying phase transitions. Using various theoretical and simulational techniques applied to realistic fluid models, we elucidate how the character of these transitions depends sensitively on both the range of fluid-fluid and substrate-fluid interactions and the temperature. Our calculations uncover previously unrecognized classes of surface phase diagram which differ from that established for simple lattice models and often assumed to be universal. The differences relate both to the topology of the phase diagram and to the nature of the transitions, with a remarkable feature being a difference between drying and wetting transitions which persists even in the approach to the bulk critical point. Most experimental and simulational studies of liquids at a substrate belong to one of these previously unrecognized classes. We predict that while there appears to be nothing particularly special about water with regard to its wetting and drying behavior, superhydrophobic behavior should be more readily observable in experiments conducted at high temperatures than at room temperature.

Effects of Packing Density and Chain Length on the Surface Hydrophobicity of Thin Films Composed of Perfluoroalkyl Acrylate Chains: A Molecular Dynamics Study

A good understanding of the surface hydrophobicity of fluorinated materials is useful for their application as coating materials. The present study investigates the surface hydrophobicity of perfluoroalkyl acrylate (PFA) thin films using molecular dynamics simulations. Surface hydrophobicity is characterized by examining the contact angle of a water droplet on PFA surfaces and the cavity formation free energy in the vicinity of the surface. It is found that the calculated microscopic contact angles are in good agreement with the experimental results and partially capture the difference in the hydrophobicity of the surface arising from the variation of packing density and side chain length of PFA. The variations of cavity formation free energy in the vicinity of the surface elucidate that the surface hydrophobicity is mainly governed by the packing density rather than the chain length of PFA. The hydrophobicity generally increases with decreasing the packing density to some extent and then turns to decrease as further reducing the packing density. At higher packing density, the surface hydrophobicity slightly decreases with increasing the chain length, while at the lower packing density, the surface hydrophobicity is increased when chain length of PFA is longer than six carbons. Furthermore, we found that the influence of packing density on the surface hydrophobicity is directly related to the variation of the surface roughness and chain flexibility, that is, the surface hydrophobicity increases with increase in the surface roughness, while the chain flexibility plays a secondary role in the enhancement by affecting the stability of water staying near the interface. The study provides a significant insight into the local hydrophobicity and microscopic structure of the PFA surfaces, which would be useful for the application of surface modification.

Binding potentials for vapour nanobubbles on surfaces using density functional theory

We calculate density profiles of a simple model fluid in contact with a planar surface using density functional theory (DFT), in particular for the case where there is a vapour layer intruding between the wall and the bulk liquid. We apply the method of Hughes et al (2015 J. Chem. Phys. 142 074702) to calculate the density profiles for varying (specified) amounts of the vapour adsorbed at the wall. This is equivalent to varying the thickness h of the vapour at the surface. From the resulting sequence of density profiles we calculate the thermodynamic grand potential as h is varied and thereby determine the binding potential as a function of h. The binding potential obtained via this coarse-graining approach allows us to determine the disjoining pressure in the film and also to predict the shape of vapour nano-bubbles on the surface. Our microscopic DFT based approach captures information from length scales much smaller than some commonly used models in continuum mechanics.

Experimental Observation of Crystal-Liquid Coexistence in Slit-Confined Nonpolar Fluids

Films of carbon tetrachloride (CCl₄) confined in slit geometry between two flat diamond substrates down to a few tens of Angstroms are studied by combining X-ray reflectivity with in-plane and out-of-plane X-ray scattering. The confined films form a heterogeneous structure with coexisting regions of liquid and crystalline phases. The liquid phase shows short-range ordering normal to the surfaces of the substrates. The experiments directly show the ability of the confinement to induce crystal objects, which is a long-discussed issue in the literature. The surface structure and morphology of the substrates may influence the actual realization of the crystalline phase in confinement.

Solvent-mediated interactions between nanostructures: From water to Lennard-Jones liquid

Solvent-mediated interactions emerge from complex mechanisms that depend on the solute structure, its wetting properties, and the nature of the liquid. While numerous studies have focused on the first two influences, here, we compare the results from water and Lennard-Jones liquid in order to reveal to what extent solvent-mediated interactions are universal with respect to the nature of the liquid. Besides the influence of the liquid, the results were obtained with classical density functional theory and brute-force molecular dynamics simulations which allow us to contrast these two numerical techniques.

Alchemical free energy calculations and umbrella sampling with local molecular field theory

Understanding the thermodynamic driving forces underlying any chemical process requires a description of the underlying free energy surface. However, computation of free energies is difficult, often requiring advanced sampling techniques. Moreover, these computations can be further complicated by the evaluation of any long-ranged interactions in the system of interest, such as Coulomb interactions in charged and polar media. Local molecular field theory is a promising approach to avoid many of the conceptual and computational difficulties associated with long-ranged interactions. We present frameworks for performing alchemical free energy calculations and non-Boltzmann sampling with local molecular field theory. We demonstrate that local molecular field theory can be used to perform these free energy calculations with accuracy comparable to traditional methodologies while eliminating the need for explicit treatment of long-ranged interactions in simulations.

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.

Stabilization and Antifouling of Polymer Films on a Planar Surface by CO2 Pressurization

In this article, we study the dewetting phenomenon of a polymer and carbon dioxide (CO₂) mixture on a planar surface by combining density functional theory and the string method. It is found that dewetting is a first-order discontinuous phase transition. When the pressure is lower than the completely dewetting pressure (P_d), CO₂ stabilizes the polymer films. The density fluctuation of the polymer decreases significantly with the inclusion of CO₂. When the pressure is above P_d, the polymer film is depleted far away from the surface, leaving a thick layer of pure CO₂ in the region near the surface. P_d is proportional to the surface energy strength. The CO₂ molecules enhance the density fluctuation of the polymer during the dewetting process. The polymer-rich phase at the triple point dewets to a CO₂-rich vapor film, as the CO₂-rich liquid film near the surface is metastable. These results have promising application in the industry of fabricating polymer films and antifouling polymers on attractive surfaces.

Modulus-pressure equation for confined fluids

Ultrasonic experiments allow one to measure the elastic modulus of bulk solid or fluid samples. Recently such experiments have been carried out on fluid-saturated nanoporous glass to probe the modulus of a confined fluid. In our previous work [G. Y. Gor et al., J. Chem. Phys., 143, 194506 (2015)], using Monte Carlo simulations we showed that the elastic modulus K of a fluid confined in a mesopore is a function of the pore size. Here we focus on the modulus-pressure dependence K(P), which is linear for bulk materials, a relation known as the Tait-Murnaghan equation. Using transition-matrix Monte Carlo simulations we calculated the elastic modulus of bulk argon as a function of pressure and argon confined in silica mesopores as a function of Laplace pressure. Our calculations show that while the elastic modulus is strongly affected by confinement and temperature, the slope of the modulus versus pressure is not. Moreover, the calculated slope is in a good agreement with the reference data for bulk argon and experimental data for confined argon derived from ultrasonic experiments. We propose to use the value of the slope of K(P) to estimate the elastic moduli of an unknown porous medium.

Critical Drying of Liquids

We report a detailed simulation and classical density functional theory study of the drying transition in a realistic model fluid at a smooth substrate. This transition (in which the contact angle θ→180°) is shown to be critical for both short-ranged and long-ranged substrate-fluid interaction potentials. In the latter case critical drying occurs at exactly zero attractive substrate strength. This observation permits the accurate elucidation of the character of the transition via a finite-size scaling analysis of the density probability function. We find that the critical exponent ν_{∥} that controls the parallel correlation length, i.e., the extent of vapor bubbles at the wall, is over twice as large as predicted by mean field and renormalization group calculations. We suggest a reason for the discrepancy. Our findings shed new light on fluctuation phenomena in fluids near hydrophobic and solvophobic interfaces.

A Relation for Nanodroplet Diffusion on Smooth Surfaces

In this work, we study the diffusion of nanodroplets on smooth surfaces through molecular dynamics (MD) simulations and theoretical analyses. Molecular dynamics simulations show that nanodroplet surface diffusion is different from that of single molecules and solid particles. The dependence of nanodroplet diffusion coefficient on temperature undergoes a transition from linear to nonlinear as the surface wettability is weakened due to the coupling of temperature and surface energy. We also develop a simple relation for the diffusion coefficient by using the contact angle and contact radius of the droplet. It works well for a wide range of surface wettabilities and different sized nanodroplets, as confirmed by MD simulations.

Understanding n-Octane Behavior near Graphene with Scaled Solvent-Solute Attractions

We employ molecular dynamics simulations of n-octane near a layered graphene surface to study the related phenomena of solvation, density fluctuations, wettability, and structure and dynamics of n-octane molecules in the inhomogeneous interfacial environment. That solvation in bulk n-octane displays a lengthscale-dependent crossover similar to that of hydrophobic solvation in water is known. Here we show that, near an extended graphene interface having attractive interactions with n-octane, lengthscale-dependent solvation is similar to that in the bulk and displays a small to large crossover. However, as the n-octane-graphene interactions are reduced to make the surface increasingly solvophobic, the crossover behavior is modulated and essentially absent near the most solvophobic surfaces, similar to that in water near hydrophobic interfaces. We show that the macroscopic measure of wettability, namely, the contact angle, characterizes n-octane-graphene coupling over a limited range of attractions. In contrast, molecular measures such as the free energy of cavity formation or the local compressibility in the interfacial region provide an effective measure of this coupling over a broader range of attractions. Finally, as n-octane-graphene attractions are increased, the n-octane liquid displays a wetting transition and corresponding change from sigmoidal to layered density profile. Analysis of the local structure shows that n-octane molecules prefer approximately linear conformations and surface-parallel orientations near the graphene surface, and their translational dynamics slow down with increasing n-octane-graphene attractions. Our study highlights molecular scale behavior of n-octane molecules that is relevant to understanding nanoparticle-solvent coupling in composite materials with enhanced mechanical or thermal properties.