Phase behavior and dynamics of fluids in mesoporous glasses

Mercury porosimetry in mesoporous glasses: a comparison of experiments with results from a molecular model

We present results from experiments and molecular modeling of mercury porosimetry into mesoporous Vycor and controlled pore glass (CPG) solid materials. The experimental intrusion/extrusion curves show a transition from a type H2 hysteresis for the Vycor glass to a type H1 hysteresis for the CPG. Mercury entrapment is observed in both materials, but we find that the amount of entrapped mercury depends on the chosen experimental relaxation time. No additional entrapment is found in a second intrusion/extrusion cycle, but hysteresis is still observed. This indicates that hysteresis and entrapment are of different origin. The experimental observations are qualitatively reproduced in theoretical calculations based on lattice models, which provide significant insights of the molecular mechanisms occurring during mercury porosimetry experiments in these porous glasses.

A computational study of the reconstruction of amorphous mesoporous materials from gas adsorption isotherms and structure factors via evolutionary optimization

A general method for the three-dimensional reconstruction of mesoporous materials by evolutionary optimization against target data is developed. The method is applied specifically in reconstruction of amorphous material models using gas adsorption data, structure factor data, or a combination of both. A recently introduced lattice-gas approach is used to model adsorption in these calculations, and a high-pass limited Fourier representation is used to facilitate evolution of large-scale structures during the optimization. Reconstructions are made of several material models which mimic real materials obtained either by phase separation and etching or by sol-gel processing. Analysis of the reconstructions provides considerable insight into the type and quantity of structural information probed by gas adsorption and small-angle scattering experiments. We find that reconstructions based only on structure factors tend to underestimate the mean pore size. We also find that in many cases excellent reconstructions can be obtained using only adsorption-branch data, and that in all cases reconstructions based jointly on both types of data are superior to those based only on one, suggesting that these measures contain "complementary" information. It is also found that in most cases the use of desorption data is not warranted, and that the use of adsorption data taken at many temperatures will not improve reconstructions. The reproducibility of the method is shown to be satisfactory. The method can be computationally expensive if gas adsorption data are used, but it is easily parallelized, and therefore results can still be obtained in reasonable time. Finally, the possible application of this approach to real systems, including templated porous materials, is discussed.

Gas adsorption in mesoporous micelle-templated silicas: MCM-41, MCM-48, and SBA-15

This paper reports a molecular simulation and experimental study on the adsorption and condensation of simple fluids in mesoporous micelle-templated silicas MCM-41, MCM-48, and SBA-15. MCM-41 is described as a regular cylindrical silica nanopore, while SBA-15 is assumed to be made up of cylindrical nanopores that are connected through lateral channels. The 3D-connected topology of MCM-48 is described using a gyroid periodic minimal surface. Argon adsorption at 77 K is calculated for the three materials using Grand Canonical Monte Carlo simulations. Qualitative comparison with experiments for nitrogen adsorption in mesoporous micelle-templated silicas is made. The adsorption isotherm for SBA-15 resembles that for MCM-41. In particular, capillary condensation and evaporation are not affected by the presence of the connecting lateral channels. In contrast, the argon adsorption isotherm for MCM-48 departs from that for MCM-41 having the same pore size. While condensation in MCM-41 is a one-step process, filling of MCM-48 involves two successive jumps in the adsorbed amounts which correspond to condensation in different domains of the porosity. The condensation pressure for MCM-48 is larger than that for MCM-41. We attribute this result to the morphology of the MCM-48 surface (made up of both concave and convex regions) that differs from that for MCM-41 (concave only). Our results suggest that the pore connectivity affects pore filling when the size of the connections is comparable to that of the nanopores.

Exploration of molecular dynamics during transient sorption of fluids in mesoporous materials

In recent years, considerable progress has been made in the development of novel porous materials with controlled architectures and pore sizes in the mesoporous range. An important feature of these materials is the phenomenon of adsorption hysteresis: for certain ranges of applied pressure, the amount of a molecular species adsorbed by the mesoporous host is higher on desorption than on adsorption, indicating a failure of the system to equilibrate. Although this phenomenon has been known for over a century, the underlying internal dynamics responsible for the hysteresis remain poorly understood. Here we present a combined experimental and theoretical study in which microscopic and macroscopic aspects of the relaxation dynamics associated with hysteresis are quantified by direct measurement and computer simulations of molecular models. Using nuclear magnetic resonance techniques and Vycor porous glass as a model mesoporous system, we have explored the relationship between molecular self-diffusion and global uptake dynamics. For states outside the hysteresis region, the relaxation process is found to be essentially diffusive in character; within the hysteresis region, the dynamics slow down dramatically and, at long times, are dominated by activated rearrangement of the adsorbate density within the host material.

