pith. sign in

arxiv: 2503.08784 · v3 · submitted 2025-03-11 · ❄️ cond-mat.str-el

Elastic Response and Instabilities of Anomalous Hall Crystals

Pith reviewed 2026-05-22 23:56 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords anomalous hall crystalselastic stiffnessmechanical instabilitychern numbermultilayer graphenewigner crystalquantum anomalous hall effecttopological phases
0
0 comments X

The pith

Anomalous Hall crystals have stiffness an order of magnitude smaller than Wigner crystals because their finite Chern number prevents tight charge localization, which triggers mechanical instabilities in graphene models.

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

This paper studies the elastic response of anomalous Hall crystals, phases that simultaneously break translation symmetry and exhibit the quantum anomalous Hall effect. In a continuum model with quadratic dispersion and uniform Berry curvature, time-dependent Hartree-Fock calculations show the AHC stiffness is far lower than that of a Wigner crystal, due to the Chern number blocking exponential localization of charge density. Modifying the dispersion to include a local minimum like that in rhombohedral pentalayer graphene allows deformations that lower kinetic energy to overwhelm the weak stiffness and produce instability. The same instability appears in a microscopic model of rhombohedral pentalayer graphene over experimentally relevant parameters. These findings indicate that the topological character of AHCs must be considered when linking them to observed quantum anomalous Hall phases in multilayer graphene.

Core claim

The stiffness of the anomalous Hall crystal is an order of magnitude smaller than that of the Wigner crystal, attributed to the finite Chern number preventing exponential localization of the charge density. By modifying the dispersion relation to include a local minimum modeled after rhombohedral pentalayer graphene, deformations away from the triangular lattice minimize the kinetic energy, which overwhelms the small stiffness and triggers a mechanical instability. A microscopic model of rhombohedral pentalayer graphene exhibits a similar instability over an experimentally relevant parameter regime.

What carries the argument

Finite Chern number of the anomalous Hall crystal, which blocks exponential localization of charge density and thereby caps the elastic stiffness relative to a topologically trivial Wigner crystal.

If this is right

  • The topologically limited stiffness makes anomalous Hall crystals susceptible to mechanical instabilities.
  • In models of rhombohedral pentalayer graphene, kinetic-energy lowering deformations dominate and drive the system away from the triangular lattice.
  • Any interpretation of integer quantum anomalous Hall observations as anomalous Hall crystals must incorporate the possibility of lattice distortions.
  • The elastic response is set primarily by the topological Chern number rather than by interaction strength alone.

Where Pith is reading between the lines

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

  • Stable anomalous Hall crystal phases may require additional mechanisms, such as substrate effects or longer-range interactions, that lie outside the models studied.
  • Similar stiffness limitations and instabilities could appear in other topological charge-density-wave phases that carry a nonzero Chern number.
  • Lattice-structure probes in multilayer graphene devices could distinguish anomalous Hall crystal explanations from competing scenarios that lack topological stiffness reduction.

Load-bearing premise

The continuum model with quadratic dispersion and uniform Berry curvature, together with the modified dispersion that adds a local minimum modeled on rhombohedral pentalayer graphene, sufficiently represents the physics of candidate systems in multilayer graphene.

What would settle it

A direct measurement showing that the elastic stiffness of an anomalous Hall crystal phase equals or exceeds that of a Wigner crystal, or an experiment finding that the lattice remains undistorted under conditions where the anomalous Hall crystal is proposed to form.

Figures

Figures reproduced from arXiv: 2503.08784 by Adarsh S. Patri, F\'elix Desrochers, Joe Huxford, Mark R. Hirsbrunner, T. Senthil, Yong Baek Kim.

