Investigating cluster formation in adsorption of CO2, CH4, and Ar in zeolites and metal organic frameworks at subcritical temperatures

Mass transfer of toluene in a series of metal-organic frameworks: molecular clusters inside the nanopores cause slow and step-like release

The mass transfer of the guest molecules in the pores is fundamental for the application of nanoporous materials like metal-organic frameworks, MOFs. In the present work, we explore the uptake and release of toluene in a series of Zr-based MOFs with different pore sizes. We find that intermolecular guest-guest interaction, sterically controlled by the pore size, has a substantial impact on the release kinetics. While the adsorption is rather fast, the desorption process is many orders of magnitude slower. Depending on the pore size, molecular clusters form, here (most likely) toluene dimers, which are rather stable and their break-up is rate-limiting during the desorption process. This results in a step-like desorption kinetics, deviating from the plain Fickian-diffusion-controlled release. Temperature-dependent experiments show that the minimum and maximum of the release rates are obtained at the same toluene loadings, independent of the temperature. Moreover, the activation energy for the release coincides with the binding energy of a toluene dimer. The work shows the importance of intermolecular guest-guest interaction, controlled by the MOF-nanoconfinement, for the uptake and release from nanoporous materials.

How Reliable Is the Ideal Adsorbed Solution Theory for the Estimation of Mixture Separation Selectivities in Microporous Crystalline Adsorbents?

Microporous crystalline adsorbents such as zeolites and metal-organic frameworks (MOFs) have potential use in a wide variety of separation applications. The adsorption selectivity S _ads is a key metric that quantifies the efficacy of any microporous adsorbent in mixture separations. The Ideal Adsorbed Solution Theory (IAST) is commonly used for estimating the value of S _ads, with unary isotherms of the constituent guests as data inputs. There are two basic tenets underlying the development of the IAST. The first tenet mandates a homogeneous distribution of adsorbates within the pore landscape. The second tenet requires the surface area occupied by a guest molecule in the mixture to be the same as that for the corresponding pure component. Configurational-bias Monte Carlo (CBMC) simulations are employed in this article to highlight several scenarios in which the IAST fails to provide a quantitatively correct description of mixture adsorption equilibrium due to a failure to conform to either of the two tenets underpinning the IAST. For CO₂ capture with cation-exchanged zeolites and MOFs with open metal sites, there is congregation of CO₂ around the cations and unsaturated metal atoms, resulting in failure of the IAST due to an inhomogeneous distribution of adsorbates in the pore space. Thermodynamic non-idealities also arise due to the preferential location of CO₂ molecules at the window regions of 8-ring zeolites such as DDR and CHA or within pockets of MOR and AFX zeolites. Thermodynamic non-idealities are evidenced for water/alcohol mixtures due to molecular clustering engendered by hydrogen bonding. It is also demonstrated that thermodynamic non-idealities can be strong enough to cause selectivity reversals, which are not anticipated by the IAST.

Ca-Embedded C2N: an efficient adsorbent for CO2 capture

Carbon dioxide as a greenhouse gas causes severe impacts on the environment, whereas it is also a necessary chemical feedstock that can be converted into carbon-based fuels via electrochemical reduction. To efficiently and reversibly capture CO₂, it is important to find novel materials for a good balance between adsorption and desorption. In this study, we performed first-principles calculations and grand canonical Monte Carlo (GCMC) simulations, to systematically study metal-embedded carbon nitride (C₂N) nanosheets for CO₂ capture. Our first-principles results indicated that Ca atoms can be uniformly trapped in the cavity center of C₂N structure, while the transition metals (Sc, Ti, V, Cr, Mn, Fe, Co) are favorably embedded in the sites off the center of the cavity. The determined maximum number of CO₂ molecules with strong physisorption showed that Ca-embedded C₂N monolayer is the most promising CO₂ adsorbent among all considered metal-embedded materials. Moreover, GCMC simulations revealed that at room temperature the gravimetric density for CO₂ adsorbed on Ca-embedded C₂N reached 50 wt% at 30 bar and 23 wt% at 1 bar, higher than other layered materials, thus providing a satisfactory system for the CO₂ capture and utilization.

Sizable dynamics in small pores: CO2 location and motion in the α-Mg formate metal-organic framework

