Enthalpy–entropy correlation for hydrogen adsorption on zeolites

Quantum chemical modeling of hydrogen binding in metal-organic frameworks: validation, insight, predictions and challenges

A detailed chemical understanding of H2 interactions with binding sites in the nanoporous crystalline structure of metal-organic frameworks (MOFs) can lay a sound basis for the design of new sorbent materials. Computational quantum chemical calculations can aid in this quest. To set the stage, we review general thermodynamic considerations that control the usable storage capacity of a sorbent. We then discuss cluster modeling of H2 ligation at MOF binding sites using state-of-the-art density functional theory (DFT) calculations, and how the binding can be understood using energy decomposition analysis (EDA). Employing these tools, we illustrate the connections between the character of the MOF binding site and the associated adsorption thermodynamics using four experimentally characterized MOFs, highlighting the role of open metal sites (OMSs) in accessing binding strengths relevant to room temperature storage. The sorbents are MOF-5, with no open metal sites, Ni2(m-dobdc), containing Lewis acidic Ni(II) sites, Cu(I)-MFU-4l, containing π basic Cu(I) sites and V2Cl2.8(btdd), also containing π-basic V(II) sites. We next explore the potential for binding multiple H2 molecules at a single metal site, with thermodynamics useful for storage at ambient temperature; a materials design goal which has not yet been experimentally demonstrated. Computations on Ca2+ or Mg2+ bound to catecholate or Ca2+ bound to porphyrin show the potential for binding up to 4 H2; there is precedent for the inclusion of both catecholate and porphyrin motifs in MOFs. Turning to transition metals, we discuss the prediction that two H2 molecules can bind at V(II)-MFU-4l, a material that has been synthesized with solvent coordinated to the V(II) site. Additional calculations demonstrate binding three equivalents of hydrogen per OMS in Sc(I) or Ti(I)-exchanged MFU-4l. Overall, the results suggest promising prospects for experimentally realizing higher capacity hydrogen storage MOFs, if nontrivial synthetic and desolvation challenges can be overcome. Coupled with the unbounded chemical diversity of MOFs, there is ample scope for additional exploration and discovery.

Characterization of polar surface groups on siliceous materials by inverse gas chromatography and the enthalpy-entropy compensation effect

Surface-modified porous silica is a well-established composite material. To improve its embedding and application behavior, adsorption studies of various probe molecules have been performed using the technique of inverse gas chromatography (IGC). For this purpose, IGC experiments were carried out in the infinite dilution mode on macro-porous micro glass spheres before and after surface modification with (3-mercaptopropyl)trimethoxysilane. To provide information about the polar interactions between probe molecules and the silica surface, in particular, eleven polar molecules have been injected. In summary, the free surface energy for pristine silica ( γ S t o t a l = 229 mJ/m²) and for (3-mercaptopropyl)trimethoxysilane-modified silica ( γ S t o t a l = 135 mJ/m²) indicates a reduced wettability after surface modification. This is due to the reduction of the polar component of the free surface energy ( γ S S P ) from 191 mJ/m² to 105 mJ/m². Simultaneously, with the reduction of surface silanol groups caused by surface modification of silica and, therefore, the decrease in polar interactions, a substantial loss of Lewis acidity was observed by various IGC approaches. Experiments with all silica materials have been conducted at temperatures in the range from 90°C to 120°C to determine the thermodynamic parameters, such as adsorption enthalpy ( Δ H a d s ) and adsorption entropy ( Δ S a d s ), using the Arrhenius regression procedure evaluating the IGC data. With the help of the enthalpy-entropy compensation, two types of adsorption complexes are assumed between polar probe molecules and the silica surface because of different isokinetic temperatures. Identical adsorption complexes with an isokinetic temperature of 370°C have been assigned to alkanes and weakly interacting polar probes such as benzene, toluene, dichloromethane, and chloroform. Polar probe molecules with typical functional groups such as OH, CO, and CN, having the ability to form hydrogen bonds to the silica surface, exhibit a lower isokinetic temperature of 60°C. Quantum chemical calculations of the probe molecules on a non-hydroxylated and hydroxylated silica cluster supported the formation of hydrogen bonds in the case of a strong polar adsorption complex with a bonding distance of 1.7 nm-1.9 nm to the silica surface.

A Theoretical Investigation of Cu+, Ni2+ and Co2+-exchanged Zeolites for Hydrogen Storage

