Kekulé Superconductivity in Twisted Magic Angle Bilayer Graphene

Motivated by recent scanning tunneling experiments, this paper proposes a microscopic theory identifying an intra-valley, finite-momentum Kekulé pair-density wave (PDW) as the mechanism for unconventional superconductivity in twisted magic-angle bilayer graphene, a state characterized by spontaneous $C_3$ symmetry breaking, triplet pairing, and a BEC-like regime consistent with experimental signatures.

The Gist

The Big Picture: A Mystery on a Twisted Carpet

Imagine you have two sheets of graphene (a material made of carbon atoms arranged in a honeycomb pattern, like chicken wire). If you stack them perfectly, nothing special happens. But, if you twist one sheet slightly on top of the other, something magical occurs: the atoms create a giant, slow-moving pattern called a Moiré pattern. It looks like the ripples you see when you hold two window screens slightly out of alignment.

Scientists discovered that at a very specific "magic" twist angle, electrons in this twisted carpet stop moving like individual cars and start behaving like a superfluid. They can flow without any resistance. This is superconductivity.

For the last decade, physicists have been trying to figure out how these electrons pair up to do this. It's like trying to figure out the secret handshake that allows a crowd of strangers to suddenly dance in perfect unison.

The New Discovery: The "Kekulé" Dance

This paper proposes a new theory for that secret handshake. The authors, Ke Wang and K. Levin, suggest that the electrons aren't just pairing up randomly; they are forming a specific, rhythmic pattern called a Kekulé order.

Think of the electrons as dancers. In most theories, the dancers pair up across the whole dance floor (inter-valley pairing). But this paper argues that the dancers are pairing up within their own specific group (intra-valley pairing), but they are doing it in a way that creates a wave of movement across the floor.

The Analogy of the Pair-Density Wave (PDW):
Imagine a crowded dance floor. Usually, everyone stands still or moves randomly. In this new theory, the dancers form pairs that are constantly moving forward together, creating a wave of motion. This is called a Pair-Density Wave (PDW). It's like a "conga line" of electron pairs that has a specific rhythm and direction.

Why is this "Kekulé" special?

The name "Kekulé" comes from a pattern of bonds in chemistry that looks like a flower. In this twisted graphene, the electrons arrange themselves so that the "dance floor" itself changes shape slightly, creating a pattern that repeats every three steps.

The paper identifies four key features of this new dance:

1. Breaking the Rules (Nematicity):
   Normally, a honeycomb dance floor looks the same no matter which way you turn it (it has 3-fold symmetry). But this new dance breaks that rule. The dancers spontaneously choose to face one specific direction, making the floor look different if you rotate it. It's like a crowd of people suddenly all deciding to face North, breaking the circular symmetry of the room.

2. The "Triplet" Spin:
   Electrons have a property called "spin" (like a tiny magnet). Usually, superconductors pair electrons with opposite spins (one up, one down), like a classic waltz. This paper suggests these electrons are pairing with spins pointing in the same direction (a "triplet"). It's more like a synchronized swim team moving in perfect unison rather than a traditional couple.

