Kekul\'e Superconductivity in Twisted Magic Angle Bilayer Graphene

Ke Wang; K. Levin

REVIEW 5 major objections 6 minor 1 cited by

Kekul\'e Superconductivity in Twisted Magic Angle Bilayer Graphene

T0 review · 5 major / 6 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper claims that twisted bilayer graphene's superconductivity is an intra-valley finite-momentum pair-density wave with intrinsic Kekulé order, unifying nematicity, triplet pairing, tunneling spectra, and short coherence lengths.

desk verdict A serious, internally coherent candidate theory for Kekulé PDW superconductivity in TBG, but the thermodynamic selection of Q=M is asserted rather than shown and the pairing glue is an assumed input. read the letter →

arxiv 2510.06451 v3 pith:XOSHCD72 submitted 2025-10-07 cond-mat.supr-con cond-mat.mtrl-scicond-mat.str-el

Ke Wang , K. Levin This is my paper

classification cond-mat.supr-concond-mat.mtrl-scicond-mat.str-el

keywords twistedbilayergrapheneKekulésuperconductivitypair-densitywaveintra-valleypairingnematicordertripletgaplessquasiparticleFermisurfaceBEC-BCScrossover

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper sets out to show that the superconductivity in twisted bilayer graphene is not the conventional intervalley-paired superconductor but an intra-valley pair-density wave: Cooper pairs form within a single valley with a finite center-of-mass momentum at the edge of the moiré Brillouin zone, which imprints a Kekulé (√3×√3 bond-order) modulation on the charge density. Starting from a continuum model of the two flat moiré bands plus a short-range nearest-neighbor attraction, the paper solves the gap equation self-consistently and finds this Kekulé pair-density wave to be the stable ground state, with spin-triplet pairing and spontaneous breaking of the lattice's threefold rotation symmetry, making the order nematic. The same state produces a V-shaped-to-U-shaped evolution of the density of states as the attraction grows, a finite zero-bias conductance tied to a small gapless quasiparticle Fermi surface, and a proximity to a BEC-like regime at modest coupling, matching the very short coherence lengths seen in experiments. If correct, it would unify several observed features—nematicity, the anomalous upper critical field, tunneling lineshapes, and short coherence lengths—as consequences of one order parameter.

What carries the argument

The load-bearing object is the particle-particle form factor Λ_mn(Q,q), which the paper calls the 'quantum texture': products of the flat-band Bloch wavefunctions that enter the intra-valley gap equation. Unlike conventional intervalley pairing, where such form factors are essentially unity, here the diagonal components deviate from one and the off-diagonal components are large, so the Bloch wavefunction texture directly dictates the momentum, band, and spin character of the order parameter. Together with the combined twofold-rotation/time-reversal symmetry, this structure suppresses diagonal triplet pairing and routes pairing through off-diagonal bands, while a threefold rotation acts as a

What would settle it

Compute the grand-canonical potential as a continuous function of the pairing momentum Q across the first mini-Brillouin zone at fixed chemical potential and interaction strength, and locate its global minimum: if any Q other than the M point has lower Ω, the proposed Kekulé pair-density-wave ground state is ruled out. The same calculation should be repeated for both singlet and triplet channels at the attraction strengths used in the paper.

Watch

Extended reading notes

Core claim

The central discovery is that, within a continuum model of twisted bilayer graphene's flat moiré bands and a nearest-neighbor intralayer attraction, the self-consistent superconducting ground state is an intra-valley pair-density wave—pairs formed within a single valley—with pairing momentum at the M point of the mini-Brillouin zone. This state intrinsically carries a Kekulé modulation, spontaneously breaks C3 rotation symmetry, and is spin-triplet, with the two layers carrying equal condensate amplitudes. A 'quantum texture'—products of the flat-band Bloch wavefunctions entering the pairing form factor—controls the gap equation: off-diagonal band elements dominate and the diagonal triplet e

Load-bearing premise

The load-bearing premise is that electrons within each layer feel a short-range nearest-neighbor attraction of order a few meV; if the real pairing force instead binds electrons from opposite Dirac cones (valleys), the intra-valley Kekulé pair-density-wave picture collapses.

Editorial extensions

If this is right

