Graphene bilayers exhibit zero-energy flat bands at a discrete series of magic twist angles. In the absence intra-sublattice inter-layer hopping, states satisfy Dirac equation with non-abelian SU(2) gauge potential that cannot be diagonalized globally. We develop semiclassical WKB approximation scheme for this by introducing dimensionless Planck's constant proportional to angle, solving lineari...
