Instabilities of interacting electrons on the honeycomb bilayer

We investigate the instabilities of interacting electrons on the honeycomb bilayer by means of the functional renormalization group. Using model parameters as determined by ab initio calculations for graphene and graphite for interactions up to the third-nearest neighbor puts the system close to the boundary between antiferromagnetic and quantum spin Hall instabilities. Importantly, the energy scales for these instabilities are large. Thus imperfections and deviations from the basic model are expected to play a major role in real bilayer graphene, where interaction effects seem to be seen only at smaller scales. We therefore analyze how reducing the critical scale and small doping of the layers affect the instabilities. For smaller on-site repulsions, we also find an instability toward a gapless charge-density wave state.