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

Deformable hard particles confined in a disordered porous matrix

With suitably designed Monte Carlo simulations, we have investigated the properties of mobile, impenetrable, yet deformable particles that are immersed into a porous matrix, the latter one realized by a frozen configuration of spherical particles. By virtue of a model put forward by Batista and Miller [Phys. Rev. Lett. 105, 088305 (2010)], the fluid particles can change in their surroundings, formed by other fluid particles or the matrix particles, their shape within the class of ellipsoids of revolution; such a change in shape is related to a change in energy, which is fed into suitably defined selection rules in the deformation "moves" of the Monte Carlo simulations. This concept represents a simple yet powerful model of realistic, deformable molecules with complex internal structures (such as dendrimers or polymers). For the evaluation of the properties of the system, we have used the well-known quenched-annealed protocol (with its characteristic double average prescription) and have analyzed the simulation data in terms of static properties (the radial distribution function and aspect ratio distribution of the ellipsoids) and dynamic features (notably the mean squared displacement). Our data provide evidence that the degree of deformability of the fluid particles has a distinct impact on the aforementioned properties of the system.

Prediction of Adsorption Isotherms and Selectivities: Comparison between Classical Density Functional Theory Based on the Perturbed-Chain Statistical Associating Fluid Theory Equation of State and Ideal Adsorbed Solution Theory

This study gives an assessment of the predictive capability of classical density functional theory (DFT) for adsorption processes of pure substances and mixtures of spherical and nonspherical molecular species. A Helmholtz energy functional based on the perturbed-chain statistical associating fluid theory (PC-SAFT) is applied to calculate isotherms and selectivities of multicomponent adsorption. In order to unambiguously assess the accuracy of the DFT model, we conduct molecular simulations. Monte Carlo (MC) simulations are performed in the grand canonical ensemble using the transition matrix. Two types of systems are studied: a model system, where fluid-fluid and solid-fluid interactions are defined as (single-site) Lennard-Jones interactions, and a more realistic methane-n-butane mixture in a graphite-like pore. Differences between a slit-shaped and a cylindrical pore geometry are examined for the model system. Adsorption isotherms and selectivities obtained from DFT calculations and MC simulations are found in very good agreement, particularly at high pressures. Capillary condensation observed along adsorption isotherms containing n-butane was accurately predicted, both, in equilibrium pressure and in density-increase. Comparisons with results from the ideal adsorbed solution theory are presented, confirming powerful predictions of the DFT approach.

Phase Equilibria of Polydisperse Square-Well Chain Fluid Confined in Random Porous Media: TPT of Wertheim and Scaled Particle Theory

Extension of Wertheim's thermodynamic perturbation theory and its combination with scaled particle theory is proposed and applied to study the liquid-gas phase behavior of polydisperse hard-sphere square-well chain fluid confined in the random porous media. Thermodynamic properties of the reference system, represented by the hard-sphere square-well fluid in the matrix, are calculated using corresponding extension of the second-order Barker-Henderson perturbation theory. We study effects of polydispersity and confinement on the phase behavior of the system. While polydispersity causes increase of the region of phase coexistence due to the critical temperature increase, confinement decreases the values of both critical temperature and critical density making the region of phase coexistence smaller. This effect is enhanced with the increase of the size ratio of the fluid and matrix particles. The increase of the average chain length at fixed values of polydispersity and matrix density shifts the critical point to a higher temperature and a slightly lower density.

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.

The Structural Properties and Diffusion of a Three-Dimensional Isotropic Core-Softened Model Fluid in Disordered Porous Media. Molecular Dynamics Simulation

