Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Superconductivity from collective excitations in magic angle twisted bilayer graphene

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read Plasmons can glue Cooper pairs in magic-angle twisted bilayer graphene.

desk verdict A timely new application of plasmon-mediated pairing to magic-angle tBG, with a sound qualitative mechanism but quantitative claims that need an honest referee to tighten. read the letter →

arxiv 1909.02574 v1 pith:6ASDP4KP submitted 2019-09-05 cond-mat.mes-hall cond-mat.supr-con

Gargee Sharma , Maxim Trushin , Oleg P. Sushkov , Giovanni Vignale , Shaffique Adam This is my paper · ORCID
classification cond-mat.mes-hallcond-mat.supr-con
keywords superconductivitytwistedbilayergraphenemagicangleplasmonscollectiveexcitationsMigdal-EliashbergtheoryvanHovesingularityscreenedCoulombinteraction
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

This paper argues that the superconductivity in magic-angle twisted bilayer graphene can arise from a purely electronic mechanism: collective oscillations of the electron fluid, such as plasmons, bind electrons into Cooper pairs, with no phonons involved. The authors solve a standard strong-coupling gap equation on a deliberately simple one-parameter lattice model and find an asymmetrical, dome-shaped dependence of the critical temperature on carrier density, centered near n = $10^{12}$ $cm^{-2}$, with Tc of order a few kelvin, matching experiments. If correct, this identifies the pairing glue in a strongly correlated system and explains why superconductivity shows up in twisted bilayer graphene but not in monolayer graphene.

What carries the argument

The machinery is the Migdal-Eliashberg gap equation for a Coulomb system with no attractive static interaction, together with a single-mode reduction of the pairing kernel. The dynamically screened interaction $V(\mathbf{q}, i\omega)=V(\mathbf{q})/\varepsilon(\mathbf{q}, i\omega)$ is momentum-averaged to give the dimensionless coupling $\lambda(i\omega)=N(E_F)\langle\langle V(i\omega)\rangle\rangle$, fitted to the Lorentzian form $\lambda_{nm}=\mu(1-\omega_b^2/[ (\omega_n-\omega_m)^2+\omega_b^2])$. Here $\mu$ is the high-frequency Coulomb pseudopotential and $\omega_b$ is the collective-mode frequency; $\lambda$ is near zero at low frequency and rises at high frequency, opposite to a phonon kernel, and that frequency dependence generates attraction without phonons. The band structure enters through a one-parameter nearest-neighbor hexagonal tight-binding model with $t_{\mathrm{eff}}=W/3$, which reproduces the tBG flat bands and the asymmetric van Hove singularity (the divergence in the density of states at a saddle point) that controls the dome shape. The gap equation has a solution only when the gap changes sign as a function of Matsubara frequency, and the paper solves it directly or with a pseudopotential method when $T_C\ll E_F$.

What would settle it

A measurement that would decide the issue: tune the twist angle and dielectric screening so that the single-mode criterion $\omega_b/(r_s^2 E_F) < 1$ is violated (e.g., lower density or smaller $r_s$); the theory predicts no $T_C$ solution, so observed superconductivity there would falsify it. Alternatively, a direct measurement showing no high-frequency attraction in the effective pairing interaction would contradict the mechanism.

Watch

Extended reading notes

Core claim

The paper's central claim is that purely electronic collective modes can mediate superconductivity in magic-angle twisted bilayer graphene. The mechanism relies on the frequency dependence of the dynamically screened Coulomb interaction: the dielectric function $\epsilon(\mathbf{q}, i\omega)=1+e^2 E_F q/(2\kappa\omega^2)$ weakens screening at high frequency, making the effective interaction attractive in that dynamic range despite a repulsive static Coulomb force. Using a one-parameter nearest-neighbor tight-binding model with hopping $t_{\mathrm{eff}}=W/3$ fixed by the flat-band bandwidth from the continuum model, the authors compute the polarization, build the screened interaction, and reduce the Migdal-Eliashberg gap equation to a momentum-independent kernel. Solving the gap equation gives an asymmetrical superconducting dome around $n \approx 10^{12}\,\mathrm{cm}^{-2}$ and $T_C=\mathcal{O}(10\,\mathrm{K})$, in agreement with experiments. Because the static interaction is repulsive, the gap function must change sign as a function of Matsubara frequency; the collective modes provide the dynamic attraction that makes this possible.

Load-bearing premise

