Monte Carlo approach to the gas-liquid transition in porous materials

Lattice model of adsorption in disordered porous materials: mean-field density functional theory and Monte Carlo simulations

We present mean-field density functional theory calculations and Monte Carlo simulations for a lattice model of a fluid confined in a disordered porous material. The model is obtained by a coarse graining of an off-lattice model of adsorption of simple molecules in silica xerogels. In some of our calculations a model of a porous glass is also considered. The lattice models exhibit behavior that is qualitatively similar to that of their off-lattice counterparts but the computations required are much more tractable and this makes it feasible to investigate the effects of porous material microstructure at longer length scales. We focus on exploring in detail the behavior in the adsorption/desorption hysteresis region for these models. In agreement with recent results for a model that uses a random distribution of solid sites on the lattice [Kierlik et al., Phys. Rev. Lett. 87, 055701 (2001)] we show that the disorder of the solid matrix induces multiple metastable states within the hysteresis region, which are evident in both the mean-field theory calculations and the Monte Carlo simulations. These multiple metastable states can be connected by scanning curves that are very similar to those seen in experimental studies of adsorption hysteresis. The results from mean-field theory predict that while there is hysteresis in the adsorption/desorption isotherms it is not possible to locate a condition of phase equilibrium that satisfies thermodynamic consistency. A wider significance of these results is discussed.

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.

Effect of templated quenched disorder on fluid phase equilibrium

Templating offers a means to direct the structure of quenched disorder. We show here that changes in phase equilibrium due to the presence of quenched disorder can themselves be altered by templating. We calculate the phase diagram of a fluid in a collection of template-directed, quenched particles by solving a set of replica Ornstein-Zernike equations within the mean spherical approximation and show templating to enhance phase behavior, that is, shift the phase envelope upward from its location for a nontemplated system of identical available volume. This enhancement is due to an augmented number of fluid-fluid interactions.

Phase diagram of a symmetric binary fluid in a porous matrix

The phase behavior of a binary symmetric fluid in thermal equilibrium with a porous matrix has been studied with the optimized random phase approximation and grand canonical Monte Carlo simulations. Depending on the matrix properties and the matrix-fluid and fluid-fluid interactions we find three types of phase diagram characterized by a tricritical point, a tricritical point with a triple point, or a critical end point. Small changes in the properties of the matrix or in the interactions are demonstrated to lead to drastic modifications of the phase diagram of the fluid, in qualitative agreement with observations in experimental studies. We show, in particular, that the change between the different types of phase diagram is triggered not only by the fluid-fluid interactions (internal parameters) but also by the properties of the matrix and of the matrix-fluid potentials (external parameters).

Computer simulations of phase equilibrium for a fluid confined in a disordered porous structure

We present calculations of the phase diagrams of a Lennard-Jones 12-6 fluid confined in a disordered porous structure made up of a dispersion of spherical particles, following up on an earlier work on the same system. In particular we present additional calculations using more realizations of the matrix and we investigate the applicability of the Gibbs-Duhem integration method to the calculation of phase equilibrium in these systems. The essential picture of disordered and inhomogeneous coexisting vapor and liquid phases, which emerged in the earlier work, is confirmed by the new calculations. However, a second phase transition associated with the wetting of the porous material by the fluid is found to be more sensitive to variations of the matrix realization. While for the present model this transition appears for particular realizations of the matrix, it does not seem to survive averaging over realizations.