Landau level spectrum for bilayer graphene in a tilted magnetic field

We consider a graphene bilayer in a constant magnetic field of arbitrary orientation, i.e., tilted with respect to the graphene plane. In the low energy approximation to the tight-binding model with Peierls substitution, we find the Landau level spectrum analytically in terms of spheroidal functions and the respective eigenvalues. We compare our result to the perpendicular and purely in-plane field cases. In the limit of perpendicular field we reproduce the known equidistant spectrum for Landau levels. In the opposite limit of large in-plane field this spectrum becomes two-fold degenerate, which is a consequence of Dirac point splitting induced by the in-plane field.