{"id":"93e173bb-5dfb-4f23-9eac-e6d1743a3148","arxiv_id":"2607.27096","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"Skyrme-RPA spectra for Ca, Ni, Zr and Sn reduce Coulomb NP g-factor uncertainties to 10–17%, agreeing with Exp+EWSR within a few half-ranges.","lead":"Microscopic Skyrme-RPA nuclear models tighten the uncertainty on the Coulomb nuclear-polarization correction to bound-electron g factors of four hydrogenlike ions to roughly 10–17%. This matters because NP uncertainty is now a limiting factor for using precision g-factor measurements to test QED or extract nuclear radii in heavy ions.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The central claim is modest and carefully scoped: a multi-Skyrme RPA scan tightens the Coulomb LO NP g-factor uncertainty to ~10–17% for four spherical ions, with averages lying close to the traditional Exp+EWSR route (Table III). The reader correctly flags that the uncertainty recipe itself is a modeling choice rather than a first-principles error bar; that is already acknowledged in §III–IV and is why CONDITIONAL is appropriate. No independent technical failure (wrong diagrams, incomplete multipoles, inconsistent radial functions, or arithmetic inconsistency in the tables) appears on a second pass. Transverse NP and deformed nuclei are explicitly left open. Therefore the verdict needs no adjustment.","tokens_in":12866,"tokens_out":369,"duration_ms":7834,"concrete_test":"Recompute one ion (e.g. 120Sn 1s) with an independent RPA implementation or a non-Skyrme EDFs (relativistic mean-field or Gogny) using the same radial form factors (16)–(17); if the new central value lies outside the Table III Skyrme half-range, the 10–12% band is under-estimated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's weakest-assumption note (that half-range Skyrme scatter plus Exp+EWSR offset may miss a common systematic bias) is already the paper's own stated limitation and is the sole reason for the CONDITIONAL verdict. Within the Coulomb-only, spherical-nucleus scope the manuscript actually claims, the numerics, multipole breakdown, and comparison to Exp+EWSR are internally consistent and do not rest on a hidden derivation error or circular construction. No stronger load-bearing flaw is present.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript evaluates the Coulomb part of the leading-order nuclear-polarization (NP) correction to the bound-electron g factor of hydrogenlike 40Ca19+, 60Ni27+, 90Zr39+, and 120Sn49+ by supplying nuclear excitation energies and reduced transition probabilities from Hartree–Fock RPA with nine Skyrme forces. These spectra enter effective self-energy diagrams in which the photon propagator carries an NP insertion (Eqs. 15–25). Results are compared multipole-by-multipole and part-by-part (irreducible/reducible/residual/pole) with the conventional combination of experimental low-lying levels plus energy-weighted-sum-rule (EWSR) giant-resonance estimates. The authors conclude that the model scatter, together with the offset to Exp+EWSR, tightens the relative theoretical uncertainty of the Coulomb NP correction to roughly 10–12% (17% for 40Ca), substantially below the 30–50% figures customarily quoted.","tokens_in":12987,"tokens_out":1163,"duration_ms":31136,"significance":"NP uncertainties are already comparable to, or larger than, residual two-loop QED and finite-nuclear-size errors in the medium- and high-Z g-factor program. A controlled reduction of that uncertainty, even if restricted to the Coulomb sector and to spherical closed-shell nuclei, directly improves the interpretability of existing and planned high-precision measurements (e.g., Sn49+) and strengthens the case for using electronic g factors as complementary probes of nuclear charge radii. The work supplies concrete numerical tables, an explicit multipole hierarchy, and a transparent comparison between microscopic Skyrme-RPA and the widely used Exp+EWSR recipe; these are reusable benchmarks for the community.","major_comments":[{"comment":"§III and Table III: the quoted 10–17% uncertainties are obtained by taking the average of nine Skyrme results as the central value and the absolute difference to the Exp+EWSR entry as the (one-sided) error bar. This prescription is not derived from a statistical or Bayesian model of nuclear-theory error; it assumes that the half-range among the selected Skyrme set plus the Exp+EWSR offset already brackets all relevant systematics. Because every Skyrme-RPA calculation and the EWSR centroids share the same mean-field/RPA framework and the same homogeneous-sphere radial form factors (Eqs. 16–17), a common bias cannot be excluded by construction. The central claim of “tighter constraints” therefore rests on an ad-hoc uncertainty rule that should be stated more explicitly as a convention, or supplemented by at least one independent nuclear approach (e.g., a relativistic mean-field RPA or a di","section":"§III, Table III"},{"comment":"The entire analysis is restricted to the Coulomb (00) component of the NP-modified photon propagator. The abstract and conclusions correctly flag the transverse contribution as open, yet the numerical claims are phrased as uncertainties on “the nuclear-polarization corrections.” Given that earlier literature (Refs. [19,36]) has reported non-negligible or even counter-intuitive transverse pieces in related systems, a short quantitative estimate—or a clear statement that the 10–17% figures apply strictly to the Coulomb sector and must be enlarged once transverse terms are known—should appear wherever the final uncertainty numbers are quoted (Abstract, end of §III, §IV).","section":"Abstract; §IV"}],"minor_comments":[{"comment":"Table III, row label “Range/25.2”: this is almost certainly a typesetting concatenation of “Range/2” with the first numerical entry 5.2. Please separate the label from the data columns.","section":"Table III"},{"comment":"Table I and II caption and body use “g-factor” with a hyphen while the title and most of the text use “g factor”; pick one convention.","section":"Tables I–II"},{"comment":"Eqs. (26)–(27): the isospin label τ is introduced without a brief reminder that τ=0 (1) denotes isoscalar (isovector); a one-line clarification would help non-nuclear readers.","section":"§II.D"},{"comment":"The radius R0 of the homogeneous sphere that defines FJ(x) is never varied. A one-sentence sensitivity check (or a pointer to prior work) would reassure the reader that this choice is sub-dominant compared with the Skyrme scatter already quoted.","section":"§II.B, Eqs. (16)–(17)"},{"comment":"Reference [11] is listed as “submitted, arXiv:2202.01668”; the arXiv number appears inconsistent with a 2025 submission date—please verify the identifier.","section":"References"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid, incremental advance that fits physics.atom-ph well. The only load-bearing softness is the informal uncertainty rule; once the authors label it clearly as a convention and restrict the abstract numbers to the Coulomb sector, the paper is acceptable. No concerns about novelty disclosure or citation manipulation."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one thing worth knowing is that they replace the usual 30–50% hand-waving on the Coulomb NP correction to the bound-electron g factor with a concrete 10–17% band for four spherical ions, obtained by feeding nine Skyrme-RPA spectra into the standard effective-self-energy diagrams. That is useful for the precision g-factor and electronic-radius community, even if it is not a foundational advance.\n\nWhat is new is the systematic nuclear-model scan, not the QED. The formalism (NP-modified photon propagator, irreducible/reducible/residual/pole pieces) is taken from Plunien, Nefiodov, and their own earlier papers; the expressions in Eqs. 18–25 look standard. Tables I–II cleanly show that the irreducible 1^{-} piece dominates and that J=4,5 still add a couple of percent. Table III is the real payload: nine Skyrme forces (the same conservative set used in their 2022 muonic work) versus Exp+EWSR, with averages and half-ranges that sit within two–three half-ranges of the older recipe. The agreement is reassuring and suggests the giant-resonance contribution is mostly an integral property.\n\nSoft spots are real but proportionate. The uncertainty recipe itself (half-range scatter plus offset to Exp+EWSR) is a modeling choice; a common Skyrme/EWSR bias could still be hiding, and they only treat the Coulomb part for spherical closed-shell nuclei. They say so. No derivation error or circular construction appears. Citations are appropriate; mild self-citation is just continuity of the same program.\n\nThis is for people who actually need the next digit on heavy H-like g factors or who want a benchmark before trusting Exp+EWSR on deformed systems. It deserves a serious referee. I would cite the Table III numbers if I were quoting NP uncertainties in the next year, and I would accept it for review without hesitation.","headline":"Solid incremental work that tightens Coulomb NP g-factor uncertainties to ~10–17% via multi-Skyrme RPA; the numbers are transparent and the main caveat is already stated by the authors.","tokens_in":13664,"tokens_out":491,"would_cite":true,"duration_ms":9413,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Skyrme nuclear models cut the uncertainty on the nuclear-polarization correction to bound-electron g factors down to about 10–17%.","keywords":["bound-electron g factor","nuclear polarization","Skyrme interaction","random-phase approximation","hydrogenlike ions","energy-weighted sum rules","QED tests"],"falsifier":"A complete experimental multipole strength distribution (or an independent microscopic calculation outside the Skyrme-RPA family) for one of the four nuclei that moves the nuclear-polarization correction outside the quoted 10–17% band relative to the Skyrme average.","tokens_in":13677,"feed_emoji":"⚛️","tokens_out":955,"duration_ms":18177,"temperature":0.7,"pith_summary":"High-precision g-factor measurements in heavy hydrogenlike ions are limited by nuclear polarization: the electron feels the nucleus’s fluctuating excitations, not just its static size. That correction is small but its uncertainty has been estimated at 30–50%, largely because experimental nuclear data are incomplete and giant resonances are treated with rough sum-rule estimates. This paper replaces that patchwork with microscopic spectra from nine Skyrme interactions computed in the Hartree–Fock random-phase approximation, for 40Ca, 60Ni, 90Zr and 120Sn. The Skyrme results cluster tightly and sit close to the usual experimental-plus-sum-rule values, so the authors argue the Coulomb part of the leading-order correction can be quoted with relative uncertainties of only 10–12% (17% for calcium). That tighter bound is offered as practical input for ongoing and future g-factor experiments that aim to test QED or extract nuclear radii.","feed_headline":"Nuclear models cut g-factor uncertainty to 10–17%","feed_subtitle":"Skyrme spectra tighten the nuclear-polarization error that limits precision tests in heavy ions","key_machinery":"Nuclear-polarization insertion in the photon propagator: the fluctuating nuclear four-current supplies a correction D_NP to the photon line, built from Skyrme-RPA excitation energies and reduced transition probabilities B(EJ); that modified propagator is then inserted into the one-loop effective self-energy diagrams that shift the bound-electron g factor.","core_discovery":"Across nine realistic Skyrme parameterizations, the Coulomb leading-order nuclear-polarization corrections to the bound-electron g factor of 60Ni27+, 90Zr39+ and 120Sn49+ agree with one another and with the conventional experimental-plus-energy-weighted-sum-rule estimates closely enough that the relative theoretical uncertainty can be reduced to about 10–12% (17% for 40Ca19+), taking the Skyrme average as the central value and the offset to the sum-rule result as a conservative absolute uncertainty.","pith_inferences":["If the transverse NP piece later proves non-negligible, the same Skyrme spectra can be reused to bound it without waiting for new nuclear data.","The slow multipole convergence (J=4,5 still add ~2.5%) suggests future work should routinely keep higher multipoles once the uncertainty floor drops below a few percent.","Agreement between fragmented RPA spectra and few-state EWSR estimates implies that integral sum-rule strength, not spectral fine structure, dominates the electronic NP shift—useful for quick estimates across the nuclear chart."],"forward_implications":["Quoted NP uncertainties for 60Ni27+, 90Zr39+ and 120Sn49+ g factors can be tightened from the traditional 30–50% level to roughly 10–12%.","Only the irreducible parts of the a/b self-energy diagrams are needed for a reliable leading-order estimate; the c-diagram is negligible at the present precision.","The simple experimental-plus-EWSR recipe is vindicated as a proxy when microscopic spectra are unavailable, including for deformed nuclei.","Tighter NP errors strengthen the case for using heavy hydrogenlike g factors to extract nuclear charge radii or to cross-check muonic-atom results."],"fun_headline_variants":["Skyrme models cut nuclear-polarization g-factor uncertainty to 10-17%","Nuclear Skyrme spectra tighten bound-electron g-factor error bounds","Microscopic Skyrme force constrains g-factor corrections to 10-12%","Nine Skyrme sets reduce nuclear-polarization uncertainty in heavy ions","Skyrme averages yield 10-17% bounds on nuclear g-factor corrections"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The spread among nine chosen Skyrme models, plus their offset from the experimental-plus-sum-rule estimate, is assumed to capture the true theoretical uncertainty rather than a shared systematic bias of all such models.","fun_headline_variants_meta":{"raw":{"variants":["Skyrme models cut nuclear-polarization g-factor uncertainty to 10-17%","Nuclear Skyrme spectra tighten bound-electron g-factor error bounds","Microscopic Skyrme force constrains g-factor corrections to 10-12%","Nine Skyrme sets reduce nuclear-polarization uncertainty in heavy ions","Skyrme averages yield 10-17% bounds on nuclear g-factor corrections"]},"model":"grok-4.5","effort":"low","cost_usd":0.00438,"raw_usage":{"total_tokens":1305,"prompt_tokens":808,"num_sources_used":0,"completion_tokens":84,"cost_in_usd_ticks":43804000,"prompt_tokens_details":{"text_tokens":808,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":413,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":808,"tokens_out":84,"duration_ms":8777,"temperature":1.0,"reasoning_tokens":413,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-30T11:14:23.876591+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A complete experimental multipole strength distribution (or an independent microscopic calculation outside the Skyrme-RPA family) for one of the four nuclei that moves the nuclear-polarization correction outside the quoted 10–17% band relative to the Skyrme average.","supporting_citations":[],"review_version":1}