High H2 Storage of Hexagonal Metal−Organic Frameworks from First-Principles-Based Grand Canonical Monte Carlo Simulations

Stimulated by the recent report by Yaghi and co-workers of hexagonal metal−organic frameworks (MOF) exhibiting reversible binding of up to 7.5 wt % at 77 K and 70 bar for MOF-177 (called here IRMOF-2-24), we have predicted additional trigonal organic linkers, including IRMOF-2-60, which we calculate to bind 9.7 wt % H2 storage at 77 K and 70 bar, the highest known value for 77 K. These calculations are based on grand canonical Monte Carlo (GCMC) simulations using force fields that match accurate quantum mechanical calculations on the binding of H2 to prototypical systems. These calculations were validated by comparison to the experimental loading curve for IRMOF-2-24 at 77K. We then used the theory to predict the effect of doping Li into the hexagonal MOFs, which leads to substantial H2 density even at ambient temperatures. For example, IRMOF-2-96-Li leads to 6.0 wt % H2 storage at 273 K and 100 bar, the first material to attain the 2010 DOE target

Improved H_2 Storage in Zeolitic Imidazolate Frameworks Using Li^+, Na^+, and K^+ Dopants, with an Emphasis on Delivery H_2 Uptake

We use grand canonical Monte Carlo simulations with first principles based force fields to show that alkali metal (Li^+, Na^+, and K^+)-doped zeolitic imidazolate frameworks (ZIFs) lead to significant improvement of H_2 uptake at room temperature. For example, at 298 K and 100 bar, Li-ZIF-70 totally binds to 3.08 wt % H_2, Na-ZIF-70 to 2.19 wt % H_2, and K-ZIF-70 to 1.62 wt % H_2, much higher than 0.74 wt % H_2 for pristine ZIF-70. Thus, the dopant effect follows the order of Li-ZIF > Na-ZIF > K-ZIF, which correlates with the H_2 binding energies to the dopants. Moreover, the total H_2 uptake is higher at lower temperatures: 243 K > 273 K > 298 K. On the other hand, delivery H_2 uptake, which is the difference between the total adsorption at the charging pressure (say 100 bar) and the discharging pressure (say 5 bar), is the important factor for practical on-board hydrogen storage in vehicles. We show that delivery H_2 uptake leads to Na-ZIF-70 (1.37 wt %) > K-ZIF-70 (1.25 wt %) > Li-ZIF-70 (1.07 wt %) > ZIF-70 (0.68 wt %), which is different from the trend from the total and excess uptake. Moreover, the delivery uptake increases with increasing temperatures (i.e., 298 K > 273 K > 243 K)! To achieve high delivery H_2 uptake at room temperature, the large free volume of ZIFs is required. We find that higher H_2 binding energy needs not always lead to higher delivery H_2 uptake

Zeolitic Imidazolate Frameworks as H_2 Adsorbents: Ab Initio Based Grand Canonical Monte Carlo Simulation

We report the H_2 uptake behavior of 10 zeolitic−imidazolate frameworks (ZIFs), based on grand canonical Monte Carlo (GCMC) simulations. The force fields (FFs) describing the interactions between H_2 and ZIF in the GCMC were based on ab initio quantum mechanical (QM) calculations (MP2) aimed at correctly describing London dispersion (van der Waals attraction). Thus these predictions of H_2 uptake are based on first principles (non empirical) and hence applicable to new framework materials for which there is no empirical data. For each of these 10 ZIFs we report the total and excess H_2 adsorption isotherms up to 100 bar at both 77 and 300 K. We report the hydrogen adsorption sites in the ZIFs and the relationships between H_2 uptake amount, isosteric heat of adsorption (Q_(st)), surface area, and free volume. Our simulation shows that various ZIFs lead to a variety of H_2 adsorption behaviors in contrast to the metal−organic frameworks (MOFs). This is because ZIFs leads to greater diversity in the adsorption sites (depending on both organic linkers and zeolite topologies) than in MOFs. In particular, the ZIFs uptake larger amounts of H_2 at low pressure because of the high H_2 adsorption energy, and ZIFs have a variety of H_2 adsorption sites. For example, ZIF-11 has an initial Q_(st) value of ~15 kJ/mol, which is higher than observed for MOFs. Moreover, the preferential H_2 adsorption site in ZIFs is onto the organic linker, not nearby the metallic joint as is the case for MOFs

