Ultrasonic study of water adsorbed in nanoporous glasses.

Thermodynamic properties of fluids confined in nanopores differ from those observed in the bulk. To investigate the effect of nanoconfinement on water compressibility, we perform water sorption experiments on two nanoporous glass samples while concomitantly measuring the speed of longitudinal and shear ultrasonic waves in these samples. These measurements yield the longitudinal and shear moduli of the water-laden nanoporous glass as a function of relative humidity that we utilize in the Gassmann theory to infer the bulk modulus of the confined water. This analysis shows that the bulk modulus (inverse of compressibility) of confined water is noticeably higher than that of the bulk water at the same temperature. Moreover, the modulus exhibits a linear dependence on the Laplace pressure. The results for water, which is a polar fluid, agree with previous experimental and numerical data reported for nonpolar fluids. This similarity suggests that irrespective of intermolecular forces, confined fluids are stiffer than bulk fluids. Accounting for fluid stiffening in nanopores may be important for accurate interpretation of wave propagation measurements in fluid-filled nanoporous media, including in petrophysics, catalysis, and other applications, such as in porous materials characterization.

Use of Boundary-Driven Nonequilibrium Molecular Dynamics for Determining Transport Diffusivities of Multicomponent Mixtures in Nanoporous Materials.

The boundary-driven molecular modeling strategy to evaluate mass transport coefficients of fluids in nanoconfined media is revisited and expanded to multicomponent mixtures. The method requires setting up a simulation with bulk fluid reservoirs upstream and downstream of a porous media. A fluid flow is induced by applying an external force at the periodic boundary between the upstream and downstream reservoirs. The relationship between the resulting flow and the density gradient of the adsorbed fluid at the entrance/exit of the porous media provides for a direct path for the calculation of the transport diffusivities. It is shown how the transport diffusivities found this way relate to the collective, Onsager, and self-diffusion coefficients, typically used in other contexts to describe fluid transport in porous media. Examples are provided by calculating the diffusion coefficients of a Lennard-Jones (LJ) fluid and mixtures of differently sized LJ particles in slit pores, a realistic model of methane in carbon-based slit pores, and binary mixtures of methane with hypothetical counterparts having different attractions to the solid. The method is seen to be robust and particularly suited for the study of study of transport of dense fluids and liquids in nanoconfined media.

A Strong Test of Atomically Detailed Models of Molecular Adsorption in Zeolites Using Multilaboratory Experimental Data for CO2 Adsorption in Ammonium ZSM-5.

A recent international interlaboratory study led by the U.S. National Institute of Standards (NIST) reported CO2 adsorption isotherms measured independently by 11 groups on reference material RM 8852, an ammonium ZSM-5 zeolite. Good reproducibility and high reliability of this experimental data provide a strong test for the ability of atomically detailed models to predict adsorption of CO2 in zeolites. We developed force fields for CO2 in ammonium zeolites based on first-principles calculations and also independently performed experiments with RM 8852 by microcalorimetry. At low pressures good agreement was obtained between predictions and experiments. At high pressures, however, deviations were observed. We show that the charge-balancing cations in the experimental material are the predominant source of the discrepancy between simulation and experiment at high pressures; the experimental sample treatment causes deammoniation. In addition, accounting for a small amount of noncrystalline mesoporosity in the zeolite brings predictions into much better agreement with experiments.

Small Pore Aluminosilicate EMM-37: Synthesis and Structure Determination Using Continuous Rotation Electron Diffraction.

A new aluminosilicate zeolite, denoted EMM-37, with a 3D small pore channel system, has been synthesized using a diquaternary ammonium molecule as the structure directing agent (SDA) and metakaolin as the aluminum source. The structures of both as-made and calcined forms of EMM-37 were solved and refined using continuous rotation electron diffraction (cRED) data. cRED is a powerful method for the collection of 3D electron diffraction data from submicron- and nanosized crystals, which allows for successful solution and refinement of complex structures in symmetry as low as P1̅.

Structures and dimensions of micelle-templated nanoporous silicas derived from swollen spherical micelles of temperature-dependent size.