This study investigates the H 2 adsorption on Cu + , Ni 2+ and Co 2+ -exchanged SSZ-13 (CHA) and SSZ-39 (AEI) using periodic DFT computations. Most stable Cu + position was found to be the 6-membered-ring window for both zeolites. Similarly, for investigated Ni 2+ and Co 2+ loadings on 6-membered-ring windows, the third nearest neighbor Al positions, i.e., Al-O-Si-O-Si-O-Al coordination, was found to be the most stable position. H 2 adsorption was investigated for all the Cu + , Ni 2+ and Co 2+ centers. AEI and CHA resulted in similar H 2 -Cu interactions for the Al and B substituted structures. H 2 adsorption on Cu + located in the 8-membered-ring gave the highest adsorption energy for both frameworks. Replacing Al with B in the framework increased the electron back donation from Cu + (3d) orbitals to H 2 antibonding orbital (s H2 * ). The H 2 adsorption energies on the Ni 2+ and Co 2+ -exchanged zeolites were found to be between -15 and -44 kJ/mol. Higher energy values were observed on the AEI framework, especially when two Al atoms have the Al-O-Si-O-Al configuration. Lesser interaction of the d-orbitals in the case of the Co 2+ and Ni 2+ cations resulted in heat of H 2 adsorption close to optimum values required for H 2 storage on porous materials.

Assessment of Performance of Density Functionals for Predicting Potential Energy Curves in Hydrogen Storage Applications

The availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H₂) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 data set consisting of equilibrium geometries and assessed the performance of 55 density functionals over this data set (Veccham, S. P.; Head-Gordon, M. J. Chem. Theory Comput. 2020, 16, 4963-4982). As it is crucial for computational tools to accurately model the entire potential energy curve (PEC), in addition to the equilibrium geometry, we extended this data set with 389 new data points to include two compressed and three elongated geometries along 78 PECs for H₂ binding, forming the H2Bind78 × 7 data set. By assessing the performance of 55 density functionals on this significantly larger and more comprehensive H2Bind78 × 7 data set, we identified the best performing density functionals for H₂ binding applications: PBE0-DH, ωB97X-V, ωB97M-V, and DSD-PBEPBE-D3(BJ). The addition of Hartree-Fock exchange improves the performance of density functionals, albeit not uniformly throughout the PEC. We recommend the usage of ωB97X-V and ωB97M-V density functionals as they offer good performance for both geometries and energies. In addition, we also identified B97M-V and B97M-rV as the best semilocal density functionals for predicting H₂ binding energy at its equilibrium geometry.

Ambient-Temperature Hydrogen Storage via Vanadium(II)-Dihydrogen Complexation in a Metal-Organic Framework

The widespread implementation of H2 as a fuel is currently hindered by the high pressures or cryogenic temperatures required to achieve reasonable storage densities. In contrast, the realization of materials that strongly and reversibly adsorb hydrogen at ambient temperatures and moderate pressures could transform the transportation sector and expand adoption of fuel cells in other applications. To date, however, no adsorbent has been identified that exhibits a binding enthalpy within the optimal range of -15 to -25 kJ/mol for ambient-temperature hydrogen storage. Here, we report the hydrogen adsorption properties of the metal-organic framework (MOF) V2Cl2.8(btdd) (H2btdd, bis(1H-1,2,3-triazolo[4,5-b],[4',5'-i])dibenzo[1,4]dioxin), which features exposed vanadium(II) sites capable of backbonding with weak π acids. Significantly, gas adsorption data reveal that this material binds H2 with an enthalpy of -21 kJ/mol. This binding energy enables usable hydrogen capacities that exceed that of compressed storage under the same operating conditions. The Kubas-type vanadium(II)-dihydrogen complexation is characterized by a combination of techniques. From powder neutron diffraction data, a V-D2(centroid) distance of 1.966(8) Å is obtained, the shortest yet reported for a MOF. Using in situ infrared spectroscopy, the H-H stretch was identified, and it displays a red shift of 242 cm-1. Electronic structure calculations show that a main contribution to bonding stems from the interaction between the vanadium dπ and H2 σ* orbital. Ultimately, the pursuit of MOFs containing high densities of weakly π-basic metal sites may enable storage capacities under ambient conditions that far surpass those accessible with compressed gas storage.

A Step Forward in Understanding the Hydrogen Adsorption and Compression on Activated Carbons

