Full, three-quarter, half and quarter Wigner crystals in Bernal bilayer graphene

Through Hartree-Fock calculations incorporating long-range Coulomb interactions and symmetry breaking, this study reveals that adjusting the displacement field and carrier density in Bernal bilayer graphene induces isospin-polarized metallic phases separated by full, three-quarter, half, and quarter Wigner crystal states that track the Van-Hove singularities.

The Gist

Imagine a bustling city built on a very specific, double-layered grid of graphene (a material made of carbon atoms, like a chicken wire fence). In this city, the "citizens" are electrons. Usually, these electrons zip around freely, like commuters rushing to work. But under certain conditions, they decide to stop rushing and start organizing themselves into a rigid, perfect formation.

This paper is about discovering four different ways these electron citizens can organize themselves into a "Wigner Crystal"—a fancy physics term for a solid crystal made entirely of moving particles.

Here is the story of how they found these crystals, explained simply:

1. The Setup: The "Traffic Jam" and the "Gate"

In this graphene city, the researchers used two main tools to control the traffic:

• The Gate (Displacement Field): Imagine a giant gate that can be opened or closed to change the height of the buildings. When they adjusted this gate, it created a "traffic jam" in the energy levels of the electrons. This is called a Van Hove Singularity. Think of it as a bottleneck where all the cars (electrons) get stuck in one spot.
• The Crowd Density (Carrier Density): They also controlled how many electrons were in the city. By adding just the right amount of "traffic," they forced the electrons to interact with each other more than they interacted with their own speed.

2. The Discovery: The "Isospin" Costumes

Electrons have a secret identity feature called isospin. Think of this as a costume they can wear. There are four types of costumes in total:

1. Spin Up + Valley A
2. Spin Up + Valley B
3. Spin Down + Valley A
4. Spin Down + Valley B

Usually, in a normal metal, the city is a "Full Metal." Everyone is wearing a different costume, and they are all mixed together in a chaotic crowd.

But in this study, the researchers found that when the electrons get stuck in that traffic jam, they start polarizing. This means they decide to stop wearing all four costumes and start wearing only some of them.

• Three-Quarter Metal: Everyone agrees to drop one costume type.
• Half Metal: Everyone drops two costume types.
• Quarter Metal: Everyone drops three costume types, leaving only one.

It's like a party where everyone suddenly agrees to only wear red hats, or only blue hats, creating distinct groups.

3. The Main Event: The Wigner Crystals

Here is the big surprise. Between these "costume parties" (the metallic phases), the electrons decided to stop moving entirely and form a Wigner Crystal.

Imagine the electrons are no longer a chaotic crowd. Instead, they arrange themselves into a perfect, repeating grid, like soldiers standing in formation or tiles on a floor. This happens because they repel each other (like magnets with the same pole facing each other), so they want to get as far apart from each other as possible.

The paper found four specific types of these crystals, depending on how many "costume types" (isospins) were involved:

• Full Wigner Crystal: The crystal is formed by electrons wearing all four costume types.
• Three-Quarter Wigner Crystal: The crystal forms with three costume types.
• Half Wigner Crystal: The crystal forms with two costume types.
• Quarter Wigner Crystal: The crystal forms with just one costume type.

4. Why Does This Matter?

The researchers used a powerful computer simulation (Hartree-Fock calculations) to map out exactly where these crystals appear. They created a "map" (a phase diagram) showing that if you tweak the gate and the crowd density just right, you can switch between a flowing metal and a rigid crystal.

The "Aha!" Moment:
Recent experiments had seen strange things in graphene: high electrical resistance and weird non-linear currents. Scientists were confused. This paper suggests that these weird behaviors are actually the fingerprints of these Wigner crystals.

Specifically, they found that a "Full Wigner Crystal" (where the electrons are organized but not polarized) sits right at the boundary between a "polarized" metal and a "non-polarized" metal. This matches a mysterious high-resistance state observed in labs that seems to be the "parent" of a superconducting state (where electricity flows with zero resistance).

The Analogy Summary

Think of the electrons as dancers in a club:

• Normal Metal: Everyone is dancing wildly, mixing with everyone else.
• Polarized Metal: Everyone agrees to dance in specific groups (e.g., only people in red shirts dance together).
• Wigner Crystal: The music stops, and everyone freezes into a perfect, rigid grid formation to avoid bumping into each other.

The paper shows us that in this specific graphene club, the dancers can freeze into four different types of grids, depending on how many groups (costumes) are present. This discovery helps explain why the club sometimes gets so "stiff" (high resistance) and hints at how it might suddenly start dancing again with zero friction (superconductivity).

Technical Summary

1. Problem Statement

Recent experiments on Bernal-stacked bilayer graphene (BBG) have revealed a rich phase diagram of highly correlated electronic states, including isospin-polarized metals and high-resistance states exhibiting non-linear transport. These high-resistance states, often observed near Van-Hove singularities induced by an out-of-plane displacement field, are hypothesized to be Wigner crystals (WCs)—periodic lattices of charge carriers formed to minimize Coulomb repulsion.

However, the nature of these phases remains debated. While translational symmetry breaking (leading to WCs) has been observed in twisted multilayer graphene and moiré systems, its existence in untwisted BBG at zero magnetic field is less established. Furthermore, the interplay between isospin polarization (spin and valley degrees of freedom) and the formation of Wigner crystals with varying degrees of filling (full, three-quarter, half, quarter) requires a theoretical framework that explicitly allows for the breaking of both translational and rotational symmetries.

2. Methodology