The superconducting order is generically nematic: even with a fully isotropic attraction, the condensate breaks the lattice's threefold rotational symmetry, so transport and scanning-tunneling measurements should reveal orientation-dependent order.
The order is spin-triplet, which directly explains the experimentally observed violation of the conventional paramagnetic limit—superconductivity surviving in-plane magnetic fields far beyond the weak-coupling bound.
Tunneling spectroscopy should show a V-shaped density of states with a finite zero-bias conductance in the weaker-attraction regime, produced by an intrinsic small gapless quasiparticle Fermi surface, and a fully gapped U-shaped spectrum without zero-bias conductance at stronger attraction.
The Kekulé modulation seen by STM in the superconducting state is a secondary charge-density-wave induced by the pair-density-wave condensates, not the pairing glue, so its measured strength should track pairing strength, as reported.
The proximity to a BEC-like regime at modest coupling means the pairing gap is large relative to the Fermi energy, giving the very short coherence lengths observed and predicting a crossover that should show up in superfluid-stiffness measurements.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Extension: if the C3-breaking mechanism is as generic as the paper's proof suggests, then observing nematic superconductivity in any moiré flat-band material becomes evidence that its pairing is intra-valley rather than intervalley.
Extension: the V-to-U criterion implies a controllable crossover—lowering the order-parameter magnitude by doping, field, or temperature should push the same sample from a gapped U-shaped spectrum into a V-shaped spectrum with finite zero-bias conductance near the gap amplitudes shown in the paper.
Extension: applying the same two-flat-band logic to twisted trilayer graphene suggests the trilayer should reach the U-shaped regime while the bilayer remains V-shaped; an explicit three-band calculation would test this expectation.
Extension: the predicted near-BEC regime implies that phase fluctuations, rather than the single-particle gap, may set the measured transition temperature, so superfluid-stiffness measurements should show a pronounced suppression as T_c is approached.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The paper is worth a real look. What is genuinely new here is the TBG-specific machinery: band-projected form factors from the BM flat bands, the singlet/triplet competition evaluated at Q=M, the explicit C3-symmetry-breaking proof, and the two-band Bogoliubov analysis that produces the V-to-U DOS crossover. These are computed, not assumed. The C3-SSB proof is particularly clean because it uses only the intra-valley nature of the order and the layer-dependent momentum boost, independent of the interaction details. The DOS and zero-bias conductance results are also concrete predictions that experiment could test. The paper is honest about its own limits—ref. [47] is cited to refute a full BEC reading, and the two-valley independence and flat-band truncation are flagged. That is integrity.

The soft spots are real but proportionate. The biggest one is exactly what the stress-test notes: the central claim that the stable pairing momentum is Q=M is asserted, not demonstrated. The text says for given μ and spin, Ω is minimal at M, but no Ω(Q) scan is shown anywhere. Fig. 3 only compares singlet vs triplet at fixed Q=M. That means the paper has established that a Q=M triplet solution exists and is locally stable at that Q, but not that it is the global minimum in the full variational space. If a singlet at a different Q, or a triplet at another Q, is lower, the main conclusion loses its footing. This is not an external-consensus objection; it is a missing internal calculation. It should be fixable, but it is load-bearing.

Second, the meV-scale nearest-neighbor attraction is an input, not derived. The paper argues the essential features are generic to short-range interactions, and the C3 proof supports that, but the actual gap scale and the quantitative triplet/singlet ordering will depend on the form of V. The agreement with experiment is also qualitative—shape-level, no direct overlays of the ZBC data from refs. [45,50]. That is fair for a theory paper, but it sets the confidence level.

The citation pattern is honest. They cite the prior Kekulé superconductivity work [21-23] and the intervalley theories they are contrasting with. No red flags there.

Who is this for? People working on the mechanism and order parameter of superconductivity in twisted graphene. It is a strong candidate to put on the table, but it is not a settled result. I would send it to peer review, and I would ask the referee to insist on the Ω(Q) scan and a clearer statement about the origin and range of validity of the short-range attraction.

Bottom line: a serious paper with one missing verification that could break the central claim. Worth engaging, but not yet a definitive resolution.

Referee Report

5 major / 6 minor

Summary. The paper develops a mean-field theory of intra-valley, finite-momentum pairing in twisted bilayer graphene. Starting from the Bistritzer\u2013MacDonald continuum model supplemented by an assumed short-range nearest-neighbor attractive interaction, the authors derive a spectral-basis gap equation with non-trivial form factors ('quantum textures'), argue that the stable solution has equal layer condensates, pairing momentum Q=M, and spin-triplet symmetry, and prove that the condensate spontaneously breaks C3. They then analyze the quasi-particle density of states and identify a U-to-V crossover as a function of gap size, together with temperature-dependent zero-bias conductance behavior, and interpret these as explaining STM, Pauli-limit, nematicity, and coherence-length observations in the twisted graphene family.

Significance. If the central claim holds, the paper provides a single order\u2014an intra-valley spin-triplet pairing at Q=M with a Kekul\u00e9 modulation\u2014that unifies several experimental puzzles: Kekul\u00e9 order in the superconducting state, Pauli-limit violation, nematicity, short coherence lengths, and V/U-shaped tunneling spectra. The formal core is coherent: the spectral-basis gap equation, the explicit singlet/triplet grand-potential comparison at fixed Q, the two-band Bogoliubov dispersion with a clear gap-closing condition, and the C3-symmetry-breaking argument are all internally consistent. The DOS and zero-bias behavior are concrete falsifiable predictions. However, two load-bearing points are not yet established: the global-minimum property of Q=M is asserted without a variational scan, and the microscopic pairing attraction is put in by hand. These currently limit confidence in the central claim.

