pith. sign in

arxiv: 2605.16218 · v1 · pith:3UY7JDSUnew · submitted 2026-05-15 · ❄️ cond-mat.str-el · cond-mat.mes-hall

Lattice Relaxation Flattens Chern Bands in Rhombohedral Graphene Stacks

Pith reviewed 2026-05-20 15:59 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.mes-hall
keywords lattice relaxationChern insulatorrhombohedral graphenemoiré potentialstrain fieldsHartree-Fockvalley polarizationtopological bands
0
0 comments X

The pith

Lattice relaxation through strain fields flattens and isolates a valley-polarized Chern band with |C|=1 in rhombohedral graphene stacks aligned to hBN.

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

The paper models the moiré potential in rhombohedral graphene on hBN as arising from layer-shear strain fields due to lattice relaxation. Although these fields decay exponentially away from the interface, they still affect electrons in distant layers enough to matter. Without a displacement field, relaxation heightens the electronic distinctions between different stackings, which persist in the moiré-distant regime and strengthen under electron interactions. The key result is that lattice relaxation is essential for flattening and isolating a valley-polarized Hartree-Fock band that carries a Chern number of 1. This work shows how lattice effects and Coulomb interactions together shape topological bands in these systems.

Core claim

In the proposed model, the moiré potential is set by the pattern of layer-shear strain fields from lattice relaxation. In the absence of a displacement field, these effects amplify stacking differences that survive at the single-particle level in the moiré-distant regime and are enhanced by interactions. Lattice relaxation is crucial for flattening and isolating a valley-polarized Hartree-Fock band with Chern number 1.

What carries the argument

The moiré potential defined by the pattern of layer-shear strain fields produced by lattice relaxation in the heterostructure.

If this is right

  • Lattice relaxation amplifies electronic differences between the two stackings with hBN even without displacement field.
  • These differences persist in the moiré-distant regime at single-electron level and are further enhanced by electron interactions.
  • Lattice relaxation is required to flatten and isolate the valley-polarized band with |C|=1.
  • This intertwines long-range Coulomb interactions and lattice relaxation, challenging conventional moiré effect views.
  • Results suggest exploring varied twist angles and displacement fields for finding topological states.

Where Pith is reading between the lines

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

  • If lattice relaxation is key, then standard models that assume rigid lattices or ignore strain may underestimate band flattening in multilayer graphene.
  • Varying the number of layers could test how the exponential decay of strain affects the survival of the Chern band.
  • The mechanism might apply to other van der Waals heterostructures where relaxation-induced strains influence topology.
  • Future calculations could include finite displacement fields to see how the Chern band evolves.

Load-bearing premise

The imprints of exponentially decaying layer-shear strain fields on electrons away from the contact layer are still non-negligible.

What would settle it

A calculation of the Hartree-Fock band structure in these stacks that neglects lattice relaxation entirely, to check if the |C|=1 band still flattens and isolates.

Figures

Figures reproduced from arXiv: 2605.16218 by H\'ector Ochoa, Luca Nashabeh.

