Band-basis decomposition of superfluid weight in magic-angle twisted bilayer graphene: Quantifying geometric and conventional contributions

This paper decomposes the superfluid weight in magic-angle twisted bilayer graphene into conventional and geometric contributions, revealing that while remote bands exclusively enhance the geometric component to approximately 55–58%, the geometric fraction peaks at 27–33% near the strongest superconducting fillings of $\nu = \pm 2$.

The Gist

Imagine you have a very special, ultra-thin sandwich made of two layers of graphene (a material as thin as a single atom) twisted at a very precise "magic" angle. Scientists have discovered that when you cool this sandwich down, it becomes a superconductor—a material that conducts electricity with zero resistance, like a frictionless slide for electrons.

The big mystery has been: How does this work?

Usually, for electricity to flow super smoothly, the electrons need to be moving fast (like cars on a highway). But in this magic-angle sandwich, the electrons are stuck in "flat" energy bands. Think of these bands as a perfectly flat, featureless parking lot. If you're in a flat parking lot, you can't accelerate; your speed is zero. So, how can these "stopped" electrons create a supercurrent?

This paper solves that mystery by breaking down the "superfluid weight" (a measure of how strong the supercurrent is) into two distinct ingredients. The authors used a mathematical microscope to separate the total strength into:

1. The "Conventional" Part (The Engine): This is the part that comes from electrons actually moving. In a normal superconductor, this is the main driver.
2. The "Geometric" Part (The Map): This is the new, weird part. It comes from the shape and twist of the electron's path through the material's quantum landscape, even if the electrons aren't moving fast.

The Analogy: The Hiker and the Mountain

Imagine you are a hiker trying to get to a campsite (the superconducting state).

• The Conventional Way: You run fast along a flat road. Your speed (velocity) gets you there. In the magic-angle graphene, the "road" is flat, so you can't run. Your speed is zero.
• The Geometric Way: Even if you stand still, the shape of the world around you helps you. Imagine the ground is a giant, invisible trampoline with a specific pattern of bumps and dips (the "Quantum Geometry"). If you wiggle in the right way, the shape of the trampoline itself pushes you forward, even if you aren't running.

What Did They Find?

The researchers did a detailed accounting of how much "push" comes from running (Conventional) versus how much comes from the shape of the trampoline (Geometric).

1. The Flat-Band Surprise (The Parking Lot)
When they looked only at the "flat" bands (the parking lot where electrons are stuck), they found that the "Geometric" push accounts for about 22% to 26% of the total supercurrent.

• Why this matters: In a normal world, if you have no speed, you have no current. Here, the shape of the material provides a quarter of the power! It's like a car that can drive just by the design of its wheels, even with the engine off.

2. The Hidden Helpers (Remote Bands)
The researchers then asked, "What if we include the electrons in the layers above and below the main flat band?" (Think of these as "remote bands").

• When they added these extra layers, the "Geometric" contribution jumped to about 55% to 58%.
• The Twist: The "Conventional" part (the running speed) didn't change much. It stayed the same. All the extra power came from the Geometric part.
• The Metaphor: It's like realizing that while you were standing still in the parking lot, the entire neighborhood was actually vibrating in a way that was pushing you forward. The "remote" electrons act like a hidden chorus that amplifies the geometric push.

3. The Sweet Spot
They found that this geometric effect is strongest right where superconductivity is actually observed in experiments (near specific electron fillings). It's like the "magic" only happens when you have just the right number of hikers on the trail.

Why Should You Care?

For years, scientists knew that "geometry" played a role in these materials, but they didn't know how much. They were guessing whether it was a tiny 1% or a huge 50%.

This paper puts a number on it. It proves that in magic-angle graphene, the shape of the quantum world is just as important as the speed of the electrons.

• The Takeaway: Superconductivity in these materials isn't just about electrons zooming around; it's about the electrons dancing to the rhythm of the material's twisted structure. The "geometry" of the material is doing half the heavy lifting.

This discovery helps scientists design better superconductors in the future. Instead of just trying to make electrons move faster, we might be able to engineer materials with the perfect "quantum shape" to make electricity flow without resistance, even at higher temperatures.

Technical Summary

1. Problem Statement

Magic-angle twisted bilayer graphene (MATBG) exhibits superconductivity in nearly flat bands. A central puzzle in this field is understanding how a large superfluid weight ($D_s$) is sustained when the band velocity (which typically drives conventional superfluidity) vanishes in flat bands.

• Theoretical Context: Previous work established that "quantum geometry" (specifically the quantum metric) provides a lower bound for $D_s$ in flat-band systems. Experimental measurements (e.g., by Tian et al. using cQED) show $D_s$ is roughly 10 times larger than conventional BCS estimates, suggesting a strong geometric enhancement.
• The Gap: Despite qualitative understanding, there has been no systematic, quantitative decomposition of $D_s$ into conventional (band-velocity) and geometric (interband coherence) components for MATBG. Existing studies lacked explicit dependence on pairing symmetry, chemical potential, gap magnitude, and the number of included bands (truncation effects).

2. Methodology

The authors employed the Bistritzer-MacDonald (BM) continuum model with relaxation-corrected tunneling parameters to simulate MATBG.