major comments (5)

[Establishing the Thermodynamically Stable Phase] The central variational claim is unsubstantiated: the text states that for fixed \mu and spin channel, \Omega is minimal at Q=M, but Fig. 3 compares singlet and triplet only at Q=M. No \Omega(Q) scan, analytic minimization, or stability analysis over Q is shown, and the text explicitly says the singlet/triplet ordering depends on Q. Therefore the conclusion that the global ground state is a triplet Q=M PDW is conditional. Provide \Omega(Q) for singlet and triplet along high-symmetry paths for representative \mu and V, and state where the global minimum lies.
[Eqs. (2)-(4), 'Microscopic Model'] The pairing mechanism is an assumed intralayer nearest-neighbor attraction of strength V\u22481\u20135 meV, with no microscopic origin. This assumption selects the intra-valley channel, fixes the gap scale, and determines the BEC-proximity claim. Since most previous TBG mechanisms are intervalley, the paper should either derive this attraction from a concrete microscopic process or show robustness across a range of short-range interaction forms and strengths. Absent that, the 'microscopic theory' is a model study; this is a correctness risk for the central claim.
[Zero-Bias Conductance, Figs. 6 and 8] The 'zero-bias conductance' is effectively the DOS at \omega=0, not a tunneling conductance. dI/dV for a tunnel junction includes matrix elements, barrier transparency, and Andreev processes, so it is not equal to DOS(0). The claimed agreement with experimental ZBC curves (Refs. [45,50]) is therefore not yet demonstrated. Please compute a conductance with a stated tunneling model, or restrict the claim to a qualitative DOS prediction.
[Establishing the Thermodynamically Stable Phase (layer branches)] The assertion that \Delta_1=\Delta_2 is the global minimum while \Delta_1=-\Delta_2 is a saddle point is made without supporting calculation. All results use the equal-layer branch; if its stability is not checked as a function of V and Q, the triplet and DOS conclusions may not be robust. Provide the grand-potential comparison or a stability analysis for both mirror branches.
[PDW ORDER AND MICROSCOPIC MODEL, Eq. (1)] The label 'finite-momentum pair-density wave' needs clarification. The text defines the center-of-mass momentum as 2Q and later finds Q=M; in the mini-BZ 2M is a moir\u00e9 reciprocal lattice vector, so the condensate has zero crystal momentum modulo G and is commensurate with the moir\u00e9 lattice. The Kekul\u00e9 modulation is at the atomic scale, not a long-wavelength gap modulation. Unless the order parameter breaks moir\u00e9 translation symmetry, the state should be described as a commensurate Kekul\u00e9 SC rather than a PDW, or the definition of PDW should be explicitly adjusted.

minor comments (6)

[Nematic Nature of Superconducting order] Typo: 'order order' should be 'order'.
[U-to-V transitions in tunneling] 'Landau Ginsberg coherence length' should be 'Landau\u2013Ginzburg coherence length'.
[Methods / Flat band Approximation] Main text says a two-flat-band approximation is adequate; Methods says 10 bands are retained in practice. Please specify which calculation (gap, DOS, ZBC) uses which truncation.
[Data availability] Code and data 'available on reasonable request' is not a stable repository; for reproducibility, an archived version would be preferable.
[Fig. 6 caption] Caption says 'Density of states and order parameter versus temperature' while text calls it zero-bias conductance; unify the terminology.
[Introduction, Ref. [24]] Footnote 24 is an extended caveat in the reference list; moving it to a main-text note or appendix would improve readability.

Circularity Check

1 steps flagged · score 4.0 of 10

Kekulé modulation restates the chosen intra-valley pairing ansatz; quantitative predictions remain self-contained.

self definitional [Introduction ('methodology for implementing superconducting Kekulé order') and 'PDW ORDER AND MICROSCOPIC MODEL', Eq. (1); cf. Abstract.]
"Motivated by this STM data, we now consider intra-valley pairing superconductivity in twisted bilayer graphene (TBG), which is a methodology for implementing superconducting Kekulé order [21–23]. ... On the atomic lattice this reduces to a modulation ∝ cos(2K·r) (since 3K is a reciprocal vector), yielding a Kekulé √3×√3 superlattice."

The theory's input is intra-valley pairing, which the paper itself calls 'a methodology for implementing superconducting Kekulé order,' with Q 'naturally expected to lie near the Dirac point K.' The headline result—'a PDW that intrinsically carries a Kekulé modulation'—is then not an independent output of the self-consistent calculation: Eq. (1) shows that the charge response of two intra-valley condensates at ±2Q is ρ∝cos(2K·r) algebraically, because 3K is a reciprocal vector. The reported 'intrinsic Kekulé modulation' is therefore the assumed intra-valley-PDW ansatz expressed in the charge channel, i.e. a restatement of the input rather than a derived prediction. The remaining results (singlet/triplet Ω comparison at fixed Q=M, C3-SSB proof from BM form factors, DOS U/V crossover) are ge

