In this section, we consider the influence of interlayer disorder on the interacting mean-field state in RMG, i.e. we preserve translational invariance within each graphene layer, but take into account random tight-binding parameters in the perpendicular-to-layer direction. This has a close analogy with studies of the SSH model [1]}, [2]} wherein the effects on the zero-energy edge states of chiral-symmetry-preserving or breaking disorder are considered [3]}, [4]}, [5]}, [6]}, [7]}. A major difference here is that we consider the influence of disorder on the interacting mean-field state in which the exchange potential has already broken chiral symmetry and gapped the spectrum.
