Wigner crystallization in Bernal bilayer graphene

Technical Summary: Wigner crystallization in Bernal bilayer graphene

Problem Statement
The paper investigates the nature and stability of the Wigner crystal (WC) state in Bernal bilayer graphene (BBG) under the influence of a perpendicular displacement field ($U$). While recent experiments in twisted and untwisted graphene multilayers have revealed rich correlated phases, the specific interplay between the unique band structure of BBG (specifically the "Mexican hat" dispersion induced by $U$), Berry curvature, and trigonal warping on the WC phase remains an open theoretical question. The authors aim to determine the phase boundaries of the WC in terms of electron density ($n$), displacement field ($U$), and temperature ($T$), and to characterize the ground state properties, including orbital magnetization and valley ordering, distinguishing them from conventional WC phases.

Methodology
The authors employ a single-particle harmonic oscillator (HO) description to model the WC state. In this framework, each electron in the triangular WC lattice is treated as a particle confined in a parabolic potential well created by the Coulomb repulsion of its neighbors. The strength of this confinement is determined by the electron density.

The study proceeds in two main theoretical regimes:

1. Rotationally Symmetric Model (No Trigonal Warping): The authors first neglect trigonal warping to treat the low-energy dispersion as rotationally symmetric around the $K$ and $K'$ points. They derive an effective single-particle Schrödinger equation in momentum space, incorporating the Berry connection ($A(p)$) as a vector potential. This allows for analytical solutions to estimate energy spectra, wavefunctions, and the Lindemann criterion for melting.
2. Trigonal Warping Model: The authors then include the effects of trigonal warping, which breaks rotational symmetry and splits the "Mexican hat" dispersion into discrete "mini-valleys" (a central pocket and three side pockets). For this case, they perform numerical diagonalization of the two-dimensional Schrödinger equation on a momentum-space grid to compute the Lindemann ratio and map the phase diagram.

The stability of the WC is assessed using the Lindemann criterion, where melting occurs when the root-mean-square width of the electron wavefunction exceeds a critical fraction ($\eta_c \approx 0.23$) of the lattice constant. The authors also calculate exchange integrals to estimate ordering temperatures and evaluate magnetization via the expectation value of the orbital angular momentum operator.

Key Contributions and Results

• Displacement Field Dependence and Dielectric Screening: The paper establishes that at zero displacement field, the WC state is precluded in BBG due to the vanishing band mass at low energy and strong interband dielectric screening. A WC phase only emerges when a sufficient displacement field is applied to open a band gap and flatten the band bottom.
• Spontaneous Orbital Magnetization (No Trigonal Warping): In the absence of trigonal warping, the authors demonstrate that Berry curvature drives a transition in the orbital angular momentum ($\ell$) of the WC electrons. As the displacement field $U$ increases, the Berry flux enclosed by the wavefunction increases. When this flux exceeds a critical value ($\tilde{\Phi} > 1/2$), the ground state transitions from $\ell=0$ to $\ell=1$ (or $\ell=-1$ in the $K'$ valley). This transition manifests as a discontinuous jump in the orbital magnetization of the WC state.
• Doubly Re-entrant Melting (With Trigonal Warping): When trigonal warping is included, the phase diagram exhibits a unique "doubly re-entrant" behavior as a function of electron density at low $U$.
  1. At very low densities, electrons crystallize in the lower-energy "side pockets."
  2. As density increases, the WC melts into a Fermi liquid (FL) at a critical density $n_{side}^c$.
  3. At a higher density $n^*$, it becomes energetically favorable for electrons to occupy the heavier-mass "central pocket," causing the system to re-freeze into a WC state.
  4. Finally, at an even higher density $n_{central}^c$, this central-pocket WC melts again.
• Minivalley Ordering: The presence of trigonal warping breaks the rotational symmetry, rendering angular momentum $\ell$ ill-defined. Instead, the ground state is characterized by a nontrivial spatial pattern of "minivalley ordering." Drawing on recent work (Ref. [102]), the authors predict that electrons polarize into specific mini-valleys in an alternating stripe pattern on the WC lattice.
• Phase Boundaries: The authors provide numerical estimates for the phase boundaries. The maximum density for the WC state in BBG is estimated to be on the order of $1.2 \times 10^{10} \text{ cm}^{-2}$, occurring at a displacement field of approximately $0.7 \text{ V/nm}$. The maximum melting temperature is estimated to be on the order of $10 \text{ K}$ (decreasing with density), with the re-entrant central pocket region having a lower maximum temperature of $\sim 1 \text{ K}$.

Significance and Claims
The paper claims to provide a comprehensive theoretical phase diagram for the WC state in BBG, highlighting how nontrivial quantum geometry (Berry curvature) and lattice symmetry breaking (trigonal warping) fundamentally alter the nature of the Wigner crystal.

• Orbital Effects: The work identifies a mechanism for spontaneous orbital magnetization in a WC driven purely by Berry curvature and displacement fields, distinct from spin-based ordering.
• Complex Phase Behavior: The prediction of "doubly re-entrant" melting is presented as a novel feature arising from the competition between different band pockets (central vs. side) with distinct effective masses and energies.
• Experimental Context: The authors note that their predicted density range for the WC in BBG ($n \lesssim 10^{10} \text{ cm}^{-2}$) is lower than densities where insulating states have been observed in recent experiments, suggesting that current experimental observations in BBG may be dominated by disorder or correspond to different regimes. They suggest that the re-entrant behavior and minivalley ordering might be more readily observable in rhombohedral-stacked graphene with higher layer numbers ($N$), where the band structure is flatter and Berry flux is larger ($N\pi$).
• Complementarity: The paper is explicitly stated to be complementary to a companion work (Ref. [1]) which focuses on spin ordering; this paper focuses on the orbital effects and the detailed phase diagram including trigonal warping.

The authors maintain a modest tone regarding the limitations of the HO approximation, acknowledging that it may overestimate the stability of the WC phase compared to more rigorous quantum Monte Carlo calculations, and that the "Einstein phonon" limit used here cannot reproduce phases with nontrivial Chern numbers (such as the quantum anomalous Hall crystal).