Pith. sign in

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 →

arxiv 2411.11359 v2 pith:5SI3KHOA submitted 2024-11-18 cond-mat.mes-hall cond-mat.str-el

classification cond-mat.mes-hallcond-mat.str-el
keywords rhombohedralmultilayergraphenetopologicalDiracnodalspiralsemimetaldrumheadsurfacestatesflatbandsNanoARPESSu-Schrieffer-Heegermodelthickness-dependenttransitionABCstacking
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper uses spatially resolved angle-resolved photoemission to watch the electronic structure of rhombohedral multilayer graphene as it thickens from about 3 layers to about 50. It claims that the gapped subbands, which resemble a three-dimensional Su-Schrieffer-Heeger chain in thin layers, progressively close their gap and merge into gapless 3D Dirac cones that wind around the K and K' points in opposite directions, while the topological surface flat band stays near the Fermi level and becomes a drumhead surface state. That would make bulk rhombohedral graphite a topological Dirac nodal spiral semimetal, connecting the few-layer SSH picture to a 3D nodal semimetal in the bulk limit. A sympathetic reader would care because the layer number is the dial that controls both band flatness and topology, which is directly relevant to the correlated quantum phases observed in this material.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 7 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [Methods/ARPES] The momentum resolution is quoted as '0.01 Å'; the unit should be Å^-1.
  3. [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.
  4. [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.
  5. [Discussion] The phrase '3D flat band in the bulk RMG' is confusing because the drumhead state is a surface (quasi-2D) state; please clarify.
  6. [Abstract and throughout] There are grammatical errors such as 'Dirac nodes that spirals'; these should be corrected.
  7. [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

0 steps flagged · score 0.0 of 10

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 1 free parameters · 4 assumptions · 0 invented entities

The central claim rests on standard domain assumptions of ARPES kz mapping, the SSH/SWM model, DFT accuracy, and sample layer-count assignment. No new physical entities are introduced. The only fitted quantity is the slope of the Dirac node angle vs photon energy, which is a characterization of the observed rotation, not an input to a derivation.

free parameters (1)
  • Linear fit slope of Dirac node angle vs photon energy = Not stated numerically; sign shown in Fig. 4d
    Fitted to the extracted angular positions of the photoemission intensity maximum. The sign of the slope is used as evidence for the helicity of the Dirac nodal spiral, but the exact slope value is not used in further derivations.
assumptions (4)
  • domain assumption Rhombohedral multilayer graphene's low-energy bands are described by the SSH/SWM model with hoppings gamma0-gamma4.
    Used throughout to interpret subbands and surface flat bands. Hopping parameters are inherited from prior literature, not fitted in this paper.
  • domain assumption Changing the photon energy in ARPES accesses different out-of-plane momentum kz values.
    Standard ARPES assumption. The paper acknowledges kz broadening (Ref 29) but still uses photon-energy cuts to infer the Dirac node spiral in Fig. 4.
  • domain assumption The N=48 sample represents the bulk limit and the N=24 sample has approximately 24 ABC-stacked layers.
    The N=24 layer count is explicitly stated as nominal with possible stacking faults (Fig. S4), and N=48 is treated as the bulk limit for DFT simulation. This enters in the N=24 and N=48 results sections.
  • domain assumption DFT with GGA and Grimme D2 vdW correction gives an accurate band structure for comparison with ARPES.
    Standard computational assumption used for all band structure comparisons; the accuracy is not independently benchmarked in the paper.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2411.11359 by the authors.

Figure 2
Figure 2. From [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

35 extracted references · 35 canonical work pages

  1. [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

  2. [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

  3. [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

  4. [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

  5. [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

  6. [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

  7. [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

  8. [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

Show all 35 references
  1. [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

  2. [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

  3. [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

  4. [12]

    Solitons in Polyacetylene

    Su WP, Schrieffer JR, Heeger AJ. Solitons in Polyacetylene. Phys Rev Lett 1979;42:1698 – 701

  5. [13]

    Electronic structure of multilayer graphene

    Min H, MacDonald AH. Electronic structure of multilayer graphene. Prog Theor Phys 2008;176:227–52

  6. [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

  7. [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

  8. [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

  9. [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

  10. [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

  11. [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

  12. [20]

    Generalized Gradient Approximation Made Simple

    Perdew JP, Burke K, Ernzerhof M. Generalized Gradient Approximation Made Simple. Phys Rev Lett 1996;77:3865–8

  13. [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

  14. [22]

    Projector augmented-wave method

    Blöchl PE. Projector augmented-wave method. Phys Rev B 1994;50:17953–79

  15. [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

  16. [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

  17. [25]

    Slonczewski JC, Weiss PR, Band Structure of Graphite. Phys. Rev. 1958;109:272

  18. [26]

    McClure W, Band structure of graphite and the de -Haas-van Alphen effect . Phys. Rev. 1957;08:612

  19. [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

  20. [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

  21. [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

  22. [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

  23. [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

  24. [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

  25. [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

  26. [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

  27. [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

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.