Capillary condensation in a disordered mesoporous medium: a grand canonical Monte Carlo study

Machine Learning Predictions of Simulated Self-Diffusion Coefficients for Bulk and Confined Pure Liquids

Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. In this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.

Fluctuation-based fracture mechanics of heterogeneous materials

We present results of a hybrid analytical-simulation investigation of the fracture resistance of heterogeneous materials. We show that bond-energy fluctuations sampled by Monte Carlo simulations in the semigrand canonical ensemble provide a means to rationalize the complexity of heterogeneous fracture processes, encompassing probability and percolation theories of fracture. For a number of random and textured model materials, we derive upper and lower bounds of fracture resistance and link bond fracture fluctuations to statistical descriptors of heterogeneity, such as two-point correlation functions, to identify the origin of toughening mechanisms. This includes a shift from short- to long-range interactions of bond fracture processes in random systems to the transition from critical to subcritical bond fracture percolation in textured materials and the activation of toughness reserves at compliant interfaces. Induced by elastic mismatch, they connect to a number of disparate experimental observations, including toughening of brittle solids by deformable polymers or organics in, e.g., gas shale, nacre; stress-induced transformational toughening in ceramics; and toughening of sparse elastic networks in hydrogels, to name a few.

Advances in the Application of Vegetable-Oil-Based Cutting Fluids to Sustainable Machining Operations—A Review

Scientists and tribologists are currently exploring sustainable and inexhaustible lubricants as a result of increased awareness of environmental and health-related issues. Vegetable oils are being investigated as a potential form of environmentally friendly cutting fluids due to their excellent renewability, biodegradability, and lubricating performance. This report provides an overview of different vegetable oils used as cutting fluids in the machining of engineering materials. The effects of virgin vegetable oils, emulsified vegetable-based oils, and vegetable-oil-based nano-cutting fluids on the cutting force, the surface finish of machined parts, the tool wear, and the temperature of the cutting area were surveyed critically. Compared to mineral-oil-based cutting fluids, studies have demonstrated that vegetable-oil-based cutting fluids meet cleaner manufacturing standards with good or better efficiency.

Nanoscale Accessible Porosity as a Key Parameter Depicting the Topological Evolution of Organic Porous Networks

A significant part of the hydrocarbons contained in source rocks remains confined within the organic matter-called kerogen-from where they are generated. Understanding the sorption and transport properties of confined hydrocarbons within the kerogens is, therefore, paramount to predict production. Specifically, knowing the impact of thermal maturation on the evolution of the organic porous network is key. Here, we propose an experimental procedure to study the interplay between the chemical evolution and the structural properties of the organic porous network at the nanometer scale. First, the organic porous networks of source rock samples, covering a significant range of natural thermal maturation experienced by the Vaca Muerta formation (Neuquén Basin, Argentina), are physically reconstructed using bright-field electron tomography. Their structural description allows us to measure crucial parameters such as the porosity, specific pore volume and surface area, aperture and cavity size distributions, and constriction. In addition, a model-free computation of the topological properties (effective porosity, connectivity, and tortuosity) is conducted. Overall, we document a general increase of the specific pore volume with thermal maturation. This controls the topological features depicting increasing accessibility to alkane molecules, sensed by the evolution of the effective porosity. Collectively, our results highlight the input of bright-field electron tomography in the study of complex disordered amorphous porous media, especially to describe the interplay between the structural features and transport properties of confined fluids.

Phase diagram of brittle fracture in the semi-grand-canonical ensemble

We present a simulation method to assess the quasistatic fracture resistance of materials. Set within a semi-grand-canonical Monte Carlo (SGCMC) simulation environment, an auxiliary field-the bond rupture potential-is introduced to generate a sufficiently large number of possible microstates in the semi-grand-canonical ensemble, and associated energy and bond fluctuations. The SGCMC approach permits identifying the full phase diagram of brittle fracture for harmonic and nonharmonic bond potentials, analogous to the gas-liquid phase diagram, with the equivalent of a liquidus line ending in a critical point. The phase diagram delineates a solid phase, a fractured phase, and a gas phase, and provides clear evidence of a first-order phase transition intrinsic to fracture. Moreover, energy and bond fluctuations generated with the SGCMC approach permit determination of the maximum energy dissipation associated with bond rupture, and hence of the fracture resistance of a widespread range of materials that can be described by bond potentials.

