Pith. sign in

REVIEW 4 major objections 5 minor 1 cited by

Using sign-problem-free quantum Monte Carlo, this paper establishes the neutrality phase diagram of twisted bilayer graphene: a continuous Dirac-semimetal to KIVC transition on approaching the magic angle, a strain-driven transition into an

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · deepseek-v4-flash

2026-08-01 09:49 UTC pith:LUMQLGIQ

load-bearing objection A solid finite-T QMC study of TBG at neutrality with a real tension between the T=0 phase diagram and the authors' own Fermi-surface instability caveat. the 4 major comments →

arxiv 2607.20619 v1 pith:LUMQLGIQ submitted 2026-07-22 cond-mat.str-el cond-mat.mes-hall

Monte Carlo Studies of Twisted Bilayer Graphene: Strain and Thermal Fluctuations

Johannes S. Hofmann , Patrick Ledwith , Ashvin Vishwanath , Jong Yeon Lee , Erez Berg This is my paper
classification cond-mat.str-el cond-mat.mes-hall PACS 71.27.+a71.30.+h
keywords twisted bilayer graphenequantum Monte Carlocharge neutralityKIVC orderentropy plateauheterostrainanisotropic semimetalMott semimetal
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Simulating the two narrow moiré bands of twisted bilayer graphene at charge neutrality with sign-problem-free determinant quantum Monte Carlo, the paper establishes a T=0 phase diagram in twist angle and uniaxial heterostrain: a Dirac semimetal gives way continuously to a gapped Kramers inter-valley coherent (KIVC) state as the twist angle approaches the magic angle, and under strain that KIVC state yields continuously to an anisotropic semimetal with gapless excitations at the moiré Brillouin zone center. The paper further shows that in the KIVC regime the entropy per moiré cell rises sharply with temperature and plateaus between about 15 and 40 K near kB ln(8 choose 4), the value expected for localized electrons with independent spin, valley, and orbital degrees of freedom — a Mott-like flavor-incoherent regime that survives even though the topological bands admit no symmetric localized Wannier basis. It also tracks the spectral function continuously from coherent K-point Dirac quasiparticles at large angle to Γ-centered gapless quasiparticles near the magic angle, and discusses how experiments can distinguish the anisotropic semimetal from this Mott semimetal. A sympathetic reader would care because these are numerically exact results in a strongly correlated regime where no expansion parameter exists, and because the predicted entropy plateau and spectral signatures are directly testable in twisting-microscope and thermodynamic measurements.

Core claim

The paper's central claim is that at charge neutrality the T→0 ground state of the projected TBG model is a gapped KIVC state in an intermediate twist-angle window, bounded on the large-angle side by a Dirac semimetal and on the small-angle side — only when uniaxial heterostrain is present — by an anisotropic semimetal whose gapless excitations sit near Γ_M rather than ±K_M; both transitions are continuous and are located by finite-size scaling of the KIVC correlation length, at θ≈1.23° and θ≈1.14° for the strain studied. At finite temperature, the paper claims the KIVC regime exhibits a broad entropy plateau, S ≈ kB ln(8 choose 4) per moiré cell for 15 K ≲ T ≲ 40 K, reflecting a Mott-like s

What carries the argument

The load-bearing object is the particle-hole-symmetric continuum model of the two flat moiré bands with projected dual-gated Coulomb interactions, formulated in momentum space to circumvent the Wannier obstruction; dropping the O(θ) rotation of the sublattice Pauli matrices makes the model exactly particle-hole symmetric and hence free of the fermion sign problem at charge neutrality, even under heterostrain. The numerical engine is determinant quantum Monte Carlo with a newly introduced hybrid Monte Carlo update that allows system sizes up to L=12; the KIVC order parameter, a momentum-space fermion bilinear with form factors, is used with correlation-length finite-size scaling to locate pha

Load-bearing premise

The whole construction rests on dropping the O(θ)≈0.02 rad rotation of the sublattice Pauli matrices to secure particle-hole symmetry and a sign-free simulation; if those terms, or the neglected remote bands and chosen screening parameters, appreciably shift the KIVC gap or the strain-induced Dirac-point motion, the predicted phase boundaries and 15–40 K entropy plateau could move or disappear.

What would settle it