Subambient temperatures are employed in Pluronic-block-copolymer-templated syntheses of many large-pore silicas: SBA-15 (2-D hexagonal with cylindrical mesopores), FDU-12 (face-centered cubic with spherical mesopores), nanotubes and hollow nanospheres. Herein, the origin of a significant temperature dependence of the unit-cell parameter and pore diameter of silicas templated by swollen micelles of Pluronics under subambient conditions was elucidated. The temperature dependence of size of swollen spherical micelles of Pluronic F127 (EO106PO70EO106) in 2 M HCl solution was studied in 12-25 °C range using dynamic light scattering and was correlated with structure types, unit-cell sizes and pore sizes of silicas synthesized at four silica-precursor/Pluronic ratios with a swelling agent (toluene, ethylbenzene). The increase in size of swollen micelles with temperature decrease was paralleled by the increase in the unit-cell size and pore diameter, even if the micelle shape changed in the process of formation of the micelle-templated silica. The decrease in the silica-precursor/Pluronic F127 ratio at constant temperature triggered a succession of phases, including SBA-15 - nanotube sequence that may involve an intermediate nanotube bundle structure, which is uncommon and potentially useful. The temperature decrease also led to a succession of phases, including FDU-12 - SBA-15, hollow nanospheres - nanotube bundles, and nanotubes - SBA-15 sequences.

Control of zeolite framework flexibility and pore topology for separation of ethane and ethylene.

The discovery of new materials for separating ethylene from ethane by adsorption, instead of using cryogenic distillation, is a key milestone for molecular separations because of the multiple and widely extended uses of these molecules in industry. This technique has the potential to provide tremendous energy savings when compared with the currently used cryogenic distillation process for ethylene produced through steam cracking. Here we describe the synthesis and structural determination of a flexible pure silica zeolite (ITQ-55). This material can kinetically separate ethylene from ethane with an unprecedented selectivity of ~100, owing to its distinctive pore topology with large heart-shaped cages and framework flexibility. Control of such properties extends the boundaries for applicability of zeolites to challenging separations.

Determination of Isosteric Heat of Adsorption by Quenched Solid Density Functional Theory.

The heat of adsorption is one of the most important parameters characterizing energetic heterogeneity of the adsorbent surface. Heats of adsorption are either determined directly by calorimetry or calculated from adsorption isotherms measured at different temperatures using the thermodynamic Clausius-Clapeyron equation. Here, we present a method for calculating the isosteric heat of adsorption that requires as input only a single adsorption isotherm measured at one temperature. The proposed method is implemented with either nonlocal (NLDFT) or quenched solid (QSDFT) density functional theory models of adsorption that are currently widely used for calculating pore size distributions in various micro- and mesoporous solids. The pore size distribution determined from the same experimental isotherm is used for predicting the isosteric heat. The QSDFT method has advantages of taking into account two factors contributing to the structural heterogeneity of adsorbents: the molecular level roughness of the surface and the pore size distribution. The method is illustrated with examples of low temperature nitrogen and argon adsorption on selected samples of carbons of different degree of graphitization and MCM-41 mesoporous silicas of different pore size. The isosteric heat predictions from the NLDFT and QSDFT methods are compared against relevant experiments and the results of Monte Carlo (MC) simulations, with good agreement found in the cases where the surface model adequately reflects the pore surface roughness. Analyses with the QSDFT method show that the isosteric heat of adsorption significantly depends of the molecular level roughness of the adsorbent surface, which is ignored in NLDFT and MC models. The proposed QSDFT method with further verification can be used for calculating the isosteric heat as an additional parameter characterizing the adsorbent surface in parallel with routine calculations of the pore size distribution from a single adsorption isotherm.

New High- and Low-Temperature Phase Changes of ZIF-7: Elucidation and Prediction of the Thermodynamics of Transitions.

We have found that the 3D zeolitic imidazolate framework ZIF-7 exhibits far more complex behavior in response to the adsorption of guest molecules and changes in temperature than previously thought. We believe that this arises from the existence of different polymorphs and different types of adsorption sites. We report that ZIF-7 undergoes a displacive, nondestructive phase change upon heating to above ∼700 °C in vacuum, or to ∼500 °C in CO2 or N2. This is the first example of a temperature-driven phase change in 3D ZIF frameworks. We predicted the occurrence of the high-temperature transition on the basis of thermodynamic arguments and analyses of the solid free-energy differences obtained from CO2 and n-butane adsorption isotherms. In addition, we found that ZIF-7 exhibits complex behavior in response to the adsorption of CO2 manifesting in double transitions on adsorption isotherms and a doubling of the adsorption capacity. We report adsorption microcalorimetry, molecular simulations, and detailed XRD investigations of the changes in the crystal structure of ZIF-7. Our results highlight mechanistic details of the phase transitions in ZIF-7 that are driven by adsorption of guest molecules at low temperature and by entropic effects at high temperature. We derived a phase diagram of CO2 in ZIF-7, which exhibits surprisingly complex re-entrant behavior and agrees with our CO2 adsorption measurements over a wide range of temperatures and pressures. We predicted phase diagrams of CH4, C3H6, and C4H10. Finally, we modeled the temperature-induced transition in ZIF-7 using molecular dynamics simulations in the isobaric-isothermal ensemble, confirming our thermodynamic arguments.