The microscopic structure and dynamic properties of an isotropic three-dimensional core-softened model fluid in disordered matrices of Lennard-Jones particles have been studied. Molecular dynamics computer simulations in Grand Canonical ensemble were used as the methodological tools. It was shown that the microscopic structure of the fluid is characterized by anomalies similar to those found in a bulk model, but that it is affected by the fluid-matrix interactions. The dynamic properties also exhibit anomalous dependence on fluid density, but the magnitude of these anomalies is suppressed in comparison to the bulk fluid model. The anomalous behaviour of the diffusion coefficient is attributed to structural changes in the first coordination shell of a given fluid particle. It seems that the anomalies can only be suppressed at matrix densities which are higher than those studied in the present work.

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.

Effect of disorder on condensation in the lattice gas model on a random graph

The lattice gas model of condensation in a heterogeneous pore system, represented by a random graph of cells, is studied using an exact analytical solution. A binary mixture of pore cells with different coordination numbers is shown to exhibit two phase transitions as a function of chemical potential in a certain temperature range. Heterogeneity in interaction strengths is demonstrated to reduce the critical temperature and, for large-enough degreeS of disorder, divides the cells into ones which are either on average occupied or unoccupied. Despite treating the pore space loops in a simplified manner, the random-graph model provides a good description of condensation in porous structures containing loops. This is illustrated by considering capillary condensation in a structural model of mesoporous silica SBA-15.

Fluids in porous media: the case of neutral walls

The bulk phase behavior of a fluid is typically altered when the fluid is brought into confinement by the walls of a random porous medium. Inside the porous medium, phase-transition points are shifted, or may disappear altogether. A crucial determinant is how the walls interact with the fluid particles. In this work, we consider the situation whereby the walls are neutral with respect to the liquid and vapor phases. In order to realize the condition of strict neutrality, we use a symmetric binary mixture inside a porous medium that interacts identically with mixture species. Monte Carlo simulations are then used to obtain the phase behavior. Our main finding is that, in the presence of the porous medium, a liquid-vapor critical point still exists. At the critical point, the distribution of the order parameter remains scale invariant, but self-averaging is violated. These findings provide further evidence that random confinement by neutral walls induces critical behavior of the random Ising model (i.e., Ising models with dilution type disorder, where the disorder couples to the energy).

Dynamic arrest of colloids in porous environments: disentangling crowding and confinement

Using numerical simulations we study the slow dynamics of a colloidal hard-sphere fluid adsorbed in a matrix of disordered hard-sphere obstacles. We calculate separately the contributions to the single-particle dynamic correlation functions due to free and trapped particles. The separation is based on a Delaunay tessellation to partition the space accessible to the centres of fluid particles into percolating and disconnected voids. We find that the trapping of particles into disconnected voids of the matrix is responsible for the appearance of a nonzero long-time plateau in the single-particle intermediate scattering functions of the full fluid. The subdiffusive exponent z, obtained from the logarithmic derivative of the mean squared displacement, is essentially unaffected by the motion of trapped particles: close to the percolation transition, we determined z approximately = 0.5 for both the full fluid and the particles moving in the percolating void. Notably, the same value of z is found in single-file diffusion and is also predicted by mode-coupling theory along the diffusion-localization line. We also reveal subtle effects of dynamic heterogeneity in both the free and the trapped component of the fluid particles, and discuss microscopic mechanisms that contribute to this phenomenon.

Finite-size scaling in Ising-like systems with quenched random fields: evidence of hyperscaling violation

In systems belonging to the universality class of the random field Ising model, the standard hyperscaling relation between critical exponents does not hold, but is replaced with a modified hyperscaling relation. As a result, standard formulations of finite-size scaling near critical points break down. In this work, the consequences of modified hyperscaling are analyzed in detail. The most striking outcome is that the free-energy cost ΔF of interface formation at the critical point is no longer a universal constant, but instead increases as a power law with system size, ΔF∝L(θ), with θ as the violation of hyperscaling critical exponent and L as the linear extension of the system. This modified behavior facilitates a number of numerical approaches that can be used to locate critical points in random field systems from finite-size simulation data. We test and confirm the approaches on two random field systems in three dimensions, namely, the random field Ising model and the demixing transition in the Widom-Rowlinson fluid with quenched obstacles.

