pith. sign in

arxiv: 2601.16851 · v2 · submitted 2026-01-23 · ❄️ cond-mat.mes-hall · cond-mat.mtrl-sci

Twisted bilayer graphene from first-principles: structural and electronic properties

Albert Zhu , Daniel Bennett , Daniel T. Larson , Mohammed M. Al Ezzi , Efstratios Manousakis , Efthimios Kaxiras This is my paper

Pith reviewed 2026-05-16 11:52 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.mtrl-sci
keywords twisted bilayer graphenedensity functional theorylattice relaxationFermi velocityelectronic band structurecontinuum elastic modelk·p modelcommensurate supercells
0
0 comments X

The pith

First-principles DFT shows twisted bilayer graphene lattice relaxation matches continuum elastic models for all angles down to 1 degree, with Fermi velocity and bandwidth following k·p trends but offset by a small angle shift.

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

The authors perform density functional theory calculations with an optimized local basis set to generate fully relaxed atomic structures and electronic bands for twisted bilayer graphene across a wide range of twist angles. They obtain commensurate supercells down to 0.987 degrees and find that the resulting lattice distortions agree closely with predictions from continuum elastic theory for every angle studied. Where direct comparison is possible, the local-basis results match plane-wave DFT calculations in both geometry and bands. The twist-angle dependence of the Fermi velocity and bandwidth reproduces the shape of results from an exact k·p continuum model but is displaced by a small constant offset in angle, while the low-energy states exhibit clear band inversion and symmetry properties.

Core claim

Using DFT with an optimized local basis in SIESTA, the authors obtain relaxed commensurate supercells for twist angles as small as 0.987° and show that the resulting lattice distortions match continuum elastic theory for every angle examined. Direct comparison with VASP plane-wave calculations at larger angles confirms agreement in both geometry and bands. The Fermi velocity and bandwidth versus twist angle follow the same qualitative trend as the k·p continuum model but exhibit a small offset in the effective twist angle, while low-energy wavefunctions display characteristic band inversion and symmetry properties.

What carries the argument

Density functional theory with an optimized local basis set applied to fully relaxed commensurate supercells of twisted bilayer graphene.

If this is right

  • The calculations supply a benchmark for continuum elastic models of structural relaxation in tBLG.
  • Local-basis and plane-wave DFT agree closely on both atomic positions and electronic bands where both methods apply.
  • The Fermi velocity and bandwidth follow the same qualitative twist-angle dependence as the k·p model but shifted by a small angle offset.
  • Low-energy wavefunctions show band inversion and specific symmetries that can be used to interpret spectroscopic probes.
  • The results provide an ab initio reference structure and band data to anchor future calculations that include many-body effects.

Where Pith is reading between the lines

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

  • If the small angle offset persists under improved functionals, continuum models may require a single adjusted parameter to achieve quantitative agreement with first-principles data.
  • The relaxation patterns could be inserted into models of strain-driven superconductivity to test whether the offset affects the location of flat bands.
  • The symmetry labels on the low-energy states suggest selection rules that could be checked in tunneling or optical experiments on devices with controlled twist angles.
  • Extending the same local-basis approach to include explicit van der Waals corrections might reduce or eliminate the reported offset without changing the structural conclusions.

Load-bearing premise

The exchange-correlation functional and local basis set chosen for the calculations accurately capture both the structural relaxation and the low-energy electronic bands of twisted bilayer graphene.

What would settle it

An experimental measurement of Fermi velocity versus twist angle that deviates from the reported DFT dependence by more than the claimed small offset would falsify the quantitative accuracy of the electronic structure results.

Figures

