Grand canonical Monte Carlo (GCMC) simulations were performed to investigate hydrogen sorption in an rht-type metal−organic framework (MOF), PCN-61. The MOF was shown to have a large hydrogen uptake, and this was studied using three different hydrogen potentials, effective for bulk hydrogen, but of varying sophistication: a model that includes only repulsion/dispersion parameters, one augmented with charge-quadrupole interactions, and one supplemented with many-body polarization interactions. Calculated hydrogen uptake isotherms and isosteric heats of adsorption, Q st , were in quantitative agreement with experiment only for the model with explicit polarization. This success in reproducing empirical measurements suggests that modeling MOFs that have open metal sites is feasible, though it is often not considered to be well described via a classical potential function; here it is shown that such systems may be accurately described by explicitly including polarization effects in an otherwise traditional empirical potential. Decomposition of energy terms for the models revealed deviations between the electrostatic and polarizable results that are unexpected due to just the augmentation of the potential surface by the addition of induction. Charge-quadrupole and induction energetics were shown to have a synergistic interaction, with inclusion of the latter resulting in a significant increase in the former. Induction interactions strongly influence the structure of the sorbed hydrogen compared to the models lacking polarizability; sorbed hydrogen is a dipolar dense fluid in the MOF. This study demonstrates that many-body polarization makes a critical contribution to gas sorption structure and must be accounted for in modeling MOFs with polar interaction sites.

Monte Carlo simulations were performed modeling hydrogen sorption in a recently synthesized metal-organic framework material (MOF) that exhibits large molecular hydrogen uptake capacity. The MOF is remarkable because at 78 K and 1.0 atm it sorbs hydrogen at a density near that of liquid hydrogen (at 20 K and 1.0 atm) when considering H2 density in the pores. Unlike most other MOFs that have been investigated for hydrogen storage, it has a highly ionic framework and many relatively small channels. The simulations demonstrate that it is both of these physical characteristics that lead to relatively strong hydrogen interactions in the MOF and ultimately large hydrogen uptake. Microscopically, hydrogen interacts with the MOF via three principle attractive potential energy contributions: Van der Waals, charge-quadrupole, and induction. Previous simulations of hydrogen storage in MOFs and other materials have not focused on the role of polarization effects, but they are demonstrated here to be the dominant contribution to hydrogen physisorption. Indeed, polarization interactions in the MOF lead to two distinct populations of dipolar hydrogen that are identified from the simulations that should be experimentally discernible using, for example, Raman spectroscopy. Since polarization interactions are significantly enhanced by the presence of a charged framework with narrow pores, MOFs are excellent hydrogen storage candidates.

An anisotropic many-body H2 potential energy function has been developed for use in heterogeneous systems. The intermolecular potential has been derived from first principles and expressed in a form that is readily transferred to exogenous systems, e.g. in modeling H2 sorption in solid-state materials. Explicit many-body polarization effects, known to be important in simulating hydrogen at high density, are incorporated. The analytic form of the potential energy function is suitable for methods of statistical physics, such as Monte Carlo or Molecular Dynamics simulation. The model has been validated on dense supercritical hydrogen and demonstrated to reproduce the experimental data with high accuracy.

Grand canonical Monte Carlo simulations of H2 sorption were performed in the metal–organic framework rht-MOF-1. The binding sites were revealed by combining simulation and inelastic neutron scattering data.