Implications of interface conventions for morphometric thermodynamics

Inside and out: surface thermodynamics from positive to negative curvature

To explore the curvature dependence of solid-fluid interfacial thermodynamics, we calculate, using Grand Canonical Monte Carlo simulation, the surface free energy γ for a 2-d hard-disk fluid confined in a circular hard container of radius R as a function of the bulk packing fraction η and wall curvature C = −1/ R. (The curvature is negative because the surface is concave) Combining this with our previous data [J. Phys. Chem. B, 124, 7938 (2020)] for the positive curvature case (a hard-disk fluid at a circular wall, C = +1/ R), we obtain a full picture of surface thermodynamics in this system over the full range of positive and negative wall curvatures. Our results show that γ is linear in C with a slope that is the same for both positive and negative wall curvatures, with deviations seen only at high negative curvatures (strong confinement) and high density. This observation indicates that the surface thermodynamics of this system is consistent with the predictions of so-called Morphometric Thermodynamics at both positive and negative curvatures. In addition, we show that classical Density Functional Theory and a generalized scaled particle theory can be constructed that give excellent agreement with the simulation data over most of the range of curvatures and densities. For extremely high curvatures, where only one or two disks can occupy the container at maximum packing, it is possible to calculate γ exactly. In this limit the simulations and DFT calculations are in remarkable agreement with the exact results.

Equation of state for confined fluids

Fluids confined in small volumes behave differently than fluids in bulk systems. For bulk systems, a compact summary of the system’s thermodynamic properties is provided by equations of state. However, there is currently a lack of successful methods to predict the thermodynamic properties of confined fluids by use of equations of state, since their thermodynamic state depends on additional parameters introduced by the enclosing surface. In this work, we present a consistent thermodynamic framework that represents an equation of state for pure, confined fluids. The total system is decomposed into a bulk phase in equilibrium with a surface phase. The equation of state is based on an existing, accurate description of the bulk fluid and uses Gibbs’ framework for surface excess properties to consistently incorporate contributions from the surface. We apply the equation of state to a Lennard-Jones spline fluid confined by a spherical surface with a Weeks–Chandler–Andersen wall-potential. The pressure and internal energy predicted from the equation of state are in good agreement with the properties obtained directly from molecular dynamics simulations. We find that when the location of the dividing surface is chosen appropriately, the properties of highly curved surfaces can be predicted from those of a planar surface. The choice of the dividing surface affects the magnitude of the surface excess properties and its curvature dependence, but the properties of the total system remain unchanged. The framework can predict the properties of confined systems with a wide range of geometries, sizes, interparticle interactions, and wall–particle interactions, and it is independent of ensemble. A targeted area of use is the prediction of thermodynamic properties in porous media, for which a possible application of the framework is elaborated.

Connect the Thermodynamics of Bulk and Confined Fluids: Confinement-Adsorption Scaling

A fluid (a gas or a liquid) adsorbed in a porous material can behave very differently from its bulk counterpart. The advent of various synthesized materials with nanopores and their wide applications have provided strong impetus for studying fluids in confinement because our current understanding is still incomplete. From a large number of Monte Carlo simulations, we found a scaling relation that allows for connecting some thermodynamic properties (chemical potential, free energy per particle, and grand potential per particle) of a confined fluid to the bulk ones. Upon rescaling the adsorbed fluid density, the adsorption isotherms for many different confining environments collapse to the corresponding bulk curve. We also reveal the intimate connection of the reported scaling relation to Gibbs theory of inhomogeneous fluids and morphological thermodynamics. The advance in our understanding of confined fluids, gained from this study, also opens attractive perspectives for circumventing experimental difficulty for directly measuring some fluid thermodynamic properties in nanoporous materials.

Local pressure for confined systems

We derive a general closed expression for the local pressure exerted onto the corrugated walls of a channel confining a fluid medium. When the fluid medium is at equilibrium, the local pressure is a functional of the shape of the walls. It is shown that, due to the intrinsic nonlocal character of the interactions among the particles forming the fluid, the applicability of approximate schemes such as the concept of a surface of tension or morphometric thermodynamics is limited to wall curvatures that are small compared to the range of particle-particle interactions.

Bending and Gaussian rigidities of confined soft spheres from second-order virial series

We use virial series to study the equilibrium properties of confined soft-spheres fluids interacting through the inverse-power potentials. The confinement is induced by hard walls with planar, spherical, and cylindrical shapes. We evaluate analytically the coefficients of order two in density of the wall-fluid surface tension γ and analyze the curvature contributions to the free energy. Emphasis is in bending and Gaussian rigidities, which are found analytically at order two in density. Their contribution to γ(R) and the accuracy of different truncation procedures to the low curvature expansion are discussed. Finally, several universal relations that apply to low-density fluids are analyzed.