Finite-size scaling study of the vapor-liquid critical properties of confined fluids: Crossover from three dimensions to two dimensions

Role of site-site interaction on the phase equilibria of multiple-site associating fluids in a functionalized slit pore

Vapor-liquid phase equilibria for multiple sites associating fluids with different associating strengths are investigated in a slit pore using grand-canonical transition matrix Monte Carlo method. The increase of critical temperature from two-site to four-site associating fluids at constant site strength is quite significant as compared to that of the one-site to two-site associating fluids, which is more pronounced at higher associating strength (ϵ* = 6). Monomer fraction and cluster size distribution are used to investigate the association of fluid particles in coexistence phases. The monomer fraction for both phases decreases with increased associating sites on the fluid particles due to more site-site interaction with neighboring fluid particles and forming a larger cluster. Therefore, the number of associating sites and their distribution play a vital role in the association of fluid particles. Moreover, the saturation chemical potential changes with the arrangement of the sites. For two-site associating fluids, we observe early vapor-liquid transition when the sites are oppositely placed, and when the sites are placed at 90°, the vapor-liquid transition is observed at the higher chemical potential. Moreover, four-site associating fluids with a square arrangement show early vapor-liquid phase transition, mainly because these arrangements of sites effectively interact with surface sites and the molecules in the next layer.

Confined Quantum Hard Spheres

We present computer simulation and theoretical results for a system of N Quantum Hard Spheres (QHS) particles of diameter σ and mass m at temperature T, confined between parallel hard walls separated by a distance Hσ, within the range 1≤H≤∞. Semiclassical Monte Carlo computer simulations were performed adapted to a confined space, considering effects in terms of the density of particles ρ*=N/V, where V is the accessible volume, the inverse length H−1 and the de Broglie’s thermal wavelength λB=h/2πmkT, where k and h are the Boltzmann’s and Planck’s constants, respectively. For the case of extreme and maximum confinement, 0.5<H−1<1 and H−1=1, respectively, analytical results can be given based on an extension for quantum systems of the Helmholtz free energies for the corresponding classical systems.

Phase Separation and Multibody Effects in Three-Dimensional Active Brownian Particles

Simulation studies of the phase diagram of repulsive active Brownian particles in three dimensions reveal that the region of motility-induced phase separation between a high and low density phase is enclosed by a region of gas-crystal phase separation. Near-critical loci and structural crossovers can additionally be identified in analogy with simple fluids. Motivated by the striking similarity to the behavior of equilibrium fluids with short-ranged pairwise attractions, we show that a direct mapping to pair potentials in the dilute limit implies interactions that are insufficiently attractive to engender phase separation. Instead, this is driven by the emergence of multibody effects associated with particle caging that occurs at sufficiently high number density. We quantify these effects via information-theoretical measures of n-body effective interactions extracted from the configurational structure.

Effects of Confinement on the Dielectric Response of Water Extends up to Mesoscale Dimensions

The extent of confinement effects on water is not clear in the literature. While some properties are affected only within a few nanometers from the wall surface, others are affected over long length scales, but the range is not clear. In this work, we have examined the dielectric response of confined water under the influence of external electric fields along with the dipolar fluctuations at equilibrium. The confinement induces a strong anisotropic effect which is evident up to 100 nm channel width, and may extend to macroscopic dimensions. The root-mean-square fluctuations of the total orientational dipole moment in the direction perpendicular to the surfaces is 1 order of magnitude smaller than the value attained in the parallel direction and is independent of the channel width. Consequently, the isotropic condition is unlikely to be recovered until the channel width reaches macroscopic dimensions. Consistent with dipole moment fluctuations, the effect of confinement on the dielectric response also persists up to channel widths considerably beyond 100 nm. When an electric field is applied in the perpendicular direction, the orientational relaxation is 3 orders of magnitude faster than the dipolar relaxation in the parallel direction and independent of temperature.

Influence of disordered porous media on the anomalous properties of a simple water model

The thermodynamic, dynamic, and structural behavior of a water-like system confined in a matrix is analyzed for increasing confining geometries. The liquid is modeled by a two-dimensional associating lattice gas model that exhibits density and diffusion anomalies, similar to the anomalies present in liquid water. The matrix is a triangular lattice in which fixed obstacles impose restrictions to the occupation of the particles. We show that obstacles shorten all lines, including the phase coexistence, the critical and the anomalous lines. The inclusion of a very dense matrix not only suppresses the anomalies but also the liquid-liquid critical point.

Computer simulation of liquid-vapor coexistence of confined quantum fluids

