Gross-Neveu-Heisenberg criticality from competing nematic and antiferromagnetic orders in bilayer graphene

Abstract

We study the phase diagram of an effective model of competing nematic and antiferromagnetic orders of interacting electrons on the Bernal-stacked honeycomb bilayer, as relevant for bilayer graphene. In the noninteracting limit, the model features a semimetallic ground state with quadratic band touching points at the Fermi level. Taking the effects of short-range interactions into account, we demonstrate the existence of an extended region in the mean-field phase diagram characterized by coexisting nematic and antiferromagnetic orders. By means of a renormalization group approach, we reveal that the quantum phase transition from nematic to coexistent nematic-antiferromagnetic orders is continuous and characterized by emergent Lorentz symmetry. It falls into the $(2+1)$-dimensional relativistic Gross-Neveu-Heisenberg quantum universality class, which has recently been much investigated in the context of interacting Dirac systems in two spatial dimensions. The coexistence-to-antiferromagnetic transition, by contrast, turns out to be weakly first order as a consequence of the absence of continuous spatial rotational symmetry on the honeycomb bilayer. Implications for experiments in bilayer graphene are discussed.