Figure 1
Figure 1. Figure 1: FIG. 1. The geometry of the system consists of five layers of [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: As expected, the displacement field un is largest in the first layer, as this is the only layer which directly experiences the hBN adhesion potential driving the re￾laxation. Subsequent layers have an exponentially sup￾pressed tendency to match this displacement, as rhombo￾hedral stacking is a locally stable stacking configuration. Variable Value Source aG 2.46 Å [40] ahBN 2.504 Å [40] d 3.33 Å [19] λ 3.25… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Results of relaxation for [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Calculated single-electron band structure for different stacking configurations and displacement fields, plotted over [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Relaxation helps isolate the filled conduction band in the Hartree band structure at [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. In the moiré proximal regime, relaxation effects create a topological band, with [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
read the original abstract

Motivated by recent observations of integer and fractional Chern insulators in rhombohedral graphene stacks aligned with hexagonal boron nitride (hBN), we propose and study a model in which the moir\'e potential is defined by the pattern of layer-shear strain fields produced by lattice relaxation in these heterostructures. Although these strain fields decrease exponentially with the number of layers, their imprints on electrons residing away from the contact layer are non-negligible. In the absence of a displacement field, lattice relaxation effects amplify the electronic differences among the two different stackings with hBN. These differences, although attenuated at the single-electron level, survive in the so-called moir\'e-distant regime and are further enhanced with the inclusion of electron interactions. We find that lattice relaxation plays a crucial role in flattening and isolating a valley-polarized Hartree-Fock electron band with $|C|=1$ Chern number. Our results challenge the conventional wisdom on moir\'e effects in these heterostructures by illustrating the intertwined effects of long-range Coulomb interactions and lattice relaxation, and opens the door to explore different regimes of twist angles and displacement fields for the search for topological states.

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 / 2 minor

Summary. The manuscript proposes a model for the moiré potential in rhombohedral graphene stacks aligned with hBN, defined by the pattern of layer-shear strain fields from lattice relaxation. Although these fields decay exponentially with layer number, the authors argue their electronic imprints remain non-negligible away from the contact layer. In the absence of a displacement field, relaxation amplifies stacking differences that persist in the moiré-distant regime and are enhanced by interactions; Hartree-Fock calculations then show lattice relaxation is crucial for flattening and isolating a valley-polarized band with |C|=1 Chern number. The work challenges conventional moiré-only pictures by emphasizing the interplay of long-range Coulomb interactions and relaxation.

Significance. If the central result holds, the paper would establish lattice relaxation as an essential ingredient for stabilizing isolated Chern bands in these heterostructures, beyond twist-angle moiré potentials alone. It would open exploration of topological states across twist angles and displacement fields by showing how strain imprints survive attenuation and are amplified by interactions, providing a concrete mechanism for the observed integer and fractional Chern insulators.

major comments (2)
  1. [Abstract and strain model] Abstract and the paragraph on layer-shear strain fields: the claim that exponentially decaying strain fields still produce non-negligible imprints on electrons in distant layers is load-bearing for the assertion that relaxation drives flattening in the moiré-distant regime. The manuscript must supply quantitative support (e.g., explicit decay length versus interlayer hopping, or layer-resolved potential amplitudes) showing the effect is not confined to the contact layer; otherwise the mechanism reduces to the conventional relaxation-free picture.
  2. [Hartree-Fock results] Hartree-Fock band structure section: the flattening and isolation of the |C|=1 valley-polarized band is attributed to the relaxation-induced potential plus interactions. A control calculation omitting the strain fields (or with artificially suppressed distant-layer coupling) is needed to demonstrate that relaxation is indeed crucial rather than merely perturbative.
minor comments (2)
  1. [Methods] Clarify the precise definition of the moiré-distant regime and how it is implemented numerically.
  2. [Figures] Ensure all figures include layer-resolved charge density or potential plots to illustrate the non-negligible distant-layer imprint.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment point by point below and have revised the manuscript to incorporate the requested quantitative support and control calculations.

read point-by-point responses
  1. Referee: [Abstract and strain model] Abstract and the paragraph on layer-shear strain fields: the claim that exponentially decaying strain fields still produce non-negligible imprints on electrons in distant layers is load-bearing for the assertion that relaxation drives flattening in the moiré-distant regime. The manuscript must supply quantitative support (e.g., explicit decay length versus interlayer hopping, or layer-resolved potential amplitudes) showing the effect is not confined to the contact layer; otherwise the mechanism reduces to the conventional relaxation-free picture.

    Authors: We agree that quantitative support is essential to substantiate the non-negligible electronic imprints from decaying strain fields. In the revised manuscript, we have added explicit layer-resolved potential amplitudes and a comparison of the strain decay length (approximately 2 layers) to the interlayer hopping scale. These calculations demonstrate that the effective potential in distant layers remains several meV, sufficient to influence band structure in the moiré-distant regime when amplified by interactions. A new supplementary figure illustrates this layer dependence. revision: yes

  2. Referee: [Hartree-Fock results] Hartree-Fock band structure section: the flattening and isolation of the |C|=1 valley-polarized band is attributed to the relaxation-induced potential plus interactions. A control calculation omitting the strain fields (or with artificially suppressed distant-layer coupling) is needed to demonstrate that relaxation is indeed crucial rather than merely perturbative.

    Authors: We acknowledge the value of a direct control calculation. We have performed additional Hartree-Fock simulations omitting the relaxation-induced strain fields while retaining the conventional moiré potential. The revised manuscript now includes these results, which show that the valley-polarized band remains more dispersive and less isolated without relaxation effects. This confirms the crucial role of lattice relaxation in achieving the observed flattening and |C|=1 Chern band. revision: yes

Circularity Check

0 steps flagged

No significant circularity; model and HF results are independent of input assumptions

full rationale

The paper defines a physical model in which the moiré potential is constructed from computed layer-shear strain fields of lattice relaxation, then performs Hartree-Fock calculations on the resulting single-particle Hamiltonian to obtain the flattened |C|=1 band. No step reduces a prediction to a fitted parameter by construction, nor does any central claim rest on a self-citation chain or imported uniqueness theorem. The exponential decay of strain is treated as an input from elasticity theory and the non-negligible imprint is a numerical outcome of the model, not a definitional tautology. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only; no explicit free parameters, axioms, or invented entities are stated. The model implicitly assumes strain fields remain relevant away from the interface and that Hartree-Fock captures the essential interaction physics.

pith-pipeline@v0.9.0 · 5739 in / 1105 out tokens · 36284 ms · 2026-05-20T15:59:55.222668+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.

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

44 extracted references · 44 canonical work pages

  1. [1]

    Regnault and B

    N. Regnault and B. A. Bernevig, Fractional chern insu- lator, Phys. Rev. X1, 021014 (2011)

  2. [2]

    Neupert, L

    T. Neupert, L. Santos, C. Chamon, and C. Mudry, Frac- tional quantum hall states at zero magnetic field, Phys. Rev. Lett.106, 236804 (2011)

  3. [3]

    D. N. Sheng, Z.-C. Gu, K. Sun, and L. Sheng, Fractional quantum hall effect in the absence of landau levels, Na- ture Communications2, 389 (2011)

  4. [4]

    K. Sun, Z. Gu, H. Katsura, and S. Das Sarma, Nearly flatbands with nontrivial topology, Phys. Rev. Lett.106, 236803 (2011)

  5. [5]

    Tang, J.-W

    E. Tang, J.-W. Mei, and X.-G. Wen, High-temperature fractional quantum hall states, Phys. Rev. Lett.106, 236802 (2011)

  6. [6]

    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)

  7. [7]

    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é- superlattices, Nature637, 1090 (2025)

  8. [8]

    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 in- sulators in rhombohedral graphene superlattices, Nature Materials24, 1042 (2025)

  9. [9]

    S. H. Aronson, T. Han, Z. Lu, Y. Yao, J. P. Butler, K. Watanabe, T. Taniguchi, L. Ju, and R. C. Ashoori, Displacement field-controlled fractional chern insulators and charge density waves in a graphene/hbn moir\’e su- perlattice, Physical Review X15, 031026 (2025)

  10. [10]

    Y. Choi, Y. Choi, M. Valentini, C. L. Patterson, L. F. W. Holleis, O. I. Sheekey, H. Stoyanov, X. Cheng, T. Taniguchi, K. Watanabe, and A. F. Young, Supercon- ductivity and quantized anomalous hall effect in rhom- bohedral graphene, Nature639, 342 (2025)

  11. [11]

    Z. Huo, W. Wang, J. Xie, Y. H. Kwan, J. Herzog- Arbeitman, Z. Zhang, Q. Yang, M. Wu, K. Watanabe, T. Taniguchi, K. Liu, N. Regnault, B. A. Bernevig, and X. Lu, Does moire matter? critical moire dependence with quantum fluctuations in graphene based integer and fractional chern insulators (2025), arXiv:2510.15309 [cond-mat.mes-hall]

  12. [12]

    J. Xie, Z. Zhang, X. Chen, Y. H. Kwan, Z. Huo, J. Herzog-Arbeitman, L. Guo, K. Watanabe, 7 T. Taniguchi, K. Liu, X. C. Xie, B. A. Bernevig, Z.-D. Song, and X. Lu, Unconventional orbital mag- netism in graphene-based fractional chern insulators (2025), arXiv:2506.01485 [cond-mat.mes-hall]

  13. [13]

    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)

  14. [14]

    Y. Zeng, Z. Xia, K. Kang, J. Zhu, P. Knüppel, C. Vaswani, K. Watanabe, T. Taniguchi, K. F. Mak, and J. Shan, Thermodynamic evidence of fractional chern in- sulator in moiré mote2, Nature622, 69 (2023)

  15. [15]

    Z. Dong, A. S. Patri, and T. Senthil, Theory of quan- tum anomalous hall phases in pentalayer rhombohedral graphene moiré structures, Phys. Rev. Lett.133, 206502 (2024)

  16. [16]

    B. Zhou, H. Yang, and Y.-H. Zhang, Fractional quan- tum anomalous hall effect in rhombohedral multilayer graphene in the moiréless limit, Phys. Rev. Lett.133, 206504 (2024)

  17. [17]

    Z. Guo, X. Lu, B. Xie, and J. Liu, Fractional Chern in- sulator states in multilayer graphene moiré superlattices, Phys. Rev. B110, 075109 (2024)

  18. [18]

    Herzog-Arbeitman, Y

    J. Herzog-Arbeitman, Y. Wang, J. Liu, P. M. Tam, Z. Qi, Y. Jia, D. K. Efetov, O. Vafek, N. Regnault, H. Weng, Q. Wu, B. A. Bernevig, and J. Yu, Moiré fractional chern insulators. ii. first-principles calculations and continuum models of rhombohedral graphene superlattices, Phys. Rev. B109, 205122 (2024)

  19. [19]

    Y. H. Kwan, J. Yu, J. Herzog-Arbeitman, D. K. Efetov, N. Regnault, and B. A. Bernevig, Moiré fractional Chern insulators. III. Hartree-Fock phase diagram, magic angle regime for Chern insulator states, role of moiré potential, and Goldstone gaps in rhombohedral graphene superlat- tices, Phys. Rev. B112, 075109 (2025)

  20. [20]

    Huang, X

    K. Huang, X. Li, S. Das Sarma, and F. Zhang, Self- consistent theory of fractional quantum anomalous hall states in rhombohedral graphene, Phys. Rev. B110, 115146 (2024)

  21. [21]

    J. Yu, J. Herzog-Arbeitman, Y. H. Kwan, N. Regnault, and B. A. Bernevig, Moiré fractional chern insulators. iv. fluctuation-driven collapse in multiband exact diago- nalization calculations on rhombohedral graphene, Phys. Rev. B112, 075110 (2025)

  22. [22]

    Huang, S

    K. Huang, S. Das Sarma, and X. Li, Fractional quan- tum anomalous hall effect in rhombohedral multilayer graphene with a strong displacement field, Phys. Rev. B111, 075130 (2025)

  23. [23]

    H. Li, B. A. Bernevig, and N. Regnault, Multiband ex- act diagonalization and an iteration approach to search forfractionalcherninsulatorsinrhombohedralmultilayer graphene, Phys. Rev. B112, 075130 (2025)

  24. [24]

    J. Dong, T. Wang, T. Wang, T. Soejima, M. P. Zaletel, A. Vishwanath, and D. E. Parker, Anomalous hall crys- tals in rhombohedral multilayer graphene. i. interaction- driven chern bands and fractional quantum hall states at zero magnetic field, Phys. Rev. Lett.133, 206503 (2024)

  25. [25]

    Soejima, J

    T. Soejima, J. Dong, T. Wang, T. Wang, M. P. Zale- tel, A. Vishwanath, and D. E. Parker, Anomalous hall crystals in rhombohedral multilayer graphene. ii. gen- eral mechanism and a minimal model, Phys. Rev. B110, 205124 (2024)

  26. [26]

    Z. Dong, A. S. Patri, and T. Senthil, Stability of anoma- lous hall crystals in multilayer rhombohedral graphene, Phys. Rev. B110, 205130 (2024)

  27. [27]

    G. Chen, A. L. Sharpe, E. J. Fox, Y.-H. Zhang, S. Wang, L. Jiang, B. Lyu, H. Li, K. Watanabe, T. Taniguchi, Z. Shi, T. Senthil, D. Goldhaber-Gordon, Y. Zhang, and F. Wang, Tunable correlated chern insulator and ferro- magnetisminamoirésuperlattice,Nature579,56(2020)

  28. [28]

    Waters, A

    D. Waters, A. Okounkova, R. Su, B. Zhou, J. Yao, K. Watanabe, T. Taniguchi, X. Xu, Y.-H. Zhang, J. Folk, and M. Yankowitz, Chern insulators at integer and frac- tional filling in moiré pentalayer graphene, Phys. Rev. X 15, 011045 (2025)

  29. [29]

    Xiang, J

    H. Xiang, J. Ding, J. Hua, N. Liu, W. Zhou, Q. Chen, K. Watanabe, T. Taniguchi, N. Xin, W. Zhu, and S. Xu, Continuously tunable anomalous hall crystals in rhom- bohedral heptalayer graphene (2025), arXiv:2502.18031 [cond-mat.mes-hall]

  30. [30]

    Z. Wang, Q. Liu, X. Han, Z. Li, W. Zhao, Z. Qu, C. Han, K. Watanabe, T. Taniguchi, Z. V. Han, S. Zhou, B. Tong, G. Liu, L. Lu, J. Liu, F. Wu, and J. Lu, Electrical switch- ing of chern insulators in moire rhombohedral heptalayer graphene (2025), arXiv:2503.00837 [cond-mat.mes-hall]

  31. [31]

    B. I. Halperin, Z. Tešanović, and F. Axel, Compatibility of crystalline order and the quantized hall effect, Phys. Rev. Lett.57, 922 (1986)

  32. [32]

    hall crystal

    Z. Tešanović, F. m. c. Axel, and B. I. Halperin, “hall crystal” versus wigner crystal, Phys. Rev. B39, 8525 (1989)

  33. [33]

    Miao and C

    W. Miao and C. Li, Various electronic crystal phases in rhombohedral graphene multilayers, Phys. Rev. B113, 155136 (2026)

  34. [34]

    J. R. Wallbank, Generic miniband structure of graphene on a hexagonal substrate, Phys. Rev. B87, 10.1103/PhysRevB.87.245408 (2013)

  35. [35]

    Mucha-Kruczyński, J

    M. Mucha-Kruczyński, J. R. Wallbank, and V. I. Fal’ko, Heterostructures of bilayer graphene and $h$-BN: Inter- play between misalignment, interlayer asymmetry, and trigonal warping, Phys. Rev. B88, 205418 (2013)

  36. [36]

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

  37. [37]

    Kindermann, B

    M. Kindermann, B. Uchoa, and D. L. Miller, Zero-energy modes and gate-tunable gap in graphene on hexagonal boron nitride, Phys. Rev. B86, 115415 (2012)

  38. [38]

    Moon and M

    P. Moon and M. Koshino, Electronic properties of graphene/hexagonal-boron-nitride moiré superlattice, Phys. Rev. B90, 155406 (2014)

  39. [39]

    San-Jose, A

    P. San-Jose, A. Gutiérrez-Rubio, M. Sturla, and F. Guinea, Electronic structure of spontaneously strained graphene on hexagonal boron nitride, Phys. Rev. B90, 115152 (2014)

  40. [40]

    San-Jose, A

    P. San-Jose, A. Gutiérrez-Rubio, M. Sturla, and F. Guinea, Spontaneous strains and gap in graphene on boron nitride, Phys. Rev. B90, 075428 (2014)

  41. [41]

    Ochoa and R

    H. Ochoa and R. M. Fernandes, Degradation of Phonons in Disordered Moiré Superlattices, Phys. Rev. Lett.128, 065901 (2022)

  42. [42]

    C.-H. Park, N. Bonini, T. Sohier, G. Samsonidze, B. Kozinsky, M. Calandra, F. Mauri, and N. Marzari, Electron–Phonon Interactions and the Intrinsic Electri- cal Resistivity of Graphene, Nano Lett.14, 1113 (2014)

  43. [43]

    Sohier, M

    T. Sohier, M. Calandra, C.-H. Park, N. Bonini, N. Marzari, and F. Mauri, Phonon-limited resistivity 8 of graphene by first-principles calculations: Electron- phonon interactions, strain-induced gauge field, and Boltzmann equation, Phys. Rev. B90, 125414 (2014)

  44. [44]

    C.Fang, M.J.Gilbert,andB.A.Bernevig,Bulktopolog- ical invariants in noninteracting point group symmetric insulators, Phys. Rev. B86, 115112 (2012). 9 SUPPLEMENTARY INFORMATION A. Details of Single-Electron Hamiltonian Lettingξ=±1be the valley index withKξ=ξ(g1−g2)/3(see definitions in Fig. 1 of the main text), andη=±1 be the ABC stacking direction index, ...