The top five occupied bands -two Dirac and three kgm bands -were fitted using the graphene-and hBN-like model Hamiltonians as shown in Figure ). The applicability of these Hamiltonians to graphene and hBN is demonstrated in Figure and Table . We set up using PythTB, the on-site elements and hopping terms were fitted to the DFTB band structure using the basin-hopping method as implemented in SciPy with the rootmean-square deviation (RMSD) between TB and DFTB bands as property to be minimized.