Pith. sign in

REVIEW 1 major objections 3 minor 22 references

General continuum model for twisted bilayer graphene and arbitrary smooth deformations

T0 review · 1 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A single continuum Hamiltonian describes twisted bilayer graphene and any smooth deformation of the two layers.

desk verdict Balents gives a clean real-space derivation of the TBG continuum model that subsumes earlier strain treatments and reproduces BM, but the claim of a 'full' Hamiltonian is undercut by the omitted scalar deformation potential. read the letter →

arxiv 1909.01545 v2 pith:BK3TYGEV submitted 2019-09-04 cond-mat.mes-hall cond-mat.str-el

Leon Balents This is my paper · ORCID
classification cond-mat.mes-hallcond-mat.str-el
keywords twistedbilayergraphenecontinuummodelEuleriandisplacementfieldmoirébandsartificialgaugeinterlayertunnelingeffectivetheoryreal-spacederivation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper derives one continuum Hamiltonian for the low-energy electrons of bilayer graphene that stays valid for any smooth deformation of the two layers, not just the rigid twist that previous continuum models were built for. The deformation of each layer is encoded in a displacement field attached to actual positions in the sheet, and the Hamiltonian keeps all terms linear in the small gradients of those fields. For a rigid twist the formula reduces exactly to the standard moiré continuum model. For arbitrary smooth deformations it also covers non-uniform strain, lattice relaxation, and the coupling of electrons to low-energy phonons, all within the same expression.

What carries the argument

The load-bearing object is the Eulerian displacement field $u_l(x)$, which specifies the actual displacement of layer $l$ at the physical point $x$, rather than at a label attached to the undeformed lattice. This choice keeps the Hamiltonian density local in space. The derivation then uses a conformal transformation of the Dirac fields under the change of coordinates, whose determinant factors cancel, an expansion of the Dirac operator to first order in displacement gradients, and a symmetry analysis of the interlayer tunneling matrix under time-reversal-times-twofold rotation, a mirror reflection, and threefold rotation. The symmetries force the tunneling to be the three-term Fourier sum over the reciprocal vectors $Q_j$ with the coefficients $T_j$ shown above, and the AA and AB stacking configurations pin down the two real parameters in $T_j$.

What would settle it

Compute the interlayer tunneling matrix between two layers with a small but non-uniform displacement by exact tight-binding overlap integrals; if the leading correction to $T(u_1-u_2)$ is of order $a\,\partial u$ (the same size as the strain and K-point shift terms the model keeps), then Eq. (39) is incomplete, and the discrepancy would show up as a strain-dependent shift of the lowest moiré bands in numerical band-structure calculations.

Watch

Extended reading notes

Core claim

The central result is the real-space Hamiltonian in Eq. (39), which combines three effects into one expression. Each layer contributes a Dirac kinetic term with Pauli matrices rotated by the displacement gradient, a K-point shift $v(K\cdot\partial_\mu u_l)\tau^\mu$ that moves the Dirac point with the deformation, and an artificial gauge field $A_l$ generated by strain. The interlayer tunneling is a sum over the three shortest reciprocal-lattice vectors $Q_j$, with coefficients $T_j = u I + w(\bar\zeta^j \tau^+ + \zeta^j \tau^-)$ fixed by the lattice symmetries and set by the AA- and AB-stacking interlayer hoppings. When the two layers are rotated by equal and opposite angles, the phase factors $e^{-i Q_j\cdot(u_1-u_2)}$ reproduce the moiré Bloch wavevectors of the standard twisted-bilayer continuum model, so Eq. (39) recovers that model exactly in that limit. Because the displacement fields are otherwise unrestricted, the same formula applies to strains, relaxation, and phonons.

Load-bearing premise

The derivation assumes the displacement gradients are small and, within that expansion, that interlayer tunneling depends only on the relative displacement of the two layers and not on how that displacement varies in space; if gradient corrections to tunneling are not subdominant, Eq. (39) misses terms of the same order as the strain and K-point shift terms.

Editorial extensions