The load-bearing premise is that the one-parameter effective lattice model with $t_\mathrm{eff}=W/3$ faithfully represents the low-energy flat bands of real magic-angle twisted bilayer graphene, especially the position and asymmetry of the van Hove singularity that sets the shape of the predicted superconducting dome.

Editorial extensions

If this is right

  • If the mechanism is right, no phonon modes are needed: the superconducting state in magic-angle twisted bilayer graphene is driven by electronic collective modes, and the calculated $T_C$ from plasmons exceeds the phonon-mediated value.
  • The superconducting dome in carrier density is asymmetrical, with its main peak when the Fermi energy crosses the van Hove point and a secondary maximum at lower twist angles; density-scan experiments should see this profile.
  • Superconductivity appears only above a threshold combination of density and coupling strength $\omega_b/(r_s^2 E_F) < 1$, so the state should be absent at very low electron concentrations.
  • As the twist angle moves away from the magic angle toward the monolayer limit, the dome narrows and disappears, consistent with the absence of superconductivity in monolayer graphene.
  • The gap function must change sign as a function of Matsubara frequency; a purely static repulsive interaction cannot produce pairing, so the frequency-dependent attraction is an observable signature.

Reading between the lines

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

  • We infer that the same effective-lattice-plus-screening construction should apply to other moiré flat-band systems: wherever a collective density mode sits within a narrow, van Hove-structured band, a dome of electronic superconductivity could appear.
  • We infer that substrate engineering of the dielectric constant should move $T_C$ through $r_s$; increasing screening should weaken pairing and eventually destroy the dome, a quantitative prediction the paper leaves implicit.
  • We infer that if the collective mode is strongly damped away from the magic angle, the dome should narrow and shift rather than simply disappear, so dome-width measurements combined with loss-function data could test the mechanism.
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.

Referee Report

3 major / 4 minor

Summary. The paper proposes a purely electronic pairing mechanism for superconductivity in magic-angle twisted bilayer graphene (tBG), mediated by collective density-fluctuation modes (plasmons) within the Migdal-Eliashberg framework. Starting from a one-parameter nearest-neighbor tight-binding model that mimics the flat bands, the authors compute the RPA dynamical polarization and screened Coulomb interaction, then project the frequency-dependent coupling onto a Lorentzian form to solve the gap equation. They report an asymmetric superconducting dome in carrier density around n ~ 10^12 cm^-2 with T_C of order a few kelvin, in agreement with experimental observations, and a non-monotonic twist-angle dependence. The central claims are that the collective electronic modes provide a stronger pairing glue than phonons and that this mechanism is tBG-specific.

Significance. If the result holds, the paper would establish a qualitatively new pairing mechanism for magic-angle tBG, with falsifiable predictions (the density dome and its asymmetry) that do not use experimental superconducting data as input. The numerical treatment goes beyond the single-mode approximation by solving the full momentum-averaged Eliashberg equations with the RPA-screened interaction, which is a strength. However, the central quantitative prediction rests on the Migdal approximation in a strongly coupled, high-frequency regime where its validity is not established; the paper's justification for neglecting vertex corrections is cursory and appears to invert the standard Migdal condition. The fitted interaction parameters are not reported, limiting reproducibility.

major comments (3)
  1. [Concluding remarks (p. 4, bottom)] The justification for neglecting vertex corrections is inadequate and appears to invert the standard Migdal expansion. With ω_b ~ 15 meV and E_F ~ 1.5 meV, the boson frequency is an order of magnitude larger than the Fermi energy, so the usual Migdal parameter (ω_b/E_F) is ~10, not small. In a strongly coupled system with r_s ~ 12, the relevant expansion parameter is effectively λ ω_b/E_F, which is of order 100, not 0.1. The paper's statement that vertex corrections are 'insignificant for processes much larger than EF' is the opposite of the standard Migdal condition, and the citation to Takada [30] does not by itself establish the required cancellation for this model. Since the quoted T_C ~ 2.6 K and the entire dome in Figs. 2-3 are computed from the truncated equations (6)-(8) with vertex corrections set to zero, the central numerical claim is not yet supported without a quantitative estimate of the leading vertex diagram.
  2. [Dynamical Coulomb interaction (p. 4) and Fig. 3] The parameters μ and ω_b are extracted by fitting the momentum-averaged RPA coupling to the Lorentzian form of Eq. (10), but the fitted values are never reported for any twist angle or carrier density. The comparison between the single-mode model of Eq. (2) and the full numerical solution in Fig. 3 cannot be assessed, and the results are not reproducible. Please report the fitted μ and ω_b for the three densities shown in Fig. 3 and for the range of twist angles, along with a measure of the quality of the fit to the numerically computed λ(iω).
  3. [Dynamical Coulomb interaction (pp. 2-3) and Figs. 2-3] The one-parameter tight-binding model of Eq. (4) is the only source of the density of states that determines the dome shape, yet no quantitative validation against the actual tBG band structure is provided. The paper asserts that the model reproduces the symmetry and the van Hove singularity of tBG, but the position of the vHs in energy and the asymmetry of the DOS around the M point are not compared with the continuum model of Ref. [4] or with the band structure used to obtain W. Because the predicted asymmetric dome (Fig. 2) and its twist-angle dependence (Fig. 3) are direct consequences of this DOS, the model's quantitative reliability needs to be demonstrated.