A Monte Carlo study of crowding effects on the self-assembly of amphiphilic molecules

In this work, lattice Monte Carlo was used to study the effects of crowding on the self-assembly of surfactants. Simulation results show that crowding strongly shifts the critical micelle concentration (CMC) of surfactants from the bulk value. Two effects originated from crowding are found to govern the CMC shift: one is the depletion effect by crowding agents and the other is the available volume for micelle formation. The depletion effects inevitably result in the enrichment of surfactants in crowding-free regions and cause the decrease in CMC. On the other hand, the appearance of crowding agents decreases the available volume for micelle formation, which reduces the conformational entropy and impedes the micelle formation. Three factors, including the radius of crowding agents, the arrangement of crowding agents, and the volume fraction of crowding agents, are considered in this work to study the crowding effects. The trends of CMC shifts are interpreted from the competition between the depletion effects and the available volume for micelle formation.

Confinement effects on phase behavior of soft matter systems

When systems that can undergo phase separation between two coexisting phases in the bulk are confined in thin film geometry between parallel walls, the phase behavior can be profoundly modified. These phenomena shall be described and exemplified by computer simulations of the Asakura-Oosawa model for colloid-polymer mixtures, but applications to other soft matter systems (e.g. confined polymer blends) will also be mentioned. Typically a wall will prefer one of the phases, and hence the composition of the system in the direction perpendicular to the walls will not be homogeneous. If both walls are of the same kind, this effect leads to a distortion of the phase diagram of the system in thin film geometry, in comparison with the bulk, analogous to the phenomenon of "capillary condensation" of simple fluids in thin capillaries. In the case of "competing walls", where both walls prefer different phases of the two phases coexisting in the bulk, a state with an interface parallel to the walls gets stabilized. The transition from the disordered phase to this "soft mode phase" is rounded by the finite thickness of the film and is not a sharp phase transition. However, a sharp transition can occur where this interface gets localized at (one of) the walls. The relation of this interface localization transition to wetting phenomena is discussed. Finally, an outlook to related phenomena is given, such as the effects of confinement in cylindrical pores on the phase behavior, and more complicated ordering phenomena (lamellar mesophases of block copolymers or nematic phases of liquid crystals under confinement).

Adsorption of fluids in a pore with chemical heterogeneities: the cooperative effect

In this work, we study the cooperative adsorption of fluids in a heterogeneous pore, in which the pore walls are composed of homogeneous substrates with chemical groups (CGs) decorating them. The adsorption caused by the homogeneous substrates alone and that by CGs do not add up to the overall adsorption, indicating the existence of a cooperative effect. The cooperative effect is the source of cooperative adsorption, and is characterized in this work by the ratio of the overall adsorption to the sum of adsorption by the substrate only and that by CGs. It is found that the cooperative adsorption does not depend monotonically on the substrate or the CGs. Two different origins of the cooperative adsorption play different roles depending on which one dominates the overall adsorption. Our simulations reveal that, when the homogeneous substrate dominates the overall adsorption, weakening of the attractive fluid-substrate interaction or alternatively strengthening of the fluid-CGs interaction leads to a stronger cooperative effect and enhances the cooperative adsorption. However, when CGs dominate the overall adsorption, weakening of the attractive fluid-CG interaction or strengthening the fluid-substrate interaction results in strong cooperative adsorption. In order to investigate the effects of the distribution of CGs on cooperative adsorption, a design-test method is generalized and used in this work. Simulation results show that the overall adsorption can be significantly affected by the CG distribution.

Colloid-polymer mixtures in the presence of quenched disorder: a theoretical and computer simulation study