full rationale

The circularity is confined to the labeling/identification of the central state. The paper begins by assuming intra-valley pairing, explicitly calling it 'a methodology for implementing superconducting Kekulé order,' and the finite-momentum pairing is taken with Q near the Dirac point. The charge modulation of such a two-valley PDW is then shown in Eq. (1) to be cos(2K·r), i.e. the Kekulé √3×√3 pattern. Consequently, the abstract's statement that the pairing 'intrinsically carries a Kekulé modulation' is a definitional consequence of the ansatz, not a derivation of Kekulé order from the microscopic interaction. However, most quantitative content is independent and self-contained: the singlet-versus-triplet grand-potential comparison at fixed Q=M, the C3-symmetry-breaking proof from BM wavefunction form factors, and the self-consistent DOS U/V crossover and zero-bias conductance curves are computed, not fitted to the target Kekulé feature. The Q=M global-minimum assertion is not supported by a shown Ω(Q) scan, but that is a missing variational verification, not circularity. One self-citation (ref. [55]) is used to argue that Kekulé order is intertwined with pairing rather than acting as glue, but this is interpretational and not load-bearing for the central derivations. Overall, the central identification is partly circular, while the ancillary predictions carry genuine independent content; a score of 4 reflects this partial, non-dominating circularity.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The 'quantum texture' is a name for the band-projected form-factor matrix Λmn, not a new entity; the paper introduces no new particles, forces, or conserved quantities. The main ledger burden is the assumed attractive interaction (ad hoc input) and the chosen chemical potential/filling, plus the asserted Δ1=Δ2 branch selection.

free parameters (3)

Attractive interaction strength V ≡ V0(q=0) = Scanned from -0.8 to -5 meV; BEC crossing at V ≈ -1.7 meV
The nearest-neighbor attraction is the model's key input; its magnitude is tuned so the self-consistent |Δ| lands in the experimentally observed range (~0.3–0.73 meV), and the U/V transition and BEC-proximity claims are presented as functions of this parameter.
Chemical potential μ = -0.59 meV (3/8 filling of the lower flat band, claimed ν = -2.5)
Chosen by hand to correspond to experimental filling ν=-2.5; the DOS curves, singlet/triplet Ω comparison, and BEC criterion all depend on this choice.
Retained band number in numerics = 10 bands
Truncation choice for the self-consistent gap equation; the paper argues the two-flat-band limit is adequate but retains 10 bands 'for precision', and grid/convergence details are not stated.

assumptions (6)

