Connecting theory and simulation with experiment for the study of diffusion in nanoporous solids

Nanoporous solids are ubiquitous in chemical, energy, and environmental processes, where controlled transport of molecules through the pores plays a crucial role. They are used as sorbents, chromatographic or membrane materials for separations, and as catalysts and catalyst supports. Defined as materials where confinement effects lead to substantial deviations from bulk diffusion, nanoporous materials include crystalline microporous zeotypes and metal–organic frameworks (MOFs), and a number of semi-crystalline and amorphous mesoporous solids, as well as hierarchically structured materials, containing both nanopores and wider meso- or macropores to facilitate transport over macroscopic distances. The ranges of pore sizes, shapes, and topologies spanned by these materials represent a considerable challenge for predicting molecular diffusivities, but fundamental understanding also provides an opportunity to guide the design of new nanoporous materials to increase the performance of transport limited processes. Remarkable progress in synthesis increasingly allows these designs to be put into practice. Molecular simulation techniques have been used in conjunction with experimental measurements to examine in detail the fundamental diffusion processes within nanoporous solids, to provide insight into the free energy landscape navigated by adsorbates, and to better understand nano-confinement effects. Pore network models, discrete particle models and synthesis-mimicking atomistic models allow to tackle diffusion in mesoporous and hierarchically structured porous materials, where multiscale approaches benefit from ever cheaper parallel computing and higher resolution imaging. Here, we discuss synergistic combinations of simulation and experiment to showcase theoretical progress and computational techniques that have been successful in predicting guest diffusion and providing insights. We also outline where new fundamental developments and experimental techniques are needed to enable more accurate predictions for complex systems

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 performed 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 utilized 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 non-polar 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

New insights into the breathing phenomenon in ZIF-4

Structural changes in ZIFs upon adsorption remain a paradigm due to the sensitivity of the adsorption mechanism to the nature of the organic ligands and gas probe molecules. Synchrotron X-ray diffraction under operando conditions clearly demonstrates for the first time that ZIF-4 exhibits a structural reorientation from a narrow-pore (np) to a new expanded-pore (ep) structure upon N2 adsorption, while it does not do so for CO2 adsorption. The existence of an expanded-pore structure of ZIF-4 has also been predicted by molecular simulations. In simulations the expanded structure was stabilized by entropy at high temperatures and by strong adsorption of N2 at low temperatures. These results are in perfect agreement with manometric adsorption measurements for N2 at 77 K that show the threshold pressure for breathing at ∼30 kPa. Inelastic neutron scattering (INS) measurements show that CO2 is also able to promote structural changes but, in this specific case, only at cryogenic temperatures (5 K).The authors would like to acknowledge financial support from the MINECO (MAT2016-80285-p), Generalitat Valenciana (PROMETEOII/2014/004), H2020 (MSCA-RISE-2016/NanoMed Project), Spanish ALBA synchrotron (Projects AV-2017021985 and IH-2018012591) and Oak Ridge beam time availability (Project IPTS-20843.1). JSA and JGL acknowledge financial support from UA (ACIE17-15) to cover all the expenses for INS measurements at Oak Ridge. JGL acknowledges GV (GRISOLIAP/2016/089) for the research contract

Adsorption hysteresis of nitrogen and argon in pore networks and characterization of novel microand 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 wormlike cylindrical channels of mode diameter of ∼95 Å, KLE silica with cagelike spheroidal pores of ca. 140 Å, KLE/IL silica with spheroidal pores of ∼140 Å connected by cylindrical channels of ∼26 Å, and, also for a comparison, on Vycor glass with a disordered network of pores of mode diameter of ∼70 Å. 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

Proton dynamics in phosphotungstic acid impregnated mesoporous silica proton exchange membrane materials

Phosphotungstic acid is an excellent proton conductor that can be incorporated into porous supports, and nanocomposite proton exchange membrane materials made from mesoporous silica impregnated with phosphotungstic acid have been suggested for use in fuels cells operating > 100 °C. In this work, quasielastic neutron scattering was used to study proton self-diffusion in mesoporous disordered and P6mm symmetry silica impregnated with two concentrations of phosphotungstic acid. Overall, the silica structure had a significantly greater effect on proton conduction and diffusion than phosphotungstic acid concentration, with higher proton conduction occurring for the P6mm symmetry silica samples. Quasielastic neutron scattering revealed two populations of protons diffusing through each sample, and that proton conduction is limited by the slower of these populations, which diffuse via a jump-diffusion mechanism. Whilst the fundamental jump-diffusion mechanism by which these slower protons moved was found to be similar for both silica supports and phosphotungstic acid concentrations, the faster diffusion occurring in P6mm structured silica arises from a lower residence time of protons moving between sites in the jump-diffusion model, suggesting a lower energy barrier

