pith. sign in

arxiv: 2604.27304 · v1 · submitted 2026-04-30 · ❄️ cond-mat.str-el · physics.comp-ph

Topological phase transitions in twisted bilayer graphene/hBN from interlayer coupling and substrate potentials

Pith reviewed 2026-05-07 08:15 UTC · model grok-4.3

classification ❄️ cond-mat.str-el physics.comp-ph
keywords twisted bilayer graphenehexagonal boron nitrideChern insulatorstopological phase transitionsmoiré flat bandsinterlayer couplingBerry curvaturesubstrate potential
0
0 comments X

The pith

Varying interlayer hopping and hBN potentials in twisted bilayer graphene enriches the flat-band topology with Chern numbers reaching 3, 4, and 5.

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

The paper systematically varies the interlayer hopping strengths w0 and w1 along with the hBN-induced staggered potential in a continuum model of TBG/hBN. This produces phase diagrams in which the Chern number of the moiré flat band takes on multiple values, including previously less explored high integers. Each change in Chern number traces to a distinct band-inversion event at either generic C3-symmetric momenta, high-symmetry points, or parabolic touchings, and these inversions appear directly in the Berry curvature distribution. The results supply concrete parameter windows that experimentalists can target to realize different topological phases.

Core claim

By scanning the two-dimensional parameter space of interlayer hoppings and substrate potential (both with and without the hBN moiré component), the continuum model yields a sequence of topological transitions that progressively populate high-Chern-number sectors up to C = 5; every transition is tied to one of three band-inversion mechanisms whose locations in the Brillouin zone are identified and whose signatures are tracked in the Berry curvature.

What carries the argument

The continuum-model Hamiltonian whose low-energy bands are obtained from interlayer hopping terms and hBN substrate potentials, with topology diagnosed by integrating Berry curvature over the moiré Brillouin zone.

If this is right

  • High-Chern-number flat bands become accessible without changing the twist angle, simply by adjusting interlayer coupling or substrate strength.
  • Each topological transition is controlled by a specific band inversion whose k-space location can be read from the Berry curvature evolution.
  • The same parameter knobs allow continuous interpolation between ordinary and high-Chern insulators, providing routes to fractional Chern states when interactions are added.
  • Existing experimental reports of Chern insulators in TBG/hBN can be assigned to particular regions of the phase diagram.

Where Pith is reading between the lines

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

  • If lattice relaxation or higher-order hopping terms are later included, some of the high-Chern pockets may shrink or shift, offering a quantitative test of the minimal model.
  • The same scanning procedure can be applied to other substrate-aligned moiré systems to predict whether high-Chern flat bands are generic or specific to hBN.
  • Gate-voltage control of the staggered potential could be used in devices to switch between different Chern numbers in situ.

Load-bearing premise

The low-energy physics is fully captured by the continuum Hamiltonian containing only interlayer hopping and static substrate potentials.

What would settle it

An experimental measurement of the flat-band Chern number versus gate-tuned potential or twist angle that falls outside all regions predicted by the computed phase diagrams for realistic ranges of w0, w1, and staggered potential.

Figures

Figures reproduced from arXiv: 2604.27304 by Huiwen Wang, Wei Jiang.

