Geometric superfluid stiffness of Kekulé superconductivity in magic-angle twisted bilayer graphene

This paper proposes that a finite-momentum pair-density-wave (PDW) state in magic-angle twisted bilayer graphene reconciles observed tunneling signatures with low-temperature superfluid stiffness suppression by generating Bogoliubov Fermi surfaces, thereby establishing a direct experimental link between zero-bias conductance and phase rigidity.

The Gist

Imagine you have a magical, ultra-thin sheet of graphene that you've twisted at a very specific "magic angle." When you do this, electrons get stuck in a slow-motion dance, forming a superconductor—a material where electricity flows with zero resistance.

For a long time, scientists have been puzzled by this material. They see two conflicting clues that don't seem to fit together:

1. The "Leaky" Clue (Tunneling): When they poke the material with a tiny probe, they see a lot of "leakage." It's as if the superconductor has a hole in its armor, letting low-energy electrons slip through easily. In a normal superconductor, this shouldn't happen; the armor should be solid.
2. The "Stiff" Clue (Superfluid Stiffness): When they try to wiggle the flow of electricity, the material fights back with surprising strength. It's very "stiff" and rigid, resisting changes in flow. Usually, if a material has those "leaky" holes (from clue #1), it should be wobbly and weak, not stiff.

The Big Question: How can the material be full of holes and incredibly strong at the same time?

The Solution: A "Dancing Wave" (The PDW)

The authors of this paper propose a new way to visualize the electrons. Instead of pairing up in a simple, uniform way (like a calm lake), they suggest the electrons are forming a Pair-Density Wave (PDW).

Think of it like this:

• Normal Superconductor: Imagine a crowd of people holding hands and walking in a perfectly straight, uniform line.
• This Magic Graphene: Imagine the crowd is walking in a line, but they are also bobbing up and down in a rhythmic wave pattern as they move. The "pairing" of electrons isn't uniform; it has a specific, wavy rhythm to it.

This wavy pattern is called a Kekulé state (named after a famous pattern in chemistry). It's like a checkerboard that is constantly shifting.

The Magic Ingredient: The "Ghost" Surface

Here is the clever part of the theory. Because of this wavy pairing, the material creates something called a Bogoliubov Fermi Surface.

Let's use an analogy: Imagine a frozen lake (the superconductor). Usually, the ice is solid everywhere. But in this magic material, the wavy pairing creates a hidden, invisible "ghost" path on the ice where the ice is actually thin and slushy.

• The Leak: Electrons can easily slide along this slushy path. This explains the "leaky" tunneling signal.
• The Stiffness: Here is the twist. Even though there is a slushy path, the geometry of the ice around it is so strange and complex that the whole system becomes incredibly rigid.

The authors discovered that the "slushy path" (the gapless electrons) and the "rigid ice" (the stiffness) are actually two sides of the same coin. The very same geometric rules that allow electrons to slip through the slush also create a unique "geometric stiffness" that holds the material together.

The "Traffic Jam" Analogy

Imagine a highway (the material).

• In a normal superconductor, the road is a perfect, empty highway. Cars (electrons) flow smoothly, and if you try to push them, they resist strongly.
• In this magic graphene, there are potholes (the slushy path).
  • Old Theory: If you have potholes, the road should be bumpy and weak.
  • New Theory: The potholes are arranged in a very specific, artistic pattern. This pattern creates a "geometric tension" in the road. Even though cars can drive through the potholes (causing the leak), the road itself is so tightly woven that it refuses to bend when you try to push it.

What This Means for the Future

The paper predicts a specific "fingerprint" for this phenomenon. If you tweak the material (by changing the density of electrons or applying an electric field):

1. The "leakiness" (zero-bias conductance) should go up.
2. The "stiffness" should go down.

Crucially, these two things should move in perfect lockstep. If you see the leak get bigger, you should see the stiffness get weaker, and they should be mathematically linked. This gives experimentalists a clear way to test if their magic graphene is actually doing this "wavy dance."

Summary

The paper solves a mystery by suggesting that the electrons in magic-angle graphene aren't just pairing up; they are dancing in a complex, wavy pattern. This dance creates "holes" that let electrons leak through, but the geometry of the dance is so unique that it simultaneously makes the material incredibly stiff. It's a beautiful example of how a "flaw" (the leak) and a "strength" (the stiffness) can actually be born from the same source.

Technical Summary

1. Problem Statement

Magic-angle twisted bilayer graphene (MATBG) exhibits superconductivity with a complex order parameter that remains unresolved. Experimental observations present a significant contradiction:

• Tunneling Spectroscopy: Experiments (STM) reveal a "V-shaped" density of states (DOS) with substantial residual zero-bias conductance (ZBC), suggesting a nodal or gapless superconducting state.
• Superfluid Stiffness: Conversely, measurements of the superfluid stiffness ($D_s$) show a robust, low-temperature power-law suppression ($\sim T^2$), which is typically associated with gapped or fully coherent states, not gapless ones.

In conventional Bardeen-Cooper-Schrieffer (BCS) theory, the presence of gapless quasiparticles (implied by the V-shaped tunneling) should strongly suppress the superfluid stiffness at low temperatures, creating a tension between the two experimental observables. Previous theoretical attempts using quantum-geometric effects have explained the stiffness but failed to reconcile it with the tunneling phenomenology.

