pith. sign in

arxiv: 2605.20677 · v1 · 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-05-21 04:11 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords twisted bilayer graphenestrain tuningKekulé spiral orderintervalley coherencequantum Monte Carloflat band modelmany-body physicsmoiré materials
0
0 comments X

The pith

Strain drives a transition from Kramers intervalley coherent to incommensurate Kekulé spiral order 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 seeks to map how mechanical strain changes the ground state of interacting electrons in twisted bilayer graphene when the system is doped to two electrons or holes per moiré unit cell. It combines adjusted quantum Monte Carlo simulations that approximately manage the sign problem, exact diagonalization on small clusters, and Hartree-Fock calculations to compare candidate ordered states. The results indicate that zero or low strain favors a time-reversal symmetric state with intervalley coherence, while moderate strain stabilizes a spiral Kekulé pattern that breaks additional symmetries. A sympathetic reader would care because this offers a non-perturbative window into phases that experiments associate with insulation and possible superconductivity, and the method can be reused for other flat-band systems where direct simulation is blocked by complexity.

Core claim

In the projected correlated flat-band model for twisted bilayer graphene at fillings ν = ±2, the interacting ground state undergoes a strain-tuned transition from the Kramers intervalley coherent (KIVC) state to the incommensurate Kekulé spiral (IKS) state, as revealed by a combined computational study using adjusted continuous-field momentum-space quantum Monte Carlo, exact diagonalization, and Hartree-Fock methods.

What carries the argument

The strain-tuned transition from the Kramers intervalley coherent (KIVC) state to the incommensurate Kekulé spiral (IKS) state, identified through quantum many-body calculations in the flat-band projection.

If this is right

  • The KIVC-to-IKS transition can be tested by applying controlled uniaxial or biaxial strain in transport or spectroscopy experiments on magic-angle samples.
  • The IKS phase may exhibit distinct low-energy excitations or broken-symmetry signatures compared with the KIVC phase.
  • The combined QMC-ED-HF protocol supplies a benchmark for ground-state properties at fillings away from charge neutrality.
  • Similar strain-driven switches may appear in other moiré flat-band platforms once the same computational approach is applied.

Where Pith is reading between the lines

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

  • If the IKS state is robust, tuning through the transition region could be used to search for nearby superconducting domes or metallic phases.
  • The method opens a route to study how the same orders compete when additional perturbations such as electric fields or heterostrain are added.
  • Connection to real-space imaging or local probes could directly visualize the spiral modulation predicted for the IKS state.

Load-bearing premise

Adjusting the continuous-field momentum-space quantum Monte Carlo method to approximately handle the sign problem still produces accurate ground-state properties for the flat-band model away from the charge-neutrality point.

What would settle it

A calculation or measurement showing that the KIVC state remains stable at all accessible strain values, or that the IKS order fails to appear once the sign-problem approximation is removed, would falsify the reported 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

1 major / 1 minor

Summary. The paper investigates the strain dependence of the ground state in twisted bilayer graphene at fillings ν=±2 away from charge neutrality. In a projected correlated flat-band model, the authors combine an adjusted continuous-field momentum-space quantum Monte Carlo (QMC) method that approximately treats the sign problem, exact diagonalization (ED), and Hartree-Fock (HF) calculations to study the transition from the Kramers intervalley coherent (KIVC) state to the incommensurate Kekulé spiral (IKS) state as strain is varied.

Significance. If the approximate QMC results prove reliable, the work offers a multi-method numerical window into the strain-tuned phase diagram of magic-angle TBG at ν=±2, where perturbative or mean-field approaches have been dominant. The combined QMC+ED+HF protocol and focus on a tunable parameter (strain) are strengths that could help interpret experimental insulating states and motivate further studies of flat-band systems.

major comments (1)
  1. [Abstract and QMC Methods] Abstract and QMC protocol description: The central claim of a strain-tuned KIVC-to-IKS transition rests on ground-state order parameters and energies obtained from the adjusted continuous-field momentum-space QMC at ν=±2. The manuscript states that the method is modified to treat the sign problem approximately but provides no explicit description of the approximation (e.g., any tunable parameter or bias introduced), no quantitative error estimates, and no benchmarks such as recovery of known results at the sign-problem-free charge-neutrality point or direct comparison of order-parameter magnitudes with ED on accessible clusters. Without these controls, systematic bias could affect the apparent location or even the identity of the transition.
