DOI: 10.5281/zenodo.7243734
title: Relaxation effects in twisted bilayer molybdenum disulfide: structure, stability, and electronic properties

description: <strong>Abstract</strong> Manipulating the interlayer twist angle is a powerful tool to tailor the properties of layered two-dimensional crystals. The twist angle has a determinant impact on these systems' atomistic structure and electronic properties. This includes the corrugation of individual layers, formation of stacking domains and other structural elements, and electronic structure changes due to the atomic reconstruction and superlattice effects. However, how these properties change with the twist angle, <em>θ</em>, is not yet well understood. Here, we monitor the change of twisted bilayer (tBL) MoS<sub>2</sub> characteristics as a function of <em>θ</em>. We identify distinct structural regimes, each with particular structural and electronic properties. We employ a hierarchical approach ranging from a reactive force field through the density-functional-based tight-binding approach and density-functional theory. To obtain a comprehensive overview, we analyzed a large number of tBLs with twist angles in the range of \(\theta=0.2^\circ\dots59.6^\circ\). Some systems include up to half a million atoms, making structure optimization and electronic property calculation challenging. For \(13^\circ \lessapprox \theta \lessapprox 47^\circ\), the structure is well-described by a moiré regime composed of two rigidly twisted monolayers. At small twist angles (\(\theta\leq3^\circ\) and \(57^\circ\leq\theta\)), a domain-soliton regime evolves, where the structure contains large triangular stacking domains, separated by a network of strain solitons and short-ranged high-energy nodes. The corrugation of the layers and the emerging superlattice of solitons and stacking domains affects the electronic structure. Emerging predominant characteristic features are Dirac cones at <em>K</em> and kagome bands. These features flatten for <em>θ</em> approaching 0<sup>∘</sup> and 60<sup>∘</sup>. Our results show at which range of <em>θ</em> the characteristic features of the reconstruction, namely extended stacking domains, the soliton network, and superlattice, emerge and give rise to exciting electronics. We expect our findings also to be relevant for other tBL systems. DOI: 10.1088/2053-1583/aceb75 <strong>Overview</strong> This repository contains calculation files, optimized structures, and visualization movies for studies of twisted-bilayer MoS<sub>2</sub>, focussing on structural properties and electronic structure. Each directory has its own README.md file with additional information, separated by what data is included and the method used. <strong>Geometry optimization</strong> Directory `calc_structure_optimization_ReaxFF`: calculation files of the structure optimization of all studied structures, done with ReaxFF. Directory `calc_structure_optimization_DFT`: validation calculation files of the ReaxFF-optimized structures using DFT optimization. <strong>Electronic structure calculations</strong> Directory `calc_electronic_properties_DFT`: calculation files of electronic structure calculations on the DFT level. Directory `calc_electronic_properties_DFTB`: calculation files of electronic structure calculations on the DFTB level <strong>Results</strong> Directory `structures_rigidly_twisted`: structure files in cif format of the rigidly twisted (flat) systems, labeled by their twist angle. Directory `structures_fully_optimized`: structure files in cif format of the fully ReaxFF-optimized systems, labeled by their twist angle. Directory `movies`: visualization of the change of the interlayer distance landscape and the strain fields with the twist angle. Additionally, the script `plot_interlayer_distance.py` is included, which was used to create the individual frames of the movie showing the interlayer distance.