Effect of Confinement on Capillary Phase Transition in Granular Aggregates

Using a 3D mean-field lattice-gas model, we analyze the effect of confinement on the nature of capillary phase transition in granular aggregates with varying disorder and their inverse porous structures obtained by interchanging particles and pores. Surprisingly, the confinement effects are found to be much less pronounced in granular aggregates as opposed to porous structures. We show that this discrepancy can be understood in terms of the surface-surface correlation length with a connected path through the fluid domain, suggesting that this length captures the true degree of confinement. We also find that the liquid-gas phase transition in these porous materials is of second order nature near capillary critical temperature, which is shown to represent a true critical temperature, i.e., independent of the degree of disorder and the nature of the solid matrix, discrete or continuous. The critical exponents estimated here from finite-size scaling analysis suggest that this transition belongs to the 3D random field Ising model universality class as hypothesized by F. Brochard and P.G. de Gennes, with the underlying random fields induced by local disorder in fluid-solid interactions.

Next-Generation Accurate, Transferable, and Polarizable Potentials for Material Simulations

PHAHST (potentials with high accuracy, high speed, and transferability) intermolecular potential energy functions have been developed from first principles for H₂, N₂, the noble gases, and a metal-organic material, HKUST-1. The potentials are designed from the outset to be transferable to heterogeneous environments including porous materials, interfaces, and material simulations. This is accomplished by theoretically justified choices for all functional forms, parameters, and mixing rules, including explicit polarization in every environment and fitting to high quality electronic structure calculations using methods that are tractable for real systems. The models have been validated in neat systems by comparison to second virial coefficients and bulk pressure-density isotherms. For inhomogeneous applications, our main target, comparisons are presented to previously published experimental studies on the metal-organic material HKUST-1 including adsorption, isosteric heats of adsorption, binding site locations, and binding site energies. A systematic prescription is provided for developing compatible potentials for additional small molecules and materials. The resulting models are recommended for use in complex heterogeneous simulations where existing potentials may be inadequate.

Capillary Stress and Structural Relaxation in Moist Granular Materials

A numerical and theoretical framework to address the poromechanical effect of capillary stress in complex mesoporous materials is proposed and exemplified for water sorption in cement. We first predict the capillary condensation/evaporation isotherm using lattice-gas simulations in a realistic nanogranular cement model. A phase-field model to calculate moisture-induced capillary stress is then introduced and applied to cement at different water contents. We show that capillary stress is an effective mechanism for eigenstress relaxation in granular heterogeneous porous media, which contributes to the durability of cement.

Molecular Simulation of Naphthalene, Phenanthrene, and Pyrene Adsorption on MCM-41

The adsorption of three typical polycyclic aromatic hydrocarbons (PAHs), naphthalene, phenanthrene, and pyrene with different ring numbers, on a common mesoporous material (MCM-41) was simulated based on a well-validated model. The adsorption equilibriums (isotherms), states (angle distributions and density profiles), and interactions (radial distribution functions) of three PAHs within the mesopores were studied in detail. The results show that the simulated isotherms agreed with previous experimental results. Each of the PAHs with flat molecules showed an adsorption configuration that was parallel to the surface of the pore, in the following order according to the degree of arrangement: pyrene (Pyr) > phenanthrene (Phe) > naphthalene (Nap). In terms of the interaction forces, there were no hydrogen bonds or other strong polar forces between the PAHs and MCM-41, and the O⁻H bond on the adsorbent surface had a unique angle in relation to the PAH molecular plane. The polarities of different H atoms on the PAHs were roughly the same, while those of the C atoms on the PAHs decreased from the molecular centers to the edges. The increasing area of the π-electron plane on the PAHs with the increasing ring number could lead to stronger adsorption interactions, and thus a shorter distance between the adsorbate and the adsorbent.

Mesoscale structure, mechanics, and transport properties of source rocks' organic pore networks