A Metal–Organic Framework with a Hierarchical System of Pores and Tetrahedral Building Blocks

No AbstractPeer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/50670/1/2528_ftp.pd

Capillary Condensation and Interface Structure of a Model Colloid-Polymer Mixture in a Porous Medium

We consider the Asakura-Oosawa model of hard sphere colloids and ideal polymers in contact with a porous matrix modeled by immobilized configurations of hard spheres. For this ternary mixture a fundamental measure density functional theory is employed, where the matrix particles are quenched and the colloids and polymers are annealed, i.e. allowed to equilibrate. We study capillary condensation of the mixture in a tiny sample of matrix as well as demixing and the fluid-fluid interface inside a bulk matrix. Density profiles normal to the interface and surface tensions are calculated and compared to the case without matrix. Two kinds of matrices are considered: (i) colloid-sized matrix particles at low packing fractions and (ii) large matrix particles at high packing fractions. These two cases show fundamentally different behavior and should both be experimentally realizable. Furthermore, we argue that capillary condensation of a colloidal suspension could be experimentally accessible. We find that in case (ii), even at high packing fractions, the main effect of the matrix is to exclude volume and, to high accuracy, the results can be mapped onto those of the same system without matrix via a simple rescaling.Comment: 12 pages, 9 figures, submitted to PR

Combining mercury thermoporometry with integrated gas sorption and mercury porosimetry to improve accuracy of pore-size distributions for disordered solids

The typical approach to analysing raw data, from common pore characterization methods such as gas sorption and mercury porosimetry, to obtain pore size distributions for disordered porous solids generally makes several critical assumptions that impact the accuracy of the void space descriptors thereby obtained. These assumptions can lead to errors in pore size of as much as 500%. In this work, we eliminated these assumptions by employing novel experiments involving fully integrated gas sorption, mercury porosimetry and mercury thermoporometry techniques. The entrapment of mercury following porosimetry allowed the isolation (for study) of a particular subset of pores within a much larger interconnected network. Hence, a degree of specificity of findings to particular pores, more commonly associated with use of templated, model porous solids, can also be achieved for disordered materials. Gas sorption experiments were conducted in series, both before and after mercury porosimetry, on the same sample, and the mercury entrapped following porosimetry was used as the probe fluid for theromporometry. Hence, even if one technique, on its own, is indirect, requiring unsubstantiated assumptions, the fully integrated combination of techniques described here permits the validation of assumptions used in one technique by another. Using controlled-pore glasses as model materials, mercury porosimetry scanning curves were used to establish the correct correspondence between the appropriate Gibbs–Thomson parameter, and the nature of the meniscus geometry in melting, for thermoporometry measurements on entrapped mercury. Mercury thermoporometry has been used to validate the pore sizes, for a series of sol–gel silica materials, obtained from mercury porosimetry data using the independently-calibrated Kloubek correlations. The pore sizes obtained for sol–gel silicas from porosimetry and thermoporometry have been shown to differ substantially from those obtained via gas sorption and NLDFT analysis. DRIFTS data for the samples studied has suggested that the cause of this discrepancy may arise from significant differences in the surface chemistries between the samples studied here and that used to calibrate the NLDFT potentials

Molecular simulation of carbon dioxide adsorption in chemically and structurally heterogeneous porous carbons

Capture of carbon dioxide from fossil fuel power plants via adsorption and sequestration of carbon dioxide in unmineable coal seams are achievable near-term methods of reducing atmospheric emissions of this greenhouse gas. To investigate the influence of surface heterogeneity upon predicted adsorption behavior in activated carbons and coal, isotherms were generated via grand canonical Monte Carlo simulation for CO 2 adsorption in slit-shaped pores with underlying graphitic structure and several variations of chemical heterogeneity (oxygen and hydrogen content), pore width, and surface functional group orientation. Adsorption generally increased with increasing surface oxygen content, although exceptions to this trend were observed on structurally heterogeneous surfaces with holes or furrows that yield strongly adsorbing preferred binding sites. Among the heterogeneous pore structures investigated, those with coal-like surfaces adsorbed carbon dioxide more strongly than planar, homogeneous graphitic slit pores of comparable width. Electrostatic adsorbate–adsorbent interactions significantly influenced adsorption onto model surfaces. © 2006 American Institute of Chemical Engineers Environ Prog, 25: 343–354, 2006Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/55883/1/10168_ftp.pd