Kolmogorov-Crespi Potential For Multilayer Transition-Metal Dichalcogenides: Capturing Structural Transformations in Moire Superlattices

We develop parameters for the interlayer Kolmogorov-Crespi (KC) potential to study structural features of four transition-metal dichalcogenides (TMDs): MoS2 WS2, MoSe2, and WSe2. We also propose a mixing rule to extend the parameters to their heterostructures. Moire super lattices of twisted bilayer TMDs have been recently shown to host shear solitons, topological point defects, and ultraflat bands close to the valence band edge. Performing structural relaxations at the density functional theory (DFT) level is a major bottleneck in the study of these systems. We show that the parametrized KC potential can be used to obtain atomic relaxations in good agreement with DFT relaxations. Furthermore, the moire superlattices relaxed using DFT and the proposed force field yield very similar electronic band structures.