Twisted Bilayer Graphene in Commensurate Angles
Pith reviewed 2026-05-23 20:43 UTC · model grok-4.3
The pith
Dirac cones exist at the Brillouin zone vertices in the exact continuum Schrödinger operator for twisted bilayer graphene at commensurate angles.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We prove the existence of Dirac cones at the vertices of the Brillouin zone for such angles in the 2D continuum model of electronic transport in twisted bilayer graphene at commensurate angles. The model uses two honeycomb potentials with the symmetries of graphene, either sharing a common origin or shifted by a half-lattice spacing, and twisted relative to each other. Quantitative bounds show that for small potentials the slope of the Dirac cones flattens at commensurate angles near incommensurate angles.
What carries the argument
The exact continuum Schrödinger operator without the Bistritzer-MacDonald approximations, for commensurate twist angles with AA or AB stacking of graphene-symmetric honeycomb potentials.
If this is right
- The existence of Dirac cones holds for a wide class of potentials in both stacking types.
- Quantitative bounds are established for the flattening of the Dirac cone slope for small potentials.
- This provides the first rigorous proof of Dirac cones in the continuum setting for TBG at commensurate angles.
Where Pith is reading between the lines
- The proof technique might extend to incommensurate angles or other 2D heterostructures.
- Experimental verification could involve measuring the density of states or ARPES near commensurate angles to observe the flattening effect.
- Similar continuum models for other twisted materials could benefit from this approach to establish topological features rigorously.
Load-bearing premise
The two honeycomb potentials possess the symmetries of graphene and the twist angles are exactly the commensurate angles defined by the model for AA or AB stacking.
What would settle it
A calculation or numerical simulation of the spectrum of the Schrödinger operator at a commensurate angle that does not exhibit linear dispersion at the Brillouin zone vertices would falsify the claim.
Figures
read the original abstract
We study a 2D continuum model of electronic transport in twisted bilayer graphene (TBG) at commensurate angles. We use two honeycomb potentials with the symmetries of graphene, either sharing a common origin (AA stacking) or shifted by a half-lattice spacing (AB stacking), and twisted relative to each other. While the electronic properties of TBG are most commonly studied via the approximate Bistritzer-MacDonald (BM) model, our approach studies the exact continuum Schr\"{o}dinger operator without these approximations. Our results hold for a wide class of potentials in both stacking types. We describe the exact angles for which the two twisted lattices are commensurate and prove the existence of Dirac cones at the vertices of the Brillouin zone for such angles. Additionally, we establish quantitative bounds showing that, for small potentials, the slope of the Dirac cones flattens at commensurate angles near incommensurate angles. This work is the first to rigorously establish the existence of Dirac cones for twisted bilayer graphene in the continuum setting, without the BM approximations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes the exact 2D continuum Schrödinger operator for twisted bilayer graphene at commensurate angles, employing two honeycomb potentials with graphene symmetries under either AA or AB stacking. It proves existence of Dirac cones at the Brillouin-zone vertices for these angles and derives quantitative bounds showing flattening of the Dirac-cone slope for small potentials near incommensurate angles. The work is presented as the first rigorous continuum proof without the Bistritzer-MacDonald approximation.
Significance. If the central existence proof and bounds hold, the result supplies the first mathematically rigorous demonstration of Dirac cones in the exact continuum model of TBG at commensurate angles. The quantitative flattening bounds constitute an additional concrete contribution that could be compared with numerical or approximate-model studies.
major comments (2)
- [Main existence result (abstract and presumed §3–4)] The abstract states that the existence proof and bounds hold for a wide class of potentials under the stated symmetries, yet the provided text gives no explicit statement of the spectral estimates or perturbation argument used to control the Schrödinger operator near the Brillouin-zone vertices; without these details the load-bearing step of the central claim cannot be verified.
- [Quantitative bounds paragraph (abstract)] The quantitative bounds on slope flattening are asserted for small potentials; the manuscript must specify whether the constants are uniform over the admissible potential class or depend on additional parameters, as this directly affects the strength of the flattening statement.
minor comments (2)
- [Model description] The definition of the exact commensurate angles (AA/AB) should be stated with an explicit formula or lattice-vector condition in the model section.
- [Introduction] A brief comparison table or paragraph contrasting the present continuum operator with the BM model would clarify the precise sense in which the approximation is avoided.
Simulated Author's Rebuttal
We thank the referee for their report and the opportunity to clarify our manuscript. We respond to each major comment below and will revise the manuscript to address the points raised.
read point-by-point responses
-
Referee: [Main existence result (abstract and presumed §3–4)] The abstract states that the existence proof and bounds hold for a wide class of potentials under the stated symmetries, yet the provided text gives no explicit statement of the spectral estimates or perturbation argument used to control the Schrödinger operator near the Brillouin-zone vertices; without these details the load-bearing step of the central claim cannot be verified.
Authors: The full manuscript contains the spectral estimates and perturbation arguments in Sections 3 and 4, where the honeycomb symmetries are used to reduce the analysis to the K-points of the Brillouin zone, followed by a controlled perturbation of the untwisted operators at commensurate angles. We acknowledge that the connection between the abstract claim and these sections could be made more explicit and will add a concise outline of the key estimates and argument structure to the introduction in the revised version. revision: yes
-
Referee: [Quantitative bounds paragraph (abstract)] The quantitative bounds on slope flattening are asserted for small potentials; the manuscript must specify whether the constants are uniform over the admissible potential class or depend on additional parameters, as this directly affects the strength of the flattening statement.
Authors: The constants appearing in the flattening bounds are uniform over the admissible class of potentials obeying the stated honeycomb symmetries; they depend only on the potential amplitude and the angular deviation from incommensurability. We will insert an explicit statement to this effect both in the abstract and in the statement of the relevant theorem. revision: yes
Circularity Check
No significant circularity identified
full rationale
The paper is a direct rigorous existence proof for Dirac cones via spectral analysis of the exact continuum Schrödinger operator under explicit symmetry assumptions on the honeycomb potentials (AA/AB stacking) and the definition of commensurate angles. The abstract and claim structure show no reduction of any prediction or result to fitted inputs, self-definitional loops, or load-bearing self-citations; the derivation is self-contained as a mathematical argument within the stated model class.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The potentials are honeycomb lattices with the symmetries of graphene (AA or AB stacking)
- domain assumption The electronic transport is governed by the exact continuum Schrödinger operator
Reference graph
Works this paper leans on
-
[1]
M. J. Ablowitz, C. W Curtis, and Y. Zhu, On tight-binding approximations in optical lattices , Studies in Applied Mathematics 129 (2012), no. 4, 362–388
work page 2012
- [2]
- [3]
- [4]
- [5]
-
[6]
S. Becker and M. Zworski, Dirac points for twisted bilayer graphene with in-plane mag netic field , Journal of Spectral Theory 14 (2024), no. 2, 479–511
work page 2024
-
[7]
G. Berkolaiko and A. Comech, Symmetry and Dirac points in graphene spectrum , Journal of Spectral Theory 8 (2018), no. 3, 1099–1147
work page 2018
-
[8]
R. Bistritzer and A. H MacDonald, Moir´ e bands in twisted double-layer graphene , Proceedings of the National Academy of Sciences 108 (2011), no. 30, 12233–12237
work page 2011
-
[9]
´E. Canc` es, L. Garrigue, and D. Gontier, Simple derivation of Moir´ e-scale continuous models for tw isted bilayer graphene, Physical Review B 107 (2023), no. 15, 155403
work page 2023
-
[10]
Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. K axiras, and P. Jarillo-Herrero, Unconventional superconductivity in magic-angle graphene superlattices , Nature 556 (2018), no. 7699, 43–50
work page 2018
-
[11]
M. Carlsson and O. Rubin, On perturbation of operators and Rayleigh-Schr¨ odinger coefficients, Complex Analysis and Operator Theory 18 (2024), no. 3, 47
work page 2024
-
[12]
G. Catarina, B. Amorim, E. V Castro, J. Lopes, and N. Peres , Twisted bilayer graphene: Low-energy physics, electronic and optical properties , Handbook of Graphene 3 (2019), 177–232
work page 2019
-
[13]
C. Fefferman, J. Lee-Thorp, and M. Weinstein, Honeycomb Schr¨ odinger operators in the strong binding reg ime, Communications on Pure and Applied Mathematics 71 (2018), no. 6, 1178–1270. 34 TAL MALINOVITCH
work page 2018
-
[14]
C. Fefferman and M. Weinstein, Honeycomb lattice potentials and Dirac points , Journal of the American Mathe- matical Society 25 (2012), no. 4, 1169–1220. 15. , Wave packets in honeycomb structures and two-dimensional D irac equations, Communications in Math- ematical Physics 326 (2014), 251–286
work page 2012
-
[15]
Z. Hennighausen and S. Kar, Twistronics: a turning point in 2d quantum materials , Electronic Structure 3 (2021), no. 1, 014004
work page 2021
-
[16]
T. Kong, D. Liu, M. Luskin, and A. Watson, Modeling of electronic dynamics in twisted bilayer graphen e, SIAM Journal on Applied Mathematics 84 (2024), no. 3, 1011–1038
work page 2024
-
[17]
M. Koshino, Band structure and topological properties of twisted doubl e bilayer graphene , Physical Review B 99 (2019), no. 23, 235406
work page 2019
-
[18]
P. Kuchment, An overview of periodic elliptic operators , Bulletin of the American Mathematical Society 53 (2016), 343–414
work page 2016
-
[19]
P. Kuchment and O. Post, On the spectra of carbon nano-structures , Communications in Mathematical Physics 33 (1973), 335–343
work page 1973
-
[20]
P. Ledwith, E. Khalaf, and A. Vishwanath, Strong coupling theory of magic-angle graphene: A pedagogi cal introduction, Annals of Physics 435 (2021), 168646
work page 2021
-
[21]
Y. Li, A. Eaton, H.A. Fertig, and B. Seradjeh, Dirac magic and Lifshitz transitions in AA-stacked twisted multilayer graphene, Physical Review Letters 128 (2022), no. 2, 026404
work page 2022
-
[22]
J. Liu, Z. Ma, J. Gao, and X. Dai, Quantum valley Hall effect, orbital magnetism, and anomalous Hall effect in twisted multilayer graphene systems , Physical Review X 9 (2019), no. 3, 031021
work page 2019
-
[23]
J. Lopes dos Santos, N. Peres, and A. Castro Neto, Continuum model of the twisted graphene bilayer , Physical Review B—Condensed Matter and Materials Physics 86 (2012), no. 15, 155449
work page 2012
-
[24]
D. Massatt, S. Carr, and M. Luskin, Electronic observables for relaxed bilayer 2d heterostruc tures in momentum space, arXiv preprint arXiv:2109.15296 (2021)
-
[25]
D. Massatt, S. Carr, M. Luskin, and C. Ortner, Incommensurate heterostructures in momentum space , Multiscale Modeling & Simulation 16 (2018), no. 1, 429–451
work page 2018
-
[26]
A. H. C. Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, an d A. K. Geim, The electronic properties of graphene , Reviews of modern physics 81 (2009), no. 1, 109
work page 2009
-
[27]
K. S. Novoselov, Nobel lecture: Graphene: Materials in the flatland , Reviews of modern physics 83 (2011), no. 3, 837
work page 2011
-
[28]
K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zha ng, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Electric field effect in atomically thin carbon films , Science 306 (2004), no. 5696, 666–669
work page 2004
-
[29]
C. Nunn, A proof of a generalization of Niven ’s theorem using algebra ic number theory , Rose-Hulman Undergrad- uate Mathematics Journal 22 (2021), no. 2, 3
work page 2021
-
[30]
M. Reed and B. Simon, Methods of modern mathematical physics - IV: Analysis of Ope rators, vol. 4, Elsevier, 1978
work page 1978
-
[31]
M. A. Rodriguez-Andrade, G. Arag´ on-Gonz´ alez, J. L. Ar ag´ on, A. G´ omez-Rodr ´ ıguez, and D. Romeu,The co- incidence site lattices in 2d hexagonal lattices using Cliffo rd algebra, Advances in Applied Clifford Algebras 25 (2015), 425–440
work page 2015
- [32]
-
[33]
J. C. Slonczewski and P. R. Weiss, Band structure of graphite , Physical Review 109 (1958), no. 2, 272
work page 1958
-
[34]
G. Tarnopolsky, A. J. Kruchkov, and A. Vishwanath, Origin of magic angles in twisted bilayer graphene , Physical review letters 122 (2019), no. 10, 106405
work page 2019
-
[35]
P. R. Wallace, The band theory of graphite , Physical review 71 (1947), 622–634
work page 1947
- [36]
-
[37]
A. Watson and M. Luskin, Existence of the first magic angle for the chiral model of bila yer graphene, Journal of Mathematical Physics 62 (2021), no. 9, 091502
work page 2021
-
[38]
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 (2019), no. 6431, 1059–1064
work page 2019
-
[39]
M. Zworski, Mathematical results on the chiral models of twisted bilaye r graphene (with an appendix by Mengxuan Yang and Zhongkai Tao) , Journal of Spectral Theory 14 (2024), no. 3, 1063–1107. TBG IN COMMENSURATE ANGLES 35 Appendix A. Notation • We will denote by ⟨·. ·⟩ the Euclidean inner product on vectors in R2 or Z2, and the size of these vector will...
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.