Chemical potentials and potential distributions of inclusion gas in quenched-annealed random porous media

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.

Fluids in porous media. IV. Quench effect on chemical potential

It appears to be a common sense to measure the crowdedness of a fluid system by the densities of the species constituting it. In the present work, we show that this ceases to be valid for confined fluids under some conditions. A quite thorough investigation is made for a hard sphere (HS) fluid adsorbed in a hard sphere matrix (a quench-annealed system) and its corresponding equilibrium binary mixture. When fluid particles are larger than matrix particles, the quench-annealed system can appear much more crowded than its corresponding equilibrium binary mixture, i.e., having a much higher fluid chemical potential, even when the density of each species is strictly the same in both systems, respectively. We believe that the insight gained from this study should be useful for the design of functionalized porous materials.

A three dimensional integral equation approach for fluids under confinement: Argon in zeolites

In this work, we explore the ability of an inhomogeneous integral equation approach to provide a full three dimensional description of simple fluids under conditions of confinement in porous media. Explicitly, we will consider the case of argon adsorbed into silicalite-1, silicalite-2, and an all-silica analogue of faujasite, with a porous structure composed of linear (and zig-zag in the case of silicalite-1) channels of 5-8 Å diameter. The equation is based on the three dimensional Ornstein-Zernike approximation proposed by Beglov and Roux [J. Chem. Phys. 103, 360 (1995)] in combination with the use of an approximate fluid-fluid direct correlation function furnished by the replica Ornstein-Zernike equation with a hypernetted chain closure. Comparison with the results of grand canonical Monte Carlo/molecular dynamics simulations evidences that the theory provides an accurate description for the three dimensional density distribution of the adsorbed fluid, both at the level of density profiles and bidimensional density maps across representative sections of the porous material. In the case of very tight confinement (silicalite-1 and silicalite-2), solutions at low temperatures could not be found due to convergence difficulties, but for faujasite, which presents substantially larger channels, temperatures as low as 77 K are accessible to the integral equation. The overall results indicate that the theoretical approximation can be an excellent tool to characterize the microscopic adsorption behavior of porous materials.

A computational study of electrolyte adsorption in a simple model for intercalated clays

A pillared interlayered clay is represented by a two-dimensional quenched charged disordered medium, in which the pillar configuration is produced by the quench of a two-dimensional electrolyte and the subsequent removal of the anions (that act as a template). The cation charge is counterbalanced by a neutralizing background that is an ideal representation of the layer's negative charge in the experimental system. In this paper we investigate the adsorption of electrolyte particles in this charged disordered medium resorting both to the use of the replica Ornstein-Zernike equation in the hypernetted chain approximation and grand canonical Monte Carlo simulations. The theoretical approach qualitatively reproduces the simulated behavior of the adsorbed fluids. Theoretical estimates of the material porosities obtained for various types of pillar distributions are in good agreement with the simulation. We investigate the influence of the matrix on correlation functions and adsorption isotherms.

Phase separation of model adsorbates in random matrices

The liquid-liquid phase behavior of binary mixtures in random pores is investigated with non-additive hard spheres using both ROZ (Replica Ornstein-Zernike) integral equations and cavity biased grand canonical Monte Carlo simulations. The critical densities of the coexistence phase envelopes are determined as function of the non-additivity parameter Delta, varying from Delta = 0.2, 0.4, 0.6, up to 0.8. The matrix is made of quenched hard spheres. Its porosity is varied to ascertain the effects of confinement, with packing densities rho(0) ranging from 0.1, 0.3, to 0.5. To obtain fiduciary results from ROZ, we use the accurate ZSEP closure relation proposed earlier with and without thermodynamic consistency. The ZSEP closure is known to enforce the zero-separation theorems via special adjustable parameters in the bridge function. Two versions of this closure are used to assess their accuracies (vis-à-vis the Monte Carlo data): first ZSEP-T, namely, the ZSEP closure with added thermodynamic consistency (the Gibbs-Duhem relation); and second purely ZSEP without adding thermodynamic consistency. It is found that both closures give correct qualitative trend, with errors of ZSEP falling within 8-9%, while ZSEP-T, being more accurate, to within 1-2%. As non-additivity is increased, both versions become more accurate. The critical density rho(c) is found to decrease with decreasing porosity. In addition, rho(c) also decreases with increasing Delta, in a non-monotone fashion.