First principles derived, transferable force fields for CO2 adsorption in Na-exchanged cationic zeolites.

The development of accurate force fields is vital for predicting adsorption in porous materials. Previously, we introduced a first principles-based transferable force field for CO2 adsorption in siliceous zeolites (Fang et al., J. Phys. Chem. C, 2012, 116, 10692). In this study, we extend our approach to CO2 adsorption in cationic zeolites which possess more complex structures. Na-exchanged zeolites are chosen for demonstrating the approach. These methods account for several structural complexities including Al distribution, cation positions and cation mobility, all of which are important for predicting adsorption. The simulation results are validated with high-resolution experimental measurements of isotherms and microcalorimetric heats of adsorption on well-characterized materials. The choice of first-principles method has a significant influence on the ability of force fields to accurately describe CO2-zeolite interactions. The PBE-D2 derived force field, which performed well for CO2 adsorption in siliceous zeolites, does not do so for Na-exchanged zeolites; the PBE-D2 method overestimates CO2 adsorption energies on multi-cation sites that are common in cationic zeolites with low Si/Al ratios. In contrast, a force field derived from the DFT/CC method performed well. Agreement was obtained between simulation and experiment not only for LTA-4A on which the force field fitting is based, but for other two common adsorbents, NaX and NaY.

Experimental and computational studies of pyridine-assisted post-synthesis modified air stable covalent-organic frameworks.

Pyridine-modified COF-10 exhibits enhanced stability in humid air relative to un-modified COF-10. Solid state NMR and computational studies were used to probe the nature of pyridine interactions with the framework. We propose two models for pyridine-framework interactions with different stabilities.

Positive curvature effects and interparticle capillary condensation during nitrogen adsorption in particulate porous materials.

The adsorption of nitrogen in a collection of spheres that touch or merge in a sintering-like manner is modeled using a Derjaguin-Broeckhof-de Boer approach. The proposed model accounts for both positive curvature effects and for capillary condensation at the contact between two spheres. A methodology is proposed to fit the P/P(0)>0.4 adsorption region with the coordination number of the spheres as the only adjustable parameter. The use of the model is illustrated on a series of silica aerogels. The suitability of various standard isotherms needed for the modeling is also discussed.

Density functional theory model of adsorption deformation.

Molecules adsorbed in pores cause elastic deformations of the solid matrix leading to either contraction or swelling of the material. Although experimental manifestation of adsorption-induced deformation in clays, coals, carbons, silicas, and other materials has been known for a long time, a rigorous theoretical description of this phenomenon is lacking. We report the nonlocal density functional theory (NLDFT) calculations that reproduce almost quantitatively the adsorption and strain isotherms of Kr and Xe on zeolite X. This system exhibits characteristic contraction at low vapor pressures and swelling at high vapor pressures. We show that the experimentally observed changes in the adsorbent volume are proportional to the solvation (disjoining) pressure caused by the adsorption stress exerted on the pore walls. The proposed NLDFT model can be used for the interpretation of adsorption measurements in micro- and mesoporous materials and for the characterization of their mechanical properties.

Density functional theory model of adsorption on amorphous and microporous silica materials.