In source rocks, natural hydrocarbons are generated from organic matter dispersed in a fine-grained mineral matrix. The potential recovery of hydrocarbons is therefore influenced by the geometry of the organic hosted porous networks. Here, the three-dimensional structures of such networks are revealed using electron tomography with a subnanometer resolution. The reconstructions are first characterized in terms of morphology and topology and then used to build a multiscale simulation tool to study the mechanics and the transport properties of confined fluids. Our results offer evidence of the prevalent role of connected nanopores, which subsequently constitutes a material limit for long-term hydrocarbon production.

Organic matter is responsible for the generation of hydrocarbons during the thermal maturation of source rock formation. This geochemical process engenders a network of organic hosted pores that governs the flow of hydrocarbons from the organic matter to fractures created during the stimulation of production wells. Therefore, it can be reasonably assumed that predictions of potentially recoverable confined hydrocarbons depend on the geometry of this pore network. Here, we analyze mesoscale structures of three organic porous networks at different thermal maturities. We use electron tomography with subnanometric resolution to characterize their morphology and topology. Our 3D reconstructions confirm the formation of nanopores and reveal increasingly tortuous and connected pore networks in the process of thermal maturation. We then turn the binarized reconstructions into lattice models including information from atomistic simulations to derive mechanical and confined fluid transport properties. Specifically, we highlight the influence of adsorbed fluids on the elastic response. The resulting elastic energy concentrations are localized at the vicinity of macropores at low maturity whereas these concentrations present more homogeneous distributions at higher thermal maturities, due to pores’ topology. The lattice models finally allow us to capture the effect of sorption on diffusion mechanisms with a sole input of network geometry. Eventually, we corroborate the dominant impact of diffusion occurring within the connected nanopores, which constitute the limiting factor of confined hydrocarbon transport in source rocks.

Stochastic analysis of capillary condensation in disordered mesopores

Most mesoporous materials of practical interest are inherently disordered, which has a significant impact on the condensation and evaporation of vapours in their pores. Traditionally, the effect of disorder is theoretically analyzed in a perturbative approach whereby slight elements of disorder (constriction, corrugation) are added to geometrically ideal pores. We propose an alternative approach, which consists of using a stochastic geometrical model to describe both the porous material and the condensate within the pores. This is done through a multiphase generalisation of the standard Gaussian random field model of disordered materials. The model parameters characterising the condensate provide a low-dimensional approximation of its configuration space, and we use a Derjaguin-Broekhoff-de Boer approximation to calculate the free-energy landscape. Our analysis notably questions the existence of vapour-like metastable states in realistically disordered mesoporous materials. Beyond capillary condensation, our general methodology is applicable to a broad array of confined phenomena.

Adhesion and friction in polymer films on solid substrates: conformal sites analysis and corresponding surface measurements

In this work, we present a statistical mechanical analysis to elucidate the molecular-level factors responsible for the static and dynamic properties of polymer films. This analysis, which we term conformal sites theory, establishes that three dimensionless parameters play important roles in determining differences from bulk behavior for thin polymer films near to surfaces: a microscopic wetting parameter, α_wx, defined as the ratio of polymer-substrate interaction to polymer-polymer interaction; a dimensionless film thickness, H*; and dimensionless temperature, T*. The parameter α_wx introduced here provides a more fundamental measure of wetting than previous metrics, since it is defined in terms of intermolecular forces and the atomic structure of the substrate, and so is valid at the nanoscale for gas, liquid or solid films. To test this theoretical analysis, we also report atomic force microscopy measurements of the friction coefficient (μ), adhesion force (F_A) and glass transition temperature (T_g) for thin films of two polymers, poly(methyl methacrylate) (PMMA) and polystyrene (PS), on two planar substrates, graphite and silica. Both the friction coefficient and the glass transition temperature are found to increase as the film thickness decreases, and this increase is more pronounced for the graphite than for the silica surface. The adhesion force is also greater for the graphite surface. The larger effects encountered for the graphite surface are attributed to the fact that the microscopic wetting parameter, α_wx, is larger for graphite than for silica, indicating stronger attraction of polymer chains to the graphite surface.

Mesoscale texture of cement hydrates

