
3
Figure 1. a) Graphene unit cell and its lattice symmetries
(top) and its reciprocal unit cell (bottom). b) Phase diagram
at ˜ν=1. It contains four phases: charge density wave (CDW),
Kekul´e distortion (KD) and the two entangled phases, the
antiferrimagnetic phase (AFI) and the canted antiferromagnet
(CAF).
therefore can account for spin-valley entanglement [18].
For ˜ν=1, we take P=F1>< F1, and for ˜ν=2
P=F1><F1+F2><F2.
Ground states for ˜ν=1.As discussed in Ref. [18], the
energy functional in this case reduces to:
E˜ν=1
HF =cos2a1{∆zη2
z+∆⊥η2
⊥},(11)
with ∆z=uH
z−uX
z, ∆⊥=uH
⊥−uX
⊥and η2
⊥=η2
x+η2
y(see
S-V-A) for further details). The resulting phase diagram
is shown in Fig.1(b) and contains four phases. These
are a charge density wave (CDW) with η=ˆz, s=ˆz
and a1=0, and a Kekul´e distortion (KD) state with
η=η⊥,s=ˆzand a1=0. Interestingly, we see that also
spin-valley entangled phases with a1=π2, appear when
∆z>0,∆⊥>0. These entangled phases are degenerate in
the absence of Zeeman fields, but in their presence they
split antiferrimagnetic phase (AFI) with η=ˆz, s=ˆz
and a1=π
2and the canted antiferromagnet (CAF) with
η=η⊥,s=ˆzand a1=π
2, as discussed in Ref. [18].
Notice that in the N=0 LL, ∆z=∆⊥=0, and therefore
all of the above states would be degenerate and with a
vanishing HF energy.
Ground states for ˜ν=2.The HF functional for ν=2 is
more difficult to minimize analytically. To make progress,
we first consider the subset of states from Eq. (10) with-
out spin-valley entanglement. These can be classified into
the valley active states [20] :
F1=η1s,F2=η2−s,(12)
in which the valley degree of freedom varies, and the spin
active states :
F1=ηs1,F2=−ηs2,(13)
in which the spin degree of freedom varies. We first
minimize the energy functional within this subspace and
States appearing at ˜ν=2
States Wavefunctions {∣F⟩1,∣F⟩2}
CDW (Charge density wave) {∣ˆz⟩ ∣s⟩,∣ˆz⟩ ∣−s⟩}
KD (Kekul´e distortion) {∣η⊥⟩ ∣s⟩,∣η⊥⟩ ∣−s⟩}
FM (Ferromagnet) {∣ˆz⟩ ∣s⟩,∣−ˆz⟩ ∣s⟩}
AF (Antiferromagnet) {∣ˆz⟩ ∣s⟩,∣−ˆz⟩ ∣−s⟩
KD-AF (Kekul´e
antiferromagnet)
{∣η⊥⟩ ∣s⟩,∣−η⊥⟩ ∣−s⟩}
Table I. Competing states at ˜ν=2 and their wavefunctions.
then perform a quadratic expansion of all possible devia-
tions of parameters that account for spin-valley entangled
states (see S-V-B), VI, VII for further details). For sim-
plicity we will also neglect the Zeeman term that is typi-
cally weak compared to the interaction terms [2, 21, 22].
In contrast to ˜ν=1, for ˜ν=2 we find that whenever a
spin-valley disentangled state is energetically favorable it
is also an exact local minima of the energy with respect to
all possible quadratic deviations that include spin-valley
entanglement. This indicates that these spin-valley dis-
entangled states are also possibly exact global minima of
the energy.
Following this procedure, we find a total of five pos-
sible ground states for ˜ν=2 that are realized as a
function of the four Hartree and exchange parameters
uH
z, uX
z, uH
⊥, uX
⊥. These possible five states are listed in
Table I (see S-V-B) for more details on these states).
To visualize the energetic competition among these five
phases, we have chosen to draw two-dimensional phase
diagrams as functions of the two Hartree parameters
˜uH
z,⊥=uH
z,⊥∆z−∆⊥for fixed values of ∆z=uH
z−uX
z,
∆⊥=uH
⊥−uX
⊥. We find that there are a total of six
different kinds of phase diagrams depending on the val-
ues and signs of ∆z,⊥. Two of these representative phase
diagrams are depicted in Fig. 2, and the remainder are
presented in S-V-B).
Interestingly, according to the model and the estimates
of Ref. [17], uH
z=uH
⊥=0 and uX
z>0, uX
⊥<0 (see S-IV
for further details). This means that graphene in the
N=1 LL would have a phase diagram like the one in
Fig. 2(a), and it would be located exactly at the origin
of this phase diagram, which we indicate by a black dot
in Fig. 2(a). Therefore, we see that the model and the
parameter estimates of Ref. [17] place graphene right at
the boundary between the CDW and the AF states. We
note that even at this boundary, these phases remain
stable against spin-valley entangled rotations (see S-VII
for further details).
One of the interesting qualitative differences that we
have found in the N=1 LL is the existence of a new
phase that features a combination of Kekul´e state and
antiferromagnet, that we term the Kekul´e- antiferromag-
net (KD-AF). In this phase one set of electrons has an
XY vector in the valley sphere with spin up while the oth-
ers occupy the opposite valley vector with spin down, as