Metal-organic frameworks (MOFs) are promising materials for carbon dioxide (CO₂) adsorption and storage; however, many details regarding CO₂ dynamics and specific adsorption site locations within MOFs remain unknown, restricting the practical uses of MOFs for CO₂ capture. The intriguing α-magnesium formate (α-Mg₃(HCOO)₆) MOF can adsorb CO₂ and features a small pore size. Using an intertwined approach of ¹³C solid-state NMR (SSNMR) spectroscopy, ¹H-¹³C cross-polarization SSNMR, and computational molecular dynamics (MD) simulations, new physical insights and a rich variety of information have been uncovered regarding CO₂ adsorption in this MOF, including the surprising suggestion that CO₂ motion is restricted at elevated temperatures. Guest CO₂ molecules undergo a combined localized rotational wobbling and non-localized twofold jumping between adsorption sites. MD simulations and SSNMR experiments accurately locate the CO₂ adsorption sites; the mechanism behind CO₂ adsorption is the distant interaction between the hydrogen atom of the MOF formate linker and a guest CO₂ oxygen atom, which are ca. 3.2 Å apart.

Nanoporous Materials Can Tune the Critical Point of a Pure Substance

Molecular simulations and NMR relaxometry experiments demonstrate that pure benzene or xylene confined in isoreticular metal–organic frameworks (IRMOFs) exhibit true vapor–liquid phase equilibria where the effective critical point may be reduced by tuning the structure of the MOF. Our results are consistent with vapor and liquid phases extending over many MOF unit cells. These results are counterintuitive since the MOF pore diameters are approximately the same length scale as the adsorbate molecules. As applications of these materials in catalysis, separations, and gas storage rely on the ability to tune the properties of adsorbed molecules, we anticipate that the ability to systematically control the critical point, thereby preparing spatially inhomogeneous local adsorbate densities, could add a new design tool for MOF applications.

Separation of benzene from mixtures with water, methanol, ethanol, and acetone: highlighting hydrogen bonding and molecular clustering influences in CuBTC

Configurational-bias Monte Carlo (CBMC) simulations are used to establish the potential of CuBTC for separation of water/benzene, methanol/benzene, ethanol/benzene, and acetone/benzene mixtures. For operations under pore saturation conditions, the separations are in favor of molecules that partner benzene; this is due to molecular packing effects that disfavor benzene. CBMC simulations for adsorption of quaternary water/methanol/ethanol/benzene mixtures show that water can be selectively adsorbed at pore saturation, making CuBTC effective in drying applications. Ideal Adsorbed Solution Theory (IAST) calculations anticipate the right hierarchy of component loadings but the quantitative agreement with CBMC mixture simulations is poor for all investigated mixtures. The failure of the IAST to provide reasonable quantitative predictions of mixture adsorption is attributable to molecular clustering effects that are induced by hydrogen bonding between water-water, methanol-methanol, and ethanol-ethanol molecule pairs. There is, however, no detectable hydrogen bonding between benzene and partner molecules in the investigated mixtures. As a consequence of molecular clustering, the activity coefficients of benzene in the mixtures is lowered below unity by one to three orders of magnitude at pore saturation; such drastic reductions cannot be adequately captured by the Wilson model, that does not explicitly account for molecular clustering. Molecular clustering effects are also shown to influence the loading dependence of the diffusivities of guest molecules.

Formation of a quasi-solid structure by intercalated noble gas atoms in pores of CuI-MFU-4l metal–organic framework

Ten crystallographically different positions for Xe and eight positions for Kr form a quasi-solid structures within the large-pore metal–organic framework Cu^I-MFU-4l.

Molybdenum disulfide as a highly efficient adsorbent for non-polar gases

Molybdenum disulfide with vacancies is predicted to be an efficient absorbent for CO₂ and CH₄ molecules, and the gravimetric adsorption can be as high as 42 wt%.

Diffusion of interacting particles in discrete geometries: Equilibrium and dynamical properties

We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.110601]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self- and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self- and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

First principle chemical kinetics in zeolites: the methanol-to-olefin process as a case study

