{"id":"3e46c273-fff7-418b-b985-461c1290b80c","arxiv_id":"2608.03795","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The ratio of pionium formation probabilities in ns-states stays at (n2/n1)^3 to O(10^-3) under first-order strong-interaction corrections, and the Deser level shift is rewritten in terms of the strong-interaction radius and the wave-function correction at the origin.","lead":"This paper argues that the strong nuclear force barely changes the ratio of how often pionium atoms form in different excited states, even though it strongly distorts the wave function at the origin, and it rewrites the standard Deser energy-shift formula using a strong-interaction radius and a wave-function correction that it hopes could be measured.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The O(10^-3) n-independence of c_n rests on unshown integrals; the flagged O(10^-2) log terms are never demonstrated to be n-independent, so the accuracy claim of Eq. (66) is not fully established.","rationale":"The paper is a derivation of an advertised cancellation: strong-interaction O(1) corrections to the wave function at the origin factor out and cancel in ns formation ratios because the relative correction c_n is n-independent to 10^-3. The Deser reformulation (Eqs. (74)-(75)) is an algebraic identity once Eq. (61) is accepted. The key risk is exactly where the reader placed it: the unshown evaluation behind Eq. (61) and the n-independence of the log-enhanced contributions the paper itself had flagged. My proposed test directly checks that gap. I find no internal inconsistency in the Zeldovich formalism as written, and the numerical checks in [25] and the explicit estimates (72)-(73) give partial support, so the honest verdict remains CONDITIONAL rather than REJECT or UNVERDICTED. I do not change the reader's verdict.","tokens_in":17433,"tokens_out":11827,"duration_ms":141366,"concrete_test":"Re-derive Eq. (61) from Eqs. (57)-(60) with the explicit Coulomb reduced functions χ_n0, and evaluate c_n = R_n0(0) for n=1,2,3,4 using a Yukawa potential U_s(r)=g e^{-r/r_s}/r with r_s=1 fm, g fixed to give the pionium scattering length a_s≈0.22/m_π, and r_B=367 fm, separating the O(a_s/r_s), O((a_s/r_B)ln(r_B/r_s)), and polynomial pieces. If any |c_n - c_1| exceeds 10^-3, or if an independent numerical solution of the radial Schrödinger equation with the same short-range plus Coulomb potential gives |ψ_n(0)|^2/(π^{-1}(μα)^3 n^{-3}) that deviates from n-independence by more than 0.1%, the O(10^-3) accuracy of Eq. (66) fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of Eq. (66) — that strong interaction leaves the ns formation ratio at (n2/n1)^3 + O(10^-3) — rests entirely on c_n = Δψ_ns(0)/ψ_ns^c(0) being n-independent to O(a_s/r_B). The decisive object is Eq. (61): R_n0(0) = -∫ r^2 e^{-2νr/n} r U_s(r)[1/r - 2ν ln(νr) + P_{2n}(νr)] dr. This formula is introduced with 'The calculations lead to the following result' and its derivation is not displayed. The following sentence asserts that the n-dependent part is O(a_s/r_B) and that the O((a_s/r_B)ln(r_B/a_s)) part is n-independent, with justification only from Eq. (62) and the first two universal terms in Eq. (65). Those two facts do not by themselves control the full integral over the short-range region: at r ~ r_s the polynomial P_{2n}(νr) can contribute at linear order in νr, and a coefficient that varies with n would enter at the O(r_s/r_B) level; a log-enhanced O(10^-2) piece with n-dependence would directly spoil the O(10^-3) claim. Section 2.1 itself flags that the plane-wave route omits terms of order (a_s/r_B)ln(r_B/r_s) ~ 10^-2 and says their accuracy 'requires further study'; the Zeldovich section is supposed to remedy this, but the needed n-independence of the log terms is an unshown calculation. This is a missing verification of the load-bearing accuracy estimate, not an observed contradiction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the effect of the strong interaction on the formation probabilities and energy shifts of hadronic (pionium) atoms. It argues that in first-order perturbation theory the strong interaction modifies the ns Coulomb wave function at the origin by a state-independent relative correction c_n, up to terms of order a_s/r_B ~ 10^-3. Consequently the ratio w_{n1s}/w_{n2s} remains (n2/n1)^3 + O(10^-3), Eq. (66), despite strong-interaction-induced distortions that are O(1). The paper also rewrites the Deser level-shift formula in terms of an effective strong-interaction radius <r>_s and c_n, Eqs. (74)-(75), and gives numerical estimates for the 1s and 2s-2p shifts.","tokens_in":17869,"tokens_out":8029,"duration_ms":79074,"significance":"The standard derivation of the Deser formula (68)-(70) is clean, and the numerical estimates in Eqs. (72)-(73) are consistent with the stated scattering lengths. If the n-independence of c_n were rigorously established, Eq. (66) would provide a useful robustness argument for the DIRAC determination of pion scattering lengths. The paper explicitly identifies the limitation of the plane-wave approximation, which is commendable. However, the central O(10^-3) accuracy claim is not yet demonstrated: it rests on integrals (Eqs. (61)-(63)) whose derivation is not shown. The advertised Deser reformulation is an algebraic identity rather than a new physical relation. The paper is promising but requires a major revision to substantiate the headline claim.","major_comments":[{"comment":"The central accuracy claim, Eq. (66), rests on the assertion that R_n0(0) in Eq. (61) is n-independent up to O(a_s/r_B). This equation is introduced with 'The calculations lead to the following result' and no derivation. The explanation after Eq. (62) is insufficient: Eq. (62) fixes only the linear term of the Coulomb wave function at r=0, while the integrand in Eq. (61) extends over r ~ r_s, where the polynomial P_{2n}(νr) contributes at O(νr) with coefficients that in principle depend on n. A logarithmically enhanced term of order (a_s/r_B) ln(r_B/r_s) ~ 10^-2 is acknowledged in Section 2.1, and the claim that this term is n-independent is not demonstrated. Please provide the derivation of Eqs. (61) and (63) or an explicit n-dependent error bound.","section":"Section 2.2, Eqs. (61)-(63)"},{"comment":"Equation (75) is an algebraic identity rather than a new relation: substituting <r>_s = a_s/c_n from Eq. (74) returns exactly the standard Deser formula (70). The paper should state explicitly that this is a parameterization, and that the only nontrivial content is the n-independence of c_n. As written, the advertised 'formulation ... in terms of the effective radius' overstates the novelty. This does not invalidate the rest of the paper, but the presentation should be corrected.","section":"Section 3, Eqs. (74)-(75)"},{"comment":"The factorized form ψ_n0(r) = R_n(r)ψ_n0^(0)(r) with R_n independent of n is stated for an arbitrary potential U_s. Equation (63) contains an O(r_s/r_B) error term, and for r comparable to r_B the n-dependence of ψ_n0^(0)(r) may enter the factorization. The numerical support cited as [25] is a preprint covering the first four states and does not establish a universal statement for all n. Please state the domain of validity in r and n, or supply a proof that the remainder is n-independent.","section":"Section 2.2, Eqs. (63)-(64)"}],"minor_comments":[{"comment":"The ratio should read (n2/n1)^3, not (n2/n2)^3.","section":"Section 1.1, Eq. (12)"},{"comment":"The sign of c_n is inconsistent with Eq. (37) and with Eq. (74): Eq. (41) should be c_n = -∫ U_s(r) r dr + O(10^-3) if U_s is defined as in Eq. (36).","section":"Section 2.1, Eq. (41)"},{"comment":"The final equality should be '= (n2/n1)^3 + O(10^-3)', not '= O(10^-3)' as written.","section":"Section 2.1, Eq. (42)"},{"comment":"The integrand appears to contain an extra factor r ('r^2 e^{-2νr/n} r U_s(r)'); please check dimensions and the intended measure.","section":"Section 2.2, Eq. (61)"},{"comment":"Use 'eV' rather than 'эВ' and a decimal point in '-0.45 eV' for consistency with the rest of the text.","section":"Section 3, Eqs. (72)-(73)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the standard parts are correct, but the advertised O(10^-3) n-independence result is not yet supported by a displayed derivation. I would suggest requesting a complete derivation of Eqs. (61)-(63) or an explicit numerical/n-dependent error analysis before publication. The reformulation in Section 3 should also be reframed as a parameterization."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version. The paper argues that the (n2/n1)^3 formation ratio for pionium ns-states survives first-order strong-interaction corrections at the 10^-3 level, and it rewrites the Deser shift in terms of an effective strong-interaction radius and the relative correction to the Coulomb wave function at the origin. The first claim is plausible and, if true, practically reassuring for DIRAC-style analyses. The second claim is true but essentially definitional: Eq. (74) defines <r>_s = a_s/c_n, so Eq. (75) is just the standard Deser formula with a_s repackaged. That's not a problem with the physics, but it's not a new dynamical relation either.\n\nWhat the paper does well: it sets up the Zeldovich perturbation method cleanly and correctly identifies that earlier plane-wave treatments dropped ~10^-2 log-enhanced terms. The numerical estimates for the 1s and 2s shifts (Eqs. (72)-(73)) check out — I recomputed and get about -3.6 eV and -0.45 eV. The paper is also upfront about model dependence and about the large uncertainties in absolute formation probabilities.\n\nThe soft spot is the load-bearing step. The n-independence of c_n to O(10^-3) rests on Eq. (61), which is introduced with \"the calculations lead to the following result\" and no derivation. The discussion around Eqs. (62)-(65) explains why the leading terms are n-independent, but it does not control the full integral. The author explicitly notes (end of Sec. 2.1) that the plane-wave route omits terms of order (a_s/r_B)ln(r_B/r_s) ~ 10^-2, and claims that the Zeldovich version eliminates them. But the n-independence of those log terms is asserted, not demonstrated. If they carried n-dependence at the 10^-2 level, the O(10^-3) accuracy of Eq. (66) would fail. Numerical support is promised from same-group refs [25, 26], but [26] is unpublished. This is a missing verification, not a contradiction, but it is exactly where a referee should push.\n\nWho this is for: hadronic-atom phenomenologists, especially the DIRAC people. They will find the reformulation handy and the robustness claim encouraging, and the paper is worth engaging. But the O(10^-3) claim should be treated as provisional until the derivation of (61) is shown.\n\nRecommendation: send it to peer review. It is not a finished proof, but it is a legitimate contribution with a clear, fixable gap. I would not cite the accuracy claim as established, but I would engage with the paper and ask for the missing calculation.","headline":"A plausible but incomplete robustness claim for the n^-3 formation ratio; the Deser reformulation is a true but essentially definitional identity.","tokens_in":18431,"tokens_out":3793,"would_cite":false,"duration_ms":37572,"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":"The paper argues that strong interactions shift but do not spoil the Coulomb n^-3 scaling of hadronic-atom formation probabilities, and recasts the Deser level-shift formula in terms of an effective strong-interaction radius and an n-indepe","keywords":["hadronic atoms","pionium","Deser formula","scattering length","strong interaction","Coulomb wave function","perturbation theory","energy-level shifts"],"falsifier":"Evaluate c_n with the full Coulomb continuum wave functions for n=1,2,3 using a Yukawa strong potential of range 1/m_ρ, and check whether c_2-c_1 remains below 10^-3; alternatively, measure the 2s/1s or 3s/1s pionium formation-probability ratio experimentally at sub-percent precision and test (n2/n1)^3. Either would settle the n-independence claim.","tokens_in":17220,"feed_emoji":"⚛️","tokens_out":9565,"duration_ms":98582,"temperature":0.7,"pith_summary":"Hadronic atoms are hydrogen-like bound states of two oppositely charged hadrons, and their strong-interaction level shifts are a precision probe of low-energy hadron scattering. This paper tries to establish that the simplest Coulomb predictions for these atoms are far more stable than the wave functions themselves: even though the strong interaction can distort the s-wave function at the origin by order one, the distortion is essentially the same for every principal quantum number n. If that holds, the ratio of formation probabilities in different ns states keeps the pure-Coulomb value (n2/n1)^3 to an accuracy of about 0.1%, so the strong-interaction uncertainties that plague absolute probabilities largely cancel in the ratios used by pionium-lifetime experiments. The paper also rewrites the standard Deser level-shift formula in terms of the mean radius of the strong-interaction region and the n-independent relative wave-function correction, giving a route to measure that radius from the 2s-2p splitting. The argument is carried by a first-order perturbation method that expresses the correction through the unperturbed wave function of the very state being corrected.","feed_headline":"Pionium's n-cubed formation law survives strong force to 0.1%","feed_subtitle":"Reformulated Deser shift turns the 2s-2p splitting into a strong-interaction radius measurement.","key_machinery":"The mechanism is a small-distance identity of the Coulomb wave functions: the logarithmic derivative of the reduced ns wave function at the origin, |dψ_n0/dr|_{r=0}/ψ_n0(0) = -μ α, is the same for all n, and the first two terms of the small-r expansion of n^{3/2}χ_n0(r) are also n-independent (Eqs. 62, 65). These coincidences make the integrals that define the first-order relative correction R_n0(0) differ from state to state only by terms of order r_s/r_B ~ 10^-3. The derivation uses a perturbation-theory variant that expresses the first-order correction to a discrete state in terms of that same state's unperturbed wave function and a second linearly independent solution, so it bypasses the","core_discovery":"The central claim, stated the way the paper's author would state it, is that in first-order perturbation theory the relative strong-interaction correction to the ns Coulomb wave function at the origin, c_n = Δψ_ns(0)/ψ_ns^c(0), is independent of n up to terms of order a_s/r_B ~ 10^-3, for any form of the strong potential. From this it follows that the formation-probability ratio w_n1s/w_n2s remains (n2/n1)^3 + O(10^-3) (Eq. 66) even though the absolute wave functions at short distances are strongly modified. The same n-independence lets the Deser formula be reformulated as ΔE_ns = -(2π/μ) <r>_s c_n |ψ_n0^c(0)|^2 (Eq. 75), where <r>_s = a_s/c_n is the mean radius of the strong-interaction reg","pith_inferences":["The n-independence rests on a special property of the Coulomb s-state family near the origin; analogous ratios for non-s states or for states where the wave function vanishes at the origin would not be protected, so the result should not be read as a general statement about all formation channels.","The same perturbative machinery should apply to other short-range modifications of the Coulomb problem, e.g. finite nuclear size or vacuum-polarization potentials; the prediction is that their leading correction to formation-probability ratios also cancels at the 10^-3 level.","If the full Coulomb-continuum calculation the paper flags as needed were carried out and showed 10^-2-level n-dependence, the central accuracy claim would collapse; this makes a numerical evaluation of c_2-c_1 with exact continuum wave functions the natural check.","A dedicated pionium experiment comparing 3s/2s or 3s/1s formation rates directly would be a sharper test of the n-independence than the current lifetime measurement, because it isolates the wave-function ratio without relying on annihilation widths."],"forward_implications":["The pure-Coulomb ratio w_{n1s}/w_{n2s} = (n2/n1)^3 survives strong interactions to O(10^-3), so pionium formation-ratio inputs to lifetime analyses can be treated as Coulomb-like even if the short-distance wave functions are badly distorted.","Measuring the 2s-2p energy splitting together with the ground-state formation probability (or c_1) determines the mean strong-interaction radius <r>_s, complementing the scattering lengths extracted from the lifetime.","The Deser formula's factorization into strong and electromagnetic factors persists in the reformulated version, now with the relative wave-function correction c_n and radius <r>_s carrying the strong part.","Free-pair to bound-state formation ratios w_ps/w_ns also retain their Coulomb value to about 10^-3 in the same approximations.","Absolute formation probabilities remain model-dependent at the tens-of-percent level; only the ratios are protected by the n-independence."],"supporting_citations":[{"why":"supplies the original Deser formula for the strong shift of the ns level in terms of the scattering length.","marker":"[8]"},{"why":"supplies the formation-probability expression and the pure-Coulomb n^-3 law that the paper generalizes.","marker":"[14]"},{"why":"gives one earlier model-potential estimate of the ground-state strong correction that omits the continuous spectrum.","marker":"[15]"},{"why":"gives the delta-potential estimate of the same correction, another baseline the paper shows to be incomplete.","marker":"[16]"},{"why":"supplies the perturbation-theory method used to express first-order wave-function corrections without summing the full spectrum.","marker":"[20]"},{"why":"shows the omitted continuous-spectrum contribution can be of order a_s/r_s, motivating the stronger correction and the plane-wave estimate.","marker":"[21]"},{"why":"provides the continuous-spectrum corrections and the free-pair-to-bound-state formation ratio used in the paper.","marker":"[22]"},{"why":"supplies numerical calculations confirming the n-independence of the relative correction and its potential-model sensitivity.","marker":"[25]"},{"why":"provides the 2s-2p splitting, scattering-length inputs, and the experimental proposal that make the reformulated Deser formula testable.","marker":"[27]"}],"fun_headline_variants":["Pionium formation law holds to 0.1% under strong force","Deser shift reformulated via strong-interaction effective radius","Strong force doesn't break pionium's n^3 formation ratios","Pionium wave-function correction is n-independent","New Deser formula ties energy shift to strong radius"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The result depends on the omitted (a_s/r_B) ln(r_B/r_s) ~ 10^-2 contributions—dropped by using plane waves for the continuous spectrum—being independent of n; if those terms varied with n at the 10^-2 level, the claimed 10^-3 accuracy of the formation-probability ratio would fail.","fun_headline_variants_meta":{"raw":{"variants":["Pionium formation law holds to 0.1% under strong force","Deser shift reformulated via strong-interaction effective radius","Strong force doesn't break pionium's n^3 formation ratios","Pionium wave-function correction is n-independent","New Deser formula ties energy shift to strong radius"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000197,"raw_usage":{"total_tokens":1179,"prompt_tokens":699,"completion_tokens":480,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":443,"completion_tokens_details":{"reasoning_tokens":396}},"tokens_in":443,"tokens_out":480,"duration_ms":5920,"temperature":1.0,"reasoning_tokens":396,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T12:22:58.584880+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate c_n with the full Coulomb continuum wave functions for n=1,2,3 using a Yukawa strong potential of range 1/m_ρ, and check whether c_2-c_1 remains below 10^-3; alternatively, measure the 2s/1s or 3s/1s pionium formation-probability ratio experimentally at sub-percent precision and test (n2/n1)^3. Either would settle the n-independence claim.","supporting_citations":[{"cited_title":"Deser, M.L","cited_arxiv_id":null,"evidence_quote":"supplies the original Deser formula for the strong shift of the ns level in terms of the scattering length."},{"cited_title":"Efimov, M","cited_arxiv_id":null,"evidence_quote":"gives one earlier model-potential estimate of the ground-state strong correction that omits the continuous spectrum."},{"cited_title":"Bel'kov, V.N","cited_arxiv_id":null,"evidence_quote":"gives the delta-potential estimate of the same correction, another baseline the paper shows to be incomplete."},{"cited_title":"Baz', Ya.B","cited_arxiv_id":null,"evidence_quote":"supplies the perturbation-theory method used to express first-order wave-function corrections without summing the full spectrum."},{"cited_title":"Kuraev, Yad","cited_arxiv_id":null,"evidence_quote":"shows the omitted continuous-spectrum contribution can be of order a_s/r_s, motivating the stronger correction and the plane-wave estimate."},{"cited_title":"Afanasyev and O","cited_arxiv_id":null,"evidence_quote":"provides the continuous-spectrum corrections and the free-pair-to-bound-state formation ratio used in the paper."},{"cited_title":"A method for research on behaviour of a dimesoatomic wave function at small distances","cited_arxiv_id":"hep-ph/9812293","evidence_quote":"supplies numerical calculations confirming the n-independence of the relative correction and its potential-model sensitivity."},{"cited_title":"HadAtom01","cited_arxiv_id":"hep-ph/0112293","evidence_quote":"provides the 2s-2p splitting, scattering-length inputs, and the experimental proposal that make the reformulated Deser formula testable."}],"review_version":1}