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.

Prediction of Multiple Hydrogen Ligation at a Vanadium(II) Site in a Metal-Organic Framework

Densifying hydrogen in a metal-organic framework (MOF) at moderate pressures can circumvent challenges associated with high-pressure compression. The highly tunable structural and chemical composition in MOFs affords vast possibilities to optimize binding interactions. At the heart of this search are the nanoscale characteristics of molecular adsorption at the binding site(s). Using density functional theory (DFT) to model binding interactions of hydrogen to the exposed metal site of cation-exchanged MFU-4l, we predict multiple hydrogen ligation of H₂ at the first coordination sphere of V(II) in V(II)-exchanged MFU-4l. We find that the strength of this binding between the metal site and H₂ molecules can be tuned by altering the halide counterion adjacent to the metal site and that the fluoride containing node affords the most favorable interactions for high-density H₂ storage. Using energy decomposition analysis, we delineate electronic contributions that enable multiple hydrogen ligation and demonstrate its benefits for hydrogen adsorption and release at modest pressures.

A Bird's-Eye View on Polymer-Based Hydrogen Carriers for Mobile Applications

Globally, reducing CO2 emissions is an urgent priority. The hydrogen economy is a system that offers long-term solutions for a secure energy future and the CO2 crisis. From hydrogen production to consumption, storing systems are the foundation of a viable hydrogen economy. Each step has been the topic of intense research for decades; however, the development of a viable, safe, and efficient strategy for the storage of hydrogen remains the most challenging one. Storing hydrogen in polymer-based carriers can realize a more compact and much safer approach that does not require high pressure and cryogenic temperature, with the potential to reach the targets determined by the United States Department of Energy. This review highlights an outline of the major polymeric material groups that are capable of storing and releasing hydrogen reversibly. According to the hydrogen storage results, there is no optimal hydrogen storage system for all stationary and automotive applications so far. Additionally, a comparison is made between different polymeric carriers and relevant solid-state hydrogen carriers to better understand the amount of hydrogen that can be stored and released realistically.

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.

Active metal single-sites based on metal–organic frameworks: construction and chemical prospects

Metal single-point is a novel and potential design strategy that has been applied for the development of metal organic frameworks.

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.

The effect of co-adsorbed solvent molecules on H2 binding to metal alkoxides

The introduction of metal alkoxides has been proposed as an attractive option to enhance hydrogen binding energies in porous materials such as metal-organic frameworks (MOFs) for room-temperature hydrogen storage applications. The presence of residual solvent molecules from MOF synthesis can, however, affect the performance of these functional groups. We perform quantum chemical calculations to predict solvent binding energies onto the metal-alkoxides and the temperatures required to drive off the solvent molecules and successfully activate porous materials with these moieties. Calculations are performed for Li, Mg, Zn, Cu, and Ni alkoxides and chloroform (CHCl3), dimethylformamide (DMF), ethanol, methanol, and water solvent molecules. We identify CHCl3 as a promising solvent that can be removed from these alkoxides at mild temperatures, whereas DMF binds strongly to the metal alkoxides and removal would require temperatures above the present upper bound of thermal stability in MOFs. As a second objective, we calculated the binding energies of hydrogen to metal alkoxide-solvent complexes to explore the effect of any solvent molecules that cannot be removed.