Electron localization, charge redistribution, and emergence of topological states at graphite junctions
Pith reviewed 2026-05-10 15:30 UTC · model grok-4.3
The pith
Junctions between Bernal and rhombohedral graphite half-crystals produce localized electronic states and flat bands.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
At junctions between half-crystals of graphite with Bernal (AB) or rhombohedral (ABC) stacking, junction-localized electronic states emerge ubiquitously. All systems but one involving a rhombohedral half-crystal support a flat-band expected to exhibit electronic instabilities and strongly-correlated states. Nascent flat-band states associated with finite rhombohedral stacking sequences extend the physics into pure Bernal systems.
What carries the argument
Junction-localized electronic states and flat bands arising from charge redistribution at stacking boundaries between Bernal and rhombohedral graphite half-crystals, computed with charge self-consistent tight-binding and embedding potentials.
If this is right
- Flat bands in most graphite junction systems are expected to exhibit electronic instabilities and strongly-correlated states.
- Topologically non-trivial states arise due to interrupted rhombohedral stacking localized at the junction edges.
- The flat-band physics extends into pure Bernal systems through finite rhombohedral stacking sequences.
- Localized states appear as a ubiquitous feature at all studied graphite stacking junctions.
Where Pith is reading between the lines
- Similar localization and flat-band effects may appear at stacking boundaries in other layered van der Waals materials.
- Local spectroscopy at artificially created stacking junctions could directly image the predicted states.
- Charge redistribution at these junctions may alter transport or optical response in graphite heterostructures.
- The approach suggests stacking-order engineering as an alternative route to flat-band systems without moiré patterns.
Load-bearing premise
The charge self-consistent tight-binding method combined with embedding potentials sufficiently captures the low-energy electronic structure and charge redistribution at the junctions without significant artifacts from the model approximations.
What would settle it
Angle-resolved photoemission or scanning tunneling spectroscopy that either detects or fails to detect flat bands at pure Bernal graphite junctions or at rhombohedral junctions.
Figures
read the original abstract
Low-energy electronic behavior in graphite crystals is highly dependent on the relative stacking arrangement of the constituent layers. Topologically non-trivial electronic states can arise due to interrupted rhombohedral (ABC) stacking, localized at the edges of the stacking region, but not in the case of Bernal (AB) stacking. Here, we study the electronic properties of junctions between half-crystals of graphite of either Bernal or rhombohedral stacking, using a charge self-consistent tight-binding method and embedding potentials to account for the influence of layers far from the junction. We find junction-localized electronic states to be a ubiquitous feature, and all systems but one involving a rhombohedral half-crystal support a flat-band expected to exhibit electronic instabilities and strongly-correlated states. Nascent flat-band states associated with finite rhombohedral stacking sequences extend the physics into pure Bernal systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates electronic properties at junctions between Bernal (AB) and rhombohedral (ABC) graphite half-crystals via charge self-consistent tight-binding calculations with embedding potentials. It claims junction-localized states are ubiquitous, that flat bands (expected to drive instabilities and correlated states) appear in all but one rhombohedral half-crystal case, and that nascent flat bands from finite rhombohedral sequences extend the physics into pure Bernal systems.
Significance. If the results are robust, the findings would extend known topological and flat-band physics from stacking edges to internal junctions in graphite, with potential relevance for correlated electron behavior in layered carbon systems. The embedding-potential approach for treating extended systems is a methodological strength.
major comments (1)
- [Methods] Methods section: The central claims rest on the low-energy spectrum from the charge self-consistent tight-binding Hamiltonian with embedding potentials. No benchmarks against full DFT (or GW) calculations for the exact junction geometries are provided, and no convergence tests of the embedding decay length versus flat-band width are reported. This leaves open whether the flat bands are physical or numerical artifacts of the model approximations.
minor comments (1)
- [Abstract] Abstract: The statement 'all systems but one' is imprecise without indicating the total number of configurations examined or identifying the exceptional case; adding this would aid readability.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for identifying a key point regarding validation of the computational approach. We address the major comment below.
read point-by-point responses
-
Referee: [Methods] Methods section: The central claims rest on the low-energy spectrum from the charge self-consistent tight-binding Hamiltonian with embedding potentials. No benchmarks against full DFT (or GW) calculations for the exact junction geometries are provided, and no convergence tests of the embedding decay length versus flat-band width are reported. This leaves open whether the flat bands are physical or numerical artifacts of the model approximations.
Authors: The tight-binding model parameters are taken from established DFT-based fits that accurately capture the low-energy physics of graphite in both Bernal and rhombohedral stackings. Full self-consistent DFT or GW calculations on the junction geometries with embedding potentials for semi-infinite leads are computationally prohibitive for the system sizes required. However, we agree that explicit convergence tests are needed. We have now performed additional calculations varying the embedding decay length and find that the flat-band features and their widths converge for decay lengths beyond approximately 10 layers; the reported results use lengths well into this converged regime. We will add a dedicated subsection and supplementary figure documenting these tests in the revised manuscript. revision: partial
- Direct benchmarks against full DFT or GW calculations for the exact junction geometries, which remain beyond current computational capabilities for the embedded semi-infinite systems studied.
Circularity Check
No circularity: direct computational exploration via established TB method
full rationale
The paper performs numerical calculations of electronic structure at graphite junctions using a charge-self-consistent tight-binding Hamiltonian with embedding potentials. Results (junction-localized states and flat bands) are outputs of the simulation for specific geometries, not derived by re-expressing fitted parameters or self-cited premises as predictions. No self-definitional loops, fitted-input predictions, or load-bearing self-citation chains appear in the method or results chain. The approach is self-contained against external benchmarks (prior TB parametrizations) and does not reduce the central claims to tautology.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Charge self-consistent tight-binding with embedding potentials accurately reproduces low-energy electronic states and charge redistribution at graphite junctions.
Reference graph
Works this paper leans on
-
[1]
J. W. McClure, Band structure of graphite and de Haas- van Alphen effect, Phys. Rev.108, 612 (1957)
work page 1957
-
[2]
J. C. Slonczewski and P. R. Weiss, Band structure of graphite, Phys. Rev.109, 272 (1958)
work page 1958
-
[3]
J. C. Charlier, X. Gonze, and J. P. Michenaud, First- principles study of the electronic properties of graphite, Phys. Rev. B43, 4579 (1991)
work page 1991
-
[4]
H. S. Lipson and A. R. Stokes, The structure of graphite, Proc. R. Soc. Lond. A181, 101 (1942)
work page 1942
-
[5]
R. R. Haering, Band structure of rhombohedral graphite, Can. J. Phys.36, 352 (1958)
work page 1958
-
[6]
C. H. Lui, Z. Li, Z. Chen, P. V. Klimov, L. E. Brus, and T. F. Heinz, Imaging stacking order in few-layer graphene, Nano Lett.11, 164 (2011)
work page 2011
-
[7]
F. Laves and Y. Baskin, On the formation of the rhom- bohedral graphite modification, Zeitschrift f¨ ur Kristallo- graphie107, 337 (1956)
work page 1956
- [8]
-
[9]
J. P. Nery, M. Calandra, and F. Mauri, Ab-initio en- ergetics of graphite and multilayer graphene: stability of Bernal versus rhombohedral stacking, 2D Mater.8, 035006 (2021)
work page 2021
-
[10]
K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Electric field effect in atomically thin carbon films, Science306, 666 (2004)
work page 2004
-
[11]
K. F. Mak, M. Y. Sfeir, J. A. Misewich, and T. F. Heinz, The evolution of electronic structure in few-layer graphene revealed by optical spectroscopy, Proc. Natl. Acad. Sci. U.S.A107, 14999 (2010)
work page 2010
-
[12]
J. Yin, S. Slizovskiy, Y. Cao, S. Hu, Y. Yang, I. Lobanova, B. A. Piot, S. K. Son, S. Ozdemir, T. Taniguchi, K. Watanabe, K. S. Novoselov, F. Guinea, A. K. Geim, V. Fal’ko, and A. Mishchenko, Dimensional reduction, quantum Hall effect and layer parity in graphite films, Nature Phys.437, 14999 (2019)
work page 2019
-
[13]
Y. Shi, S. Xu, Y. Yang, S. Slizovskiy, S. V. Morozov, S.- K. Son, S. Ozdemir, C. Mullan, J. Barrier, J. Yin, A. I. Berdyugin, B. A. Piot, T. Taniguchi, K. Watanabe, V. I. Fal’ko, K. S. Novoselov, A. K. Geim, and A. Mishchenko, Electronic phase separation in multilayer rhombohedral graphite, Nature584, 210 (2020)
work page 2020
-
[14]
Y. Yang, Y.-C. Zou, C. R. Woods, Y. Shi, J. Yin, S. Xu, S. Ozdemir, T. Taniguchi, K. Watanabe, A. K. Geim, K. S. Novoselov, S. J. Haigh, and A. Mishchenko, Stack- ing order in graphite films controlled by van der Waals technology, Nano Lett.19, 8526 (2019)
work page 2019
-
[15]
F. Pan, K. Ni, Y. Ma, H. Wu, X. Tang, J. Xiong, Y. Yang, C. Ye, H. Yuan, M.-L. Lin, J. Dai, M. Zhu, P.-H. Tan, Y. Zhu, and K. S. Novoselov, Phase-changing in graphite assisted by interface charge injection, Nano Lett.21, 5648 (2021)
work page 2021
-
[16]
H. Wu, X. Yu, M. Zhu, Z. Zhu, J. Zhang, S. Zhang, S. Qin, G. Wang, G. Peng, J. Dai, and K. S. Novoselov, Direct visualization and manipulation of stacking or- ders in few-layer graphene by dynamic atomic force mi- croscopy, J. Phys. Chem. Lett.12, 7328 (2021)
work page 2021
-
[17]
Y. Yeo, Y. Sharaby, N. Roy, N. Raab, K. Watanabe, T. Taniguchi, and M. B. Shalom, Polytype switching by super-lubricant van der Waals cavity arrays, Nature638, 389 (2025)
work page 2025
-
[18]
C. Hu, T. Wu, X. Huang, Y. Dong, J. Chen, Z. Zhang, B. Lyu, S. Ma, K. Watanabe, T. Taniguchi, G. Xie, X. Li, Q. Liang, and Z. Shi, In-situ twistable bilayer graphene, Sci. Rep.12, 204 (2022)
work page 2022
- [19]
-
[20]
D. P. Arovas and F. Guinea, Stacking faults, bound states, and quantum Hall plateaus in crystalline graphite, Phys. Rev. B78, 245416 (2008)
work page 2008
-
[21]
M. Taut, K. Koepernik, and M. Richter, Electronic struc- ture of stacking faults in hexagonal graphite, Phys. Rev. B88, 205411 (2013)
work page 2013
-
[22]
M. Koshino and E. McCann, Multilayer graphenes with mixed stacking structure: Interplay of Bernal and rhom- bohedral stacking, Phys. Rev. B87, 045420 (2013)
work page 2013
-
[23]
M. Taut, K. Koepernik, and M. Richter, Electronic struc- ture of stacking faults in rhombohedral graphite, Phys. Rev. B90, 085312 (2014)
work page 2014
-
[24]
M. Taut and K. Koepernik, Electronic structure of in- terfaces between hexagonal and rhombohedral graphite, Phys. Rev. B94, 035446 (2016)
work page 2016
-
[25]
A. Garc´ ıa-Ruiz, S. Slizovskiy, M. Mucha-Kruczy´ nski, and V. I. Fal’ko, Spectroscopic signatures of electronic excita- tions in Raman scattering in thin films of rhombohedral graphite, Nano Lett.19, 6152 (2019)
work page 2019
-
[26]
J. H. Muten, A. J. Copeland, and E. McCann, Exchange interaction, disorder, and stacking faults in rhombohe- dral graphene multilayers, Phys. Rev. B104, 035404 (2021)
work page 2021
-
[27]
A. Garcia-Ruiz, S. Slizovskiy, and V. I. Fal’ko, Flat bands for electrons in rhombohedral graphene multilayers with a twin boundary, Adv. Mater. Interfaces10, 2202221 (2023)
work page 2023
-
[28]
P. J. Sarsfield, S. Slizovskiy, M. Koshino, and V. Fal’ko, Electronic properties of stacking faults in Bernal graphite, npj Comput. Mater.11, 142 (2025)
work page 2025
- [29]
-
[30]
E. Cisternas and J. D. Correa, Theoretical reproduction of superstructures revealed by STM on bilayer graphene, Chem. Phys.409, 74 (2012). 12
work page 2012
-
[31]
A. Dey, S. A. Chowdhury, T. Pe˜ na, S. Singh, S. M. Wu, and H. Askari, An Atomistic Insight into Moir´ e Recon- struction in Twisted Bilayer Graphene beyond the Magic Angle, ACS Appl. Eng. Mater.1, 970 (2023)
work page 2023
-
[32]
M. S. Dresselhaus and G. Dresselhaus, Intercalation com- pounds of graphite, Adv. Phys.30, 139 (1981)
work page 1981
-
[33]
B. Partoens and F. M. Peeters, From graphene to graphite: Electronic structure around theKpoint, Phys. Rev. B74, 075404 (2006)
work page 2006
-
[34]
D. D. L. Chung, Review Graphite, J. Mater. Sci.37, 1475 (2002)
work page 2002
-
[35]
S. Slizovskiy, E. McCann, M. Koshino, and V. I. Fal’ko, Films of rhombohedral graphite as two-dimensional topo- logical semimetals, Commun. Phys.2, 164 (2019)
work page 2019
-
[36]
T. Cea, P. A. Pantale´ on, V. T. Phong, and F. Guinea, Superconductivity from repulsive interactions in rhombo- hedral trilayer graphene: A Kohn-Luttinger-like mecha- nism, Phys. Rev. B105, 075432 (2022)
work page 2022
-
[37]
V. Kaladzhyan, S. Pinon, F. Joucken, Z. Ge, E. A. Quezada-Lopez, T. Taniguchi, K. Watanabe, J. Velasco, Jr., and C. Bena, Surface states and quasiparticle inter- ference in Bernal and rhombohedral graphite with and without trigonal warping, Phys. Rev. B104, 155418 (2021)
work page 2021
-
[38]
E. N. Economou,Green ’s Functions in Quantum Physics (Springer-Verlag Berlin Heidelberg, 1979)
work page 1979
-
[39]
M. Lannoo and P. Friedel,Atomic and Electronic Struc- ture of Surfaces(Springer-Verlag Berlin Heidelberg, 1991)
work page 1991
-
[40]
See Supplemental Material at [insert URL] for an analytic derivation of the surface Green’s functions in the mini- mal model, full numerical results from the self-consistent charge redistribution procedure at all 12 junctions, full numerical results for the local density of states on layers J−1,J,J+ 1, andJ+ 2 at each junction, analytical expressions for ...
-
[41]
N. V. Tepliakov, Q. Wu, and O. V. Yazyev, Crystal field effect and electric field screening in multilayer graphene with and without twist, Nano Lett.21, 4636 (2021)
work page 2021
-
[42]
S. Slizovskiy, A. Garcia-Ruiz, A. I. Berdyugin, N. Xin, T. Taniguchi, K. Watanabe, A. K. Geim, N. D. Drum- mond, and V. I. Fal’ko, Out-of-plane dielectric suscepti- bility of graphene in twistronic and Bernal bilayers, Nano Lett.21, 6678 (2021)
work page 2021
-
[43]
Y. Wu, A. Ilie, and S. Crampin, Self-consistent charge and dipole density functional tight binding method and application to carbon-based systems, Comput. Mater. Sci.134, 206 (2017)
work page 2017
-
[44]
J. C. Slater, Atomic shielding constants, Phys. Rev.36, 57 (1930)
work page 1930
-
[45]
W. P. Su, J. R. Schrieffer, and A. J. Heeger, Solitons in polyacetylene, Phys. Rev. Lett.42, 1698 (1979)
work page 1979
-
[46]
J. K. Asb´ oth, L. Oroszl´ any, and A. P´ alyi, The Su- Schrieffer-Heeger (SSH) Model, inA Short Course on Topological Insulators: Band Structure and Edge States in One and Two Dimensions(Springer International Publishing, 2016) pp. 1–22
work page 2016
- [47]
-
[48]
J. L. Ma˜ nes, F. Guinea, and M. A. H. Vozmediano, Ex- istence and topological stability of Fermi points in mul- tilayered graphene, Phys. Rev. B75, 155424 (2007)
work page 2007
-
[49]
M. Koshino and E. McCann, Trigonal warping and Berry’s phaseN πin ABC-stacked multilayer graphene, Phys. Rev. B80, 165409 (2009)
work page 2009
-
[50]
D. Pierucci, H. Sediri, M. Hajlaoui, J.-C. Girard, T. Brumme, M. Calandra, E. Velez-Fort, G. Patri- arche, M. G. Silly, G. Ferro, V. Souli` ere, M. Marangolo, F. Sirotti, F. Mauri, and A. Ouerghi, Evidence for flat bands near the Fermi level in epitaxial rhombohedral multilayer graphene, ACS Nano9, 5432 (2015)
work page 2015
-
[51]
I. Hagym´ asi, M. S. M. Isa, Z. Tajkov, K. M´ arity, L. Oroszl´ any, J. Koltai, A. Alassaf, P. Kun, K. Kandrai, A. P´ alink´ as, P. Vancs´ o, L. Tapaszt´ o, and P. Nemes-Incze, Observation of competing, correlated ground states in the flat band of rhombohedral graphite, Sci. Adv.8, eabo6879 (2022)
work page 2022
- [52]
- [53]
-
[54]
D.-H. Xu, J. Yuan2, Z.-J. Yao, Y. Zhou, J.-H. Gao, and F.-C. Zhang, Stacking order, interaction, and weak sur- face magnetism in layered graphene sheets, Phys. Rev. B 86, 201404 (2012)
work page 2012
-
[55]
W. A. Mu˜ noz, L. Covaci, and F. M. Peeters, Tight- binding description of intrinsic superconducting correla- tions in multilayer graphene, Phys. Rev. B87, 134509 (2013)
work page 2013
-
[56]
N. B. Kopnin, M. Ij¨ as, A. Harju, and T. T. Heikkil¨ a, High-temperature surface superconductivity in rhombo- hedral graphite, Phys. Rev. B87, 140503 (2013)
work page 2013
-
[57]
H. Henck, J. Avila, Z. B. Aziza, D. Pierucci, J. Baima, B. Pamuk, J. Chaste, D. Utt, M. Bartos, K. Nogajew- ski, B. A. Piot, M. Orlita, M. Potemski, M. Calandra, M. C. Asensio, F. Mauri, C. Faugeras, and A. Ouerghi, Flat electronic bands in long sequences of rhombohedral- stacked graphene, Phys. Rev. B97, 245421 (2018)
work page 2018
-
[58]
O. A. Awoga, T. L¨ othman, and A. M. Black-Schaffer, Superconductivity and magnetism in the surface states of ABC-stacked multilayer graphene, Phys. Rev. B108, 144504 (2023)
work page 2023
-
[59]
C. Coletti, S. Forti, A. Principi, K. V. Emtsev, A. A. Zakharov, K. M. Daniels, B. K. Daas, M. V. S. Chan- drashekhar, T. Ouisse, D. Chaussende, A. H. MacDon- ald, M. Polini, and U. Starke, Revealing the electronic band structure of trilayer graphene on SiC: An angle- resolved photoemission study, Phys. Rev. B88, 155439 (2013)
work page 2013
-
[60]
C. Bao, W. Yao, E. Wang, C. Chen, J. Avila, M. C.Asensio, and S. Zhou, Stacking-Dependent Elec- tronic Structure of Trilayer Graphene Resolved by Nanospot Angle-Resolved Photoemission Spectroscopy, Nano Lett.17, 1564 (2017)
work page 2017
-
[61]
W. Wang, Y. Shi, A. A. Zakharov, M. Syv¨ aj¨ arvi, R. Yaki- mova, R. I. G. Uhrberg, and J. Sun, Flat-Band Electronic Structure and Interlayer Spacing Influence in Rhombohe- dral Four-Layer Graphene, Nano Lett.18, 5862 (2018)
work page 2018
-
[62]
H. Zhou, T. Xie, T. Taniguchi, K. Watanabe, and A. F. Young, Superconductivity in rhombohedral trilayer graphene, Nature598, 434 (2021). 13
work page 2021
-
[63]
P. A. Pantale´ on, A. Jimeno-Pozo, H. Sainz-Cruz, V. T. Phong, T. Cea, and F. Guinea, Superconductivity and correlated phases in non-twisted bilayer and trilayer graphene, Nat. Rev. Phys.5, 304 (2023)
work page 2023
-
[64]
E. V. Bostr¨ om, A. Fischer, J. B. Profe, J. Zhang, D. M. Kennes, and A. Rubio, Phonon-mediated unconventional superconductivity in rhombohedral stacked multilayer graphene, npj Comput. Mater.10, 163 (2024)
work page 2024
-
[65]
T. Han, Z. Lu, Z. Hadjri, L. Shi, Z. Wu, W. Xu, Y. Yao, A. A. Cotten, O. S. Sedeh, H. Weldeyesus, J. Yang, J. Seo, S. Ye, M. Zhou, H. Liu, G. Shi, Z. Hua, K. Watan- abe, T. Taniguchi, P. Xiong, D. M. Zumb¨ uhl, L. Fu, and L. Ju, Signatures of chiral superconductivity in rhombo- hedral graphene, Nature643, 654 (2025). Supplemental Material: Electron locali...
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.