Interlayer Potential for Graphene/h-BN Heterostructures.

We present a new force-field potential that describes the interlayer interactions in heterojunctions based on graphene and hexagonal boron nitride (h-BN). The potential consists of a long-range attractive term and a short-range anisotropic repulsive term. Its parameters are calibrated against reference binding and sliding energy profiles for a set of finite dimer systems and the periodic graphene/h-BN bilayer, obtained from density functional theory using a screened-exchange hybrid functional augmented by a many-body dispersion treatment of long-range correlation.