We apply a multiscale modeling approach to study lattice reconstruction in marginally twisted bilayers of transition metal dichalcogenides (TMD). For this, we develop density functional theory parametrized interpolation ...