We use theory and computer simulation to study the structure and phase behavior of colloid-polymer mixtures in the presence of quenched disorder. The Asakura-Oosawa model (AO) (Asakura and Oosawa 1954 J. Chem. Phys. 22 1255) is used to describe the colloid-colloid, colloid-polymer, and polymer-polymer pair interactions. We then investigate the behavior of this model in the presence of frozen-in (quenched) obstacles. The obstacles will be placed according to two different scenarios, both of which are experimentally feasible. In the first scenario, polymers are distributed at positions drawn from an ideal gas configuration. In the second scenario, colloidal particles are distributed at positions drawn from an equilibrium hard sphere configuration. We investigate how the unmixing transition of the AO model is affected by the type of quenched disorder. The theoretical formalism is based on the replica method of Given and Stell (1994 Physica A 209 495). Our foremost aim is to test the accuracy of three common closures to the replica Ornstein-Zernike equations, namely the hypernetted chain, the Percus-Yevick, and the Martinov-Sarkisov equations. The accuracy is determined by comparison with grand canonical Monte Carlo simulations. We find that, for quenched polymer disorder, all three closures perform remarkably well. However, when quenched colloid disorder is considered, i.e. the second mentioned scenario, the predictions of all three closures worsen dramatically.

Thermodynamics of reversible gas adsorption on alkali-metal exchanged zeolites—the interplay of infrared spectroscopy and theoretical calculations

Detailed understanding of weak solid-gas interactions giving rise to reversible gas adsorption on zeolites and related materials is relevant to both, fundamental studies on gas adsorption and potential improvement on a number of (adsorption based) technological processes. Combination of variable-temperature infrared spectroscopy with theoretical calculations constitutes a fruitful approach towards both of these aims. Such an approach is demonstrated here (mainly) by reviewing recent studies on hydrogen and carbon monoxide adsorption (at a low temperature) on alkali-metal exchanged ferrierite. However, the methodology discussed, which involves the interplay of experimental measurements and theoretical calculations at the periodic DFT level, should be equally valid for many other gas-solid systems. Specific aspects considered are the identification of gas adsorption complexes and thermodynamic studies related to standard adsorption enthalpy and entropy.

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.

Critical behavior of colloid-polymer mixtures in random porous media

We show that the critical behavior of a colloid-polymer mixture inside a random porous matrix of quenched hard spheres belongs to the universality class of the random-field Ising model. We also demonstrate that random-field effects in colloid-polymer mixtures are surprisingly strong. This makes these systems attractive candidates to study random-field behavior experimentally.

Capillary phase transitions of linear and branched alkanes in carbon nanotubes from molecular simulation

Capillary phase transitions of linear (from C(1) to C(12)) and branched (C(5) isomers) alkanes in single-walled carbon nanotubes have been investigated using the gauge-cell Monte Carlo simulation. The isotherm at a supercritical temperature increases monotonically with chemical potential and coincides with that from the traditional grand canonical Monte Carlo simulation, whereas the isotherm at a subcritical temperature exhibits a sigmoid van der Waals loop including stable, metastable, and unstable regions. Along this loop, the coexisting phases are determined using an Maxwell equal-area construction. A generic confinement effect is found that reduces the saturation chemical potential, lowers the critical temperature, increases the critical density, and shrinks the phase envelope. The effect is greater in a smaller diameter nanotube and is greater in a nanotube than in a nanoslit.

Phase coexistence in heterogeneous porous media: A new extension to Gibbs ensemble Monte Carlo simulation method