minor comments (4)
  1. [Abstract and Introduction] The abstract calls this a 'one parameter' model, but the calculation uses several additional inputs (μ, ω_b, k_c, M, v_F, r_s). Please clarify which parameter is the single free one or rephrase the description.
  2. [Equation (10) and following text] The inline formula for Z_n has an unbalanced square bracket: 'Z_n = 1 + μ(ω_b/ω_n) arctan{ω_n E_F / [(ω_n^2 + ω_b(E_F+ω_b)]}' should read '... / (ω_n^2 + ω_b(E_F+ω_b))}'. Please fix.
  3. [Figure 1c caption] The dotted lines indicating the threshold density are not defined or labeled on the figure, and the axis labels are unclear. Please specify what is plotted and what the dotted lines denote.
  4. [Equation (2) and surrounding text] The text states M ≫ 1 but then identifies M ~ ω_b/E_F, which is about 10 for the quoted parameters. Please explain why this choice of M is adequate and how the results depend on M.

Circularity Check

0 steps flagged · score 2.0 of 10

No substantive circularity: the fitted parameters mu and omega_b enter the model's own interaction, while T_C and the dome shape are genuine outputs of the Eliashberg gap equation; self-citations are non-load-bearing.

full rationale

The derivation chain is: Eq. (4) one-parameter moiré tight-binding model (teff = W/3, W from the continuum model) -> RPA polarization and dielectric function -> momentum-averaged screened Coulomb interaction -> Lorentzian fit of lambda_nm with parameters mu and omega_b -> Eliashberg gap equation -> T_C and dome shape (Figs. 2-3). The parameters mu and omega_b are obtained by an internal fit to the model's own RPA-screened interaction ('The parameters mu and omega_b are then extracted by fitting the actual coupling lambda(iomega) to Eq. (10)'), and no experimental superconductivity datum is used to set them. T_C is therefore not an input renamed as a prediction; the agreement with experiments is an external benchmark rather than a fitted constraint. The dome shape is explicitly traced to the model's density of states ('the shape of the domes is set by the density of states of the non-interacting bands'), so it is an output of the stated model, not a re-statement of the experimental dome. The self-citations (Refs. 17 and 21) occur only in background literature surveys and are not load-bearing. The manuscript's own limitation notes ('vertex corrections are neglected... may have a quantitative impact on the calculated TC', in the Concluding Remarks, and 'Other competing states... may obviously affect the shape of the dome', in the Superconductivity section) identify robustness risks, not circularity: they concern the accuracy of the truncated Migdal-Eliashberg approximation and the omission of competing orders, not the re-injection of the target result as an input. Hence no circular step meeting the quoted-reduction standard is present; at most one or two non-load-bearing self-citations justify a low nonzero score.

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

The central claim rests on the fidelity of the effective lattice model and on the RPA plus Migdal-Eliashberg framework. The only numerical parameters fit inside the paper are mu and omega_b, which are fit to the model's own screened interaction rather than to the experimental T_C. The stated v_F contains an internal inconsistency that changes the existence of a superconducting solution by two orders of magnitude. No new entities are introduced; the collective mode is a known plasmon or density-fluctuation excitation.

free parameters (6)
  • mu (high-frequency dimensionless coupling)
    Fit to the numerical momentum-averaged coupling lambda(i*omega) via Eq. (10); sets the strength of the repulsive kernel and therefore T_C.
  • omega_b (collective-mode frequency) = ~15 meV in single-mode estimate
    Fit to lambda(i*omega) in the full calculation; input in the single-mode model. T_C depends sensitively on omega_b and E_F through Eq. (3).
  • k_c (momentum-averaging cutoff) = 2 k_F
    Chosen following Grabowski-Sham; the averaging window in Eq. (9) affects the fitted mu and omega_b, yet no sensitivity check is given.
  • M (integer Matsubara cutoff in single-mode model) = ~10
    Set by M ~ omega_b/E_F; used in Eq. (2) to reduce the infinite set of gap equations, no convergence study shown.
  • renormalized Fermi velocity v_F = 1.5e4 m/s as stated
    Used to convert density to E_F; with the stated value E_F = 0.015 meV, not 1.5 meV, so the quoted T_C is not reproduced.
  • r_s (strong-coupling parameter) = ~12 (with kappa ~ 12)
    Derived from kappa and v_F following ref. [25]; controls the existence threshold in Eq. (3), but its value is an estimate from the model, not a fit to T_C.
