pith. sign in

arxiv: 2606.13434 · v1 · pith:AUEHNMPInew · submitted 2026-06-11 · 🧮 math.NA · cond-mat.mtrl-sci· cs.NA

Momentum Space Algorithm for Electronic Structure of Double-Incommensurate Trilayer Graphene

Pith reviewed 2026-06-27 05:51 UTC · model grok-4.3

classification 🧮 math.NA cond-mat.mtrl-scics.NA
keywords momentum space algorithmdouble-incommensurate trilayer graphenetight-binding modeldensity of statestruncation schememagic angleselectronic structure
0
0 comments X

The pith

A momentum space algorithm with four-dimensional lattice truncation computes electronic structure for double-incommensurate trilayer graphene with improved convergence.

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

The paper introduces a momentum space framework that exactly transforms the tight-binding model for double-incommensurate trilayer graphene into a momenta description. It adds an efficient truncation scheme for the resulting four-dimensional lattice to achieve faster numerical convergence of the density of states and momentum local density of states. The method is applied to an ab initio model of twisted trilayer graphene, with convergence validated and altered band behavior near flat bands at magic angles observed. This approach is higher order than continuum models because it preserves the exact tight-binding structure.

Core claim

The momentum space algorithm for double-incommensurate trilayers is an exact transformation of the tight-binding model into a four-dimensional lattice space; the introduced truncation scheme improves convergence of the density of states and momentum local density of states without significant errors, and numerical tests on twisted trilayer graphene show it captures altered band behavior near flat bands at magic angles.

What carries the argument

The efficient truncation scheme of the four-dimensional lattice, which approximates the infinite momentum lattice while preserving essential physics of the tight-binding model.

If this is right

  • Density of states and momentum local density of states converge faster than with prior continuum descriptions.
  • The method reveals altered band behavior near flat bands at magic angles due to its higher-order accuracy.
  • Convergence estimates hold for the implemented ab initio model of twisted trilayer graphene.
  • The parallel structure to classical band structure enables direct comparison with experimental observables.

Where Pith is reading between the lines

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

  • The truncation approach could be adapted to compute other observables such as conductivity or optical responses in similar incommensurate systems.
  • Higher accuracy near magic angles may help model interaction-driven effects like superconductivity when combined with many-body methods.
  • The four-dimensional lattice representation might allow systematic error bounds that are harder to obtain in real-space tight-binding calculations.

Load-bearing premise

The truncation scheme of the four-dimensional lattice improves convergence without introducing significant errors or losing essential physics of the original tight-binding model.

What would settle it

A computation with successively finer truncations that shows the density of states or local density of states near the magic angles changing by more than the reported convergence threshold would falsify the accuracy claim.

Figures

Figures reproduced from arXiv: 2606.13434 by Daniel Massatt, Ken Beard.

