{"id":"9e46ad61-face-42f5-b0de-97e2acb73d7c","arxiv_id":"2509.01311","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Hadronic vacuum polarization shifts in muonic atoms across the periodic table remain near 67% of the muonic vacuum polarization shift, with few-percent nuclear model dependence.","lead":"The authors calculate hadronic vacuum polarization corrections to energy levels in muonic atoms from Z=20 to Z=100 using a semi-empirical potential fitted to electron-positron data. They find the correction stays close to 67% of the muonic vacuum polarization shift and depends only weakly on the nuclear charge model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Missing specification of the Dirac wavefunctions in Eq. (7): point-nucleus vs finite-size eigenstates is load-bearing for the high-Z HVP shifts and the claimed ratio stability.","rationale":"The reader's weakest assumption—unstated finite-size treatment of the Dirac wavefunctions in Eq. (7)—is exactly the concern I would raise. The paper is otherwise a clean, internally consistent first-order perturbation calculation: the analytical potentials for the homogeneous sphere are a useful contribution, the Fermi-model comparison is a sensible model-dependence check, and the ratio to µVP is a meaningful diagnostic. However, because the high-Z HVP shift is a short-range overlap integral, the choice of wavefunction Hamiltonian is load-bearing. The text is ambiguous: 'fully relativistic energy shift calculations, including finite nuclear size effects' could describe only the perturbing potential, not the wavefunction generation. This is an addressable but unresolved detail, so the conditional verdict is appropriate. I do not see a reason to reject the paper on this basis, nor to accept it without the requested clarification/recalculation. Hence the verdict remains unchanged from the reader's CONDITIONAL.","tokens_in":5411,"tokens_out":5880,"duration_ms":80182,"concrete_test":"Recompute ΔE_HVP and ΔE_µVP for 1s1/2 at Z = 60, 82, 92, 100 using Eq. (7) with two wavefunction sets: (a) Dirac eigenstates of the Fermi finite-size nuclear potential with the same parameters as in the paper; (b) point-Coulomb Dirac eigenstates. Keep V_HVP fixed (Fermi model) in both cases. Report ΔE_HVP(b)/ΔE_HVP(a) and ΔE_HVP/ΔE_µVP for both wavefunction sets. If the wavefunction choice changes the 1s HVP shift by more than ~2% at Z = 100, or shifts the ratio outside the stated 0.60–0.70 band, the missing specification is decisive. As a secondary check, repeat for 2s1/2 at Z = 90, where the sphere/Fermi model difference is already 5–6%.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All numerical results pass through Eq. (7), the first-order expectation value of V_HVP(r) in Dirac radial wavefunctions G and F. The paper never states which Dirac Hamiltonian produced those wavefunctions—specifically, whether the Coulomb potential used to generate the bound states contains the same finite nuclear charge distribution used in V_HVP. This is not a pedantic detail. In muonic atoms at Z ≳ 60, the 1s wavefunction is strongly pulled into the nuclear volume; finite nuclear size changes the small-r normalization and the energy. V_HVP is short-ranged and largest inside or right at the nuclear surface (Eqs. (3), (5), (6)), so the overlap integral in Eq. (7) is sensitive to the near-origin behavior of G and F. If point-nucleus Dirac wavefunctions were used, the high-Z HVP shifts in Fig. 1 would differ materially, and the claimed 'irregularities that correlate with nuclear radii' would arise only from the perturbing potential, not from a consistent solution of the bound muon in the same finite-size field. The phrase 'fully relativistic energy shift calculations, including finite nuclear size effects' is ambiguous: it may mean only the perturbing potential includes the finite charge distribution, leaving the wavefunctions point-nucleus. This unstated modeling choice determines whether the central numerical claim survives at high Z.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes first-order hadronic vacuum polarization (HVP) energy shifts for 1s1/2, 2s1/2, 2p1/2, and 2p3/2 states of muonic atoms with Z = 20–100, using a semi-empirical HVP potential fitted to e+e− → hadrons data. Two nuclear charge distributions are compared: a two-parameter Fermi model (numerical integration) and a homogeneously charged sphere model (partly analytical). The authors report that the HVP shift grows with Z, shows high-Z irregularities correlated with nuclear charge radii, and that the HVP-to-muonic-VP ratio stays near 0.60–0.70, within 10% of the nonrelativistic 67% benchmark. They conclude that their results improve on the earlier Sundaresan-Watson treatment.","tokens_in":5730,"tokens_out":1432,"duration_ms":19963,"significance":"If the central numerical claims are correct, the paper provides a systematic reference dataset for HVP corrections in medium- and high-Z muonic atoms, extending earlier work that focused on low-Z systems. The two-model comparison and the stability analysis of the HVP/muonic-VP ratio are useful for ongoing muonic-atom spectroscopy programs. The approach is transparent and the analytic reduction for the sphere model is a practical asset. The main strength is the breadth of the survey and the explicit correlation of high-Z irregularities with measured nuclear radii. However, the paper does not establish the reliability of the underlying Dirac wavefunctions, gives no uncertainty estimates, and does not quantitatively benchmark against the earlier result it claims to improve.","major_comments":[{"comment":"The manuscript does not state how the Dirac radial wavefunctions G and F entering Eq. (7) are generated. This is load-bearing: at Z ≳ 60 the muon 1s wavefunction is strongly modified by finite nuclear size, and V_HVP is short-ranged and peaked near or inside the nuclear surface. If point-nucleus Dirac wavefunctions were used, the high-Z shifts in Fig. 1 would differ materially, and the claimed correlation of the shifts with nuclear radii would be an artifact of the perturbing potential only. Please specify the Dirac Hamiltonian (including the nuclear charge distribution) used to obtain G and F, and ideally demonstrate the sensitivity by comparing point-nucleus and finite-size wavefunctions for a few high-Z cases.","section":"Theory, Eq. (7)"},{"comment":"No uncertainties are reported for any numerical result. The final shifts are quotes to many digits (e.g., −250 eV scale in Fig. 1), yet the input HVP fit parameters B1 and C1 carry fit errors from the e+e− analysis, and the Fermi-model integration and radii data introduce additional uncertainty. A quantitative error budget is needed before these results can serve as a reference for high-precision muonic spectroscopy.","section":"Results and Discussion, Fig. 1"},{"comment":"The paper claims 'improved theoretical predictions for low-lying energy levels in muonic systems compared to [18]', but no numerical comparison with Sundaresan-Watson is shown. The reader cannot verify whether the differences are due to updated HVP parameters, different wavefunctions, or different nuclear models. A table comparing representative Z values with Ref. [18] would make the improvement claim concrete and falsifiable.","section":"Introduction / Conclusions"}],"minor_comments":[{"comment":"The derivation from Eq. (2) to Eq. (3) is only sketched. The integral over q is not shown, and the sign convention for the exponential integral E_n may confuse readers. A short derivation in an appendix would improve reproducibility.","section":"Theory, Eq. (3)"},{"comment":"The caption of Fig. 1 defines δ but the main text references it only implicitly. Please state explicitly that δ is the relative difference between Fermi and sphere models and clarify its sign convention in the main text.","section":"Results and Discussion, Eq. (8)"},{"comment":"The conclusion states that HVP reaches up to 68% of muonic VP, but Fig. 1 shows ratios near 0.675–0.68 for s states and lower values for p states. The sentence should distinguish the state-dependence, otherwise it overstates a single representative value.","section":"Conclusions"},{"comment":"Equation (1) and the parameter values A1, B1, C1 are attributed to Refs. [17,19], but the reader must infer that A1=0 from Ref. [19]. Please indicate explicitly that these are the low-energy Burkhardt-Pietrzyk parameters and specify their fit uncertainty.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The missing specification of the Dirac wavefunctions is the most consequential issue and should be resolved before publication. The paper otherwise appears technically sound within its stated scope, but the absence of any uncertainty estimate and the lack of a quantitative comparison with the earlier Sundaresan-Watson work weaken the 'improved predictions' claim. I recommend major revision rather than rejection because the central ratio-stability result is likely robust to the wavefunction choice for s-states, but the high-Z absolute values and the claimed nuclear-radius correlation require the wavefunction clarification and sensitivity check."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a clean, practical calculation: HVP shifts for 1s, 2s, 2p1/2, 2p3/2 in muonic atoms from Z=20 to 100, using the Burkhardt–Pietrzyk semiphenomenological hadronic polarization function and two nuclear charge models. The analytic sphere formulas are a nice touch, and the Fermi-model numerical integration is a sensible cross-check. The main new result is systematic: the HVP/muVP ratio stays within 0.60–0.70 across the whole Z range, with s-states near 0.68 and p-states climbing from 0.62 toward the nonrelativistic 67%—that is genuinely new and useful for muonic-atom spectroscopy, even if it doesn't change the theoretical framework.\n\nThe soft spot the stress-test flagged is real. Equation (7) evaluates the HVP potential in Dirac wavefunctions G and F, but the paper never says what Hamiltonian produced those wavefunctions. The abstract claims the calculation includes finite nuclear size, and the perturbing potential is built from a finite-size charge distribution, but that's not the same as using Dirac eigenstates of a finite-size Coulomb potential. At Z around 100 the muon 1s wavefunction is strongly pulled into the nucleus; the HVP potential is short-ranged and peaks right there. If those wavefunctions are point-nucleus, the high-Z shifts and the 'irregularities' in Fig. 1 would change materially. The difference matters, and the paper needs to state it explicitly.\n\nOther issues are milder but still worth fixing. There are no error bars anywhere—the B1, C1 fit parameters carry uncertainties, the experimental charge radii carry uncertainties, and those propagate directly into the shifts. The comparison with Sundaresan and Watson (1975) is asserted but not quantified, so 'improved' is hard to evaluate. No code or tabulated data is provided, which limits reuse.\n\nNone of this sinks the paper. The calculation is straightforward, the two nuclear models agree to a few percent, and the ratio stability is a solid empirical observation. It just needs a revision that pins down the wavefunction prescription, adds uncertainties, and gives the earlier comparison concretely.\n\nThis deserves a serious referee. It's not a breakthrough, but it's a reference-quality calculation for a community that will use these numbers.","headline":"Useful systematic scan of HVP shifts in muonic atoms, but the unstated choice for the Dirac wavefunctions—point nucleus or finite-size—needs to be pinned down before the high-Z numbers can be trusted.","tokens_in":6197,"tokens_out":1336,"would_cite":true,"duration_ms":19664,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"For muonic atoms from Z=20 to Z=100, hadronic vacuum polarization shifts grow with nuclear charge, track measured nuclear-radius fluctuations, and stay within 10% of the 67% non-relativistic ratio to the muonic vacuum polarization.","keywords":["hadronic vacuum polarization","muonic atoms","vacuum polarization ratio","nuclear charge radii","Dirac equation","finite nuclear size","energy levels","QED corrections"],"falsifier":"Recompute the 1s HVP shift for Z=90-100 with Dirac wavefunctions obtained from a Fermi-distribution nucleus versus a point nucleus; if the two differ by more than about 30% at Z=100, the quoted values are not robust to the wavefunction model. In parallel, measure the 2s-2p Lamb shift of a heavy muonic atom (e.g., Z about 80) to enough precision that the HVP contribution is resolved; if the implied HVP-to-muVP ratio deviates by more than 10% from 0.67, the stability claim fails.","tokens_in":5341,"feed_emoji":"⚛️","tokens_out":6517,"duration_ms":68470,"temperature":0.7,"pith_summary":"This paper computes hadronic vacuum polarization (HVP) energy shifts for the 1s, 2s, 2p1/2, and 2p3/2 states of muonic atoms with nuclear charges Z=20 to 100, using a semi-empirical potential fitted to e+e- annihilation data and two nuclear charge distribution models. It finds that HVP shifts grow with Z but show irregular spikes that disappear when a smooth radius formula replaces measured nuclear charge radii, so the spikes trace the non-monotonic experimental radii. The paper then compares HVP with the muonic vacuum polarization (muVP) shift, finding the ratio stays between 0.60 and 0.70 for all studied states and Z values, within 10% of the non-relativistic 67% benchmark. The two nuclear models agree to a few percent for most nuclei, with differences up to 5-6% for the heaviest 2s states.","feed_headline":"Hadronic vacuum polarization ratio in muonic atoms holds near 67%","feed_subtitle":"Relativistic calculation finds HVP reaches 68% of muon-VP at high Z; spikes trace nuclear radii.","key_machinery":"The machinery is the semi-empirical hadronic Uehling potential, Eq. (2)/(3): a one-dimensional integral over the hadronic polarization function, parameterized piecewise by constants A, B, C fitted to e+e- to hadrons data and restricted here to the low-energy (0-0.7 GeV) domain, folded with the nuclear charge density rho(x). For a homogeneously charged sphere the integral reduces to closed exponential-integral expressions, Eqs. (5)-(6); for the Fermi distribution it is evaluated numerically. First-order energy shifts come from the radial overlap of this potential with Dirac wavefunctions, Eq. (7). The ratio to the muonic vacuum polarization (Uehling) shift provides the benchmark comparison.","core_discovery":"The central claim is that HVP corrections to energy levels of muonic atoms can be reliably computed across the whole periodic table with a dispersion-relation-based potential, and that two quantitative results hold. First, the Z-dependence of the HVP shift tracks the measured nuclear charge radii, not a smooth Z^(1/3) law: the observed spikes in the high-Z region correlate directly with the non-monotonic behavior of experimental nuclear charge radii, and replacing the radii by a smooth formula restores a monotonic profile. Second, the ratio of HVP to muonic vacuum polarization remains remarkably stable, within 10% of the established non-relativistic 67% value, despite strong relativistic eff","pith_inferences":["The paper leaves unspecified whether the Dirac wavefunctions in Eq. (7) are eigenstates of a Hamiltonian that includes the finite nuclear charge distribution; redoing the calculation with wavefunctions computed in a finite-size Fermi potential versus a point nucleus would likely change the high-Z 1s numbers, so this should be pinned down before using the values in fits.","The stability of the HVP-to-muVP ratio suggests the two corrections share the same short-distance weighting; if future higher-order hadronic corrections preserve this ratio, HVP uncertainties could be absorbed into an effective vacuum-polarization coupling in heavy muonic atoms.","The correlation with nuclear radii could be turned around: precise muonic-atom HVP measurements at high Z, combined with muVP and other QED corrections, might serve as a probe of the nuclear charge radius independent of electron scattering.","The method is restricted to q^2 below 0.7 GeV, so higher-energy hadronic contributions are dropped; if those matter at the percent level, the quoted few-percent model agreement may not cover the full hadronic uncertainty."],"forward_implications":["For light and medium muonic atoms, using either the Fermi or homogeneous-sphere nuclear model changes HVP shifts by less than a few percent, so the simpler sphere model is adequate there.","At high Z the HVP shift is isotope-sensitive: measured nuclear radius fluctuations imprint directly on the shift, so precise calculations need the actual radius of the isotope under study rather than a smooth global formula.","The HVP-to-muVP ratio stays within 10% of 0.67 across Z=20-100 for the four lowest states, so the non-relativistic 67% benchmark remains a reasonable planning baseline, but p-states at low Z sit about 5-7% lower and require state-specific values.","For heavy elements the HVP contribution reaches up to 68% of the muVP shift, making HVP a non-negligible part of the Lamb-shift budget in precision muonic x-ray spectroscopy."],"supporting_citations":[{"why":"Supplies the semi-empirical dispersion-relation method and potential parametrization for hadronic vacuum polarization that this paper extends to a systematic Z range.","marker":"[16]"},{"why":"Provides the e+e- to hadrons cross-section data that the hadronic polarization function is fitted to.","marker":"[17]"},{"why":"Gives the low-energy fit parameters (A1=0, B1=0.0023092, C1=3.9925370 GeV^-2) used in the adopted potential.","marker":"[19]"},{"why":"Supplies the table of experimental nuclear ground-state charge radii whose non-monotonic Z-behavior explains the observed spikes in HVP shifts.","marker":"[20]"},{"why":"Establishes the non-relativistic 67% HVP-to-muonic-VP benchmark ratio that this paper tests at high Z.","marker":"[21]"},{"why":"Prior calculation of HVP corrections in muonic atoms that the present work updates and improves upon.","marker":"[18]"}],"fun_headline_variants":["Muonic atom HVP shifts mirror nuclear radii spikes","HVP-to-muon VP ratio holds near 67% across Z","Hadronic VP in muonic atoms: radii explain Z spikes","Relativistic HVP ratio stable despite high-Z effects","Nuclear radii drive hadronic VP irregularities in muonic atoms"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is that the Dirac radial wavefunctions used in Eq. (7) adequately represent the muon in the actual finite-size nuclear potential; the paper does not state whether they include the finite nuclear charge distribution, and for Z near 100 the 1s wavefunction is strongly modified by nuclear size, so the numerical HVP shifts depend on this silent choice.","fun_headline_variants_meta":{"raw":{"variants":["Muonic atom HVP shifts mirror nuclear radii spikes","HVP-to-muon VP ratio holds near 67% across Z","Hadronic VP in muonic atoms: radii explain Z spikes","Relativistic HVP ratio stable despite high-Z effects","Nuclear radii drive hadronic VP irregularities in muonic atoms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000133,"raw_usage":{"total_tokens":924,"prompt_tokens":647,"completion_tokens":277,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":391,"completion_tokens_details":{"reasoning_tokens":206}},"tokens_in":391,"tokens_out":277,"duration_ms":3592,"temperature":1.0,"reasoning_tokens":206,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T12:38:49.195161+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the 1s HVP shift for Z=90-100 with Dirac wavefunctions obtained from a Fermi-distribution nucleus versus a point nucleus; if the two differ by more than about 30% at Z=100, the quoted values are not robust to the wavefunction model. In parallel, measure the 2s-2p Lamb shift of a heavy muonic atom (e.g., Z about 80) to enough precision that the HVP contribution is resolved; if the implied HVP-to-muVP ratio deviates by more than 10% from 0.67, the stability claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the semi-empirical dispersion-relation method and potential parametrization for hadronic vacuum polarization that this paper extends to a systematic Z range."},{"cited_title":"Borie and G","cited_arxiv_id":null,"evidence_quote":"Provides the e+e- to hadrons cross-section data that the hadronic polarization function is fitted to."},{"cited_title":"Burkhardt, F","cited_arxiv_id":null,"evidence_quote":"Gives the low-energy fit parameters (A1=0, B1=0.0023092, C1=3.9925370 GeV^-2) used in the adopted potential."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the table of experimental nuclear ground-state charge radii whose non-monotonic Z-behavior explains the observed spikes in HVP shifts."},{"cited_title":"Burkhardt and B","cited_arxiv_id":null,"evidence_quote":"Establishes the non-relativistic 67% HVP-to-muonic-VP benchmark ratio that this paper tests at high Z."},{"cited_title":"Breidenbach, E","cited_arxiv_id":null,"evidence_quote":"Prior calculation of HVP corrections in muonic atoms that the present work updates and improves upon."}],"review_version":1}