The effect of confinement on phase behavior of simple fluids is still an area of intensive research. In between experiment and theory, molecular simulation is a powerful tool to study the effect of confinement in realistic porous materials, containing some disorder. Previous simulation works aiming at establishing the phase diagram of a confined Lennard-Jones-type fluid, concentrated on simple pore geometries (slits or cylinders). The development of the Gibbs ensemble Monte Carlo technique by Panagiotopoulos [Mol. Phys. 61, 813 (1987)], greatly favored the study of such simple geometries for two reasons. First, the technique is very efficient to calculate the phase diagram, since each run (at a given temperature) converges directly to an equilibrium between a gaslike and a liquidlike phase. Second, due to volume exchange procedure between the two phases, at least one invariant direction of space is required for applicability of this method, which is the case for slits or cylinders. Generally, the introduction of some disorder in such simple pores breaks the initial invariance in one of the space directions and prevents to work in the Gibbs ensemble. The simulation techniques for such disordered systems are numerous (grand canonical Monte Carlo, molecular dynamics, histogram reweighting, N-P-T+test method, Gibbs-Duhem integration procedure, etc.). However, the Gibbs ensemble technique, which gives directly the coexistence between phases, was never generalized to such systems. In this work, we focus on two weakly disordered pores for which a modified Gibbs ensemble Monte Carlo technique can be applied. One of the pores is geometrically undulated, whereas the second is cylindrical but presents a chemical variation which gives rise to a modulation of the wall potential. In the first case almost no change in the phase diagram is observed, whereas in the second strong modifications are reported.

Thermodynamic Characterization of Fluids Confined in Heterogeneous Pores by Monte Carlo Simulations in the Grand Canonical and the Isobaric−Isothermal Ensembles

Materials presenting nanoscale porosity are able to condense gases in their structure. This "capillary condensation" phenomenon has been studied for more than one century. Theoretical models help to understand experimental results but fail in explaining all experimental features. Most of the time, the difficulties in making quantitative or even qualitative predictions are due to the geometric complexity of the porous materials, such as large pore size distribution, chemical heterogeneities, or pore interconnections. Numerical calculations (lattice gas models or molecular simulations) are of considerable interest to calculate the adsorption properties of a fluid confined in a porous model with characteristic sizes up to several tens of nanometers. For instance, the grand canonical Monte Carlo method allows one to compute the average amount of fluid adsorbed in the porous model as a function of the temperature and the chemical potential of the fluid. However, the grand potential, necessary for a complete characterization of the system, is not a direct output of the algorithm. It is shown in this paper that the use of the isobaric-isothermal (NPT) ensemble allows one to circumvent this problem; that is, it is possible to get in one single Monte Carlo run the absolute grand potential for any given thermodynamic state of the fluid. A simplified thermodynamic integration scheme is then used to evaluate the grand potential over the whole isotherm branch passing through this initially given point. Since the usual NPT technique is a priori limited to homogeneous pores, it is proposed, for the first time, to generalize this procedure to a pore presenting a chemical heterogeneity along its axis. The new method gives the same results as the previous for homogeneous pores and allows new predictions for chemically heterogeneous pores. Comparison with the full integration scheme shows that the proposed direct calculation is faster since it avoids multiple Monte Carlo runs and more precise because it avoids the possible cumulative errors of the integration procedure.

Capillary Condensation in a Geometrically and a Chemically Heterogeneous Pore: A Molecular Simulation Study

A computer simulation study has been carried out, using an extended Gibbs ensemble Monte Carlo technique, to examine the influence of so-called geometric and chemical disorder on the thermodynamic behavior of simple fluids confined in porous media. The technique allows the equilibrium coexistence of gas and liquid phases to be calculated in a single run. The phase diagram of Lennard-Jones fluid has been calculated in a perfectly cylindrical pore as a reference. Some disorder is then introduced in the porous material, first by spatially modifying the external potential of the initially cylindrical pore, to imitate the geometric disorder of a more realistic pore (undulation, constrictions, etc.) and second by modulating the amplitude of the same initially cylindrical potential to reproduce the energetic disorder of realistic pores due to chemical variations along it. It is shown that the chemical disorder has a much stronger effect on the phase diagram of the confined fluid. The complete adsorption/desorption isotherms are also calculated to help in understanding the large effects of chemical disorder.

Phase diagram of a binary symmetric hard-core Yukawa mixture

