pith. sign in

arxiv: 2605.20677 · v2 · pith:D54Z6FREnew · submitted 2026-05-20 · ❄️ cond-mat.str-el

Strain-Tuned Incommensurate Kekul\'e Spiral Order in Twisted Bilayer Graphene: a Quantum Many-Body Study

Cheng Huang , Yves H. Kwan , Maksim Ulybyshev , Fakher F. Assaad , Laura Classen , Zi Yang Meng This is my paper

Pith reviewed 2026-06-30 17:38 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords twisted bilayer grapheneKramers intervalley coherentincommensurate Kekulé spiralquantum Monte Carlostrain tuningflat bandmany-body statesphase transition
0
0 comments X

The pith

At filling factors ±2, strain drives a transition from the Kramers intervalley coherent state to the incommensurate Kekulé spiral state in twisted bilayer graphene.

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

The paper establishes that the interacting ground state at ν = ±2 undergoes a strain-tuned transition from the KIVC state to the IKS state. This is determined through a combination of quantum Monte Carlo simulations adjusted for the sign problem, exact diagonalization, and Hartree-Fock calculations in a projected flat-band model. A sympathetic reader would care because it provides a non-perturbative understanding of the phases away from charge neutrality, where insulating and superconducting behaviors are observed experimentally. The work demonstrates how strain can be used to control the many-body order in these systems.

Core claim

The study shows that increasing strain at ν = ±2 induces a transition from the Kramers intervalley coherent (KIVC) state to the incommensurate Kekulé spiral (IKS) state, as revealed by the combined QMC, ED, and HF protocol in the correlated flat-band setting.

What carries the argument

The strain dependence of the ground state in the projected correlated flat-band model at ν=±2, computed via adjusted continuous-field momentum-space quantum Monte Carlo together with ED and HF.

If this is right

  • The KIVC state is stable at low strain while IKS dominates at higher strain.
  • The phase boundary between KIVC and IKS can be located using the computational protocol.
  • This applies to understanding the rich phase diagram of twisted bilayer graphene at fillings away from neutrality.
  • Similar transitions may occur in other strongly-correlated flat-band systems.

Where Pith is reading between the lines

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

  • The transition might be observable in transport or spectroscopy experiments by varying strain.
  • Extending the method could allow study of superconductivity in these phases.
  • Strain tuning could be a general control knob for intervalley coherent orders in moiré systems.

Load-bearing premise

The approximate treatment of the sign problem in the quantum Monte Carlo method still produces reliable ground-state properties and phase boundaries at ν = ±2.

What would settle it

Experimental measurement of the strain value at which the dominant order switches from KIVC to IKS at filling factor 2 or -2 would confirm or refute the transition.

Figures

Figures reproduced from arXiv: 2605.20677 by Cheng Huang, Fakher F. Assaad, Laura Classen, Maksim Ulybyshev, Yves H. Kwan, Zi Yang Meng.