Figure 1
Figure 1. Figure 1: FIG. 1. Representative real space charge density modula [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) A depiction of the dilation deformation of the triangular lattice employed to compute the stiffness. The stiffness [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The stiffness of the WC with [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Continuum R5G dispersion in a strong displace [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original abstract

Anomalous Hall crystals (AHCs) are exotic phases of matter that simultaneously break continuous translation symmetry and exhibit the quantum anomalous Hall effect. AHCs have recently been proposed to explain the observation of an integer quantum anomalous Hall phase in a multilayer graphene system. Despite intense theoretical and experimental interest, little is known about the mechanical properties of AHCs. We study the elastic properties of AHCs first by using a continuum model with quadratic dispersion and uniform Berry curvature. We find using time-dependent Hartree-Fock that the stiffness of the AHC is an order of magnitude smaller than that of the WC, which we attribute to the finite Chern number of the AHC preventing exponential localization of the charge density. By modifying the dispersion relation to include a local minimum modeled after that of rhombohedral pentalayer graphene (R5G), we find that deformations away from the triangular lattice minimize the AHC's kinetic energy, which overwhelms the small stiffness and triggers a mechanical instability. Using a microscopic model of R5G, we observe a similar mechanical instability over an experimentally relevant parameter regime. We conclude that the topologically limited stiffness of AHCs makes them susceptible to mechanical instabilities, an important consideration when interpreting experiments in terms of AHCs.

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

Summary. The paper studies the elastic response of anomalous Hall crystals (AHCs) proposed to explain integer quantum anomalous Hall phases in multilayer graphene. Using time-dependent Hartree-Fock on a continuum model with quadratic dispersion and uniform Berry curvature, it reports that AHC stiffness is an order of magnitude smaller than that of a Wigner crystal, attributed to the finite Chern number preventing exponential charge localization. Modifying the dispersion to include a local minimum modeled after rhombohedral pentalayer graphene (R5G) produces mechanical instabilities because deformations lower kinetic energy enough to overcome the reduced stiffness; a microscopic R5G model reproduces similar instabilities in an experimentally relevant regime. The central conclusion is that topologically limited stiffness renders AHCs susceptible to mechanical instabilities, which must be considered when interpreting experiments.

Significance. If the results hold, the work identifies a concrete mechanical limitation on AHCs that is directly relevant to ongoing experiments in multilayer graphene. The combination of continuum time-dependent Hartree-Fock calculations with a microscopic R5G model provides a useful bridge between idealized topological effects and realistic band structures. The finding that finite Chern number caps stiffness via delocalization is a clear, falsifiable claim that could guide future searches for stable AHC phases.

major comments (3)
  1. [Continuum model section] Continuum-model section (quadratic dispersion, uniform Berry curvature): the attribution of the order-of-magnitude stiffness reduction specifically to the finite Chern number (preventing exponential localization) is load-bearing for the topological claim, yet the text does not report a control calculation with zero Chern number but otherwise identical dispersion and interaction; without this isolation, the topological origin remains an interpretation rather than a demonstrated necessity.
  2. [Modified dispersion section] Modified-dispersion section: the instability conclusion rests on the statement that kinetic-energy lowering from lattice deformations overwhelms the small elastic stiffness. No explicit numerical comparison of the elastic energy scale versus the kinetic-energy gain (e.g., via the second-derivative matrix or total-energy differences at small strain) is provided, leaving the mechanism qualitative even though the microscopic R5G calculation is invoked to confirm it.
  3. [Microscopic R5G model] Microscopic R5G model: the instability is reported only inside an 'experimentally relevant' parameter regime, but the manuscript does not tabulate the precise values of interaction strength, displacement-field range, or filling factors used, nor does it show how the instability boundary shifts when the actual (non-uniform) Berry curvature of R5G is retained instead of the uniform approximation.
minor comments (2)
  1. [Abstract] The abstract states the stiffness is 'an order of magnitude smaller' but supplies no numerical values, error estimates, or figure references; adding a short table or quoted numbers from the time-dependent Hartree-Fock runs would make the quantitative claim immediately verifiable.
  2. [Modified dispersion section] Notation for the modified dispersion (local minimum 'modeled after' R5G) should be defined with an explicit functional form or equation number so that readers can reproduce the kinetic-energy minimization argument.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment below.

read point-by-point responses
  1. Referee: [Continuum model section] Continuum-model section (quadratic dispersion, uniform Berry curvature): the attribution of the order-of-magnitude stiffness reduction specifically to the finite Chern number (preventing exponential localization) is load-bearing for the topological claim, yet the text does not report a control calculation with zero Chern number but otherwise identical dispersion and interaction; without this isolation, the topological origin remains an interpretation rather than a demonstrated necessity.

    Authors: We thank the referee for this observation. The Wigner crystal serves as our primary control, exhibiting exponential localization and an order-of-magnitude higher stiffness in the absence of finite Chern number. Constructing an otherwise identical model with zero Chern number while retaining uniform Berry curvature is not straightforward, since a uniform Berry curvature in a single gapped band implies nonzero Chern number. We will revise the text to more explicitly frame the WC comparison as supporting the topological origin while acknowledging the limitation of the available controls. revision: partial

  2. Referee: [Modified dispersion section] Modified-dispersion section: the instability conclusion rests on the statement that kinetic-energy lowering from lattice deformations overwhelms the small elastic stiffness. No explicit numerical comparison of the elastic energy scale versus the kinetic-energy gain (e.g., via the second-derivative matrix or total-energy differences at small strain) is provided, leaving the mechanism qualitative even though the microscopic R5G calculation is invoked to confirm it.

    Authors: We agree that an explicit numerical comparison would make the mechanism more quantitative. In the revised manuscript we will add direct comparisons, including the elastic energy cost obtained from the stiffness tensor at small strains and the corresponding kinetic-energy lowering from the modified dispersion, to demonstrate that the kinetic term dominates. revision: yes

  3. Referee: [Microscopic R5G model] Microscopic R5G model: the instability is reported only inside an 'experimentally relevant' parameter regime, but the manuscript does not tabulate the precise values of interaction strength, displacement-field range, or filling factors used, nor does it show how the instability boundary shifts when the actual (non-uniform) Berry curvature of R5G is retained instead of the uniform approximation.

    Authors: We will add an explicit table or list of the precise parameter values (interaction strength, displacement-field range, and filling factors) employed in the microscopic R5G calculations. The microscopic model already incorporates the actual non-uniform Berry curvature of R5G; the fact that the instability persists in this setting indicates robustness beyond the uniform approximation used in the continuum model. We will add clarifying text to this effect. revision: partial

Circularity Check

0 steps flagged

No circularity: derivation uses explicit models and computations without reduction to inputs by construction

full rationale

The paper computes elastic stiffness via time-dependent Hartree-Fock on a continuum model with quadratic dispersion and uniform Berry curvature, attributes the reduction relative to Wigner crystal to the finite Chern number via the resulting charge density (not by definitional equivalence), then explicitly modifies the dispersion to add a local minimum modeled after R5G and separately applies a microscopic R5G model to show instability. No steps fit parameters to data subsets then rename the output as prediction, invoke self-citations as load-bearing uniqueness theorems, or smuggle ansatzes via prior work. All results follow from the stated model Hamiltonians and numerical methods without the central claims reducing to the inputs by the paper's own equations.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Based on abstract only; paper relies on standard condensed-matter approximations without introducing new free parameters or entities visible here.

axioms (2)
  • domain assumption Time-dependent Hartree-Fock approximation is adequate for computing elastic properties of AHCs
    Invoked for the stiffness calculation in the continuum model.
  • domain assumption The modified dispersion with local minimum accurately represents rhombohedral pentalayer graphene
    Used to trigger the instability result.

pith-pipeline@v0.9.0 · 5777 in / 1257 out tokens · 46000 ms · 2026-05-22T23:56:37.768764+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.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Emergence of Topological Electron Crystals in Bilayer Graphene--Mott Insulator Heterostructures

    cond-mat.mes-hall 2025-12 unverdicted novelty 6.0

    Interlayer attraction in bilayer graphene-Mott insulator heterostructures stabilizes topological electron crystals with triangular, honeycomb, and kagome geometries via self-consistent Hartree-Fock calculations.

  2. Various electronic crystal phases in rhombohedral graphene multilayers

    cond-mat.mes-hall 2025-12 unverdicted novelty 5.0

    Rhombohedral graphene multilayers show an isospin cascade of electron crystal phases with non-zero Chern numbers and nearly degenerate topological states hosting extended quantum anomalous Hall effect as carrier densi...

Reference graph

Works this paper leans on

71 extracted references · 71 canonical work pages · cited by 2 Pith papers · 1 internal anchor

  1. [1]

    J. Xie, Z. Huo, X. Lu, Z. Feng, Z. Zhang, W. Wang, Q. Yang, K. Watanabe, T. Taniguchi, K. Liu, et al. , Tunable fractional Chern insulators in rhombohedral graphene superlattices, arXiv preprint arXiv:2405.16944 (2024)

  2. [2]

    Waters, A

    D. Waters, A. Okounkova, R. Su, B. Zhou, J. Yao, K. Watanabe, T. Taniguchi, X. Xu, Y.-H. Zhang, J. Folk, et al. , Interplay of electronic crystals with integer and fractional Chern insulators in moiré pentalayer graphene, arXiv preprint arXiv:2408.10133 (2024)

  3. [3]

    S. H. Aronson, T. Han, Z. Lu, Y. Yao, 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é superlattice, arXiv preprint arXiv:2408.11220 (2024)

  4. [4]

    Y. Choi, Y. Choi, M. Valentini, C. L. Patterson, L. F. Holleis, O. I. Sheekey, H. Stoyanov, X. Cheng, T. Taniguchi, K. Watanabe, et al., Electric field control of superconductivity and quantized anomalous Hall effects in rhombohedral tetralayer graphene, arXiv preprint arXiv:2408.12584 (2024)

  5. [5]

    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. B 109, 205122 (2024)

  6. [6]

    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. B 110, 115146 (2024)

  7. [7]

    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. B 111, 075130 (2025)

  8. [8]

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

  9. [9]

    Bonsall and A

    L. Bonsall and A. A. Maradudin, Some static and dy- namical properties of a two-dimensional Wigner crystal, Phys. Rev. B 15, 1959 (1977)

  10. [10]

    Côté and A

    R. Côté and A. H. MacDonald, Collective modes of the two-dimensional Wigner crystal in a strong magnetic field, Phys. Rev. B 44, 8759 (1991)

  11. [11]

    Chitra, T

    R. Chitra, T. Giamarchi, and P. Le Doussal, Pinned Wigner crystals, Phys. Rev. B 65, 035312 (2001)

  12. [12]

    M. M. Fogler and D. A. Huse, Dynamical response of a pinned two-dimensional Wigner crystal, Phys. Rev. B 62, 7553 (2000)

  13. [13]

    Côté, M.-A

    R. Côté, M.-A. Lemonde, C. B. Doiron, and A. M. Et- touhami, Dynamical matrix of two-dimensional electron crystals, Phys. Rev. B 77, 115303 (2008)

  14. [14]

    Marzari, A

    N. Marzari, A. A. Mostofi, J. R. Yates, I. Souza, and D. Vanderbilt, Maximally localized Wannier functions: Theory and applications, Rev. Mod. Phys. 84, 1419 (2012)

  15. [15]

    Brouder, G

    C. Brouder, G. Panati, M. Calandra, C. Mourougane, and N. Marzari, Exponential localization of Wannier functions in insulators, Phys. Rev. Lett. 98, 046402 (2007)

  16. [16]

    Q. Li, J. Dong, P. J. Ledwith, and E. Khalaf, Con- straints on real space representations of Chern bands, arXiv preprint arXiv:2407.02561 (2024)

  17. [17]

    S. A. Parameswaran, R. Roy, and S. L. Sondhi, Fractional quantum Hall physics in topological flat bands, Comptes Rendus Physique 14, 816 (2013)

  18. [18]

    J. Yu, B. A. Bernevig, R. Queiroz, E. Rossi, P. Törmä, and B.-J. Yang, Quantum geometry in quantum materi- als, arXiv preprint arXiv:2501.00098 (2024)

  19. [19]

    See Supplemental Material

  20. [21]

    R. J. Hunt, M. Szyniszewski, G. I. Prayogo, R. Maezono, and N. D. Drummond, Quantum monte carlo calculations of energy gaps from first principles, Phys. Rev. B 98, 075122 (2018)

  21. [22]

    2(e)) at large interaction strengths, but this is an artifact arising from the finite momentum cutoff

    We note that the square lattice appears lower in energy than the triangular lattice (Fig. 2(e)) at large interaction strengths, but this is an artifact arising from the finite momentum cutoff. Indeed, the perturbative calculation predicts that the triangular lattice always remains stable and lower in energy than the square lattice for large Vc. Computing th...

  22. [23]

    Abouelkomsan, K

    A. Abouelkomsan, K. Yang, and E. J. Bergholtz, Quan- tum metric induced phases in moiré materials, Phys. Rev. Res. 5, L012015 (2023)

  23. [24]

    Claassen, C

    M. Claassen, C. H. Lee, R. Thomale, X.-L. Qi, and T. P. Devereaux, Position-momentum duality and fractional quantum Hall effect in Chern insulators, Phys. Rev. Lett. 114, 236802 (2015)

  24. [25]

    J. Wang, J. Cano, A. J. Millis, Z. Liu, and B. Yang, Exact Landau level description of geometry and interaction in a flatband, Phys. Rev. Lett. 127, 246403 (2021)

  25. [26]

    P. J. Ledwith, A. Vishwanath, and D. E. Parker, Vor- texability: A unifying criterion for ideal fractional Chern insulators, Phys. Rev. B 108, 205144 (2023)

  26. [27]

    Estienne, N

    B. Estienne, N. Regnault, and V. Crépel, Ideal Chern bands as Landau levels in curved space, Phys. Rev. Res. 5, L032048 (2023)

  27. [28]

    Qi, Generic wave-function description of fractional quantum anomalous Hall states and fractional topologi- cal insulators, Phys

    X.-L. Qi, Generic wave-function description of fractional quantum anomalous Hall states and fractional topologi- cal insulators, Phys. Rev. Lett. 107, 126803 (2011). 8

  28. [29]

    S. A. Parameswaran, R. Roy, and S. L. Sondhi, Fractional Chern insulators and the W∞ algebra, Phys. Rev. B 85, 241308 (2012)

  29. [30]

    Roy, Band geometry of fractional topological insula- tors, Phys

    R. Roy, Band geometry of fractional topological insula- tors, Phys. Rev. B 90, 165139 (2014)

  30. [31]

    Tanatar and D

    B. Tanatar and D. M. Ceperley, Ground state of the two- dimensional electron gas, Phys. Rev. B 39, 5005 (1989)

  31. [32]

    Rapisarda and G

    F. Rapisarda and G. Senatore, Diffusion monte carlo study of electrons in two-dimensional layers, Australian journal of physics 49, 161 (1996)

  32. [33]

    Xie and S

    M. Xie and S. Das Sarma, Integer and fractional quan- tum anomalous Hall effects in pentalayer graphene, Phys. Rev. B 109, L241115 (2024)

  33. [34]

    Impurity-induced thermal crossover in fractional Chern insulators

    K. Huang, S. D. Sarma, and X. Li, Impurity-induced thermal crossover in fractional Chern insulators, arXiv preprint arXiv:2409.04349 (2024)

  34. [35]

    Das Sarma and M

    S. Das Sarma and M. Xie, Thermal crossover from a Chern insulator to a fractional Chern insulator in penta- layer graphene, Phys. Rev. B 110, 155148 (2024)

  35. [36]

    Shavit, Entropy-enhanced fractional quantum anoma- lous Hall effect, Phys

    G. Shavit, Entropy-enhanced fractional quantum anoma- lous Hall effect, Phys. Rev. B 110, L201406 (2024)

  36. [37]

    Elastic Response and Instabilities of Anomalous Hall Crystals

    Z. Wei, A.-K. Wu, M. Gonçalves, and S.-Z. Lin, Edge- driven transition between extended quantum anomalous Hall crystal and fractional Chern insulator in rhombo- hedral graphene multilayers, Physical Review B 111, 035116 (2025). Supplemental Material for “Elastic Response and Instabilities of Anomalous Hall Crystals” Félix Desrochers ,1, ∗ Mark R. Hirsbrun...

  37. [38]

    Moiré rhombohedral graphene Hamiltonian 20

  38. [39]

    Hartree-Fock calculations 22 C

    Band basis 21 B. Hartree-Fock calculations 22 C. Existence of the AHC 22 D. Lattice deformations 23 E. Optimal orientation of the AHC 25 F. Elastic properties of the AHC 25 References 27 I. MECHANICAL RESPONSE In this section, we comment on how the shear and dilation stiffnesses defined in the main text are related to the usual elastic coefficients that appea...

  39. [40]

    At n1 = 12 , the system size used to compute the stiffnesses in the main text, we see we have captured roughly 80% of the correlation energy stiffness, so our results are qualitatively unchanged by the finite system size and our conclusions hold. IV. V ARIA TIONAL ANSA TZ: PER TURBA TIVE ENERGY CALCULA TION Here we employ a variational ansatz for the parent ...

  40. [41]

    42 (for a relative contribution of order 0

    1) or ξ > 0. 42 (for a relative contribution of order 0. 01). Given that there will be more terms contributing at higher 13 order, a larger ξ may be required for good convergence. Note that this estimate for the convergence is in terms of the variational parameter ξ rather than a physical parameter such as Vc. We discuss how the regime of applicability de...

  41. [42]

    (S65b) The associated real space basis vectors respect Ai · Gj = 2πδ ij

    Moiré lattice The initial graphene reciprocal lattice basis vectors are G1 = 4π √ 3aG (0, 1) (S65a) G2 = 4π√ 3aG ( − √ 3, 1 ) . (S65b) The associated real space basis vectors respect Ai · Gj = 2πδ ij. The real space basis vectors of the hBN substrate are obtained as A′ j =MR [θ ]Aj; M = 1 1 +εI, (S66) whereR[θ ] is a counter-clockwise rotation matrix and ...

  42. [43]

    77◦ reported in experiments [18], leads to a moiré lattice constant of λ M = |AmBZ i | ≈ aG/ √ ε 2 +θ 2 ≈ 11. 4 nm

  43. [44]

    Moiré rhombohedral graphene Hamiltonian The moiré rhombohedral graphene Hamiltonian is H =HR5G +HM +HC, (S68) where HR5G is the continuum kinetic term, HM the moiré potential from the hBN substrate that is acting on the bottom graphene layer, and HC the Coulomb potential. The continuum kinetic term is obtained by expanding the rhombohedral pentalayer grap...

  44. [45]

    charge neutrality scheme

    Band basis Due to the large number of bands (coming from the sublattices, layers, and reciprocal lattice vectors within the momentum cutoff considered), the Hartree-Fock calculation will be done by only considering the nbands lowest con- duction bands [17, 23–26]. To do so, it is first convenient to work in a band basis obtained by diagonalizing the 22 quad...

  45. [46]

    L. D. Landau, L. Pitaevskii, A. M. Kosevich, and E. M. Lifshitz, Theory of elasticity: volume 7 , Vol. 7 (Elsevier, 2012)

  46. [47]

    Côté, M.-A

    R. Côté, M.-A. Lemonde, C. B. Doiron, and A. M. Ettouhami, Dynamical matrix of two-dimensional electron crystals, Phys. Rev. B 77, 115303 (2008)

  47. [48]

    Tan and T

    T. Tan and T. Devakul, Parent Berry curvature and the ideal anomalous Hall crystal, Phys. Rev. X 14, 041040 (2024)

  48. [49]

    T. Tan, J. May-Mann, and T. Devakul, Wavefunction approach to the fractional anomalous Hall crystal, arXiv preprint arXiv:2409.06775 (2024)

  49. [50]

    Pulay, Convergence acceleration of iterative sequences

    P. Pulay, Convergence acceleration of iterative sequences. the case of SCF iteration, Chemical Physics Letters 73, 393 (1980)

  50. [51]

    Pulay, Improved SCF convergence acceleration, Journal of Computational Chemistry 3, 556 (1982)

    P. Pulay, Improved SCF convergence acceleration, Journal of Computational Chemistry 3, 556 (1982)

  51. [52]

    Kresse and J

    G. Kresse and J. Furthmüller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Physical review B 54, 11169 (1996)

  52. [53]

    Rohwedder and R

    T. Rohwedder and R. Schneider, An analysis for the DIIS acceleration method used in quantum chemistry calculations, Journal of mathematical chemistry 49, 1889 (2011)

  53. [54]

    Z. Dong, A. S. Patri, and T. Senthil, Stability of anomalous Hall crystals in multilayer rhombohedral graphene, Phys. Rev. B 110, 205130 (2024)

  54. [55]

    D. J. Rowe, Equations-of-Motion Method and the Extended Shell Model, Reviews of Modern Physics 40, 153 (1968)

  55. [56]

    Khalaf, N

    E. Khalaf, N. Bultinck, A. Vishwanath, and M. P. Zaletel, Soft modes in magic angle twisted bilayer graphene (2020), arXiv:2009.14827 [cond-mat]

  56. [57]

    Y. H. Kwan, Y. Hu, S. H. Simon, and S. A. Parameswaran, Exciton Band Topology in Spontaneous Quantum Anomalous Hall Insulators: Applications to Twisted Bilayer Graphene, Physical Review Letters 126, 137601 (2021)

  57. [59]

    T. M. R. Wolf and C. Huang, Quasi-boson approximation yields accurate correlation energy in the 2D electron gas, Physical Review Research 6, 033296 (2024), arXiv:2311.11449 [cond-mat]

  58. [60]

    Jain and C

    A. Jain and C. Huang, Elementary excitations, melting temperature and correlation energy in wigner crystal, arXiv preprint arXiv:2504.04538 (2025)

  59. [61]

    Bateman and B

    H. Bateman and B. M. Project, Higher Transcendental Functions [Volumes I-III] (McGraw-Hill Book Company, 1953) funding by Office of Naval Research (ONR)

  60. [62]

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

  61. [63]

    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, Nature 626 (2024)

  62. [64]

    Zhang, B

    F. Zhang, B. Sahu, H. Min, and A. H. MacDonald, Band structure of abc-stacked graphene trilayers, Phys. Rev. B 82, 035409 (2010)

  63. [65]

    Y. Park, Y. Kim, B. L. Chittari, and J. Jung, Topological flat bands in rhombohedral tetralayer and multilayer graphene on hexagonal boron nitride moiré superlattices, Phys. Rev. B 108, 155406 (2023)

  64. [66]

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

  65. [67]

    J. Jung, A. M. DaSilva, A. H. MacDonald, and S. Adam, Origin of band gaps in graphene on hexagonal boron nitride, Nature communications 6, 6308 (2015)

  66. [68]

    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, the role of the moiré potential and Goldstone gaps in rhombohedral graphene superlattices, arXiv preprint arXiv:2312.11617 (2023)

  67. [69]

    J. Dong, T. Wang, T. Wang, T. Soejima, M. P. Zaletel, A. Vishwanath, and D. E. Parker, Anomalous Hall crystals in rhombohedral multilayer graphene. I. Interaction-driven Chern bands and fractional quantum Hall states at zero magnetic field, Phys. Rev. Lett. 133, 206503 (2024)

  68. [70]

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

  69. [71]

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

  70. [72]

    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. B 110, 115146 (2024)

  71. [73]

    Huang, S

    K. Huang, S. Das Sarma, and X. Li, Fractional quantum anomalous Hall effect in rhombohedral multilayer graphene with a strong displacement field, Phys. Rev. B 111, 075130 (2025)