Capillary Condensation in Pores with Energetically Heterogeneous Walls: Density Functional versus Monte Carlo Calculations

2D-NLDFT adsorption models for porous oxides with corrugated cylindrical pores

In this work, we develop two-dimensional models based on the non-local density functional theory (2D-NLDFT) for the analysis of pore size distribution (PSD) of oxide materials with cylindrical pores with rough and heterogeneous walls. The existing standard NLDFT models for porous oxides assume the smooth energetically uniform surface of the pore walls. Due to this assumption, the calculated theoretical isotherms show typical layering transitions, which are not consistent with the experimental adsorption isotherms measured on real oxide materials. As a result, the fits of standard NLDFT models to N₂ or Ar adsorption isotherms show deviations from the experimental points in association with artifacts observed on the calculated PSD plots. The 2D-NLDFT framework allows us to improve the standard model by introducing the corrugation and energetic heterogeneity to the surface of cylindrical pores. The surface roughness and energetic heterogeneity are known characteristics of the oxide surfaces. With these assumptions we develop a comprehensive approach in which both branches of the adsorption isotherm may be used for the PSD analysis of mesoporous oxide materials. We validate this approach by using Ar data measured at 87 K on the reference set of MCM-41 silica samples (Kruk and Jaroniec, 2000). The generated kernels are smooth, do not show layering transitions and fit accurately the reference data.

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 system size on the properties of a fluid adsorbed in a nanopore: Physical manifestations and methodological consequences

Hysteresis and discontinuities in the isotherms of a fluid adsorbed in a nanopore in general hamper the determination of equilibrium thermodynamic properties, even in computer simulations. A way around this has been to consider both a reservoir of small size and a pore of small extent in order to restrict the fluctuations of density and approach a classical van der Waals loop. We assess this suggestion by thoroughly studying through Monte Carlo simulations and density functional theory the influence of system size on the equilibrium configurations of the adsorbed fluid and on the resulting isotherms. We stress the importance of pore-symmetry-breaking states that even for modest pore sizes lead to discontinuous isotherms and we discuss the physical relevance of these states and the methodological consequences for computing thermodynamic quantities.

Influence of surface chemical heterogeneities on adsorption/desorption hysteresis and coexistence diagram of metastable states within cylindrical pores

Grand canonical Monte Carlo simulations are performed to determine the adsorption/desorption isotherms at different temperatures of a Lennard-Jones fluid confined within a simple model of cylindrical pores presenting chemical heterogeneities. A complex hysteresis loop is observed, showing hysteresis subloops (scanning curves). This is shown to be consistent with the existence of several metastable states (local minima in the system free energy). A recent extension to the Gibbs ensemble technique is then used to calculate the complete coexistence diagram of these local minima.

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.

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.

Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons

A new approach based on the nonlocal density functional theory to determine pore size distribution (PSD) of activated carbons and energetic heterogeneity of the pore wall is proposed. The energetic heterogeneity is modeled with an energy distribution function (EDF), describing the distribution of solid-fluid potential well depth (this distribution is a Dirac delta function for an energetic homogeneous surface). The approach allows simultaneous determining of the PSD (assuming slit shape) and EDF from nitrogen or argon isotherms at their respective boiling points by using a set of local isotherms calculated for a range of pore widths and solid-fluid potential well depths. It is found that the structure of the pore wall surface significantly differs from that ofgraphitized carbon black. This could be attributed to defects in the crystalline structure of the surface, active oxide centers, finite size of the pore walls (in either wall thickness or pore length), and so forth. Those factors depend on the precursor and the process of carbonization and activation and hence provide a fingerprint for each adsorbent. The approach allows very accurate correlation of the experimental adsorption isotherm and leads to PSDs that are simpler and more realistic than those obtained with the original nonlocal density functional theory.