Stacking domains and dislocation networks in marginally twisted bilayers of transition metal dichalcogenides

We apply a multiscale modeling approach to study lattice reconstruction in marginally twisted bilayers of transition metal dichalcogenides (TMD). For this, we develop DFT-parametrized interpolation formulae for interlayer adhesion energies of MoSe2, WSe2, MoS2, and WS2, combine those with elasticity theory, and analyze the bilayer lattice relaxation into mesoscale domain structures. Paying particular attention to the inversion asymmetry of TMD monolayers, we show that 3R and 2H stacking domains, separated by a network of dislocations develop for twist angles $^{}