Figures reproduced from arXiv: 2601.16851 by Albert Zhu, Daniel Bennett, Daniel T. Larson, Efstratios Manousakis, Efthimios Kaxiras, Mohammed M. Al Ezzi.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the atomic structure and forces from DFT calculations using [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (b) shows the layer separation d, defined as the dif￾ference in the out-of-plane coordinates between the top and bottom layers after geometry relaxation. Most of the corruga￾tion occurs around the AA sites. For comparison to the DFT results, we performed relaxation calculations using a contin￾uum model employing the concept of the generalized stacking fault energy (GSFE) [62]. For completeness, a derivatio… view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The absolute value [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Electronic band structure of tBLG for commensurate twist angles near the magic angle. The black lines show the bands from DFT [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Probability density of the [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Electronic band structures and wavefunctions as a function of twist angle [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
read the original abstract

We present a comprehensive first-principles study of twisted bilayer graphene (tBLG) for a wide range of twist angles, with a focus on structural and electronic properties. By employing density functional theory (DFT) with an optimized local basis set, we simulate tBLG, obtaining fully relaxed commensurate structures for twist angles down to 0.987{\deg}. For all angles the lattice relaxation agrees well with continuum elastic models. For angles accessible to plane-wave DFT (VASP), we provide a detailed comparison with our local basis DFT (SIESTA) calculations, demonstrating excellent agreement in both the atomic and electronic structure. The dependence of the Fermi velocity and band width on the twist angle shows qualitative agreement with results from an `exact' $\mathbf{k \cdot p}$ continuum model, but reveals a small twist angle offset. Additionally, we provide details of the low-energy wavefunction character, band inversion and symmetries. Our results provide an ab initio reference point for the microscopic structure and electronic properties of tBLG which will serve as the foundation for future studies incorporating many-body effects.

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

0 major / 2 minor

Summary. The manuscript presents a first-principles DFT study of twisted bilayer graphene using SIESTA with an optimized local basis set. It obtains fully relaxed commensurate structures down to 0.987°, demonstrates agreement of lattice relaxations with continuum elastic models for all angles, shows excellent numerical agreement with VASP calculations on both structure and bands for angles accessible to plane-wave methods, and reports qualitative agreement with an exact k·p continuum model for the twist-angle dependence of Fermi velocity and bandwidth, accompanied by a small observed offset. Additional results cover low-energy wavefunction character, band inversion, and symmetries, positioning the work as an ab initio reference for future many-body studies.

Significance. If the results hold, the work supplies a validated ab initio reference dataset for tBLG structural and electronic properties across a wide angle range, with direct cross-code validation between independent DFT implementations (SIESTA and VASP) that strengthens in the local-basis approach for large supercells. Explicit agreement with continuum elastic models for relaxations down to ~1° and the presentation of the k·p offset as an observation rather than a fitted parameter add to its utility as a benchmark.

minor comments (2)
  1. [§2] §2 (Methods): The optimization procedure for the local basis set is described only at a high level; adding a short table or paragraph with the final basis parameters (e.g., cutoff radii, polarization functions) would improve reproducibility.
  2. [Figure 4] Figure 4 caption: The definition of the plotted bandwidth (full width or half-width) is not stated explicitly; a one-sentence clarification would remove ambiguity when comparing to the k·p results.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive and constructive report, which accurately summarizes the scope and contributions of our work. We are pleased that the manuscript is recommended for acceptance.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper reports independent first-principles DFT results (SIESTA relaxations and bands for tBLG down to 0.987°, cross-validated against VASP where feasible) and direct numerical comparisons to continuum elastic models and a k·p model. No central quantity is obtained by fitting a parameter to the same data being reported, no self-definitional equations appear, and no load-bearing premise reduces to a self-citation chain. The noted small twist-angle offset versus the continuum model is an observed discrepancy, not a fitted input. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central results rest on standard DFT approximations whose accuracy for tBLG is assumed rather than re-derived; no new entities are introduced and the only adjustable elements are the usual convergence parameters of the chosen basis and functional.

axioms (1)
  • domain assumption Standard DFT exchange-correlation functional and pseudopotentials are sufficiently accurate for structural relaxation and low-energy bands in tBLG
    Invoked throughout the comparisons with continuum models and between SIESTA and VASP

pith-pipeline@v0.9.0 · 5515 in / 1263 out tokens · 20524 ms · 2026-05-16T11:52:10.724134+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

67 extracted references · 67 canonical work pages

  1. [1]

    Y . Cao, V . Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxi- ras, and P. Jarillo-Herrero, Nature556, 43 (2018)

    work page 2018

  2. [2]

    X. Lu, P. Stepanov, W. Yang, M. Xie, M. A. Aamir, I. Das, C. Urgell, K. Watanabe, T. Taniguchi, G. Zhang,et al., Nature 574, 653 (2019)

    work page 2019

  3. [3]

    Yankowitz, S

    M. Yankowitz, S. Chen, H. Polshyn, Y . Zhang, K. Watanabe, T. Taniguchi, D. Graf, A. F. Young, and C. R. Dean, Science 363, 1059 (2019)

    work page 2019

  4. [4]

    M. Oh, K. P. Nuckolls, D. Wong, R. L. Lee, X. Liu, K. Watan- abe, T. Taniguchi, and A. Yazdani, Nature600, 240 (2021)

    work page 2021

  5. [5]

    Tanaka, J

    M. Tanaka, J. ˆI.-j. Wang, T. H. Dinh, D. Rodan-Legrain, S. Za- man, M. Hays, A. Almanakly, B. Kannan, D. K. Kim, B. M. Niedzielski,et al., Nature638, 99 (2025)

    work page 2025

  6. [6]

    Banerjee, Z

    A. Banerjee, Z. Hao, M. Kreidel, P. Ledwith, I. Phinney, J. M. Park, A. Zimmerman, M. E. Wesson, K. Watanabe, T. Taniguchi,et al., Nature638, 93 (2025)

    work page 2025

  7. [7]

    Y . Cao, V . Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y . Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxi- ras,et al., Nature556, 80 (2018)

    work page 2018

  8. [8]

    Codecido, Q

    E. Codecido, Q. Wang, R. Koester, S. Che, H. Tian, R. Lv, S. Tran, K. Watanabe, T. Taniguchi, F. Zhang,et al., Sci. Adv. 5, eaaw9770 (2019)

    work page 2019

  9. [9]

    Inbar, J

    A. Inbar, J. Birkbeck, J. Xiao, T. Taniguchi, K. Watanabe, B. Yan, Y . Oreg, A. Stern, E. Berg, and S. Ilani, Nature614, 11 682 (2023)

    work page 2023

  10. [10]

    K. P. Nuckolls, R. L. Lee, M. Oh, D. Wong, T. Soejima, J. P. Hong, D. Cualuguaru, J. Herzog-Arbeitman, B. A. Bernevig, K. Watanabe,et al., Nature620, 525 (2023)

    work page 2023

  11. [11]

    H. Tian, E. Codecido, D. Mao, K. Zhang, S. Che, K. Watanabe, T. Taniguchi, D. Smirnov, E.-A. Kim, M. Bockrath,et al., Nat. Phys.20, 1407 (2024)

    work page 2024

  12. [12]

    Krishna Kumar, G

    R. Krishna Kumar, G. Li, R. Bertini, S. Chaudhary, K. Nowakowski, J. M. Park, S. Castilla, Z. Zhan, P. A. Pan- tale´on, H. Agarwal,et al., Nat. Mater. , 1 (2025)

    work page 2025

  13. [13]

    A. L. Sharpe, E. J. Fox, A. W. Barnard, J. Finney, K. Watanabe, T. Taniguchi, M. Kastner, and D. Goldhaber-Gordon, Science 365, 605 (2019)

    work page 2019

  14. [14]

    Serlin, C

    M. Serlin, C. Tschirhart, H. Polshyn, Y . Zhang, J. Zhu, K. Watanabe, T. Taniguchi, L. Balents, and A. Young, Science 367, 900 (2020)

    work page 2020

  15. [15]

    K. P. Nuckolls, M. Oh, D. Wong, B. Lian, K. Watanabe, T. Taniguchi, B. A. Bernevig, and A. Yazdani, Nature588, 610 (2020)

    work page 2020

  16. [16]

    S. Wu, Z. Zhang, K. Watanabe, T. Taniguchi, and E. Y . Andrei, Nat. Mater.20, 488 (2021)

    work page 2021

  17. [17]

    Y . Choi, H. Kim, Y . Peng, A. Thomson, C. Lewandowski, R. Polski, Y . Zhang, H. S. Arora, K. Watanabe, T. Taniguchi, et al., Nature589, 536 (2021)

    work page 2021

  18. [18]

    Stepanov, M

    P. Stepanov, M. Xie, T. Taniguchi, K. Watanabe, X. Lu, A. H. MacDonald, B. A. Bernevig, and D. K. Efetov, Phys. Rev. Lett. 127, 197701 (2021)

    work page 2021

  19. [19]

    Y . Xie, A. T. Pierce, J. M. Park, D. E. Parker, E. Khalaf, P. Led- with, Y . Cao, S. H. Lee, S. Chen, P. R. Forrester,et al., Nature 600, 439 (2021)

    work page 2021

  20. [20]

    J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Cas- tro Neto, Phys. Rev. Lett.99, 256802 (2007)

    work page 2007

  21. [21]

    Bistritzer and A

    R. Bistritzer and A. H. MacDonald, Proc. Natl. Acad. Sci.108, 12233 (2011)

    work page 2011

  22. [22]

    E. J. Mele, Phys. Rev. B84, 235439 (2011)

    work page 2011

  23. [23]

    J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Cas- tro Neto, Phys. Rev. B86, 155449 (2012)

    work page 2012

  24. [24]

    S. Carr, S. Fang, Z. Zhu, and E. Kaxiras, Phys. Rev. Res.1, 013001 (2019)

    work page 2019

  25. [25]

    S. Fang, S. Carr, Z. Zhu, D. Massatt, and E. Kaxiras, arXiv:1908.00058 (2019)

    work page arXiv 1908

  26. [26]

    Guinea and N

    F. Guinea and N. R. Walet, Phys. Rev. B99, 205134 (2019)

    work page 2019

  27. [27]

    Kang and O

    J. Kang and O. Vafek, Phys. Rev. B107, 075408 (2023)

    work page 2023

  28. [28]

    Vafek and J

    O. Vafek and J. Kang, Phys. Rev. B107, 075123 (2023)

    work page 2023

  29. [29]

    W. Miao, C. Li, X. Han, D. Pan, and X. Dai, Phys. Rev. B107, 125112 (2023)

    work page 2023

  30. [30]

    M. M. A. Ezzi, G. N. Pallewela, C. De Beule, E. J. Mele, and S. Adam, Phys. Rev. Lett.133, 266201 (2024)

    work page 2024

  31. [31]

    Trambly de Laissardi `ere, D

    G. Trambly de Laissardi `ere, D. Mayou, and L. Magaud, Nano Lett.10, 804 (2010)

    work page 2010

  32. [32]

    Jung and A

    J. Jung and A. H. MacDonald, Phys. Rev. B89, 035405 (2014)

    work page 2014

  33. [33]

    Fang and E

    S. Fang and E. Kaxiras, Phys. Rev. B93, 235153 (2016)

    work page 2016

  34. [34]

    S. Carr, D. Massatt, S. Fang, P. Cazeaux, M. Luskin, and E. Kaxiras, Phys. Rev. B95, 075420 (2017)

    work page 2017

  35. [35]

    L. Zou, H. C. Po, A. Vishwanath, and T. Senthil, Phys. Rev. B 98, 085435 (2018)

    work page 2018

  36. [36]

    H. C. Po, L. Zou, A. Vishwanath, and T. Senthil, Phys. Rev. X 8, 031089 (2018)

    work page 2018

  37. [37]

    Koshino, N

    M. Koshino, N. F. Q. Yuan, T. Koretsune, M. Ochi, K. Kuroki, and L. Fu, Phys. Rev. X8, 031087 (2018)

    work page 2018

  38. [38]

    Kang and O

    J. Kang and O. Vafek, Phys. Rev. X8, 031088 (2018)

    work page 2018

  39. [39]

    H. C. Po, L. Zou, T. Senthil, and A. Vishwanath, Phys. Rev. B 99, 195455 (2019)

    work page 2019

  40. [40]

    Tarnopolsky, A

    G. Tarnopolsky, A. J. Kruchkov, and A. Vishwanath, Phys. Rev. Lett.122, 106405 (2019)

    work page 2019

  41. [41]

    S. Carr, S. Fang, H. C. Po, A. Vishwanath, and E. Kaxiras, Phys. Rev. Res.1, 033072 (2019)

    work page 2019

  42. [42]

    B. A. Bernevig, Z.-D. Song, N. Regnault, and B. Lian, Phys. Rev. B103, 205411 (2021)

    work page 2021

  43. [43]

    Song and B

    Z.-D. Song and B. A. Bernevig, Phys. Rev. Lett.129, 047601 (2022)

    work page 2022

  44. [44]

    Bennett, D

    D. Bennett, D. T. Larson, L. Sharma, S. Carr, and E. Kaxiras, Phys. Rev. B109, 155422 (2024)

    work page 2024

  45. [45]

    Uchida, S

    K. Uchida, S. Furuya, J.-I. Iwata, and A. Oshiyama, Phys. Rev. B90, 155451 (2014)

    work page 2014

  46. [46]

    Lucignano, D

    P. Lucignano, D. Alf`e, V . Cataudella, D. Ninno, and G. Cantele, Phys. Rev. B99, 195419 (2019)

    work page 2019

  47. [47]

    Yananose, G

    K. Yananose, G. Cantele, P. Lucignano, S.-W. Cheong, J. Yu, and A. Stroppa, Phys. Rev. B104, 075407 (2021)

    work page 2021

  48. [48]

    M. M. A. Ezzi, A. Zhu, D. Bennett, D. T. Larson, and E. Kaxi- ras, arXiv:2505.19485 (2025)

    work page arXiv 2025

  49. [49]

    Bennett, M

    D. Bennett, M. Pizzochero, J. Junquera, and E. Kaxiras, Phys. Rev. B111, 125123 (2025)

    work page 2025

  50. [50]

    A. Zhu, D. Bennett, D. T. Larson, M. M. A. Ezzi, E. Manousakis, and E. Kaxiras, Phys. Rev. B113, L041112 (2026)

    work page 2026

  51. [51]

    J. M. Soler, E. Artacho, J. D. Gale, A. Garc ´ıa, J. Junquera, P. Ordej´on, and D. S ´anchez-Portal, J. Phys. Condens. Matter 14, 2745 (2002)

    work page 2002

  52. [52]

    Garc ´ıa, N

    A. Garc ´ıa, N. Papior, A. Akhtar, E. Artacho, V . Blum, E. Bosoni, P. Brandimarte, M. Brandbyge, J. I. Cerd ´a, F. Corsetti,et al., J. Chem. Phys.152(2020)

    work page 2020

  53. [53]

    Kresse and J

    G. Kresse and J. Furthm ¨uller, Phys. Rev. B54, 11169 (1996)

    work page 1996

  54. [54]

    Kresse and J

    G. Kresse and J. Furthm ¨uller, Comput Mater. Sci.6, 15 (1996)

    work page 1996

  55. [55]

    S. Carr, D. Massatt, S. B. Torrisi, P. Cazeaux, M. Luskin, and E. Kaxiras, Phys. Rev. B98, 224102 (2018)

    work page 2018

  56. [56]

    Hamann, Phys

    D. Hamann, Phys. Rev. B88, 085117 (2013)

    work page 2013

  57. [57]

    Garc ´ıa, M

    A. Garc ´ıa, M. J. Verstraete, Y . Pouillon, and J. Junquera, Com- put. Phys. Commun.227, 51 (2018)

    work page 2018

  58. [58]

    Van Setten, M

    M. Van Setten, M. Giantomassi, E. Bousquet, M. J. Verstraete, D. R. Hamann, X. Gonze, and G.-M. Rignanese, Comput. Phys. Commun.226, 39 (2018)

    work page 2018

  59. [59]

    Junquera, O

    J. Junquera, O. Paz, D. S ´anchez-Portal, and E. Artacho, Phys. Rev. B64, 235111 (2001)

    work page 2001

  60. [60]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)

    work page 1996

  61. [61]

    Grimme, J

    S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, J. Chem. Phys. 132(2010)

    work page 2010

  62. [62]

    Kaxiras and M

    E. Kaxiras and M. S. Duesbery, Phys. Rev. Lett.70, 3752 (1993)

    work page 1993

  63. [63]

    Bennett and B

    D. Bennett and B. Remez, npj 2D Mater. Appl.6, 1 (2022)

    work page 2022

  64. [64]

    Bennett, Phys

    D. Bennett, Phys. Rev. B105, 235445 (2022)

    work page 2022

  65. [65]

    Ramos-Alonso, B

    A. Ramos-Alonso, B. Remez, D. Bennett, R. M. Fernandes, and H. Ochoa, Phys. Rev. Lett.134, 026501 (2025)

    work page 2025

  66. [66]

    J. Ahn, S. Park, and B.-J. Yang, Phys. Rev. X9, 021013 (2019)

    work page 2019

  67. [67]

    Z. Song, Z. Wang, W. Shi, G. Li, C. Fang, and B. A. Bernevig, Phys. Rev. Lett.123, 036401 (2019)

    work page 2019