When one stacks a couple of graphene layers on top of each other with a twist by an angle \(\theta \) , one needs to take into account interlayer couplings, which extends modifies the above monoloayer Hamiltonian (REF ). The low energy continuum model capturing these features and robustly accounting for dispersionless states is the so-called `chiral' one [1]}:
\(\tilde{H} = v_F \int d^2 \mathbf {r} \, \Phi ^{\dagger } (\mathbf {r},t) D \Phi (\mathbf {r},t).\) 