To optimally design next generation catalysts a thorough understanding of the chemical phenomena at the molecular scale is a prerequisite. Apart from qualitative knowledge on the reaction mechanism, it is also essential to be able to predict accurate rate constants. Molecular modeling has become a ubiquitous tool within the field of heterogeneous catalysis. Herein, we review current computational procedures to determine chemical kinetics from first principles, thus by using no experimental input and by modeling the catalyst and reacting species at the molecular level. Therefore, we use the methanol-to-olefin (MTO) process as a case study to illustrate the various theoretical concepts. This process is a showcase example where rational design of the catalyst was for a long time performed on the basis of trial and error, due to insufficient knowledge of the mechanism. For theoreticians the MTO process is particularly challenging as the catalyst has an inherent supramolecular nature, for which not only the Brønsted acidic site is important but also organic species, trapped in the zeolite pores, must be essentially present during active catalyst operation. All these aspects give rise to specific challenges for theoretical modeling. It is shown that present computational techniques have matured to a level where accurate enthalpy barriers and rate constants can be predicted for reactions occurring at a single active site. The comparison with experimental data such as apparent kinetic data for well-defined elementary reactions has become feasible as current computational techniques also allow predicting adsorption enthalpies with reasonable accuracy. Real catalysts are truly heterogeneous in a space- and time-like manner. Future theory developments should focus on extending our view towards phenomena occurring at longer length and time scales and integrating information from various scales towards a unified understanding of the catalyst. Within this respect molecular dynamics methods complemented with additional techniques to simulate rare events are now gradually making their entrance within zeolite catalysis. Recent applications have already given a flavor of the benefit of such techniques to simulate chemical reactions in complex molecular environments.

Diffusion of interacting particles in discrete geometries

We evaluate the self-diffusion and transport diffusion of interacting particles in a discrete geometry consisting of a linear chain of cavities, with interactions within a cavity described by a free-energy function. Exact analytical expressions are obtained in the absence of correlations, showing that the self-diffusion can exceed the transport diffusion if the free-energy function is concave. The effect of correlations is elucidated by comparison with numerical results. Quantitative agreement is obtained with recent experimental data for diffusion in a nanoporous zeolitic imidazolate framework material, ZIF-8.

Influence of adsorption thermodynamics on guest diffusivities in nanoporous crystalline materials

Published experimental data, underpinned by molecular simulations, are used to highlight the strong influence of adsorption thermodynamics on diffusivities of guest molecules inside ordered nanoporous crystalline materials such as zeolites, metal-organic frameworks (MOFs), and zeolitic imidazolate frameworks (ZIFs). For cage-type structures (e.g. LTA, CHA, DDR, and ZIF-8), the variation of the free energy barrier for inter-cage hopping across the narrow windows, -δFi, provides a rationalization of the observed strong influence of pore concentrations, ci, on diffusivities. In open structures with large pore volumes (e.g. FAU, IRMOF-1, CuBTC) and within channels (MFI, BEA, MgMOF-74, MIL-47, MIL-53), the pore concentration (ci) dependence of the self- (Di,self), Maxwell-Stefan (Đi), and Fick (Di) diffusivities are often strongly dictated by the inverse thermodynamic correction factor, 1/Γi≡∂ln ci/∂ln pi; the magnitudes of the diffusivities are dictated by the binding energies for adsorption. For many guest-host combinations Đi-ci dependence is directly related to the 1/Γivs. ci variation. When molecular clustering occurs, we get 1/Γi > 1, causing unusual Đivs. ci dependencies. The match, or mis-match, between the periodicity of the pore landscape and the conformations of adsorbed chain molecules often leads to non-monotonic variation of diffusivities with chain lengths.

X-Nuclei NMR Self-Diffusion Studies in Mesoporous Silica Foam and Microporous MOF CuBTC

A standard X-observe NMR probe was equipped with a z-gradient coil to enable high-sensitivity pulsed field gradient NMR diffusion studies of Li+ and Cs+ cations of aqueous salt solutions in a high-porosity mesocellular silica foam (MCF) and of CO2 adsorbed in metal-organic frameworks (MOF). The coil design and the necessary probe modifications, which yield pulsed field gradients of up to ±16.2Tm-1, are introduced. The system was calibrated at 2H resonance frequency and successfully applied for diffusion studies at 7Li, 23Na, 13C and 133Cs frequencies. Significant reductions of the diffusivities of the cations in LiClac and CsClac solution introduced into MCFs are observed. By comparison of the diffusion behavior with the bulk solutions, a tortuosity of the silica foam of 4.5±0.6 was derived. Single component self-diffusion of CO2 and CH4 (measured by 1H NMR) as well as self-diffusion of the individual components in CO2/CH4 mixtures was studied in the MOF CuBTC. The experimental results confirm high mobilities of the adsorbed gases and trends for diffusion separation factors predicted by MD simulations.

An experimental and molecular simulation study of the adsorption of carbon dioxide and methane in nanoporous carbons in the presence of water

