Modelization of the H 2 adsorption on graphene and molecular dynamics simulation

In the search for ef cient molecular dynamics simulation models both simplicity and acceptable accuracy matter. In the present study, a model of the graphene-H2 physisorption system is used to explore its performance and limitations under canonical NVT and microcanonical NVE simulation conditions. The model implies several simpli - cations that can be summarized in (a) a single ideal planar frozen graphene-like layer of C atoms, (b) rigid rotor H2 molecules and (c) interaction potentials written as C–H2 and H2–H2 site–site Improved Lennard-Jones potentials parameterized to reproduce DFT calculations. This model can be used in a variety of molecular dynamics simulation conditions, both in NVT and NVE ensembles. Such simula- tions lead to the formation of a single layer of adsorbed H2 molecules in dynamically stable equilibrium with a uid- phase region. In addition, the incipient formation of sec- ondary layers for high-density conditions is also observed. Some properties as average pressure, temperatures and fluid-phase densities are discussed as well as possible improvements of the model.