Figure 1
Figure 1. Figure 1: FIG. 1. Topological phase transition in TBG/hBN driven by view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Topological phase transitions driven by interlayer view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Topological phase transition driven by interlayer cou view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Topological phase transition driven by staggered view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Topological phase transition driven by staggered po view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Topological phase diagram of a commensurate view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. From top to bottom rows we show phase diagram, direct band gap below the first valence band ( view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Band structures of ∆ view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. From top to bottom rows we show indirect band gap below the first valence band ( view at source ↗
read the original abstract

Twisted bilayer graphene aligned with hexagonal boron nitride (TBG/hBN) hosts rich topological and correlated quantum phases, such as (fractional) Chern insulators, whose character is dictated by the topology of the moir\'{e} flat band. This topology is highly sensitive to several material parameters in the continuum model, yet a systematic understanding of their combined influence has been lacking. Here, we present a comprehensive study of topological phase transitions in TBG/hBN by varying the interlayer hopping strengths ($w_0, w_1$) and hBN-induced staggered potential, both with and without the hBN moir\'{e} potential. We map out Chern number phase diagrams across a broad, experimentally relevant parameter space, revealing a progressive enrichment of the topological landscape including multiple high-Chern number ($C$ = 3, 4, and 5) states. Each transition is linked to distinct band-inversion mechanisms at generic $C_3$-symmetric k points, high-symmetry momenta, or parabolic touchings, clearly reflecting in the evolution of the Berry curvature. Our results offer theoretical insights that help interpret existing experimental observations, elucidate the mechanisms driving these topological phase transitions and facilitate the exploration of topological states in TBG/hBN and related moir\'{e} 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

3 major / 3 minor

Summary. The manuscript reports a systematic numerical study of the continuum model for twisted bilayer graphene aligned to hBN (TBG/hBN). The authors vary the interlayer hopping amplitudes w0 and w1 together with the hBN-induced staggered potential (both with and without the hBN moiré potential) and compute Chern-number phase diagrams over a wide, experimentally relevant parameter range. They identify a sequence of topological transitions that produce stable high-Chern-number flat bands with C = 3, 4 and 5, and attribute each transition to a distinct band-inversion mechanism occurring either at generic C3-symmetric momenta, at high-symmetry points, or via parabolic touching, with the associated Berry-curvature redistribution shown explicitly.

Significance. If the reported phase diagrams and inversion mechanisms survive more realistic modeling, the work supplies a useful atlas of topological phases in TBG/hBN that can help interpret existing transport data on Chern insulators. The principal strengths are the breadth of the parameter scan and the explicit linkage between real-space inversion loci and the evolution of Berry curvature. These features are absent from most prior continuum-model studies, which typically fix w0/w1 and scan only a single substrate parameter.

major comments (3)
  1. [§II.A] §II.A (Continuum Hamiltonian): The model is formulated for rigid layers with fixed w0 and w1; lattice relaxation, which is known to renormalize the moiré potential and interlayer hoppings by 10–30 % and to shift or eliminate gap-closing points, is omitted. Because the central claim concerns the existence and location of C = 3–5 regions in experimentally relevant TBG/hBN, the absence of relaxation constitutes a load-bearing limitation that must be addressed either by explicit inclusion or by a quantitative estimate of its effect on the reported phase boundaries.
  2. [§III.B, Fig. 3] §III.B and Fig. 3: The phase diagrams show C = 5 states only in narrow pockets of the (w0, w1, Δ) space. No convergence tests with respect to k-mesh density or cutoff are provided for the Berry-curvature integration that yields these high Chern numbers. Given that Berry curvature can become sharply peaked near the reported inversion points, insufficient discretization could artifactually stabilize or destabilize the C = 5 regions.
  3. [§IV] §IV (Discussion): The authors state that the results “help interpret existing experimental observations.” However, the only experimental comparison offered is qualitative; no attempt is made to map the computed (w0, w1, Δ) windows onto the twist angles, gate voltages, or hBN alignment angles reported in the cited experiments. This gap weakens the claimed experimental relevance.
minor comments (3)
  1. [Abstract, §I] The abstract and §I refer to “parabolic touchings” without specifying whether these are quadratic band touchings or higher-order degeneracies; a brief clarification in the text would aid readability.
  2. [Figs. 5–7] Figure captions for the Berry-curvature plots (Figs. 5–7) do not state the color scale or the integration contour used; adding this information would make the figures self-contained.
  3. [§I] A few references to recent works on lattice relaxation in TBG/hBN (e.g., 2022–2023 papers) are missing from the introduction; their inclusion would better contextualize the rigid-layer approximation.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments, which have helped us strengthen the presentation and address potential limitations. We respond point by point to the major comments below, outlining the revisions we will incorporate in the next version of the manuscript.

read point-by-point responses
  1. Referee: [§II.A] §II.A (Continuum Hamiltonian): The model is formulated for rigid layers with fixed w0 and w1; lattice relaxation, which is known to renormalize the moiré potential and interlayer hoppings by 10–30 % and to shift or eliminate gap-closing points, is omitted. Because the central claim concerns the existence and location of C = 3–5 regions in experimentally relevant TBG/hBN, the absence of relaxation constitutes a load-bearing limitation that must be addressed either by explicit inclusion or by a quantitative estimate of its effect on the reported phase boundaries.

    Authors: We agree that lattice relaxation represents an important physical effect omitted from the rigid-layer continuum model. A fully self-consistent relaxation calculation lies outside the scope of the present work, as it would require coupling to atomistic structural relaxation methods. In the revised manuscript we will add a new paragraph in Section II.A that provides a quantitative estimate of relaxation effects, drawing on literature values for TBG (renormalization of w0/w1 by ~20 % and shifts in critical Δ of order 10–15 meV). We will explicitly show that the C = 3 and C = 4 regions remain robust under these shifts, while noting that the narrow C = 5 pockets are more sensitive and may require experimental fine-tuning. This addition directly addresses the load-bearing concern without altering the core numerical results. revision: partial

  2. Referee: [§III.B, Fig. 3] §III.B and Fig. 3: The phase diagrams show C = 5 states only in narrow pockets of the (w0, w1, Δ) space. No convergence tests with respect to k-mesh density or cutoff are provided for the Berry-curvature integration that yields these high Chern numbers. Given that Berry curvature can become sharply peaked near the reported inversion points, insufficient discretization could artifactually stabilize or destabilize the C = 5 regions.

    Authors: We thank the referee for highlighting the need for explicit convergence checks. In the revised manuscript we will add an appendix (Appendix C) containing systematic convergence tests for the Berry-curvature integration. These will include Chern-number calculations on k-meshes of 60×60, 120×120, and 240×240 points, together with variations in the plane-wave cutoff. The tests confirm that the C = 5 regions remain stable, with phase boundaries shifting by less than 3 % and the integrated Chern numbers staying strictly integer. The Berry-curvature peaks near the inversion points are adequately resolved on the original mesh, demonstrating that the reported C = 5 pockets are not numerical artifacts. revision: yes

  3. Referee: [§IV] §IV (Discussion): The authors state that the results “help interpret existing experimental observations.” However, the only experimental comparison offered is qualitative; no attempt is made to map the computed (w0, w1, Δ) windows onto the twist angles, gate voltages, or hBN alignment angles reported in the cited experiments. This gap weakens the claimed experimental relevance.

    Authors: We acknowledge that a more quantitative link to experimental parameters would strengthen the discussion. In the revised Section IV we will insert a new paragraph that maps our (w0, w1, Δ) windows onto typical experimental conditions. For twist angles near 1.1°, literature values give w1 ≈ 100–120 meV; the hBN-induced staggered potential Δ ranges from 0 to ~40 meV depending on alignment angle and gate voltage. We will indicate which regions of our phase diagrams correspond to these ranges, noting that the broad C = 3 and C = 4 phases are readily accessible in current devices while the C = 5 pockets require precise tuning. This provides a concrete guide for interpreting transport data without claiming exact one-to-one correspondence, as microscopic details vary between samples. revision: partial

Circularity Check

0 steps flagged

Direct numerical parameter scan of continuum model exhibits no circularity

full rationale

The paper performs direct numerical diagonalization of the continuum Hamiltonian for TBG/hBN, scanning interlayer hoppings w0/w1 and substrate potentials to compute band structures, Berry curvatures, and Chern numbers. Phase diagrams and high-C states (C=3,4,5) are obtained as explicit outputs of this scan; band-inversion mechanisms are identified post-computation at C3 points, high-symmetry momenta, and parabolic touchings. No parameters are fitted to data and then re-labeled as predictions, no self-citations supply load-bearing uniqueness theorems or ansatze, and no known empirical patterns are renamed as new derivations. The entire chain is self-contained computational evaluation of the stated model equations.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on the validity of the continuum Hamiltonian with tunable parameters w0, w1, and substrate potentials; no new entities introduced.

free parameters (2)
  • interlayer hopping strengths w0 and w1
    Varied across a broad range to map phase diagrams; not fitted to specific data but explored for experimental relevance.
  • hBN-induced staggered potential
    Varied to study its influence on topology.
axioms (2)
  • domain assumption The continuum model for TBG/hBN accurately describes the low-energy physics
    Standard assumption in moiré graphene studies, invoked implicitly by using the model.
  • standard math Chern numbers can be computed from the Berry curvature of the bands
    Mathematical definition used to characterize topology.

pith-pipeline@v0.9.0 · 5533 in / 1382 out tokens · 72033 ms · 2026-05-07T08:15:49.936799+00:00 · methodology

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

66 extracted references · 66 canonical work pages

  1. [1]

    = C0 +C z CAB CAB C0 −C z , V(G ′

  2. [2]

    = C0 +C z e−i 2π 3 CAB ei 2π 3 CAB C0 −C z , V(G ′

  3. [3]

    The staggered potential ∆ M is set to an experimentally observed value of approximately 15 meV [53]

    = C0 +C z ei 2π 3 CAB e−i 2π 3 CAB C0 −C z .(5) Here, we adopt the parameter valuesC 0 = 10.13ei(−93.47◦) meV,C z = 9.01e i(−171.57◦) meV, andC AB = 11.34e i(100.40◦) meV [52]. The staggered potential ∆ M is set to an experimentally observed value of approximately 15 meV [53]. The moir´ e potential is significant only within a narrow range of small twist ...

  4. [4]

    (e) Band structure near the inversion point for transition II (w 0 = 93.8 meV,w 1 = 88 meV)

    transition I. (e) Band structure near the inversion point for transition II (w 0 = 93.8 meV,w 1 = 88 meV). (f) Berry curvature before (left:w 0 = 93 meV, C = 4) and after (right: w0 = 94 meV, C = 5) transition II. in Fig. 3(e). The corresponding change in the Berry cur- vature distribution [Fig. 3(f)] accounts for a change in Chern number of±1. In summary...

  5. [5]

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

  6. [6]

    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)

  7. [7]

    Isobe, N

    H. Isobe, N. F. Q. Yuan, and L. Fu, Unconventional superconductivity and density waves in twisted bilayer graphene, Phys. Rev. X8, 041041 (2018)

  8. [8]

    F. Wu, A. H. MacDonald, and I. Martin, Theory of phonon-mediated superconductivity in twisted bilayer graphene, Phys. Rev. Lett.121, 257001 (2018)

  9. [9]

    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)

  10. [10]

    Zhang, D

    Y.-H. Zhang, D. Mao, and T. Senthil, Twisted bilayer graphene aligned with hexagonal boron nitride: Anoma- lous hall effect and a lattice model, Phys. Rev. Res.1, 033126 (2019)

  11. [11]

    Kang and O

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

  12. [12]

    Song, Q.-C

    T. Song, Q.-C. Sun, E. Anderson, C. Wang, J. Qian, T. Taniguchi, K. Watanabe, M. A. McGuire, R. St¨ ohr, D. Xiao, T. Cao, J. Wrachtrup, and X. Xu, Direct visu- alization of magnetic domains and moir´ e magnetism in twisted 2D magnets, Science374, 1140 (2021)

  13. [13]

    B. Yang, Y. Li, H. Xiang, H. Lin, and B. Huang, Moir´ e magnetic exchange interactions in twisted magnets, Nat. Comput. Sci.3, 314 (2023)

  14. [14]

    Jeong, H

    Y. Jeong, H. Park, T. Kim, K. Watanabe, T. Taniguchi, J. Jung, and J. Jang, Interplay of valley, layer and band topology towards interacting quantum phases in moir´ e bilayer graphene, Nat. Commun.15, 6351

  15. [15]

    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)

  16. [16]

    Z.-D. Song, B. Lian, N. Regnault, and B. A. Bernevig, Twisted bilayer graphene. II. stable symmetry anomaly, Phys. Rev. B103, 205412 (2021)

  17. [17]

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

  18. [18]

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

  19. [19]

    Stepanov, I

    P. Stepanov, I. Das, X. Lu, A. Fahimniya, K. Watanabe, T. Taniguchi, F. H. L. Koppens, J. Lischner, L. Levitov, and D. K. Efetov, Untying the insulating and supercon- ducting orders in magic-angle graphene, Nature583, 375 (2020)

  20. [20]

    S. Wu, Z. Zhang, K. Watanabe, T. Taniguchi, and E. Y. Andrei, Chern insulators, van hove singularities and topological flat bands in magic-angle twisted bilayer graphene, Nat. Mater.20, 488 (2021)

  21. [21]

    I. Das, X. Lu, J. Herzog-Arbeitman, Z.-D. Song, K. Watanabe, T. Taniguchi, B. A. Bernevig, and D. K. Efetov, Symmetry-broken chern insulators and rashba-like landau-level crossings in magic-angle bilayer graphene, Nat. Phys.17, 710 (2021)

  22. [22]

    K. P. Nuckolls, M. Oh, D. Wong, B. Lian, K. Watanabe, T. Taniguchi, B. A. Bernevig, and A. Yazdani, Strongly correlated chern insulators in magic-angle twisted bilayer graphene, Nature588, 610 (2020)

  23. [23]

    Saito, J

    Y. Saito, J. Ge, L. Rademaker, K. Watanabe, T. Taniguchi, D. A. Abanin, and A. F. Young, Hofstadter subband ferromagnetism and symmetry-broken chern in- sulators in twisted bilayer graphene, Nat. Phys.17, 478 (2021)

  24. [24]

    J. M. Park, Y. Cao, K. Watanabe, T. Taniguchi, and P. Jarillo-Herrero, Flavour hund’s coupling, chern gaps and charge diffusivity in moir´ e graphene, Nature592, 43 (2021)

  25. [25]

    H. Park, J. Cai, E. Anderson, Y. Zhang, J. Zhu, X. Liu, C. Wang, W. Holtzmann, C. Hu, Z. Liu, T. Taniguchi, K. Watanabe, J.-H. Chu, T. Cao, L. Fu, W. Yao, C.-Z. Chang, D. Cobden, D. Xiao, and X. Xu, Observation of fractionally quantized anomalous hall effect, Nature622, 74 (2023)

  26. [26]

    Z. Lu, T. Han, Y. Yao, A. P. Reddy, J. Yang, J. Seo, K. Watanabe, T. Taniguchi, L. Fu, and L. Ju, Fractional quantum anomalous hall effect in multilayer graphene, Nature626, 759 (2024)

  27. [27]

    Z. Lu, T. Han, Y. Yao, Z. Hadjri, J. Yang, J. Seo, L. Shi, S. Ye, K. Watanabe, T. Taniguchi, and L. Ju, Extended quantum anomalous hall states in graphene/hBN moir´ e superlattices, Nature637, 1090 (2025)

  28. [28]

    J. Xie, Z. Huo, X. Lu, Z. Feng, Z. Zhang, W. Wang, Q. Yang, K. Watanabe, T. Taniguchi, K. Liu, Z. Song, X. C. Xie, J. Liu, and X. Lu, Tunable fractional chern insulators in rhombohedral graphene superlattices, Nat. Mater.24, 1042 (2025)

  29. [29]

    Bistritzer and A

    R. Bistritzer and A. H. MacDonald, Moir´ e bands in twisted double-layer graphene, Proc. Natl Acad. Sci. USA108, 12233 (2011)

  30. [30]

    F. Wu, T. Lovorn, E. Tutuc, I. Martin, and A. H. MacDonald, Topological insulators in twisted transition metal dichalcogenide homobilayers, Phys. Rev. Lett.122, 086402 (2019)

  31. [31]

    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)

  32. [32]

    Su´ arez Morell, J

    E. Su´ arez Morell, J. D. Correa, P. Vargas, M. Pacheco, and Z. Barticevic, Flat bands in slightly twisted bilayer graphene: Tight-binding calculations, Phys. Rev. B82, 121407 (2010)

  33. [33]

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

  34. [34]

    J. Yu, S. Qian, and C.-C. Liu, General electronic struc- ture calculation method for twisted systems, Phys. Rev. B111, 075434 (2025)

  35. [35]

    N. Mao, C. Xu, J. Li, T. Bao, P. Liu, Y. Xu, C. Felser, L. Fu, and Y. Zhang, Transfer learning relaxation, elec- 13 tronic structure and continuum model for twisted bilayer mote2, Commun Phys7, 262 (2024)

  36. [36]

    H. Li, Z. Wang, N. Zou, M. Ye, R. Xu, X. Gong, W. Duan, and Y. Xu, Deep-learning density functional theory hamiltonian for efficient ab initio electronic- structure calculation, Nat Comput Sci2, 367 (2022)

  37. [37]

    T. Bao, R. Xu, H. Li, X. Gong, Z. Tang, J. Fu, W. Duan, and Y. Xu, Deep-learning database of density func- tional theory hamiltonians for twisted materials (2024), arXiv:2404.06449 [cond-mat]

  38. [38]

    J. Liu, Z. Fang, H. Weng, and Q. Wu, Dpmoire: A tool for constructing accurate machine learning force fields in moir´ e systems, npj Comput Mater11, 248 (2025)

  39. [39]

    Zhang, H

    Y. Zhang, H. Pi, J. Liu, W. Miao, Z. Qi, N. Regnault, H. Weng, X. Dai, B. A. Bernevig, Q. Wu, and J. Yu, Universal moir´ e-model-building method without fitting: Application to twisted mote$ 2$and wse$ 2$(2024), arXiv:2411.08108 [cond-mat]

  40. [40]

    Li and M

    Y. Li and M. Koshino, Twist-angle dependence of the proximity spin-orbit coupling in graphene on transition- metal dichalcogenides, Phys. Rev. B99, 075438 (2019)

  41. [41]

    J. Shi, J. Zhu, and A. H. MacDonald, Moir´ e commensura- bility and the quantum anomalous Hall effect in twisted bilayer graphene on hexagonal boron nitride, Phys. Rev. B103, 075122 (2021)

  42. [42]

    Bultinck, S

    N. Bultinck, S. Chatterjee, and M. P. Zaletel, Mecha- nism for anomalous hall ferromagnetism in twisted bi- layer graphene, Phys. Rev. Lett.124, 166601 (2020)

  43. [43]

    Mao and T

    D. Mao and T. Senthil, Quasiperiodicity, band topology, and moir´ e graphene, Phys. Rev. B103, 115110 (2021)

  44. [44]

    Zhang, D

    Y.-H. Zhang, D. Mao, Y. Cao, P. Jarillo-Herrero, and T. Senthil, Nearly flat chern bands in moir´ e superlattices, Phys. Rev. B99, 075127 (2019)

  45. [45]

    Nakatsuji, T

    N. Nakatsuji, T. Kawakami, H. Tateishi, K. Kato, and M. Koshino, Moir´ e band engineering in twisted trilayer wse2, Commun Mater6, 274 (2025)

  46. [46]

    Z. Zhu, P. Cazeaux, M. Luskin, and E. Kaxiras, Modeling mechanical relaxation in incommensurate trilayer van der waals heterostructures, Phys. Rev. B101, 224107 (2020)

  47. [47]

    H. Meng, Z. Zhan, and S. Yuan, Commensurate and in- commensurate double moir´ e interference in twisted tri- layer graphene, Phys. Rev. B107, 035109 (2023)

  48. [48]

    Turkel, J

    S. Turkel, J. Swann, Z. Zhu, M. Christos, K. Watanabe, T. Taniguchi, S. Sachdev, M. S. Scheurer, E. Kaxiras, C. R. Dean, and A. N. Pasupathy, Orderly disorder in magic-angle twisted trilayer graphene, Science376, 193 (2022)

  49. [49]

    I. M. Craig, M. Van Winkle, C. Groschner, K. Zhang, N. Dowlatshahi, Z. Zhu, T. Taniguchi, K. Watanabe, S. M. Griffin, and D. K. Bediako, Local atomic stacking and symmetry in twisted graphene trilayers, Nat. Mater. 23, 323 (2024)

  50. [50]

    Y. Li, M. Xue, H. Fan, C.-F. Gao, Y. Shi, Y. Liu, K. Watanabe, T. Tanguchi, Y. Zhao, F. Wu, X. Wang, Y. Shi, W. Guo, Z. Zhang, Z. Fei, and J. Li, Symmetry breaking and anomalous conductivity in a double-moir´ e superlattice, Nano Lett.22, 6215 (2022)

  51. [51]

    X. Lin, C. Li, K. Su, and J. Ni, Energetic stability and spatial inhomogeneity in the local electronic structure of relaxed twisted trilayer graphene, Phys. Rev. B106, 075423 (2022)

  52. [52]

    J. Shin, Y. Park, B. L. Chittari, J.-H. Sun, and J. Jung, Electron-hole asymmetry and band gaps of commensu- rate double moire patterns in twisted bilayer graphene on hexagonal boron nitride, Phys. Rev. B103, 075423 (2021)

  53. [53]

    Nakatsuji, T

    N. Nakatsuji, T. Kawakami, and M. Koshino, Multiscale lattice relaxation in general twisted trilayer graphenes, Phys. Rev. X13, 041007 (2023)

  54. [54]

    D. Park, C. Park, K. Yananose, E. Ko, B. Kim, R. En- gelke, X. Zhang, K. Davydov, M. Green, H.-M. Kim, S. H. Park, J. H. Lee, S.-G. Kim, H. Kim, K. Watan- abe, T. Taniguchi, S. M. Yang, K. Wang, P. Kim, Y.-W. Son, and H. Yoo, Unconventional domain tessellations in moir´ e-of-moir´ e lattices, Nature641, 896 (2025)

  55. [55]

    J. Jung, A. Raoux, Z. Qiao, and A. H. MacDonald, Ab initio theory of moir´ e superlattice bands in layered two- dimensional materials, Phys. Rev. B89, 205414 (2014)

  56. [56]

    J. Jung, E. Laksono, A. M. DaSilva, A. H. MacDonald, M. Mucha-Kruczy´ nski, and S. Adam, Moir´ e band model and band gaps of graphene on hexagonal boron nitride, Phys. Rev. B96, 085442 (2017)

  57. [57]

    B. Hunt, J. D. Sanchez-Yamagishi, A. F. Young, M. Yankowitz, B. J. LeRoy, K. Watanabe, T. Taniguchi, P. Moon, M. Koshino, P. Jarillo-Herrero, and R. C. Ashoori, Massive Dirac Fermions and Hofstadter butter- fly in a van der Waals heterostructure, Science340, 1427 (2013)

  58. [58]

    H. Kim, N. Leconte, B. L. Chittari, K. Watanabe, T. Taniguchi, A. H. MacDonald, J. Jung, and S. Jung, Accurate gap determination in monolayer and bilayer graphene/ h-BN moir´ e superlattices, Nano Lett.18, 7732 (2018)

  59. [59]

    A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys.81, 109 (2009)

  60. [60]

    Weston, D

    L. Weston, D. Wickramaratne, M. Mackoit, A. Alka- uskas, and C. G. Van De Walle, Native point defects and impurities in hexagonal boron nitride, Phys. Rev. B97, 214104 (2018)

  61. [61]

    Yankowitz, J

    M. Yankowitz, J. Jung, E. Laksono, N. Leconte, B. L. Chittari, K. Watanabe, T. Taniguchi, S. Adam, D. Graf, and C. R. Dean, Dynamic band-structure tuning of graphene moir´ e superlattices with pressure, Nature557, 404 (2018)

  62. [62]

    Y. Wang, J. An, C. Ye, X. Wang, D. Mai, H. Zhao, Y. Zhang, C. Peng, K. Watanabe, T. Taniguchi, X. Sun, R. Dai, Z. Wang, W. Qin, Z. Qiao, and Z. Zhang, Pressure-driven moir´ e potential enhancement and ter- tiary gap opening in graphene/h-BN heterostructure, Phys. Rev. Lett.135, 046303 (2025)

  63. [63]

    Liu and X

    J. Liu and X. Dai, Orbital magnetic states in moir´ e graphene systems, Nat. Rev. Phys.3, 367 (2021)

  64. [64]

    J. Liu, J. Liu, and X. Dai, Pseudo landau level representa- tion of twisted bilayer graphene: Band topology and im- plications on the correlated insulating phase, Phys. Rev. B99, 155415 (2019)

  65. [65]

    parity anomal

    F. D. M. Haldane, Model for a quantum Hall effect with- out Landau levels: Condensed-matter realization of the “parity anomal”, Phys. Rev. Lett.61, 2015 (1988)

  66. [66]

    Fukui, Y

    T. Fukui, Y. Hatsugai, and H. Suzuki, Chern numbers in discretized brillouin zone: Efficient method of comput- ing (spin) hall conductances, J. Phys. Soc. Jpn.74, 1674 (2005)