We assess the accuracy of the self-consistent Ornstein-Zernike approximation for a binary symmetric hard-core Yukawa mixture by comparison with Monte Carlo simulations of the phase diagrams obtained for different choices of the ratio alpha of the unlike-to-like interactions. In particular, from the results obtained at alpha=0.75 we find evidence for a critical endpoint in contrast to recent studies based on integral equation and hierarchical reference theories. The variation of the phase diagrams with range of the Yukawa potential is investigated.

Isotropic-nematic transition of hard rods immersed in random sphere matrices

Using replica density functional theory and Monte Carlo computer simulations we investigate a system of annealed hard spherocylinders adsorbed in a matrix of quenched hard spheres. Theoretical predictions for the partition coefficient, defined as the ratio of density of rods in the matrix and that in a reservoir, agree well with simulation results. Theory predicts the isotropic-nematic transition to remain first order upon increasing sphere packing fraction, and to shift towards lower rod densities. This scenario is consistent with our simulation results that clearly show a jump in the nematic order parameter upon increasing the rod density at constant matrix packing fraction, corresponding to the isotropic-nematic transition, even for sphere matrix packing fractions < or approximately equal to 0.3.

Grand Potential, Helmholtz Free Energy, and Entropy Calculation in Heterogeneous Cylindrical Pores by the Grand Canonical Monte Carlo Simulation Method

The adsorption of fluids in porous media is still an open area of research, since no model is able to explain all experimental features. The difficulties rise from the complexity of the real porous materials which present surface heterogeneities, large pore size distributions, and complex networks of interconnected pores. In parallel to experimental efforts trying to produce more ordered porous materials, theoreticians try to introduce more disorder in their models, with the help of molecular simulation for instance. This grand canonical Monte Carlo simulation study concentrates on the adsorption of a simple Lennard-Jones fluid in three porous substrates, to compare the effect of purely geometric heterogeneity (spatial deformation of the external potential) as opposed to purely chemical heterogeneity (amplitude variations of the external potential). This separation is unrealistic, since geometric fluctuations of a real pore diameter along its axis generally induce variations in the amplitude of the external potential created by the pore. However it enables one to compare both effects. In this paper, a thermodynamic integration scheme is applied to a complete set of adsorption/desorption isotherms. The grand potential, free energy, and entropy are calculated, which allows one to discuss the features of the phase diagrams. It is shown that a purely geometric deformation (undulation) of the external potential does not affect the thermodynamic characteristics of the confined fluid. On the other hand, amplitude modulation of the external potential (chemical heterogeneity) strongly distorts the phase diagram. This heterogeneity is actually able to stabilize a "bridgelike" phase which corresponds to an accumulation of molecules in the most attractive region of the pore.

Computer simulation of the phase diagram for a fluid confined in a fractal and disordered porous material

We present a grand canonical Monte Carlo simulation study of the phase diagram of a Lennard-Jones fluid adsorbed in a fractal and highly porous aerogel. The gel environment is generated from an off-lattice diffusion limited cluster-cluster aggregation process. Simulations have been performed with the multicanonical ensemble sampling technique. The biased sampling function has been obtained by histogram reweighting calculations. Comparing the confined and the bulk system liquid-vapor coexistence curves we observe a decrease of both the critical temperature and density in qualitative agreement with experiments and other Monte Carlo studies on Lennard-Jones fluids confined in random matrices of spheres. At variance with these numerical studies we do not observe upon confinement a peak on the liquid side of the coexistence curve associated with a liquid-liquid phase coexistence. In our case only a shouldering of the coexistence curve appears upon confinement. This shoulder can be associated with high density fluctuations in the liquid phase. The coexisting vapor and liquid phases in our system show a high degree of spatial disorder and inhomogeneity.

Vapor-liquid transitions of dipolar fluids in disordered porous media: Performance of angle-averaged potentials