Hydrogen adsorption on activated carbons (ACs) is a promising alternative to compression and liquefaction for storing hydrogen. Herein, we have studied hydrogen adsorption on six commercial ACs (CACs) with surface areas ranging from 996 to 2216 m² g^-1 in a temperature range of 77 to 273 K and pressures up to 15 MPa. Excess hydrogen adsorption capacities of 2.3 to 5.8 wt % were obtained at 77 K and 4 MPa. We demonstrated that, contrary to what is normally done, hydrogen capacity is more accurately predicted by the surface area determined by the nonlocal density functional theory method applied to N₂ and CO₂ adsorption data than by the Brunauer-Emmett-Teller (BET) area. The modified Dubinin-Astakhov (MDA) equation was used to fit the experimental adsorption data, and the relationship between the MDA parameters (n_max, V_a, α, and β) and the textural properties of the CACs was determined for the first time. We concluded that the n_max and V_a parameters are related to the BET area, while the α and β parameters are related to the average micropore size and total pore volume, respectively. α and β were used to evaluate the enthalpy and entropy of adsorption and we show that these parameters can be used to assess the best carbon for hydrogen storage or compression.

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H₂ interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

Non-Linear Enthalpy-Entropy Correlation for Nitrogen Adsorption in Zeolites

The thermodynamics of dinitrogen adsorption in faujasite-type zeolites, Na-Y, Ca-Y and Sr-Y, were investigated by means of variable-temperature infrared spectroscopy, a technique that affords determination of the standard adsorption enthalpy (ΔH⁰) and entropy (ΔS⁰) from an analysis of the IR spectra recorded over a range of temperatures. The results obtained, taken together with previously reported values for N₂ adsorption on protonic zeolites, revealed a non-linear correlation between ΔH⁰ and ΔS⁰. Implications of such a correlation for gas separation and purification by adsorption in porous solids are highlighted.

Catalysis-tunable Heusler alloys in selective hydrogenation of alkynes: A new potential for old materials

