Therefore, for very small twist angles, where the moiré cells contain far more than 5,000 atoms, we parameterize a TB model, fit the hopping and on-site parameters as function of the moiré cell's lattice constant, and then extrapolate the TB results to very small twist angles with respect to the two high-symmetry stackings 𝑅 ℎ 𝑀 and 𝐻 ℎ ℎ .