{"id":"37a53ae6-c4d3-4c07-bc44-ef8f6ee72f0a","arxiv_id":"2608.12747","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"If the helium ionization anomaly comes from a new scalar boson, positronium 1S-2S transitions should shift by 0.25 to 0.85 MHz, providing an independent test.","lead":"This paper calculates how a hypothetical scalar force, proposed to explain a 9-sigma discrepancy in helium ionization energies, would shift positronium spectral lines. It predicts shifts of 0.25 to 0.85 MHz in a transition that current experiments could eventually measure.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Predicted Ps shifts inherit any error in the numerical helium matrix element C_He(λ), which is used with no stated uncertainty or benchmark; this is the main load-bearing gap.","rationale":"The analytic core of the paper is internally consistent: the hydrogenic integrals in Eqs. (4)–(6) are dimensionally correct and the elimination of g_e via Eq. (9) is sound. I found no fatal internal inconsistency, and the paper is appropriately transparent that the whole program is conditional on the scalar-boson interpretation of the helium anomaly. The genuinely load-bearing weakness is the unquantified numerical helium matrix element C_He(λ). Because the predicted positronium shifts are inversely proportional to it, an error in C_He translates directly into every headline number, including the mass dependence in Fig. 2. The paper gives no uncertainty, no convergence data, and no independent validation for this operator, despite describing the calculation as 'benchmarked' in the abstract. This is the same weakest assumption identified by the reader, and I agree with it. A single independent high-precision recomputation of C_He(λ) at representative mediator masses would settle whether the concern is substantive. If the independent calculation agrees to a few percent, the paper's predictions are robust and the conditional acceptance is justified; if it does not, the error budget of the central claim is incomplete.","tokens_in":8944,"tokens_out":17152,"duration_ms":168429,"concrete_test":"Compute C_He(λ) for mediator masses 0, 100, 400, and 800 eV using an independent high-precision Hylleraas wavefunction for the He 2^3S state (e.g., a large Drake-type expansion with tabulated expectation values), and compare against the values obtained from Ref [24]. Propagate any fractional difference through Eq. (9) for the 1^3S1→2^3S1 and 2S-ionization shifts. If the difference exceeds roughly 10% at any λ, the quoted 0.250–0.850 MHz and 0.14–0.21 MHz windows are not yet reliable; if it is a few percent or less, the concern is resolved and the conditional prediction can stand.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central conversion in Eq. (9), Δν_Ps = (ΔE_He/h) F_nl(λ)/C_He(λ), is only as reliable as the numerical matrix element C_He(λ) = ⟨e^{-r12/λ}/r12⟩_{2^3S} from Eq. (8). The paper states that 'C_He(λ) is evaluated numerically using the He wave function of Ref [24]' and later calls the He state 'approximate', but gives no uncertainty estimate, no convergence study, and no independent check for the Yukawa-weighted 1/r12 operator. The abstract claims a 'benchmarked correlated calculation', yet no benchmark data appear in the manuscript. Because C_He appears in the denominator of Eq. (9), any fractional error in C_He maps one-to-one onto every quoted Ps frequency shift, and any error in its λ-dependence distorts the shape of the curves in Fig. 2. The paper propagates only the ΔE_He uncertainty, so the cited ranges 0.250–0.850 MHz and 0.14–0.21 MHz have no error budget from the wavefunction. This is a genuine correctness risk rather than a disagreement with external consensus: if the Ref [24] wavefunction is accurate for this operator, the concern is resolved; if it is not, the headline predictions shift correspondingly.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper takes the recently reported discrepancy between measured and calculated 2^3S_1 ionization energies of metastable helium as motivation for the scalar-boson interpretation proposed in Ref. [6]. Starting from a finite-range Yukawa potential (Eq. (1)), the authors compute first-order energy shifts for positronium n=1 and n=2 states, using analytic hydrogenic matrix elements (Eqs. (4)-(6)) together with a numerical helium matrix element C_He(λ) (Eq. (8)). Combining the helium anomaly with the ratio F_nl/C_He (Eq. (9)), they predict a 1^3S_1→2^3S_1 shift of 0.250–0.850 MHz and a 2^3S_1 ionization shift of 0.14–0.21 MHz for mediator masses below 800 eV. The paper then assesses the experimental feasibility of testing these shifts, including 1S–2S two-photon spectroscopy, Rydberg-ionization measurements, and n=2 fine-structure transitions, concluding that the 1S–2S interval is the most promising probe.","tokens_in":9194,"tokens_out":11954,"duration_ms":113556,"significance":"If the calculation is correct, the paper provides a concrete, falsifiable bridge between a proposed new-physics explanation of the helium anomaly and an independent leptonic system. The analytic positronium matrix elements are simple to verify, the conversion in Eq. (9) does not involve fitting to positronium data, and the experimental discussion is realistic and well grounded in the literature. The main weakness is that the numerical helium matrix element C_He(λ) is presented without a quantitative uncertainty estimate or benchmark, despite the abstract's claim of a 'benchmarked correlated calculation.' Because C_He appears in the denominator of Eq. (9), its uncertainty maps directly onto all quoted shifts. With that gap filled, the paper would be a valuable contribution to the search for new leptonic forces; in its present form, the central numerical input is not yet adequately supported.","major_comments":[{"comment":"The numerical matrix element C_He(λ) is the only input connecting the helium anomaly to the predicted positronium shifts, but the manuscript provides no uncertainty estimate, no convergence study, and no benchmark for this quantity. The abstract and introduction describe a 'benchmarked correlated calculation,' yet no benchmark data appear in the body of the paper. Since C_He appears in the denominator of Eq. (9), any fractional error in C_He propagates one-to-one into every quoted shift, and any error in its λ-dependence distorts the curves in Fig. 2. The authors should provide a quantitative uncertainty on C_He(λ), for example by comparing the wavefunction of Ref. [24] with an independent high-accuracy correlated wavefunction, and by checking the λ→∞ limit against a known high-precision value of ⟨1/r12⟩ for the He 2^3S state. This uncertainty must be propagated into the shaded bands of Fig. 2 and into the headline ranges 0.250–0.850 MHz and 0.14–0.21 MHz.","section":"He-positronium conversion (Eqs. (7)-(9), Fig. 2)"},{"comment":"The statement that C_He(λ) 'is evaluated numerically using the He wave function of Ref [24]' is insufficient for reproducibility. The authors should specify the form of the wavefunction (basis type and number of terms), the integration method, and the numerical precision. This matters because the operator e^{-r12/λ}/r12 emphasizes short electron-electron separations, and the accuracy of the Ref. [24] wavefunction for this operator is not established by the original paper's intended use. Without this information, Eq. (9) cannot be independently checked, and the reliability of the central prediction remains unverified.","section":"Numerical evaluation of C_He (Eq. (8))"}],"minor_comments":[{"comment":"There is a typographical error in the integrand: 'r,dr' should be replaced by the correct integration measure, presumably e^{-r/λ} r dr (with all variables defined).","section":"Eq. (3)"},{"comment":"The 'projected sensitivity' curves should state the confidence level or statistical criterion used, so that the reader can interpret the exclusion power of a future 100 kHz measurement.","section":"Fig. 1 and text"},{"comment":"The quantity referred to as the '2^3S_1 ionization shift' should be defined precisely, for example as the shift of the 2^3S_1 binding energy relative to the free e+e- threshold, and the theoretical uncertainty assigned to this observable should be stated explicitly.","section":"Fig. 2 and text"},{"comment":"The propagation of the asymmetric ΔE_He uncertainty used for the shaded bands is described only as 'independently for its upper and lower bounds'; please state the exact formula used so that the bands can be reproduced.","section":"Eq. (9) and Fig. 2"},{"comment":"The citation [24] is a paper on exotic spin-dependent interactions; please clarify whether it contains the wavefunction used here or whether it should be supplemented by a more direct citation to the original helium wavefunction calculation.","section":"References"},{"comment":"The '0–800 eV' mediator-mass range should be more carefully motivated; the authors should state that this range is inherited from Ref. [6] and briefly list the criteria defining its upper limit.","section":"Abstract and text"}],"recommendation":"major_revision","confidential_remarks":"The central concern is the unquantified numerical matrix element C_He(λ); the editor should require an uncertainty estimate and a benchmark before acceptance. If the authors can supply that, the paper would be a solid contribution to the field. The manuscript is in scope for the journal and does not suffer from circular reasoning, since the positronium predictions are derived from the helium anomaly rather than fitted to positronium data."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper is worth your time. It takes the scalar-boson interpretation of the metastable helium ionization-energy anomaly and works out, for the first time, what that implies for positronium transitions. The analytic core is standard hydrogenic integrals, but the ratio F_nl/C_He combined with the measured ΔE_He gives concrete predicted shifts (0.25–0.85 MHz for 1^3S1–2^3S1) that are large enough to be testable with planned Ps spectroscopy. That is a genuinely new bridge between two experimental programs, and it is built honestly: there is no fitting to Ps data, the conditional status is stated clearly, and the uncertainties that are propagated are handled transparently.\n\nWhat the paper does well: the integrals (4)–(6) are correct, Eq. (9) follows cleanly, and the discussion of which Ps observables are most sensitive (and why 2S–2P fine structure is a dead end at low mediator mass) is clear and useful. The self-citations are appropriate; they are context, not load-bearing. The paper is careful to note that the whole program is conditional on the scalar interpretation.\n\nThe soft spot is exactly where the reader put it: C_He(λ) = ⟨e^{-r12/λ}/r12⟩_{2^3S} is taken from the wavefunction of Ref. [24] with no stated uncertainty, no convergence study, and no independent check. The abstract advertises a 'benchmarked correlated calculation,' but no benchmark appears in the manuscript, and the text itself calls the He state 'approximate.' Because C_He appears in the denominator of Eq. (9), a fractional error in C_He lands directly on every quoted Ps shift, and an error in its λ-dependence distorts the shape of Fig. 2. The paper propagates only the ΔE_He uncertainty, so the quoted ranges have no error budget from the wavefunction. This is a genuine correctness risk, not a stylistic quibble. It is also easily fixable: the authors should supply the numerical values, a convergence study in the basis size, and a benchmark against an independent He wavefunction or an analytic limit (e.g., the λ→∞ limit of C_He reduces to the known ⟨1/r12⟩ expectation value for He 2^3S, which is available in the literature). Without that, the central numbers are not independently checkable.\n\nMy overall take: this is a smart, useful paper that deserves a serious referee. The conditional framing is honest, the analytic part is solid, and the predicted shifts are big enough to matter. I would not cite it for the helium matrix element itself, but I would cite it for the Ps predictions and the framework. The authors should be asked to quantify C_He before publication; that is a revision, not a rejection.\n\nRecommendation: send to peer review. The one required change is an uncertainty or benchmark for C_He, plus the numerical values or code.\n\nBest,\n\n[Name]","headline":"A clean conditional prediction that converts the He anomaly into a Ps spectroscopy target; the one soft spot that matters is the unquantified helium matrix element.","tokens_in":9752,"tokens_out":2280,"would_cite":true,"duration_ms":19623,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"If a new scalar force explains the helium anomaly, positronium's 1S–2S interval must shift by 0.25–0.85 MHz, a detectable signal.","keywords":["positronium","helium ionization energy anomaly","scalar boson","Yukawa potential","precision spectroscopy","quantum electrodynamics","bound-state QED","new physics"],"falsifier":"Measure the positronium $1\\,^3S_1\\to2\\,^3S_1$ interval with total uncertainty at or below 100 kHz (requiring a next-generation QED calculation) and compare to the updated QED prediction: if the measured value agrees with QED to within 100 kHz across the mediator-mass range 0–800 eV, the predicted 0.250–0.850 MHz shift is ruled out, disfavoring the scalar-boson interpretation of the helium anomaly.","tokens_in":8718,"feed_emoji":"⚛️","tokens_out":9883,"duration_ms":75160,"temperature":0.7,"pith_summary":"This paper argues that if the recently reported 9σ discrepancy between the measured and calculated ionization energy of metastable helium is caused by a new scalar boson interacting with electrons, then positronium must show computable frequency shifts. Using analytic hydrogenic wavefunctions for positronium and a numerical helium matrix element, the authors derive that the $1\\,^3S_1\\to2\\,^3S_1$ transition shifts by 0.250–0.850 MHz and the 2S ionization energy by 0.14–0.21 MHz over the viable mediator-mass range. These shifts are large enough that planned precision positronium spectroscopy could detect them, making positronium a direct test of the scalar-boson explanation. The paper identifies which measurement is most sensitive and what a null or matching result would mean.","feed_headline":"Helium anomaly predicts positronium shift of 0.25-0.85 MHz","feed_subtitle":"A 0.25–0.85 MHz shift in positronium's 1S-2S line would confirm the new scalar force; its absence would refute it.","key_machinery":"The machinery is the first-order perturbation formula $\\Delta E_{n\\ell}=-g_e^2/(\\hbar c)\\,4\\pi F_{n\\ell}(\\lambda)$ for a Yukawa potential $V(r)=-g_e^2/(\\hbar c)\\,4\\pi e^{-r/\\lambda}/r$, together with the ratio identity $\\Delta\\nu_{\\mathrm{Ps}} = \\Delta E_{\\mathrm{He}}/h \\cdot F_{n\\ell}(\\lambda)/C_{\\mathrm{He}}(\\lambda)$. Here $F_{n\\ell}(\\lambda)$ is the radial integral $\\int_0^\\infty |R_{n\\ell}(r)|^2 e^{-r/\\lambda} r\\,dr$, evaluated analytically for hydrogenic positronium (with $a_{\\mathrm{Ps}}=2a_0$), and $C_{\\mathrm{He}}(\\lambda)=\\langle e^{-r_{12}/\\lambda}/r_{12}\\rangle_{2^3S}$ is the corresponding matrix element in the metastable helium $2\\,^3S$ state, computed numerically. This identity is what carries the argument: it converts the measured helium anomaly into a specific, mass-dependent shift for every positronium transition, with no free parameters beyond the mediator mass.","core_discovery":"The central discovery is a quantitative link between the helium ionization-energy anomaly and observable positronium spectra: the frequency shift of a positronium transition is exactly the helium anomaly times the ratio of a positronium radial integral $F_{n\\ell}(\\lambda)$ to the helium matrix element $C_{\\mathrm{He}}(\\lambda)$. The positronium integrals are analytic for a Yukawa perturbation, and the helium matrix element is computed numerically, so the ratio turns the measured helium shift into a definite prediction for positronium for every mediator mass. Over the mediator masses that can explain the anomaly (up to 800 eV), the $1\\,^3S_1\\to2\\,^3S_1$ shift is 0.250–0.850 MHz, and the $2\\,^3S_1$ ionization shift is 0.14–0.21 MHz. The paper concludes that the 1S–2S interval is the most sensitive probe, that the 2S-ionization measurement can independently constrain the radial form of the interaction, and that the 2S–2P fine-structure transitions are insensitive to light mediators.","pith_inferences":["The same ratio method could be applied to muonium, whose different reduced mass would give different radial overlaps; comparing positronium and muonium shifts could separate electron-only couplings from flavor-dependent couplings.","The exact analytic $F_{n\\ell}(\\lambda)$ expressions are generic response functions for any short-range perturbation of a hydrogenic system, so they could be reused to forecast shifts from other hypothetical forces (e.g., dark-photon or axion-like exchanges) in positronium and hydrogen.","A paired measurement of the 1S–2S and 2S-ionization shifts would not only test the scalar hypothesis but, if both are seen, would determine the mediator mass from their ratio, since that ratio is a known function of $\\lambda$.","The existing 5.7 MHz offset between the current best positronium 1S–2S measurement and QED theory is much larger than the predicted scalar shift; resolving that offset is a prerequisite for using positronium to test the helium anomaly, a point the paper leaves implicit."],"forward_implications":["If the scalar-boson interpretation is correct, a positronium $1\\,^3S_1\\to2\\,^3S_1$ measurement with ~100 kHz total uncertainty will see a shift between 0.250 and 0.850 MHz for mediator masses below 800 eV.","A null positronium result—no shift beyond QED uncertainty—would strongly disfavor the scalar-boson explanation of the helium anomaly.","A measurement of the $2\\,^3S_1$ ionization energy probes a different combination of states, so its 0.14–0.21 MHz shift would independently test the radial form of the interaction.","The $2\\,^3S_1\\to2\\,^3P_J$ fine-structure transitions are insensitive to very light scalar mediators because the $2S$ and $2P$ shifts cancel in the $\\lambda\\to\\infty$ limit, so they are not a viable test.","Reducing the theoretical QED uncertainty for positronium from 580 kHz toward 100 kHz is necessary to test the full mediator-mass range and would make the projected experimental sensitivity competitive with or stronger than the helium-anomaly constraint."],"supporting_citations":[{"why":"Provides the experimental ionization energy of metastable $^3$He that defines one arm of the anomaly.","marker":"[3]"},{"why":"Supplies the QED theoretical ionization energy of helium whose difference from experiment yields $\\Delta E_{\\mathrm{He}}$.","marker":"[4]"},{"why":"Provides the earlier $^4$He measurement showing the same anomaly, confirming isotope independence.","marker":"[5]"},{"why":"Establishes that only a scalar boson interaction is consistent with the helium observations; the paper's starting point.","marker":"[6]"},{"why":"Gives the current QED theoretical uncertainty (0.58 MHz) for the positronium 1S–2S interval used in sensitivity projections.","marker":"[21]"},{"why":"Provides the helium wavefunction used to compute the numerical matrix element $C_{\\mathrm{He}}(\\lambda)$ that links the anomaly to positronium.","marker":"[24]"},{"why":"Reports the most precise existing measurement of the positronium 1S–2S interval, the experimental baseline the new tests must surpass.","marker":"[28]"}],"fun_headline_variants":["Helium anomaly predicts positronium 1S–2S shift of 0.25–0.85 MHz","Positronium shift of 0.25–0.85 MHz tests helium anomaly","New scalar force would shift positronium by 0.25–0.85 MHz","Helium anomaly implies positronium 0.25–0.85 MHz shift"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions assume the numerical helium matrix element $C_{\\mathrm{He}}(\\lambda)$ — computed with the wavefunction of Ref. [24] — is accurate, and the paper assigns it no uncertainty; because the predicted positronium shifts scale inversely with $C_{\\mathrm{He}}$, any error in that matrix element changes all predicted shifts proportionally.","fun_headline_variants_meta":{"raw":{"variants":["Helium anomaly predicts positronium 1S–2S shift of 0.25–0.85 MHz","Positronium shift of 0.25–0.85 MHz tests helium anomaly","New scalar force would shift positronium by 0.25–0.85 MHz","Helium anomaly implies positronium 0.25–0.85 MHz shift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001308,"raw_usage":{"total_tokens":5347,"prompt_tokens":977,"completion_tokens":4370,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":4270}},"tokens_in":593,"tokens_out":4370,"duration_ms":27732,"temperature":1.0,"reasoning_tokens":4270,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:06:40.597882+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the positronium $1\\,^3S_1\\to2\\,^3S_1$ interval with total uncertainty at or below 100 kHz (requiring a next-generation QED calculation) and compare to the updated QED prediction: if the measured value agrees with QED to within 100 kHz across the mediator-mass range 0–800 eV, the predicted 0.250–0.850 MHz shift is ruled out, disfavoring the scalar-boson interpretation of the helium anomaly.","supporting_citations":[{"cited_title":"Ionization energy of metastable 3he (2 3s1) and the alpha- and helion-particle charge-radius difference from precision spectroscopy of thenprydberg series,","cited_arxiv_id":null,"evidence_quote":"Provides the experimental ionization energy of metastable $^3$He that defines one arm of the anomaly."},{"cited_title":"Completeα 7mlamb shift of helium triplet states,","cited_arxiv_id":null,"evidence_quote":"Supplies the QED theoretical ionization energy of helium whose difference from experiment yields $\\Delta E_{\\mathrm{He}}$."},{"cited_title":"The ionization energy of metastable triplet helium,","cited_arxiv_id":null,"evidence_quote":"Provides the earlier $^4$He measurement showing the same anomaly, confirming isotope independence."},{"cited_title":"PositroniumS-state spectrum: Analytic results ato(mα 6),","cited_arxiv_id":null,"evidence_quote":"Gives the current QED theoretical uncertainty (0.58 MHz) for the positronium 1S–2S interval used in sensitivity projections."},{"cited_title":"Constraints on exotic spin-dependent interac- tions between electrons from helium fine-structure spec- troscopy,","cited_arxiv_id":null,"evidence_quote":"Provides the helium wavefunction used to compute the numerical matrix element $C_{\\mathrm{He}}(\\lambda)$ that links the anomaly to positronium."},{"cited_title":"Measurement of the positronium 1 3s1–2 3s1 interval by continuous-wave two-photon excitation,","cited_arxiv_id":null,"evidence_quote":"Reports the most precise existing measurement of the positronium 1S–2S interval, the experimental baseline the new tests must surpass."}],"review_version":1}