{"id":"ce797bee-e3cc-436c-8dc3-277800cfd22e","arxiv_id":"2607.05917","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"Exact form-factor convolution of the Yukawa potential yields model-independent Lamb-shift bounds on α(λ) for λ in [10^{-4}, 10^6] fm, with a 2s–2p cancellation singularity near 10^5 fm.","lead":"An exact momentum-space method turns the six-dimensional Yukawa energy-shift integral into a one-dimensional convolution with atomic form factors, then uses hydrogen and deuterium Lamb-shift data to bound the fifth-force coupling. The bounds reveal a cancellation near 10^5 fm where the 2s–2p shift vanishes, and show that present spectroscopic precision is still too coarse to see such forces.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"Central bounds rest on treating the Lamb-shift residual as pure first-order Yukawa; short-range degeneracy with nuclear-radius uncertainty can invalidate the claimed model-independent α(λ) limits.","rationale":"The Reader correctly identified that only the abstract was available and that the residual-to-bound step is the weakest assumption; LOW confidence and UNVERDICTED are appropriate. I sharpen the same concern: finite-size degeneracy at short λ is the concrete place the assumption fails most clearly, and the title/abstract conflict (muonic vs electronic) blocks even a consistency check of the quoted resonance at λ ∼ 10^5 fm (which matches the electronic Bohr radius, not the muonic one). The momentum-space reduction itself is standard and not in doubt; the 2s–2p cancellation is physically expected. No stronger verdict is possible until the full manuscript’s uncertainty budget and the identity of the atomic system are inspected. The Reader’s UNVERDICTED therefore stands; the concrete test is the minimal calculation that would confirm or refute model independence once the text is available.","tokens_in":2082,"tokens_out":628,"duration_ms":78857,"concrete_test":"Obtain the full text and extract the residual/uncertainty ΔE used to set |δE_Yukawa(α,λ)| < ΔE at λ = 1 fm and at λ ∼ 10^5 fm. Re-evaluate the α bound at λ = 1 fm after promoting the rms charge radius to a free nuisance parameter with its experimental prior and profiling it out; if the resulting |α| limit weakens by more than a factor of ∼3, the short-range model-independence claim fails and the published constraints are overstated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim converts non-observation of an anomalous 2s–2p residual into upper limits on α(λ) via an exact momentum-space formula folding the atomic form factor with the Yukawa potential. That conversion requires 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; any portion of the residual currently absorbed into the fitted proton or deuteron radius is therefore double-counted or unavailable for new-physics bounds. The abstract asserts model independence and that induced frequency shifts lie well below spectroscopic resolution, yet does not indicate whether the nuclear radius is floated as a nuisance when α is constrained. Without that partition the short-distance exclusion can be artificially tight. A secondary obstacle is the title–abstract mismatch (muonic vs ordinary H/D), which changes the Bohr radius, the 2s–2p cancellation location, and the applicable uncertainty budget, blocking verification from the abstract alone.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","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.","tokens_in":2311,"tokens_out":909,"duration_ms":16610,"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":[{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"minor_comments":[{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"Figures showing α(λ) should indicate the location of the 2s–2p cancellation and the transition between attractive and repulsive regimes more prominently.","section":null},{"comment":"A short comparison with earlier approximate treatments of the same problem would help the reader gauge the improvement obtained by the exact convolution.","section":null}],"recommendation":"major_revision","confidential_remarks":"The title–abstract mismatch is sufficiently glaring that it raises a question of whether the manuscript was hastily retitled or whether two related analyses were conflated. I would ask the authors to clarify this point before the paper is sent to further referees. The short-distance degeneracy with the nuclear radius is a standard issue in this literature; if the authors have already addressed it carefully in the full text, the major comment can be downgraded after revision."},"author_rebuttal":{"model":"grok-4.5","summary":"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.","responses":[{"response":"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_made":"yes","referee_comment":"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."},{"response":"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_made":"yes","referee_comment":"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."},{"response":"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_made":"yes","referee_comment":"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."}],"tokens_in":1869,"tokens_out":1014,"duration_ms":28463,"standing_objections":[]},"desk_editor":{"model":"grok-4.5","letter":"The one thing worth knowing is the method: they rewrite the sixfold coordinate-space energy shift as an exact one-dimensional momentum-space convolution of the atomic form factor with the Yukawa potential. That is cleaner than the usual approximate treatments, and it lets them fold in empirical nuclear form factors directly. The second thing is the cancellation: at λ ∼ 10^5 fm the 2s and 2p shifts cancel, the energy difference vanishes, and the bound on α blows up. That singularity is a real feature of the 2s–2p difference and people quoting fifth-force limits from Lamb shifts should know it is there.\n\nWhat the paper does well is the reduction itself and the claim that the resulting α(λ) curves are independent of the nuclear-charge model once the form factor is fixed. Using published H and D spectroscopy data over six orders of magnitude in range is the right application. The statement that the induced frequency shifts sit well below current resolution is also honest; it correctly identifies spectroscopic precision, not nuclear modeling, as the bottleneck for most of the range.\n\nThe soft spots are real but localized. First, the title says muonic hydrogen and deuterium while the abstract and the stated data are ordinary H and D. That is not cosmetic: the Bohr radius, the cancellation location, and the uncertainty budget all change. Second, and more load-bearing, at λ ≲ a few fm a Yukawa is nearly degenerate with a shift of the rms charge radius already present in the Standard-Model finite-size correction. If the empirical form factor is frozen and the entire residual (or its uncertainty) is attributed to a first-order Yukawa, the short-distance exclusion can look tighter than it is. The abstract asserts model independence and does not say whether the radius is floated as a nuisance when α is constrained. That partition needs to be shown explicitly; without it the short-range claim is overstated. The long-range and intermediate-range results, and the cancellation window itself, do not suffer from this.\n\nThis is for people who extract light-mediator limits from precision atomic spectroscopy. A referee who can check the integral transform and demand a clear radius-versus-α treatment at short distance will get value from it. The central argument for the method and the cancellation holds up; the short-range bounds need tightening, not discarding.\n\nSend it to peer review. It deserves a serious referee.","headline":"Useful exact form-factor convolution and a real 2s–2p cancellation window; short-range α bounds are the soft spot because of radius degeneracy.","tokens_in":2938,"tokens_out":591,"would_cite":false,"duration_ms":28325,"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":"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.","keywords":["Yukawa fifth force","Lamb shift","atomic form factor","momentum-space convolution","hydrogen spectroscopy","deuterium spectroscopy","new-physics constraints","fifth-force coupling"],"falsifier":"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.","tokens_in":2966,"feed_emoji":"⚛️","tokens_out":958,"duration_ms":41658,"temperature":0.7,"pith_summary":"The paper shows that the energy shift induced by a Yukawa-type fifth force in the Lamb transition can be rewritten exactly as a one-dimensional momentum-space convolution of the atomic form factor with the Fourier transform of the Yukawa potential. This eliminates the six-fold coordinate-space integral and lets measured nuclear charge distributions enter the calculation directly, without intermediate model assumptions. When the formula is evaluated on the latest hydrogen and deuterium spectroscopic residuals, it produces upper limits on the coupling α for interaction ranges λ spanning 10^{-4} to 10^6 fm. At a critical range near 10^5 fm the 2s and 2p orbital shifts cancel, the net energy difference vanishes, and the bound on α diverges, separating regimes of attractive and repulsive effective forces. The frequency shifts that would be produced by any still-allowed fifth force lie well below present experimental resolution, so spectroscopic precision itself is the limiting factor for further progress.","feed_headline":"Lamb-shift residuals bound Yukawa forces across six decades of range","feed_subtitle":"Exact form-factor convolution yields model-free limits; a 10^5 fm cancellation opens a blind spot","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Exact form-factor convolution bounds Yukawa forces via Lamb residuals","Model-free Yukawa limits from muonic hydrogen Lamb-shift data","Lamb-shift cancellation at 10^5 fm erases net Yukawa energy shift","Atomic form factors yield nuclear-model-independent Yukawa bounds","Spectroscopy precision bottlenecks Yukawa probes in Lamb residuals"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Exact form-factor convolution bounds Yukawa forces via Lamb residuals","Model-free Yukawa limits from muonic hydrogen Lamb-shift data","Lamb-shift cancellation at 10^5 fm erases net Yukawa energy shift","Atomic form factors yield nuclear-model-independent Yukawa bounds","Spectroscopy precision bottlenecks Yukawa probes in Lamb residuals"]},"model":"grok-4.5","cost_usd":0.018846,"raw_usage":{"total_tokens":3743,"prompt_tokens":807,"num_sources_used":0,"completion_tokens":76,"cost_in_usd_ticks":188460000,"prompt_tokens_details":{"text_tokens":807,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2860,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":807,"tokens_out":76,"duration_ms":155627,"temperature":1.0,"reasoning_tokens":2860,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T19:29:25.347143+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"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.","supporting_citations":[],"review_version":1}