▎ 摘　　要

In the search for efficient molecular dynamics simulation models both simplicity and acceptable accuracy matter. In the present study, a model of the graphene-H-2 physisorption system is used to explore its performance and limitations under canonical NVT and microcanonical NVE simulation conditions. The model implies several simplifications that can be summarized in (a) a single ideal planar frozen graphene-like layer of C atoms, (b) rigid rotor H-2 molecules and (c) interaction potentials written as C-H-2 and H-2-H-2 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 simulations lead to the formation of a single layer of adsorbed H-2 molecules in dynamically stable equilibrium with a fluid-phase region. In addition, the incipient formation of secondary 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.