assumptions (6)
  • domain assumption The one-parameter nearest-neighbor tight-binding model (Eq. 4) with t_eff = W/3 faithfully represents the low-energy flat bands, bandwidth, symmetry, and van Hove singularity of magic-angle twisted bilayer graphene.
    Invoked in 'Dynamical Coulomb interaction'; if the model density of states is wrong, the predicted dome shape and T_C density dependence in Figs. 2-3 are not valid.
  • domain assumption The polarization function expanded to second order in q is sufficient to capture collective modes over 0 < omega < 2W.
    Stated after Eq. (5) citing ref. [25]; the RPA dielectric function and fitted omega_b depend on this truncation.
  • domain assumption The RPA dielectric function epsilon(q,i*omega) = 1 - V_q Pi(q,i*omega) is valid at strong coupling r_s ~ 12.
    Used to define the screened Coulomb interaction; no beyond-RPA checks are presented.
  • domain assumption Migdal-Eliashberg theory with vertex corrections neglected is applicable because the boson frequency omega_b is much larger than E_F.
    Acknowledged in Concluding remarks and justified by ref. [30]; this is not obviously valid for a flat band with r_s ~ 12.
  • domain assumption The momentum-averaged interaction can be mapped to the two-parameter Lorentzian form of Eq. (10) with fitted mu and omega_b.
    This mapping supports the analytic gap equation; no fit quality metrics are reported.
  • domain assumption A single collective-mode frequency scale omega_b characterizes the relevant density-fluctuation spectrum.
    The single-mode approximation and Eq. (3) assume a sharp mode; the text says undamped plasmons are not strictly required but the formalism still relies on one boson scale.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Superconductivity from collective excitations in magic angle twisted bilayer graphene." pith.science (2026). https://pith.science/paper/6ASDP4KP

@misc{pith2026190902574,
  author       = {Pith},
  title        = {Pith review of: Superconductivity from collective excitations in magic angle twisted bilayer graphene},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6ASDP4KP}},
  note         = {Machine review of arXiv:1909.02574}
}
read the original abstract

A purely electronic mechanism is proposed for the unconventional superconductivity recently observed in twisted bilayer graphene (tBG) close to the magic angle. Using the Migdal-Eliashberg framework on a one parameter effective lattice model for tBG we show that a superconducting state can be achieved by means of collective electronic modes in tBG. We posit robust features of the theory, including an asymmetrical superconducting dome and the magnitude of the critical temperature that are in agreement with experiments.

Figures

Figures reproduced from arXiv: 1909.02574 by the authors.

