Strain-Induced Landau Levels in Arbitrary Dimensions with an Exact Spectrum

Certain nonuniform strain applied to graphene flakes has been shown to induce pseudo-Landau levels in the single-particle spectrum, which can be rationalized in terms of a pseudomagnetic field for electrons near the Dirac points. However, this Landau level structure is, in general, approximate and restricted to low energies. Here, we introduce a family of strained bipartite tight-binding models in arbitrary spatial dimension d and analytically prove that their entire spectrum consists of perfectly degenerate pseudo-Landau levels. This construction generalizes the case of triaxial strain on graphene's honeycomb lattice to arbitrary d; in d=3, our model corresponds to tetraxial strain on the diamond lattice. We discuss general aspects of pseudo-Landau levels in arbitrary d.