Pith. sign in

REVIEW 3 major objections 6 minor 5 references

Drumhead Surface States of Rhombohedral Graphite with Near Ideal Quantum Geometry Condition

T0 review · 3 major / 6 minor · reviewed 2026-07-31 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Thick rhombohedral graphite’s drumhead surface states curve by tens of meV and meet the ideal quantum-geometry condition only on their inner rim.

desk verdict Solid first-principles baseline: DSS curvature matches ARPES and IQG is rim=1 / center≈0.96 for thick RG; ordinary single-particle limits, not a load-bearing flaw. read the letter →

arxiv 2607.28491 v1 pith:NEXKLJ6Q submitted 2026-07-30 cond-mat.mtrl-sci

Lin-Lin Wang This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords rhombohedralgraphitedrumheadsurfacestatesidealquantumgeometryBerrycurvaturemetricweaktopologicalinsulatorflatbandsARPES
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 shows that the nearly flat drumhead surface states on thick and semi-infinite rhombohedral graphite are not perfectly flat: they have a convex curvature of roughly 34–38 meV that matches recent ARPES data. Using density-functional theory and Wannier functions it further maps the quantum-geometry ratio of Berry curvature to quantum metric across these states. That ratio reaches a minimum of about 0.96 at the K-point center yet equals exactly 1 on the inner rim of the drumhead region. Bulk rhombohedral graphite is identified as a weak topological insulator once spin-orbit coupling is included. The single-particle results are offered as a realistic starting point for later many-body theories of superconductivity or fractional Chern states that rely on ideal quantum geometry.

What carries the argument

The ratio |Ω|/TrG evaluated from Kubo formulas on DFT-derived Wannier Hamiltonians of finite and semi-infinite slabs; this ratio is the concrete diagnostic of how close the pristine drumhead states sit to the ideal quantum-geometry condition.

What would settle it

A higher-resolution ARPES map or a quantum-geometry measurement that shows either a flat (zero-curvature) drumhead or a |Ω|/TrG profile that is uniformly 1 (or uniformly far from 1) across the entire drumhead region would directly contradict the central claim.

Watch

Extended reading notes

Core claim

First-principles calculations establish that the drumhead surface states of thick (L = 60) and semi-infinite rhombohedral-graphite slabs possess a sizable convex curvature whose depth agrees with ARPES, and that the ideal-quantum-geometry ratio |Ω|/TrG is only approximately 1 at the K-point center while strictly equaling 1 on the inner rim of the drumhead region.

Load-bearing premise

The single-particle DFT plus Wannier band structure with fixed experimental lattice constants is assumed accurate enough that many-body theories can safely adopt its curvature and quantum-geometry profile as their non-interacting starting point.

Editorial extensions

If this is right

  • Many-body models of surface superconductivity in thick rhombohedral graphite should start from a convex, ~35 meV deep drumhead rather than a strictly flat band.
  • The ideal-quantum-geometry condition is realized on an annular rim, not at the K-point center, so annular Fermi-surface physics is favored.
  • Bulk rhombohedral graphite must be treated as a weak topological insulator once spin-orbit coupling is retained.
  • Finite-slab calculations converge to the same convex |Ω|/TrG profile already by L ≈ 18–60 layers, giving a practical thickness target for theory and experiment.

Reading between the lines

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

  • Displacement-field or moiré tuning that moves the chemical potential onto the inner rim could lock the system into the strict IQG condition more effectively than doping the K-point center.
  • The residual 4 % deviation from ideal geometry at the center may set a quantitative upper bound on the superfluid stiffness or FCI gap that pure flat-band models overestimate.
  • Side-surface Dirac cones protected by the weak-TI index offer an independent transport signature that could be checked on cleaved edges.
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

3 major / 6 minor