Measure the specific heat (entropy) per moiré cell of a clean, unstrained near-magic-angle device at charge neutrality from 2 to 60 K: the claim requires a sharp entropy rise near 10 K and a plateau near kB ln(8 choose 4) up to ≈40 K; absence of the plateau would falsify the Mott-regime claim.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • The exact ground state at neutrality is KIVC between a Dirac semimetal and (under strain) an anisotropic semimetal; strain is therefore a control knob that can destroy the insulator without doping.
  • The entropy plateau at kB ln(8 choose 4) between 15 and 40 K is a distinctive thermodynamic fingerprint: measurements of the entropy per moiré cell can detect the Mott-like flavor-incoherent regime even where KIVC order is suppressed.
  • The spectral function's continuous evolution means that a gapped spectrum at K_M and a Γ_M-centered gap closing are both expected as θ approaches the magic angle; quantum twisting microscope data at low temperature can locate the phase boundaries.
  • Cooling the Mott semimetal should open a gap, whereas the anisotropic semimetal stays gapless to much lower T; a Zeeman field distinguishes them by spin-splitting vs. opening a gap ∝ U B/T.
  • The phase boundaries computed with and without the Hartree-Fock subtraction scheme coincide when strain is expressed as the energy splitting E_str, indicating the diagram is robust to the renormalization convention used for the single-particle dispersion.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The paper's finite-T results suggest a practical diagnostic: a specific-heat or entropy measurement in the 15–40 K window on a nominally unstrained device could detect local strain variations, since strained (anisotropic semimetal) regions have much lower entropy at the same temperature — effectively making the entropy plateau a strain microscope.
  • The reported reentrant enhancement of KIVC correlations near θ≈1.06°, where the strain-meandering Dirac cones meet at Γ_M, hints that tuning strain alone could drive an insulator-to-semimetal-to-insulator sequence at fixed twist angle; a dedicated finite-size scaling scan in that corner of the phase diagram would settle whether that is a genuine phase or a finite-size effect.
  • The coexistence of a gapped single-particle spectrum with a high-entropy local-moment plateau suggests that the entropy is carried by flavor (spin/valley/orbital) degrees of freedom that are invisible to single-particle probes; measurements of the spin susceptibility or magnetic entropy would isolate this contribution.
  • If the small-s² picture is right, the same concentrated-Berry-curvature mechanism that protects the Wannier obstruction also produces the nonlocal trion-like carriers responsible for the Γ_M spectral weight; this implies the 'missing' electron spectral weight in tunneling experiments could be recovered in three-particle (trion) correlation functions.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The manuscript reports sign-problem-free determinant quantum Monte Carlo simulations of a projected, particle-hole-symmetric continuum model of twisted bilayer graphene at charge neutrality, including uniaxial heterostrain and dual-gated Coulomb interactions. The central claims are: (i) at T=0 and zero strain there is a continuous Dirac-semimetal to Kramers inter-valley coherent (KIVC) transition as the twist angle decreases; (ii) under finite strain the KIVC phase gives way at smaller angles to an anisotropic semimetal with gapless excitations near the moiré Brillouin zone center; (iii) in the KIVC regime the entropy rises with temperature and plateaus near k_B ln C(8,4) between about 15 K and 40 K, indicating a flavor-incoherent Mott-like regime; and (iv) the single-particle spectral function evolves continuously with twist angle and at intermediate temperatures resembles either an anisotropic semimetal or a Mott semimetal. The paper includes a new hybrid Monte Carlo variant, details of the KIVC correlation-length analysis, entropy calculation, and an analytic small-s Appendix G that yields a Mott semimetal self-energy.

Significance. If the central results hold, this is a valuable numerically exact benchmark for TBG at charge neutrality: it treats interactions in a projected model without a Wannier basis, extends to L=12 via a new HMC approach, and compares renormalized and unrenormalized band-structure schemes. The entropy plateau compared with the parameter-free combinatorial value ln C(8,4) is a clean and falsifiable prediction, and the spectral-function comparison to an analytic self-energy is a useful cross-check. However, the paper's own Discussion contains a statement that directly undermines the T=0 phase diagram: the anisotropic semimetal is said to be expected to be unstable at T=0 toward excitonic order, which contradicts the abstract and Fig. 1(a). This, together with missing error bars in the central entropy calculation and unresolved finite-size issues in the reentrant KIVC region, means the paper requires substantive revision before its main claims are established.

major comments (4)
  1. [Discussion (final paragraph); Fig. 1(a); End Matter A.1] The Discussion says the anisotropic semimetal 'exhibits (small) electron and hole Fermi surfaces; it is thus expected to be unstable at T=0 [103] towards an excitonic order, most likely some form of IVC.' This is in direct tension with the abstract and Fig. 1(a), where the anisotropic semimetal is presented as a T=0 ground state and the KIVC-to-anisotropic-semimetal boundary near θ≈1.14° is called a continuous transition. If Ref. [103] applies to this model, the T=0 ground state in that region should be ordered, not semimetallic, and the transition would be between two ordered states. The phase identification in Fig. 4(b) rests on the absence of KIVC order at T=1.8 K for L≤12, not on a measurement of the competing excitonic order parameter. This must be resolved: either the phase diagram should be explicitly reframed as a finite-T (1.8 K) phase diagram, or the excitonic instability shoul
  2. [End Matter A.1, Appendix F (Fig. 8), Appendix D] The reentrant KIVC correlations at θ=1.06° are acknowledged in footnote [106] to require finer θ and larger L to establish, and Appendix F admits that 'we cannot reach system sizes large enough to faithfully resolve those Fermi surfaces' of the anisotropic semimetal. Since the KIVC-to-anisotropic-semimetal boundary is extracted from finite-size crossings at L≤12, the possibility of a reentrant or competing ordered phase in the small-θ strained region is not excluded. The phase boundary drawn in Fig. 1(a)/Fig. 4(a) is therefore less secure than the central claim requires. The authors should either provide a quantitative finite-size analysis of the reentrant region or explicitly state the resulting uncertainty in the phase diagram.
  3. [Appendix F, Eq. (F3); Fig. 2(b-d); Fig. 8(b-d)] The entropy curves are obtained from Eq. (F3), which integrates the internal energy ⟨H⟩_β over inverse temperature starting from the infinite-temperature limit. This is numerically delicate because the integrand is a difference of large energies at low T, yet no error bars, jackknife estimates, or sensitivity checks (e.g., dependence on β range, Trotter step, or HMC trajectory length) are shown in any entropy panel. Since the ln C(8,4) plateau is one of the three headline results, the existence, width, and value of the plateau need to be supported by error bars or an equivalent uncertainty analysis.
  4. [Appendix G, Eqs. (G11)-(G16); Fig. 3] The 'Mott semimetal' identification at T=30 K rests on comparing DQMC spectra with the analytic self-energy Eq. (G16), derived from the small-s wavefunction ansatz Eq. (G2) with s≈0.25 and the classical Hamiltonian Eq. (G11). At s≈0.25 the expansion parameter s²≈0.06 is not extremely small, so the agreement in Fig. 3 is a consistency check rather than a rigorous confirmation. Furthermore, the statement that the Γ-point gap 'closely tracks' the single-particle dispersion may be partly built into the analytic Green's function (G17), whose quasiparticle branch approaches ω=E_BM(1-|λ_k|²) as k→0. The authors should provide a quantitative measure of agreement between the DQMC spectrum and Eq. (G17), or soften the claim that the QMC spectra identify a Mott semimetal.