minor comments (1)
  1. [Abstract] The abstract and introduction would benefit from a brief statement clarifying the range of strain values explored and the precise definition of the incommensurate wavevector in the IKS state.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for the positive assessment of its potential significance. We address the major comment below and will revise the manuscript to improve clarity on the QMC protocol.

read point-by-point responses

  1. Referee: [Abstract and QMC Methods] Abstract and QMC protocol description: The central claim of a strain-tuned KIVC-to-IKS transition rests on ground-state order parameters and energies obtained from the adjusted continuous-field momentum-space QMC at ν=±2. The manuscript states that the method is modified to treat the sign problem approximately but provides no explicit description of the approximation (e.g., any tunable parameter or bias introduced), no quantitative error estimates, and no benchmarks such as recovery of known results at the sign-problem-free charge-neutrality point or direct comparison of order-parameter magnitudes with ED on accessible clusters. Without these controls, systematic bias could affect the apparent location or even the identity of the transition.

    Authors: We thank the referee for this constructive comment. The manuscript refers to our prior work for the technical details of the continuous-field momentum-space QMC, but we agree that a self-contained description of the sign-problem approximation is needed here. In the revised manuscript we will expand the Methods section to explicitly describe the approximation (including the form of the bias and any tunable parameters), provide quantitative error estimates from the auxiliary-field sampling, and add benchmarks: recovery of the known sign-problem-free results at charge neutrality (ν=0) together with direct order-parameter comparisons against ED on the smallest accessible clusters. These additions will allow readers to assess possible systematic effects on the reported KIVC–IKS transition. revision: yes

Circularity Check

0 steps flagged

No significant circularity: results generated from Hamiltonian via numerical methods

full rationale

The paper applies continuous-field momentum-space QMC (with approximate sign-problem handling), ED, and HF to the projected correlated flat-band model at ν=±2 under strain. Ground-state order parameters and the KIVC-IKS transition are obtained by direct simulation of the many-body Hamiltonian rather than by any self-definitional mapping, fitted parameter renamed as prediction, or load-bearing self-citation chain that reduces the reported phase boundary to an input. The protocol is self-contained against external benchmarks because the outputs remain falsifiable through cluster-size scaling, method cross-checks, and comparison to known limits at charge neutrality. No quoted step equates a claimed result to its own definition or prior fit.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The study assumes a projected flat-band model and an approximate sign-problem treatment whose validity is taken as given for the reported transition.

free parameters (1)
  • strain magnitude
    Varied continuously to locate the KIVC-IKS transition point.
axioms (1)
  • domain assumption The projected correlated flat-band Hamiltonian captures the essential low-energy physics away from the chiral limit.
    Invoked when restricting the calculation to the flat-band subspace.

pith-pipeline@v0.9.0 · 5796 in / 1278 out tokens · 40413 ms · 2026-05-21T04:11:01.621150+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We adjust our continuous field momentum-space quantum Monte Carlo (QMC) method to treat the sign problem approximately... 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é spiral state (IKS).

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We make an approximation by taking the absolute value of the original sampling weight as the new sampling weight... we evaluate ⟨Ô⟩ according to ∫ dC O_C W_C / Re(W_C) |Re(W_C)| / ∫ dC |Re(W_C)|.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

57 extracted references · 57 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 2021

  28. [28]

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

    work page 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)

    work page 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)

    work page 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)

    work page 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

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 2020

  37. [38]

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

    work page 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

    work page

  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)

    work page 2021

  40. [41]

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

    work page 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)

    work page 2025

  44. [45]

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

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 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)

    work page 2011

  53. [54]

    HPC2021,InformationTechnologyServices,TheUniversityof Hong Kong

    work page

  54. [55]

    Beijing PARATERA Tech CO.,Ltd

    work page

  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)

    work page 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...

    work page