Strength and other mechanical properties of cement and concrete rely upon the formation of calcium-silicate-hydrates (C-S-H) during cement hydration. Controlling structure and properties of the C-S-H phase is a challenge, due to the complexity of this hydration product and of the mechanisms that drive its precipitation from the ionic solution upon dissolution of cement grains in water. Departing from traditional models mostly focused on length scales above the micrometer, recent research addressed the molecular structure of C-S-H. However, small-angle neutron scattering, electron-microscopy imaging, and nanoindentation experiments suggest that its mesoscale organization, extending over hundreds of nanometers, may be more important. Here we unveil the C-S-H mesoscale texture, a crucial step to connect the fundamental scales to the macroscale of engineering properties. We use simulations that combine information of the nanoscale building units of C-S-H and their effective interactions, obtained from atomistic simulations and experiments, into a statistical physics framework for aggregating nanoparticles. We compute small-angle scattering intensities, pore size distributions, specific surface area, local densities, indentation modulus, and hardness of the material, providing quantitative understanding of different experimental investigations. Our results provide insight into how the heterogeneities developed during the early stages of hydration persist in the structure of C-S-H and impact the mechanical performance of the hardened cement paste. Unraveling such links in cement hydrates can be groundbreaking and controlling them can be the key to smarter mix designs of cementitious materials.

Tuning the adsorption and separation properties of noble gases and N2 in CuBTC by ligand functionalization

The Grand Canonical Monte Carlo method was used to investigate adsorption and separation properties of noble gas and N₂ mixtures on CuBTCs, functionalized with different groups, including amino, hydroxyl and fluorine groups.

Cold-induced spreading of water drops on hydrophobic surfaces

Superhydrophobic surfaces are characterized by their peculiarities, such as water-repellent, anti-icing, and freezing-delay properties. Wetting dynamics of deposited water drops on cooling hydrophobic surfaces, which directly affects the aforementioned properties, has not been studied thoroughly. Here, water drops are cooled on different hydrophobic surfaces in a controlled environment. During the cooling process, a significant increase in the drop footprint and decrease in the apparent contact angle are observed because of premature and capillary condensation, followed by thin water film formation adjacent to the solid-liquid-gas line. The water thin film propagates on the hydrophobic substrates radially away from the trijunction, followed by spreading of the drop on the film, which was experimentally validated through high-speed visualization. In addition, the roles of physical variables, such as the substrate temperature, humidity of surrounding air, types of hydrophobic surfaces, surface roughness, and drop volume, on post-spreading shape are investigated experimentally.

Molecular simulation of adsorption and transport in hierarchical porous materials

Adsorption and transport in hierarchical porous solids with micro- (~1 nm) and mesoporosities (>2 nm) are investigated by molecular simulation. Two models of hierarchical solids are considered: microporous materials in which mesopores are carved out (model A) and mesoporous materials in which microporous nanoparticles are inserted (model B). Adsorption isotherms for model A can be described as a linear combination of the adsorption isotherms for pure mesoporous and microporous solids. In contrast, adsorption in model B departs from adsorption in pure microporous and mesoporous solids; the inserted microporous particles act as defects, which help nucleate the liquid phase within the mesopore and shift capillary condensation toward lower pressures. As far as transport under a pressure gradient is concerned, the flux in hierarchical materials consisting of microporous solids in which mesopores are carved out obeys the Navier-Stokes equation so that Darcy's law is verified within the mesopore. Moreover, the flow in such materials is larger than in a single mesopore, due to the transfer between micropores and mesopores. This nonzero velocity at the mesopore surface implies that transport in such hierarchical materials involves slippage at the mesopore surface, although the adsorbate has a strong affinity for the surface. In contrast to model A, flux in model B is smaller than in a single mesopore, as the nanoparticles act as constrictions that hinder transport. By a subtle effect arising from fast transport in the mesopores, the presence of mesopores increases the number of molecules in the microporosity in hierarchical materials and, hence, decreases the flow in the micropores (due to mass conservation). As a result, we do not observe faster diffusion in the micropores of hierarchical materials upon flow but slower diffusion, which increases the contact time between the adsorbate and the surface of the microporosity.

Atomistic model of micelle-templated mesoporous silicas: structural, morphological, and adsorption properties