Freezing of fluids confined in a disordered nanoporous structure

Freezing of a simple fluid in a disordered nanoporous carbon is studied using molecular simulations. Only partial crystallization occurs, and the confined phase is composed of crystalline and amorphous nanodomains. This freezing behavior departs strongly from that for nanopores of simple geometry. We present a method for analyzing the freezing in such disordered materials in terms of a transition in the average size and number of crystalline clusters. The results provide a basis for the interpretation of experiments on freezing in such materials, particularly 1H-NMR and scattering experiments.

Gas adsorption and desorption in silica aerogels: a theoretical study of scattering properties

We present a numerical study of the structural correlations associated with gas adsorption and desorption in silica aerogels in order to provide a theoretical interpretation of scattering experiments. Following our earlier work, we use a coarse-grained lattice-gas description and determine the nonequilibrium behavior of the adsorbed gas within a local mean-field analysis. We focus on the differences between the adsorption and desorption mechanisms and their signature in the fluid-fluid and gel-fluid structure factors as a function of temperature. At low temperature, but still in the regime where the isotherms are continuous, we find that the adsorbed fluid density, during both filling and draining, is correlated over distances that may be much larger than the gel correlation length. In particular, extended fractal correlations may occur during desorption, indicating the existence of a ramified cluster of vapor filled cavities. This also induces an important increase of the scattering intensity at small wave vectors. The similarity and differences with the scattering of fluids in other porous solids such as Vycor are discussed.

Wetting of rings on a nanopatterned surface: a lattice model study

We perform mean-field density functional theory calculations on a lattice model to study the wetting of a solid substrate decorated with a ring pattern of nanoscale dimensions. We have found three different liquid morphologies on the substrate: a ring morphology where the liquid covers the pattern, a bulge morphology where a droplet is forming on one side of the ring, and a morphology where the liquid forms a cap spanning the nonwetting disk inside the pattern. We investigate the relative stability of these morphologies as a function of the ring size, wall-fluid interaction, and temperature. The results found are in very good agreement with experiments and calculations performed on similar systems at a micrometer length scale. The bulge morphology has also been observed in Monte Carlo simulations of the lattice model. Our results show that (i) morphologies of wetting patterns previously observed on a much larger (microm) scale can also form on a nm length scale, (ii) whether or not this happens depends crucially on the size of the wettable pattern, and (iii) the wettable ring may only be partially wet by the bulge morphology of the fluid. This morphology is a result of a spontaneously broken symmetry in the system.

Mean-field theory of liquid droplets on roughened solid surfaces: application to superhydrophobicity

We present calculations of the density distributions and contact angles of liquid droplets on roughened solid surfaces for a lattice gas model solved in a mean-field approximation. For the case of a smooth surface, this approach yields contact angles that are well described by Young's equation. We consider rough surfaces created by placing an ordered array of pillars on a surface, modeling so-called superhydrophobic surfaces, and we have made calculations for a range of pillar heights. The apparent contact angle follows two regimes as the pillar height increases. In the first regime, the liquid penetrates the interpillar volume, and the contact angle increases with pillar height before reaching a constant value. This behavior is similar to that described by the Wenzel equation for contact angles on rough surfaces, although the contact angles are underestimated. In the second regime, the liquid does not penetrate the interpillar volume substantially, and the contact angle is independent of the pillar height. This situation is similar to that envisaged in the Cassie-Baxter equation for contact angles on heterogeneous surfaces, but the contact angles are overestimated by this equation. For larger pillar heights, two states of the droplet can be observed, one Wenzel-like and the other Cassie-like.

Soft core fluid in a quenched matrix of soft core particles: a mobile mixture in a model gel

We present a density-functional study of a binary phase-separating mixture of soft core particles immersed in a random matrix of quenched soft core particles of larger size. This is a model for a binary polymer mixture immersed in a cross-linked rigid polymer network. Using the replica "trick" for quenched-annealed mixtures we derive an explicit density functional theory that treats the quenched species on the level of its one-body density distribution. The relation to a set of effective external potentials acting on the annealed components is discussed. We relate matrix-induced condensation in bulk to the behavior of the mixture around a single large particle. The interfacial properties of the binary mixture at a surface of the quenched matrix display a rich interplay between capillary condensation inside the bulk matrix and wetting phenomena at the matrix surface.