3. The Shape-Shifting Sound (V vs. U):
   When scientists measure the energy of these electrons, they get a graph.

   • The U-Shape: When the attraction between electrons is strong, the graph looks like a deep "U". This means there is a clear gap; no electrons can exist at the lowest energy levels. It's a solid, gapped state.
   • The V-Shape: As the attraction gets weaker, the bottom of the "U" gets sharper and turns into a "V". This is a big deal because a "V" shape usually means the material is still conducting a little bit of electricity even at the lowest energy (it's "gapless").
   • The Twist: The authors explain that this "V" shape isn't a mistake or a defect. It's caused by a Bogoliubov Fermi Surface. Imagine a tiny, hidden island of electrons that exists inside the gap. This explains why experiments see a "V" shape and a specific type of electrical conductance that other theories couldn't explain.

4. The "Short" Dance:
   The paper notes that these electron pairs are very tightly bound and don't travel far before they break up (short coherence length). This suggests the "glue" holding them together is very strong and local, consistent with the idea that they are dancing in a tight, rhythmic pattern rather than a long, flowing wave.

Why Does This Matter?

For years, scientists have been arguing about whether the "glue" holding these electrons together is vibrations in the lattice (phonons) or magnetic fluctuations.

This paper argues that the "Kekulé" pattern isn't just a side effect; it is the engine of the superconductivity.

• The "Glue" vs. The "Machinery": The authors say the Kekulé order isn't the glue itself, but the machinery that organizes the glue. It's the difference between the music (the glue) and the choreography (the machinery). The choreography here is the Kekulé PDW.

The Bottom Line

The authors have built a microscopic model that explains several confusing experimental results:

• Why the material looks different when rotated (Nematicity).
• Why the energy spectrum changes from a "U" to a "V" shape.
• Why there is a specific electrical signal at zero voltage.

They conclude that the electrons in twisted graphene are likely forming a Kekulé Pair-Density Wave. It's a state where electrons pair up with the same spin, break the rotational symmetry of the crystal, and create a hidden "island" of conducting electrons inside the superconducting gap.

In simple terms: They found the choreography. The electrons aren't just pairing up; they are performing a specific, rhythmic, direction-breaking dance that explains all the weird signals scientists have been seeing in their labs. This theory fits the data better than previous ideas and offers a clear path for future experiments to confirm it.

Technical Summary

1. Problem Statement

The nature of superconductivity in twisted graphene systems (specifically magic-angle twisted bilayer graphene, TBG) remains a central unsolved problem in condensed matter physics. While experimental progress has been made in characterizing transport and thermodynamic properties, the microscopic mechanism is debated.

• The Conflict: Recent Scanning Tunneling Microscopy (STM) experiments have observed Kekulé ordering (a $\sqrt{3} \times \sqrt{3}$ reconstruction of the unit cell) within the superconducting and pseudogap phases. However, most existing theoretical frameworks focus on inter-valley pairing (Cooper pairs formed between electrons in different valleys, $K$ and $K'$).
• The Gap: There is a lack of a microscopic theory that explains how intra-valley pairing (pairs within the same valley) can generate the observed Kekulé patterns, nematicity, and specific tunneling signatures (U-to-V transitions) without invoking a Kekulé "glue" mechanism. The authors argue that the Kekulé order observed in STM is a consequence of the pairing machinery itself, not the pairing mechanism.

2. Methodology

The authors develop a microscopic theory based on the Bistritzer-MacDonald (BM) continuum model for twisted bilayer graphene, incorporating a generic short-range attractive interaction.

• Model Framework:
  • Hamiltonian: Uses the BM model with parameters $w_{AA}$ and $w_{AB}$ (interlayer tunneling) and a twist angle $\theta \approx 0.94^\circ$.
  • Interaction: Assumes a short-range, non-local attractive interaction between sites on the honeycomb lattice. The authors argue that screening renders the repulsive Coulomb core short-ranged, allowing for effective attraction.
  • Pairing Channel: Focuses exclusively on intra-valley pairing. The pairing momentum $Q$ is treated as a variational parameter constrained to the mini-Brillouin zone (mBZ).
• Key Theoretical Tools:
  • Quantum Textures: The authors introduce the concept of "quantum textures," defined by momentum- and band-dependent form factors ($\Lambda_{mn}$) derived from the Bloch wavefunctions of the flat bands. Unlike conventional BCS theory where form factors are often unity, these textures are crucial for intra-valley pairing and dictate the symmetry and strength of the order.
  • Symmetry Analysis: The system is analyzed under the $C_3$ rotational symmetry of the lattice. The interaction potential is decomposed into irreducible representations ($A, E_1, E_2$).
  • Thermodynamic Stability: The grand-canonical thermodynamic potential ($\Omega$) is minimized to determine the stable ground state, comparing singlet vs. triplet channels and different pairing momenta.
  • Numerical Implementation: The gap equation is solved self-consistently using a multi-band approach (including remote bands up to $N=20$) to ensure robustness.

3. Key Contributions

1. Identification of Intra-Valley Kekulé PDW: The paper proposes that the superconducting state is an intra-valley, finite-momentum Pair Density Wave (PDW). This state intrinsically carries a Kekulé modulation.
2. Mechanism for Kekulé Order: The authors demonstrate that the Kekulé pattern observed in STM is not a separate order parameter driving pairing, but a secondary order induced by the PDW. Specifically, the bilinear coupling of condensates from the two valleys ($\Delta_{2Q} \Delta^*_{-2Q}$) generates a charge density wave (CDW) with momentum $4Q$, which reduces to a Kekulé $\sqrt{3} \times \sqrt{3}$ modulation on the atomic lattice.
3. Role of Quantum Textures: The study highlights that the "quantum texture" (the structure of the form factors $\Lambda$) is the dominant factor in intra-valley pairing. It strongly suppresses diagonal (intra-band) pairing in the triplet channel while enhancing off-diagonal (inter-band) pairing, leading to a specific thermodynamic preference.
4. Resolution of Experimental Signatures: The theory provides a unified explanation for four distinct experimental features that were previously difficult to reconcile.

4. Key Results

The proposed intra-valley Kekulé PDW state exhibits four salient features that align with experimental data:

• (i) Spontaneous Nematicity ($C_3$ Breaking):
  The ground state condenses at the M point of the mini-Brillouin zone. Since there are three symmetry-related M points, selecting one spontaneously breaks the three-fold rotational symmetry ($C_3$), resulting in nematic order. This contrasts with inter-valley theories where nematicity requires specific combinations of $E$-channels.
• (ii) Triplet Pairing:
  Thermodynamic calculations show that the unitary spin-triplet state has a lower free energy ($\Omega$) than the singlet state at the M point. This is driven by the strong off-diagonal form factors in the triplet channel, consistent with experimental indications of non-singlet pairing (e.g., Pauli limit violation).
• (iii) U-to-V Transition in Tunneling DOS:
  The theory predicts a transition in the quasiparticle density of states (DOS) as the attraction strength varies:
  • Strong Attraction (U-shape): The spectrum is fully gapped.
  • Weaker Attraction (V-shape): The spectrum develops a finite zero-bias conductance.
  • Mechanism: This transition is driven by the emergence of a Bogoliubov Fermi Surface (BFS). When the pairing gap is small, the condition $|\Delta \phi| \leq |\delta \epsilon|$ is met for some momenta, creating gapless regions. This is distinct from nodal superconductors where the DOS vanishes at zero energy; here, the BFS yields a finite zero-bias conductance.
• (iv) Temperature-Dependent Zero-Bias Conductance:
  The calculated zero-bias conductance shows a systematic temperature dependence consistent with experimental observations. The finite conductance is intrinsic to the BFS nature of the state, not merely an artifact of lifetime broadening (Dynes parameter).

Additional Findings:

• The state is found to be near a BEC-like phase (Bose-Einstein Condensation) even with modest interaction strengths, consistent with the extremely short coherence lengths observed experimentally.
• The required interaction energy scale is $\sim 1$ meV, which is experimentally accessible and significantly lower than estimates for inter-valley pairing mechanisms.

5. Significance

• Unification of Phenomenology: This work provides a single microscopic framework that explains the coexistence of Kekulé ordering, nematicity, triplet pairing, and specific tunneling spectra in twisted graphene.
• Paradigm Shift: It challenges the prevailing view that Kekulé order acts as a "pairing glue." Instead, it posits that Kekulé order is a consequence of the intra-valley PDW pairing machinery.
• Experimental Predictions:
  • The theory predicts a finite-q charge modulation at the M point observable in STM, distinct from the $q=0$ signal of inter-valley coherent orders.
  • It offers a clear discriminator between nodal superconductivity (zero zero-bias conductance) and the proposed BFS mechanism (finite zero-bias conductance).
• Robustness: The qualitative results (nematicity, triplet pairing, BFS) are shown to be robust against variations in twist angle, bandwidth, and interaction range, suggesting they are generic features of intra-valley pairing in moiré systems.

In conclusion, the paper identifies the intra-valley Kekulé PDW as the compelling candidate for the unconventional superconductivity in twisted graphene, resolving long-standing discrepancies between theory and the rich phenomenology observed in STM experiments.