Abstract
We use the recently developed van der Waals density functional (vdW-DF) along with the C09x exchange functional to explore the interactions of molecular hydrogen with a single sheet of graphene. Our calculations demonstrate the importance of the inclusion of dispersion interactions for evaluating the adsorption capacity of graphitic, sparsely packed materials. In particular we show that, although the commonly used local density approximation (LDA) exchange-correlation functional gives reasonable H$_2$-graphene separation distances, these interactions are in fact too strong and decrease too rapidly as the H$_2$-graphene distance increases. On the other hand, calculations employing the generalized gradient approximation (GGA) for exchange and correlation exhibit very little binding. The vdW-DF$^{\rm C09x}$ functional, however, gives binding very similar to benchmark second-order M$\o$ller-Plesset (MP2) theory results for the H$_2$-coronene interaction and has an appropriate long range interaction with the graphene sheet. This shortcoming of the LDA and GGA functionals is shown to have significant consequences on the overall adsorption densities of H$_2$ near the graphene, thus emphasizing the fundamental importance of properly characterizing the interactions of adsorbed molecules in porous media.