Water Condensation in Hydrophobic Silicalite-1 Zeolite: A Molecular Simulation Study

We report grand canonical Monte Carlo simulations of the gas and liquid phase adsorption of water in silicalite-1 zeolite. Simple but effective models and simulation methods, found useful for studying gas adsorption in nanoporous materials, have been extended to describe the intrusion/extrusion cycle of water in this hydrophobic solid. The picture of water confined to hydrophobic spaces of nanoscopic dimensions that emerges from this study is one of a strongly depleted and highly inhomogeneous fluid.

Does water condense in carbon pores?

Using grand canonical Monte Carlo (GCMC) simulations of molecular models, we investigate the nature of water adsorption and desorption in slit pores with graphitelike surfaces. Special emphasis is placed on the question of whether water exhibits capillary condensation (i.e., condensation when the external pressure is below the bulk vapor pressure). Three models of water have been considered. These are the SPC and SPC/E models and a model where the hydrogen bonding is described by tetrahedrally coordinated square-well association sites. The water-carbon interaction was described by the Steele 10-4-3 potential. In addition to determining adsorption/desorption isotherms, we also locate the states where vapor-liquid equilibrium occurs for both the bulk and confined states of the models. We find that for wider pores (widths >1 nm), condensation does not occur in the GCMC simulations until the pressure is higher than the bulk vapor pressure, P0. This is consistent with a physical picture where a lack of hydrogen bonding with the graphite surface destabilizes dense water phases relative to the bulk. For narrow pores where the slit width is comparable to the molecular diameter, strong dispersion interactions with both carbon surfaces can stabilize dense water phases relative to the bulk so that pore condensation can occur for P < P0 in some cases. For the narrowest pores studied--a pore width of 0.6 nm--pore condensation is again shifted to P > P0. The phase-equilibrium calculations indicate vapor-liquid coexistence in the slit pores for P < P0 for all but the narrowest pores. We discuss the implications of our results for interpreting water adsorption/desorption isotherms in porous carbons.

Lattice density functional for colloid-polymer mixtures: comparison of two fundamental measure theories

We consider a binary mixture of colloid and polymer particles with positions on a simple cubic lattice. Colloids exclude both colloids and polymers from nearest neighbor sites. Polymers are treated as effective particles that are mutually noninteracting, but exclude colloids from neighboring sites; this is a discrete version of the (continuum) Asakura-Oosawa-Vrij model. Two alternative density functionals are proposed and compared in detail. The first is based on multioccupancy in the zero-dimensional limit of the bare model, analogous to the corresponding continuum theory that reproduces the bulk fluid free energy of free volume theory. The second is based on mapping the polymers onto a multicomponent mixture of polymer clusters that are shown to behave as hard cores; the corresponding property of the extended model in strong confinement permits direct treatment with lattice fundamental measure theory. Both theories predict the same topology for the phase diagram with a continuous fluid-fcc freezing transition at low polymer fugacity and, upon crossing a tricritical point, a first-order freezing transition for high polymer fugacities with rapidly broadening density jump.

Application of the Bethe-Peierls approximation to a lattice-gas model of adsorption on mesoporous materials

We calculate adsorption and desorption isotherms in models of several classes of porous materials using a lattice-gas model solved in the Bethe-Peierls (quasichemical) approximation. Isotherms and fluid density profiles from the Bethe-Peierls and Bragg-Williams approximations are compared with grand-canonical Monte Carlo simulation results. The Bethe-Peierls approximation produces both more accurate adsorption and desorption isotherms and more realistic fluid density profiles than the Bragg-Williams approximation. Details of the application of the Bethe-Peierls approximation applied to a three-dimensionally inhomogeneous system are given. We show that the numerical solution of this theory can be accomplished using a self-consistent iterator very similar to that currently used in studies employing the Bragg-Williams approximation. This iterative scheme is substantially more efficient than the numerical optimization method used in many previous studies of lattice-gas models in the quasichemical approximation. We find that use of the Bethe-Peierls approximation is only slightly more computationally demanding than the Bragg-Williams approximation, and thus recommend it for use in future work on this class of models.

Dynamic aspects of mercury porosimetry: a lattice model study