2. Methodology

The authors propose a Pair-Density-Wave (PDW) state, specifically a Kekulé-type superconducting order, as the unifying framework. Their methodology involves:

• Theoretical Model: They utilize the Bistritzer-MacDonald (BM) continuum model for MATBG, supplemented by a short-range attractive interaction that couples A and B sublattices.
• Symmetry Analysis: They focus on an intravalley, unitary spin-triplet PDW state with ordering wave vector $Q = M$ (the M-point of the mini-Brillouin zone). This state spontaneously breaks $C_3$ rotational symmetry but preserves mirror symmetry. Crucially, the pairing occurs within a single valley where time-reversal ($T$) and inversion ($C_2$) symmetries are individually broken, though the global state preserves them.
• Form Factor Analysis: They analyze the particle-particle form factor $\Lambda$, which characterizes the Cooper pair wavefunction. Due to the broken symmetries in the pairing sector, $\Lambda$ becomes highly non-local in band space and momentum-dependent.
• Linear Response Theory: The superfluid stiffness ($D_s$) is calculated explicitly using linear-response theory (London relation) and the curvature of the grand-canonical thermodynamic potential.
• Geometric Decomposition: The stiffness is decomposed into contributions from "active" (flat) bands and "inactive" (remote) bands, incorporating non-Abelian Berry connections and geometric matrix elements.

3. Key Contributions

The paper makes several novel theoretical contributions:

1. Resolution of the Tension: The authors demonstrate that the PDW/Kekulé state naturally reconciles the coexistence of gapless quasiparticles (V-shaped tunneling) and robust superfluid stiffness.
2. Bogoliubov Fermi Surface (BFS): The PDW order induces a Bogoliubov Fermi Surface (BFS)—a codimension-one manifold of gapless quasiparticles in the Brillouin zone. This explains the finite zero-bias conductance and V-shaped spectra.
3. Geometric Origin of Stiffness: A central discovery is that the BFS contributes directly to the superfluid stiffness through PDW-specific quantum-geometric matrix elements. Unlike conventional nodal superconductors where gapless states suppress stiffness, here the geometric structure (involving inter-band pairing between flat and remote bands) allows the gapless states to contribute positively to the stiffness even at $T=0$.
4. Active-Active Channel Contribution: The authors identify a new geometric contribution arising from the "active-active" channel (interactions between flat bands). This term, driven by the non-local form factor $\Lambda$ and the Berry connection, remains significant even in the flat-band limit and dominates the low-temperature behavior.
5. Unified Prediction: They predict a direct, geometric correlation between tunneling spectroscopy and phase rigidity: enhanced residual zero-bias conductance should accompany reduced low-temperature stiffness under density or displacement-field tuning.

4. Key Results

• Bogoliubov Fermi Surface: Numerical calculations confirm the existence of a robust BFS in the PDW state. The gapless dispersion yields a V-shaped DOS, consistent with STM experiments.
• Temperature Dependence of Stiffness: The calculated superfluid stiffness $D_s(T)$ exhibits a suppression proportional to $T^2$ (specifically $D_s(T) \approx D_s(0) - aT^n$ with $n \approx 2.1$) at low temperatures. This matches experimental data from MATBG.
• Mechanism of $T^2$ Suppression: The $T^2$ behavior arises from the geometric contribution of the BFS. The coefficient $a$ is determined by the geometrically weighted density of states $W(E)$ at the Fermi level.
• Non-Linear Order Parameter Dependence: The zero-temperature stiffness $D_{s,0}$ shows a strong non-linear dependence on the order parameter $\Delta$. In the "V-regime" (associated with the BFS), the stiffness grows rapidly with $\Delta$, whereas in the "U-regime" (conventional), it is more linear.
• Correlation between ZBC and Stiffness: The paper establishes that the same low-energy gapless quasiparticles responsible for the finite zero-bias conductance (ZBC) govern the low-temperature suppression of $D_s$. This correlation persists even in the flat-band limit, distinguishing it from conventional nodal superconductivity.
• BKT Transition: The estimated Berezinskii-Kosterlitz-Thouless (BKT) transition temperature $T_{BKT}$ shows an approximately linear relationship with the zero-temperature stiffness, consistent with experimental trends.

5. Significance

• Unification of Observables: This work provides the first unified microscopic mechanism that simultaneously explains the nodal-like tunneling spectra and the robust, $T^2$-suppressed superfluid stiffness in MATBG.
• Role of Geometry: It highlights the critical role of quantum geometry (Berry connections and non-local form factors) in determining the electromagnetic response of superconductors with finite-momentum pairing.
• Experimental Verification: The paper offers a falsifiable prediction: tuning the carrier density or displacement field should simultaneously modulate the zero-bias conductance and the superfluid stiffness in a correlated manner. This provides a direct experimental link between tunneling spectroscopy and phase rigidity.
• New State of Matter: It solidifies the PDW/Kekulé state as a viable candidate for the superconducting order in MATBG, suggesting that the moiré superlattice promotes finite-momentum pairing that fundamentally alters the topological and geometric properties of the superconducting state.

In summary, the paper argues that the "paradox" of gapless tunneling coexisting with robust stiffness is resolved by the unique geometric properties of the Bogoliubov Fermi surface in a Kekulé PDW state, where the gapless modes are intrinsically linked to the superfluid rigidity via quantum geometry.