The liquid-vapor coexistence (LV) of bulk and confined quantum fluids has been studied by Monte Carlo computer simulation for particles interacting via a semiclassical effective pair potential Veff(r) = VLJ + VQ, where VLJ is the Lennard-Jones 12-6 potential (LJ) and VQ is the first-order Wigner-Kirkwood (WK-1) quantum potential, that depends on β = 1∕kT and de Boer's quantumness parameter Λ=h/σ√mε, where k and h are the Boltzmann's and Planck's constants, respectively, m is the particle's mass, T is the temperature of the system, and σ and ε are the LJ potential parameters. The non-conformal properties of the system of particles interacting via the effective pair potential Veff(r) are due to Λ, since the LV phase diagram is modified by varying Λ. We found that the WK-1 system gives an accurate description of the LV coexistence for bulk phases of several quantum fluids, obtained by the Gibbs Ensemble Monte Carlo method (GEMC). Confinement effects were introduced using the Canonical Ensemble (NVT) to simulate quantum fluids contained within parallel hard walls separated by a distance Lp, within the range 2σ ≤ Lp ≤ 6σ. The critical temperature of the system is reduced by decreasing Lp and increasing Λ, and the liquid-vapor transition is not longer observed for Lp∕σ < 2, in contrast to what has been observed for the classical system.

Critical behavior of a water monolayer under hydrophobic confinement

The properties of water can have a strong dependence on the confinement. Here, we consider a water monolayer nanoconfined between hydrophobic parallel walls under conditions that prevent its crystallization. We investigate, by simulations of a many-body coarse-grained water model, how the properties of the liquid are affected by the confinement. We show, by studying the response functions and the correlation length and by performing finite-size scaling of the appropriate order parameter, that at low temperature the monolayer undergoes a liquid-liquid phase transition ending in a critical point in the universality class of the two-dimensional (2D) Ising model. Surprisingly, by reducing the linear size L of the walls, keeping the walls separation h constant, we find a 2D-3D crossover for the universality class of the liquid-liquid critical point for L=h ≃ 50, i.e. for a monolayer thickness that is small compared to its extension. This result is drastically different from what is reported for simple liquids, where the crossover occurs for L=h ≃ 5, and is consistent with experimental results and atomistic simulations. We shed light on these findings showing that they are a consequence of the strong cooperativity and the low coordination number of the hydrogen bond network that characterizes water.

Fluids in extreme confinement

For extremely confined fluids with a two-dimensional density n in slit geometry of an accessible width L, we prove that in the limit L → 0, the lateral and transversal degrees of freedom decouple, and the latter become ideal-gas-like. For a small wall separation, the transverse degrees of freedom can be integrated out and renormalize the interaction potential. We identify nL(2) as the hidden smallness parameter of the confinement problem and evaluate the effective two-body potential analytically, which allows calculating the leading correction to the free energy exactly. Explicitly, we map a fluid of hard spheres in extreme confinement onto a 2D fluid of disks with an effective hard-core diameter and a soft boundary layer. Two-dimensional phase transitions are robust and the transition point experiences a shift O(nL(2)).

Dimensional crossover in fluids under nanometer-scale confinement

Several earlier studies have shown signatures of crossover in various static and dynamics properties of a confined fluid when the confining dimension decreases to about a nanometer. The density fluctuations govern the majority of such properties of a fluid. Here, we illustrate the crossover in density fluctuation in a confined fluid, to provide a generic understanding of confinement-induced crossover of fluid properties, using computer simulations. The crossover can be understood as a manifestation of changes in the long-wavelength behavior of fluctuation in density due to geometrical constraints. We further show that the confining potential significantly affects the crossover behavior.

Widom line and noise-power spectral analysis of a supercritical fluid

We have performed extensive molecular dynamics simulations to study noise-power spectra of density and potential energy fluctuations of a Lennard-Jones model of a fluid in the supercritical region. Emanating from the liquid-vapor critical point, there is a locus of isobaric specific heat maxima, called the Widom line, which is often regarded as an extension of the liquid-vapor coexistence line. Our simulation results show that the noise-power spectrum of the density fluctuations on the Widom line of the liquid-vapor transition exhibits three distinct 1/f^{γ} behaviors with exponents γ=0, 1.2, and 2, depending on the frequency f. We find that the intermediate frequency region with an exponent γ∼ 1 appears as the temperature approaches the Widom temperature from above or below. On the other hand, we do not find three distinct regions of 1/f^{γ} in the power spectrum of the potential energy fluctuations on the Widom line. Furthermore, we find that the power spectra of both the density and potential energy fluctuations at low frequency have a maximum on the Widom line, suggesting that the noise power can provide an alternative signature of the Widom line.