{"id":"3ca3a78f-e7a6-4b06-b0a5-455428e8fe8f","arxiv_id":"2411.11359","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"NanoARPES measurements trace rhombohedral graphene from 3 layers to bulk, showing gapped subbands and flat bands evolve into a Dirac nodal spiral semimetal with drumhead surface states.","lead":"Using high-resolution angle-resolved photoemission, this paper follows how rhombohedral multilayer graphene's electronic bands change from a few layers to the bulk: gapped bands and flat surface bands evolve into a three-dimensional Dirac nodal spiral with a drumhead surface state.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The nodal-spiral evidence rests on equating photoemission intensity maxima with Dirac-node positions, an assumption that is unvalidated and could bias the reported helicity.","rationale":"The reader's weakest-assumption analysis pinpoints the method used to locate Dirac nodes: the maximum of photoemission intensity in constant-energy maps. I agree that this is the most load-bearing experimental step. The paper's measured dispersions and flat bands are internally consistent and align with the SSH/SWM model and DFT slab calculations, which gives independent support for the qualitative band evolution. However, the specific assertion of a helicity-resolved Dirac nodal spiral depends on the assumption that the brightest FS spot tracks the node, and the manuscript provides no quantitative test of this assumption. The self-reported limitations, such as the nominal N=24 thickness and the post-submission note about a similar paper (Ref. 35), were also weighed. The N=24 issue weakens the intermediate point in the thickness evolution but does not undermine the bulk-limit claim. The Ref. 35 note is a provenance issue rather than a correctness issue and does not change the burden of proof. A direct re-evaluation of the raw maps with alternative node estimators and simulated intensities is the cleanest way to settle the concern. Since the reader already recommended CONDITIONAL and my concern does not move the verdict, I recommend no adjustment.","tokens_in":7993,"tokens_out":2557,"duration_ms":30201,"concrete_test":"Re-analyze the raw Fermi-surface maps in Fig. 4b,c with at least two independent node estimators: (i) fit the angular intensity profile around the dashed circle in each map after normalizing the maps; (ii) extract node positions from the band dispersion cuts at the same photon energies. In parallel, simulate the ARPES intensity from the DFT Bloch states with a slowly varying matrix element and kz lifetime broadening, and compare the intensity maximum with the true Dirac-node position. If the simulated or alternative-estimator node angles differ from the reported angles by more than the scatter in Fig. 4d, the spiral evidence needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of a topological Dirac nodal spiral semimetal in the N≈50 sample rests on Fig. 4b,c,d, where the Dirac-node angle is extracted by tracing the evolution of the strongest spectral intensity around the K/K' Fermi-surface contours. This procedure assumes that matrix-element effects, kz broadening, and the coexisting flat surface band do not shift the bright spot away from the true node. The paper provides no check of this assumption, for example by simulating ARPES intensities from the DFT bands or by comparing FS-maximum positions with dispersion minima. In addition, the conversion from photon energy to kz is not specified, so the plotted angle-versus-kz helical dependence is actually angle-versus-photon-energy unless an inner potential is assumed. If the intensity maximum is displaced by matrix elements or spectral-weight modulations, the clockwise/counterclockwise helicity and the kz periodicity would not be quantitatively established, although the coexistence of flat bands near EF at all photon energies would remain. This concern does not attack the DFT/slab prediction or the qualitative band evolution, but it limits the strength of the experimental proof of the spiral.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8120,"tokens_out":9030,"duration_ms":92940,"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":[{"comment":"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.","section":"Fig. 4b-d (Results)"},{"comment":"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.","section":"Fig. 4d (Results) and Methods/ARPES"}],"minor_comments":[{"comment":"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.","section":"Results, last paragraph"},{"comment":"The momentum resolution is quoted as '0.01 Å'; the unit should be Å^-1.","section":"Methods/ARPES"},{"comment":"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.","section":"Fig. 2 and Table S1"},{"comment":"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.","section":"Fig. 2d-e and Fig. S4"},{"comment":"The phrase '3D flat band in the bulk RMG' is confusing because the drumhead state is a surface (quasi-2D) state; please clarify.","section":"Discussion"},{"comment":"There are grammatical errors such as 'Dirac nodes that spirals'; these should be corrected.","section":"Abstract and throughout"},{"comment":"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.","section":"Note added in proof / Ref. 35"}],"recommendation":"major_revision","confidential_remarks":"The qualitative evolution of subbands and flat bands is well supported and likely solid. The one part I would not accept as is is the experimental extraction of the Dirac nodal spiral from Fermi-surface intensity maxima, together with the missing kz conversion. If the authors can provide an ARPES intensity simulation or a dispersion-minimum cross-check, the paper would be much stronger. I also suggest the editor ask the authors to integrate the comparison with Ref. 35 into the main text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the good news: this is the first systematic NanoARPES study of rhombohedral multilayer graphene from N=3 to N~48, and the qualitative evolution is compelling. The subband gap closes, the subband number grows, the flat band persists, and the photon-energy-independent flat band in the bulk sample is strong evidence for a surface (drumhead) state. The DFT slab calculations match the observed dispersions nicely, and the SSH/SWM framework provides a clean conceptual skeleton. The authors are also honest about the competing paper (Ref 35).\n\nThe soft spot is exactly where the reader's report puts it: the Dirac nodal spiral. The identification of the Dirac node with the maximum of FS intensity in Fig. 4 is not validated. Matrix-element effects, kz broadening, or the coexisting flat band could shift the bright spot. The paper does not show a simulation of ARPES intensity from the DFT bands, nor a comparison of FS-maximum positions to the dispersion minima. In addition, the photon-energy to kz conversion is not stated, so the plotted angle versus kz is actually angle versus photon energy unless an inner potential is assumed. That limits the quantitative claim about helicity and periodicity. These are fixable with a straightforward intensity simulation or a tight-binding forward model.\n\nMinor issues: the N=24 thickness is nominal with possible stacking faults; gaps and bandwidths have no error bars; the post-submission note deserves a short comparison to Ref 35 to clarify what is new.\n\nOverall, the central qualitative conclusion—thickness-driven evolution from 3D SSH gapped subbands to a nodal spiral semimetal with drumhead surface states—is well supported in its broad strokes, even if the spiral geometry is the least secure part. This is a serious experimental paper with real data and a clear theoretical comparison. It merits peer review; a good referee should ask for the intensity simulation and a kz calibration, but not for a rewrite from scratch.\n\nWho is it for? Groups working on rhombohedral graphene and on topological flat bands generally. I would bring it to a reading group and cite it in any context that needs the thickness-dependent ARPES baseline.","headline":"Solid systematic ARPES with one under-supported spiral claim that needs a forward-model check.","tokens_in":8801,"tokens_out":2345,"would_cite":false,"duration_ms":24202,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Rhombohedral multilayer graphene evolves from a gapped SSH-like stack into a topological Dirac nodal spiral semimetal as its layer count grows.","keywords":["rhombohedral multilayer graphene","topological Dirac nodal spiral semimetal","drumhead surface states","flat bands","NanoARPES","Su-Schrieffer-Heeger model","thickness-dependent topological transition","ABC stacking"],"falsifier":"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.","tokens_in":7743,"feed_emoji":"🌀","tokens_out":9000,"duration_ms":83889,"temperature":0.7,"pith_summary":"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.","feed_headline":"Stacking graphene turns flat bands into a nodal spiral semimetal","feed_subtitle":"Thicker rhombohedral graphene closes its band gap and winds Dirac nodes into spirals while surface flat bands survive.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Predicts that ABC-stacked graphene evolves from bulk subbands to a 3D Dirac cone structure in rhombohedral graphite, the specific transition this paper tests.","marker":"[18]"},{"why":"Provides the Slonczewski-Weiss-McClure band-structure model for ABC trilayer graphene used to describe the subbands and flat bands.","marker":"[9]"},{"why":"Establishes the N-pi Berry phase and topological character of ABC-stacked multilayer graphene, grounding the claim of nontrivial topology.","marker":"[8]"},{"why":"Introduces the Su-Schrieffer-Heeger chain whose 3D generalization explains the gapped subbands and topological surface flat bands.","marker":"[12]"},{"why":"Density-functional study of rhombohedral stacks identifying topological surface states and their bulk limit, used to support the drumhead surface states.","marker":"[27]"},{"why":"Shows interlayer interaction and kz dispersion in multilayer graphene, used to interpret the photon-energy dependence of subband spectral weight.","marker":"[28]"}],"fun_headline_variants":["Thicker rhombohedral graphene turns gap into a Dirac spiral","Flat bands persist as graphene's gap closes into a nodal spiral","Stacking rhombohedral graphene: gapped subbands become spiral nodes","Thickness-driven graphene transition: gap to Dirac nodal spiral semimetal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Thicker rhombohedral graphene turns gap into a Dirac spiral","Flat bands persist as graphene's gap closes into a nodal spiral","Stacking rhombohedral graphene: gapped subbands become spiral nodes","Thickness-driven graphene transition: gap to Dirac nodal spiral semimetal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000769,"raw_usage":{"total_tokens":3408,"prompt_tokens":945,"completion_tokens":2463,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":2388}},"tokens_in":561,"tokens_out":2463,"duration_ms":16947,"temperature":1.0,"reasoning_tokens":2388,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:36:32.946866+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Evolution and dimensional crossover from the bulk subbands in ABC -stacked graphene to a three -dimensional Dirac cone structure in rhombohedral graphite","cited_arxiv_id":null,"evidence_quote":"Predicts that ABC-stacked graphene evolves from bulk subbands to a 3D Dirac cone structure in rhombohedral graphite, the specific transition this paper tests."},{"cited_title":"Band structure of ABC -stacked graphene trilayers","cited_arxiv_id":null,"evidence_quote":"Provides the Slonczewski-Weiss-McClure band-structure model for ABC trilayer graphene used to describe the subbands and flat bands."},{"cited_title":"Trigonal warping and Berry’s phase N π in ABC-stacked multilayer graphene","cited_arxiv_id":null,"evidence_quote":"Establishes the N-pi Berry phase and topological character of ABC-stacked multilayer graphene, grounding the claim of nontrivial topology."},{"cited_title":"Solitons in Polyacetylene","cited_arxiv_id":null,"evidence_quote":"Introduces the Su-Schrieffer-Heeger chain whose 3D generalization explains the gapped subbands and topological surface flat bands."},{"cited_title":"Density functional investigation of rhombohedral stacks of graphene: Topological surface states, nonlinear dielectric response, and bulk limit","cited_arxiv_id":null,"evidence_quote":"Density-functional study of rhombohedral stacks identifying topological surface states and their bulk limit, used to support the drumhead surface states."},{"cited_title":"Phys Rev Lett 2007;98:206802","cited_arxiv_id":null,"evidence_quote":"Shows interlayer interaction and kz dispersion in multilayer graphene, used to interpret the photon-energy dependence of subband spectral weight."}],"review_version":1}