Chiral finite-momentum superconductivity in the tetralayer graphene

Motivated by recent experiments on rhombohedral tetralayer graphene, this study utilizes the random-phase approximation to demonstrate that density-density interactions drive chiral $p$-wave pairing, successfully predicting the observed superconducting range and distinguishing between finite-momentum chiral states in SC1 and SC2 and zero-momentum spin-singlet superconductivity in SC4.

The Gist

Imagine a microscopic dance floor made of four layers of graphene (a material as thin as a single atom but incredibly strong). Usually, electrons on this floor just bounce around chaotically. But under very specific conditions—like a gentle electric push and a low crowd density—they suddenly decide to hold hands and dance in perfect unison. This is superconductivity: a state where electricity flows with zero resistance.

This paper is a theoretical investigation into how and why this dance happens in a specific type of four-layer graphene, and it uncovers some very surprising, exotic steps.

Here is the breakdown using simple analogies:

1. The Setting: A Flat, Crowded Dance Floor

The researchers are studying a system where the electrons are moving on a "flat" energy landscape. Imagine a dance floor that is perfectly level. On a flat floor, it's hard to run fast; you tend to shuffle slowly.

• The Problem: Because the electrons are moving so slowly and the crowd is thin (low density), they interact strongly with each other. It's like a crowded room where everyone is standing very close; you can't move without bumping into someone.
• The Question: Can these slow, bumping electrons still form a superconducting dance? And if they do, what kind of dance is it?

2. The Discovery: Two Different Types of Dances

The paper finds that the electrons don't just do one dance; they do different ones depending on how crowded the floor is and how hard the electric "push" is.

Dance Type A: The "Chiral Finite-Momentum" Waltz (SC1 & SC2)

In the low-density regions (not too many electrons), the electrons form a very strange pair.

• The Analogy: Imagine two dancers who are twins (same spin and valley). Instead of standing still and holding hands, they decide to run in a circle together while holding hands. They are moving forward as a pair, even though they are paired up.
• The Twist: This is called "chiral" (meaning it has a specific direction, like a left-handed or right-handed turn).
• The Catch: Because the dance floor is so flat and the electrons are moving so slowly, they are very "wobbly." It's like trying to dance a perfect waltz on a boat in rough seas. The paper suggests that in these regions, the "wobble" (phase fluctuations) is so strong that it might actually stop the superconductivity from happening, even if the electrons want to pair up. This explains why the superconductivity in these areas is fragile and hard to detect.

Dance Type B: The "Zero-Momentum" Tango (SC4)

In the high-density regions (a fuller crowd), the dance changes completely.

• The Analogy: Here, the electrons pair up with their opposites (different valleys). They stand still, holding hands tightly in the center of the floor. They aren't running in circles; they are a stationary, stable unit.
• The Result: This dance is much more stable and doesn't suffer from the "wobbly boat" problem as much.

3. The "Wobbly Boat" Problem (Phase Fluctuations)

One of the paper's biggest contributions is explaining why the superconductivity is so weak in some areas.

• The Theory: Standard physics theories often assume that if electrons want to pair up, they will immediately become superconductors.
• The Reality: The authors show that in this graphene, the "wobble" (phase fluctuations) is a major obstacle.
• The Metaphor: Imagine a group of people trying to march in perfect lockstep. If the ground is shaky (strong fluctuations), they might trip and fall apart before they can march. The paper calculates that in the low-density areas, the ground is so shaky that the "marching" (superconductivity) gets cancelled out, even though the pairing force is strong. This explains why the experimental superconducting temperatures are so low (near absolute zero) compared to what simple math predicts.

4. The Map of the Dance Floor

The authors created a theoretical "map" (a phase diagram) that matches the experimental map almost perfectly.

• SC1 & SC2: The "Wobbly Circle Runners" (Low density, chiral, moving pairs).
• SC3: A transition zone, likely a mix of the two styles.
• SC4: The "Stable Stationary Tango" (High density, spin-singlet, stationary pairs).

Why Does This Matter?