Recent advances on simulation and theory of hydrogen storage in metal–organic frameworks and covalent organic frameworks

This critical review covers the application of computer simulations, including quantum calculations (ab initio and DFT), grand canonical Monte-Carlo simulations, and molecular dynamics simulations, to the burgeoning area of the hydrogen storage by metal–organic frameworks and covalent-organic frameworks. This review begins with an overview of the theoretical methods obtained from previous studies. Then strategies for the improvement of hydrogen storage in the porous materials are discussed in detail. The strategies include appropriate pore size, impregnation, catenation, open metal sites in metal oxide parts and within organic linker parts, doping of alkali elements onto organic linkers, substitution of metal oxide with lighter metals, functionalized organic linkers, and hydrogen spillover (186 references)

High H_2 Uptake in Li-, Na-, K-Metalated Covalent Organic Frameworks and Metal Organic Frameworks at 298 K

The Yaghi laboratory has developed porous covalent organic frameworks (COFs), COF102, COF103, and COF202, and metal–organic frameworks (MOFs), MOF177, MOF180, MOF200, MOF205, and MOF210, with ultrahigh porosity and outstanding H2 storage properties at 77 K. Using grand canonical Monte Carlo (GCMC) simulations with our recently developed first principles based force field (FF) from accurate quantum mechanics (QM), we calculated the molecular hydrogen (H2) uptake at 298 K for these systems, including the uptake for Li-, Na-, and K-metalated systems. We report the total, delivery and excess amount in gravimetric and volumetric units for all these compounds. For the gravimetric delivery amount from 1 to 100 bar, we find that eleven of these compounds reach the 2010 DOE target of 4.5 wt % at 298 K. The best of these compounds are MOF200-Li (6.34) and MOF200-Na (5.94), both reaching the 2015 DOE target of 5.5 wt % at 298 K. Among the undoped systems, we find that MOF200 gives a delivery amount as high as 3.24 wt % while MOF210 gives 2.90 wt % both from 1 to 100 bar and 298 K. However, none of these compounds reach the volumetric 2010 DOE target of 28 g H_2/L. The best volumetric performance is for COF102-Na (24.9), COF102-Li (23.8), COF103-Na (22.8), and COF103-Li (21.7), all using delivery g H_2/L units for 1–100 bar. These are the highest volumetric molecular hydrogen uptakes for a porous material under these thermodynamic conditions. Thus, one can obtain outstanding H_2 uptakes with Li, Na, and K doping of simple frameworks constructed from simple, cheap organic linkers. We present suggestions for strategies for synthesis of alkali metal-doped MOFs or COFs

Lithium-Doped Metal-Organic Frameworks for Reversible H_2 Storage at Ambient Temperature

To maximize reversible H_2 storage near room temperature and modest pressures, we propose Li doping of the metal-organic framework (MOFs) structures developed by the Yaghi group at UCLA (constructed using octahedral Zn−O−C clusters with aromatic carbon ring linkers). We tested this design using grand canonical Monte Carlo simulations with first-principles derived force fields and predict that at −30 °C and 100 bar the Li−MOF-C30 leads gravimetric H_2 uptake of 6.0 wt %, reaching the 2010 Department of Energy target (6.0 wt % in the temperature ranges of −30 to 80 °C and pressures ≤100 bar). Thus this promising material for practical hydrogen storage is worthy of developing experimental procedures for synthesis and characterization

Adsorption Mechanism and Uptake of Methane in Covalent Organic Frameworks: Theory and Experiment

