{"id":"5e5433e1-0ab3-4fac-89f8-8269440f8412","arxiv_id":"2607.07658","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":1,"one_line_summary":"The hVP contribution to the HFS in muonic hydrogen is 2.153(11) µeV, deviating from previous evaluations by ~10x the anticipated experimental precision, due to corrected recoil and finite-size interplay.","lead":"This paper recalculates the hadronic vacuum polarization (hVP) contribution to the Lamb shift and hyperfine splitting in hydrogen-like atoms and ions, finding significant deviations from previous HFS evaluations. A smart generalist reads this because upcoming precision spectroscopy experiments (CREMA, FAMU) will test these predictions, potentially impacting proton radius and fundamental constant determinations.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"FF model dependence is the right concern but is not load-bearing for the central claim; the deviation from previous work (~2.6 µeV) vastly exceeds the 4% FF model uncertainty (~0.086 µeV). Verdict unchanged.","rationale":"The reader's verdict of ACCEPT with HIGH confidence is appropriate. The central claim — that the hVP contribution to µH HFS is 2.153(11) µeV, deviating from previous evaluations by roughly ten times the anticipated experimental precision — is well-supported by: (1) a clean dispersive formalism with multiple limiting-case checks, (2) agreement with the literature for the Lamb shift and Mu HFS, (3) concrete identification of errors in prior work [9], and (4) transparent discussion of systematic uncertainties even though they are not all included in the tabulated errors. The weakest assumption (proton FF model dependence at 4%) is acknowledged and does not threaten the central finding because the deviation from previous work is ~30× larger than the FF model uncertainty. The fully inflated uncertainty (~0.11 µeV including e+e− data scatter and FF model dependence) remains much smaller than the ~2.6 µeV deviation from Borie's result. No internal inconsistency or critical red flag was identified. The paper's physical argument for recoil suppression by nuclear FFs is clearly articulated and quantitatively supported by Fig. 2 and Table A1.","tokens_in":20466,"tokens_out":4002,"duration_ms":391482,"concrete_test":"Recompute the µH HFS result using the dispersion-theoretical proton FFs of Lin et al. [Ref. 62] in Eq. (10) with the same DHMZ R-ratio input. The paper states this yields a ~3.7% larger correction. If the shift is indeed ~3.7% (i.e., ~0.080 µeV), the 4% model uncertainty assignment is validated and the central claim of ~10× experimental precision deviation is unaffected. If the shift exceeds 10% (>0.2 µeV), the FF model dependence would become comparable to the experimental precision and would need to be included in the tabulated error budget.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the proton FF model dependence as the weakest assumption. The paper uses Borah et al. [Ref. 26] FFs (constrained by the µH Lamb shift charge radius) and acknowledges a 4% model uncertainty versus the dispersion-theoretical FFs of Lin et al. [Ref. 62], which is not included in the tabulated errors. However, this concern does not land as load-bearing for the central claim of significant deviation from previous evaluations. The deviation from Borie's result [Ref. 66] is |4.8 − 2.153| ≈ 2.6 µeV, while the 4% FF model uncertainty corresponds to ~0.086 µeV — a factor of ~30 smaller. Even combining all untabulated systematic uncertainties (e+e− data scatter inflation factor ~6.1 giving ~0.067 µeV, plus the 4% FF model uncertainty ~0.086 µeV), the total inflated uncertainty is ~0.11 µeV, still far below the ~2.6 µeV deviation. The master formula Eq. (10) is supported by multiple consistency checks: it reduces to known results in the pointwise [Eq. (13)], structureless [Eq. (14)], and non-recoil [Eq. (12)] limits; the Lamb shift analog [Eq. (23a)] agrees with previous calculations; the Mu HFS result agrees with the literature (Table I); and the identification of specific calculational errors in [9] (extra factor 2 on F₂, factors 1/2 and M/2m on Eqs. 25 and 26) is concrete and independently verifiable. The most substantive residual concern is that the quoted uncertainty 0.011 µeV substantially understates the full systematic budget, but this is a presentation issue rather than a threat to the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript evaluates hadronic vacuum polarization (hVP) contributions to the Lamb shift and hyperfine splitting (HFS) in ordinary and muonic hydrogen and hydrogen-like helium-3 ions, using the dispersive data-driven approach with the DHMZ parametrization of the R-ratio. The central physical result is that nuclear elastic form factors suppress the recoil corrections to the hVP-HFS that are large in muonium (where both constituents are pointlike), rendering them negligible in composite-nucleus systems. The authors obtain 2.153(11) µeV for the hVP contribution to the µH ground-state HFS, deviating from previous evaluations by roughly ten times the anticipated CREMA/FAMU experimental precision. They attribute the discrepancies to (i) an incorrect hVP-to-µVP rescaling used in prior work and (ii) specific calculational errors in Ref. [9] (Faustov and Martynenko). Lamb shift results are shown to agree with the literature, and a first evaluation of the subleading O(Z^5 α^6) hVP–finite-size correction is presented.","tokens_in":21301,"tokens_out":1331,"duration_ms":232411,"significance":"The timing is excellent: CREMA and FAMU are pursuing the µH ground-state HFS at ~1 ppm precision, and a factor-of-ten discrepancy in the hVP contribution relative to prior evaluations is directly relevant to the interpretation of those measurements. The formalism is laid out clearly, with the master formula Eq. (10) supported by multiple consistency checks: it reduces to known results in the pointlike [Eq. (13)], structureless [Eq. (14)], and non-recoil [Eq. (12)] limits. The Mu HFS result (Table I) agrees with independent evaluations, providing a validation of the data-driven pipeline. The identification of specific errors in Ref. [9] (extra factor of 2 on F_2, missing factors of 1/2 and M/2m in Eqs. 25 and 26) is concrete and independently verifiable. The first evaluation of the O(Z^5 α^6) hVP–finite-size correction for the Lamb shift, particularly the 6.7(2) µeV effect in µ³He⁺, is a genuine addition. No free parameters are fitted; the results depend on empirical R-ratio data and external nuclear form factor parametrizations.","major_comments":[{"comment":"Sec. V, Table II: The quoted uncertainty on the µH HFS result, 0.011 µeV, does not include the 4% proton FF model uncertainty (~0.086 µeV) that the authors themselves identify in Sec. V. The tabulated error thus substantially understates the full systematic budget. While the central claim of significant deviation from previous work remains valid—the ~2.6 µeV discrepancy with Borie's result far exceeds even the inflated uncertainty—the authors should either include the FF model uncertainty in the tabulated error or, at minimum, state the total inflated uncertainty explicitly in the table caption. As written, a reader taking Table II at face value would underestimate the uncertainty by nearly an order of magnitude.","section":null},{"comment":"Sec. V: The treatment of the e⁺e⁻ data scatter and the a_µ discrepancy inflation factor (2.44) is described in prose but not propagated into the tabulated uncertainties. The text states that accounting for experimental scatter gives an extra factor of 2.51 for (µ)H and 1.06 for (µ)³He⁺, followed by an additional inflation factor of 2.42 and 1.25, respectively. It would strengthen the paper to provide a compact summary table or footnote giving the final inflated uncertainty for each system, so that the full systematic budget is transparent and not buried in narrative arithmetic from the prose.","section":null}],"minor_comments":[{"comment":"Sec. III: The criterion for when recoil corrections may be neglected is stated qualitatively ('the VP spectral function extends to scales comparable with the heavier mass, and no FF cutoff intervenes at a lower scale'). A more quantitative version, e.g., specifying the ratio √t₀/M or the FF cutoff scale Λ_FF, would be helpful.","section":null},{"comment":"Fig. 2(b): The curves for W(t) and W_non-recoil(t) are stated to be 'on top of each other,' but the figure as presented makes it difficult to assess the actual size of the residual difference. A ratio plot or an inset showing the fractional difference would improve clarity.","section":null},{"comment":"Table A1: The caption could clarify that the 'finite size' and 'pointlike' rows correspond to Eqs. (10)/(12) and Eqs. (13)/(14)/(16), respectively, to help the reader navigate the various limits.","section":null},{"comment":"Sec. V, discussion of Ref. [9]: The statement that 'Eq. (12) of that reference appears to contain a mistake' is specific and useful. It would help the reader if the corresponding corrected expressions were written out explicitly, perhaps in a short appendix, so that the identification of errors is self-contained.","section":null},{"comment":"The abstract states the µH result 'differs from previous evaluations by roughly ten times the experimental precision anticipated by the upcoming CREMA and FAMU measurements.' This is accurate, but the abstract could also mention that the quoted uncertainty itself is subject to additional systematic inflation, as discussed in Sec. V.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is well-motivated and technically sound. The main physics result—suppression of recoil corrections by nuclear form factors—is well-argued and supported by the limiting cases and numerical checks. The two major comments concern uncertainty presentation rather than correctness of the central derivation. The FF model dependence concern raised in the stress-test note does not undermine the central claim: the ~2.6 µeV deviation from Borie's result vastly exceeds the ~0.086 µeV FF model uncertainty. The identification of errors in Ref. [9] is specific and verifiable. I recommend minor revision to address the uncertainty presentation and the minor points above."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful reading and for the constructive suggestions regarding the transparency of our uncertainty budget. Both major comments are well-taken and will be addressed in the revised manuscript.","responses":[{"response":"The referee is correct. The 4% proton FF model uncertainty is discussed in Sec. V but is explicitly excluded from the numbers in Table II, as stated in the text ('This uncertainty is not included in the tables below'). We agree that this creates a risk of misinterpretation: a reader consulting Table II without reading the surrounding prose would underestimate the full systematic budget by nearly an order of magnitude. In the revised manuscript, we will add a footnote to Table II stating the total uncertainty including the FF model dependence explicitly. For µH, the 4% FF model uncertainty corresponds to approximately 0.086 µeV, which, added in quadrature with the quoted 0.011 µeV, gives a total of approximately 0.087 µeV. For H, the corresponding total is approximately 0.0039 kHz. We will also add the analogous statement for the (µ)H entries. We note that the central conclusion—that our result deviates from previous evaluations by far more than the anticipated experimental precision—remains valid under either uncertainty estimate, as the referee also acknowledges.","revision_made":"yes","referee_comment":"Sec. V, Table II: The quoted uncertainty on the µH HFS result, 0.011 µeV, does not include the 4% proton FF model uncertainty (~0.086 µeV) that the authors themselves identify in Sec. V. The tabulated error thus substantially understates the full systematic budget. While the central claim of significant deviation from previous work remains valid—the ~2.6 µeV discrepancy with Borie's result far exceeds even the inflated uncertainty—the authors should either include the FF model uncertainty in the tabulated error or, at minimum, state the total inflated uncertainty explicitly in the table caption."},{"response":"We agree that the current presentation buries the full uncertainty budget in narrative arithmetic, making it difficult for the reader to reconstruct the total inflated uncertainty for each system. In the revised manuscript, we will add a compact summary—either as a footnote to Table II or as a small additional table—listing, for each system, the quoted (statistical + R-ratio systematic) uncertainty, the scatter inflation factor, the a_µ-discrepancy inflation factor, and the final inflated uncertainty. Concretely, for µH the final inflated uncertainty including both the e⁺e⁻ data scatter and the a_µ discrepancy is approximately 0.011 × 2.51 × 2.42 ≈ 0.067 µeV; adding the 4% FF model uncertainty in quadrature gives approximately 0.11 µeV. For µ³He⁺, the corresponding inflated uncertainty is approximately 0.057 × 1.06 × 1.25 ≈ 0.075 µeV, which is still dominated by the helion FF scatter. We will present these numbers transparently so that the full systematic budget is immediately accessible.","revision_made":"yes","referee_comment":"Sec. V: The treatment of the e⁺e⁻ data scatter and the a_µ discrepancy inflation factor (2.44) is described in prose but not propagated into the tabulated uncertainties. The text states that accounting for experimental scatter gives an extra factor of 2.51 for (µ)H and 1.06 for (µ)³He⁺, followed by an additional inflation factor of 2.42 and 1.25, respectively. It would strengthen the paper to provide a compact summary table or footnote giving the final inflated uncertainty for each system, so that the full systematic budget is transparent and not buried in narrative arithmetic from the prose."}],"tokens_in":20429,"tokens_out":1069,"duration_ms":67144,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The headline: this paper produces new numerical evaluations of hadronic vacuum polarization (hVP) contributions to the hyperfine splitting (HFS) in muonic hydrogen and muonic helium-3, and the results deviate from all prior literature by amounts that matter for upcoming experiments. The µH HFS result of 2.153(11) µeV differs from Borie's 4.8 µeV by roughly 2.6 µeV — about ten times the CREMA/FAMU target precision. The paper also presents the first evaluation of the O(Z⁵α⁶) hVP–finite-size correction to the Lamb shift in µ³He⁺, finding 6.7(2) µeV, which is non-negligible relative to the pointlike hVP uncertainty there. The Lamb shift pointlike results agree with existing literature, which is a useful consistency check. The Mu HFS result (0.2333(11) kHz) also agrees with independent evaluations, confirming the formalism in a regime where recoil matters and finite-size effects are absent. The identification of specific calculational errors in Faustov and Martynenko's earlier work (extra factor of 2 on F₂, missing factors of 1/2 and M/2m in their Eqs. 25 and 26) is concrete and checkable — their master formula reduces to this paper's Eq. (10) once corrected. The physical argument for why recoil corrections are suppressed by nuclear form factors is well-made and supported by Fig. 2: the FF cutoff acts below the nuclear mass scale, so the logarithmic enhancement never develops. The formalism has multiple consistency checks built in — pointlike, structureless, and non-recoil limits all reproduce known results. The soft spot is the uncertainty budget. The quoted 0.011 µeV error on the µH HFS reflects only the R-ratio statistical uncertainty and FF parameter covariance. The paper acknowledges a 4% proton FF model uncertainty (from comparing Borah et al. vs. Lin et al. dispersion-theoretical FFs) but does not fold it into the tabulated errors. The e⁺e⁻ data scatter inflation factor (~2.5 for µH) is discussed but also not included in the tables. Even combining all these, the total inflated uncertainty is ~0.11 µeV, still a factor of ~25 below the 2.6 µeV deviation from prior work. So this is a presentation issue, not a threat to the central claim. The reader and stress-test note both correctly identify this as the weakest assumption and correctly conclude it is not load-bearing. This paper is for precision atomic spectroscopy theorists and experimentalists interpreting CREMA/FAMU data. It deserves a serious referee — the formalism is clean, the new results are consequential, and the error identifications in prior work need independent verification. I'd ask the referee to check whether the 4% FF model uncertainty and the data-scatter inflation should be folded into the quoted errors before publication, but the physics is sound.","headline":"New hVP HFS evaluations for muonic hydrogen and helium-3 deviate significantly from prior work; errors identified in earlier calculations.","tokens_in":21601,"tokens_out":708,"would_cite":true,"duration_ms":142259,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.30.jf","32.10.Fn","36.10.Dr","13.40.Gp"],"model":"glm-5.2","headline":"Nuclear form factors suppress recoil corrections in hadronic vacuum polarization of light atoms","keywords":["hadronic vacuum polarization","hyperfine splitting","muonic hydrogen","muonic helium-3","nuclear form factors","recoil corrections","finite-size effects","Lamb shift"],"falsifier":"If the CREMA or FAMU measurement of the muonic hydrogen ground-state HFS lands at a value inconsistent with the theory prediction that uses 2.153(11) μeV for the hVP piece, the discrepancy would point either to a problem in the hVP evaluation or to missing physics elsewhere in the theory compilation.","tokens_in":20598,"feed_emoji":"⚛️","tokens_out":1392,"duration_ms":78829,"temperature":0.7,"pith_summary":"This paper evaluates the hadronic vacuum polarization (hVP) contribution to the Lamb shift and hyperfine splitting in ordinary and muonic hydrogen and helium-3 ions. The central finding is a mechanism: in systems with composite nuclei, the nuclear elastic form factors act as a soft cutoff on the loop integrals, suppressing the recoil corrections that would otherwise be logarithmically enhanced by the lepton-to-nucleus mass ratio. In muonium, where both constituents are pointlike, recoil corrections reduce the hVP contribution by more than a factor of five; in muonic hydrogen and helium, the form factors cut off the integral well below the nuclear mass scale, so the recoil correction becomes negligible. Using this insight and a data-driven dispersive evaluation of the R ratio, the authors obtain hVP contributions to the ground-state hyperfine splitting of 2.153(11) μeV in muonic hydrogen and −15.19(57) μeV in muonic helium-3. These results deviate from all prior evaluations, with the muonic hydrogen value differing by roughly ten times the precision targeted by upcoming CREMA and FAMU experiments. The paper attributes the discrepancies to two sources: a rescaling of muonic vacuum polarization that misstates the hVP-to-μVP ratio in the HFS context, and calculational errors in an earlier combined hVP–finite-size treatment. For the Lamb shift, the results agree with existing literature, and the paper provides a first evaluation of a subleading hVP–finite-size correction that is non-negligible in muonic helium-3.","feed_headline":"Hadronic vacuum polarization in muonic hydrogen rewritten by form-factor cutoff","feed_subtitle":"Nuclear elastic form factors suppress recoil corrections that prior work overcounted, shifting the hVP contribution to the hyperfine split 2","key_machinery":"hadronic vacuum polarization (hVP): the non-perturbative QCD contribution to the photon self-energy, evaluated here via a dispersive integral over the empirical R ratio (the e+e− annihilation cross section into hadrons). The elastic electromagnetic form factors (Sachs G_E and G_M) of the proton and helion encode the nuclear finite-size effects and provide the cutoff mechanism.","core_discovery":"The paper identifies a quantitative criterion for when recoil corrections to a vacuum-polarization contribution matter: they are logarithmically enhanced only if the VP spectral function extends to scales comparable with the heavier constituent mass and no form-factor cutoff intervenes at a lower scale. For composite nuclei, the elastic form factors always provide the lower cutoff, so the recoil logarithm never develops. This mechanism explains why the full and non-recoil weighting functions are nearly indistinguishable in muonic hydrogen and helium, while in pointlike muonium they differ by more than a factor of five. Applying this with a dispersive data-driven evaluation yields hVP–HFS数值值在","pith_inferences":[],"forward_implications":["Upcoming CREMA and FAMU measurements of the muonic hydrogen ground-state hyperfine splitting at 1 ppm precision will directly test the paper's revised hVP value of 2.153(11) μeV against prior estimates that differ by ~2 μeV.","The recoil-suppression criterion can be applied to other hydrogen-like systems with composite nuclei—including muonic deuterium and heavier muonic ions—to decide a priori whether recoil corrections to any VP contribution can be neglected.","The first evaluation of the O(Z⁵α⁶) hVP–finite-size correction in muonic helium-3 (6.7(2) μeV) sets a benchmark for theory compilations of n=2 levels in that system.","Correlations between hVP contributions across different observables (Lamb shift, HFS, muon g−2) arise from shared R-ratio input, and the paper's weighting-function analysis shows these observables probe similar spectral regions, opening paths for joint consistency checks."],"fun_headline_variants":["Form factors suppress recoil corrections in muonic hydrogen HFS","Hadronic vacuum polarization in muonic hydrogen bounded by nuclear size","Finite nuclear size suppresses recoil hVP corrections in muonic atoms","Recoil corrections to vacuum polarization curbed by nuclear form factors","Form factor cutoff shifts hVP hyperfine splitting in muonic hydrogen"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The proton elastic form factors used for the finite-size calculation are taken from a specific fit that imposes the muonic-hydrogen Lamb-shift charge radius; an alternative dispersion-theoretical parametrization yields a 3.7% larger HFS correction, and the paper assigns a 4% model uncertainty that it does not fold into its tabulated errors.","fun_headline_variants_meta":{"raw":{"variants":["Form factors suppress recoil corrections in muonic hydrogen HFS","Hadronic vacuum polarization in muonic hydrogen bounded by nuclear size","Finite nuclear size suppresses recoil hVP corrections in muonic atoms","Recoil corrections to vacuum polarization curbed by nuclear form factors","Form factor cutoff shifts hVP hyperfine splitting in muonic hydrogen"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1273,"prompt_tokens":677,"completion_tokens":596,"prompt_tokens_details":null},"tokens_in":677,"tokens_out":596,"duration_ms":32394,"temperature":1.0,"reasoning_tokens":587,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T03:28:38.138797+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the CREMA or FAMU measurement of the muonic hydrogen ground-state HFS lands at a value inconsistent with the theory prediction that uses 2.153(11) μeV for the hVP piece, the discrepancy would point either to a problem in the hVP evaluation or to missing physics elsewhere in the theory compilation.","supporting_citations":[],"review_version":1}