The structural, morphological, and adsorption properties of MCM-41 porous silicas are investigated using a realistic numerical model obtained by means of ab initio calculations [Ugliengo, P.; et al. Adv. Mater.2008, 20, 1]. Simulated X-ray diffraction, small angle neutron scattering, and electronic microscopy for the atomistic model are in good agreement with experimental data. The morphological features are also assessed from chord length distributions and porous volume and specific geometrical surface calculations, etc. The N(2), CO(2), and H(2)O adsorption isotherms in the atomistic model of MCM-41 are also in reasonable agreement with their experimental counterpart. An important finding of the present work is that water forms a film adsorbed on specific hydrophilic regions of the surface while the rest of the surface is depleted in water molecules. This result suggests that the surface of MCM-41 materials is heterogeneous, as it is made up of both hydrophilic and hydrophobic patches. While adsorption and irreversible capillary condensation can be described using the thermodynamical approach by Derjaguin (also known as the Derjaguin-Broekhoff-De Boer model), the Freundlich equation fits nicely the data for reversible and continuous filling in small pores.

Adsorption, capillary bridge formation, and cavitation in SBA-15 corrugated mesopores: a Derjaguin-Broekhoff-de Boer analysis

A Derjaguin-Broekhoff-de Boer analysis of adsorption and desorption in SBA-15 mesoporous silica is presented, using realistic geometrical models that account for the pore corrugation in these materials. The model parameters are derived from independent electron tomography and small-angle scattering characterization. A geometrical characteristic of the pore that is found to be important for adsorption is the corrugation length, l(C), which describes the longitudinal size of the geometrical defects along a given pore. Capillary bridges are possible only for large values of l(C). The results are explained in terms of two spinodal and two equilibrium pressures, characterizing the wide and the narrow sections of the pores. Simplified analytical expressions are obtained, which provide necessary conditions for bridge formation and for cavitation in terms of the radii of the narrow and wide sections of the pores, as well as of l(C). Quite generally, the results show that the deviation of the pore shape from that of ideal cylinders is key to understanding adsorption and desorption in corrugated mesopores, notably in SBA-15.

Enhanced CO2 solubility in hybrid MCM-41: molecular simulations and experiments

Grand canonical Monte Carlo simulations are performed in a hybrid adsorbent model in order to interpret the CO(2) solubility behavior. The hybrid adsorbent is prepared by confining a physical solvent (OMCTS) into the pores of a mimetic MCM-41 solid support. As a result, simulated adsorption isotherms of CO(2) nicely match the experimental data for three distinctive systems: bulk solvent, raw MCM-41, and hybrid MCM-41. The microscopic mechanisms underlying the apparition of enhanced solubility are then clearly identified. In fact, the presence of solvent molecules favors the layering of CO(2) molecules within the pores; therefore, the CO(2) solubility in the hybrid adsorbent markedly increases in comparison to that found in the raw adsorbent as well as in the bulk solvent. In addition, a good understanding of confined solvents' properties and solid surface structures is essential to fully evaluate the efficiency of hybrid adsorbents in capturing CO(2). The sorbent-solid interactions along with the solvent molecular size's impact on CO(2) solubility are therefore investigated in this study. We found that an ideal hybrid system should possess a weak solvent-solid interaction but a strong solvent-CO(2) interaction. Besides, an optimal solvent size is obtained for the enhanced CO(2) solubility in the hybrid system. According to the simulation results, the solvent layer builds pseudomicropores inside the mesoporous MCM-41, enabling more CO(2) molecules to be absorbed under the greater influence of spatial confinement and surface interaction. In addition, the molecular sieving effect is clearly observed in the case of larger solvent molecular sizes.

Comparison between Adsorption in Pores of a Simple Geometry and Realistic Models of Porous Materials

We study, by means of Grand Canonical Monte Carlo simulations, the effect of disorder in pore morphology and topology on adsorption of simple fluids at 77 K. We consider the case of microporous carbons and mesoporous silicas. For both systems, we found that morphological and topological disorder affects the adsorption isotherm. Consequently, we show that complex porous networks cannot be modeled as a collection of unconnected pores of simple geometry. In the case of microporous carbons, the isosteric heat of adsorption for disordered systems differs significantly from that obtained for an assembly of unconnected slit pores. In contrast, the isosteric heat of adsorption for the disordered mesoporous silica sample compares reasonably well with that observed for a silica material having regular cylindrical pores.

