In a basis of \(p_z\)  orbitals on \(A_1\) , \(B_1\) , \(A_2\) , \(B_2\) , ..., \(A_N\) , \(B_N\)  sites, the non-interacting Hamiltonian of RMG with \(N\)  layers  [1]}, [2]}, [3]} may be written near each valley as
\(H =\begin{pmatrix}D & V & 0 & 0 & \cdots \\V^{\dagger } & D & V & 0 & \cdots \\0 & V^\dagger & D & V & \cdots \\0 & 0 & V^\dagger & D & \cdots \\\vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix},\) 
