pith. sign in

arxiv: 2605.17805 · v1 · pith:LULZ27LBnew · submitted 2026-05-18 · ❄️ cond-mat.mes-hall

Lattice Relaxation in Moir\'e Heterobilayers

Christophe De Beule , Yiyang Lai , Liangtao Peng , Daniel Bennett , Shaffique Adam This is my paper

Pith reviewed 2026-05-20 01:33 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall
keywords lattice relaxationmoiré heterobilayerscontinuum elasticitybuckling instabilityin-plane displacementgraphene on hBNTMD heterobilayersperturbative expressions
0
0 comments X

The pith

An analytical theory from continuum elasticity provides perturbative expressions for lattice relaxation displacements in moiré heterobilayers that agree with numerical solutions and predict a buckling instability near alignment.

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

The paper develops an analytical framework for lattice relaxation in twisted moiré heterobilayers that includes lattice mismatch, twist, external strain, and differing elastic constants. Starting from continuum elasticity, it derives self-consistent equations for the in-plane displacement fields and obtains simple perturbative expressions for the layer-resolved displacements. These expressions match full numerical solutions across experimentally relevant parameters in systems such as graphene on hBN and transition metal dichalcogenide heterobilayers. The theory further demonstrates that compressive in-plane strain from moiré relaxation can drive a buckling instability near perfect alignment.

Core claim

Starting from continuum elasticity, the self-consistent equations for the in-plane displacement fields yield simple perturbative expressions for the layer-resolved in-plane displacement fields induced by lattice relaxation. These analytical results agree very well with full numerical solutions over experimentally relevant parameters. Heterobilayers can exhibit a buckling instability near alignment driven by compressive in-plane strain due to moiré relaxation.

What carries the argument

The self-consistent equations for the in-plane displacement fields derived from continuum elasticity, which enable perturbative solutions for layer-resolved displacements accounting for mismatch, twist, and heterostrain.

If this is right

  • Analytical expressions accurately describe lattice relaxation in graphene on hBN and representative 2H TMD heterobilayers like MoTe2/WSe2 and WSe2/WS2.
  • The framework provides a simple way to incorporate lattice relaxation effects into models of realistic moiré heterostructures.
  • Heterobilayers exhibit a buckling instability near alignment due to compressive in-plane strain from moiré relaxation.

Where Pith is reading between the lines

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

  • This perturbative approach could enable faster computation of relaxed structures in large-scale moiré systems without requiring full numerical relaxation each time.
  • Similar methods might apply to modeling relaxation in other van der Waals heterostructures beyond the examples considered.
  • High-resolution STM or AFM imaging of atomic positions could directly test the predicted layer-resolved displacements and buckling heights.

Load-bearing premise

The continuum elasticity description remains valid at the moiré length scale and the perturbative expansion holds for the range of twists, mismatches, and strains in experiments.

What would settle it

If atomistic simulations or direct measurements of in-plane displacements in a heterobilayer such as graphene on hBN at small twist angles show significant deviations from the predicted perturbative expressions, the analytical theory would be falsified.

Figures

Figures reproduced from arXiv: 2605.17805 by Christophe De Beule, Daniel Bennett, Liangtao Peng, Shaffique Adam, Yiyang Lai.