We determined the methane (CH_4) uptake (at 298 K and 1 to 100 bar pressure) for a variety of covalent organic frameworks (COFs), including both two-dimensional (COF-1, COF-5, COF-6, COF-8, and COF-10) and three-dimensional (COF-102, COF-103, COF-105, and COF-108) systems. For all COFs, the CH_4 uptake was predicted from grand canonical Monte Carlo (GCMC) simulations based on force fields (FF) developed to fit accurate quantum mechanics (QM) [second order Møller−Plesset (MP2) perturbation theory using doubly polarized quadruple-ζ (QZVPP) basis sets]. This FF was validated by comparison with the equation of state for CH_4 and by comparison with the experimental uptake isotherms at 298 K (reported here for COF-5 and COF-8), which agrees well (within 2% for 1−100 bar) with the GCMC simulations. From our simulations we have been able to observe, for the first time, multilayer formation coexisting with a pore filling mechanism. The best COF in terms of total volume of CH_4 per unit volume COF absorbent is COF-1, which can store 195 v/v at 298 K and 30 bar, exceeding the U.S. Department of Energy target for CH_4 storage of 180 v/v at 298 K and 35 bar. The best COFs on a delivery amount basis (volume adsorbed from 5 to 100 bar) are COF-102 and COF-103 with values of 230 and 234 v(STP: 298 K, 1.01 bar)/v, respectively, making these promising materials for practical methane storage

Liquefaction of H2 molecules upon exterior surfaces of carbon nanotube bundles

We have used molecular dynamics simulations to investigate interaction of H2 molecules on the exterior surfaces of carbon nanotubes (CNTs): single and bundle types. At 80 K and 10 MPa, it is found that charge transfer occurs from a low curvature region to a high curvature region of the deformed CNT bundle, which develops charge polarization only on the deformed structure. The long-range electrostatic interactions of polarized charges on the deformed CNT bundle with hydrogen molecules are observed to induce a high local-ordering of H2 gas that results in hydrogen liquefaction. Our predicted heat of hydrogen liquefaction on the CNT bundle is 97.6 kcal kg^-1. On the other hand, hydrogen liquefaction is not observed in the CNT of a single type. This is because charge polarization is not developed on the single CNT as it is symmetrically deformed under the same pressure. Consequently, the hydrogen storage capacity on the CNT bundle is much higher due to liquefaction than that on the single CNT. Additionally, our results indicate that it would also be possible to liquefy H2 gas on a more strongly polarized CNT bundle at temperatures higher than 80 K

The theoretical study on interaction of hydrogen with single-walled boron nitride nanotubes. I. The reactive force field ReaxFFHBN development

We present a new reactive force field ReaxFFHBN derived to accurately model large molecular and condensed phase systems of H, B, and N atoms. ReaxFFHBN has been tested against quantum calculation data for B–H, B–B, and B–N bond dissociations and for H–B–H, B–N–B, and N–B–N bond angle strain energies of various molecular clusters. The accuracy of the developed ReaxFFHBN for B–N–H systems is also tested for (i) H–B and H–B bond energies as a function of out of plane in H–B(NH2)3 and H–N(BH2)3, respectively, (ii) the reaction energy for the B3N3H6+H2-->B3N3H8, and (iii) crystal properties such as lattice parameters and equations of states for the hexagonal type (h-BN) with a graphite structure and for the cubic type (c-BN) with a zinc-blende structure. For all these systems, ReaxFFHBN gives reliable results consistent with those from quantum calculations as it describes well bond breaking and formation in chemical processes and physical properties. Consequently, the molecular-dynamics simulation based on ReaxFFHBN is expected to give a good description of large systems (>2000 atoms even on the one-CPU machine) with hydrogen, boron, and nitrogen atoms

Theoretical study on interaction of hydrogen with single-walled boron nitride nanotubes. II. Collision, storage, and adsorption

Collision and adsorption of hydrogen with high incident kinetic energies on a single-walled boron nitride (BN) nanotube have been investigated. Molecular-dynamics (MD) simulations indicate that at incident energies below 14 eV hydrogen bounces off the BN nanotube wall. On the other hand, at incident energies between 14 and 22 eV each hydrogen molecule is dissociated at the exterior wall to form two hydrogen atoms, but only one of them goes through the wall. However, at the incident energies between 23 and 26 eV all of the hydrogen atoms dissociated at the exterior wall are found to be capable of going inside the nanotube and then to recombine to form hydrogen molecules inside the nanotube. Consequently, it is determined that hydrogen should have the incident energy >22 eV to go inside the nanotube. On the other hand, we find that the collisions using the incident energies >26 eV could result in damaging the nanotube structures. In addition our MD simulations find that hydrogen atoms dissociated at the wall cannot bind to either boron or nitrogen atoms in the interior wall of the nanotube