• Current Operator Splitting: They applied a band-basis decomposition to the current operator. By diagonalizing the normal-state Hamiltonian $H(k)$, they split the velocity operator into:
  • Intra-band ($v^{intra}$): Diagonal terms representing band velocities ($\partial \epsilon_n / \partial k$).
  • Inter-band ($v^{inter}$): Off-diagonal terms proportional to the Berry connection ($A_{nm}$) and energy differences.
• Decomposition of $D_s$: Substituting these into the Kubo formula for superfluid weight yielded three distinct terms:
  $$D_s = D_s^{conv} + D_s^{geom} + D_s^{cross}$$
  • $D_s^{conv}$: Conventional contribution (intra-band).
  • $D_s^{geom}$: Geometric contribution (inter-band coherence).
  • $D_s^{cross}$: Interference terms between intra- and inter-band currents.
• Computational Strategy:
  • Model: BM model with $N_{shell}=3$ reciprocal lattice shells ($196 \times 196$ matrix).
  • Pairing Symmetries: Investigated uniform $s$-wave, sublattice $s$-wave, and nematic $d$-wave pairings.
  • Band Truncation: To handle the ultraviolet (UV) divergence inherent in continuum models (where $D_s^{conv}$ diverges as more bands are added), they compared:
    1. Flat-band projection ($n_{keep}=2$): Retaining only the two flat bands.
    2. Extended truncation ($n_{keep}=2 \to 6$): Including remote bands to test convergence.
  • Integration: Used a centered $14 \times 14$ Monkhorst-Pack mesh to avoid divergences at the $\Gamma_M$ point.

3. Key Contributions

• First Quantitative Decomposition: Provided the first Brillouin-zone averaged, multi-band convergence study separating geometric and conventional contributions in MATBG.
• Resolution of UV Divergence: Demonstrated that while the total $D_s$ diverges in the continuum model with added bands, the conventional part ($D_s^{conv}$) converges rapidly (within 2%), while the geometric part ($D_s^{geom}$) absorbs all additional weight from remote bands.
• Symmetry Analysis: Proved that for symmetric gap functions (uniform $s$-wave, nematic $d$-wave), the cross term $D_s^{cross}$ vanishes to machine precision due to the orthogonality of Pauli matrices in Nambu space.

4. Key Results

A. Flat-Band Limit ($n_{keep}=2$)

At charge neutrality ($\mu=0$) with a gap $\Delta_0 = 1$ meV:

• Geometric Fraction: Quantum geometry accounts for 22–26% of the total superfluid weight, depending on pairing symmetry.
  • Uniform $s$-wave: 21.5%
  • Sublattice $s$-wave: 26.2% (Highest, due to sublattice polarization).
  • Nematic $d$-wave: 24.4%
• Absolute Contribution: The geometric term contributes 13–16 eV Ų, a substantial value that would be zero in a dispersive-band superconductor with the same density of states.
• Cross Terms: Vanish ($<10^{-13}$) for uniform and nematic pairings; negligible ($<0.01\%$) for sublattice $s$-wave.

B. Effect of Remote Bands ($n_{keep} \to 6$)

• Convergence: $D_s^{conv}$ converges rapidly to $\approx 54$ eV Ų regardless of the number of remote bands included.
• Geometric Enhancement: Including remote bands increases the total $D_s$ significantly, with the geometric fraction rising from 21.5% to 58.3%.
• Physical Insight: Remote bands contribute exclusively through interband coherence (geometric terms), not conventional band velocity. This implies that simplified flat-band-only models systematically underestimate the geometric contribution.

C. Dependence on Filling and Gap

• Filling ($\nu$): The geometric fraction is not constant. It peaks at ~33% near the superconducting domes ($\nu \approx \pm 2$, specifically $\mu \approx \mp 4$ meV) where the Fermi level intersects regions of high quantum metric ($\Gamma_M, K_M$ points). It drops to 9–15% at high doping where flat bands are fully occupied/empty.
• Gap Magnitude ($\Delta_0$): In the experimentally relevant range ($0.3 - 1.0$ meV), the geometric fraction is relatively insensitive to $\Delta_0$ (decreasing slightly from 32% to 22% as $\Delta_0$ increases). The ratio remains roughly constant as long as $\Delta_0 \ll W$ (bandwidth).

D. Pairing Symmetry

• $D_s^{geom}$ is relatively stable across pairing symmetries (determined by the normal-state band geometry).
• $D_s^{conv}$ varies significantly (~34%) depending on how the gap modulates the quasiparticle spectrum (e.g., uniform $s$-wave yields the highest conventional weight).

5. Significance and Implications

• Quantifying the "Anomaly": The study quantifies that geometric effects contribute a factor of 1.27–1.35 to the superfluid weight in the flat-band limit, rising to ~2.4 when remote bands are included. While this explains part of the experimental observation ($D_s/D_s^{BCS} \approx 10$), the remaining enhancement is likely due to interaction-renormalized quantum metrics, vertex corrections, and strong-coupling effects.
• Model Validity: The results validate the use of flat-band projections for capturing intrinsic physics but warn that they underestimate the total geometric contribution by a factor of ~2.5 if remote bands are ignored.
• Future Directions: The paper establishes a quantitative baseline, motivating the need for UV-complete tight-binding calculations with self-consistent pairing to fully resolve the magnitude of the geometric enhancement in MATBG.

In summary, this work rigorously dissects the origin of superfluid stiffness in MATBG, confirming that while conventional band velocity provides the baseline, quantum geometry is a dominant, non-negligible contributor that is further amplified by interband coherence with remote bands.