We have generalized the mean-field Hartree Fock description [1]}, [2]}, [3]}, [4]}, [5]} to provide a comprehensive qualitative description of broken symmetry ground states in RMG, including the effects of defects including random disorder and stacking faults.
The non-trivial topology of the low-energy bands is reflected in large Berry curvature and Chern numbers per spin-valley flavor.
An obvious generalization is to a number of stacking faults separating rhombohedral sections with different numbers of layers \(m\) , each section contributing a pair of low-energy flat bands and a peak in the Berry curvature at a characteristic \(m\) -dependent wave vector.
In RMG, each stacking fault contributes a pair of low-energy flat bands because they are more complicated than the domain walls usually considered in the SSH model which consist of isolated monomers or trimers [6]}. The Bernal fault corresponds to two monomers (and a 4-mer), the twin boundary fault is a monomer plus a trimer.
