Competing nematic, antiferromagnetic, and spin-flux orders in the ground state of bilayer graphene

We analyze the phase diagram of bilayer graphene (BLG) at zero temperature and zero doping. Assuming that at high energies the electronic system of BLG can be described within a weak-coupling theory (consistent with the experimental evidence), we systematically study the evolution of the couplings with going from high to low energies. The divergences of the couplings at some energies indicate the tendency towards certain symmetry breakings. Carrying out this program, we found that the phase diagram is determined by microscopic couplings defined on the short distances (initial conditions). We explored all plausible space of these initial conditions and found that the three states have the largest phase volume of the initial couplings: nematic, antiferromagnetic, and spin flux (a.k.a. quantum spin Hall). In addition, ferroelectric and two superconducting phases appear only near the very limits of the applicability of the weak-coupling approach. The paper also contains the derivation and analysis of the renormalization group equations and the group theory classification of all the possible phases which might arise from the symmetry breakings of the lattice, spin rotation, and gauge symmetries of graphene.