ad hoc to paper A short-range nearest-neighbor attractive interaction of strength ~1 meV exists within each TBG layer and drives the pairing.
Introduced to enable the intra-valley pairing channel; no microscopic origin (glue) is specified or derived. The paper is explicit that pairing 'machinery' rather than 'mechanism' is treated (PDW Order and Microscopic Model section; footnote 24).
domain assumption The two valleys contribute independently; a single-moiré-valley analysis suffices.
'That the two valleys contribute independently is well-justified for small twist angles' (PDW Order and Microscopic Model section).
domain assumption The A-channel (l=0) dominates the interaction because chiral form factors vanish at small momentum transfer; f0(q)→1.
Used to project V(q−q') onto V0(q)=g0 f0*(q) and restrict to AB-bonding order; relies on the smallness of |q·δj| relative to the atomic scale (footnotes 39–40).
domain assumption Two isolated flat bands capture the pairing physics; remote bands can be integrated out (Δ/E_band ≪ W/Δ).
Flat Band Approximation section; the two-band reduction underlies Eq. (10) and the U/V analysis. 10 bands are retained numerically 'for precision'.
ad hoc to paper The physical ground state has equal layer condensates Δ1=Δ2; the Δ1=−Δ2 branch is a saddle point.
Stated in 'Establishing the Thermodynamically Stable Phase' without a supporting calculation; this branch selection enters the subsequent singlet/triplet and DOS analysis.
standard math C3 and C2T symmetries of the BM Hamiltonian classify the pairing; the full state is time-reversal-symmetric via the two-valley combination [Δ2Q, Δ−2Q].
Standard BM model symmetries (Methods section); used for the C3-SSB proof and for Λ^t_nn(Q,0)=0 (footnote 44).

how reviews work

0 comments

Cite this review

Pith. "Pith review of Kekul\'e Superconductivity in Twisted Magic Angle Bilayer Graphene." pith.science (2026). https://pith.science/paper/XOSHCD72

@misc{pith2026251006451,
  author       = {Pith},
  title        = {Pith review of: Kekul\'e Superconductivity in Twisted Magic Angle Bilayer Graphene},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XOSHCD72}},
  note         = {Machine review of arXiv:2510.06451}
}

abstract

While it has been one of the most important new physics discoveries in the last decade, the nature of superconductivity in the twisted graphene family remains an unsolved problem. Motivated by recent scanning tunneling experiments that report Kekul\'e ordering in moir\'e graphene superconductors, we develop a microscopic theory of this superconductivity for the twisted bilayer system. The pairing we find is an intra-valley, finite-momentum pair-density wave (PDW) that intrinsically carries a Kekul\'e modulation. This state exhibits four salient features: (i) spontaneous breaking of $C_3$ rotation symmetry, producing nematic order (ii)with triplet pairing; and (iii) a quasiparticle density of states that evolves from a V-shaped profile to a fully gapped, U-shaped spectrum as the attraction increases which is accompanied by (iv) systematic behavior of the temperature dependent zero bias conductance. These features align with key experimental signatures. We find, as well, that with only modest interaction strengths, the state is near to a BEC-like phase, consistent with the observed extremely short coherence lengths. Taken together, these results identify a microscopic intra-valley Kekul\'e PDW as a compelling candidate for unconventional superconductivity in the twisted graphene family.

Figures

Figures reproduced from arXiv: 2510.06451 by the authors.

**Figure 2.** FIG. 2. Flat-band dispersion and form factors from the Bistritzer– [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗

**Figure 3.** FIG. 3. The figure shows the grand-canonical thermodynamic [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

Figures from the paper (5 more)

**Figure 4.** Figure 4: FIG. 4. Order-parameter amplitude [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]

**Figure 5.** Figure 5: FIG. 5. (a) Density of states (DOS) versus energy [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]

**Figure 6.** Figure 6: FIG. 6. Density of states and order parameter versus temperature, [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]

**Figure 7.** Figure 7: FIG. 7. U-shaped density of states (DOS) from two-flat-band super [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]

**Figure 8.** Figure 8: FIG. 8. Quasi-particle density of states (reflecting zero bias con [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]

Discussion (0). Sign in to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

d-Wave pair density wave superconductivity in a two-orbital model
cond-mat.str-el 2026-04 unverdicted novelty 6.0 of 10

Two-orbital models with interorbital t-J or t-V interactions exhibit incommensurate d_xy pair density wave superconductivity driven by interband pairing, confirmed by RPA and strong-coupling hard-core pair calculation...

Reference graph

Works this paper leans on

57 extracted references · 4 linked inside Pith · cited by 1 Pith paper

[47]

Van Loon and C

S. Van Loon and C. S´a de Melo, Topological two-band electron- hole superconductors with d-wave symmetry: Absence of dirac quasiparticle annihilation in magic-angle twisted trilayer graphene, Physical Review B111, 064515 (2025)

2025
[1]

Y . Cao, V . Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxi- ras, and P. Jarillo-Herrero, Unconventional superconductivity in magic-angle graphene superlattices, Nature556, 43 (2018)

2018
[2]

Bistritzer and A

R. Bistritzer and A. H. MacDonald, Moir ´e bands in twisted double-layer graphene, Proceedings of the National Academy of Sciences108, 12233 (2011)

2011
[3]

K. P. Nuckolls and A. Yazdani, A microscopic perspective on moir´e materials, Nature Reviews Materials9, 460 (2024)

2024
[4]

Polshyn, M

H. Polshyn, M. Yankowitz, S. Chen, Y . Zhang, K. Watan- abe, T. Taniguchi, C. R. Dean, and A. F. Young, Large linear- in-temperature resistivity in twisted bilayer graphene, Nature Physics15, 1011 (2019)

2019
[5]

Y . Cao, D. Chowdhury, D. Rodan-Legrain, O. Rubies-Bigorda, K. Watanabe, T. Taniguchi, T. Senthil, and P. Jarillo-Herrero, Strange metal in magic-angle graphene with near planckian dis- sipation, Phys. Rev. Lett.124, 076801 (2020)

2020
[6]

Banerjee, Z

A. Banerjee, Z. Hao, M. Kreidel, P. Ledwith, I. Phinney, J. M. Park, A. Zimmerman, M. E. Wesson, K. Watanabe, T. Taniguchi,et al., Superfluid stiffness of twisted trilayer graphene superconductors, Nature638, 93 (2025)

2025
[7]

Tanaka, J

M. Tanaka, J. ˆI.-j. Wang, T. H. Dinh, D. Rodan-Legrain, S. Za- man, M. Hays, A. Almanakly, B. Kannan, D. K. Kim, B. M. Niedzielski,et al., Superfluid stiffness of magic-angle twisted bilayer graphene, Nature638, 99 (2025)

2025

Show all 57 references

[8]

Jiang, X

Y . Jiang, X. Lai, K. Watanabe, T. Taniguchi, K. Haule, J. Mao, and E. Y . Andrei, Charge order and broken rotational symmetry in magic-angle twisted bilayer graphene, Nature573, 91 (2019)

2019
[9]

M. Oh, K. P. Nuckolls, D. Wong, R. L. Lee, X. Liu, K. Watan- abe, T. Taniguchi, and A. Yazdani, Evidence for unconventional superconductivity in twisted bilayer graphene, Nature600, 240 (2021)

2021
[10]

Y . Xie, B. Lian, B. J ¨ack, X. Liu, C.-L. Chiu, K. Watanabe, T. Taniguchi, B. A. Bernevig, and A. Yazdani, Spectroscopic signatures of many-body correlations in magic-angle twisted bi- layer graphene, Nature572, 101 (2019)

2019
[11]

J. M. Park, Y . Cao, K. Watanabe, T. Taniguchi, and P. Jarillo- Herrero, Tunable strongly coupled superconductivity in magic- angle twisted trilayer graphene, Nature590, 249 (2021)

2021
[12]

Zondiner, A

U. Zondiner, A. Rozen, D. Rodan-Legrain, Y . Cao, R. Queiroz, T. Taniguchi, K. Watanabe, Y . Oreg, F. von Oppen, A. Stern, et al., Cascade of phase transitions and dirac revivals in magic- angle graphene, Nature582, 203 (2020)

2020
[13]

Yankowitz, S

M. Yankowitz, S. W. Chen, H. Polshyn, Y . X. Zhang, K. Watan- abe, T. Taniguchi, D. Graf, A. F. Young, and C. R. Dean, Tun- ing superconductivity in twisted bilayer graphene, Science363, 1059 (2019)

2019
[14]

A. L. Sharpe, E. J. Fox, A. W. Barnard, J. Finney, K. Watan- abe, T. Taniguchi, M. A. Kastner, and D. Goldhaber-Gordon, Emergent ferromagnetism near three-quarters filling in twisted bilayer graphene, Science365, 605 (2019)

2019
[15]

Stepanov, I

P. Stepanov, I. Das, X. Lu, A. Fahimniya, K. Watanabe, T. Taniguchi, F. H. Koppens, J. Lischner, L. Levitov, and D. K. Efetov, The interplay of insulating and superconducting orders in magic-angle graphene bilayers, arXiv:1911.09198 (2019)

1911 arXiv
[16]

Saito, J

Y . Saito, J. Ge, K. Watanabe, T. Taniguchi, and A. F. Young, De- coupling superconductivity and correlated insulators in twisted bilayer graphene, arXiv:1911.13302 (2019)

1911 arXiv
[17]

K. P. Nuckolls, R. L. Lee, M. Oh, D. Wong, T. Soejima, J. P. Hong, D. C ˘alug˘aru, J. Herzog-Arbeitman, B. A. Bernevig, K. Watanabe,et al., Quantum textures of the many-body wave- functions in magic-angle graphene, Nature620, 525 (2023)

2023
[18]

H. Kim, Y . Choi, ´E. Lantagne-Hurtubise, C. Lewandowski, A. Thomson, L. Kong, H. Zhou, E. Baum, Y . Zhang, L. Holleis, et al., Imaging inter-valley coherent order in magic-angle twisted trilayer graphene, Nature623, 942 (2023)

2023
[19]

Wagner, Y

G. Wagner, Y . H. Kwan, N. Bultinck, S. H. Simon, and S. Parameswaran, Global phase diagram of the normal state of twisted bilayer graphene, Physical review letters128, 156401 (2022)

2022
[20]

Y . H. Kwan, G. Wagner, T. Soejima, M. P. Zaletel, S. H. Simon, S. A. Parameswaran, and N. Bultinck, Kekul´e spiral order at all nonzero integer fillings in twisted bilayer graphene, Physical Review X11, 041063 (2021)

2021
[21]

Roy and I

B. Roy and I. F. Herbut, Unconventional superconductivity on honeycomb lattice: Theory of kekule order parameter, Phys. Rev. B82, 035429 (2010)

2010
[22]

Tsuchiya, J

S. Tsuchiya, J. Goryo, E. Arahata, and M. Sigrist, Cooperon condensation and intravalley pairing states in honeycomb dirac systems, Physical Review B94, 104508 (2016)

2016
[23]

Barlas, F

Y . Barlas, F. Zhang, and E. Rossi, Quantum geometry in- duced kekul ´e superconductivity in haldane phases (2025), arXiv:2508.21791 [cond-mat.supr-con]

2025
[25]

You and A

Y .-Z. You and A. Vishwanath, Superconductivity from valley fluctuations and approximate so (4) symmetry in a weak cou- pling theory of twisted bilayer graphene, npj Quantum Materi- als4, 16 (2019)

2019
[26]

F. Wu, A. H. MacDonald, and I. Martin, Theory of phonon- mediated superconductivity in twisted bilayer graphene, Physi- cal review letters121, 257001 (2018)

2018
[27]

B. Lian, Z. Wang, and B. A. Bernevig, Twisted bilayer graphene: a phonon-driven superconductor, Physical review letters122, 257002 (2019)

2019
[28]

Cea and F

T. Cea and F. Guinea, Coulomb interaction, phonons, and su- perconductivity in twisted bilayer graphene, Proceedings of the National Academy of Sciences118, e2107874118 (2021)

2021
[29]

X. Hu, T. Hyart, D. I. Pikulin, and E. Rossi, Geometric and conventional contribution to the superfluid weight in twisted bi- layer graphene, Physical Review Letters123, 237002 (2019)

2019
[30]

Julku, T

A. Julku, T. J. Peltonen, L. Liang, T. T. Heikkil¨a, and P. T¨orm¨a, Superfluid weight and berezinskii-kosterlitz-thouless transition temperature of twisted bilayer graphene, Physical Review B 101, 060505 (2020). 10

2020
[31]

Christos, S

M. Christos, S. Sachdev, and M. S. Scheurer, Nodal band-off- diagonal superconductivity in twisted graphene superlattices, Nature Communications14, 7134 (2023)

2023
[32]

Khalaf, S

E. Khalaf, S. Chatterjee, N. Bultinck, M. P. Zaletel, and A. Vishwanath, Charged skyrmions and topological origin of superconductivity in magic-angle graphene, Science advances 7, eabf5299 (2021)

2021
[33]

Gonzalez and T

J. Gonzalez and T. Stauber, Kohn-luttinger superconductivity in twisted bilayer graphene, Physical review letters122, 026801 (2019)

2019
[34]

H. C. Po, L. Zou, A. Vishwanath, and T. Senthil, Origin of mott insulating behavior and superconductivity in twisted bi- layer graphene, Physical Review X8, 031089 (2018)

2018
[35]

D. F. Agterberg, J. S. Davis, S. D. Edkins, E. Fradkin, D. J. Van Harlingen, S. A. Kivelson, P. A. Lee, L. Radzihovsky, J. M. Tranquada, and Y . Wang, The physics of pair-density waves: cuprate superconductors and beyond, Annual Review of Con- densed Matter Physics11, 231 (2020)

2020
[36]

This will then lead to an observed Kekule pattern in the charge channel which can be observed in STM

A coexisting intervalley condensate with zero center-of-mass momentum,∆ 0, can also couple to the PDW order to induce a secondary CDW: The charge modulation from this second mechanism has a momentum of2Qand, like the first, behaves ascos(2K·r)on the atomic scale. This will the...
[37]

Z. Bi, N. F. Q. Yuan, and L. Fu, Designing flat bands by strain, Phys. Rev. B100, 035448 (2019)

2019
[38]

Hejazi, C

K. Hejazi, C. Liu, H. Shapourian, X. Chen, and L. Balents, Mul- tiple topological transitions in twisted bilayer graphene near the first magic angle, Phys. Rev. B99, 035111 (2019)

2019
[39]

In the Bistritzer-MacDonald model,|Gj| can be up to three or four times the magnitude of the primitive moir´e reciprocal lattice vectors

Rigorously, the momentum transfer involved isq+G j, where qlies within the first mini-Brillouin Zone andG j is a moir´e re- ciprocal lattice vector. In the Bistritzer-MacDonald model,|Gj| can be up to three or four times the magnitude of the primitive moir´e reciprocal lattice...
[40]

This conclusion is general and applies to a large class of short- range interaction, as the argument holds if the set of nearest- neighbor vectors,δ j, is replaced by any short-range atomic vec- tor and itsC 3 rotational partners
[41]

This differs from the known FF state, whereQmust be solved self-consistently to satisfy the zero-current condition

The zero-current condition is automatically satisfied once the time-reversal symmetric partner at the other valley is included. This differs from the known FF state, whereQmust be solved self-consistently to satisfy the zero-current condition. Important here also is that a con...
[42]

Y . Cao, J. M. Park, K. Watanabe, T. Taniguchi, and P. Jarillo- Herrero, Pauli-limit violation and re-entrant superconductivity in moir´e graphene, Nature595, 526 (2021)

2021
[43]

Since the Hamiltonian preserves full spin rotation symmetry, we select the triplet channel withS z = 0without loss of gen- erality
[44]

TheC 2Tsymmetry enforces the conditionΛ t nn(Q,q=0) = 0
[45]

H. Kim, Y . Choi, C. Lewandowski, A. Thomson, Y . Zhang, R. Polski, K. Watanabe, T. Taniguchi, J. Alicea, and S. Nadj- Perge, Evidence for unconventional superconductivity in twisted trilayer graphene, Nature606, 494 (2022)

2022
[46]

Q. Chen, Z. Wang, R. Boyack, S. Yang, and K. Levin, When superconductivity crosses over: From bcs to bec, Reviews of Modern Physics96, 025002 (2024)

2024
[48]

Although their Bogoliubov Fermi surfaces are larger, there is some similarity here to work from Christos et al, Nat. Commun. 14, 7134 and Agterberg et al, Phys. Rev. Lett. 118, 127001
[49]

2, we observe that the off-diagonal components of Λare significantly larger than the diagonal ones

From Fig. 2, we observe that the off-diagonal components of Λare significantly larger than the diagonal ones. Consequently, we retain only the off-diagonal pairing
[50]

J. M. Park, S. Sun, K. Watanabe, T. Taniguchi, and P. Jarillo- Herrero, Simultaneous transport and tunneling spectroscopy of moir\’e graphene: Distinct observation of the superconducting gap and signatures of nodal superconductivity, arXiv preprint arXiv:2503.16410 (2025)

2025 arXiv
[51]

Pushp, C

A. Pushp, C. V . Parker, A. N. Pasupathy, K. K. Gomes, S. Ono, J. Wen, Z. Xu, G. Gu, and A. Yazdani, Extend- ing universal nodal excitations optimizes superconductivity in bi2sr2cacu2o8+δ, Science324, 1689 (2009)

2009
[52]

Reber, N

T. Reber, N. Plumb, Y . Cao, Z. Sun, Q. Wang, K. McElroy, H. Iwasawa, M. Arita, J. Wen, Z. Xu,et al., Prepairing and the filling of the gap in the cuprates from the tomographic density of states, Physical Review B87, 060506 (2013)

2013
[53]

Private Communications with Hyunjin Kim
[54]

Y . Cao, D. Rodan-Legrain, J. M. Park, N. F. Q. Yuan, K. Watan- abe, T. Taniguchi, R. M. Fernandes, L. Fu, and P. Jarillo- Herrero, Nematicity and competing orders in superconducting magic-angle graphene, Science372, 264 (2021)

2021
[55]

Z. Wang, K. Wang, and K. Levin, Superconductivity on the edge of vanishing magnetic order (2025), arXiv:2506.18996 [cond-mat.supr-con]

2025 arXiv
[56]

Jian, Y .-F

S.-K. Jian, Y .-F. Jiang, and H. Yao, Emergent spacetime su- persymmetry in 3d weyl semimetals and 2d dirac semimetals, Phys. Rev. Lett.114, 237001 (2015)

2015
[57]

K. Wang, Q. Chen, R. Boyack, and K. Levin, Anomalous su- perfluid density in pair-density-wave superconductors (2025), arXiv:2506.13631 [cond-mat.supr-con]

2025
[58]

The proof is similar

An equivalent approach is to show that the condensate is not invariant underC 3. The proof is similar

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.

[1] [47]

Van Loon and C

S. Van Loon and C. S´a de Melo, Topological two-band electron- hole superconductors with d-wave symmetry: Absence of dirac quasiparticle annihilation in magic-angle twisted trilayer graphene, Physical Review B111, 064515 (2025)

2025

[2] [1]

Y . Cao, V . Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxi- ras, and P. Jarillo-Herrero, Unconventional superconductivity in magic-angle graphene superlattices, Nature556, 43 (2018)

2018

[3] [2]

Bistritzer and A

R. Bistritzer and A. H. MacDonald, Moir ´e bands in twisted double-layer graphene, Proceedings of the National Academy of Sciences108, 12233 (2011)

2011

[4] [3]

K. P. Nuckolls and A. Yazdani, A microscopic perspective on moir´e materials, Nature Reviews Materials9, 460 (2024)

2024

[5] [4]

Polshyn, M

H. Polshyn, M. Yankowitz, S. Chen, Y . Zhang, K. Watan- abe, T. Taniguchi, C. R. Dean, and A. F. Young, Large linear- in-temperature resistivity in twisted bilayer graphene, Nature Physics15, 1011 (2019)

2019

[6] [5]

Y . Cao, D. Chowdhury, D. Rodan-Legrain, O. Rubies-Bigorda, K. Watanabe, T. Taniguchi, T. Senthil, and P. Jarillo-Herrero, Strange metal in magic-angle graphene with near planckian dis- sipation, Phys. Rev. Lett.124, 076801 (2020)

2020

[7] [6]

Banerjee, Z

A. Banerjee, Z. Hao, M. Kreidel, P. Ledwith, I. Phinney, J. M. Park, A. Zimmerman, M. E. Wesson, K. Watanabe, T. Taniguchi,et al., Superfluid stiffness of twisted trilayer graphene superconductors, Nature638, 93 (2025)

2025

[8] [7]

Tanaka, J

M. Tanaka, J. ˆI.-j. Wang, T. H. Dinh, D. Rodan-Legrain, S. Za- man, M. Hays, A. Almanakly, B. Kannan, D. K. Kim, B. M. Niedzielski,et al., Superfluid stiffness of magic-angle twisted bilayer graphene, Nature638, 99 (2025)

2025