We present a novel quenched solid density functional theory (QSDFT) model of adsorption on heterogeneous surfaces and porous solids, which accounts for the effects of surface roughness and microporosity. Within QSDFT, solid atoms are considered as quenched component(s) of the solid-fluid system with given density distribution(s). Solid-fluid intermolecular interactions are split into hard-sphere repulsive and mean-field attractive parts. The former are treated with the multicomponent fundamental measure density functional. Capabilities of QSDFT are demonstrated by drawing on the example of adsorption on amorphous silica materials. We show that, using established intermolecular potentials and a realistic model for silica surfaces, QSDFT quantitatively describes adsorption/desorption isotherms of Ar and Kr on reference MCM-41, SBA-15, and LiChrosphere materials in a wide range of relative pressures. QSDFT offers a systematic approach to the practical problems of characterization of microporous, mesoporous, and amorphous silica materials, including an assessment of microporosity, surface roughness, and adsorption deformation. Predictions for the pore diameter and the extent of pore surface roughness in MCM-41 and SBA-15 materials are in very good agreement with recent X-ray diffraction studies.

Characterization of micro-mesoporous materials from nitrogen and toluene adsorption: experiment and modeling.

Universal mechanisms of adsorption and capillary condensation of toluene and nitrogen on ordered MCM-41 and PHTS materials are studied by means of high-resolution experiments and Monte Carlo molecular simulations. A molecular simulation model of toluene adsorption in silica nanopores, which accounts for surface heterogeneity, and a hybrid molecular-macsroscopic method for pore size distribution (PSD) calculations have been developed. For a range of reference materials, the PSD results obtained from toluene isotherms are consistent with the results of nitrogen adsorption using the nonlocal density functional theory method.

Adsorption hysteresis of nitrogen and argon in pore networks and characterization of novel micro- and mesoporous silicas.

We report results of nitrogen and argon adsorption experiments performed at 77.4 and 87.3 K on novel micro/mesoporous silica materials with morphologically different networks of mesopores embedded into microporous matrixes: SE3030 silica with worm-like cylindrical channels of mode diameter of approximately 95 angstroms, KLE silica with cage-like spheroidal pores of ca. 140 angstroms, KLE/IL silica with spheroidal pores of approximately 140 angstroms connected by cylindrical channels of approximately 26 angstroms, and, also for a comparison, on Vycor glass with a disordered network of pores of mode diameter of approximately 70 angstroms. We show that the type of hysteresis loop formed by adsorption/desorption isotherms is determined by different mechanisms of condensation and evaporation and depends upon the shape and size of pores. We demonstrate that adsorption experiments performed with different adsorptives allow for detecting and separating the effects of pore blocking/percolation and cavitation in the course of evaporation. The results confirm that cavitation-controlled evaporation occurs in ink-bottle pores with the neck size smaller than a certain critical value. In this case, the pressure of evaporation does not depend upon the neck size. In pores with larger necks, percolation-controlled evaporation occurs, as observed for nitrogen (at 77.4 K) and argon (at 87.3 K) on porous Vycor glass. We elaborate a novel hybrid nonlocal density functional theory (NLDFT) method for calculations of pore size distributions from adsorption isotherms in the entire range of micro- and mesopores. The NLDFT method, applied to the adsorption branch of the isotherm, takes into account the effect of delayed capillary condensation in pores of different geometries. The pore size data obtained by the NLDFT method for SE3030, KLE, and KLE/IL silicas agree with the data of SANS/SAXS techniques.

Using nitrogen and carbon dioxide molecules to probe arsenic(V) bioaccessibility in soils.

Highly specialized personnel and high cost are typically required for in vivo risk assessment of arsenic (As) exposure to humans in As-contaminated soils. Arsenic bioaccessibility in soils, as determined with the aid of in vitro tests, is quite variable, and its magnitude depends upon unidentified soil properties. Use of soil chemical properties is a common practice for construction of As(V) sorption and bioaccessibility models with relative success. We propose a novel As(V) bioaccessibility model, which was tested on 17 soils. The model includes only two parameters characterizing surface properties of soils that are readily determined from N2- and CO2-based specific surface areas (SSAs), and total organic carbon (OC) content. We found that N2 and CO2 molecules act as As(V) "surrogates", probing easily accessible and relatively difficult to access soil porosity, respectively. Three interrelated linear models were constructed using two terms (CO2/N2-based SSAs and OC) that were significant (p <0.001) in explaining 51 and 95% of the variability observed in As(V) sorption and bioaccessibility, respectively. The proposed models successfully predicted bioaccessible As concentrations for 4 out of the 5 soils that were not included in the bioaccessibility models, reaching RMSE values of < or =10%.

Nitrogen and carbon dioxide adsorption by soils.