minor comments (5)
  1. [Abstract] Typographical errors: 'The later can be a anisotropic' should be 'The latter can be an anisotropic'; 'distiguish' should be 'distinguish'; 'heterstrain' appears elsewhere and should be 'heterostrain'.
  2. [Fig. 1 caption] The caption says 'green symbols indicate the parameters of panels (b-d)', but the green symbols are not clearly visible or labeled in the rendered figure. Please ensure the markers are identifiable.
  3. [Appendix D, Eq. (D2)] The moiré lattice constant a_M is used in Eq. (D2) but not explicitly defined before this equation. Please define it in the text or in the caption.
  4. [Fig. 2 caption] The caption refers to 'dashed horizontal lines' representing entropy from the flat-band-limit ground-state degeneracy, but does not explain which line corresponds to which system size or how the quoted U(4) degeneracy is computed. Please clarify.
  5. [End Matter A.2 / Appendix B] The distinction between 'renormalized' and 'unrenormalized' bands is central to Fig. 1(a) and Appendix A.2, but the main text could make clearer which scheme is used for the later quantitative claims (entropy, spectral functions). A one-sentence summary in the main text would help.

Circularity Check

0 steps flagged

No significant circularity: central phase diagram and entropy plateau are direct QMC measurements; analytic comparison curves are not fitted to QMC data.

full rationale

The paper's central claims are extracted directly from sign-problem-free DQMC: the phase diagram comes from finite-size scaling of the KIVC correlation length ξ_KIVC (Appendix D, Eq. D2), the entropy is obtained by integrating ⟨H⟩ over β (Eq. F3), and the plateau is compared to the parameter-free value k_B ln C(8,4), not to a fitted model. The atomic-limit Hubbard U=33/35 meV is explicitly fit to the entropy for a comparison curve only ('the black curve represents the high temperature behavior obtained from a single-site SU(8) Hubbard model with U=33meV and U=35meV'), and is not used to generate a prediction. The analytic Mott self-energy Eq. (G16) is adopted from self-cited Refs [71,86], but the small parameter s≈0.25 was fixed by matching the single-particle TBG Bloch wavefunctions, not by fitting QMC observables; the QMC spectral functions are then compared to these analytic curves, so the comparison is a genuine check rather than a reduction-by-construction. Self-citations to [69] for the sign-problem-free property are supported by the explicit particle-hole symmetry transformation stated in Appendix B, and the QMC algorithm is described self-containedly in Appendix C. The only notable tension is a correctness/consistency issue, not circularity: the Discussion states the anisotropic semimetal is 'expected to be unstable at T=0 [103] towards an excitonic order, most likely some form of IVC,' which appears to conflict with the T=0 phase diagram labeling that region as a stable semimetal. This is flagged as a scientific risk (the boundary is inferred from absence of KIVC order at T=1.8 K and limited L), but it is not an equation-level circular reduction. Model approximations such as dropping θ-order terms (Appendix B) and the HF subtraction scheme are stated assumptions, not circular steps.

Axiom & Free-Parameter Ledger

6 free parameters · 7 axioms · 1 invented entities

The central QMC results do not rely on fitted parameters for the phase boundaries. The listed free parameters are comparison/interpretation aids (Hubbard U, small-s ansatz) or model choices (screening, strain coupling, 0.05° shift). The main load-bearing assumptions are the projected flat-band Hilbert space, the particle-hole-symmetric single-particle Hamiltonian, and finite-size/low-temperature scaling.