If this is right

  • The same Hamiltonian can be applied to non-uniform strain and twist-angle inhomogeneity, because only gradients of the displacement fields need to be small, not the displacements themselves.
  • Adding dynamics to the displacement fields gives a model of electron-phonon coupling to the original acoustic phonons of the two layers, a starting point for studies of phonon-mediated superconductivity and transport.
  • Twist, strain, and relaxation can be turned on together by writing $u_l = (3-2l)\frac{\theta}{2}\hat z\times x + \hat u_l$, with $\hat u_l$ the strain or phonon part.
  • The K-point shift and the artificial gauge field appear at the same order in displacement gradients, so both must be kept for any deformation that goes beyond a rigid twist.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • An extension the paper does not spell out: the same real-space construction should carry over to other hexagonal van der Waals bilayers by replacing the Dirac valley structure and the symmetry group, leaving the overall form of Eq. (39) unchanged.
  • If gradient corrections to interlayer tunneling turn out not to be subdominant, the model's first failing should appear as a strain-dependent renormalization of the interlayer coupling; a tight-binding benchmark under uniform strain could measure that correction.
  • Treating the displacement fields as dynamical variables suggests a direct parallel with charge-density-wave phase dynamics, in which slow spatial variation of the twist angle would be described by phase-like equations of motion; that connection is not developed in the paper.
  • Keeping the second valley and second-order displacement gradients would provide controlled corrections for intervalley scattering and larger deformations, both outside the stated domain of Eq. (39).
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 3 minor

Summary. This paper presents a real-space derivation of a continuum Hamiltonian for bilayer graphene with arbitrary smooth lattice deformations, parametrized by layer displacement fields with small gradients. The derivation proceeds in three steps: a coordinate transformation to the global (Eulerian) frame, intralayer strain effects encoded through an artificial gauge field, and interlayer tunneling constrained by the remaining discrete symmetries. The final result is Eq. (39). For the special case of a rigid twist, Eq. (39) is shown to reduce to the Bistritzer-MacDonald model in Sec. 3.4. The paper claims that Eq. (39) is the full continuum band Hamiltonian for arbitrary small displacement gradients and that it goes beyond the BM model by describing uniform and non-uniform strains as well as phonon coupling.

Significance. If the completeness claim is corrected, the paper is a useful conceptual and pedagogical contribution. The Eulerian-coordinate formulation is a genuine simplification, and the symmetry-based construction of the interlayer tunneling matrix is transparent and instructive. The explicit recovery of the Bistritzer-MacDonald model in Sec. 3.4 is a nontrivial consistency check that gives confidence in the main mechanism. The assumptions -- small displacement gradients, locality, and weak tunneling -- are stated clearly. However, the central claim that Eq. (39) is the full continuum band Hamiltonian for arbitrary smooth deformations is currently stronger than the derivation supports, because a symmetry-allowed first-order strain term is omitted. The rigid-twist application is unaffected, but the advertised generalizations to strains and phonons are directly affected.

major comments (1)
  1. [Sec. 3.2 and Eq. (39)] The intralayer part of Eq. (39) omits the scalar deformation potential. In a single graphene valley, the most general first-order strain Hamiltonian contains a term g Tr(epsilon_l) psi_l^dagger psi_l with g of order t, in addition to the pseudo-gauge field A_l of Eq. (12) and the K-point shift of Eq. (10). This term is invariant under all the symmetries used in Sec. 3.3, cannot be absorbed into A_l or the K-point shift, and is not subleading compared with the retained terms v(K . d_mu u_l + A_l) tau^mu, whose coefficients are also fixed by the Dirac velocity and lattice scale. For the rigid twist Tr(epsilon)=0 and the BM limit is unaffected, but for the stated applications to non-uniform strains and phonon coupling this term yields a position-dependent scalar potential and therefore changes physical predictions. The Conclusion's phrase 'full continuum band Hamiltonian' overstates what the derivation establishes. The authors should either add this term, with a parameter and a discussion of its magnitude, or explicitly state that Eq. (39) is a minimal model containing only the K-point shift and pseudo-gauge-field couplings, and soften the title-level 'general' claim accordingly.