This isn't just about graphene; it's about understanding the rules of the universe for how matter behaves when it's squeezed and cooled down.

1. Unconventional Superconductors: Most superconductors are like the "Stationary Tango" (simple, stable). This paper shows that nature can also do the "Wobbly Circle Run," which is much more exotic and fragile.
2. The Role of Fluctuations: It teaches us that sometimes, the instability of the system is just as important as the strength of the pairing. You can have strong glue holding things together, but if the table is shaking too hard, the glue won't work.
3. Future Tech: Understanding these fragile, exotic states could help us design new materials for quantum computers or ultra-efficient power grids, provided we can find a way to stabilize the "wobbly" dancers.

In a nutshell: The paper explains that in four-layer graphene, electrons can form a superconducting state, but it's a delicate balancing act. In some areas, they try to dance in a moving, spinning circle, but the floor is too shaky, making the dance unstable. In other areas, they dance a stable, stationary tango. The authors successfully mapped out where these different dances happen, matching real-world experiments perfectly.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the theoretical mechanism behind the recently discovered superconductivity in rhombohedral tetralayer graphene (experimental reference [11]). While experiments have revealed exotic properties such as an anomalous Hall effect, time-dependent resistivity fluctuations, and magnetic hysteresis, the underlying pairing mechanism remains unclear.

Key challenges and questions addressed include:

• Low Carrier Density: The system operates at extremely low electron densities ($n \approx 0.5 \times 10^{12} \text{ cm}^{-2}$) with a flat band, suggesting strong electron-electron interactions and potential suppression of superconductivity due to phase fluctuations.
• Pairing Symmetry: It is unclear whether the pairing is spin-singlet or spin-triplet, and whether it occurs within the same valley (intra-valley) or between opposite valleys (inter-valley).
• Finite Momentum: Experimental hints suggest the possibility of finite-momentum pairing (Cooper pairs with non-zero center-of-mass momentum), which is unconventional.
• Phase Diagram Discrepancy: Experiments show disconnected superconducting regions (SC1–SC4) and a high-resistivity region at low densities. Theoretical models need to explain why superconductivity vanishes at low densities despite strong pairing interactions (likely due to phase fluctuations).

2. Methodology and Model

The authors employ a self-consistent mean-field calculation based on the Random Phase Approximation (RPA) to treat density-density interactions.

• Band Structure Model:

  • A tight-binding Hamiltonian ($H_0$) is used to describe the 4-layer graphene system, incorporating 4 layers and 2 sublattices (A/B) per layer.
  • The model includes parameters for interlayer hopping ($\gamma_1, \gamma_2, \gamma_3, \gamma_4$) and the electric displacement field ($u$).
  • The Fermi energy is located in the first conduction band ($\mu_0 = 5$).

