Anisotropic Interlayer Force Field for Group-VI Transition Metal Dichalcogenides

Abstract

An anisotropic interlayer force field that describes the interlayer interactions in homogeneous and heterogeneous interfaces of group-VI transition metal dichalcogenides (MX2 where M = Mo, W and X = S, Se) is presented. The force field is benchmarked against density functional theory calculations for bilayer systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmented by a nonlocal many-body dispersion treatment of long-range correlation. The parametrization yields good agreement with reference calculations of binding energy curves and sliding potential energy surfaces. It is found to be transferable to TMD junctions outside the training set that contain the same atom types. Calculated bulk moduli agree with most previous dispersion corrected DFT predictions, which underestimate available experimental values. Calculated phonon spectra of the various junctions under consideration demonstrate the importance of appropriately treating the anisotropic nature of layered interfaces. Considering our previous parameterization for MoS2, the interlayer potential enables accurate and efficient large-scale simulations of the dynamical, tribological, and thermal transport properties of a large set of homogeneous and heterogeneous TMD interfaces.

References in corpus (5)

• Observation of Moiré Excitons in WSe2/WS2 Heterostructure Superlattices
• Pressure induced metallization with absence of structural transition in layered MoSe2
• Robust Superlubricity in Graphene/h-BN Heterojunctions
• Negative Friction Coefficients in Superlubric Graphite-Hexagonal Boron Nitride Heterojunctions
• Mechanical and Tribological Properties of Layered Materials Under High Pressure: Assessing the Importance of Many-Body Dispersion Effects

Cited by in corpus (5)

• Moiré-Driven Interfacial Thermal Transport in Twisted Transition Metal Dichalcogenides
• PYSED: A tool for extracting kinetic-energy-weighted phonon dispersion and lifetime from molecular dynamics simulations
• Universal Moiré Buckling of Freestanding 2D Bilayers
• Modular hybrid machine learning and physics-based potentials for scalable modeling of van der Waals heterostructures
• Two-dimensional moiré phonon polaritons