Kolmogorov-Crespi Potential For Multilayer Transition Metal Dichalcogenides: Capturing Structural Transformations In Moir\'e Superlattices

We develop parameters for the interlayer Kolmogorov-Crespi (KC) potential to study structural features of four transition metal dichalcogenides (TMDs): MoS$_2$, WS$_2$, MoSe$_2$ and WSe$_2$. We also propose a mixing rule to extend the parameters to their heterostructures. Moir\'e superlattices of twisted bilayer TMDs have been recently shown to host shear solitons, topological point defects and ultraflatbands close to the valence band edge. Performing structural relaxations at the 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 moir\'e superlattices relaxed using DFT and the proposed forcefield yield very similar electronic band structures.