Heusler alloys (X 2 YZ) are well-established intermetallic compound materials in various fields because their function can be precisely adjusted by elemental substitution (e.g., X 2 YZ 1-x Z' x ). Although intermetallic compound catalysts started attracting attention recently, catalysis researchers are not familiar with Heusler alloys. We report their potential as novel catalysts focusing on the selective hydrogenation of alkynes. We found that Co2MnGe and Co2FeGe alloys have great alkene selectivity. Mutual substitution of Mn and Fe (Co2Mn x Fe1-x Ge) enhanced the reaction rate without changing selectivity. The substitution of Ga for Ge decreased the selectivity but increased the reaction rate monotonically with Ga composition. Elucidation of these mechanisms revealed that the fine tuning of catalytic properties is possible in Heusler alloys by separately using ligand and ensemble effects of elemental substitution.

Controlling Cooperative CO2 Adsorption in Diamine-Appended Mg2(dobpdc) Metal-Organic Frameworks

In the transition to a clean-energy future, CO2 separations will play a critical role in mitigating current greenhouse gas emissions and facilitating conversion to cleaner-burning and renewable fuels. New materials with high selectivities for CO2 adsorption, large CO2 removal capacities, and low regeneration energies are needed to achieve these separations efficiently at scale. Here, we present a detailed investigation of nine diamine-appended variants of the metal-organic framework Mg2(dobpdc) (dobpdc4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) that feature step-shaped CO2 adsorption isotherms resulting from cooperative and reversible insertion of CO2 into metal-amine bonds to form ammonium carbamate chains. Small modifications to the diamine structure are found to shift the threshold pressure for cooperative CO2 adsorption by over 4 orders of magnitude at a given temperature, and the observed trends are rationalized on the basis of crystal structures of the isostructural zinc frameworks obtained from in situ single-crystal X-ray diffraction experiments. The structure-activity relationships derived from these results can be leveraged to tailor adsorbents to the conditions of a given CO2 separation process. The unparalleled versatility of these materials, coupled with their high CO2 capacities and low projected energy costs, highlights their potential as next-generation adsorbents for a wide array of CO2 separations.

Entropy prediction for H2 adsorption in metal-organic frameworks

Entropy is an important thermodynamic property and serves as a bridge connecting equilibrium and non-equilibrium systems, which provides a basic understanding of various practical phenomena. In this study, classical density functional theory was introduced to efficiently predict entropy. The theory was applied to a high-throughput prediction of entropy and excess entropy for H2 adsorption in metal-organic frameworks. It seems that the entropy screening and uptake screening are generally equivalent at high temperature. Based on the entropy screening, the best hydrogen storage materials have been identified. The correlations between entropy and thermodynamic properties, such as uptake, isosteric heat and adsorption degree, were examined and are explained. The results imply that among the tested thermodynamic properties, the correlation between entropy and isosteric heat is the strongest.

Direct inference of site strength in basic solids upon CO2 adsorption: enthalpy-entropy compensation effects

The adsorption of CO2 coupled to calorimetry is a state-of-the-art technique for characterizing the basic properties of solids. In this paper, we show that the differential heat and entropy curves measured upon CO2 adsorption on a basic solid can be reasonably estimated from a single CO2 isotherm with no need for any independent heat (calorimetric) measurement. Our method relies on two important observations: (1) formulation of generalized F-H-TS thermodynamic isotherms, the former (F) being directly generated from the raw CO2 isotherms, and (2) the presence of unexpected enthalpy-entropy compensation effects upon CO2 adsorption linking the integral enthalpy and entropy of adsorption until saturation for different solids. Our thermodynamic method has been validated using a broad library of basic solids with variable site strength and heterogeneity. Finally, a new scale of basicity is proposed using the parameters fitted from the thermodynamic isotherm (free energy basis) as descriptors of basic strength. This method opens an avenue to the inference of site strength of basic solids without the need for expensive calorimeters.

Probing zeolites by vibrational spectroscopies

This review addresses the most relevant aspects of vibrational spectroscopies (IR, Raman and INS) applied to zeolites and zeotype materials. Surface Brønsted and Lewis acidity and surface basicity are treated in detail. The role of probe molecules and the relevance of tuning both the proton affinity and the steric hindrance of the probe to fully understand and map the complex site population present inside microporous materials are critically discussed. A detailed description of the methods needed to precisely determine the IR absorption coefficients is given, making IR a quantitative technique. The thermodynamic parameters of the adsorption process that can be extracted from a variable-temperature IR study are described. Finally, cutting-edge space- and time-resolved experiments are reviewed. All aspects are discussed by reporting relevant examples. When available, the theoretical literature related to the reviewed experimental results is reported to support the interpretation of the vibrational spectra on an atomic level.

H2 interaction with divalent cations in isostructural MOFs: a key study for variable temperature infrared spectroscopy

Systematic studies of H2 adsorption by variable temperature infrared (VTIR) spectroscopy have added value in the characterization of hydrogen storage materials. As a key study to describe the potential of the method, here we report VTIR spectroscopy results of H2 adsorption at isostructural MOFs CPO-27-M (M = Mg, Mn, Co, Ni, Zn). The strongest perturbation of H2 vibrational frequency is due to the interaction with an open metal site. Although ionic radius is an empirical value, the direct correlation between ionic radii of the metal cation and H2 interaction energy is found in MOFs of the same topology. The highest enthalpy of hydrogen adsorption 15 ± 1 kJ mol(-1) was found for Ni(2+). VTIR results of H2 adsorption at isostructural MOFs CPO-27-M (M = Mg, Mn, Co, Ni, Zn) were compared with data obtained from analogous studies performed on a large variety of microporous materials (MOFs and zeolites), underlining the relevance of the approach to get reliable energetic and entropic (ΔH(0) and ΔS(0)) values to be compared with computational data and isosteric heats.

Optimizing nanoporous materials for gas storage

In this work, we address the question of which thermodynamic factors determine the deliverable capacity of methane in nanoporous materials. The deliverable capacity is one of the key factors that determines the performance of a material for methane storage in automotive fuel tanks. To obtain insights into how the molecular characteristics of a material are related to the deliverable capacity, we developed several statistical thermodynamic models. The predictions of these models are compared with the classical thermodynamics approach of Bhatia and Myers [Bhatia and Myers, Langmuir, 2005, 22, 1688] and with the results of molecular simulations in which we screen the International Zeolite Association (IZA) structure database and a hypothetical zeolite database of over 100,000 structures. Both the simulations and our models do not support the rule of thumb that, for methane storage, one should aim for an optimal heat of adsorption of 18.8 kJ mol(-1). Instead, our models show that one can identify an optimal heat of adsorption, but that this optimal heat of adsorption depends on the structure of the material and can range from 8 to 23 kJ mol(-1). The different models we have developed are aimed to determine how this optimal heat of adsorption is related to the molecular structure of the material.

Hydrogen storage in metal-organic frameworks

New materials capable of storing hydrogen at high gravimetric and volumetric densities are required if hydrogen is to be widely employed as a clean alternative to hydrocarbon fuels in cars and other mobile applications. With exceptionally high surface areas and chemically-tunable structures, microporous metal-organic frameworks have recently emerged as some of the most promising candidate materials. In this critical review we provide an overview of the current status of hydrogen storage within such compounds. Particular emphasis is given to the relationships between structural features and the enthalpy of hydrogen adsorption, spectroscopic methods for probing framework-H(2) interactions, and strategies for improving storage capacity (188 references).