REVIEW 2 major objections 1 minor 1 cited by
At filling factors ±2, strain drives a transition from the Kramers intervalley coherent state to the incommensurate Kekulé spiral state in twisted bilayer graphene.
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 →
T0 review · grok-4.3
2026-06-30 17:38 UTC pith:D54Z6FRE
load-bearing objection The reported strain-tuned KIVC-to-IKS transition at ν=±2 rests on an approximate sign-problem fix in QMC whose reliability is not demonstrated. the 2 major comments →
Strain-Tuned Incommensurate Kekul\'e Spiral Order in Twisted Bilayer Graphene: a Quantum Many-Body Study
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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.
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.
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
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [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.
- [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)
- [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
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
-
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
-
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
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
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.
Figures
Forward citations
Cited by 1 Pith paper
-
Trion Excitations in Twisted Bilayer Graphene: A Quantum Monte Carlo Study
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
-
[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]
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]
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]
-
[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]
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]
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]
- [9]
-
[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]
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]
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]
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]
-
[15]
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]
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]
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]
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]
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]
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]
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]
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]
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]
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]
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]
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]
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]
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]
- [30]
- [31]
- [32]
-
[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
-
[35]
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
-
[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
-
[37]
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
-
[38]
Z. Bi, N. F. Q. Yuan, and L. Fu, Designing flat bands by strain, Phys. Rev. B100, 035448 (2019)
work page 2019
-
[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
-
[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
-
[41]
P.Potasz,M.Xie,andA.H.MacDonald,Exactdiagonalization for magic-angle twisted bilayer graphene, Phys. Rev. Lett.127, 147203 (2021)
work page 2021
-
[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
- [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
-
[45]
F.GuineaandN.R.Walet,Electrostaticeffects,banddistortions, andsuperconductivityintwistedgraphenebilayers,Proceedings of the National Academy of Sciences115, 13174 (2018)
work page 2018
-
[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
-
[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
- [48]
- [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
-
[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
-
[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
-
[53]
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
-
[54]
HPC2021,InformationTechnologyServices,TheUniversityof Hong Kong
-
[55]
Beijing PARATERA Tech CO.,Ltd
-
[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
-
[57]
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...
-
[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...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.