Twisted bilayer graphene I. Matrix elements, approximations, perturbation theory and a k· p 2-Band model
Abstract
We investigate the Twisted Bilayer Graphene (TBG) model to obtain an analytic understanding of its energetics and wavefunctions needed for many-body calculations. We provide an approximation scheme which first elucidates why the BM KM-point centered calculation containing only 4 plane-waves provides a good analytical value for the first magic angle. The approximation scheme also elucidates why most many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a M-point centered model of 6 plane-waves, we analytically understand the small M-point gap between the active and passive bands in the isotropic limit w0=w1. Furthermore, we analytically calculate the group velocities of passive bands in the isotropic limit, and show that they are almost doubly degenerate, while no symmetry forces them to be. Furthermore, away from M and KM points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite it is defined only by vanishing KM-point Dirac velocity. We derive analytically a connected "magic manifold" w1=21+w02-2+3w02, on which the bands remain extremely flat as w0 is tuned between the isotropic (w0=w1) and chiral (w0=0) limits. We analytically show why going away from the isotropic limit by making w0 less (but not larger) than w1 increases the M- point gap between active and passive bands. Finally, perturbatively, we provide an analytic M point k· p 2-band model that reproduces the TBG band structure and eigenstates in a certain w0,w1 parameter range. Further refinement of this model suggests a possible faithful 2-band M point k· p model in the full w0, w1 parameter range.