REVIEW 2 major objections 7 minor 35 references
Thickness-dependent Topological Phases and Flat Bands in Rhombohedral Multilayer Graphene
T0 review · 2 major / 7 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Rhombohedral multilayer graphene evolves from a gapped SSH-like stack into a topological Dirac nodal spiral semimetal as its layer count grows.
desk verdict Solid systematic ARPES with one under-supported spiral claim that needs a forward-model check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the 3D generalization of the 1D Su-Schrieffer-Heeger chain, realized by ABC stacking: the nearest-neighbor interlayer hopping forms a dimerized chain whose topological end states are the surface flat bands. The Slonczewski-Weiss-McClure model, including hoppings from gamma0 through gamma4, produces N flat-band crossing points and 2N-2 gapped bulk subbands; in the bulk limit the crossings line up along a helical spiral around K/K', forming the Dirac nodal spiral, and the flat band becomes the drumhead surface state. Experimentally, the load-bearing machinery is NanoARPES with photon-energy tuning: changing photon energy changes kz, so Fermi-surface maps at different photon energies reveal the angular rotation of the Dirac node, while the flat band's constant presence at all photon energies identifies the surface state.
What would settle it
A calculation of the ARPES spectral function for the reported nodal spiral that includes dipole matrix elements and kz broadening, compared pixel-by-pixel with the measured Fermi-surface maps, would settle whether the extracted bright-spot angles coincide with the true Dirac nodes; a mismatch larger than the fitting uncertainty would falsify the reported helicity.
Extended reading notes
Core claim
The paper reports NanoARPES measurements on N=3, N=24, and N≈48 rhombohedral graphene. In N=3, two gapped hole-type subbands sit below the surface flat band with a gap of about 276 meV; in N=24 the subband count increases and the gap shrinks to about 83 meV; in the bulk-like N=48 sample the subbands merge into the bulk continuum, the gap closes, and the flat band remains pinned near the Fermi level. Photon-energy-dependent maps around K and K' trace the Dirac node's angular position and show that it rotates with kz, with opposite helicity at the two valleys, which the paper identifies as the Dirac nodal spiral. DFT slab calculations for N=3, 24, and 48 reproduce the measured dispersions, and the flat band's persistence at all photon energies—absent from bulk DFT but present in slab calculations—marks it as the drumhead surface state.
Load-bearing premise
The claim's load-bearing assumption is that the brightest spot in each constant-energy photoemission map marks the true Dirac node; if matrix-element effects or kz broadening shift that maximum, the extracted spiral angle and helicity would not be reliable.
Editorial extensions
If this is right
- Bulk rhombohedral graphite is a topological Dirac nodal spiral semimetal, meaning its gapped subbands merge into gapless 3D Dirac cones that wind around K and K' with opposite helicity.
- The layer number N directly tunes the flat-band bandwidth and subband gap, so choosing N controls where correlation-driven instabilities such as superconductivity and fractional Chern states are most likely to appear.
- The surface flat band in the bulk limit is a drumhead surface state protected by bulk-boundary correspondence, so it should survive as a 2D conducting surface even when the bulk is a semimetal.
- As the flat band becomes increasingly three-dimensional with larger N, its Bloch-wave quantum geometry could support large superfluid stiffness in all three directions, extending the case for 3D flat-band superconductivity.
- Photon-energy-dependent ARPES can serve as a direct probe of Dirac nodal spiral semimetals, since changing photon energy moves kz and traces the node's angular position.
Reading between the lines
- A natural testable extension is to compute the ARPES spectral function for the nodal spiral including full dipole matrix elements and kz broadening; if the intensity maximum coincides with the true Dirac node, the angular tracking is robust, and if not, the reported helicity would need re-analysis.
- The same SSH-to-DNSS evolution should occur in other ABC-stacked honeycomb multilayers and could be tuned mechanically or by twist, suggesting a design rule for topological nodal-spiral semimetals beyond graphene.
- The measured bandwidth-versus-N curve provides input for Hubbard-style models, and a 'sweet spot' layer number with maximal flatness and minimal gap could be predicted and then tested by transport experiments on 3-, 4-, 5-, and thicker samples.
- Because the drumhead state is a surface state, scanning tunneling microscopy on a cleaved rhombohedral graphite surface could image its real-space shape and compare it with the k-space spiral predicted here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports spatially resolved ARPES measurements of rhombohedral multilayer graphene with nominally N = 3, 24, and ~48 layers, and uses them to claim a thickness-driven topological transition from a 3D SSH-like gapped phase to a Dirac nodal spiral semimetal. The authors compare measured band dispersions, subband gaps, and Fermi-surface maps to DFT slab calculations and to the Slonczewski-Weiss-McClure/SSH model, and interpret the photon-energy-independent flat band at EF in the bulk sample as the drumhead surface state. The central claim is that for N ~ 50 the gapped subbands become 3D Dirac cones that spiral around K/K' with opposite helicity, while the surface flat bands become drumhead states.
Significance. If established, the result would be an important systematic spectroscopic demonstration of the layer-number-controlled topological phase transition in rhombohedral graphene, a system of current interest for correlated and topological flat-band physics. The paper has several genuine strengths: the DFT and SSH/SWM comparisons are independent of the ARPES data with no fitted parameters; the surface nature of the flat band is supported by its photon-energy independence and by slab calculations; and the TEM/STEM characterization anchors the stacking. The main weakness is experimental: the Dirac nodal spiral is inferred from intensity maxima in Fermi-surface maps without validating that matrix-element and broadening effects do not shift the maxima from the true nodes, and the photon-energy-to-kz conversion is not specified. These issues are local to the spiral-extraction analysis but are load-bearing for the strongest version of the central claim.
major comments (2)
- [Fig. 4b-d (Results)] The identification of the Dirac-node position with the strongest spectral intensity on the dashed circle is load-bearing for the claimed nodal spiral, but it is not validated. Matrix-element modulations, kz broadening, and the coexisting surface flat band can all change the intensity distribution on a constant-energy map without moving the actual band crossing; the flat band near EF is particularly relevant because it sits in the same momentum region. Please provide a quantitative check, for example by simulating the ARPES intensity from the DFT bands with matrix elements and broadening, or by cross-correlating the intensity maxima with the minima of the dispersions in the photon-energy-dependent cuts of Fig. 3b. Until this is done, the clockwise/counterclockwise helicity and the spiral angle in Fig. 4d are not quantitatively established.
- [Fig. 4d (Results) and Methods/ARPES] The manuscript states that the extracted angles are plotted as a function of photon energy and kz, but the conversion from photon energy to kz is never specified. In particular, no inner potential or final-state dispersion is given, so the abscissa of Fig. 4d is only a photon-energy axis. If the claim is that the nodes spiral in momentum space with a particular kz periodicity or pitch, the authors must provide the kz conversion and its uncertainty; otherwise the linear fit in Fig. 4d demonstrates only a trend versus photon energy, which could be affected by matrix-element variations with photon energy.
minor comments (7)
- [Results, last paragraph] The sentence 'The angle positions ... plotted ... in Fig. 4c' should refer to Fig. 4d; the caption of Fig. 4d should also define the angle convention used for the linear fit.
- [Methods/ARPES] The momentum resolution is quoted as '0.01 Å'; the unit should be Å^-1.
- [Fig. 2 and Table S1] The subband gap is defined as the energy from EF to the top of the hole-type subbands; the calibration of EF should be stated explicitly.
- [Fig. 2d-e and Fig. S4] For the N=24 sample, the text notes that the actual number of layers may be smaller due to stacking faults; please quantify this uncertainty and discuss its effect on the comparison with the N=24 DFT calculation.
- [Discussion] The phrase '3D flat band in the bulk RMG' is confusing because the drumhead state is a surface (quasi-2D) state; please clarify.
- [Abstract and throughout] There are grammatical errors such as 'Dirac nodes that spirals'; these should be corrected.
- [Note added in proof / Ref. 35] The appended note states that a similar ARPES paper has been published; the main text should cite and explicitly compare with Ref. 35 so that the reader can assess the incremental contribution of this work.
Circularity Check
No significant circularity: experimental data and DFT/SSH theory are independent, and no fitted parameter is renamed as a prediction.
full rationale
The paper's central claim is that NanoARPES measurements of rhombohedral multilayer graphene show a thickness-driven evolution from gapped subbands and surface flat bands (thin layers) to a Dirac nodal spiral semimetal with drumhead surface states (bulk limit, N~50). The comparison chain is: (i) DFT slab and bulk calculations are performed with standard GGA functionals and fixed interlayer parameters; (ii) the SSH/SWM model is an independent analytical description with parameters from prior literature; (iii) the ARPES data are measured on exfoliated samples and are not used to set any model parameter. The paper reports subband gaps, bandwidths, and the angular positions of Fermi-surface intensity maxima, and compares these to DFT and model calculations. No equation in the paper defines a measured quantity in terms of the theory output, nor does any fitted constant enter the predicted dispersion. The only self-citation of note is Ref. [23], which describes the NanoARPES endstation; this is an instrumental description and does not support the central physics claim. The identification of Dirac-node positions with intensity maxima in Fig. 4 is a possible validity concern about matrix-element and kz-broadening effects, but that is an experimental systematic limitation, not a circularity in the derivation: the intensity maxima are not constructed from the theory being tested. The DFT calculations and the SSH/SWM prediction of the nodal spiral stand independently of the measured intensity distribution. Consequently, the central observation is self-contained against external benchmarks, and no load-bearing step reduces to its own input by definition or by self-citation.
Assumptions & free parameters
free parameters (1)
- Linear fit slope of Dirac node angle vs photon energy =
Not stated numerically; sign shown in Fig. 4d
assumptions (4)
- domain assumption Rhombohedral multilayer graphene's low-energy bands are described by the SSH/SWM model with hoppings gamma0-gamma4.
- domain assumption Changing the photon energy in ARPES accesses different out-of-plane momentum kz values.
- domain assumption The N=48 sample represents the bulk limit and the N=24 sample has approximately 24 ABC-stacked layers.
- domain assumption DFT with GGA and Grimme D2 vdW correction gives an accurate band structure for comparison with ARPES.
Cite this review
Pith. "Pith review of Thickness-dependent Topological Phases and Flat Bands in Rhombohedral Multilayer Graphene." pith.science (2026). https://pith.science/paper/5SI3KHOA
@misc{pith2026241111359,
author = {Pith},
title = {Pith review of: Thickness-dependent Topological Phases and Flat Bands in Rhombohedral Multilayer Graphene},
year = {2026},
howpublished = {\url{https://pith.science/paper/5SI3KHOA}},
note = {Machine review of arXiv:2411.11359}
}
read the original abstract
Rhombohedral multilayer graphene has emerged as an extraordinary platform for investigating exotic quantum states, such as superconductivity and fractional quantum anomalous Hall effects, mainly due to the existence of topological surface flatbands. Despite extensive research efforts, a systematic spectroscopic investigation on the evolution of its electronic structure from thin layers to bulk remains elusive. Using state-of-the-art angle-resolved photoemission spectroscopy with submicron spatial resolution, we directly probe and trace the thickness evolution of the topological electronic structures of rhombohedral multilayer graphene. As the layer number increases, the gapped subbands transform into the 3D Dirac nodes that spirals in the momentum space; while the flatbands are constantly observed around Fermi level, and eventually evolve into the topological drumhead surface states. This unique thickness-dependent topological phase transition can be well captured by the 3D generalization of 1D Su-Schrieffer-Heeger chain in thin layers, to the topological Dirac nodal spiral semimetal in the bulk limit. Our findings establish a solid foundation for exploring the exotic quantum phases with nontrivial topology and correlation effects in rhombohedral multilayer graphene.
Figures
Reference graph
Works this paper leans on
-
[1]
Superconductivity in rhombohedral trilayer graphene
Zhou H, Xie T, Taniguchi T, Watanabe K, Young AF. Superconductivity in rhombohedral trilayer graphene. Nature 2021;598:434–8
work page 2021
-
[2]
Spontaneous Quantum Hall States in Chirally Stacked Few-Layer Graphene Systems
Zhang F, Jung J, Fiete GA, Niu Q, MacDonald AH. Spontaneous Quantum Hall States in Chirally Stacked Few-Layer Graphene Systems. Phys Rev Lett 2011;106:156801
work page 2011
-
[3]
Spontaneous broken-symmetry insulator and metals in tetralayer rhombohedral graphene
Liu K, Zheng J, Sha Y , Lyu B, Li F, Park Y , et al. Spontaneous broken-symmetry insulator and metals in tetralayer rhombohedral graphene. Nat Nanotechnol 2024;19:188–95
work page 2024
-
[4]
Large quantum anomalous Hall effect in spin-orbit proximitized rhombohedral graphene
Han T, Lu Z, Yao Y , Yang J, Seo J, Yoon C, et al. Large quantum anomalous Hall effect in spin-orbit proximitized rhombohedral graphene. Science 2024;384:647–51
work page 2024
-
[5]
Fractional quantum anomalous hall effect in multilayer graphene
Lu Z, Han T, Yao Y , Reddy AP, Yang J, Seo J, et al. Fractional quantum anomalous hall effect in multilayer graphene. Nature 2024;626:759–64
work page 2024
-
[6]
Observation of a chern insulator in crystalline ABCA-tetralayer graphene with spin-orbit coupling
Sha Y , Zheng J, Liu K, Du H, Watanabe K, Taniguchi T, et al. Observation of a chern insulator in crystalline ABCA-tetralayer graphene with spin-orbit coupling. Science 2024;384:414–9
work page 2024
-
[7]
Correlated insulator and Chern insulators in pentalayer rhombohedral-stacked graphene
Han T, Lu Z, Scuri G, Sung J, Wang J, Han T, et al. Correlated insulator and Chern insulators in pentalayer rhombohedral-stacked graphene. Nat Nanotechnol 2024;19:181–7
work page 2024
-
[8]
Trigonal warping and Berry’s phase N π in ABC-stacked multilayer graphene
Koshino M, McCann E. Trigonal warping and Berry’s phase N π in ABC-stacked multilayer graphene. Phys Rev B 2009;80:165409
work page 2009
Show all 35 references
-
[9]
Band structure of ABC -stacked graphene trilayers
Zhang F, Sahu B, Min H, MacDonald AH. Band structure of ABC -stacked graphene trilayers. Phys Rev B 2010;82:35409
2010
-
[10]
Observation of competing, correlated ground states in the flat band of rhombohedral graphite
Hagymási I, Mohd Isa MS, Tajkov Z, Márity K, Oroszlány L, Koltai J, et al. Observation of competing, correlated ground states in the flat band of rhombohedral graphite. Science Advances 2022;8:eabo6879
2022
-
[11]
Moiréless correlations in ABCA graphene
Kerelsky A, Rubio-Verdú C, Xian L, Kennes DM, Halbertal D, Finney N, et al. Moiréless correlations in ABCA graphene. Proc Natl Acad Sci 2021;118:e2017366118
2021
-
[12]
Solitons in Polyacetylene
Su WP, Schrieffer JR, Heeger AJ. Solitons in Polyacetylene. Phys Rev Lett 1979;42:1698 – 701
1979
-
[13]
Electronic structure of multilayer graphene
Min H, MacDonald AH. Electronic structure of multilayer graphene. Prog Theor Phys 2008;176:227–52
2008
-
[14]
Flat electronic bands in long sequences of rhombohedral-stacked graphene
Henck H, Avila J, Ben Aziza Z, Pierucci D, Baima J, Pamuk B, et al. Flat electronic bands in long sequences of rhombohedral-stacked graphene. Phys Rev B 2018;97:245421
2018
-
[15]
Flat -Band Electronic Structure and Interlayer Spacing Influence in Rhombohedral Four-Layer Graphene
Wang W, Shi Y , Zakharov AA, Syväjärvi M, Yakimova R, Uhrberg RIG, et al. Flat -Band Electronic Structure and Interlayer Spacing Influence in Rhombohedral Four-Layer Graphene. Nano Lett 2018;18:5862–6
2018
-
[16]
Half- and quarter-metals in rhombohedral trilayer graphene
Zhou H, Xie T, Ghazaryan A, Holder T, Ehrets JR, Spanton EM, et al. Half- and quarter-metals in rhombohedral trilayer graphene. Nature 2021;598:429–33
2021
-
[17]
Stacking-Dependent Electronic Structure of Trilayer Graphene Resolved by Nanospot Angle -Resolved Photoemission Spectroscopy
Bao C, Yao W, Wang E, Chen C, Avila J, Asensio MC, et al. Stacking-Dependent Electronic Structure of Trilayer Graphene Resolved by Nanospot Angle -Resolved Photoemission Spectroscopy. Nano Lett 2017;17:1564–8
2017
-
[18]
Evolution and dimensional crossover from the bulk subbands in ABC -stacked graphene to a three -dimensional Dirac cone structure in rhombohedral graphite
Ho C-H, Chang C-P, Lin M-F. Evolution and dimensional crossover from the bulk subbands in ABC -stacked graphene to a three -dimensional Dirac cone structure in rhombohedral graphite. Phys Rev B 2016;93:075437
2016
-
[19]
Continuous epitaxy of single - crystal graphite films by isothermal carbon diffusion through nickel
Zhang Z, Ding M, Cheng T, Qiao R, Zhao M, Luo M, et al. Continuous epitaxy of single - crystal graphite films by isothermal carbon diffusion through nickel. Nat Nanotechnol 2022;17:1258–64. 16 / 16
2022
-
[20]
Generalized Gradient Approximation Made Simple
Perdew JP, Burke K, Ernzerhof M. Generalized Gradient Approximation Made Simple. Phys Rev Lett 1996;77:3865–8
1996
-
[21]
Ab initio molecular dynamics for liquid metals
Kresse G, Hafner J. Ab initio molecular dynamics for liquid metals. Phys Rev B 1993;47:558– 61
1993
-
[22]
Projector augmented-wave method
Blöchl PE. Projector augmented-wave method. Phys Rev B 1994;50:17953–79
1994
-
[23]
Nano-ARPES Endstation at BL07U of Shanghai Synchrotron Radiation Facility
Gao H, Xiao H, Wang F, Zhu F, Wang M, Liu Z, et al. Nano-ARPES Endstation at BL07U of Shanghai Synchrotron Radiation Facility. Synchrotron Radiat News 2024;37:12–7
2024
-
[24]
A Short Course on Topological Insulators
Asbóth JK, Oroszlány L, Pályi A. A Short Course on Topological Insulators. vol. 919. Cham: Springer International Publishing; 2016
2016
-
[25]
Slonczewski JC, Weiss PR, Band Structure of Graphite. Phys. Rev. 1958;109:272
1958
-
[26]
McClure W, Band structure of graphite and the de -Haas-van Alphen effect . Phys. Rev. 1957;08:612
1957
-
[27]
Density functional investigation of rhombohedral stacks of graphene: Topological surface states, nonlinear dielectric response, and bulk limit
Xiao R, Tasnádi F, Koepernik K, V enderbos JWF, Richter M, Taut M. Density functional investigation of rhombohedral stacks of graphene: Topological surface states, nonlinear dielectric response, and bulk limit. Phys Rev B 2011;84:165404
2011
-
[28]
Phys Rev Lett 2007;98:206802
Ohta T, Bostwick A, McChesney JL, Seyller T, Horn K, Rotenberg E, Interlayer Interaction and Electronic Screening in Multilayer Graphene Investigated with Angle -Resolved Photoemission Spectroscopy. Phys Rev Lett 2007;98:206802
2007
-
[29]
Journal of Electron Spectroscopy and Related Phenomena 2003;130:65-78
Strocov VN, Intrinsic accuracy in 3 -dimensional photoemission band mapping . Journal of Electron Spectroscopy and Related Phenomena 2003;130:65-78
2003
-
[30]
Angle-resolved photoemission studies of quantum materials
Sobota JA, He Y , Shen Z-X. Angle-resolved photoemission studies of quantum materials. Rev Mod Phys 2021;93:025006
2021
-
[31]
Rhombohedral graphene goes correlated at four or five layers
Cao Y . Rhombohedral graphene goes correlated at four or five layers. Nat Nanotechnol 2024;19:139–40
2024
-
[32]
Electronic phase separation in multilayer rhombohedral graphite
Shi Y , Xu S, Yang Y , Slizovskiy S, Morozov SV , Son S-K, et al. Electronic phase separation in multilayer rhombohedral graphite. Nature 2020;584:210–4
2020
-
[33]
Layer-dependent evolution of electronic structures and correlations in rhombohedral multilayer graphene arXiv:2312.13637
Zhou Y-Y , Zhang Y , Zhang S, Cai H, Tong L-H, Tian Y , et al. Layer-dependent evolution of electronic structures and correlations in rhombohedral multilayer graphene arXiv:2312.13637
-
[34]
Designing Three-Dimensional Flat Bands in Nodal-Line Semimetals
Lau A, Hyart T, Autieri C, Chen A, Pikulin DI. Designing Three-Dimensional Flat Bands in Nodal-Line Semimetals. Phys Rev X 2021;11:031017
2021
-
[35]
Correlated topological flat bands in rhombohedral graphite, PNAS 2024; 121: e2410714121
Zhang HY , Li Q, Scheer MG, Wang RQ, Tuo CY , Zou NL, et al. Correlated topological flat bands in rhombohedral graphite, PNAS 2024; 121: e2410714121
2024
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.