free parameters (6)
  • Hubbard U (atomic-limit comparison) = 33 meV (θ=1.1°), 35 meV (θ=1.15°)
    Chosen so the single-site SU(8) Hubbard entropy matches the high-temperature QMC entropy (Fig. 2b-c); used only for comparison, not for the plateau value.
  • Small-s wavefunction parameter = s≈0.25
    Set in Ref [71] by matching the analytic ansatz to TBG Bloch wavefunctions to within ~5%; used here in Appendix G and Fig. 3 to plot the Mott-semimetal spectral function. Not fitted to QMC data.
  • Angular shift Δθ between renormalized and unrenormalized schemes = 0.05°
    Applied to the non-renormalized phase boundary in Fig. 1(a) to account for magic-angle renormalization; hand-chosen so the two schemes' boundaries line up.
  • Screening parameters (ε_r, d0) = ε_r=10, d0=20 nm
    Hand-chosen dual-gate screening model; sets interaction scale V0. Qualitative phase diagram likely robust, quantitative boundaries may shift.
  • Strain hopping parameter γ = γ=2
    Chosen from the graphene literature for hopping variation under bond-length change; enters the pseudo-gauge potential Eq. (2) and thus the strain phase boundary.
  • Poisson ratio ν = 0.16
    Material parameter for the uniaxial strain tensor in Eq. (1); standard value for graphene.
axioms (7)
  • domain assumption The Bistritzer-MacDonald continuum model projected to the two narrow moiré bands per spin/valley captures the physics of TBG near the magic angle.
    All QMC results are computed in this projected Hilbert space; remote bands and higher moiré bands are neglected. Invoked in 'Model and Method' and Appendix B.
  • domain assumption After dropping θ-order terms, the Hamiltonian has exact particle-hole symmetry C (Eq. B3), which guarantees absence of the sign problem at charge neutrality.
    The θ-dependent rotation of Pauli matrices is dropped, 'amounting to dropping terms of order θ≈1.1°≈0.02 rad' (Appendix B). If these terms matter, the sign-problem-free property and the resulting phase boundaries are approximate.
  • domain assumption T=1.8 K is sufficiently low, and finite-size scaling of ξ_KIVC/L with L up to 12 correctly identifies the T=0 continuous transitions.
    Phase boundaries are located by correlation-length crossings (Appendix D, Fig. 4b). This assumes the quantum-critical scaling regime is reached at T=1.8 K and L≤12.
  • standard math The entropy can be obtained by integrating the temperature derivative of ⟨H⟩ from the infinite-temperature limit S(β0)=8k_B L^2 ln 2 (Appendix F).
    Uses the thermodynamic identity S=-∂Ω/∂T and integration of specific-heat-like data; the authors note ln Z cannot be measured directly and the derivative is numerically challenging.
  • domain assumption Maximum-entropy analytic continuation of imaginary-time QMC data provides reliable real-frequency spectral functions at the displayed energy scales.
    Spectral functions in Figs. 1, 3, 4, and 7 are extracted via stochastic maximum entropy (Beach method, Ref [78]); no validation against exact spectra is shown.
  • ad hoc to paper The small-s wavefunction ansatz (Eq. G2, s≈0.25) and the resulting classical Hamiltonian (Eq. G11) describe the flavor-incoherent state and yield the Mott semimetal self-energy (Eq. G16).
    Borrowed from Refs [71,86]; s is fixed by matching TBG wavefunctions to ~5%, not by QMC data. It underpins the interpretation of the Γ-centered gapless branch and the experimental discrimination proposals.
  • domain assumption The Hartree-Fock subtraction scheme (Eq. B5) using the decoupled-layer HF potential avoids double counting without biasing the correlated phase diagram.
    The renormalized band structure changes E_str(ε) significantly; the phase diagram is presented for this scheme, with the unrenormalized scheme shifted by 0.05°.
invented entities (1)
  • Dirac trion excitation no independent evidence
    purpose: Explains the zero-quasiparticle-residue spectral branch near ΓM in the Mott semimetal, giving a three-particle operator F_k that couples to the electron (Eq. G18).
    Introduced in cited Ref [86], not this paper; no three-particle spectral function is computed here, and the paper states resolving it is future work.

pith-pipeline@v1.3.0-alltime-deepseek · 26119 in / 20799 out tokens · 162615 ms · 2026-08-01T09:49:06.222311+00:00 · methodology

0 comments
read the original abstract

We study the phase diagram of twisted bilayer graphene at charge neutrality as a function of twist angle $\theta$, uniaxial heterostrain $\varepsilon$, and temperature $T$ using sign-problem-free quantum Monte Carlo simulations. At $T=0$ and zero strain, we find a continuous transition from a Dirac semimetal to a gapped Kramers inter-valley coherent (KIVC) phase as $\theta$ decreases toward the magic angle. With finite strain, the KIVC phase undergoes a further continuous transition at smaller $\theta$ into an anisotropic semimetal with gapless excitations near the center of the moir\'{e} Brillouin zone. In the KIVC regime, the entropy rises sharply with temperature and plateaus at $15\,\text{K} \lesssim T \lesssim 40\,\text{K}$ near the value expected from a Mott-like regime of localized electrons with nearly uncorrelated spin, valley, and orbital degrees of freedom, despite the topological obstruction preventing a localized tight-binding description of the active bands. The spectral function evolves continuously with $\theta$: at low $T$, a gap opens at the $K$ points and the minimal gap shifts to $\Gamma$ as $\theta$ decreases; at intermediate $T$, the spectral function smoothly interpolates between a Dirac semimetal spectrum with coherent $K$-point quasiparticles and a spectrum with gapless $\Gamma$-centered quasiparticles near the magic angle. The later can be a anisotropic or a Mott semimetal and we discuss how to distiguish them in experiment.