Figure 1
Figure 1. The single-mode model for the dimensionless pair [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Critical temperature TC as a function of electron density in tBG (close to the magic angle). The presence of an asymmetrical superconducting dome around n = 1012cm−2 and TC = O(10K) are the main predictions of our theory [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Critical temperature TC as a function of twist angle in tBG for three different carrier densities. The solid lines show the TC obtained from Eq. 2 using the parameters ωb, EF , and µ(rs) from the numerical solution. The inset shows the normalized density of states (for any angle) at the Fermi energy as the Fermi surface moves up in energy and intersects various symmetry points in the Brillouin zone K → M → Γ. It has… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Plasmonic Cooper pairing in single layer graphene

    cond-mat.supr-con 2019-09 conditional novelty 6.0 of 10

    The authors derive a dielectric-function-method gap equation for a Dirac cone and compute plasmon-mediated superconducting critical temperatures in graphene that rise to the millikelvin range with density.

Reference graph

Works this paper leans on

30 extracted references · 23 canonical work pages · cited by 1 Pith paper

  1. [30]

    Takada, Journal of the Physical Society of Japan 61, 3849 (1992)

    Y. Takada, Journal of the Physical Society of Japan 61, 3849 (1992)

    work page 1992

  2. [4]

    Bistritzer and A

    R. Bistritzer and A. H. MacDonald, Proceedings of the National Academy of Sciences 108, 12233 (2011)

    2011

  3. [1]

    Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, Nature 556, 43 (2018)

    2018

  4. [2]

    Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, et al. , Nature 556, 80 (2018)

    work page 2018

  5. [3]

    Yankowitz, S

    M. Yankowitz, S. Chen, H. Polshyn, Y. Zhang, K. Watan- abe, T. Taniguchi, D. Graf, A. F. Young, and C. R. Dean, Science , eaav1910 (2019)

    work page 2019

  6. [5]

    Xu and L

    C. Xu and L. Balents, Phys. Rev. Lett. 121, 087001 (2018)

    work page 2018

  7. [6]

    Unconventional superconductivity in nearly flat bands in twisted bilayer graphene

    B. Roy and V. Juricic, arXiv preprint arXiv:1803.11190 (2018)

    work page Pith review arXiv 2018

  8. [7]

    H. C. Po, L. Zou, A. Vishwanath, and T. Senthil, Phys. Rev. X 8, 031089 (2018)

    2018

Show all 30 references
  1. [8]

    Koshino, N

    M. Koshino, N. F. Q. Yuan, T. Koretsune, M. Ochi, K. Kuroki, and L. Fu, Phys. Rev. X 8, 031087 (2018)

    2018

  2. [9]

    Kang and O

    J. Kang and O. Vafek, Phys. Rev. X 8, 031088 (2018)

    2018

  3. [10]

    Padhi, C

    B. Padhi, C. Setty, and P. W. Phillips, Nano Lett. 18, 6175 (2018)

    2018

  4. [11]

    H. Guo, X. Zhu, S. Feng, and R. T. Scalettar, Phys. Rev. B 97, 235453 (2018)

    2018

  5. [12]

    Liu, L.-D

    C.-C. Liu, L.-D. Zhang, W.-Q. Chen, and F. Yang, Phys. Rev. Lett. 121, 217001 (2018)

    2018

  6. [13]

    Isobe, N

    H. Isobe, N. F. Q. Yuan, and L. Fu, Phys. Rev. X 8, 041041 (2018)

    2018

  7. [14]

    You and A

    Y.-Z. You and A. Vishwanath, arXiv preprint arXiv:1805.06867 (2018)

    2018 arXiv

  8. [15]

    Gonz´ alez and T

    J. Gonz´ alez and T. Stauber, Phys. Rev. Lett.122, 026801 (2019)

    2019

  9. [16]

    Xie and A

    M. Xie and A. H. MacDonald, arXiv preprint arXiv:1812.04213 (2018)

    2018 arXiv

  10. [17]

    Laksono, J

    E. Laksono, J. N. Leaw, A. Reaves, M. Singh, X. Wang, S. Adam, and X. Gu, Solid State Communications 282, 38 (2018)

    2018

  11. [18]

    B. Lian, Z. Wang, and B. A. Bernevig, Phys. Rev. Lett. 122, 257002 (2019)

    2019

  12. [19]

    Y. W. Choi and H. J. Choi, Phys. Rev. B 98, 241412 (2018)

    2018

  13. [20]

    F. Wu, A. H. MacDonald, and I. Martin, Phys. Rev. Lett. 121, 257001 (2018)

    2018

  14. [21]

    Yudhistira, N

    I. Yudhistira, N. Chakraborty, G. Sharma, D. Y. Ho, E. Laksono, O. P. Sushkov, G. Vignale, and S. Adam, Phys. Rev. B 99, 140302(R) (2019)

    2019

  15. [22]

    F. Wu, E. Hwang, and S. Das Sarma, arXiv preprint arXiv:1811.04920 (2018)

    2018 arXiv

  16. [23]

    Grigorenko, M

    A. Grigorenko, M. Polini, and K. Novoselov, Nature Pho- tonics 6, 749 (2012)

    2012

  17. [24]

    Grabowski and L

    M. Grabowski and L. Sham, Phys. Rev. B 29, 6132 (1984)

    1984

  18. [25]

    Lewandowski and L

    C. Lewandowski and L. Levitov, arXiv preprint arXiv:1905.13088 (2019)

    2019 arXiv

  19. [26]

    Wunsch, T

    B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New Journal of Physics 8, 318 (2006)

    2006

  20. [27]

    E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 (2007)

    2007

  21. [28]

    Canright and G

    G. Canright and G. Vignale, Phys. Rev. B 39, 2740 (1989)

    1989

  22. [29]

    Eliashberg, Sov

    G. Eliashberg, Sov. Phys. JETP 11, 696 (1960)

    1960

Pith tools

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