Molecular simulation of nitrogen adsorption in nanoporous silica

This article reports on a molecular simulation study of nitrogen adsorption and condensation at 77 K in atomistic silica cylindrical nanopores (MCM-41). Two models are considered for the nitrogen molecule and its interaction with the silica substrate. In the "pea" model, the nitrogen molecule is described as a single Lennard-Jones sphere and only Lennard-Jones interactions between the nitrogen molecule and the oxygens atoms of the silica substrate are taken into account. In the "bean" model (TraPPE force field), the nitrogen molecule is composed of two Lennard-Jones sites and a linear array of three charges on the atomic positions and at the center of the nitrogen-nitrogen bond. In the bean model, the interactions between the sites on the nitrogen molecule and the Si, O, and H atoms of the substrate are the sum of the Coulombic and dispersion interactions with a repulsive short-range contribution. The data obtained with the pea and bean models in silica nanopores conform to the typical behavior observed in the experiments for adsorption/condensation in cylindrical MCM-41 nanopores; the adsorbed amount increases continuously in the multilayer adsorption regime until an irreversible jump occurs because of capillary condensation and evaporation of the fluid within the pore. Our results suggest that the pea model can be used for characterization purposes where one is interested in capturing the global experimental behavior upon adsorption and desorption in silica nanopores. However, the bean model is more suitable to investigating the details of the interaction with the surface because this model, which accounts for the partial charges located on the nitrogen atoms of the molecule (quadrupole), allows a description of the specific interactions between this adsorbate and silica surfaces (silanol groups and siloxane bridges) or grafted silica surfaces. In particular, the bean model provides a more realistic picture of nitrogen adsorption in the vicinity of silica surfaces or confined in silica nanopores, where the isosteric heat of adsorption curves show that the nitrogen molecule in this model is sensitive to the surface heterogeneity.

Adsorption and structure of benzene on silica surfaces and in nanopores

Grand canonical Monte Carlo simulations are used to study the adsorption of benzene on atomistic silica surfaces and in cylindrical nanopores. The effect of temperature and surface chemistry is addressed by studying the adsorption of benzene at 293 and 323 K on both fully and partially hydroxylated silica surfaces or nanopores. We also consider the adsorption of benzene in a cylindrical nanopore of diameter D=3.6 nm and compare our results with those obtained for planar surfaces. The structure of benzene in the vicinity of the planar surface and confined in the cylindrical nanopore is determined by calculating orientational order parameters and examining positional pair correlation functions. The density profiles of adsorbed benzene reveal the strong layering of the adsorbate, which decays with the distance from the silica surface. At a given temperature and at low pressures, the film adsorbed at the fully hydroxylated silica surface is larger than that for the partially hydroxylated silica surface. This result is due to an increase in the density of silanol groups that induces an increase in the polarity of the silica surface, which becomes more attractive for the adsorbate. Our results also suggest that the benzene molecules prefer an orientation in which their ring is nearly perpendicular to the surface when fully hydroxylated surfaces are considered. When partially hydroxylated surfaces are considered, a second preferential orientation is observed where the benzene ring forms an angle of approximately 50 degrees with the silica surface. In this case, the average orientation of the benzene molecules appears disordered as in the bulk phase. These results suggest that determining the experimental orientation of benzene in the vicinity of a silica surface is a difficult task even when the surface chemistry is known. Capillary condensation in the nanopores involves a transition from a partially filled pore (a thin film adsorbed at the pore surface) to a completely filled pore configuration where the confined liquid coexists at equilibrium with the external gas phase. The disordered orientation of the adsorbed benzene molecules in the case of the partially hydroxylated surface favors the condensation of benzene molecules (the condensation pressure for this substrate is lower than that for the fully hydroxylated surface). Finally, these results are consistent with the structural analysis, showing that (1) benzene tends to relax its liquid structure a little in order to optimize its molecular arrangement near the pore wall and (2) the disordering of the liquid structure induced by the surface becomes stronger as the interaction with the pore wall increases.