Figures

Figures reproduced from arXiv: 2607.20619 by Ashvin Vishwanath, Erez Berg, Johannes S. Hofmann, Jong Yeon Lee, Patrick Ledwith.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (b-c) show the entropy per moir´e unit cell ver￾sus temperature for twist angles θ = 1.1 ◦ , 1.15◦ , and 1.3 ◦ . The first two angles are within the KIVC regime, whereas the third angle is in the Dirac semimetal phase. For θ = 1.1 ◦ , 1.15◦ and below TKIVC ≈ 10K the entropy is small, and close to the value expected from the ground state degeneracy due to the approximate U(4) symme￾try for the corresponding… view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Meandering pair ( [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Additional data for the KIVC correlation length [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Correlation functions at [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. (a) Finite temperature phase diagram at non-zero strain, [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

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

  1. Trion Excitations in Twisted Bilayer Graphene: A Quantum Monte Carlo Study

    cond-mat.str-el 2026-07 conditional novelty 7.0

    Finite-temperature QMC of twisted bilayer graphene finds gapless 'Dirac trion' three-particle excitations in the normal state, exactly orthogonal to electrons at the Γ point.

Reference graph

Works this paper leans on

122 extracted references · 1 canonical work pages · cited by 1 Pith paper

  1. [1]

    G. Li, A. Luican, J. Lopes dos Santos, A. Castro Neto, A. Reina, J. Kong, and E. Andrei, Observation of Van Hove singularities in twisted graphene layers, Nature physics6, 109 (2010)

    2010

  2. [2]

    Bistritzer and A

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

    2011

  3. [3]

    J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, Continuum model of the twisted graphene bilayer, Phys. Rev. B86, 155449 (2012)

    2012

  4. [4]

    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, Nature556, 80 (2018)

    2018

  5. [5]

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

    2018

  6. [6]

    Balents, C

    L. Balents, C. R. Dean, D. K. Efetov, and A. F. Young, Superconductivity and strong correlations in moir´ e flat bands, Nature Physics16, 725 (2020)

    2020

  7. [7]

    H. Yang, L. Liu, H. Yang, Y. Zhang, X. Wu, Y. Huang, H.-J. Gao, and Y. Wang, Advance in two-dimensional twisted moir´ e materials: Fabrication, properties, and applications, Nano Research16, 2579 (2023)

    2023

  8. [8]

    B. A. Bernevig and D. K. Efetov, Twisted bilayer graphene’s gallery of phases, Physics Today77, 38 (2024)

    2024

  9. [9]

    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, Nature574, 653 (2019)

    2019

  10. [10]

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

    2019

  11. [11]

    Y. Choi, J. Kemmer, Y. Peng, A. Thomson, H. Arora, R. Polski, Y. Zhang, H. Ren, J. Alicea, G. Refael, F. von Oppen, K. Watanabe, T. Taniguchi, and S. Nadj-Perge, Electronic correlations in twisted bilayer graphene near the magic angle, Nature Physics 10.1038/s41567-019- 0606-5 (2019)

    doi:10.1038/s41567-019- 2019

  12. [12]

    Jiang, X

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

    2019

  13. [13]

    Serlin, C

    M. Serlin, C. Tschirhart, H. Polshyn, Y. Zhang, J. Zhu, K. Watanabe, T. Taniguchi, L. Balents, and A. Young, Intrinsic quantized anomalous hall effect in a moir´ e het- erostructure, Science367, 900 (2020)

    2020

  14. [14]

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

    2019

  15. [15]

    E. Y. Andrei and A. H. MacDonald, Graphene bilayers with a twist, Nature materials19, 1265 (2020)

    2020

  16. [16]

    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

  17. [17]

    D. Wong, K. P. Nuckolls, M. Oh, B. Lian, Y. Xie, S. Jeon, K. Watanabe, T. Taniguchi, B. A. Bernevig, and A. Yazdani, Cascade of electronic transitions in magic-angle twisted bilayer graphene, Nature582, 198 (2020)

    2020

  18. [18]

    Y. Xie, A. T. Pierce, J. M. Park, D. E. Parker, E. Kha- laf, P. Ledwith, Y. Cao, S. H. Lee, S. Chen, P. R. For- rester,et al., Fractional chern insulators in magic-angle twisted bilayer graphene, Nature600, 439 (2021)

    2021

  19. [19]

    J. Yu, B. A. Foutty, Z. Han, M. E. Barber, Y. Schattner, K. Watanabe, T. Taniguchi, P. Phillips, Z.-X. Shen, S. A. Kivelson,et al., Correlated hofstadter spectrum and flavour phase diagram in magic-angle twisted bi- layer graphene, Nature Physics18, 825 (2022)

    2022

  20. [20]

    Grover, M

    S. Grover, M. Bocarsly, A. Uri, P. Stepanov, G. Di Bat- tista, I. Roy, J. Xiao, A. Y. Meltzer, Y. Myasoedov, K. Pareek,et al., Chern mosaic and berry-curvature magnetism in magic-angle graphene, Nature physics18, 885 (2022)

    2022

  21. [21]

    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 wavefunctions in magic-angle graphene, Na- ture620, 525 (2023)

    2023

  22. [22]

    Guinea and N

    F. Guinea and N. R. Walet, Electrostatic effects, band distortions, and superconductivity in twisted graphene bilayers, Proceedings of the National Academy of Sci- ences115, 13174 (2018)

    2018

  23. [23]

    H. C. Po, L. Zou, A. Vishwanath, and T. Senthil, Ori- gin of Mott Insulating Behavior and Superconductivity in Twisted Bilayer Graphene, Phys. Rev. X8, 031089 (2018)

    2018

  24. [24]

    Z. Song, Z. Wang, W. Shi, G. Li, C. Fang, and B. A. Bernevig, All magic angles in twisted bilayer graphene are topological, Phys. Rev. Lett.123, 036401 (2019)

    2019

  25. [26]

    Koshino, N

    M. Koshino, N. F. Q. Yuan, T. Koretsune, M. Ochi, K. Kuroki, and L. Fu, Maximally localized wannier or- bitals and the extended hubbard model for twisted bi- layer graphene, Phys. Rev. X8, 031087 (2018)

    2018

  26. [27]

    N. F. Q. Yuan and L. Fu, Model for the metal-insulator transition in graphene superlattices and beyond, Phys. Rev. B98, 045103 (2018)

    2018

  27. [28]

    J. F. Dodaro, S. A. Kivelson, Y. Schattner, X. Q. Sun, and C. Wang, Phases of a phenomenological model of twisted bilayer graphene, Phys. Rev. B98, 075154 (2018)

    2018

  28. [29]

    Padhi, C

    B. Padhi, C. Setty, and P. W. Phillips, Doped twisted bi- layer graphene near magic angles: Proximity to wigner crystallization, not mott insulation, Nano letters18, 6175 (2018)

    2018

  29. [30]

    Tarnopolsky, A

    G. Tarnopolsky, A. J. Kruchkov, and A. Vishwanath, Origin of magic angles in twisted bilayer graphene, Phys. Rev. Lett.122, 106405 (2019)

    2019

  30. [31]

    S. Carr, S. Fang, Z. Zhu, and E. Kaxiras, Exact contin- uum model for low-energy electronic states of twisted bilayer graphene, Phys. Rev. Res.1, 013001 (2019)

    2019

  31. [32]

    S. Carr, S. Fang, H. C. Po, A. Vishwanath, and E. Kaxi- ras, Derivation of wannier orbitals and minimal-basis tight-binding hamiltonians for twisted bilayer graphene: First-principles approach, Phys. Rev. Res.1, 033072 (2019)

    2019

  32. [33]

    H. C. Po, L. Zou, T. Senthil, and A. Vishwanath, Faith- ful tight-binding models and fragile topology of magic- angle bilayer graphene, Phys. Rev. B99, 195455 (2019)

    2019

  33. [34]

    Rademaker, D

    L. Rademaker, D. A. Abanin, and P. Mellado, Charge smoothening and band flattening due to hartree correc- tions in twisted bilayer graphene, Phys. Rev. B100, 205114 (2019)

    2019

  34. [35]

    T. Cea, N. R. Walet, and F. Guinea, Twists and the electronic structure of graphitic materials,Nano Letters, Nano Letters19, 8683 (2019)

    2019

  35. [36]

    K. Seo, V. N. Kotov, and B. Uchoa, Ferromagnetic mott state in twisted graphene bilayers at the magic angle, Phys. Rev. Lett.122, 246402 (2019)

    2019

  36. [37]

    A. O. Sboychakov, A. V. Rozhkov, A. L. Rakhmanov, and F. Nori, Many-body effects in twisted bilayer graphene at low twist angles, Phys. Rev. B100, 045111 (2019)

    2019

  37. [38]

    Thomson, S

    A. Thomson, S. Chatterjee, S. Sachdev, and M. S. Scheurer, Triangular antiferromagnetism on the honey- comb lattice of twisted bilayer graphene, Phys. Rev. B 98, 075109 (2018)

    2018

  38. [39]

    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 dissipation, Phys. Rev. Lett.124, 076801 (2020)

    2020

  39. [40]

    Xie and A

    M. Xie and A. H. MacDonald, Nature of the correlated insulator states in twisted bilayer graphene, Phys. Rev. Lett.124, 097601 (2020)

    2020

  40. [41]

    Wu and S

    F. Wu and S. Das Sarma, Collective excitations of quantum anomalous hall ferromagnets in twisted bilayer graphene, Phys. Rev. Lett.124, 046403 (2020)

    2020

  41. [42]

    B. A. Bernevig, Z.-D. Song, N. Regnault, and B. Lian, 7 Twisted bilayer graphene. i. matrix elements, approx- imations, perturbation theory, and ak·ptwo-band model, Phys. Rev. B103, 205411 (2021)

    2021

  42. [43]

    P. J. Ledwith, E. Khalaf, and A. Vishwanath, Strong coupling theory of magic-angle graphene: A pedagogical introduction, Annals of Physics435, 168646 (2021)

    2021

  43. [44]

    S. Liu, E. Khalaf, J. Y. Lee, and A. Vishwanath, Ne- matic topological semimetal and insulator in magic- angle bilayer graphene at charge neutrality, Phys. Rev. Res.3, 013033 (2021)

    2021

  44. [45]

    D. E. Parker, T. Soejima, J. Hauschild, M. P. Zaletel, and N. Bultinck, Strain-induced quantum phase tran- sitions in magic-angle graphene, Phys. Rev. Lett.127, 027601 (2021)

    2021

  45. [46]

    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 bi- layer graphene, Phys. Rev. X11, 041063 (2021)

    2021

  46. [47]

    Shavit, E

    G. Shavit, E. Berg, A. Stern, and Y. Oreg, Theory of correlated insulators and superconductivity in twisted bilayer graphene, Phys. Rev. Lett.127, 247703 (2021)

    2021

  47. [48]

    Liu and X

    J. Liu and X. Dai, Theories for the correlated insulating states and quantum anomalous hall effect phenomena in twisted bilayer graphene, Phys. Rev. B103, 035427 (2021)

    2021

  48. [49]

    Wagner, Y

    G. Wagner, Y. H. Kwan, N. Bultinck, S. H. Simon, and S. A. Parameswaran, Global phase diagram of the nor- mal state of twisted bilayer graphene, Phys. Rev. Lett. 128, 156401 (2022)

    2022

  49. [50]

    Vafek and J

    O. Vafek and J. Kang, Continuum effective hamiltonian for graphene bilayers for an arbitrary smooth lattice de- formation from microscopic theories, Phys. Rev. B107, 075123 (2023)

    2023

  50. [51]

    C. N. Breiø and B. M. Andersen, Chern insulator phases and spontaneous spin and valley order in a moir´ e lattice model for magic-angle twisted bilayer graphene, Phys. Rev. B107, 165114 (2023)

    2023

  51. [52]

    Sheffer, R

    Y. Sheffer, R. Queiroz, and A. Stern, Symmetries as the guiding principle for flattening bands of dirac fermions, Phys. Rev. X13, 021012 (2023)

    2023

  52. [53]

    J. Zhu, I. Torre, M. Polini, and A. H. MacDonald, Weak- coupling theory of magic-angle twisted bilayer graphene, Phys. Rev. B110, L121117 (2024)

    2024

  53. [54]

    K. D. Stubbs, M. Ragone, A. H. MacDonald, and L. Lin, The many-body ground state manifold of flat band in- teracting hamiltonian for magic angle twisted bilayer graphene (2025), arXiv:2503.20060 [math-ph]

    Pith/arXiv arXiv 2025

  54. [55]

    Haule, E

    M. Haule, E. Y. Andrei, and K. Haule, The Mott- semiconducting state in the magic angle bilayer graphene (2019), arXiv:1901.09852 [cond-mat.str-el]

    Pith/arXiv arXiv 2019

  55. [56]

    Kang and O

    J. Kang and O. Vafek, Non-abelian dirac node braiding and near-degeneracy of correlated phases at odd integer filling in magic-angle twisted bilayer graphene, Phys. Rev. B102, 035161 (2020)

    2020

  56. [57]

    Potasz, M

    P. Potasz, M. Xie, and A. H. MacDonald, Exact di- agonalization for magic-angle twisted bilayer graphene, Phys. Rev. Lett.127, 147203 (2021)

    2021

  57. [58]

    Da Liao, J

    Y. Da Liao, J. Kang, C. N. Breiø, X. Y. Xu, H.-Q. Wu, B. M. Andersen, R. M. Fernandes, and Z. Y. Meng, Correlation-induced insulating topological phases at charge neutrality in twisted bilayer graphene, Phys. Rev. X11, 011014 (2021)

    2021

  58. [59]

    F. Xie, A. Cowsik, Z.-D. Song, B. Lian, B. A. Bernevig, and N. Regnault, Twisted bilayer graphene. vi. an exact diagonalization study at nonzero integer filling, Phys. Rev. B103, 205416 (2021)

    2021

  59. [60]

    Pahlevanzadeh, P

    B. Pahlevanzadeh, P. Sahebsara, and D. S´ en´ echal, Charge order and antiferromagnetism in twisted bilayer graphene from the variational cluster approximation, SciPost Phys.13, 040 (2022)

    2022

  60. [61]

    F. Xie, J. Kang, B. A. Bernevig, O. Vafek, and N. Reg- nault, Phase diagram of twisted bilayer graphene at fill- ing factorν=±3, Phys. Rev. B107, 075156 (2023)

    2023

  61. [62]

    T. Wang, D. E. Parker, T. Soejima, J. Hauschild, S. Anand, N. Bultinck, and M. P. Zaletel, Ground-state order in magic-angle graphene at fillingν=−3: A full- scale density matrix renormalization group study, Phys. Rev. B108, 235128 (2023)

    2023

  62. [63]

    Zhang, G

    X. Zhang, G. Pan, B.-B. Chen, H. Li, K. Sun, and Z. Y. Meng, Polynomial sign problem and topological mott insulator in twisted bilayer graphene, Phys. Rev. B107, L241105 (2023)

    2023

  63. [64]

    F. M. Faulstich, K. D. Stubbs, Q. Zhu, T. Soejima, R. Dilip, H. Zhai, R. Kim, M. P. Zaletel, G. K.-L. Chan, and L. Lin, Interacting models for twisted bilayer graphene: A quantum chemistry approach, Phys. Rev. B107, 235123 (2023)

    2023

  64. [65]

    Datta, M

    A. Datta, M. J. Calderon, A. Camjayi, and E. Bascones, Heavy quasiparticles and cascades without symmetry breaking in twisted bilayer graphene, Nature Commu- nications14, 5036 (2023)

    2023

  65. [66]

    G. Rai, L. Crippa, D. C˘ alug˘ aru, H. Hu, F. Paoletti, L. de’ Medici, A. Georges, B. A. Bernevig, R. Valent ´ ı, G. Sangiovanni, and T. Wehling, Dynamical correla- tions and order in magic-angle twisted bilayer graphene, Phys. Rev. X14, 031045 (2024)

    2024

  66. [67]

    M. J. Calderon, A. Camjayi, A. Datta, and E. Bas- cones, Cascades in transport and optical conductivity of twisted bilayer graphene (2025), arXiv:2412.20855 [cond-mat.str-el]

    Pith/arXiv arXiv 2025

  67. [68]

    Xiao and S

    Z.-Y. Xiao and S. Zhang, Correlation effects in magic- angle twisted bilayer graphene: An auxiliary-field quan- tum monte carlo study, Phys. Rev. Res.7, 013103 (2025)

    2025

  68. [69]

    J. S. Hofmann, E. Khalaf, A. Vishwanath, E. Berg, and J. Y. Lee, Fermionic Monte Carlo Study of a Realistic Model of Twisted Bilayer Graphene, Phys. Rev. X12, 011061 (2022)

    2022

  69. [70]

    Zhang, G

    X. Zhang, G. Pan, Y. Zhang, J. Kang, and Z. Y. Meng, Momentum space quantum Monte Carlo on twisted bilayer graphene, Chinese Physics Letters38, 077305 (2021)

    2021

  70. [71]

    P. J. Ledwith, J. Dong, A. Vishwanath, and E. Khalaf, Nonlocal Moments and Mott Semimetal in the Chern Bands of Twisted Bilayer Graphene, Phys. Rev. X15, 021087 (2025)

    2025

  71. [72]

    Hu, Z.-D

    H. Hu, Z.-D. Song, and B. A. Bernevig, Projected and Solvable: Topological Heavy Fermion Model of Twisted Bilayer Graphene (2025), arXiv:2502.14039 [cond-mat.str-el]

    arXiv 2025

  72. [73]

    Such particle-hole symmetry breaking terms contain those originating from the rotation of the Dirac Pauli matrices, proportional toθ

  73. [74]

    Z. Wang, M. P. Zaletel, R. S. K. Mong, and F. F. As- saad, Phases of the (2+1) Dimensional SO(5) Nonlinear Sigma Model with Topological Term, Phys. Rev. Lett. 126, 045701 (2021)

    2021

  74. [75]

    Balents, General continuum model for twisted bilayer 8 graphene and arbitrary smooth deformations, SciPost Physics7, 048 (2019)

    L. Balents, General continuum model for twisted bilayer 8 graphene and arbitrary smooth deformations, SciPost Physics7, 048 (2019)

    2019

  75. [76]

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

    2019

  76. [77]

    Antebi, A

    O. Antebi, A. Stern, and E. Berg, In-plane orbital mag- netization as a probe for symmetry breaking in strained twisted bilayer graphene, Phys. Rev. B105, 104423 (2022)

    2022

  77. [78]

    K. S. D. Beach, Identifying the maximum entropy method as a special limit of stochastic analytic con- tinuation, arXiv e-prints , cond-mat/0403055 (2004), arXiv:cond-mat/0403055 [cond-mat.str-el]

    Pith/arXiv arXiv 2004

  78. [79]

    Bultinck, E

    N. Bultinck, E. Khalaf, S. Liu, S. Chatterjee, A. Vish- wanath, and M. P. Zaletel, Ground state and hidden symmetry of magic-angle graphene at even integer fill- ing, Phys. Rev. X10, 031034 (2020)

    2020

  79. [80]

    In realistic samples, weak inter-valley Hund’s coupling breaks this symmetry and lifts the degeneracy, se- lecting either a singlet KIVC—which then exhibits a finite-temperature Berezinskii–Kosterlitz–Thouless transition—or a triplet KIVC, which remains disordered for allT >0

  80. [81]

    M. J. Calder´ on and E. Bascones, Interactions in the 8- orbital model for twisted bilayer graphene, Phys. Rev. B102, 155149 (2020)

    2020

Showing first 80 references.