Summary. The manuscript presents a DFT+Wannier study of bulk rhombohedral graphite (RG) and its drumhead surface states (DSS). Without SOC, bulk RG is treated as a chiral semimetal with a spiral nodal line that projects to a trigonal nodal loop; with SOC it is identified as a weak TI (Wilson loops on kz=0 and 0.5 planes, Z2w3=1, Z4=0). For semi-infinite surfaces and thick slabs (up to L=60), the DSS is shown to be convex rather than strictly flat, with a depth of ~34–38 meV that matches recent ARPES (~32 meV). Quantum geometry is evaluated via Kubo formulas on the slabs: |Ω|/TrG exhibits a convex profile that reaches a minimum ≈0.96 at the K-point center while equaling 1 on the inner rim of the DSS region. Results are cross-checked with PBE+D3 and r2SCAN+rVV10 at fixed experimental lattice constants. The authors position the work as a single-particle baseline for future many-body studies of superconductivity and fractional Chern physics in thick RG.

Significance. If the reported DSS curvature and rim-versus-center IQG profile hold, the paper supplies a concrete, first-principles single-particle reference that many-body theories of surface superconductivity and related correlated states in thick RG can adopt in place of idealized nearest-neighbor tight-binding models. The dual XC-functional check, L-convergence of the |Ω|/TrG convex shape, and quantitative agreement with ARPES curvature are genuine strengths and make the baseline more credible than a pure model calculation. The weak-TI classification with SOC is a useful clarification of the bulk SPT context, even if secondary to the quantum-geometry results. The work is incremental rather than transformative, but it fills a documented gap between TB proposals (Ref. 49) and experiment (Ref. 52) at the single-particle level.

major comments (3)
  1. [§II Eqs. (1)–(3); §III Fig. 3] §II Eqs. (1)–(3) and §III Fig. 3: The Kubo sums for Ω and TrG are written over valence (n) and conduction (m) bands of the finite slab. For the DSS the valence–conduction gap is tiny (especially with normal SOC) and the surface bands hybridize with bulk continuum at the rim. The manuscript does not state which bands are retained, whether a surface-projection or energy window is applied, or how near-degeneracies and residual hybridization at L=60 are handled. Because the central claim is the quantitative rim-versus-center profile (|Ω|/TrG → 1 on the inner rim, ≈0.96 at K), this choice is load-bearing and should be specified (and, if possible, shown to be stable under reasonable window variations).
  2. [§III; Fig. 2(e–f)] §III and Fig. 2: Surface spectral functions are computed with SOC enhanced by a factor of 100 “to increase the surface band gap but with no change of band order or connectivity.” While the topology statement is unaffected, the DSS depth quoted for the semi-infinite case (38 meV) and the visual comparison to ARPES are taken from these enhanced-SOC plots. A brief check that the curvature depth (and, for thick slabs, the |Ω|/TrG profile) is unchanged at physical SOC would remove any residual concern that the enhancement distorts the single-particle baseline being advertised.
  3. [Abstract; §III] Abstract vs. §III: The abstract states that the K-point minimum is “approximately |Ω|/TrG≈1,” while the body and Fig. 3 report a converged value of 0.960 (PBE) / 0.958 (r2SCAN). The distinction between “near-ideal at center” and “strictly ideal on the inner rim” is the paper’s main quantitative message; the abstract should quote the actual minimum (≈0.96) so that the claim is not overstated relative to the figures.
minor comments (6)
  1. [Abstract; §I; §III] Throughout: “angular resolved phono-emission” / “AREPS” should be “photoemission” / “ARPES.”
  2. [§I] §I: “2-dimesional” → “two-dimensional.”
  3. [Fig. 1; §III] Fig. 1 caption and text: “K1,” “KM,” “ΓN,” etc. are typeset inconsistently; a uniform notation for bulk vs. surface high-symmetry points would help.
  4. [§II] §II: The Wannier energy window (EF±3 eV) and C sp projection are stated, but a short note on the spread or the fidelity of the Wannier interpolation near the nodal line would strengthen reproducibility.
  5. [Fig. 3] Fig. 3(b,d): The horizontal axis units and the precise cut direction ([110] vs. Γ–K–M) should be labeled on the panels themselves, not only in the caption.
  6. [Data Availability] Data Availability statement is present but empty of a repository link or DOI; a concrete deposition (structures, Wannier Hamiltonians, or QGT grids) would match the claim that the data are “openly available.”

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: DSS curvature and |Ω|/TrG profiles are computed DFT+Wannier outputs, not forced by construction or self-citation