High-resolution nitrogen (77 K) and carbon dioxide (273 K) adsorption at subatmospheric pressures has been studied for a range of model soils of various origins with different organic matter (OM) contents. It is demonstrated that N2 and CO2 molecules probe different regions of soil particles. Nitrogen is adsorbed primarily on the outer surface of soil particles, while CO2 has a higher affinity to OM domains. Low-pressure nitrogen adsorption reveals that soil particle surfaces consist of clay/mineral domains with discrete patches of OM. A linear correlation has been found between the CO2 uptake and the amount of organic carbon reduced per unit of the external surface area. A new method for discriminating the microporosity of soil particles and accessibility of OM has been proposed.

"Humic coverage index" as a determining factor governing strain-specific hydrocarbon availability to contaminant-degrading bacteria in soils.

We report development of a novel parameter for quantifying the amount of humic and fulvic acids per unit surface area in a particular soil. This quantity, the "humic coverage index" (HCI), provides a measurement of the relative spatial extents and/or thicknesses of the humic/fulvic overlayers in different soils, and, therefore, can be used in modeling various soils' behavior in sequestration processes in which humic materials are involved. HCI is herein applied to modeling biodegradation of aromatic and aliphatic hydrocarbons (phenanthrene, pyrene, and hexadecane) by several bacterial strains. Results indicate that, for the cases studies here, contaminant biodegradation is highest at a particular HCI and decreases if the coverage density of humic material is lower or higher than this optimum value. The HCI value at which maximal degradation was observed varied across different strains (indicating strain-specific differences in ability to degrade contaminants sorbed to humic materials) and, to a lesser extent, across different contaminants. The HCI concept is also demonstrated to be useful in explaining soil-, strain-, and contaminant-specific variations in the ability of fulvic acid supplementation to enhance contaminant biodegradation. Finally, we show that, in general, strains which are comparatively better at degrading contaminants in high-HCI soils also show enhanced contaminant mineralization in vitro in the presence of humic acids, such as when hydrocarbons are adsorbed onto these materials.

Sorption hysteresis of benzene in charcoal particles.

Charcoal is found in water, soil, and sediment where it may act as a sorbent of organic pollutants. The sorption of organic compounds to natural solids often shows hysteresis. The purpose of this study was to determine the source of pronounced hysteresis that we found in the sorption of a hydrophobic compound (benzene) in water to a maple-wood charcoal prepared by oxygen-limited pyrolysis at 673 K. Gas adsorption (N2, Ar, CO2), 13C NMR, and FTIR show the charcoal to be a microporous solid composed primarily of elemental (aromatic) C and secondarily of carboxyl and phenolic C. Nonlocal density functional theory (N2, Ar) and Monte Carlo (CO2) calculations reveal a porosity of 0.15 cm3/g, specific surface area of 400 m2/g, and appreciable porosity in ultramicropores < 10 A. Benzene sorption-desorption conditions were chosen to eliminate artificial causes of hysteresis (rate-limiting diffusion, degradation, colloids effect). Charcoal sorbed up to its own weight of benzene at approximately 69% of benzene water solubility. Sorption was highly irreversible over most of the range tested (10(-4)-10(3) microg/mL). A dimensionless irreversibility index (/i) (0 < or = /i < or = 1) based on local slopes of adsorption and desorption branches was evaluated at numerous places along the isotherm. /i decreases as C increases, from 0.9-1 at low concentration to approximately 0 (approximately fully reversible) at the highest concentrations. Using sedimentation and volumetric displacement measurements, benzene is observed to cause pronounced swelling (up to > 2-fold) of the charcoal particles. It is proposed that hysteresis is due to pore deformation by the solute, which results in the pathway of sorption being different than the pathway of desorption and which leads to entrapment of some adsorbate as the polyaromatic scaffold collapses during desorption. It is suggested that intra-charcoal mass transport may be influenced by structural rearrangement of the solid, in addition to molecular diffusion.

Inside the hysteresis loop: multiplicity of internal states in confined fluids.

We study the equilibrium and stability of metastable states and capillary condensation hysteresis of a Lennard-Jones fluid in cylindrical pores by means of the canonical ensemble density functional theory and gauge cell Monte Carlo simulations. We demonstrate a possibility for the existence of multiple laterally uniform internal states of equal density inside the hysteresis loop. The region of multiple states is bounded by the states of zero compressibility. The internal states can be stabilized in Monte Carlo simulations constraining the density fluctuations.