Grand canonical Monte Carlo simulations using both Glauber dynamics and Kawasaki dynamics have been carried out for a recently developed lattice model of a nonwetting fluid confined in a porous material. The calculations are aimed at investigating the molecular scale mechanisms leading to mercury retention encountered during mercury porosimetry experiments. We first describe a set of simulations on slit and ink-bottle pores. We have studied the influence of the pore width parameter on the intrusion/extrusion curve shapes and investigated the corresponding mechanisms. Entrapment appears during Kawasaki dynamics simulations of extrusion performed on ink-bottle pores when the system is studied for short relaxation times. We then consider the more realistic and complex case of a Vycor glass building on recent work on the dynamics of adsorption of wetting fluids (Woo, H. J.; Monson, P. A. Phys. Rev. E 2003, 67, 041207). Our results suggest that mercury entrapment is caused by a decrease in the rate of mass transfer associated with the fragmentation of the liquid during extrusion.

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.

Mechanisms for gas adsorption and desorption in silica aerogels: the effect of temperature

We present a theoretical study of the adsorption and desorption mechanisms of fluids in silica aerogels, focusing on the effect of temperature. We adopt a coarse-grained lattice description in which the gel structure is generated by a diffusion-limited cluster-cluster aggregation algorithm and the fluid configurations are computed using local mean-field (i.e., density functional) theory. Our calculations reproduce qualitatively the changes in the shape of the hysteresis loops observed with (4)He in gels of varying porosity. We study in detail the morphology of the condensation and evaporation events that correspond to the irreversible processes (avalanches) which are at the origin of the hysteresis. Depending on porosity and temperature, these avalanches may be localized, involve regions that extend beyond the gel correlation length, or even span the entire sample. This makes difficult the characterization of aerogels based on analyzing sorption isotherms.

Modeling mercury porosimetry using statistical mechanics

We consider mercury porosimetry from the perspective of the statistical thermodynamics of penetration of a nonwetting liquid into a porous material under an external pressure. We apply density functional theory to a lattice gas model of the system and use this to compute intrusion/extrusion curves. We focus on the specific example of a Vycor glass and show that essential features of mercury porosimetry experiments can be modeled in this way. The lattice model exhibits a symmetry that provides a direct relationship between intrusion/extrusion curves for a nonwetting fluid and adsorption/desorption isotherms for a wetting fluid. This relationship clarifies the status of methods that are used for transforming mercury intrusion/extrusion curves into gas adsorption/desorption isotherms. We also use Monte Carlo simulations to investigate the nature of the intrusion and extrusion processes.

Modeling desorption of fluids from disordered mesoporous materials

The desorption mechanism of fluids in disordered mesoporous glasses is studied by Monte Carlo simulations of a coarse-grained lattice model with realistic matrix configurations representative of Vycor. Two methods of simulation are considered: grand canonical ensemble Monte Carlo simulations and dynamic Monte Carlo simulations which mimic the diffusion of the fluid in and out of the material using Kawasaki dynamics. In the grand canonical simulations, cavitation via nucleation of bubbles inside the pores plays the dominant role in determining the fluid configurations along the desorption isotherm. The Kawasaki dynamics simulations indicate that such configurations are achieved dynamically via the gradual advancement of macroscopic front interfaces toward the interior. This is made possible by the bubble nucleation mechanism operating on a length scale that is determined by both the typical pore size and the strength of the solid-fluid interaction.

Adsorption/desorption hysteresis in inkbottle pores: a density functional theory and Monte Carlo simulation study

The mechanisms of adsorption and desorption in inkbottle-shaped pores are considered for lattice models using grand canonical mean field density functional theory and Monte Carlo simulation. We find that they depend significantly on the particular pore geometry, the nature of the fluid-solid interaction, and the temperature. We find two mechanisms for desorption. One mechanism involves the emptying of the main cavity even as the density of fluid in the necks remains high, a mechanism observed recently in studies of an off-lattice model of an inkbottle. The other is a simultaneous desorption from the entire pore space, behavior that is more closely related to the traditional picture of pore blocking in the inkbottle system.

Local mean-field study of capillary condensation in silica aerogels

We apply local mean-field (i.e., density functional) theory to a lattice model of a fluid in contact with a dilute, disordered gel network. The gel structure is described by a diffusion-limited cluster aggregation model. We focus on the influence of porosity on both the hysteretic and the equilibrium behavior of the fluid as one varies the chemical potential at low temperature. We show that the shape of the hysteresis loop changes from smooth to rectangular as the porosity increases and that this change is associated with disorder-induced out-of-equilibrium phase transitions that differ in adsorption and in desorption. Our results provide insight in the behavior of 4He in silica aerogels.