Fluids in porous media. I. A hard sponge model

The morphology of many porous materials is spongelike. Despite the abundance of such materials, simple models which allow for a theoretical description of these materials are still lacking. Here, we propose a hard sponge model which is made by digging spherical cavities in a solid continuum. We found an analytical expression for describing the interaction potential between fluid particles and the spongelike porous matrix. The diagrammatic expansions of different correlation functions are derived as well as that of grand potential. We derived also the Ornstein-Zernike (OZ) equations for this model. In contrast to Madden-Glandt model of random porous media [W. G. Madden and E. D. Glandt, J. Stat. Phys. 51, 537 (1988)], the OZ equations for a fluid confined in our hard sponge model have some similarity to the OZ equations of a three-component fluid mixture. We show also how the replica method can be extended to study our sponge model and that the same OZ equations can be derived also from the extended replica method.

Chemical potentials and phase equilibria of Lennard-Jones mixtures: a self-consistent integral equation approach

We explore the vapor-liquid phase behavior of binary mixtures of Lennard-Jones-type molecules where one component is supercritical, given the system temperature. We apply the self-consistency approach to the Ornstein-Zernike integral equations to obtain the correlation functions. The consistency checks include not only thermodynamic consistencies (pressure consistency and Gibbs-Duhem consistency), but also pointwise consistencies, such as the zero-separation theorems on the cavity functions. The consistencies are enforced via the bridge functions in the closure which contain adjustable parameters. The full solution requires the values of not only the monomer chemical potentials, but also the dimer chemical potentials present in the zero-separation theorems. These are evaluated by the direct chemical-potential formula [L. L. Lee, J. Chem. Phys. 97, 8606 (1992)] that does not require temperature nor density integration. In order to assess the integral equation accuracy, molecular-dynamics simulations are carried out alongside the states studied. The integral equation results compare well with simulation data. In phase calculations, it is important to have pressure consistency and valid chemical potentials, since the matching of phase boundaries requires the equality of the pressures and chemical potentials of both the liquid and vapor phases. The mixtures studied are methane-type and pentane-type molecules, both characterized by effective Lennard-Jones potentials. Calculations on one isotherm show that the integral equation approach yields valid answers as compared with the experimental data of Sage and Lacey. To study vapor-liquid phase behavior, it is necessary to use consistent theories; any inconsistencies, especially in pressure, will vitiate the phase boundary calculations.

Replica Ornstein-Zernike self-consistent theory for mixtures in random pores

We present a self-consistent integral equation theory for a binary liquid in equilibrium with a disordered medium, based on the formalism of the replica Ornstein-Zernike (ROZ) equations. Specifically, we derive direct formulas for the chemical potentials and the zero-separation theorems (the latter provide a connection between the chemical potentials and the fluid cavity distribution functions). Next we solve a modified-Verlet closure to ROZ equations, which has built-in parameters that can be adjusted to satisfy the zero-separation theorems. The degree of thermodynamic consistency of the theory is also kept under control. We model the binary fluid in random pores as a symmetrical binary mixture of nonadditive hard spheres in a disordered hard-sphere matrix and consider two different values of the nonadditivity parameter and of the quenched matrix packing fraction, at different mixture concentrations. We compare the theoretical structural properties as obtained through the present approach with Percus-Yevick and Martinov-Sarkisov integral equation theories, and assess both structural and thermodynamic properties by performing canonical standard and biased grand canonical Monte Carlo simulations. Our theory appears superior to the other integral equation schemes here examined and provides reliable estimates of the chemical potentials. This feature should be useful in studying the fluid phase behavior of model adsorbates in random pores in general.

Adsorption of a diatomic molecular fluid into random porous media

Structural and thermodynamic properties of a homonuclear hard dumbbell fluid adsorbed into a disordered hard sphere matrix are studied by means of integral equation techniques and computer simulation. In particular, we have rewritten the replica Ornstein-Zernike equations to deal with orientational degrees of freedom and we have solved them in two different approaches: the hypernetted chain equation and a semiempirical extension of Verlet's approximation. We have also derived direct expressions to calculate the chemical potential in these approximations. Comparison with grand canonical Monte Carlo results shows that both theoretical treatments describe adequately the physical behavior of the system, Verlet's approach being, however, clearly superior in accordance with previous findings for equilibrated hard core mixtures.