full rationale

The paper’s load-bearing chain is standard single-particle computation: experimental lattice constants → DFT (PBE and r2SCAN+rVV10) → C-sp Wannier model → Wilson loops / surface Green’s functions for topology and DSS, and Kubo formulas (Eqs. 1–3) on finite slabs (L up to 60) for |Ω| and TrG. The reported convex DSS depth (~34–38 meV) and convex |Ω|/TrG profile (center ≈0.96, inner rim =1) are numerical outputs of that pipeline, not quantities fitted to ARPES or to the IQG=1 target; the paper explicitly reports a deviation from ideal at K rather than forcing equality. Benchmarks (ARPES Ref. 52; TB IQG proposal Ref. 49) are external. The sole self-citation (Ref. 5) is unrelated surface-energy work and is not load-bearing. No equation equates input to claimed prediction by construction, and no uniqueness theorem or ansatz is imported from overlapping authors to forbid alternatives. Residual XC/finite-slab limitations are ordinary genre caveats, not circularity.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The central claims rest on standard DFT/Wannier practice plus a few domain choices (fixed experimental lattice, C sp projection, Kubo QGT on valence/conduction manifolds, finite slabs as semi-infinite proxies). No new physical entities are postulated. Free parameters are limited to conventional computational settings and one visualization SOC scale factor that the authors state does not alter topology.

free parameters (3)
  • SOC enhancement factor (×100) for surface spectral functions = 100
    Used only to enlarge the surface gap for visualization in Fig. 2; authors assert band order and topological index are unchanged. Still a hand-chosen scale.
  • Wannier energy window and orbital set = EF±3 eV, C sp
    C sp orbitals fitted to reproduce bulk bands in EF±3 eV; standard but choice-dependent for surface and QGT accuracy.
  • Slab thickness L used as semi-infinite proxy = L=60
    IQG convergence claimed by comparing L=9,18,60; L=60 taken as representative of thick/semi-infinite limit without continuum extrapolation.
assumptions (5)
  • domain assumption Kohn-Sham DFT with PBE+D3 or r2SCAN+rVV10 adequately describes the single-particle bands and velocity matrix elements of RG near EF.
    Invoked throughout §II–III; meta-GGA check addresses gap underestimation but not many-body renormalization of DSS.
  • domain assumption Experimental bulk lattice constants (a=2.461 Å, c=10.061 Å) without ionic relaxation are the correct geometry for surface and slab calculations.
    Stated in §II; interlayer distance enters DSS curvature and hopping range.
  • standard math Quantum geometry tensor of the DSS is given by the Kubo formulas (Eqs. 1–3) summing over DFT valence and conduction bands from the Wannier Hamiltonian.
    Standard single-particle QGT definition; used for all |Ω| and TrG plots in §III.
  • standard math Wilson-loop odd winding on kz=0 and kz=0.5 planes implies weak TI indices Z2w3=1, Z4=0 for bulk RG with SOC.
    §III and Fig. 2(a–b); standard Fu-Kane/Wilson-loop bulk-boundary correspondence.
  • domain assumption Longer-range hoppings automatically included in Wannier models (beyond NN/NNN) are responsible for DSS curvature relative to idealized TB flatness.
    Stated in §III citing Ref. 49; underpins why first-principles DSS is not strictly flat.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Drumhead Surface States of Rhombohedral Graphite with Near Ideal Quantum Geometry Condition." pith.science (2026). https://pith.science/paper/NEXKLJ6Q

@misc{pith2026260728491,
  author       = {Pith},
  title        = {Pith review of: Drumhead Surface States of Rhombohedral Graphite with Near Ideal Quantum Geometry Condition},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NEXKLJ6Q}},
  note         = {Machine review of arXiv:2607.28491}
}
abstract