The adsorption of carbon dioxide and methane in nanoporous carbons in the presence of water is studied using experiments and molecular simulations. For all amounts of adsorbed water molecules, the adsorption isotherms for carbon dioxide and methane resemble those obtained for pure fluids. The pore filling mechanism does not seem to be affected by the presence of the water molecules. Moreover, the pressure at which the maximum adsorbed amount of methane or carbon dioxide is reached is nearly insensitive to the loading of preadsorbed water molecules. In contrast, the adsorbed amount of methane or carbon dioxide decreases linearly with the number of guest water molecules. Typical molecular configurations obtained using molecular simulation indicate that the water molecules form isolated clusters within the host porous carbon due to the nonfavorable interaction between carbon dioxide or methane and water.

Metal-organic frameworks: a rapidly growing class of versatile nanoporous materials

Metal-organic frameworks (MOFs) represent a new class of hybrid organic-inorganic supramolecular materials comprised of ordered networks formed from organic electron donor linkers and metal cations. They can exhibit extremely high surface areas, as well as tunable pore size and functionality, and can act as hosts for a variety of guest molecules. Since their discovery, MOFs have enjoyed extensive exploration, with applications ranging from gas storage to drug delivery to sensing. This review covers advances in the MOF field from the past three years, focusing on applications, including gas separation, catalysis, drug delivery, optical and electronic applications, and sensing. We also summarize recent work on methods for MOF synthesis and computational modeling.

Hydrogen bonding effects in adsorption of water-alcohol mixtures in zeolites and the consequences for the characteristics of the Maxwell-Stefan diffusivities

This work highlights a variety of peculiar characteristics of adsorption and diffusion of polar molecules such as water, methanol and ethanol in zeolites. These peculiarities are investigated with the aid of configurational-bias Monte Carlo (CBMC) simulations of adsorption isotherms, and molecular dynamics (MD) simulations of diffusivities in FAU, MFI, DDR, and LTA zeolites. Because of strong hydrogen bonding, significant clustering of the guest molecules occurs in all investigated structures. Because of molecular clustering, the inverse thermodynamic factor 1/Gamma(i) identical with (d[ln c(i)])/(d[ln f(i)]) exceeds unity for a range molar concentrations c(i) within the micropores. The degree of clustering is lowered as the temperature is increased. For the concentration ranges for which 1/Gamma(i) > 1, the Fick diffusivity, D(i), for unary diffusion is often lower than both the Maxwell-Stefan, D(i), and the self-diffusivity, D(i,self). For water-alcohol mixtures, the hydrogen bonding between water and alcohol molecules is much more predominant than for water-water, and alcohol-alcohol molecule pairs. Consequently, the adsorption of water-alcohol mixtures shows significant deviations from the predictions of the ideal adsorbed solution theory (IAST). The water-alcohol bonding also leaves its imprint on the mixture diffusion characteristics. The Maxwell-Stefan diffusivity, D(i), of either component in water-alcohol mixtures is lower than the corresponding values of the pure components; this behavior is distinctly different from that for mixtures of nonpolar guest molecules. The binary exchange coefficient D(12) for water-alcohol mixtures is also significantly lower than either self-exchange coefficients D(11) and D(22) of the constituent species. This implies that correlation effects are significantly stronger in water-alcohol mixtures than for the constituent species. Correlation effects are found to be significant for water-alcohol mixture diffusion in DDR and LTA zeolites, even though such effects are negligible for the pure constituents. The major conclusion to emerge from this investigation is that, unlike mixtures of nonpolar molecules, it is not possible to estimate water-alcohol mixture adsorption and diffusion characteristics on the basis of pure component data.

Cobalt doping of the MOF-5 framework and its effect on gas-adsorption properties

Partial isomorphic substitution of Zn in IRMOF metal clusters by cobalt ions is described for the first time. Specifically, different numbers of Co(2+) ions have been incorporated during solvothermal crystallization into the Zn-based MOF-5 (IRMOF-1) framework, which is one of the most studied MOF materials. The amount of Zn that can be substituted seems to be limited, being no more than 25% of total metal content, that is, no more than one Co atom inside every metal cluster formed by four transition-metal ions, on average. Several characterization techniques, including X-ray diffraction, DR UV-visible spectroscopy, N(2) adsorption isotherms, and thermogravimetrical analysis, strongly support the effective incorporation of Co into the material framework. As-synthesized CoMOF-5 has cobalt ions in octahedral coordination, changing to tetrahedral by simple evacuation, presumably by the removal of two diethylformamide molecules per Co ion. Moreover, the H(2), CH(4), and CO(2) uptake of MOF-5 materials systematically increases with the Co content, particularly at high pressure. Such an increase is moderate anyway, considering that Co is incorporated into unexposed metal sites that are less accessible to gas molecules.