minor comments (3)
  1. [Sec. 2] The phrase 'expressable in an expansion' is awkward; consider 'expressible as an expansion'.
  2. [Sec. 3.4, Eq. (37)] The sign conventions in the rotated Pauli matrices tau^mu(theta) and the layer-dependent K-point shift terms are correct, but the notation would be easier to follow if the relation to the standard BM conventions for the two layers were spelled out explicitly, since Eq. (38) is used without a derivation.
  3. [Sec. 3.3, after Eq. (23)] The iteration of Eq. (24) starting from Q=0 generates T_{-Q1} and T_{-Q2}, but the text says 'we iterate this relation starting with Q = 0 and generate two further Fourier coefficients before the iteration closes.' It would be clearer to state explicitly that the three vectors -Q_j, j=0,1,2, are the ones kept and that all other Fourier coefficients are approximated as zero.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the continuum model is derived from coordinate transformations, symmetry constraints, and independently fixed AA/AB limits, with the BM limit as a check.

full rationale

The derivation chain is self-contained. Section 3.1 transforms the single-layer Dirac Hamiltonian into global coordinates, so Eq. (9) is a pure coordinate-change result. Section 3.2 imports the standard strain-induced artificial gauge field from the external literature [21, 22], not from the author's own prior work. Section 3.3 constructs the interlayer tunneling by a symmetry analysis: T(u) is expanded in reciprocal-lattice harmonics (Eq. (16)), the coefficients are constrained by C2T, Ry, and C3 (Eqs. (21)-(25)), and the two real parameters u and w are fixed by requiring the Hamiltonian to match the independent AA and AB stacking limits (Eqs. (28)-(31)). No prediction is defined in terms of the target result: the rigid-twist limit is then a derived consequence, Eq. (37), stated to be in perfect agreement with the external BM benchmark. Reproducing BM is a check, not an input. The possible objection that a scalar deformation potential is omitted, so Eq. (39) is not literally the most general first-order strain Hamiltonian, concerns completeness and correctness, not circularity; the paper explicitly states its gradient-neglect assumption in Sec. 3.3, so the derivation does not smuggle in its conclusion. There is no load-bearing self-citation chain, no fitted parameter renamed as a prediction, and no uniqueness theorem imported from the author's prior work. The central Hamiltonian is not equivalent to its inputs by construction.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The model introduces no new entities. The free parameters are the two interlayer tunneling amplitudes, which are standard inputs from microscopic hopping; they are fixed by matching AA and AB stackings. The main axioms are the effective field theory assumptions standard for continuum models.

free parameters (2)
  • u = t'_AA / 3
    Interlayer tunneling strength for the identity channel, fixed by matching the AA-stacked bilayer limit in Sec. 3.3.
  • w = t'_AB / 3
    Interlayer tunneling strength for the off-diagonal channel, fixed by matching the AB-stacked limit.
assumptions (4)
  • standard math Graphene's low-energy electrons are described by massless Dirac fermions with velocity v.
    Used as the starting point in Sec. 3.1, Eq. (4), following standard results from the literature (Ref. [22]).
  • domain assumption Effective field theory locality and analyticity: the Hamiltonian density is local and analytic in the displacement fields and their gradients.
    Assumed in Sec. 2 and used throughout to justify expanding in small gradients and neglecting gradient dependence in tunneling.
  • domain assumption Small displacement gradients and small interlayer hopping: |∂_μ u| << 1 and t' << t.
    Stated in Sec. 4; the derivation keeps only terms linear in these small parameters.
  • domain assumption Valley decoupling: the two valleys do not mix under smooth deformations.
    Assumed in Sec. 3.1, valid when displacement gradients are too small to generate the momentum needed for inter-valley scattering.

how reviews work

0 comments
Cite this review

Pith. "Pith review of General continuum model for twisted bilayer graphene and arbitrary smooth deformations." pith.science (2026). https://pith.science/paper/BK3TYGEV

