Stacking faults have been considered previously in graphene multilayer systems [1]}, [2]} and in RMG in particular [3]}, [4]}, [5]}. Single stacking faults, e.g. a Bernal fault or a twin boundary fault, Fig. REF , within a large RMG system are interesting because they introduce additional flat bands into the energy spectrum. Here we focus on the properties of the interacting ground state at half filling.