Figure 1
Figure 1. Figure 1: Atomic configuration of twisted trilayer graphene for [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Fixing 𝐿 = 20 and 𝜃 = [− √ 3 2 , 0, √ 2 2 ], increasing 𝑊 results in movement along the 𝑧-axis of the monolayer Dirac cones. In (a) we give a side-view of all three monolayer Dirac cones as 𝑊 increases from 0.15 (green) to 0.35 (orange). In (b) we give a top-down view of all three monolayer dirac cones as 𝑊 increases from 0.15 (green) to 0.35 (orange). Next, we define the projection operator, which will pr… view at source ↗
Figure 3
Figure 3. Figure 3: For fixed 𝑊 = 0.45 and 𝜃 = [−1.4, 0, 2.8], we demonstrate the filling of the moiré mapped 𝒲∗ 𝑗 (𝑞, 𝑊), that is, (𝐵𝑘 − 𝐵𝑗 )𝐵 −1 𝑘 𝐺𝑘 + (𝐵𝑙 − 𝐵𝑗 )𝐵 −1 𝑙 𝐺𝑙 for all 𝐺 ∈ 𝒲∗ (𝑞, 𝑊) ∩ ℒ∗ (𝐿). With regard to the relative error introduced by the 𝐿-truncation, we obtain the following lemma. Lemma 4.4. For double-incommensurate TTG, the following bound holds |𝒟ˆ (𝑁,𝜏,𝑊) 𝜖 (𝐸) − 𝒟ˆ (𝑁,𝜏,𝑊,𝐿) 𝜖 (𝐸)| ≲ 𝜖−3𝛾 −1 𝜖 𝐿𝑒−𝛾𝜖𝐿… view at source ↗
Figure 4
Figure 4. Figure 4: For 𝜃 = [− √ 3 2 , 0, √ 2 2 ] degrees, 𝑊 = 0.3771, 𝐿 = 30.5555 and 𝑃 = 8000, we compare the continuum (a) and tight-binding models (b). We observe that concentrations of the local density of states in the continuum model are lost in the tight-binding model. We now present several momentum LDoS plots for both our tight-binding model and the continuum model (BM) proposed in [28] with all energy units given i… view at source ↗
Figure 5
Figure 5. Figure 5: For 𝜃 = [−1.4, 0, 2.8] degrees, 𝑊 = 0.47422, 𝐿 = 19.40479, 𝑃 = 8000, we compare the continuum (a) and tight-binding models (b). We observe that the mostly flat band near the Fermi energy in the continuum model, has significantly changed in the tight-binding model. While this may still be sufficient for a superconducting phase, it suggests that the flat band is not the core mechanism. interest. We begin by … view at source ↗
Figure 6
Figure 6. Figure 6: Semilogy of the relative 𝑊-error for a range of Chebyshev order 𝑃: (a) 𝜃 = [− √ 3 2 , 0, √ 2 2 ] at 𝐿 = 35.2831, (b) 𝜃 = [−1.4, 0, 2.8] at 𝐿 = 51.5415. Semilogy of the relative 𝐿-error for a range of Chebyshev orders 𝑃: (c) 𝜃 = [− √ 3 2 , 0, √ 2 2 ] at 𝑊 = 0.3771, (d) 𝜃 = [−1.4, 0, 2.8] at 𝑊 = 0.47422. Exponential convergence is observed for both angles. 54 [PITH_FULL_IMAGE:figures/full_fig_p054_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Two semilogy plots of the relative 𝑁-error for 𝑃 = 4000: (a) 𝜃 = [− √ 3 2 , 0, √ 2 2 ] at 𝑊 = 0.1743 and 𝐿 = 73.7463, (b) 𝜃 = [−1.4, 0, 2.8] at 𝑊 = 0.4324 and 𝐿 = 73.7463. Exponential convergence is observed for both angles. (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p055_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Two DoS plots for 𝑃 = 4000 and 𝑁 = 30: (a) 𝜃 = [− √ 3 2 , 0, √ 2 2 ] at 𝑊 = 0.1743 and 𝐿 = 73.7463, (b) 𝜃 = [−1.4, 0, 2.8] at 𝑊 = 0.4324 and 𝐿 = 73.7463. 55 [PITH_FULL_IMAGE:figures/full_fig_p055_8.png] view at source ↗
read the original abstract

Numerical algorithms for computing electronic structure of incommensurate 2D materials using ab initio models is critical for predicting material properties and guiding experiment. For bilayers, momentum space and continuum models have been introduced to approximate observables of ab initio tight-binding models using a momenta description despite the lack of periodicity in the tight-binding model required for Bloch theory. A similar structure has been introduced for double-incommensurate trilayers using a continuum model, where the three lattices are all mutually incommensurate. However, this description leads to a four-dimensional lattice space, and numerical convergence of the density of states was observed to have poor convergence. In this work, we introduce a momentum space framework for double incommensurate trilayer graphene, and introduce an efficient truncation scheme of the four-dimensional lattice to drastically improve convergence of the density of states and momentum local density of states (a parallel object to classical band structure). We implement this algorithm on an ab initio model of twisted trilayer graphene and validate convergence estimates. We further verify numerically that the momentum space algorithm, inherently higher order than the continuum model as it is an exact transformation of the tight-binding model, captures altered band behavior near the flat bands at magic angles.

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

2 major / 0 minor

Summary. The paper introduces a momentum-space framework for the electronic structure of double-incommensurate trilayer graphene. It replaces the continuum model with an exact transformation of the ab initio tight-binding Hamiltonian, yielding a four-dimensional momentum lattice, and proposes an efficient truncation scheme that is asserted to restore rapid convergence of the density of states and momentum-resolved local density of states while preserving the essential physics near magic-angle flat bands.

Significance. If the truncation is shown to be faithful, the algorithm would supply a systematically improvable, higher-order alternative to continuum models for incommensurate trilayers, enabling reliable numerical access to band reconstructions that are currently inaccessible.

major comments (2)
  1. [Abstract] Abstract: the assertion that the truncation 'drastically improves convergence … without introducing significant errors or losing essential physics' rests solely on observed numerical convergence of integrated quantities (DOS, MLDOS). No a-priori bound is supplied on the neglected interlayer matrix elements, which decay only polynomially; convergence of an integrated observable does not guarantee that eigenvalue rearrangements inside the flat-band window are preserved.
  2. [Abstract] The central claim that the momentum-space algorithm is 'inherently higher order than the continuum model' and 'an exact transformation' is load-bearing, yet the manuscript provides no explicit error estimate or comparison that quantifies the order improvement once the four-dimensional lattice is truncated.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments. We address the two major comments point by point below, indicating the revisions we intend to make to the manuscript.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the assertion that the truncation 'drastically improves convergence … without introducing significant errors or losing essential physics' rests solely on observed numerical convergence of integrated quantities (DOS, MLDOS). No a-priori bound is supplied on the neglected interlayer matrix elements, which decay only polynomially; convergence of an integrated observable does not guarantee that eigenvalue rearrangements inside the flat-band window are preserved.

    Authors: We agree that the manuscript relies on numerical evidence of convergence for the DOS and MLDOS (including momentum-resolved quantities near the flat bands) rather than an a-priori bound on the truncation error for the polynomially decaying interlayer couplings. Convergence of integrated observables does not rigorously guarantee preservation of all eigenvalue rearrangements. We will revise the abstract to moderate the language, removing the claim of 'without introducing significant errors or losing essential physics' and instead stating that the truncation is validated by observed numerical convergence. We will also add a short paragraph in the main text acknowledging this limitation and noting that a rigorous a-priori error analysis for the flat-band window remains future work. revision: yes

  2. Referee: [Abstract] The central claim that the momentum-space algorithm is 'inherently higher order than the continuum model' and 'an exact transformation' is load-bearing, yet the manuscript provides no explicit error estimate or comparison that quantifies the order improvement once the four-dimensional lattice is truncated.

    Authors: The phrase 'exact transformation' in the abstract refers specifically to the untruncated momentum-space formulation, which is derived directly from the ab initio tight-binding Hamiltonian without the approximations inherent to continuum models (such as dispersion linearization). This establishes the base algorithm as higher-order in fidelity to the tight-binding model. The truncation is presented as a separate numerical approximation whose practical effect is demonstrated through convergence studies. We acknowledge that the manuscript does not supply an explicit error estimate or quantitative comparison of the truncated four-dimensional lattice versus continuum models. We will revise the abstract to clarify this distinction and avoid any implication of a quantified order improvement for the truncated version. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation self-contained via exact transformation and numerical validation

full rationale

The paper defines its momentum-space algorithm as an exact transformation of the tight-binding model and validates the truncation scheme through direct numerical convergence of DOS and MLDOS on an ab initio model of twisted trilayer graphene. No quoted step reduces a claimed prediction or uniqueness result to a fitted parameter, self-citation chain, or ansatz defined by the authors' own prior work; the central claims rest on the stated transformation property and external numerical checks rather than internal redefinition.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; the central claim rests on the unelaborated assertion that the momentum-space map is exact and that the truncation preserves observables.

pith-pipeline@v0.9.1-grok · 5752 in / 1122 out tokens · 32340 ms · 2026-06-27T05:51:18.940759+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

31 extracted references · 25 canonical work pages · 2 internal anchors

  1. [1]

    Electronic spectral properties of incommensurate twisted trilayer graphene

    Amorim, B., Castro, E.V., 2018. Electronic spectral prop- erties of incommensurate twisted trilayer graphene. URL: https://arxiv.org/abs/1807.11909, doi:10.48550/ARXIV.1807.11909

  2. [2]

    Exact ground state of interact- ing electrons in magic angle graphene

    Becker, S., Lin, L., Stubbs, K.D., 2025. Exact ground state of interact- ing electrons in magic angle graphene. Communications in Mathematical Physics 406, 148. URL:https://doi.org/10.1007/s00220-025-05300-x, doi:10.1007/s00220-025-05300-x

  3. [3]

    Bistritzer \ and\ author A

    Bistritzer, R., MacDonald, A.H., 2011. Moiré bands in twisted double-layer graphene. Proceedings of the National Academy of Sciences of the United States of America 108, 12233–12237. doi:10.1073/pnas.1108174108

  4. [4]

    Generalized Kubo formulas for the transport properties of incommensurate 2D atomic heterostructures

    Cancès, E., Cazeaux, P., Luskin, M., 2017. Generalized Kubo formulas for the transport properties of incommensurate 2D atomic heterostructures. Journal of Mathematical Physics 58, 063–502. doi:10.1063/1.4984041. 50

  5. [5]

    Semiclassical analysis of two- scale electronic hamiltonians for twisted bilayer graphene

    Cancès, E., Meng, L., 2023. Semiclassical analysis of two- scale electronic hamiltonians for twisted bilayer graphene. URL: https://arxiv.org/abs/2311.14011,arXiv:2311.14011

  6. [6]

    Cao , author V

    Cao, Y., Fatemi, V., Fang, S., Watanabe, K., Taniguchi, T., Kaxi- ras, E., Jarillo-Herrero, P., 2018. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50. URL: https://doi.org/10.1038/nature26160, doi:10.1038/nature26160

  7. [7]

    QM/MM Methods for Crystalline Defects

    Chen, H., Ortner, C., 2016. QM/MM Methods for Crystalline Defects. Part 1: Locality of the Tight Binding Model. Multiscale Modeling & Simulation 14, 232–264. URL: http://epubs.siam.org/doi/10.1137/15M1022628, doi:10.1137/15M1022628

  8. [8]

    Electronic structure theory of weakly interacting bilayers

    Fang, S., Kaxiras, E., 2016. Electronic structure theory of weakly interacting bilayers. Physical Review B 93, 235153. URL: https://link.aps.org/doi/10.1103/PhysRevB.93.235153, doi:10.1103/PhysRevB.93.235153

  9. [9]

    Guo , author J

    Guo, Y., Pack, J., Swann, J., Holtzman, L., Cothrine, M., Watanabe, K., Taniguchi, T., Mandrus, D.G., Barmak, K., Hone, J., Millis, A.J., Pasupathy, A., Dean, C.R., 2025. Superconductivity in 5.0◦twisted bilayer wse2. Nature 637, 839–845. URL: https://doi.org/10.1038/s41586-024-08381-1, doi:10.1038/s41586-024-08381-1

  10. [10]

    Tight-binding model for graphene𝜋-bands from maximally localized Wannier functions

    Jung, J., MacDonald, A.H., 2013. Tight-binding model for graphene𝜋-bands from maximally localized Wannier functions. Physical Review B 87, 195450. URL: https://link.aps.org/doi/10.1103/PhysRevB.87.195450, doi:10.1103/PhysRevB.87.195450

  11. [11]

    Evidence for unconventional superconductivity in twisted trilayer graphene

    Kim, H., Choi, Y., Lewandowski, C., Thomson, A., Zhang, Y., Polski, R., Watanabe, K., Taniguchi, T., Alicea, J., Nadj-Perge, S., 2022. Evidence for unconventional superconductivity in twisted trilayer graphene. Nature 606, 494–500. URL: https://doi.org/10.1038/s41586-022-04715-z, doi:10.1038/s41586-022-04715-z

  12. [12]

    Numerical Analysis

    Kress, R., 1998. Numerical Analysis. volume 181 ofGrad- uate Texts in Mathematics. Springer New York, New York, NY. URL: http://link.springer.com/10.1007/978-1-4612-0599-9, doi:10.1007/978-1-4612-0599-9

  13. [13]

    Fractional chern insulator states in twisted bilayer graphene: An analytical approach

    Ledwith, P.J., Tarnopolsky, G., Khalaf, E., Vishwanath, A., 2020. Fractional chern insulator states in twisted bilayer graphene: An analytical approach. Phys. Rev. Res. 2, 023237. URL: https://link.aps.org/doi/10.1103/PhysRevResearch.2.023237, doi:10.1103/PhysRevResearch.2.023237

  14. [14]

    Quantum phases in twisted homobilayer transition metal dichalcogenides

    Li, B., Qiu, W.X., Wu, F., MacDonald, A.H., 2025. Quantum phases in twisted homobilayer transition metal dichalcogenides. URL: https://arxiv.org/abs/2509.07360,arXiv:2509.07360. 51

  15. [15]

    Electronic Observables for Relaxed Bilayer Two-Dimensional Heterostructures in Momentum Space

    Massatt, D., Carr, S., Luskin, M., 2023. Electronic Observables for Relaxed Bilayer Two-Dimensional Heterostructures in Momentum Space. Multiscale Modeling & Simulation 21, 1344–1378. doi:10.1137/21M1451208

  16. [16]

    Incommensurate Heterostructures in Momentum Space

    Massatt, D., Carr, S., Luskin, M., Ortner, C., 2018. Incommensurate Heterostructures in Momentum Space. Multiscale Modeling & Simulation 16, 429–451. URL:https://epubs.siam.org/doi/10.1137/17M1141035, doi:10.1137/17M1141035

  17. [17]

    Commensurate and incommensurate double moiré interference in twisted trilayer graphene

    Meng, H., Zhan, Z., Yuan, S., 2023. Commensurate and incommensurate double moiré interference in twisted trilayer graphene. Phys. Rev. B 107, 035109. URL: https://link.aps.org/doi/10.1103/PhysRevB.107.035109, doi:10.1103/PhysRevB.107.035109

  18. [18]

    Flatbands and per- fect metal in trilayer moiré graphene

    Mora, C., Regnault, N., Bernevig, B.A., 2019. Flatbands and per- fect metal in trilayer moiré graphene. Phys. Rev. Lett. 123, 026402. URL: https://link.aps.org/doi/10.1103/PhysRevLett.123.026402, doi:10.1103/PhysRevLett.123.026402

  19. [19]

    Construction and accuracy of electronic continuum models of incommensurate bi- layer 2d materials

    Quan, X., Watson, A.B., Massatt, D., 2025. Construction and accuracy of electronic continuum models of incommensurate bi- layer 2d materials. Electronic Structure 7, 035005. URL: https://doi.org/10.1088/2516-1075/ae0545, doi:10.1088/2516- 1075/ae0545

  20. [20]

    Re- view: moiré-of-moiré superlattice in twisted trilayer graphene

    Ren, W., Zhu, Z., Zhang, X., Luskin, M., Wang, K., 2025. Re- view: moiré-of-moiré superlattice in twisted trilayer graphene. Journal of Physics: Condensed Matter 37, 353001. URL: https://doi.org/10.1088/1361-648X/adf6f9, doi:10.1088/1361- 648X/adf6f9

  21. [21]

    Fourier Analysis, Self-Adjointness

    Simon, B., Reed, M., 2007. Fourier Analysis, Self-Adjointness. volume 2 of Methods of Modern Mathematical Physics. 1 ed., Acad. Pr, San Diego

  22. [22]

    \ Song \ and\ author B

    Song, Z.D., Bernevig, B.A., 2022. Magic-angle twisted bilayer graphene as a topological heavy fermion problem. Phys. Rev. Lett. 129, 047601. URL: https://link.aps.org/doi/10.1103/PhysRevLett.129.047601, doi:10.1103/PhysRevLett.129.047601

  23. [23]

    Watson, A.B., Kong, T., MacDonald, A.H., Luskin, M.,

  24. [24]

    Journal of Mathematical Physics 64, 031502

    Bistritzer–MacDonald dynamics in twisted bilayer graphene. Journal of Mathematical Physics 64, 031502. URL: https://doi.org/10.1063/5.0115771, doi:10.1063/5.0115771, arXiv:https://pubs.aip.org/aip/jmp/article-pdf/doi/10.1063/5.0115771/16792366/031502_1_online.pdf

  25. [25]

    The kernel polynomial method

    Weiße, A., Wellein, G., Alvermann, A., Fehske, H., 2006. The kernel polynomial method. Reviews of Modern Physics 78, 275–

  26. [26]

    URL: https://link.aps.org/doi/10.1103/RevModPhys.78.275, doi:10.1103/RevModPhys.78.275. 52

  27. [27]

    Weak-field hall resistivity and spin-valley flavor symmetry breaking in magic-angle twisted bilayer graphene

    Xie, M., Macdonald, A.H., 2021. Weak-field hall resistivity and spin-valley flavor symmetry breaking in magic-angle twisted bilayer graphene. Physical Review Letters 127. doi:10.1103/PhysRevLett.127.196401

  28. [28]

    A high-order regularized delta-Chebyshev method for com- puting spectral densities

    Yi, J., Massatt, D., Horning, A., Luskin, M., Pixley, J.H., Kaye, J., 2025. A high-order regularized delta-Chebyshev method for com- puting spectral densities. URL: https://arxiv.org/abs/2512.03149, doi:10.48550/ARXIV.2512.03149

  29. [29]

    Twisted trilayer graphene, quasiperiodic superconductor

    Zhang, X., Zhu, Z., Wilson, J.H., Foster, M.S., 2025. Twisted trilayer graphene, quasiperiodic superconductor. URL: https://arxiv.org/abs/2512.22340,arXiv:2512.22340

  30. [30]

    Twisted trilayer graphene: A precisely tunable plat- form for correlated electrons

    Zhu, Z., Carr, S., Massatt, D., Luskin, M., Kaxiras, E., 2020a. Twisted trilayer graphene: A precisely tunable plat- form for correlated electrons. Phys. Rev. Lett. 125, 116404. URL: https://link.aps.org/doi/10.1103/PhysRevLett.125.116404, doi:10.1103/PhysRevLett.125.116404

  31. [31]

    Mod- eling mechanical relaxation in incommensurate trilayer van der waals heterostructures

    Zhu, Z., Cazeaux, P., Luskin, M., Kaxiras, E., 2020b. Mod- eling mechanical relaxation in incommensurate trilayer van der waals heterostructures. Phys. Rev. B 101, 224107. URL: https://link.aps.org/doi/10.1103/PhysRevB.101.224107, doi:10.1103/PhysRevB.101.224107. 53 (a) (b) (c) (d) Figure 6: Semilogy of the relative𝑊-error for a range of Chebyshev order𝑃:...