Electron-electron interactions are included within a mean-field Hartree-Fock approximation [1]}, [2]}, [3]} and, in particular, we follow the methodology applied to Bernal-stacked multilayer graphene in Refs. [4]}, [5]}. The total Hamiltonian is \({\hat{H}}_{\mathrm {tot}} = {\hat{H}} + {\hat{V}}_{\mathrm {MF}}\) 
where
\({\hat{H}} &=& \sum _{\mathbf {k} \sigma X \!X^{\prime }} H_{\mathbf {k} X\!X^{\prime }} c_{\mathbf {k} \sigma X}^\dagger c_{\mathbf {k} \sigma X^{\prime }} ,  \\{\hat{V}}_{\mathrm {MF}} &=& \sum _{\mathbf {k} \sigma X \!X^{\prime }} \left[ U_{X}^{(\mathrm {H})} \delta _{X\!X^{\prime }} + W_{\mathbf {k} \sigma X \!X^{\prime }} \right] c_{\mathbf {k} \sigma X}^\dagger c_{\mathbf {k} \sigma X^{\prime }}  .\) 
