Abstract
Graphene bilayers exhibit zero-energy flatbands at a discrete series of magic twist angles. In the absence of intrasublattice interlayer hopping, zero-energy states satisfy a Dirac equation with a non-Abelian SU(2) gauge potential that cannot be diagonalized globally. We develop a semiclassical WKB approximation scheme for this Dirac equation by introducing a dimensionless Planck's constant proportional to the twist angle, solving the linearized Dirac equation around AB and BA turning points, and connecting Airy function solutions via bulk WKB wave functions. We find zero-energy solutions at a discrete set of values of the dimensionless Planck's constant, which we obtain analytically. Our analytic flatband twist angles correspond closely to those determined numerically in previous work.
Sign-in/Register to access full text options 
Free) 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call