Modeling micelle-templated mesoporous material SBA-15: atomistic model and gas adsorption studies

We report the development of a realistic molecular model for mesoporous silica SBA-15, which includes both the large cylindrical mesopores and the smaller micropores in the pore walls. The methodology for modeling the SBA-15 structure involves molecular and mesoscale simulations combined with geometrical interpolation techniques. First, a mesoscale model is prepared by mimicking the synthesis process using lattice Monte Carlo simulations. The main physical features of this mesoscale pore model are then carved out of an atomistic silica block; both the mesopores and the micropores are incorporated from the mimetic simulations. The calculated pore size distribution, surface area, and simulated TEM images of the model structure are in good agreement with those obtained from experimental samples of SBA-15. We then investigate the adsorption of argon in this structure using Grand Canonical Monte Carlo (GCMC) simulations. The adsorption results for our SBA-15 model are compared with those for a similar model that does not include the micropores; we also compare with results obtained in a regular cylindrical pore. The simulated adsorption isotherm for the SBA-15 model shows semiquantitative agreement with the experimental isotherm for a SBA-15 sample having a similar pore size. We observe that the presence of the micropores leads to increased adsorption at low pressure compared to the case of a model without micropores in the pore walls. At higher pressures, for all models, the filling proceeds via the monolayer-multilayer adsorption on the mesopore surface followed by capillary condensation, which is mainly controlled by the mesopore diameter and is not influenced by the presence of the micropores.

Simple phenomenological model for phase transitions in confined geometry. 2. Capillary condensation/evaporation in cylindrical mesopores

A simple phenomenological model that describes capillary condensation and evaporation of pure fluids confined in cylindrical mesopores is presented. Following the work of Celestini (Celestini, F. Phys. Lett. A 1997, 228, 84), the free energy density of the system is derived using interfacial tensions and a corrective term that accounts for the interaction coupling between the vapor/adsorbed liquid and the adsorbed liquid/adsorbent interfaces. This corrective term is shown to be consistent with the Gibbs adsorption isotherm and assessed by standard adsorption tests. This model reveals that capillary condensation and evaporation are metastable and equilibrium processes, respectively, hence exhibiting the existence of a hysteresis loop inadsorption/desorption isotherm that is well-known in experiment. We extend the phenomenological model of Celestini to give a quantitative description of adsorption on the pore wall and hysteresis width evolution with temperature and confinement. Direct quantitative comparison is made with experimental data for confined argon. Used as a characterizing tool, this integrated model allows in a single fit of an experimental adsorption/desorption isotherm assessing essential characterization data such as the specific surface area, pore volume, and mean pore size.

Adsorption of simple fluid on silica surface and nanopore: effect of surface chemistry and pore shape

This paper reports a molecular simulation study on the adsorption of simple fluids (argon at 77 K) on hydroxylated silica surfaces and nanopores. The effect of surface chemistry is addressed by considering substrates with either partially or fully hydroxylated surfaces. We also investigate the effect of pore shape on adsorption and capillary condensation by comparing the results for cylindrical and hexagonal nanopores having equivalent sections (i.e., equal section areas). Due to the increase in the polarity of the surface with the density of OH groups, the adsorbed amounts for fully hydroxylated surfaces are found to be larger than those for partially hydroxylated surfaces. Both the adsorption isotherms for the cylindrical and hexagonal pores conform to the typical behavior observed in the experiments for adsorption/condensation in cylindrical nanopores MCM-41. Capillary condensation occurs through an irreversible discontinuous transition between the partially filled and the completely filled configurations, while evaporation occurs through the displacement at equilibrium of a hemispherical meniscus along the pore axis. Our data are also used to discuss the effect of surface chemistry and pore shape on the BET method. The BET surface for fully hydroxylated surfaces is much larger (by 10-20%) than the true geometrical surface. In contrast, the BET surface significantly underestimates the true surface when partially hydroxylated surfaces are considered. These results suggest that the surface chemistry and the choice of the system adsorbate/adsorbent is crucial in determining the surface area of solids using the BET method.

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.