REVIEW 3 major objections 4 minor 15 references
Reviewed by Pith at T0; open to challenge.
T0 means a machine referee read the full paper against a public rubric. The mark states how deep the mechanical check went, never who wrote it. the ladder, T0–T4 →
T0 review · grok-4.5
An exact one-dimensional convolution of atomic form factors with the Yukawa potential converts hydrogen and deuterium Lamb-shift residuals into model-independent bounds on fifth-force coupling strength across six decades of range.
2026-07-08 19:29 UTC pith:K4RP5C6N
load-bearing objection Useful exact form-factor convolution and a real 2s–2p cancellation window; short-range α bounds are the soft spot because of radius degeneracy. the 3 major comments →
Constraints on Yukawa-type New Forces from the Lamb Shift in Muonic Hydrogen and Deuterium
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The first-order Lamb-shift energy difference caused by a Yukawa potential reduces exactly to a one-dimensional integral over the product of the atomic form factor and the Yukawa propagator. Inserting empirical proton and deuteron form factors and comparing with residual discrepancies in measured 2s–2p intervals yields model-independent upper bounds on the fifth-force coupling α over λ ∈ [10^{-4}, 10^6] fm, with a resonance-like cancellation at λ ∼ 10^5 fm that drives the net shift to zero and makes the constraint diverge.
What carries the argument
The exact one-dimensional momentum-space convolution of the atomic form factor with the Fourier-transformed Yukawa potential, which replaces the original six-fold spatial integral while retaining the full empirical nuclear charge distribution.
Load-bearing premise
Any residual mismatch between measured Lamb shifts and Standard-Model theory can be attributed solely to a first-order Yukawa perturbation whose nuclear dependence is completely captured by the empirical form factor.
What would settle it
A new high-precision determination of the hydrogen or deuterium 2s–2p interval whose residual (after all known QED and nuclear corrections) exceeds the energy shift given by the convolution formula for any previously allowed (α, λ) pair would falsify the derived bound.
If this is right
- Upper limits on α are obtained across λ from 10^{-4} to 10^6 fm that do not depend on any particular model of the nuclear charge distribution.
- At λ ≈ 10^5 fm the 2s–2p signature of the fifth force vanishes, so Lamb-shift spectroscopy loses all sensitivity at that special range regardless of precision.
- Induced transition-frequency shifts for still-allowed couplings remain far below current spectroscopic resolution, identifying experimental precision as the primary bottleneck.
- The same convolution expression can be re-evaluated with future higher-precision data to tighten the bounds without reformulating the nuclear physics.
Where Pith is reading between the lines
- Spectroscopy of additional fine-structure intervals (higher n or different principal shells) could lift the 2s–2p cancellation and restore sensitivity near the blind-spot range λ ∼ 10^5 fm.
- The identical form-factor convolution can be applied at once to muonic atoms or light ions once their residuals and empirical form factors are available, extending the same model-independent reach.
- Even if the present Lamb-shift residuals are later reassigned to some non-Yukawa source, the quoted curves remain conservative upper limits on any additional Yukawa component of that strength.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops an analytical framework that converts the six-dimensional coordinate-space integral for a Yukawa-type energy shift of the Lamb interval into an exact one-dimensional momentum-space convolution of the atomic form factor with the Yukawa potential. Empirical nuclear form factors are inserted directly, and the resulting expression is applied to published hydrogen and deuterium Lamb-shift data to extract upper limits on the fifth-force coupling α over the range λ ∈ [10^{-4}, 10^6] fm. A resonance-like cancellation between the 2s and 2p orbital shifts is identified near λ ∼ 10^5 fm, where the energy difference vanishes and the bound on α diverges. The authors claim that the constraints are independent of the nuclear-charge model and that the induced frequency shifts lie well below present spectroscopic resolution.
Significance. If the derivation and the subsequent comparison with data are free of load-bearing errors, the work supplies a clean, essentially parameter-free method for translating Lamb-shift residuals into fifth-force bounds over many orders of magnitude in range. The exact reduction of the six-fold integral and the explicit demonstration of model independence with respect to the nuclear form factor are genuine technical advances over earlier approximate treatments. The identification of the 2s–2p cancellation and the clear statement that spectroscopic precision is the limiting factor are useful for the community. The result is therefore of moderate but real interest to precision atomic physics and light-force searches.
major comments (3)
- The conversion of a non-observation of a residual into an upper bound on α presupposes that every residual (or the full experimental/theoretical uncertainty) may be attributed to a first-order Yukawa shift while the nuclear charge distribution is held fixed at its empirical form factor. At λ ≲ few fm the Yukawa interaction is nearly degenerate with a shift of the rms charge radius already present in Standard-Model finite-size corrections. The manuscript does not demonstrate that the nuclear radius is floated as a nuisance parameter when α is constrained, nor does it quantify the fraction of the residual currently absorbed into the fitted proton/deuteron radius. Without that partition the short-distance exclusion can be artificially tight and the claim of model independence is incomplete.
- The title refers to “Muonic Hydrogen and Deuterium” while the abstract and the body of the paper treat ordinary (electronic) hydrogen and deuterium. The Bohr radius, the location of the 2s–2p cancellation, the relevant uncertainty budget, and the numerical values of the bounds all change by orders of magnitude between the two systems. This mismatch must be resolved; either the title is corrected or the analysis is extended to the muonic case with a clear statement of which data set is used.
- The abstract asserts that “the constraints are shown to be independent of the nuclear charge distribution model.” The body must contain an explicit numerical comparison (different form-factor parameterizations, or a pure point-like nucleus versus a realistic form factor) that demonstrates the residual dependence of α(λ) is negligible across the quoted range. If that comparison is only qualitative, the claim should be weakened.
minor comments (4)
- Notation for the Yukawa coupling α should be clearly distinguished from the fine-structure constant; a different symbol (e.g., α_Y or α_5) would avoid confusion.
- The precise numerical values of the experimental and theoretical Lamb-shift residuals adopted for H and D should be tabulated, together with their uncertainties and the references from which they are taken.
- Figures showing α(λ) should indicate the location of the 2s–2p cancellation and the transition between attractive and repulsive regimes more prominently.
- A short comparison with earlier approximate treatments of the same problem would help the reader gauge the improvement obtained by the exact convolution.
Simulated Author's Rebuttal
We thank the referee for a careful and constructive report. The three major comments identify genuine issues of interpretation (short-range radius degeneracy), presentation (title–body mismatch), and documentation of a central claim (model independence). We address each point below and will revise the manuscript accordingly. None of the comments invalidates the exact convolution framework or the intermediate- and long-range bounds; they do require clearer statements of assumptions and, in one case, a corrected title.
read point-by-point responses
-
Referee: The conversion of a non-observation of a residual into an upper bound on α presupposes that every residual may be attributed to a first-order Yukawa shift while the nuclear charge distribution is held fixed. At λ ≲ few fm the Yukawa interaction is nearly degenerate with a shift of the rms charge radius. The manuscript does not float the nuclear radius as a nuisance parameter nor quantify the fraction of residual absorbed into the fitted radius; without that partition the short-distance exclusion can be artificially tight and model independence is incomplete.
Authors: We agree that this is a load-bearing limitation of the present analysis and that the manuscript does not adequately address it. Our bounds are obtained by attributing the full experimental/theoretical residual (with the nuclear form factor held fixed at its empirical value) to a first-order Yukawa shift. At λ ≲ few fm the Yukawa potential is nearly degenerate with a change in the rms charge radius already present in Standard-Model finite-size corrections; floating the radius as a nuisance parameter would therefore absorb part of any short-range signal and loosen the bound. We will revise the text to (i) state explicitly that the nuclear radius/form factor is held fixed, (ii) discuss the degeneracy with the rms radius for λ ≲ few fm, (iii) mark the short-distance portion of the exclusion as conditional on that assumption, and (iv) note that a joint fit of α and the radius lies beyond the scope of this work. The intermediate- and long-range constraints (λ ≳ tens of fm), where the Yukawa shape is distinguishable from a pure radius shift, remain unaffected. We will also soften the abstract claim of model independence in the short-distance regime. revision: yes
-
Referee: The title refers to “Muonic Hydrogen and Deuterium” while the abstract and the body of the paper treat ordinary (electronic) hydrogen and deuterium. The Bohr radius, the location of the 2s–2p cancellation, the relevant uncertainty budget, and the numerical values of the bounds all change by orders of magnitude between the two systems. This mismatch must be resolved.
Authors: The referee is correct. The analysis, data sets, Bohr-radius scale, and the 2s–2p cancellation near λ ∼ 10^5 fm all refer to ordinary (electronic) hydrogen and deuterium. The word “Muonic” in the title is an error and will be removed. The revised title will read “Constraints on Yukawa-type New Forces from the Lamb Shift in Hydrogen and Deuterium,” matching the abstract and body. We will also check the manuscript for any residual wording that could suggest a muonic analysis and correct it. Extending the same framework to muonic hydrogen and deuterium is of interest but is outside the present work; we will not claim muonic results. revision: yes
-
Referee: The abstract asserts that “the constraints are shown to be independent of the nuclear charge distribution model.” The body must contain an explicit numerical comparison (different form-factor parameterizations, or a pure point-like nucleus versus a realistic form factor) that demonstrates the residual dependence of α(λ) is negligible across the quoted range. If that comparison is only qualitative, the claim should be weakened.
Authors: We accept the criticism. The present manuscript argues model independence largely on the structure of the momentum-space convolution (the atomic form factor multiplies the nuclear form factor, and empirical form factors are inserted directly), but it does not display a quantitative side-by-side comparison of α(λ) for different nuclear models. In the revision we will add an explicit numerical comparison—e.g., dipole versus Gaussian versus a point-like nucleus versus the empirical form factors used in the main results—showing the residual variation of the bound as a function of λ. Where the variation is negligible we will retain a carefully worded independence statement; where it is not (principally at the shortest ranges, consistent with Comment 1) we will weaken the claim. The abstract sentence will be revised to match the body after that comparison is in place. revision: yes
Circularity Check
No significant circularity: external form factors and spectroscopy data constrain α(λ) via an exact first-order energy-shift formula; nothing is forced by definition or self-fit.
full rationale
The claimed chain is a standard external-data constraint analysis. Empirical nuclear form factors and published H/D Lamb-shift spectroscopy enter as independent inputs; the six-fold coordinate-space energy shift is rewritten as a one-dimensional momentum-space convolution of the atomic form factor with the Yukawa potential (a technical identity, not a fit); the residual or uncertainty relative to Standard-Model QED is then converted into upper bounds on the free parameters α and λ. Neither α nor λ is adjusted to a subset of the same data and then re-presented as a prediction, nor is any uniqueness theorem or ansatz imported from the authors’ prior work to force the result. The resonance-like 2s–2p cancellation at λ ∼ 10^5 fm is a derived feature of the same first-order formula, not an input. Model-independence of the nuclear charge distribution is asserted after the form-factor convolution, not assumed by renaming. Potential systematic concerns (short-range degeneracy with nuclear-radius uncertainty, possible prior absorption of fifth-force effects into the spectroscopic residual, title–abstract muonic/ordinary mismatch) affect correctness or interpretation of the bounds, not circularity of the derivation. No self-definitional loop, fitted-input-as-prediction, load-bearing self-citation, or renaming of a known empirical pattern is exhibited by the abstract or the stated chain. Score 0 is therefore the honest finding.
Axiom & Free-Parameter Ledger
axioms (3)
- domain assumption The new interaction is a static Yukawa potential V(r) ∝ α exp(−r/λ)/r that can be treated in first-order perturbation theory on the atomic states.
- domain assumption Nuclear finite-size effects factorize completely into the measured nuclear charge form factor, so the sixfold integral reduces exactly to a one-dimensional momentum-space convolution.
- domain assumption Standard-Model QED Lamb-shift theory (including all known higher-order corrections) is accurate enough that any residual can be interpreted as a bound on a new Yukawa contribution.
read the original abstract
We develop an analytical framework to constrain Yukawa-type fifth-force interactions using the Lamb shift in hydrogen and deuterium. The sixfold coordinate-space integral for the energy shift is transformed into a one-dimensional momentum-space convolution of the atomic form factor and the Yukawa potential, eliminating the need for approximations. This formulation directly incorporates empirical nuclear charge distributions through their form factors, making the implementation of experimental results straightforward. Applying this method to the latest hydrogen and deuterium spectroscopy data, we derive constraints on the coupling constant $\alpha$ in the interaction range $\lambda \in [10^{-4}, 10^{6}]$~fm. A resonance-like cancellation between the $2s$ and $2p$ orbital shifts is identified at a critical range $\lambda \sim 10^{5}$~fm, where the energy shift vanishes and the constraint on $\alpha$ diverges, marking a transition between attractive and repulsive regimes. The constraints are shown to be independent of the nuclear charge distribution model. Our results indicate that the transition frequency shifts induced by such new forces lie well below current spectroscopic resolution, establishing spectroscopic precision as the primary bottleneck for probing these interactions.
Figures
Reference graph
Works this paper leans on
-
[1]
Chung-Chien Huang and Li-Bang Wang. Constraints on extra-dimension gravity by muonic hydrogen and deuterium spectroscopy.Modern Physics Letters A, 0(0):2650135, 0
-
[2]
Proton radius puzzle and large extra dimensions.Modern Physics Letters A, 28(20):1350094, 2013
LI-BANG WANG and WEI-TOU NI. Proton radius puzzle and large extra dimensions.Modern Physics Letters A, 28(20):1350094, 2013
work page 2013
-
[3]
Bounds to strong gravity from electron-positron annihilations near thresh- old.Phys
Nathaniel Alden and Roberto Onofrio. Bounds to strong gravity from electron-positron annihilations near thresh- old.Phys. Rev. D, 109:115021, Jun 2024
work page 2024
- [4]
-
[5]
Zachary Bogorad, Peter W. Graham, and Giorgio Gratta. Detecting nanometer-scale new forces with coherent neu- tron scattering.Phys. Rev. D, 108:055005, Sep 2023
work page 2023
-
[6]
Alfio M. Bonanno and Emiliano M. Glaviano. Fifth-force constraints from uv-complete scalar-tensor gravity.Phys. Rev. Lett., 136:201501, May 2026
work page 2026
-
[7]
E.J. Salumbides, W. Ubachs, and V.I. Korobov. Bounds on fifth forces at the sub-˚ alength scale.Journal of Molec- ular Spectroscopy, 300:65–69, 2014. Spectroscopic Tests of Fundamental Physics
work page 2014
-
[8]
Damiano F. G. Fiorillo, Alessandro Lella, Ciaran A. J. O’Hare, and Edoardo Vitagliano. Leading bounds on mi- crometer to picometer fifth forces from neutron star cool- ing.Phys. Rev. Lett., 135:211003, Nov 2025. Chung-Chien Huang, Li-Bang Wang: Constraints on Yukawa-type New Forces 9
work page 2025
-
[9]
ROBERTO ONOFRIO. Muonic hydrogen as a quantum gravimeter.International Journal of Modern Physics D, 23(01):1450005, 2014
work page 2014
-
[10]
Proton radius puzzle and quantum grav- ity at the fermi scale.EPL, 104(2):20002, 2013
Onofrio, Roberto. Proton radius puzzle and quantum grav- ity at the fermi scale.EPL, 104(2):20002, 2013
work page 2013
-
[11]
N. Bezginov, T. Valdez, M. Horbatsch, A. Marsman, A. C. Vutha, and E. A. Hessels. A measurement of the atomic hydrogen lamb shift and the proton charge radius.Science, 365(6457):1007–1012, 2019
work page 2019
-
[12]
Amaro, Fran¸ cois Biraben, Jo˜ ao M
Aldo Antognini, Fran¸ cois Nez, Karsten Schuhmann, Fer- nando D. Amaro, Fran¸ cois Biraben, Jo˜ ao M. R. Cardoso, Daniel S. Covita, Andreas Dax, Satish Dhawan, Marc Diepold, Luis M. P. Fernandes, Adolf Giesen, Andrea L. Gouvea, Thomas Graf, Theodor W. H¨ ansch, Paul Indeli- cato, Lucile Julien, Cheng-Yang Kao, Paul Knowles, Franz Kottmann, Eric-Olivier Le...
work page 2013
-
[13]
Amaro, Fran¸ cois Biraben, Jo˜ ao M
Randolf Pohl, Aldo Antognini, Fran¸ cois Nez, Fernando D. Amaro, Fran¸ cois Biraben, Jo˜ ao M. R. Cardoso, Daniel S. Covita, Andreas Dax, Satish Dhawan, Luis M. P. Fernan- des, Adolf Giesen, Thomas Graf, Theodor W. H¨ ansch, Paul Indelicato, Lucile Julien, Cheng-Yang Kao, Paul Knowles, Eric-Olivier Le Bigot, Yi-Wei Liu, Jos´ eA. M. Lopes, Livia Ludhova, C...
work page 2010
-
[14]
Deuteron charge radius and ryd- berg constant from spectroscopy data in atomic deuterium
Randolf Pohl, Fran¸ cois Nez, Thomas Udem, Aldo An- tognini, Axel Beyer, H´ el` ene Fleurbaey, Alexey Grinin, Theodor W H¨ ansch, Lucile Julien, Franz Kottmann, Ju- lian J Krauth, Lothar Maisenbacher, Arthur Matveev, and Fran¸ cois Biraben. Deuteron charge radius and ryd- berg constant from spectroscopy data in atomic deuterium. Metrologia, 54(2):L1, mar 2017
work page 2017
-
[15]
Randolf Pohl, Fran¸ cois Nez, Luis M. P. Fernandes, Fer- nando D. Amaro, Fran¸ cois Biraben, Jo˜ ao M. R. Cardoso, Daniel S. Covita, Andreas Dax, Satish Dhawan, Marc Diepold, Adolf Giesen, Andrea L. Gouvea, Thomas Graf, Theodor W. H¨ ansch, Paul Indelicato, Lucile Julien, Paul Knowles, Franz Kottmann, Eric-Olivier Le Bigot, Yi-Wei Liu, Jos´ e A. M. Lopes,...
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.