Searching for DFT-based methods that include dispersion interactions to calculate the physisorption of H2 on benzene and graphene

Simulations of the hydrogen storage capacities of nanoporous carbons require an accurate treatment of the interaction of the hydrogen molecule with the graphite-like surfaces of the carbon pores, which is dominated by the dispersion forces. These interactions are described accurately by high level quantum chemistry methods, like the Coupled Cluster method with single and double excitations and a non-iterative correction for triple excitations (CCSD(T)), but those methods are computationally very expensive for large systems and for massive simulations. Density functional theory (DFT)-based methods that include dispersion interactions at different levels of complexity are less accurate, but computationally less expensive. In order to find DFT-methods that include dispersion interactions to calculate the physisorption of H2 on benzene and graphene, with a reasonable compromise between accuracy and computational cost, CCSD(T), Moller-Plesset second-order perturbation theory method, and several DFT-methods hav...