Figure 1
Figure 1. Figure 1: Strained Brillouin zones, shifted Dirac points, and occupa￾tion number comparison between ED and approximated QMC. (a) Schematic Brillouin zones (green and orange hexagons) of two graphene layers with interlayer twist 𝜃 and uniaxial heterostrain of ±𝜖s/2 for each layer. The corresponding strained moiré Brillouin zone (mBZ) in the 𝜂 = + valley can be represented by the blue hexagon or black rhombus, with 𝚪 … view at source ↗
Figure 2
Figure 2. Figure 2: Occupation number 𝑛 𝜂 (k) and corresponding nesting diagnostics from QMC and HF at 𝜈 = −2 with 𝑁k = 18 × 18 and strain strength 𝜖s = 0.6%. The left panels are the QMC data, where (a) and (c) show the occupation factor 𝑛 𝜂 (k). The magenta arrow in (c) points from the minimum in (a) to the maximum in (c). The magneta arrow in (e) [(g)] corresponds to the maximum [minimum] of 𝑂 ′ (q) [𝑂 ′′(q)]. The right pan… view at source ↗
Figure 4
Figure 4. Figure 4: Intervalley structure factor 𝑆IKS (Q) from HF and QMC at 𝜈 = −2 with 𝑁k = 18×18 and strain strength 𝜖s = 0.6%. (a) The HF solution has Bragg peaks corresponding to qIKS = (7/18, 7/18). The magenta arrow points to the peak and is the same as those in [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: Occupation number 𝑛 𝜂 k and corresponding nesting di￾agnostics from QMC and HF at 𝜈 = −2 with 𝑁k = 18 × 18 and strain strength 𝜖s = 0.3%. The left panels are the QMC data, where (a) and (c) show the occupation factor 𝑛 𝜂 (k). The magenta arrow in (c) points from the minimum in (a) to the maximum in (c). The magneta arrow in (e) [(g)] corresponds to the maximum [minimum] of 𝑂 ′ (q) [𝑂 ′′(q)]. The right pane… view at source ↗
Figure 5
Figure 5. Figure 5: Occupation number 𝑛 𝜂 (k) and corresponding nesting diagnostics from QMC and HF at 𝜈 = −2 with 𝑁k = 18 × 18 and strain strength 𝜖s = 0. The left panels are the QMC data, where (a) and (c) show the occupation factor 𝑛 𝜂 (k). Panels (e) and (g) correspond to 𝑂 ′ (q) and 𝑂 ′′(q) respectively. The right panels are the corresponding data for the HF ground state which corresponds to a spin-polarized KIVC with q … view at source ↗
read the original abstract

The understanding of quantum many-body states in twisted bilayer graphene at the magic angle has been greatly improved both in experiment and in theory. However, away from the exactly solvable chiral limit and the sign-problem-free charge neutrality point, the calculation of the ground state properties and the identification of the phase diagram are challenging due to the exponential increase in the complexity, which has rendered explanations of experimentally observed insulating and superconducting phases restricted largely to the perturbative level. Here we focus on the filling factors $\nu = \pm2$ away from charge neutrality and address the question of the strain dependence of the interacting ground state. We adjust our continuous field momentum-space quantum Monte Carlo (QMC) method to treat the sign problem approximately, and perform a quantum many-body study together with exact diagonalization (ED) and Hartree-Fock (HF) mean field. Leveraging this combined protocol of QMC, ED, and HF, we investigate the strain-tuned transition from the Kramers intervalley coherent (KIVC) state to the incommensurate Kekul\'e spiral state (IKS). Our computational protocol sheds light on the KIVC-IKS transition in a projected correlated flat-band setting, and opens the door for further understanding of the rich phase diagram of twisted bilayer graphene and other strongly-correlated flat-band systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript investigates the strain dependence of the interacting ground state in twisted bilayer graphene at filling factors ν=±2 using a combination of continuous-field momentum-space quantum Monte Carlo (QMC) with an approximate treatment of the sign problem, exact diagonalization (ED), and Hartree-Fock (HF). It reports a strain-tuned transition from the Kramers intervalley coherent (KIVC) state to the incommensurate Kekulé spiral (IKS) state in a projected correlated flat-band setting.

Significance. If the approximate QMC treatment is shown to be reliable, the work would offer a non-perturbative many-body perspective on the KIVC-IKS transition at ν=±2, complementing existing perturbative approaches and providing a protocol applicable to other flat-band systems. The multi-method cross-check is a positive feature, though its scope is limited by the unvalidated approximation.

major comments (2)
  1. [Abstract] Abstract (method paragraph): The claim that the adjusted continuous-field momentum-space QMC yields reliable ground-state properties and the KIVC-IKS phase boundary rests on an approximate treatment of the sign problem whose accuracy is not quantified; no error bars, convergence tests against exact limits, or benchmarks at the relevant strains and fillings are described, leaving the reported transition line vulnerable to uncontrolled bias.
  2. [Abstract] The combined protocol is said to investigate the strain-tuned transition, but ED and HF cross-checks are noted to apply only in limited regimes and do not directly certify the QMC approximation at ν=±2; this makes the central phase-boundary result dependent on an unvalidated step.
minor comments (1)
  1. [Abstract] The abstract refers to 'our continuous field momentum-space quantum Monte Carlo (QMC) method' without a citation to the prior work defining the base method or the specific adjustment.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments regarding the validation of our approximate QMC treatment. We address each point below and will revise the manuscript accordingly to provide a more balanced presentation.

read point-by-point responses

  1. Referee: [Abstract] Abstract (method paragraph): The claim that the adjusted continuous-field momentum-space QMC yields reliable ground-state properties and the KIVC-IKS phase boundary rests on an approximate treatment of the sign problem whose accuracy is not quantified; no error bars, convergence tests against exact limits, or benchmarks at the relevant strains and fillings are described, leaving the reported transition line vulnerable to uncontrolled bias.

    Authors: We agree that the current abstract does not sufficiently qualify the approximate nature of the sign-problem treatment or provide quantitative validation. In the revised manuscript we will modify the abstract to state the approximate character explicitly and add a dedicated paragraph (or appendix) reporting benchmarks against ED on small clusters where the sign problem is controllable, together with convergence tests with respect to the approximation parameter. These additions will allow readers to assess the expected accuracy of the reported transition line. revision: yes

  2. Referee: [Abstract] The combined protocol is said to investigate the strain-tuned transition, but ED and HF cross-checks are noted to apply only in limited regimes and do not directly certify the QMC approximation at ν=±2; this makes the central phase-boundary result dependent on an unvalidated step.

    Authors: The referee is correct that ED is restricted to small system sizes and HF is mean-field; neither directly validates the QMC approximation at the fillings and strains of primary interest. We will revise the abstract and main text to clarify that the central KIVC–IKS boundary is obtained from the adjusted QMC, with ED and HF serving only as consistency checks in overlapping regimes. We will also discuss the expected systematic bias of the sign-problem approximation and note that full error quantification remains an open technical challenge for this method. revision: partial

Circularity Check

0 steps flagged

No circularity: numerical protocol uses independent methods without self-referential reductions

full rationale

The paper presents a combined QMC+ED+HF computational study of the strain-tuned KIVC-to-IKS transition at ν=±2. The abstract describes adjusting the continuous-field momentum-space QMC to treat the sign problem approximately and then using the protocol to investigate the transition. No equations, fitted parameters, or results are shown to reduce by construction to their own inputs; the methods are externally cross-validated against each other in limited regimes and benchmarked as standard numerical tools. The derivation chain is therefore self-contained and does not exhibit any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract supplies insufficient technical detail to enumerate concrete free parameters, axioms, or invented entities; the approximate sign-problem treatment in QMC is the most obvious unstated modeling choice whose validity is not independently justified here.

pith-pipeline@v0.9.1-grok · 5796 in / 1203 out tokens · 35680 ms · 2026-06-30T17:38:44.645918+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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

Reference graph

Works this paper leans on

57 extracted references · 5 canonical work pages

  1. [1]

    Y.Cao,V.Fatemi,A.Demir,S.Fang,S.L.Tomarken,J.Y.Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxi- ras, R. C. Ashoori, and P. Jarillo-Herrero, Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature556, 80–84 (2018)

    2018

  2. [2]

    Y.Cao,V.Fatemi,S.Fang,K.Watanabe,T.Taniguchi,E.Kaxi- ras,andP.Jarillo-Herrero,Unconventionalsuperconductivityin magic-angle graphene superlattices, Nature556, 43–50 (2018)

    2018

  3. [3]

    Kerelsky, L

    A. Kerelsky, L. J. McGilly, D. M. Kennes, L. Xian, M. Yankowitz, S. Chen, K. Watanabe, T. Taniguchi, J. Hone, C. Dean,et al., Maximized electron interactions at the magic angle in twisted bilayer graphene, Nature572, 95 (2019)

    2019

  4. [4]

    Zhang, H

    Y.Choi,J.Kemmer,Y.Peng,A.Thomson,H.Arora,R.Polski, Y. Zhang, H. Ren, J. Alicea, G. Refael,et al., Electronic corre- lationsintwistedbilayergraphenenearthemagicangle,Nature physics15, 1174 (2019)

    2019

  5. [5]

    Y. Xie, B. Lian, B. Jäck, 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 bilayer graphene, Nature572, 101 (2019)

    2019

  6. [6]

    D. Wong, K. P. Nuckolls, M. Oh, B. Lian, Y. Xie, S. Jeon, K. Watanabe, T. Taniguchi, B. A. Bernevig, and A. Yazdani, Cascadeofelectronictransitionsinmagic-angletwistedbilayer graphene, Nature582, 198 (2020)

    2020

  7. [7]

    N. P. Kazmierczak, M. Van Winkle, C. Ophus, K. C. Bustillo, S.Carr,H.G.Brown,J.Ciston,T.Taniguchi,K.Watanabe,and D. K. Bediako, Strain fields in twisted bilayer graphene, Nature materials20, 956 (2021)

    2021

  8. [8]

    Huder, A

    L. Huder, A. Artaud, T. Le Quang, G. T. de Laissardière, A. G. M. Jansen, G. Lapertot, C. Chapelier, and V. T. Renard, Electronic spectrum of twisted graphene layers under heteros- train, Phys. Rev. Lett.120, 156405 (2018)

    2018

  9. [9]

    Mesple, A

    F. Mesple, A. Missaoui, T. Cea, L. Huder, F. Guinea, G. Tram- blydeLaissardière,C.Chapelier,andV.T.Renard,Heterostrain determines flat bands in magic-angle twisted graphene layers, Phys. Rev. Lett.127, 126405 (2021)

    2021

  10. [10]

    D.E.Parker,T.Soejima,J.Hauschild,M.P.Zaletel,andN.Bult- inck, Strain-induced quantum phase transitions in magic-angle graphene, Phys. Rev. Lett.127, 027601 (2021)

    2021

  11. [11]

    K. P. Nuckolls, R. L. Lee, M. Oh, D. Wong, T. Soejima, J. P. Hong, D. Călugăru, 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

  12. [12]

    H. Kim, Y. Choi, E. Lantagne-Hurtubise, C. Lewandowski, A.Thomson,L.Kong,H.Zhou,E.Baum,Y.Zhang,L.Holleis, K. Watanabe, T. Taniguchi, A. F. Young, J. Alicea, and S. Nadj-Perge, Imaging inter-valley coherent order in magic- angle twisted trilayer graphene, Nature623, 942 (2023)

    2023

  13. [13]

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

    2021

  14. [14]

    Wagner, Y

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

    2022

  15. [15]

    Bultinck, E

    N. Bultinck, E. Khalaf, S. Liu, S. Chatterjee, A. Vishwanath, andM.P.Zaletel,Groundstateandhiddensymmetryofmagic- angle graphene at even integer filling, Phys. Rev. X10, 031034 (2020)

    2020

  16. [16]

    B.Lian,Z.-D.Song,N.Regnault,D.K.Efetov,A.Yazdani,and B. A. Bernevig, Twisted bilayer graphene. iv. exact insulator ground states and phase diagram, Phys. Rev. B103, 205414 (2021)

    2021

  17. [17]

    Kang and O

    J. Kang and O. Vafek, Strong coupling phases of partially filled 12 twisted bilayer graphene narrow bands, Phys. Rev. Lett.122, 246401 (2019)

    2019

  18. [18]

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

    2021

  19. [19]

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

    2021

  20. [20]

    Herzog-Arbeitman, J

    J. Herzog-Arbeitman, J. Yu, D. Călugăru, H. Hu, N. Regnault, O. Vafek, J. Kang, and B. A. Bernevig, Topological heavy fermion model as an efficient representation of atomistic strain and relaxation in twisted bilayer graphene, Phys. Rev. B112, 125128 (2025)

    2025

  21. [21]

    Herzog-Arbeitman, D

    J. Herzog-Arbeitman, D. Călugăru, H. Hu, J. Yu, N. Regnault, J.Kang,B.A.Bernevig,andO.Vafek,Kekuléspiralorderfrom strained topological heavy fermions, Phys. Rev. B112, 125129 (2025)

    2025

  22. [22]

    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

  23. [23]

    J. Ahn, S. Park, and B.-J. Yang, Failure of nielsen-ninomiya theorem and fragile topology in two-dimensional systems with space-time inversion symmetry: Application to twisted bilayer graphene at magic angle, Phys. Rev. X9, 021013 (2019)

    2019

  24. [24]

    Y. H. Kwan, Z. Wang, G. Wagner, S. H. Simon, S. A. Parameswaran, and N. Bultinck, Textured exciton insulators, Phys. Rev. B112, 035129 (2025)

    2025

  25. [25]

    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

  26. [26]

    Classen, and Z

    C.Huang,N.Parthenios,M.Ulybyshev,X.Zhang,F.F.Assaad, L. Classen, and Z. Y. Meng, Angle-tuned gross-neveu quantum criticality in twisted bilayer graphene, Nature Communications 16, 7176 (2025)

    2025

  27. [27]

    X.Zhang, G.Pan,Y.Zhang,J.Kang,andZ.Y.Meng,Momen- tum space quantum monte carlo on twisted bilayer graphene, Chinese Physics Letters38, 077305 (2021)

    2021

  28. [28]

    J. S. Hofmann, E. Khalaf, A. Vishwanath, E. Berg, and J. Y. Lee,Fermionicmontecarlostudyofarealisticmodeloftwisted bilayer graphene, Phys. Rev. X12, 011061 (2022)

    2022

  29. [29]

    Zhang, G

    X. Zhang, G. Pan, X. Y. Xu, and Z. Y. Meng, Fermion sign boundstheoryinquantummontecarlosimulation,Phys.Rev.B 106, 035121 (2022)

    2022

  30. [30]

    Zhang, G

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

    2023

  31. [31]

    Huang, X

    C. Huang, X. Zhang, G. Pan, H. Li, K. Sun, X. Dai, and Z. Y. Meng, Evolution from quantum anomalous hall insula- tor to heavy-fermion semimetal in magic-angle twisted bilayer graphene, Phys. Rev. B109, 125404 (2024)

    2024

  32. [32]

    Pan and Z

    G. Pan and Z. Y. Meng, The sign prob- lem in quantum monte carlo simulations, in Encyclopedia of Condensed Matter Physics (Second Edition), edited by T. Chakraborty (Academic Press, Oxford, 2024) second edition ed., pp. 879–893

    2024

  33. [34]

    S.Carr,S.Fang,Z.Zhu,andE.Kaxiras,Exactcontinuummodel for low-energy electronic states of twisted bilayer graphene, Phys. Rev. Res.1, 013001 (2019)

    2019

  34. [35]

    Koshino, N

    M. Koshino, N. F. Q. Yuan, T. Koretsune, M. Ochi, K. Kuroki, and L. Fu, Maximally localized wannier orbitals and the ex- tended hubbard model for twisted bilayer graphene, Phys. Rev. X8, 031087 (2018)

    2018

  35. [36]

    N. N. T. Nam and M. Koshino, Lattice relaxation and energy band modulation in twisted bilayer graphene, Phys. Rev. B96, 075311 (2017)

    2017

  36. [37]

    Vafek and J

    O. Vafek and J. Kang, Renormalization group study of hidden symmetry in twisted bilayer graphene with coulomb interac- tions, Phys. Rev. Lett.125, 257602 (2020)

    2020

  37. [38]

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

    2019

  38. [39]

    The free dispersions under strain, the comparison of structure factors between QMC and HF, the HF energy as a function of qIKS, single-particle Green’s function in QMC, and the proper- ties of the approximated QMC are given in this Supplemental Information

  39. [40]

    F. Xie, A. Cowsik, Z.-D. Song, B. Lian, B. A. Bernevig, and N.Regnault,Twistedbilayergraphene.vi.anexactdiagonaliza- tion study at nonzero integer filling, Phys. Rev. B103, 205416 (2021)

    2021

  40. [41]

    P.Potasz,M.Xie,andA.H.MacDonald,Exactdiagonalization for magic-angle twisted bilayer graphene, Phys. Rev. Lett.127, 147203 (2021)

    2021

  41. [42]

    Y. H. Kwan, Z. Wang, G. Wagner, N. Bultinck, S. H. Simon, and S. A. Parameswaran, Mean-field modeling of moiré materials: a user’s guide with selected applications to twisted bilayer graphene, Advances in Physics0, 1 (2025), https://doi.org/10.1080/00018732.2025.2600658

    work page doi:10.1080/00018732.2025.2600658 2025

  42. [43]

    Z.Wang,G.Wagner,Y.H.Kwan,N.Bultinck,S.H.Simon,and S.A.Parameswaran,Puttinganewspinontheincommensurate kekulé spiral: from spin-valley locking and collective modes to fermiology and implications for superconductivity (2025), arXiv:2509.12320 [cond-mat.str-el]

    work page arXiv 2025

  43. [44]

    Z. Wang, Y. H. Kwan, G. Wagner, S. H. Simon, N. Bultinck, andS.A.Parameswaran,Chern-texturedexcitoninsulatorswith valleyspiralorderinmoirématerials,Phys.Rev.B112,035130 (2025)

    2025

  44. [45]

    F.GuineaandN.R.Walet,Electrostaticeffects,banddistortions, andsuperconductivityintwistedgraphenebilayers,Proceedings of the National Academy of Sciences115, 13174 (2018)

    2018

  45. [46]

    L.Rademaker,D.A.Abanin,andP.Mellado,Chargesmoothen- ing and band flattening due to hartree corrections in twisted bilayer graphene, Phys. Rev. B100, 205114 (2019)

    2019

  46. [47]

    T. Cea, N. R. Walet, and F. Guinea, Electronic band structure and pinning of fermi energy to van hove singularities in twisted bilayergraphene: Aself-consistentapproach,Phys.Rev.B100, 205113 (2019)

    2019

  47. [48]

    Z.A.H.Goodwin,V.Vitale,X.Liang,A.A.Mostofi,andJ.Lis- chner, Hartree theory calculations of quasiparticle properties in twisted bilayer graphene (2020), arXiv:2004.14784 [cond- mat.mes-hall]

    work page arXiv 2020

  48. [49]

    J. Kang, B. A. Bernevig, and O. Vafek, Cascades between light and heavy fermions in the normal state of magic angle twisted bilayer graphene (2021), arXiv:2104.01145 [cond-mat.str-el]

    work page arXiv 2021

  49. [50]

    X.Zhang,K.Sun,H.Li,G.Pan,andZ.Y.Meng,Superconduc- tivity and bosonic fluid emerging from moiré flat bands, Phys. Rev. B106, 184517 (2022)

    2022

  50. [51]

    G. Pan, X. Zhang, H. Lu, H. Li, B.-B. Chen, K. Sun, and Z. Y. Meng, Thermodynamic characteristic for a correlated flat-band system with a quantum anomalous hall ground state, Phys. Rev. Lett.130, 016401 (2023)

    2023

  51. [52]

    X.Lin,B.-B.Chen,W.Li,Z.Y.Meng,andT.Shi,Excitonpro- liferation and fate of the topological mott insulator in a twisted 13 bilayer graphene lattice model, Phys. Rev. Lett.128, 157201 (2022)

    2022

  52. [53]

    Bistritzer and A

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

    2011

  53. [54]

    HPC2021,InformationTechnologyServices,TheUniversityof Hong Kong

  54. [55]

    Beijing PARATERA Tech CO.,Ltd

  55. [56]

    Grover, Entanglement of interacting fermions in quantum monte carlo calculations, Phys

    T. Grover, Entanglement of interacting fermions in quantum monte carlo calculations, Phys. Rev. Lett.111, 130402 (2013)

    2013

  56. [57]

    Strain-Tuned Incommensurate Kekulé Spiral (IKS) Order in T wisted Bilayer Graphene: a Quantum Many-Body Investigation

    M. Ulybyshev and F. F. Assaad, Beyond the instanton gas approach: dominant thimbles approximation for the Hubbard model, arXiv e-prints , arXiv:2407.09452 (2024), arXiv:2407.09452 [cond-mat.str-el]. Acknowledgments We acknowledge discussions with Nikolaos Parthenios and Jeyong Park on similar topics. C.H. and Z.Y.M. acknowledge thesupportfromtheResearchGr...

    work page arXiv 2024

  57. [58]

    Different dispersion colors are used to indicate distinct bands. S2. INTERV ALLEY COHERENCE STRUCTURE FACTORS In this section we show the comparisons of the KIVC structure factor𝑆KIVC and IKS structure factor𝑆 IKS from both approximated QMC and HF at𝜖s =0 and𝜖 s =0.6% in Supplementary Fig. S2. The definitions of these structure factors are given in the ma...