Superconductivity from Quasiparticle Pairing of Intervalley Coherent State in Rhombohedral Trilayer Graphene
Pith reviewed 2026-05-24 01:59 UTC · model grok-4.3
The pith
Superconductivity in rhombohedral trilayer graphene arises from pairing quasiparticles of the adjacent intervalley coherent state.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The rhombohedral trilayer graphene superconducting phase arises from the pairing of quasiparticles of the adjacent inter-valley coherent state. Using gapped Dirac cones with chemical potential close to the valence band edge yields Tc proportional to ε_D exp(-2/ρ_qp U) and ξ ∼ v/√(μ Tc) that align with measured values.
What carries the argument
Quasiparticle pairing within gapped Dirac cones of the intervalley coherent state, whose suppressed density of states sets the scale of Tc.
If this is right
- Tc follows Tc ∝ ε_D exp(-2/ρ_qp U) with ρ_qp the suppressed density of states of the intervalley coherent quasiparticles.
- Coherence length follows ξ ∼ v/√(μ Tc) instead of the conventional BCS relation.
- The same microscopic model reproduces both Tc and ξ across the experimental doping window without parameter adjustment.
Where Pith is reading between the lines
- The same quasiparticle-pairing route may operate in other multilayer graphene systems that host intervalley coherent order.
- Direct measurement of the quasiparticle density of states near the band edge would provide an independent test of the required suppression.
Load-bearing premise
The chemical potential sits close to the valence band edge inside the gapped Dirac cones of the intervalley coherent state.
What would settle it
Spectroscopic or transport data showing that the chemical potential lies far from the valence band edge, or that the quasiparticle density of states is not correspondingly suppressed, would remove the mechanism that produces the observed Tc and ξ.
Figures
read the original abstract
Superconductivity is observed in rhombohedral trilayer graphene in a narrow regime between the flavor-symmetric state and the symmetry breaking phase, which cannot be described by the conventional Bardeen-Cooper-Schrieffer theory. The measured coherence length, for instance, is roughly two orders of magnitude shorter than the value predicted by the Bardeen-Cooper-Schrieffer relation based on the large fermi velocity and an extremely low charge carrier density of the flavor-symmetric phase. To resolve the discrepancies, we propose that the rhombohedral trilayer graphene superconducting phase arises from the pairing of quasiparticles of the adjacent inter-valley coherent state. We illustrate the superconducting phenomenology using gapped Dirac cones with the chemical potential $\mu$ close to the valence band's edge. Our findings indicate that the transition temperature $T_c$ obeys $T_c\propto \epsilon_D\exp(-2/\rho_\mathrm{qp}U)$ with the density of states $\rho_\mathrm{qp}$ of intervalley coherent state quasiparticles, which is much suppressed compared to predictions from the Bardeen-Cooper-Schrieffer theory. The coherence length $\xi$ we predict behaves according to $\xi\sim v/\sqrt{\mu T_c}$ with $v$ being the velocity of Dirac cone. Applying our assumption to a microscopic model, our predictions align well with experimental data and effectively capture key measurable quantities such as the transition temperature $T_c$ and the coherence length $\xi$ without parameter fine-tuning.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that superconductivity in rhombohedral trilayer graphene arises from quasiparticle pairing within the adjacent intervalley coherent (IVC) state rather than the flavor-symmetric metallic phase. Modeling the IVC state with gapped Dirac cones and placing the chemical potential μ near the valence-band edge suppresses the quasiparticle density of states ρ_qp, yielding Tc ∝ ε_D exp(−2/ρ_qp U) and ξ ∼ v/√(μ Tc). The authors state that applying this assumption to a microscopic model reproduces experimental Tc and ξ values without parameter fine-tuning.
Significance. If the central assumption is microscopically justified, the proposal supplies a concrete mechanism that reconciles the two-order-of-magnitude mismatch between the observed short coherence length and conventional BCS estimates based on the large Fermi velocity and low carrier density of the flavor-symmetric state. The explicit, falsifiable expressions for Tc and ξ constitute a strength that can be tested against further doping- or field-dependent measurements.
major comments (2)
- [Abstract] Abstract (paragraph describing the model): the claim that μ lies close to the valence-band edge is introduced as an assumption required to suppress ρ_qp sufficiently; no self-consistent determination of μ from the total carrier density of the IVC state or from the microscopic Hamiltonian is shown. If the solved microscopic model instead places μ near mid-gap, ρ_qp is no longer parametrically small and both the Tc formula and the ξ scaling revert to conventional BCS expectations, removing the resolution of the coherence-length discrepancy.
- [Abstract] Abstract: the statement that predictions “align well with experimental data … without parameter fine-tuning” is not accompanied by an explicit fitting procedure, error analysis, or comparison table showing how the microscopic IVC solution determines ρ_qp, ε_D, and U independently of the measured Tc and ξ.
minor comments (1)
- Notation for the Dirac velocity v and the cutoff ε_D should be defined once at first use and used consistently.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments highlight important points about clarity in the abstract and the need for explicit documentation of the microscopic inputs. We address each below and will revise the manuscript to improve transparency while preserving the central claim that the IVC quasiparticle pairing mechanism accounts for the observed Tc and short coherence length.
read point-by-point responses
-
Referee: [Abstract] Abstract (paragraph describing the model): the claim that μ lies close to the valence-band edge is introduced as an assumption required to suppress ρ_qp sufficiently; no self-consistent determination of μ from the total carrier density of the IVC state or from the microscopic Hamiltonian is shown. If the solved microscopic model instead places μ near mid-gap, ρ_qp is no longer parametrically small and both the Tc formula and the ξ scaling revert to conventional BCS expectations, removing the resolution of the coherence-length discrepancy.
Authors: The abstract presents the μ placement as an illustrative assumption to derive the suppressed ρ_qp and resulting Tc, ξ scalings. In the body of the manuscript we apply the same framework to a microscopic model of the IVC state (with gap parameters and total filling fixed by the Hamiltonian), which self-consistently places μ near the valence-band edge for the relevant carrier densities. We will revise the abstract to state explicitly that the microscopic solution determines μ near the edge rather than mid-gap, thereby keeping ρ_qp parametrically small. Should the microscopic solution instead yield mid-gap μ, the mechanism would indeed revert to conventional BCS scaling and fail to resolve the coherence-length issue; our calculations indicate this does not occur. revision: yes
-
Referee: [Abstract] Abstract: the statement that predictions “align well with experimental data … without parameter fine-tuning” is not accompanied by an explicit fitting procedure, error analysis, or comparison table showing how the microscopic IVC solution determines ρ_qp, ε_D, and U independently of the measured Tc and ξ.
Authors: The quoted statement refers to using gap size ε_D, velocity v, and interaction strength U obtained directly from the microscopic IVC calculation (independent of measured Tc or ξ) to evaluate the analytic expressions for Tc and ξ. We will add a table (or supplementary section) that lists the microscopic inputs (ρ_qp, ε_D, U), the resulting predicted Tc and ξ, and direct numerical comparison to the experimental values, together with a brief discussion of uncertainties arising from the model parameters. revision: yes
Circularity Check
No significant circularity; derivation rests on explicit modeling assumption applied to standard BCS formulas
full rationale
The paper's central derivation begins with an explicitly stated modeling choice: gapped Dirac cones for the intervalley coherent state with chemical potential placed near the valence band edge to obtain a suppressed quasiparticle DOS. From this, the expressions Tc ∝ ε_D exp(-2/ρ_qp U) and ξ ∼ v/√(μ Tc) follow directly as adaptations of conventional BCS relations to the resulting ρ_qp; these are not shown to be equivalent to the inputs by construction, nor do they rely on self-citations, fitted parameters renamed as predictions, or uniqueness theorems. The claim that the assumption applied to a microscopic model reproduces experimental Tc and ξ is presented as an outcome check rather than a tautology. No load-bearing step reduces to a self-referential definition or citation chain.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption BCS pairing theory structure remains valid when applied to quasiparticles of the intervalley coherent state
invented entities (1)
-
quasiparticle pairing mechanism from intervalley coherent state
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
TMF ∝ ϵD exp(−2/ρU′) … ξcon ≈ v/(4√(μ TMF)) … quantum metric contribution
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
gapped Dirac cones with chemical potential close to the valence band’s edge
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 1 Pith paper
-
Pair density wave in quarter metals from a repulsive fermionic interaction in graphene heterostructures: A renormalization group study
Leading-order RG analysis shows repulsive interactions stabilize a chiral odd-parity pair density wave in quarter metals of chirally stacked graphene heterostructures.
Reference graph
Works this paper leans on
- [1]
-
[2]
Case 2: TMF ≫ m 23 E. Discussion on pairing and time-reversal symmetry in IVC 24 I. INTRODUCTION Graphene heterostructures, in particular twisted bi- layer graphene, have been observed experimentally to host various strongly correlated phases [1–11]. How- ever, many of these heterostructures are difficult to re- alize in experiments due to structural inst...
work page internal anchor Pith review Pith/arXiv arXiv 2024
-
[3]
Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, R. C. Ashoori, and P. Jarillo- Herrero, Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature (London) 556, 80 (2018), arXiv:1802.00553 [cond-mat.mes-hall]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [4]
- [5]
-
[6]
H. S. Arora, R. Polski, Y. Zhang, A. Thomson, Y. Choi, H. Kim, Z. Lin, I. Z. Wilson, X. Xu, J.-H. Chu, K. Watan- abe, T. Taniguchi, J. Alicea, and S. Nadj-Perge, Super- conductivity in metallic twisted bilayer graphene stabi- lized by WSe 2, Nature (London) 583, 379 (2020)
work page 2020
-
[7]
Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, Unconventional super- conductivity in magic-angle graphene superlattices, Na- ture (London) 556, 43 (2018), arXiv:1803.02342 [cond- mat.mes-hall]. 9
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[8]
X. Lu, P. Stepanov, W. Yang, M. Xie, M. A. Aamir, I. Das, C. Urgell, K. Watanabe, T. Taniguchi, G. Zhang, A. Bachtold, A. H. MacDonald, and D. K. Efetov, Su- perconductors, orbital magnets and correlated states in magic-angle bilayer graphene, Nature (London) 574, 653 (2019), arXiv:1903.06513 [cond-mat.str-el]
-
[10]
Tuning superconductivity in twisted bilayer graphene
M. Yankowitz, S. Chen, H. Polshyn, Y. Zhang, K. Watan- abe, T. Taniguchi, D. Graf, A. F. Young, and C. R. Dean, Tuning superconductivity in twisted bilayer graphene, Science 363, 1059 (2019), arXiv:1808.07865 [cond- mat.mes-hall]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[11]
A. Kerelsky, L. J. McGilly, D. M. Kennes, L. Xian, M. Yankowitz, S. Chen, K. Watanabe, T. Taniguchi, J. Hone, C. Dean, A. Rubio, and A. N. Pasupathy, Maxi- mized electron interactions at the magic angle in twisted bilayer graphene, Nature (London) 572, 95 (2019)
work page 2019
- [12]
- [13]
-
[14]
X. Lin, D. Liu, and D. Tom´ anek, Shear instability in twisted bilayer graphene, Phys. Rev. B 98, 195432 (2018)
work page 2018
- [15]
- [16]
- [17]
-
[18]
T. Arp, O. Sheekey, H. Zhou, C. L. Tschirhart, C. L. Patterson, H. M. Yoo, L. Holleis, E. Redekop, G. Babikyan, T. Xie, J. Xiao, Y. Vituri, T. Holder, T. Taniguchi, K. Watanabe, M. E. Huber, E. Berg, and A. F. Young, Intervalley coherence and intrinsic spin- orbit coupling in rhombohedral trilayer graphene, arXiv e-prints , arXiv:2310.03781 (2023), arXiv:...
- [19]
-
[20]
Trigonal warping and Berry's phase N pi in ABC-stacked multilayer graphene
M. Koshino and E. McCann, Trigonal warping and Berry’s phase N π in ABC-stacked multilayer graphene, Phys. Rev. B 80, 165409 (2009), arXiv:0906.4634 [cond- mat.mes-hall]
work page internal anchor Pith review Pith/arXiv arXiv 2009
- [21]
-
[22]
Y. Zhumagulov, D. Kochan, and J. Fabian, Emer- gent correlated phases in rhombohedral trilayer graphene induced by proximity spin-orbit and ex- change coupling, arXiv e-prints , arXiv:2305.14277 (2023), arXiv:2305.14277 [cond-mat.str-el]
-
[23]
J. M. Koh, J. Alicea, and E. Lantagne-Hurtubise, Corre- lated phases in spin-orbit-coupled rhombohedral trilayer graphene, Phys. Rev. B 109, 035113 (2024)
work page 2024
- [24]
-
[25]
C. H. Lui, Z. Li, K. F. Mak, E. Cappelluti, and T. F. Heinz, Observation of an electrically tunable band gap in trilayer graphene, Nature Physics 7, 944 (2011), arXiv:1105.4658 [cond-mat.mes-hall]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[26]
Y.-Z. Chou, F. Wu, J. D. Sau, and S. Das Sarma, Acoustic-phonon-mediated superconductivity in moir´ eless graphene multilayers, Phys. Rev. B 106, 024507 (2022)
work page 2022
-
[27]
E. Vi˜ nas Bostr¨ om, A. Fischer, J. B. Hauck, J. Zhang, D. M. Kennes, and A. Rubio, Phonon-mediated un- conventional s- and f-wave pairing superconductivity in rhombohedral stacked multilayer graphene, arXiv e-prints , arXiv:2311.02494 (2023), arXiv:2311.02494 [cond-mat.supr-con]
- [28]
-
[29]
A. Ghazaryan, T. Holder, M. Serbyn, and E. Berg, Un- conventional Superconductivity in Systems with Annu- lar Fermi Surfaces: Application to Rhombohedral Tri- layer Graphene, Phys. Rev. Lett. 127, 247001 (2021), arXiv:2109.00011 [cond-mat.supr-con]
- [30]
-
[31]
A. Jimeno-Pozo, H. Sainz-Cruz, T. Cea, P. A. Pan- tale´ on, and F. Guinea, Superconductivity from electronic interactions and spin-orbit enhancement in bilayer and trilayer graphene, Phys. Rev. B 107, L161106 (2023), arXiv:2210.02915 [cond-mat.mes-hall]
- [32]
- [33]
- [34]
-
[35]
Z. Dong and L. Levitov, Superconductivity in the vicinity of an isospin-polarized state in a cubic Dirac band, arXiv e-prints , arXiv:2109.01133 (2021), arXiv:2109.01133 10 [cond-mat.supr-con]
-
[36]
S. Chatterjee, T. Wang, E. Berg, and M. P. Za- letel, Inter-valley coherent order and isospin fluctua- tion mediated superconductivity in rhombohedral tri- layer graphene, Nature Communications13, 6013 (2022), arXiv:2109.00002 [cond-mat.supr-con]
- [37]
-
[38]
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of Superconductivity, Physical Review 108, 1175 (1957)
work page 1957
-
[39]
See Supplemental Material
-
[40]
Electric Field Control of Soliton Motion and Stacking in Trilayer Graphene
M. Yankowitz, J. I. J. Wang, A. G. Birdwell, Y.-A. Chen, K. Watanabe, T. Taniguchi, P. Jacquod, P. San-Jose, P. Jarillo-Herrero, and B. J. Leroy, Electric field con- trol of soliton motion and stacking in trilayer graphene, Nature Materials 13, 786 (2014), arXiv:1401.7663 [cond- mat.mtrl-sci]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[41]
S. A. Chen and K. T. Law, Ginzburg-landau theory of flat-band superconductors with quantum metric, Phys. Rev. Lett. 132, 026002 (2024)
work page 2024
- [42]
-
[43]
X. Han, Q. Liu, Y. Wang, R. Niu, Z. Qu, Z. Wang, Z. Li, C. Han, K. Watanabe, T. Taniguchi, Z. Song, J. Mao, Z. V. Han, Z. Gan, and J. Lu, Chemical Potential Char- acterization of Symmetry-Breaking Phases in a Rhombo- hedral Trilayer Graphene, Nano Letters 23, 6875 (2023)
work page 2023
-
[44]
F. Winterer, F. R. Geisenhof, N. Fernandez, A. M. Seiler, F. Zhang, and R. T. Weitz, Ferroelectric and anomalous quantum Hall states in bare rhombohedral tri- layer graphene, arXiv e-prints , arXiv:2305.04950 (2023), arXiv:2305.04950 [cond-mat.mes-hall]
-
[45]
Stacking-Dependent Band Gap and Quantum Transport in Trilayer Graphene
W. Bao, L. Jing, J. Velasco, Y. Lee, G. Liu, D. Tran, B. Standley, M. Aykol, S. B. Cronin, D. Smirnov, M. Koshino, E. McCann, M. Bockrath, and C. N. Lau, Stacking-dependent band gap and quantum trans- port in trilayer graphene, Nature Physics 7, 948 (2011), arXiv:1103.6088 [cond-mat.mes-hall]. 11 Appendix for “Origin of Superconductivity in Rhombohedral T...
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[46]
With all information, we can solve for δF . Due to stability of the mean-field solution, a linear correction has to vanish, and the lowest order term is thus second order in δ∆(k), namely F2 the Gaussian fluctuation of the free energy, given by: F2 = Z d2k (2π)2 U |δ∆(k)|2 − T 2 * Z β 0 dτ Z d2k (2π)2 δL !2+ = Z d2k (2π)2 U |δ∆(k)|2 − T 2 X n A Z d2q (2π)...
-
[47]
In particular, we will set it to µ = m + κϵD, and with 0 ≤ κ ≤ 1
Case 1: m ≫ TMF In continuous model, we can treat the chemical potential as a free parameter. In particular, we will set it to µ = m + κϵD, and with 0 ≤ κ ≤ 1. To begin with we should write down our integral over momentum space in terms 18 of energy: Z d2q (2π)2 = Z m+(1+ν)ϵD m ϵ0dϵ0 2πv2 = Z ϵD −κϵD (ε + µ)dε 2πv2 . (D9) Considering the contribution from...
-
[48]
Case 2: TMF ≫ m In this scenario, most of the calculation is identical to the previous case, with the exception of χqm when µ ≈ m: 1 4 k2A Z ϵD 0 dε 2π 1 2 tanh β 2 ε 2ε (ε + µ) = 1 32π k2A Z ϵD 0 dε 1 µ 1 ε − 1 ε + µ tanh β 2 ε = 1 32π k2A Z ϵD 0 dε 1 µ 1 ε tanh β 2 ε − 1 16π k2A Z ϵD 0 dε 1 µ 1 ε + µ tanh β 2 ε ≈ 1 32π k2A 1 m Z ϵD 0 dε 1 ε tanh β 2 ε −...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.