Figure 1
Figure 1. Figure 1: FIG. 1. First-star analytical theory for (a) the relative in-plane displacement field [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Acoustic displacement fields induced by atomic re [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Buckling transition in moiré heterobilayers. Shown [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Stacking energy landscape (up to an overall constant) [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Fourier components of the relative displacement field [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Fourier components of the relative displacement field [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
read the original abstract

We develop an analytical theory for lattice relaxation in twisted moir\'e heterobilayers, accounting for lattice mismatch, twist, external biaxial heterostrain, and different elastic constants. Starting from continuum elasticity, we derive the self-consistent equations for the in-plane displacement fields and obtain simple perturbative expressions for the layer-resolved in-plane displacement fields induced by lattice relaxation. We apply our theory to graphene on hBN and representative 2H transition metal dichalcogenide heterobilayers, including MoTe$_2$/WSe$_2$ and WSe$_2$/WS$_2$. Our analytical results agree very well with full numerical solutions over experimentally relevant parameters. We further show that heterobilayers can exhibit a buckling instability near alignment, driven by compressive in-plane strain due to moir\'e relaxation. Our results provide a simple theoretical framework for incorporating lattice relaxation in realistic moir\'e heterostructures.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript develops an analytical theory for lattice relaxation in twisted moiré heterobilayers that incorporates lattice mismatch, twist, external biaxial heterostrain, and differing elastic constants. Starting from continuum elasticity, it derives self-consistent equations for the in-plane displacement fields and obtains simple perturbative expressions for the layer-resolved displacements. The theory is applied to graphene/hBN and representative TMD heterobilayers (MoTe₂/WSe₂, WSe₂/WS₂), with analytical results shown to agree well with full numerical solutions over experimentally relevant parameters. The work further predicts a buckling instability near alignment driven by compressive in-plane strain from moiré relaxation.

Significance. If the perturbative solutions remain accurate and the buckling threshold is correctly identified within the linear regime, the paper supplies a simple, parameter-free framework for incorporating lattice relaxation into models of moiré heterostructures. This is valuable for the field, as strain effects strongly influence electronic properties in these systems. The explicit comparison to numerical solutions over a range of twists and mismatches provides concrete validation and strengthens the utility of the analytic expressions.

major comments (1)
  1. [buckling instability analysis] The buckling instability claim (near the end of the manuscript) assumes linear elasticity remains valid when in-plane compression reaches its maximum at zero twist. In this limit the moiré period diverges, the perturbative ordering for the displacements u(r) is least secure, and nonlinear or substrate-adhesion corrections may intervene before the critical compressive strain is reached. The manuscript reports agreement only for “experimentally relevant parameters” without isolating the zero-twist case or demonstrating that the buckling threshold precedes breakdown of the linear approximation.
minor comments (2)
  1. [Abstract] The abstract states that the analytical results “agree very well” with numerics; the main text should quantify this agreement (e.g., maximum relative error or R² values) for each material system and parameter range.
  2. Notation for the layer-resolved displacement fields should be introduced once in the introduction and used consistently thereafter; occasional switches between vector and component notation obscure the perturbative ordering.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their detailed and constructive feedback on our manuscript. We are pleased that the referee recognizes the value of our analytical framework for lattice relaxation in moiré heterobilayers. We address the major comment below.

read point-by-point responses

  1. Referee: [buckling instability analysis] The buckling instability claim (near the end of the manuscript) assumes linear elasticity remains valid when in-plane compression reaches its maximum at zero twist. In this limit the moiré period diverges, the perturbative ordering for the displacements u(r) is least secure, and nonlinear or substrate-adhesion corrections may intervene before the critical compressive strain is reached. The manuscript reports agreement only for “experimentally relevant parameters” without isolating the zero-twist case or demonstrating that the buckling threshold precedes breakdown of the linear approximation.

    Authors: We thank the referee for pointing out this subtlety in the buckling analysis. The maximum compression indeed occurs as the twist angle approaches zero, but the buckling threshold is crossed at a small finite twist where the moiré period is still finite. Our numerical validations include small twist angles, and the perturbative expressions are derived under the assumption of small displacements relative to the moiré period. In the revised manuscript, we will add a dedicated paragraph or subsection that isolates the behavior near zero twist, computes the critical twist angle for the onset of buckling, and verifies that at this point the linear elasticity and perturbative ordering are still valid. We will also discuss the potential impact of nonlinear corrections and substrate adhesion, providing estimates showing that for the material parameters considered, buckling occurs prior to significant deviations from the linear regime. revision: yes

Circularity Check

0 steps flagged

Derivation from continuum elasticity is self-contained with no circular reductions

full rationale

The paper begins from standard continuum elasticity to derive self-consistent equations for in-plane displacement fields u(r), then obtains perturbative expressions for layer-resolved displacements induced by lattice relaxation. These steps constitute a forward derivation from established elastic theory rather than any self-definitional loop, fitted input renamed as prediction, or load-bearing self-citation. The reported agreement with full numerical solutions over experimentally relevant parameters supplies an independent check outside the analytic expressions themselves. No uniqueness theorem, ansatz smuggling, or renaming of known results is invoked in the derivation chain to force the target quantities.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of continuum elasticity at moiré scales and the validity of a perturbative treatment of small displacements; no new free parameters or postulated entities are introduced beyond standard physical inputs such as elastic constants.

axioms (1)
  • domain assumption Continuum elasticity theory applies to the atomic lattices in moiré heterobilayers at the relevant length scales.
    The paper states that it starts from continuum elasticity to derive the self-consistent equations for displacement fields.

pith-pipeline@v0.9.0 · 5694 in / 1423 out tokens · 73296 ms · 2026-05-20T01:33:02.743471+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

72 extracted references · 72 canonical work pages · 1 internal anchor

  1. [1]

    Lattice Relaxation in Moir\'e Heterobilayers

    observed Mott and generalized Wigner crystal states in WSe2/WS2 moiré superlattices. And theoretically, the formation of quantum dots localized at the highly strained nodes of domain wall networks in MoX2/WX2 heterostructures have been proposed [33]. The goal of this paper is to generalize the analytical framework [22, 25] to heterobilayers where lattice ...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  2. [2]

    In addition, we also apply the analytical theory to several 2H TMD heterobilayers in Fig

    For graphene on hBN, we findθ ⊥ ≈0.76 ◦. In addition, we also apply the analytical theory to several 2H TMD heterobilayers in Fig. 2. Now the stacking energy contains non-negligible contributions from the second and third star. For example, up to the second star,u l ≃u l,1 +u l,2 with ul,2(r) = a 2π 3X n=1 u∥ l,2ˆg′ n +u ⊥ l,2ˆz׈g′ n sin(g′ n ·r),(29) wh...

    work page 1922

  3. [3]

    (e) Buckling profile corresponding to (d) using the results from (a) and (c)

    using the one-shot analytical result. (e) Buckling profile corresponding to (d) using the results from (a) and (c). layerl= 1,2, together with Eq. (31), we find a nontrivial buckling solution for ϵ (ϵ2 +θ 2)2 > 32π2 1−3ν κ 3a2V1 ,(33) wheren= 0for the first layer andn= 1for the second layer witha 1 > a 2. Here, we definedν= 1/(1 + 2µ/λ) the 2D Poisson’s r...

    work page

  4. [4]

    This is motivated by the fact that the latter is dominant, while the former is expected to be much smaller

    Perturbation theory It is convenient to write the original layer-resolved fields in terms of the layer center-of-mass (U) and rel- ative displacements (u). This is motivated by the fact that the latter is dominant, while the former is expected to be much smaller. We define these auxiliary fields as u1,2 = U±u 2 .(B16) We can rewrite the elastic energy in ...

    work page

  5. [5]

    +u ⊥ 1 ˆz׈gn sin(gn ·r+α ⊥ 1 ) i ,(B23) U(r) = √ 3a 2π 3X n=1 h U ∥ 1 ˆgn sin(gn ·r+β ∥

    work page

  6. [6]

    For example, we defined 11 u∥ 1,0 =u ∥ 1eiα∥ 1,u ⊥ 1,0 =u ⊥ 1 eiα⊥ 1 , and similarly for the center-of-mass displacement

    +U ⊥ 1 ˆz׈gn sin(gn ·r+β ⊥ 1 ) i ,(B24) where the subscript refers the first star of reciprocal vectors instead of the layer index. For example, we defined 11 u∥ 1,0 =u ∥ 1eiα∥ 1,u ⊥ 1,0 =u ⊥ 1 eiα⊥ 1 , and similarly for the center-of-mass displacement. If we plug these expression in the elastic and adhesion energy, we find helas = a2 L2 n3µ 2 (U ⊥ 1 )2 ...

    work page

  7. [7]

    Elastic energy Including the out-of-plane bending term, the elastic energy of a single layer can be written as [61] Helas = 1 2 Z d2r h λuiiuii + 2µuijuji +κ ∇2h 2i , (C1) with 2D Lamé parametersλandµ, and whereκis the out-of-plane bending rigidity. The strain tensor now also includes the quadratic out-of-plane part: uij = 1 2 ∂uj ∂rj + ∂ui ∂rj + ∂h ∂ri ∂...

    work page

  8. [8]

    Buckling equation In order to find the out-of-plane displacement field un- der periodic boundary conditions for a given in-plane dis- placementfield(specifiedbyu ∥ g andu ⊥ g), weminimizethe elastic energy of a single layer with respect toh−g. To this end, we first compute ∂fij,g′ ∂h−g = gi(gj +g ′ j) +g j(gi +g ′ i) hg+g ′,(C14) ∂Fg′ ∂h−g = 2 (g×g ′)2 g′...

    work page

  9. [9]

    Balents, C

    L. Balents, C. R. Dean, D. K. Efetov, and A. F. Young, Superconductivity and strong correlations in moiré flat bands, Nature Physics16, 725 (2020)

    work page 2020

  10. [10]

    E. Y. Andrei and A. H. MacDonald, Graphene bilayers with a twist, Nature Materials19, 1265 (2020)

    work page 2020

  11. [11]

    K. F. Mak and J. Shan, Semiconductor moiré materials, Nature Nanotechnology17, 686 (2022)

    work page 2022

  12. [12]

    K. S. Novoselov, A. Mishchenko, A. Carvalho, and A. H. Castro Neto, 2D materials and van der Waals het- erostructures, Science353, aac9439 (2016)

    work page 2016

  13. [13]

    E. Y. Andrei, D. K. Efetov, P. Jarillo-Herrero, A. H. MacDonald, K.F.Mak, T.Senthil, E.Tutuc, A.Yazdani, and A. F. Young, The marvels of moiré materials, Nature Reviews Materials6, 201 (2021)

    work page 2021

  14. [14]

    D. M. Kennes, M. Claassen, L. Xian, A. Georges, A. J. Millis, J. Hone, C. R. Dean, D. N. Basov, A. N. Pa- supathy, and A. Rubio, Moiré heterostructures as a condensed-matter quantum simulator, Nature Physics 17, 155 (2021)

    work page 2021

  15. [15]

    C. R. Woods, L. Britnell, A. Eckmann, R. S. Ma, J. C. Lu, H. M. Guo, X. Lin, G. L. Yu, Y. Cao, R. V. Gor- bachev, A. V. Kretinin, J. Park, L. A. Ponomarenko, 13 M. I. Katsnelson, Y. N. Gornostyrev, K. Watanabe, T. Taniguchi, C. Casiraghi, H.-J. Gao, A. K. Geim, and K. S. Novoselov, Commensurate–incommensurate tran- sition in graphene on hexagonal boron ni...

    work page 2014

  16. [16]

    J. Jung, A. M. DaSilva, A. H. MacDonald, and S. Adam, Origin of band gaps in graphene on hexagonal boron ni- tride, Nature Communications6, 6308 (2015)

    work page 2015

  17. [17]

    J. Jung, E. Laksono, A. M. DaSilva, A. H. MacDonald, M. Mucha-Kruczyński, and S. Adam, Moiré band model and band gaps of graphene on hexagonal boron nitride, Physical Review B96, 085442 (2017)

    work page 2017

  18. [18]

    N. N. T. Nam and M. Koshino, Lattice relaxation and en- ergy band modulation in twisted bilayer graphene, Phys- ical Review B96, 075311 (2017)

    work page 2017

  19. [19]

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

    work page 2018

  20. [20]

    Goldhaber-Gordon, Insulating Behavior at the Neu- trality Point in Single-Layer Graphene, Physical Review Letters110, 216601 (2013)

    F.Amet, J.R.Williams, K.Watanabe, T.Taniguchi,and D. Goldhaber-Gordon, Insulating Behavior at the Neu- trality Point in Single-Layer Graphene, Physical Review Letters110, 216601 (2013)

    work page 2013

  21. [21]

    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 But- terfly in a van der Waals Heterostructure, Science340, 1427 (2013)

    work page 2013

  22. [22]

    Bultinck, S

    N. Bultinck, S. Chatterjee, and M. P. Zaletel, Mechanism for Anomalous Hall Ferromagnetism in Twisted Bilayer Graphene, Physical Review Letters124, 166601 (2020)

    work page 2020

  23. [23]

    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)

    work page 2019

  24. [24]

    Watanabe, T

    M.Serlin, C.L.Tschirhart, H.Polshyn, Y.Zhang, J.Zhu, K. Watanabe, T. Taniguchi, L. Balents, and A. F. Young, Intrinsic quantized anomalous Hall effect in a moiré het- erostructure, Science367, 900 (2020)

    work page 2020

  25. [25]

    S. Carr, D. Massatt, S. B. Torrisi, P. Cazeaux, M. Luskin, and E. Kaxiras, Relaxation and domain formation in in- commensurate two-dimensional heterostructures, Physi- cal Review B98, 224102 (2018)

    work page 2018

  26. [26]

    Leconte, S

    N. Leconte, S. Javvaji, J. An, A. Samudrala, and J. Jung, Relaxation effects in twisted bilayer graphene: A multi- scale approach, Physical Review B106, 115410 (2022)

    work page 2022

  27. [27]

    M. H. Naik and M. Jain, Ultraflatbands and Shear Soli- tons in Moiré Patterns of Twisted Bilayer Transition Metal Dichalcogenides, Physical Review Letters121, 266401 (2018)

    work page 2018

  28. [28]

    Cantele, D

    G. Cantele, D. Alfè, F. Conte, V. Cataudella, D. Ninno, and P. Lucignano, Structural relaxation and low-energy properties of twisted bilayer graphene, Physical Review Research2, 043127 (2020)

    work page 2020

  29. [29]

    Lucignano, D

    P. Lucignano, D. Alfè, V. Cataudella, D. Ninno, and G. Cantele, Crucial role of atomic corrugation on the flat bands and energy gaps of twisted bilayer graphene at the magic angle $\theta \approx 1.08^\circ$, Physical Re- view B99, 195419 (2019)

    work page 2019

  30. [30]

    M. M. A. Ezzi, G. N. Pallewela, C. De Beule, E. Mele, and S. Adam, Analytical Model for Atomic Relaxation in Twisted Moiré Materials, Physical Review Letters133, 266201 (2024)

    work page 2024

  31. [31]

    Kang and O

    J. Kang and O. Vafek, Analytical solution for the relaxed atomic configuration of twisted bilayer graphene includ- ing heterostrain, Physical Review B112, 125138 (2025)

    work page 2025

  32. [32]

    J. Yu, B. Wang, and C.-C. Liu, Relaxation and Its Effects on Electronic Structure in Twisted Systems: An Analyt- ical Perspective (2025), arXiv:2509.13114 [cond-mat]

    work page arXiv 2025

  33. [33]

    De Beule, G

    C. De Beule, G. N. Pallewela, M. M. A. Ezzi, L. Peng, E. J. Mele, and S. Adam, Theory for Lattice Re- laxation in Marginal Twist Moirés, arXiv:2503.19162 10.48550/arXiv.2503.19162 (2025), arXiv:2503.19162 [cond-mat] version: 1

    work page doi:10.48550/arxiv.2503.19162 2025

  34. [34]

    Zhang, C

    X.-W. Zhang, C. Wang, X. Liu, Y. Fan, T. Cao, and D. Xiao, Polarization-driven band topology evolution in twisted MoTe2 and WSe2, Nature Communications15, 4223 (2024)

    work page 2024

  35. [35]

    Bennett, Theory of polar domains in moiré het- erostructures, Physical Review B105, 235445 (2022)

    D. Bennett, Theory of polar domains in moiré het- erostructures, Physical Review B105, 235445 (2022)

    work page 2022

  36. [36]

    C. R. Dean, L. Wang, P. Maher, C. Forsythe, F. Ghahari, Y. Gao, J. Katoch, M. Ishigami, P. Moon, M. Koshino, T. Taniguchi, K. Watanabe, K. L. Shepard, J. Hone, and P. Kim, Hofstadter’s butterfly and the fractal quantum Halleffectinmoirésuperlattices,Nature497,598(2013)

    work page 2013

  37. [37]

    L. A. Ponomarenko, R. V. Gorbachev, G. L. Yu, D. C. Elias, R. Jalil, A. A. Patel, A. Mishchenko, A. S. Mayorov, C. R. Woods, J. R. Wallbank, M. Mucha- Kruczynski, B. A. Piot, M. Potemski, I. V. Grigorieva, K. S. Novoselov, F. Guinea, V. I. Fal’ko, and A. K. Geim, Cloning of Dirac fermions in graphene superlattices, Na- ture497, 594 (2013)

    work page 2013

  38. [38]

    D. R. Klein, U. Zondiner, A. Keren, J. Birkbeck, A. In- bar, J. Xiao, Y. Zamir, M. Sidorova, M. M. Al Ezzi, L. Peng, K. Watanabe, T. Taniguchi, S. Adam, and S. Ilani, Imaging the sub-moiré potential using an atomic single electron transistor, Nature650, 875 (2026)

    work page 2026

  39. [39]

    T. Li, S. Jiang, B. Shen, Y. Zhang, L. Li, Z. Tao, T. De- vakul, K. Watanabe, T. Taniguchi, L. Fu, J. Shan, and K. F. Mak, Quantum anomalous Hall effect from inter- twined moiré bands, Nature600, 641 (2021)

    work page 2021

  40. [40]

    E. C. Regan, D. Wang, C. Jin, M. I. Bakti Utama, B. Gao, X. Wei, S. Zhao, W. Zhao, Z. Zhang, K. Yu- migeta, M. Blei, J. D. Carlström, K. Watanabe, T. Taniguchi, S. Tongay, M. Crommie, A. Zettl, and F. Wang, Mott and generalized Wigner crystal states in WSe2/WS2 moiré superlattices, Nature579, 359 (2020)

    work page 2020

  41. [41]

    Soltero, M

    I. Soltero, M. A. Kaliteevski, J. G. McHugh, V. Enaldiev, and V. I. Fal’ko, Competition of Moiré Network Sites to Form Electronic Quantum Dots in Reconstructed MoX2 /WX2 Heterostructures, Nano Letters24, 1996 (2024)

    work page 1996

  42. [42]

    L. D. Landau, L. P. Pitaevskii, A. M. Kosevich, and E. M. Lifshitz,Theory of Elasticity: Volume 7, 3rd ed. (Butterworth-Heinemann, Amsterdam Heidelberg, 1986)

    work page 1986

  43. [43]

    Kim and A

    E.-A. Kim and A. H. C. Neto, Graphene as an electronic membrane, Europhysics Letters84, 57007 (2008)

    work page 2008

  44. [44]

    Escudero, A

    F. Escudero, A. Sinner, Z. Zhan, P. A. Pantaleón, and F. Guinea, Designing Moiré Patterns by Strain (2023), arXiv:2309.08671

    work page arXiv 2023

  45. [45]

    J. Shi, G. Chaudhary, A. H. MacDonald, and I. Martin, Spontaneous Twirls and Structural Frustration in Moiré Materials, Physical Review Letters136, 026101 (2026)

    work page 2026

  46. [46]

    Z.Han, Y.Xia, Z.Xia, W.Zhao, Y.Zhang, K.Watanabe, T.Taniguchi, J.Shan,andK.F.Mak,TopologicalKondo insulator in MoTe2/WSe2 moiré bilayers, Nature Physics 14 22, 396 (2026)

    work page 2026

  47. [47]

    W. Zhao, B. Shen, Z. Tao, Z. Han, K. Kang, K. Watan- abe, T. Taniguchi, K. F. Mak, and J. Shan, Gate-tunable heavy fermions in a moiré Kondo lattice, Nature616, 61 (2023)

    work page 2023

  48. [48]

    B. Gao, D. G. Suárez-Forero, S. Sarkar, T.-S. Huang, D. Session, M. J. Mehrabad, R. Ni, M. Xie, P. Upad- hyay, J. Vannucci, S. Mittal, K. Watanabe, T. Taniguchi, A. Imamoglu, Y. Zhou, and M. Hafezi, Excitonic Mott insulator in a Bose-Fermi-Hubbard system of moiré WS2/WSe2 heterobilayer, Nature Communications15, 2305 (2024)

    work page 2024

  49. [49]

    L. M. Devenica, Z. Hadjri, J. Kumlin, D. Suárez-Forero, R. Li, K. Domi, B. Lyu, W. Li, L. Fausten, V. Vento, N. Ubrig, S. Liu, J. Hone, K. Watanabe, T. Taniguchi, T. Pohl, and A. Srivastava, Collective photon emission and ferroelectric exciton ordering near Mott insulating state in WSe2/WS2 heterobilayers, Nature Materials , 1 (2026)

    work page 2026

  50. [50]

    Pulay, Improved SCF convergence acceleration, Jour- nal of Computational Chemistry3, 556 (1982)

    P. Pulay, Improved SCF convergence acceleration, Jour- nal of Computational Chemistry3, 556 (1982)

    work page 1982

  51. [51]

    Ceferino and F

    A. Ceferino and F. Guinea, Pseudomagnetic fields in fully relaxed twisted bilayer and trilayer graphene, 2D Mate- rials11, 035015 (2024)

    work page 2024

  52. [52]

    W.Bao, F.Miao, Z.Chen, H.Zhang, W.Jang, C.Dames, and C. N. Lau, Controlled ripple texturing of suspended graphene and ultrathin graphite membranes, Nature Nanotechnology4, 562 (2009)

    work page 2009

  53. [53]

    Cai, D.Breid, A.J.Crosby, Z.Suo,andJ.W

    S. Cai, D.Breid, A.J.Crosby, Z.Suo,andJ.W. Hutchin- son, Periodic patterns and energy states of buckled films on compliant substrates, Journal of the Mechanics and Physics of Solids59, 1094 (2011)

    work page 2011

  54. [54]

    Cerda and L

    E. Cerda and L. Mahadevan, Geometry and Physics of Wrinkling, Physical Review Letters90, 074302 (2003)

    work page 2003

  55. [55]

    Watanabe, T

    J.Mao, S.P.Milovanović, M.Anđelković, X.Lai, Y.Cao, K. Watanabe, T. Taniguchi, L. Covaci, F. M. Peeters, A. K. Geim, Y. Jiang, and E. Y. Andrei, Evidence of flat bands and correlated states in buckled graphene super- lattices, Nature584, 215 (2020)

    work page 2020

  56. [56]

    Wang and E

    J. Wang and E. Tosatti, Universal moiré buckling of freestanding 2D bilayers, Proceedings of the National Academy of Sciences121, e2418390121 (2024)

    work page 2024

  57. [57]

    V. T. Phong and E. Mele, Boundary Modes from Peri- odic Magnetic and Pseudomagnetic Fields in Graphene, Physical Review Letters128, 176406 (2022)

    work page 2022

  58. [58]

    De Beule, V

    C. De Beule, V. T. Phong, and E. J. Mele, Roses in the nonperturbative current response of artificial crys- tals, Proceedings of the National Academy of Sciences 120, e2306384120 (2023)

    work page 2023

  59. [59]

    De Beule, R

    C. De Beule, R. Smeyers, W. N. Luna, E. Mele, and L. Covaci, Elastic Screening of Pseudogauge Fields in Graphene, Physical Review Letters134, 046404 (2025)

    work page 2025

  60. [60]

    Enaldiev, V

    V. Enaldiev, V. Zólyomi, C. Yelgel, S. Magorrian, and V. Fal’ko, Stacking Domains and Dislocation Net- works in Marginally Twisted Bilayers of Transition Metal Dichalcogenides, Physical Review Letters124, 206101 (2020)

    work page 2020

  61. [61]

    De Beule and L

    C. De Beule and L. Peng, moirelax (2026)

    work page 2026

  62. [62]

    Gonze, B

    X. Gonze, B. Amadon, P.-M. Anglade, J.-M. Beuken, F.Bottin, P.Boulanger, F.Bruneval, D.Caliste, R.Cara- cas, M. Côté, et al., Abinit: First-principles approach to material and nanosystem properties, Comput. Phys. Commun.180, 2582 (2009)

    work page 2009

  63. [63]

    Gonze and et al., The abinitproject: Impact, environ- ment and recent developments, Comput

    X. Gonze and et al., The abinitproject: Impact, environ- ment and recent developments, Comput. Phys. Commun. 248, 107042 (2020)

    work page 2020

  64. [64]

    M. J. Verstraete, J. Abreu, G. E. Allemand, B. Amadon, G. Antonius, M. Azizi, L. Baguet, C. Barat, L. Bas- togne, R. Béjaud, et al., Abinit 2025: New capabilities for the predictive modeling of solids and nanomaterials, J. Chem. Phys.163, 10.1063/5.0288278 (2025)

    work page doi:10.1063/5.0288278 2025

  65. [65]

    Torrent, F

    M. Torrent, F. Jollet, F. Bottin, G. Zérah, and X. Gonze, Implementation of the projector augmented- wave method in the abinit code: Application to the study of iron under pressure, Comput. Mater. Sci.42, 337 (2008)

    work page 2008

  66. [66]

    Van Setten, M

    M. Van Setten, M. Giantomassi, E. Bousquet, M. J. Ver- straete, D. R. Hamann, X. Gonze, and G.-M. Rignanese, The PseudoDojo: Training and grading a 85 element op- timized norm-conserving pseudopotential table, Comput. Phys. Commun.226, 39 (2018)

    work page 2018

  67. [67]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

    work page 1996

  68. [68]

    H. J. Monkhorst and J. D. Pack, Special points for Brillouin-zone integrations, Phys. Rev. B13, 5188 (1976)

    work page 1976

  69. [69]

    Nelson, David, Tsvi, Piran, and Weinberg, Steven, Statistical Mechanics of Membranes and Surfaces (World Scientific, 2004)

    work page 2004

  70. [70]

    Guinea, B

    F. Guinea, B. Horovitz, and P. Le Doussal, Gauge field induced by ripples in graphene, Physical Review B77, 205421 (2008)

    work page 2008

  71. [71]

    Guinea, B

    F. Guinea, B. Horovitz, and P. Le Doussal, Gauge fields, ripples and wrinkles in graphene layers, Solid State Com- munications Recent Progress in Graphene Studies,149, 1140 (2009)

    work page 2009

  72. [72]

    T. O. Wehling, A. V. Balatsky, A. M. Tsvelik, M. I. Kat- snelson, and A. I. Lichtenstein, Midgap states in corru- gated graphene: Ab initio calculations and effective field theory, Europhysics Letters84, 17003 (2008)

    work page 2008