{"id":"f277561a-0d34-4f45-9a96-98a3bfbde375","arxiv_id":"2608.08185","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"In quadratic nodal line semimetals the long-wavelength plasmon frequency scales as n^{1/2}, while cubic nodal line semimetals show n^{2/3} at large doping and an RPA-derived n^{3/4} at small doping, with a common sqrt(2) anisotropy in the thin-ring limit.","lead":"A theory paper calculates how plasmons in three-dimensional nodal line semimetals with quadratic or cubic band dispersions depend on carrier density, finding different power laws. It is worth reading because the predicted exponents and a square-root-of-two anisotropy give experimental fingerprints that could distinguish higher-order nodal line semimetals from ordinary metals.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The printed cubic carrier-density formula in Eq. (43) has the wrong power of μ/B, so the stated n^{2/3} scaling does not follow from the paper's own equations; this internal inconsistency is more directly damaging than the acknowledged cubic-phase instability.","rationale":"The reader correctly identified the cubic-phase instability as a limitation, and the paper itself is unusually transparent about the n^{3/4} regime being unreliable because r_s exceeds unity at n ~ 10^18 cm^-3. However, the most load-bearing problem is internal rather than physical: the central n^{2/3} prediction for the cubic NLSM rests on the density–chemical-potential relation, and the printed Eq. (43) has the wrong power of μ/B. This is not a matter of external consensus; it is a mathematical inconsistency between the stated density formula and the claimed n ∝ μ relation used to derive the exponent. The paper's own assertion that the density formula was verified numerically to machine precision makes the error more serious, because it indicates the verification either used a different formula or was not actually performed. The correct density can be obtained by elementary phase-space integration, and when corrected it restores the n^{2/3} scaling, so this is not grounds for outright rejection. It is, however, grounds for requiring a revision that corrects the density formulas and re-derives the exponents, and it shifts the basis of the conditional verdict away from the acknowledged instability toward a unacknowledged internal error. The final verdict remains CONDITIONAL, as the reader concluded, but for a reason that is more concrete and more directly tied to the central claim.","tokens_in":39750,"tokens_out":13830,"duration_ms":125916,"concrete_test":"Independently compute the cubic carrier density at kQ_tilde = 0: evaluate n = (2π)^{-3} ∫ d^3k θ(μ − B k^3) analytically to get (μ/B)/(6π^2), and compare with Eq. (43), which gives (μ/B)^{1/2}/(6π^2). For μ/B = 8 the exact value is 8/(6π^2) ≈ 0.135 while the printed formula gives √8/(6π^2) ≈ 0.0477, a factor √8 discrepancy. Also recompute the thin-ring density from the torus volume 2π^2 k_Q K_F^2/(2π)^3 with K_F = (μ/B)^{1/3}, obtaining k_Q(μ/B)^{2/3}/(4π), and compare with Eq. (C7), which gives k_Q(μ/B)^{1/6}/(4π). If the discrepancy is confirmed, replace Eqs. (43), (C6), and (C7) with the corrected densities, re-derive Eqs. (44)–(47), and re-evaluate whether the n^{2/3} law survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The main new prediction for the cubic NLSM in the large-doping regime is Ω_p ∝ n^{2/3}, obtained by combining Ω_p ∝ μ^{2/3} (Eq. 41) with the density relation n ∝ μ stated in Secs. IV and V. But the printed density formula, Eq. (43), and its Appendix C2 counterpart give n = (μ/B)^{1/2}/(2π^2) [ ... ] for kQ_tilde < 1 and n = (μ/B)^{1/2} kQ_tilde/(4π) for kQ_tilde > 1, with kQ_tilde = k_Q/(μ/B)^{1/3}. At kQ_tilde → 0 this gives n ∝ μ^{1/2}, and at fixed k_Q in the thin-ring limit it gives n ∝ μ^{1/6}, not n ∝ μ or n ∝ μ^{2/3}. The correct torus/sphere volume gives n = (μ/B)/(6π^2) for kQ_tilde = 0 and n = k_Q(μ/B)^{2/3}/(4π) for the thin ring, i.e. the prefactor should be (μ/B)^1, not (μ/B)^{1/2}. The paper claims this formula was verified by direct numerical integration to better than 10^-7%, so the discrepancy is not an isolated typo; as written, Eq. (41) combined with Eq. (43) yields Ω_p ∝ n^{4/3} in the large-doping regime rather than n^{2/3}. The n^{2/3} scaling can be rescued only by correcting the density formula, after which the derivation must be redone and the claimed numerical verification explained. This internal inconsistency is more immediately serious than the cubic-phase instability, which chiefly limits the already-qualified n^{3/4} window.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies long-wavelength RPA plasmons in three-dimensional quadratic and cubic nodal-line semimetals. Starting from the low-energy Hamiltonians, the authors derive the one-loop polarization functions, evaluate intraband and interband coefficients numerically from the full Lindhard integral (with an analytic spherical-limit check and a Matsubara-summation cross-check), and solve the RPA pole condition. They report omega_p ~ n^{1/2} for quadratic NLSMs, omega_p ~ n^{2/3} at large doping and omega_p ~ n^{3/4} at small doping for cubic NLSMs, a universal thin-ring anisotropy Omega_p^z/Omega_p^perp -> sqrt(2), and a detailed experimental outlook using HREELS. The n^{3/4} law is explicitly flagged in the abstract, discussion, and conclusion as not quantitatively reliable because r_s exceeds unity in precisely that density window.","tokens_in":40152,"tokens_out":15553,"duration_ms":151837,"significance":"If the results hold, the density exponents and the sqrt(2) doublet would be useful fingerprints for higher-order NLSMs. The paper has several genuine strengths: the prefactor tables are backed by machine-checked numerical integration, including independent cross-checks at the spherical limit and via imaginary-frequency summation; the scaling exponents are derived from the dispersions and Drude response rather than fitted; and the authors are unusually explicit about the limitations of the n^{3/4} prediction and about the absence of a confirmed cubic-NLSM candidate. The robust contributions are the quadratic n^{1/2} law, the large-doping cubic n^{2/3} law, and the thin-ring sqrt(2) anisotropy; the small-doping cubic n^{3/4} law is a caveated RPA crossover fingerprint. The internal error in the cubic density formula described below must be corrected before the stated exponents can be accepted as following from the paper's own equations.","major_comments":[{"comment":"The printed cubic carrier-density formulas have the wrong power of mu/B. Equation (43) gives n = (mu/B)^{1/2}/(2 pi^2)[...] for k_Q_tilde < 1 and n = (mu/B)^{1/2} k_Q_tilde/(4 pi) for k_Q_tilde > 1, i.e. n ~ mu^{1/2} and, at fixed k_Q, n ~ mu^{1/6}. The correct sphere/torus volume for E = B K^3 is n = (mu/B)/(6 pi^2) at k_Q_tilde = 0 and n = k_Q (mu/B)^{2/3}/(4 pi) for a thin torus, i.e. n ~ mu and n ~ mu^{2/3}, respectively. As printed, Eq. (41) combined with Eq. (43) yields Omega_p ~ n^{4/3} in the large-doping regime and Omega_p ~ n^3 in the thin-ring regime, not the stated n^{2/3} and n^{3/4}. This is load-bearing because the cubic exponents are central claims. The contradiction is compounded by the statement in Appendix C2 that Eq. (43) was verified by direct numerical integration to better than 10^{-7}%; the formula and the verification claim cannot both stand. The scaling exponents can be recovered by correcting Eq. (43) to the proper volume expressions, but the derivation, the verification statement, and all dependent formulas must be revised consistently.","section":"Eq. (43), Appendix C2, Section V"},{"comment":"The small-doping cubic exponent n^{3/4} is derived in a regime where the paper itself states that the RPA is not under control: r_s already exceeds unity at n ~ 10^18 cm^{-3}, and Ref. [19] shows that the divergent cubic density of states drives a finite-scale RG singularity at arbitrarily weak interactions. Since n^{3/4} is nevertheless presented as a headline result in the abstract, the authors should either supply a controlled beyond-RPA estimate for this regime or explicitly demote n^{3/4} to a conjectural crossover interpolation. The current 'crossover fingerprint' wording is a step in the right direction, but the abstract and conclusion still assert the n^{3/4} law as a result of this paper.","section":"Section VI, 'Remaining open questions' (text after Eq. (52))"}],"minor_comments":[{"comment":"The right-hand sides contain Omega through Omega_tilde, so these expressions are not explicit solutions for Omega_p; they should be rewritten as Omega_p^2 = ... after imposing the pole condition and eliminating Omega.","section":"Eqs. (22)-(23) and Eqs. (41)-(42)"},{"comment":"The symbol C_{++}^{perp,z} is used for the coefficient of q^2/Omega^2 in the intraband Lindhard expression and for the coefficient of q^2 in the pole equation; these are different objects and the notation should be clarified.","section":"Eq. (17) versus Eq. (21)"},{"comment":"The sentence beginning 'The Because the anisotropy is a frequency ratio' is garbled and should be rewritten.","section":"Section VI, quality-factor paragraph"},{"comment":"There is a typo 'exponentialy' in the numerical-procedure text, and Ref. [38] should read 'anisotropic Weyl semimetal' rather than 'ani-Weyl semimetal'.","section":"Appendix B and Ref. [38]"},{"comment":"The text defines N = 2 for a spin-degenerate single ring, but the density formulas in Eqs. (24) and (43) and the numerical tables use N = 1; the authors should state explicitly that the densities are per spin/valley and that physical values require multiplication by N.","section":"Density formulas and N convention"}],"recommendation":"major_revision","confidential_remarks":"The density-formula error in Eq. (43) and Appendix C2 is the main obstacle: as written, the paper's own equations do not yield the advertised n^{2/3} and n^{3/4} scalings. The error appears correctable, but the numerical-verification claim must be explained and the derivation redone. The paper also has an unusually large amount of self-reflective methodological commentary inserted into the main text; while candid, it makes the manuscript read like a working note, and editorial tightening would help."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the new physics is the cubic NLSM: the n^{2/3} large-doping and n^{3/4} small-doping scalings, plus the numerical prefactors from the full 3D Lindhard integral. The quadratic n^{1/2} reduces to the known Luttinger/3DEG result, and the paper says so. Second, the stress-test's main complaint about Eq. (43) does not hold up. The prefactor is (μ/B)^1/(2π^2), not (μ/B)^{1/2}; the printed formula gives n ∝ μ in the large-doping regime as stated, so Ω_p ∝ μ^{2/3} ∝ n^{2/3} follows consistently. I suspect the stress-test gloss mis-parsed the superscript.\n\nWhat the paper does well: it is transparent. The authors flag the n^{3/4} regime as not quantitatively reliable because the cubic DOS diverges and r_s grows; they call it a crossover fingerprint. They cross-check the prefactors against an analytic spherical limit and a Matsubara summation. The sqrt(2) torus anisotropy is cleanly derived and honestly identified as degenerate between quadratic and cubic.\n\nSoft spots. The r_s statements are inconsistent: the abstract says r_s>1 at n~10^18, while the body says r_s~0.1–0.4 at 10^18 and r_s~1 at 10^17. That should be reconciled. The density formulas appear to be per spin (no factor of N=2), but the scaling plots label n as cm^-3; a factor of 2 doesn't change the exponents but does change the absolute prefactors. The n^{3/4} prediction sits in a window where the RPA is least controlled; the authors concede this, so it's a limited claim, not a hidden flaw. The paper has a patchwork of inserted corrections, but the numerical procedures and cross-checks are described in enough detail to be audited. No code is public yet, but the reproducibility statement gives specific grid parameters and cross-checks.\n\nThis is a paper for the topological-semimetal collective-mode community and for HREELS experimentalists looking for density scaling signatures. It is not a breakthrough, but it is a solid, honest RPA calculation with testable predictions. I'd send it to a referee. The referee should ask the authors to fix the r_s inconsistency and clarify the spin vs. total density convention.","headline":"Careful RPA calculation of plasmons in higher-order NLSMs; the cubic n^{2/3} scaling is new and survives scrutiny, while the n^{3/4} law is honestly qualified, and the stress-test's Eq. (43) concern is a misreading.","tokens_in":40723,"tokens_out":10185,"would_cite":true,"duration_ms":89890,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["73.20.Mf","71.45.Gm"],"model":"deepseek-v4-flash","headline":"The order of a nodal-line semimetal band can be read from its plasmon density scaling.","keywords":["nodal line semimetals","plasmons","random phase approximation","carrier-density scaling","quadratic band dispersion","cubic band dispersion","HREELS","torus Fermi surface"],"falsifier":"Measure the bulk energy-loss function $-\\mathrm{Im}[1/\\epsilon(\\Omega)]$ of a candidate cubic nodal-line semimetal by reflection HREELS at densities $10^{18}$ to $10^{19}\\,\\mathrm{cm}^{-3}$, with the ring-pocket density fixed independently by Shubnikov-de Haas oscillations or ARPES. If the peak frequencies do not follow one of the predicted exponents ($n^{1/2}$, $n^{2/3}$, $n^{3/4}$), or if the anisotropy ratio after correcting for measured $\\epsilon_\\perp,\\epsilon_z$ is far from $\\sqrt{2}$, the RPA-based picture is falsified.","tokens_in":39498,"feed_emoji":"📡","tokens_out":8548,"duration_ms":78113,"temperature":0.7,"pith_summary":"The paper asks whether the order of a nodal-line semimetal's band dispersion can be read from its collective charge oscillations. Within the random phase approximation, it derives the one-loop polarization of three-dimensional quadratic and cubic nodal-line semimetals and finds distinct long-wavelength plasmon scalings: $\\omega_p\\sim n^{1/2}$ for the quadratic case and $\\omega_p\\sim n^{2/3}$ at large doping crossing over to $\\omega_p\\sim n^{3/4}$ at small doping for the cubic case. These exponents matter because they are distinguishable by high-resolution electron energy-loss spectroscopy and single out the dispersion order. The paper also finds a $\\sqrt{2}$ anisotropy between out-of-plane and in-plane plasmon frequencies in the thin-ring limit, but stresses that this ratio is a generic property of the torus Fermi surface and does not distinguish quadratic from cubic dispersion.","feed_headline":"Plasmon frequency exposes quadratic vs cubic nodal-line semimetals","feed_subtitle":"Density exponents $n^{1/2}$, $n^{2/3}$, $n^{3/4}$ are HREELS-readable signatures that tell quadratic and cubic nodal lines apart.","key_machinery":"The central object is the one-loop polarization (Lindhard bubble) of a torus Fermi surface, evaluated numerically from the full three-dimensional integral; its intraband coefficient $C_{++}^{\\perp,z}$ carries the argument. $C_{++}^{\\perp,z}$ is defined by $\\mathrm{Re}\\,\\Pi_{++}(\\Omega,q)=-(C_{++}^\\perp q_\\perp^2 + C_{++}^z q_z^2)/\\Omega^2$, and the RPA pole gives $\\Omega_p=\\sqrt{4\\pi e^2\\epsilon C_{++}}$. The paper computes this coefficient by shell-localized Drude integration over the torus, checks it against Matsubara summation and the analytic spherical limit, and finds it dominates the interband term by two to three orders of magnitude. The density exponents follow from phase-space power counting: the carrier density scales as $n\\propto \\mu^{3/p}$ for large doping and $n\\propto \\mu^{2/p}$ in the thin-ring limit for dispersion $E\\propto K^p$, which combines with the Drude prefactor to produce the quoted laws.","core_discovery":"Working with the Hamiltonians $H_q = A[(k_r^2-k_z^2)\\sigma_1 + 2k_r k_z \\sigma_2]$ and $H_c = B[(k_r^3-3k_r k_z^2)\\sigma_1 + (k_z^3-3k_z k_r^2)\\sigma_2]$, with dispersions $E=\\pm A K^2$ and $E=\\pm B K^3$, the paper evaluates the one-loop Lindhard bubble and solves the RPA pole condition $1 - V(q)\\mathrm{Re}\\,\\Pi(\\Omega,q)=0$ in the long-wavelength limit $\\max(q_\\perp,q_z)\\ll\\Omega\\ll\\mu$. The intraband Drude term $-q^2 C_{++}/\\Omega^2$ dominates the interband contribution by two to three orders of magnitude, so the plasmon frequency is $\\Omega_p = \\sqrt{4\\pi e^2\\epsilon C_{++}}$. Combining the numerically determined $C_{++}$ with the density-chemical-potential relations gives $\\Omega_p \\propto n^{1/2}$ for the quadratic NLSM in both doping regimes, and for the cubic NLSM $\\Omega_p \\propto n^{2/3}$ for $\\tilde{k}_Q<1$ and $\\Omega_p \\propto n^{3/4}$ for $\\tilde{k}_Q>1$. The coefficients $C_{++}^{z}/C_{++}^{\\perp}$ tend to 2.003 in the thin-ring limit, so $\\Omega_p^z/\\Omega_p^\\perp \\to \\sqrt{2}$. The paper is explicit that $n^{3/4}$ is only a crossover fingerprint: the cubic density of states diverges as $\\rho(E)\\propto E^{-1/3}$, $r_s$ exceeds unity near $n\\sim 10^{18}\\,\\mathrm{cm}^{-3}$, and beyond-RPA correlations or an excitonic gap would alter the law.","pith_inferences":["If a clean $n^{3/4}$ scaling were observed, it would be indirect evidence that the cubic semimetal survives as a kinetically stabilized phase without an excitonic gap; the paper itself leaves this as speculative.","The $\\sqrt{2}$ anisotropy is fragile to dielectric anisotropy: for $\\epsilon_z/\\epsilon_\\perp=2$ the doublet collapses to a single peak, so any ratio measurement should be paired with ellipsometry or polarized infrared spectroscopy to fix the dielectric tensor.","Because $n^{1/2}$ also describes conventional three-dimensional electron gases and Luttinger semimetals, a quadratic-NLSM identification requires an independent measurement of the density of states or effective mass on the ring pocket, not just the density exponent.","The same density-exponent logic could be tested in classical wave analogues such as photonic or acoustic nodal-line crystals, where the carrier density maps to a tunable frequency parameter, before a clean electronic higher-order candidate is identified."],"forward_implications":["In any candidate quadratic or cubic nodal-line semimetal, the carrier-density exponent of the long-wavelength plasmon frequency is the discriminating observable: $n^{1/2}$, $n^{2/3}$, or $n^{3/4}$.","A high-resolution electron energy-loss measurement should resolve the $\\sqrt{2}$ doublet (around 190 and 269 meV for representative parameters) when the density is above about $10^{18}\\,\\mathrm{cm}^{-3}$ and the crystal is clean enough that the impurity broadening is below half the chemical potential.","The cubic crossover from $n^{2/3}$ to $n^{3/4}$ occurs near $n^*\\sim 10^{19}\\,\\mathrm{cm}^{-3}$, which is reachable by ionic-liquid gating or chemical doping; the finite-momentum dispersion is roughly quadratic, $\\Omega_p(q)^2 = \\Omega_p(0)^2[1+\\beta(q/k_F)^2]$, so the ratio test is approximately $q$-independent in the HREELS window.","For a general $p$-th order nodal-line semimetal in three dimensions, the same power counting predicts $\\omega_p\\propto n^{p/4}$ in the thin-ring limit and $\\omega_p\\propto n^{(p+1)/6}$ at large doping, interpolating from the linear case to the cubic case."],"supporting_citations":[{"why":"Supplies the symmetry-stabilized quadratic and cubic nodal-line Hamiltonians that the whole calculation starts from.","marker":"[17]"},{"why":"Establishes the divergent cubic density of states and the RG instability that limit the quantitative reliability of the $n^{3/4}$ law.","marker":"[19]"},{"why":"Gives the linear-NLSM plasmon scaling $n^{1/4}$ that the higher-order results extend.","marker":"[39]"},{"why":"Provides the earlier anisotropic density-response treatment of linear nodal-line semimetals used as reference.","marker":"[40]"},{"why":"Supplies the Luttinger-semimetal $n^{1/2}$ benchmark against which the quadratic prediction is read.","marker":"[43]"},{"why":"Provides the HREELS observation of nodal-line plasmons in ZrSiS that motivates the experimental protocol.","marker":"[50]"},{"why":"Documents the temperature-dependent nodal-line plasmon response in the ZrSiX family used as the experimental analogy.","marker":"[51]"}],"fun_headline_variants":["Plasmon scaling laws distinguish cubic from quadratic nodal-line semimetals","Cubic nodal-line plasmons show distinct density exponents","Plasmon frequency reveals quadratic vs cubic nodal-line semimetals","Density-scaling law distinguishes cubic from quadratic nodal-line plasmons","Plasmon density exponents distinguish cubic and quadratic nodal-line semimetals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The cubic nodal-line semimetal remains a semimetal in the $10^{18}$ to $10^{19}\\,\\mathrm{cm}^{-3}$ density window; if an excitonic gap or another interaction-driven instability opens instead, the predicted $n^{3/4}$ scaling no longer applies.","fun_headline_variants_meta":{"raw":{"variants":["Plasmon scaling laws distinguish cubic from quadratic nodal-line semimetals","Cubic nodal-line plasmons show distinct density exponents","Plasmon frequency reveals quadratic vs cubic nodal-line semimetals","Density-scaling law distinguishes cubic from quadratic nodal-line plasmons","Plasmon density exponents distinguish cubic and quadratic nodal-line semimetals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001115,"raw_usage":{"total_tokens":4859,"prompt_tokens":1376,"completion_tokens":3483,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":992,"completion_tokens_details":{"reasoning_tokens":3393}},"tokens_in":992,"tokens_out":3483,"duration_ms":23711,"temperature":1.0,"reasoning_tokens":3393,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:18:15.483374+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the bulk energy-loss function $-\\mathrm{Im}[1/\\epsilon(\\Omega)]$ of a candidate cubic nodal-line semimetal by reflection HREELS at densities $10^{18}$ to $10^{19}\\,\\mathrm{cm}^{-3}$, with the ring-pocket density fixed independently by Shubnikov-de Haas oscillations or ARPES. If the peak frequencies do not follow one of the predicted exponents ($n^{1/2}$, $n^{2/3}$, $n^{3/4}$), or if the anisotropy ratio after correcting for measured $\\epsilon_\\perp,\\epsilon_z$ is far from $\\sqrt{2}$, the RPA-based picture is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the symmetry-stabilized quadratic and cubic nodal-line Hamiltonians that the whole calculation starts from."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the divergent cubic density of states and the RG instability that limit the quantitative reliability of the $n^{3/4}$ law."},{"cited_title":"Panfilov, A","cited_arxiv_id":null,"evidence_quote":"Gives the linear-NLSM plasmon scaling $n^{1/4}$ that the higher-order results extend."},{"cited_title":"Zhou, H.-R","cited_arxiv_id":null,"evidence_quote":"Provides the earlier anisotropic density-response treatment of linear nodal-line semimetals used as reference."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Luttinger-semimetal $n^{1/2}$ benchmark against which the quadratic prediction is read."},{"cited_title":"Mauri and M","cited_arxiv_id":null,"evidence_quote":"Provides the HREELS observation of nodal-line plasmons in ZrSiS that motivates the experimental protocol."},{"cited_title":"Mandal, Search for plasmons in isotropic Luttinger semimetals, Ann","cited_arxiv_id":null,"evidence_quote":"Documents the temperature-dependent nodal-line plasmon response in the ZrSiX family used as the experimental analogy."}],"review_version":1}