The ideal quantum geometry (IQG) condition of equal magnitude between quantum metric ($G(k)$) and Berry curvature ($\Omega(k)$) as in $|\Omega|$/Tr$G$=1 has been associated with realizing fractional Chern insulators for the flat band in few layers of rhombohedral graphene (RG). More recently the IQG condition has also been proposed for superconductivity in the flat band for thick RG layers. Using density functional theory and Wannier functions, we study the symmetry-protected topology of bulk RG, drumhead surface state (DSS) of semi-infinite RG surface, and the IQG condition for thick RG slabs. We find that bulk RG is a weak topological insulator with spin-orbit coupling (SOC), besides being effectively a chiral semimetal with chiral nodal line without SOC. We also find that the DSS flat band of semi-infinite and thick RG slabs have a sizable convex curvature with depth agreeing with the recent angular resolved phono-emission spectroscopy experiment. The calculated IQG also shows a convex shape with the K point at the center having a minimum with approximately $|\Omega|$/Tr$G$$\approx$1. But the inner rim of the DSS region shows the strict IQG condition of $|\Omega|$/Tr$G$=1. These results on semi-infinite and thick RG slabs from first-principles calculations provide useful information for the pristine DSS at the single particle level for future studies to consider when many-body interactions and strong correlations will be included.

Figures

Figures reproduced from arXiv: 2607.28491 by the authors.

Figure 1
Figure 1. Band structure of rhombohedral graphite (RG) without spin-orbit coupling (SOC). (a) The crystal structure and primitive unit cell of RG with ABC stacking of honeycomb layers. (b) The 3D Brillouin zone (BZ) and (001) surface 2D BZ of RG with the high￾symmetry k-points related via projection. (c) Band structure of RG along high-symmetry directions. The highest valence (lowest conduction) band is in blue (red). The K1 … view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

5 extracted references · 1 linked inside Pith

  1. [14]

    G. Chen, A. L. Sharpe, P. Gallagher, I. T. Rosen, E. J. Fox, L. Jiang, B. Lyu, H. Li, K. Watanabe, T. Taniguchi, J. Jung, Z. Shi, D. Goldhaber-Gordon, Y. Zhang, F. Wang, Signatures of tunable superconductivity in a trilayer graphene moiré superlattice. Nature 572, 215-219 (2019). 15. H. Zhou, T. Xie, T. Taniguchi, K. Watanabe, A. F. Young, Superconductivi...

    2019

  2. [27]

    C. L. Kane, E. J. Mele, Quantum spin Hall effect in graphene. Phys Rev Lett 95, 226801 (2005). 28. F. Tang, H. C. Po, A. Vishwanath, X. Wan, Topological materials discovery by large-order symmetry indicators. Science Advances 5, eaau8725 (2019). 29. J. P. Provost, G. Vallee, Riemannian structure on manifolds of quantum states. Communications in Mathematic...

    arXiv doi:10.48550/arxiv.2508.00469 2005

  3. [44]

    Z. Dong, A. S. Patri, T. Senthil, Theory of Quantum Anomalous Hall Phases in Pentalayer Rhombohedral Graphene Moir\'e Structures. Phys Rev Lett 133, 206502 (2024). 45. J. Dong, T. Wang, T. Wang, T. Soejima, M. P. Zaletel, A. Vishwanath, D. E. Parker, Anomalous Hall Crystals in Rhombohedral Multilayer Graphene. I. Interaction-Driven Chern Bands and Fractio...

    2024

  4. [56]

    Huang, T

    C. Huang, T. M. R. Wolf, W. Qin, N. Wei, I. V. Blinov, A. H. MacDonald, Spin and orbital metallic magnetism in rhombohedral trilayer graphene. Phys Rev B 107, L121405 (2023). 57. W. Qin, C. Huang, T. Wolf, N. Wei, I. Blinov, A. H. MacDonald, Functional Renormalization Group Study of Superconductivity in Rhombohedral Trilayer Graphene. Phys Rev Lett 130, 1...

    2023

  5. [72]

    Q. Wu, S. Zhang, H.-F. Song, M. Troyer, A. A. Soluyanov, WannierTools: An open-source software package for novel topological materials. Computer Physics Communications 224, 405-416 (2018). 73. A. Bosak, M. Krisch, M. Mohr, J. Maultzsch, C. Thomsen, Elasticity of single-crystalline graphite: Inelastic x-ray scattering study. Phys Rev B 75, 153408 (2007). 7...

    2018

Pith tools

Reviewed July 31, 2026 · model on record in the stance chip above.