@misc{pith2026190901545,
  author       = {Pith},
  title        = {Pith review of: General continuum model for twisted bilayer graphene and arbitrary smooth deformations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BK3TYGEV}},
  note         = {Machine review of arXiv:1909.01545}
}
read the original abstract

We present a simple derivation of a continuum Hamiltonian for bilayer graphene with an arbitrary smooth lattice deformation -- technically in a fashion parametrized by displacement fields with small gradients. We show that this subsumes the continuum model of Bistritzer and Macdonald for twisted bilayer grapheneas well as many generalizations and extensions of it. The derivation is carried out entirely in real space.

Figures

Figures reproduced from arXiv: 1909.01545 by the authors.

Figure 1
Figure 1. Left panel: conventions for graphene lattice. Two [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

22 extracted references · 6 canonical work pages

  1. [1]

    Bistritzer and A

    R. Bistritzer and A. H. MacDonald, Moir´ e bands in twisted double-layer graphene, Proceedings of the National Academy of Sciences 108(30), 12233 (2011), doi:10.1073/pnas.1108174108

    doi:10.1073/pnas.1108174108 2011

  2. [2]

    Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. L uo, J. D. Sanchez- Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras et al. , Correlated insulator be- haviour at half-filling in magic-angle graphene superlatti ces, Nature 556(7699), 80 (2018), doi:10.1038/nature26154

    doi:10.1038/nature26154 2018

  3. [3]

    Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. K axiras and P. Jarillo- Herrero, Unconventional superconductivity in magic-angle graphen e superlattices, Nature 556(7699), 43 (2018), doi:10.1038/nature26160

    doi:10.1038/nature26160 2018

  4. [4]

    Yankowitz, S

    M. Yankowitz, S. Chen, H. Polshyn, Y. Zhang, K. Watanabe, T. Taniguchi, D. Graf, A. F. Young and C. R. Dean, Tuning superconductivity in twisted bilayer graphene , Science 363(6431), 1059 (2019), doi:10.1126/science.aav1910

    doi:10.1126/science.aav1910 2019

  5. [5]

    Y. Choi, J. Kemmer, Y. Peng, A. Thomson, H. Arora, R. Polsk i, Y. Zhang, H. Ren, J. Alicea, G. Refael, F. von Oppen, K. Watanabe et al. , Electronic cor- relations in twisted bilayer graphene near the magic angle , Nature Physics (2019), doi:10.1038/s41567-019-0606-5

    doi:10.1038/s41567-019-0606-5 2019

  6. [6]

    A. L. Sharpe, E. J. Fox, A. W. Barnard, J. Finney, K. Watana be, T. Taniguchi, M. A. Kastner and D. Goldhaber-Gordon, Emergent ferromagnetism near three-quarters filling in twisted bilayer graphene , Science 365(6453), 605 (2019), doi:10.1126/science.aaw3780, https://science.sciencemag.org/content/365/6453/605.full.pdf

    doi:10.1126/science.aaw3780 2019

  7. [7]

    H. Yoo, R. Engelke, S. Carr, S. Fang, K. Zhang, P. Cazeaux, S. H. Sung, R. Hovden, A. W. Tsen, T. Taniguchi, K. Watanabe, G.-C. Yi et al. , Atomic and electronic reconstruction at the van der waals interface in twisted bilayer graphene , Nature Materials 18(5), 448 (2019), doi:10.1038/s41563-019-0346-z

    doi:10.1038/s41563-019-0346-z 2019

  8. [8]

    X. Lu, P. Stepanov, W. Yang, M. Xie, M. A. Aamir, I. Das, C. U rgell, K. Watanabe, T. Taniguchi, G. Zhang, A. Bachtold, A. H. MacDonald et al. , Superconductors, or- bital magnets, and correlated states in magic angle bilayer graphene, arXiv preprint arXiv:1903.06513 (2019)

    arXiv 2019

Show all 22 references
  1. [9]

    Kerelsky, L

    A. Kerelsky, L. J. McGilly, D. M. Kennes, L. Xian, M. Yanko witz, S. Chen, K. Watanabe, T. Taniguchi, J. Hone, C. Dean, A. Rubio and A. N. Pasupathy, Maximized electron interactions at the magic angle in twisted bilayer graphene , Nature 572(7767), 95 (2019), doi:10.1038/s4158...

    2019 doi

  2. [10]

    K. Kim, A. DaSilva, S. Huang, B. Fallahazad, S. Larentis , T. Taniguchi, K. Watan- abe, B. J. LeRoy, A. H. MacDonald and E. Tutuc, Tunable moir´ e bands and strong correlations in small-twist-angle bilayer graphen e, Proceedings of the Na- tional Academy of Sciences 114(13), 3...

    2017 doi

  3. [11]

    N. N. Nam and M. Koshino, Lattice relaxation and energy band modula- tion in twisted bilayer graphene , Physical Review B 96(7), 075311 (2017), doi:10.1103/PhysRevB.96.075311

    2017 doi

  4. [12]

    A. Uri, S. Grover, Y. Cao, J. Crosse, K. Bagani, D. Rodan- Legrain, Y. Myasoedov, K. Watanabe, T. Taniguchi, P. Moon et al. , Mapping the twist angle and unconventional landau levels in magic angle graphene , arXiv preprint arXiv:1908.04595 (2019)

    2019 arXiv

  5. [13]

    S. Xu, A. Berdyugin, P. Kumaravadivel, F. Guinea, R. K. K umar, D. Bandurin, S. Mo- rozov, W. Kuang, B. Tsim, S. Liu et al. , Giant oscillations in a triangular network of one-dimensional states in marginally twisted graphene , arXiv preprint arXiv:1905.12984 (2019)

    2019 arXiv

  6. [14]

    Z. Bi, N. F. Q. Yuan and L. Fu, Designing flat bands by strain , Phys. Rev. B 100, 035448 (2019), doi:10.1103/PhysRevB.100.035448

    2019 doi

  7. [15]

    B. Lian, Z. Wang and B. A. Bernevig, Twisted bilayer graphene: A phonon-driven super- conductor, Phys. Rev. Lett. 122, 257002 (2019), doi:10.1103/PhysRevLett.122.257002

    2019 doi

  8. [16]

    F. Wu, A. H. MacDonald and I. Martin, Theory of phonon-mediated super- conductivity in twisted bilayer graphene , Phys. Rev. Lett. 121, 257001 (2018), doi:10.1103/PhysRevLett.121.257001

    2018 doi

  9. [17]

    Y. W. Choi and H. J. Choi, Strong electron-phonon coupling, electron-hole asymmetr y, and nonadiabaticity in magic-angle twisted bilayer graphe ne, Phys. Rev. B 98, 241412 (2018), doi:10.1103/PhysRevB.98.241412

    2018 doi

  10. [18]

    F. Wu, E. Hwang and S. Das Sarma, Phonon-induced giant linear-in-t resistivity in magic angle twisted bilayer graphene: Ordinary strangeness and e xotic superconductivity, Phys. Rev. B 99, 165112 (2019), doi:10.1103/PhysRevB.99.165112

    2019 doi

  11. [19]

    P. A. Lee and T. M. Rice, Electric field depinning of charge density waves , Phys. Rev. B 19, 3970 (1979), doi:10.1103/PhysRevB.19.3970

    1979 doi

  12. [20]

    Fukuyama and P

    H. Fukuyama and P. A. Lee, Dynamics of the charge-density wave. i. impurity pinning in a single chain , Phys. Rev. B 17, 535 (1978), doi:10.1103/PhysRevB.17.535

    1978 doi

  13. [21]

    Suzuura and T

    H. Suzuura and T. Ando, Phonons and electron-phonon scattering in carbon nanotube s, Physical review B 65(23), 235412 (2002), doi:10.1103/PhysRevB.65.235412

    2002 doi

  14. [22]

    A. C. Neto, F. Guinea, N. M. Peres, K. S. Novoselov and A. K . Geim, The electronic properties of graphene , Reviews of modern physics 81(1), 109 (2009), doi:10.1103/RevModPhys.81.109. 12

    2009 doi

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.