The authors perform self-consistent Hartree-Fock (HF) calculations to investigate the ground state of BBG under varying carrier densities and displacement fields.

• Hamiltonian: The baseline is a nearest-neighbor tight-binding model tailored for BBG around the $K$ and $K'$ points, including parameters for interlayer hopping ($t_1, t_3, t_4$) and an out-of-plane displacement field ($D$).
• Interactions: Long-range Coulomb interactions are included via a double-gate screened potential. To avoid double-counting effects inherent in the tight-binding parameters (which are benchmarked against DFT), a reference projector (representing a completely filled valence band) is subtracted from the density matrix.
• Symmetry Breaking: Unlike previous studies that imposed translational invariance, this work explicitly allows for translational symmetry breaking.
  • The momentum space is folded into a reduced Brillouin zone defined by reciprocal lattice vectors $\mathbf{G}_{WC}$ corresponding to a potential Wigner crystal.
  • The calculation explores ground states formed by completely filling 1 to 4 bands in this reduced zone. This corresponds to holes with 1 to 4 isospin flavors (spin $\times$ valley).
• Numerical Implementation:
  • A $13 \times 13$ momentum mesh is used, folded to include 9 bands per valley/spin (totaling 1521 points).
  • The Optimal Damping Algorithm (ODA) is employed to ensure stable convergence of the self-consistent cycle.
  • The study considers spin collinearity (spin-diagonal states) to reduce computational cost, assuming spin textures are not the primary driver of the observed phases.

3. Key Contributions

• Discovery of Fractional Wigner Crystals: The paper identifies and characterizes four distinct types of Wigner crystals based on the number of occupied isospin flavors: Full (F-WC), Three-quarter (TQ-WC), Half (H-WC), and Quarter (Q-WC).
• Phase Diagram Construction: The authors map the displacement field vs. carrier density phase diagram, revealing regions where translational symmetry is broken, interspersed between isospin-polarized metallic phases.
• Correlation with Isospin Polarization: The study demonstrates that the isospin polarization of the Wigner crystal tracks the polarization of the adjacent metallic phases. For instance, a "half" WC forms in regions where the adjacent metal is a "half" metal.
• Mechanism for High-Resistance States: The work provides a theoretical candidate for the experimentally observed high-resistance "dome" in BBG, suggesting it arises from a Wigner crystal located at the interface between isospin-polarized and unpolarized metallic phases.

4. Key Results

• Phase Diagram:
  • Metallic Phases: The system exhibits Full, Three-quarter, Half, and Quarter metals. The Half metal region is degenerate (spin-polarized, valley-polarized, or spin-valley-locked) due to $SU(2) \times SU(2)$ symmetry.
  • Wigner Crystal Phases: Between these metallic regions, the ground state transitions into Wigner crystals.
    • Q-WC (Quarter WC): Occurs in regions with 1 isospin flavor. It has a large energy gap and breaks translational symmetry via valley and spin polarization channels.
    • H-WC (Half WC) & F-WC (Full WC): Occur with 2 and 4 isospin flavors respectively. They exhibit large energy gaps and break symmetry via Inter-Valley Coherence (IVC) channels.
    • TQ-WC (Three-quarter WC): Unique in that it is gapless. In this phase, holes with two of the three isospin flavors form a WC in a single spin channel, while the third flavor remains metallic.
• Energy Gaps: The calculated energy gaps for the gapped WCs (F-WC, H-WC, Q-WC) are approximately 12 meV, 7 meV, and 4 meV, respectively. These values are significantly larger than thermal energies at typical experimental temperatures (a few Kelvin), suggesting robustness.
• Real Space Profiles: The Wigner crystals exhibit periodicities of a few tens of nanometers. The F-WC profile breaks $C_3$ rotational symmetry, resembling a charge density wave.
• Comparison with Experiments:
  • The F-WC is proposed as the parent state for the high-resistance dome observed in Ref. [1], which precedes superconductivity under an in-plane magnetic field.
  • The H-WC aligns with experimental reports of spin-polarized WCs located between half-metal and quarter-metal phases (Refs. [31, 32]).
• Nematic vs. WC Competition: When translational symmetry is enforced, a nematic phase (broken $C_3$ rotational symmetry) appears near charge neutrality. However, when translational symmetry breaking is allowed, this nematic phase is superseded by a Wigner crystal that preserves $C_3$ symmetry, though the energy competition remains close.

5. Significance

• Theoretical Validation: This work provides strong theoretical evidence that Wigner crystals can form in Bernal bilayer graphene at zero magnetic field, driven by the flatness of the bands near Van-Hove singularities and strong Coulomb interactions.
• Superconductivity Connection: By identifying the high-resistance state as a Wigner crystal, the paper supports the hypothesis that superconductivity in BBG may emerge from a polarized Wigner crystal parent state, similar to mechanisms proposed for rhombohedral tetralayer graphene.
• Experimental Guidance: The prediction of specific energy gaps and real-space periodicities offers concrete targets for experimental verification via scanning tunneling microscopy (STM) and transport measurements.
• Generalization: The framework of "fractional" Wigner crystals (full, 3/4, 1/2, 1/4) based on isospin flavor counting provides a new classification scheme for correlated phases in multicomponent electron systems.

In summary, the paper establishes a comprehensive Hartree-Fock phase diagram for BBG, demonstrating that the interplay between isospin polarization and translational symmetry breaking leads to a hierarchy of Wigner crystal phases that likely underpin the exotic transport and superconducting phenomena observed in recent experiments.