• Interaction Treatment:

  • The interaction is modeled as a screened Coulomb potential ($V_{0,q}$) between electrons.
  • The interaction is projected onto the relevant conduction band.
  • RPA Renormalization: The bare interaction is renormalized by the static susceptibility ($\chi_{0,q}$) to obtain the effective interaction: $V_q = V_{0,q} / (1 + \chi_{0,q}V_{0,q})$.
  • The pairing interaction kernel ($V_{P, k'-k}$) includes complex projection factors arising from the orbital wavefunctions.

• Calculations:

  • Gap Equation: Solved self-consistently for the gap function $\Delta_{k,Q}$ with a finite center-of-mass momentum $Q$.
  • Stability Analysis: Comparison of condensation energies ($E_c$) for different pairing symmetries (intra-valley vs. inter-valley, $Q=2K$ vs. $Q=0$) and different spin configurations.
  • Phase Fluctuations: The Berezinskii-Kosterlitz-Thouless (BKT) theory is used to estimate the phase coherence temperature ($T_{BKT}$) based on the superfluid density ($\rho_s$), addressing the limitation of mean-field theory which overestimates $T_c$ by ignoring fluctuations.

3. Key Contributions and Results

A. Identification of Pairing Mechanisms

The study identifies distinct superconducting phases corresponding to the experimental regions SC1–SC4:

• SC1 and SC2 (Low to Intermediate Density):
  • Mechanism: Dominated by chiral $p$-wave pairing ($p_x \pm i p_y$).
  • Configuration: Intra-valley pairing with parallel spins (spin-triplet).
  • Momentum: Finite-momentum pairing with center-of-mass momentum $Q = 2K$ (where $K$ is the Dirac point).
  • Stability: The $p_x \pm i p_y$ state has the lowest condensation energy compared to nematic ($p_x, p_y$) or other chiral states.
• SC3 (Intermediate Density):
  • Mechanism: Chiral inter-valley pairing.
  • Configuration: Occurs between opposite valleys with parallel spins.
  • Momentum: Zero-momentum ($Q=0$).
• SC4 (High Density, $n \approx 1.9 \times 10^{12} \text{ cm}^{-2}$):
  • Mechanism: Spin-singlet pairing.
  • Configuration: Inter-valley pairing between opposite valleys.
  • Momentum: Zero-momentum ($Q=0$).

B. Role of Phase Fluctuations

• Mean-Field vs. Reality: Standard mean-field theory predicts a maximum $T_c \approx 10$ K, significantly higher than the experimental value of $\approx 0.3$ K.
• Suppression Mechanism: The authors attribute this discrepancy to strong phase fluctuations at low electron densities.
• BKT Analysis: Calculations of the superfluid density ($\rho_s$) show that at low densities, the phase coherence temperature ($T_{BKT}$) is drastically reduced. The condition $k_B T_c \leq \frac{\pi}{2} \tilde{D}$ (where $\tilde{D}$ is related to superfluid stiffness) suggests that phase fluctuations can completely suppress superconductivity in the low-density regime, explaining the "high-resistivity" region observed in experiments.

C. Phase Diagram Reconstruction

• The theoretical phase diagram (Fig. 1b) qualitatively reproduces the experimental observations (Fig. 1a).
• It successfully maps the regions of SC1, SC2, SC3, and SC4 as functions of electron density ($n$) and electric displacement field ($u$).
• The study highlights that the transition from SC3 to SC2 as density increases corresponds to a switch from inter-valley ($Q=0$) to intra-valley ($Q=2K$) pairing, driven by the evolution of valley polarization and interaction strength.

D. Complex Projection Factors

• The inclusion of complex projection factors in the pairing interaction ($V_{P, k'-k}$) preserves the chiral $p$-wave symmetry but reduces the pairing amplitude.
• This reduction lowers the mean-field $T_c$ by approximately half, bringing theoretical estimates closer to experimental values, though phase fluctuations remain the dominant suppression factor.

4. Significance and Implications

• Mechanism Resolution: The paper provides a unified theoretical framework explaining the diverse superconducting phases in rhombohedral tetralayer graphene, linking them to specific spin, valley, and momentum configurations.
• Finite-Momentum Superconductivity: It offers strong theoretical evidence for chiral finite-momentum superconductivity (a form of Fulde-Ferrell-Larkin-Ovchinnikov state) in a graphene-based system, a rare phenomenon in condensed matter physics.
• Role of Fluctuations: It underscores the critical role of phase fluctuations in low-density, flat-band systems, explaining why strong pairing interactions do not always lead to high-temperature superconductivity.
• Time-Reversal Symmetry Breaking: The proposed chiral finite-momentum pairing in SC1/SC2, combined with the spin-valley polarized metallic state, suggests a mechanism for time-reversal symmetry breaking in the normal state, consistent with experimental anomalous Hall signals.
• Guidance for Future Research: The identification of the $r_s$ parameter (ratio of interaction to kinetic energy) suggests that in the $r_s > 40$ regime, strong correlation effects (like Wigner crystallization) may compete with superconductivity, offering a roadmap for future experimental tuning.

In summary, Qin and Wu successfully model the complex phase diagram of tetralayer graphene, attributing the observed phenomena to a competition between chiral finite-momentum pairing, valley polarization, and strong phase fluctuations.