Using replica integral equations in the reference hypernetted-chain (RHNC) approximation we calculate vapor-liquid spinodals, chemical potentials, and compressibilities of fluids with angle-averaged dipolar interactions adsorbed to various disordered porous media. Comparison with previous RHNC results for systems with true angle-dependent Stockmayer (dipolar plus Lennard-Jones) interactions indicate that, for a dilute hard sphere matrix, the angle-averaged fluid-fluid (ff) potential is a reasonable alternative for reduced fluid dipole moments m( *2)=mu(2)/(epsilon(0)sigma(3))< or =2.0. This range is comparable to that estimated in bulk fluids, for which RHNC results are presented as well. Finally, results for weakly polar matrices suggest that angle-averaged fluid-matrix (fm) interactions can reproduce main features observed for true dipolar (fm) interactions such as the shift of the vapor-liquid spinodals towards lower temperatures and higher densities. However, the effective attraction induced by dipolar (fm) interaction is underestimated rather than overestimated as in the case of angle-averaged ff interactions.

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.

Water in nanopores. I. Coexistence curves from Gibbs ensemble Monte Carlo simulations

Coexistence curves of water in cylindrical and slitlike nanopores of different size and water-substrate interaction strength were simulated in the Gibbs ensemble. The two-phase coexistence regions cover a wide range of pore filling level and temperature, including ambient temperature. Five different kinds of two-phase coexistence are observed. A single liquid-vapor coexistence is observed in hydrophobic and moderately hydrophilic pores. Surface transitions split from the main liquid-vapor coexistence region, when the water-substrate interaction becomes comparable or stronger than the water-water pair interaction. In this case prewetting, one and two layering transitions were observed. The critical temperature of the first layering transition decreases with strengthening water-substrate interaction towards the critical temperature expected for two-dimensional systems and is not sensitive to the variation of pore size and shape. Liquid-vapor phase transition in a pore with a wall which is already covered with two water layers is most typical for hydrophilic pores. The critical temperature of this transition is very sensitive to the pore size, in contrast to the liquid-vapor critical temperature in hydrophobic pores. The observed rich phase behavior of water in pores evidences that the knowledge of coexistence curves is of crucial importance for the analysis of experimental results and a prerequiste of meaningful simulations.

Model fluid in a porous medium: results for a Bethe lattice

We consider a lattice gas with quenched impurities or "quenched-annealed binary mixture" on the Bethe lattice. The quenched part represents a porous matrix in which the (annealed) lattice gas resides. This model features the three main factors of fluids in random porous media: wetting, randomness, and confinement. The recursive character of the Bethe lattice enables an exact treatment, whose key ingredient is an integral equation yielding the one-particle effective field distribution. Our analysis shows that this distribution consists of two essentially different parts. The first one is a continuous spectrum and corresponds to the macroscopic volume accessible to the fluid, the second is discrete and comes from finite closed cavities in the porous medium. Those closed cavities are in equilibrium with the bulk fluid within the grand canonical ensemble we use, but are inaccessible in real experimental situations. Fortunately, we are able to isolate their contributions. Separation of the discrete spectrum facilitates also the numerical solution of the main equation. The numerical calculations show that the continuous spectrum becomes more and more rough as the temperature decreases, and this limits the accuracy of the solution at low temperatures.

Molecular simulation of the vapor-liquid phase behavior of Lennard-Jones mixtures in porous solids

We present vapor-liquid phase coexistence curves for binary fluid mixtures in a disordered porous solid. The porous material is modeled as a collection of randomly dispersed hard spheres. A variant of the Monte Carlo Gibbs ensemble method [J. K. Brennan and W. Dong, J. Chem. Phys. 116, 8948 (2002)] is used to simulate Lennard-Jones fluid mixtures at several porosities: 0.9, 0.95, and 0.975. Effects based on the size and the energetics of the mixture components are studied. Pressure-composition and pressure-density phase diagrams at reduced temperatures of 0.75 and 1.0 are reported. Compared to the bulk fluid behavior, dramatic shifts in the phase envelope were found for even highly porous structures. Both the Lennard-Jones size and energy mixture parameters were found to strongly influence the resulting shape of the phase envelope.