{"id":"baac8512-14a6-471b-bedc-d53f50d3ee96","arxiv_id":"1908.05331","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A closed-form Glauber amplitude for ground-state elementary-atom Coulomb breakup is derived and compared with Born and dipole approximations for DIRAC-type atomic pair spectra.","lead":"This preprint derives analytic formulas for the momentum and angular spectra of charged meson pairs produced when pionium or kaonic atoms are ionized in the Coulomb field of a target. It compares a simple dipole approximation with first-Born and all-order Glauber calculations and finds the simple dipole shape is sufficient for approximate estimates.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The cross-section split at q0 in Eq (7) lacks a matching or interference justification; without a q0-stability check, the all-order Glauber result may not give the true breakup cross section.","rationale":"I read the paper as making two connected claims: first, a closed all-order Glauber amplitude for ionization of a ground-state elementary atom, and second, a total breakup cross section expressed as the sum of a screened Born term and an unscreened Glauber term. The amplitude derivation is long and partly delegated to an appendix, but the single most load-bearing step for the cross-section claim is Eq (7), because without it the unscreened all-order amplitude cannot be turned into a realistic cross section for the DIRAC target. The reader's weakest_assumption is exactly this point: q0 is chosen by magnitude, no matching or interference term is provided, and the authors themselves say the Born-Glauber differences are not experimentally checkable. I agree that this is the weakest link. The conditionality of the verdict is appropriate: a numerical q0-stability test could settle whether the split is harmless. I do not see a more fundamental internal inconsistency in the amplitude derivation itself, although the omitted algebra in the appendix would also deserve independent verification. The recommendation is therefore unchanged from the reader's CONDITIONAL verdict.","tokens_in":10846,"tokens_out":5232,"duration_ms":55441,"concrete_test":"Vary the matching momentum q0 in Eq (7) by factors 0.5, 2, and 5 for the Ni-target 1S pionium breakup, recompute dsigma/dp and the total cross section, and inspect the two integrands at q0 to see whether one is already negligible there. If the total cross section or spectrum changes by more than a few percent, the q0 split is not well-defined; if the result is flat across q0 choices, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result is the all-order Glauber amplitude in Eqs (27) and (34), but its practical use relies on Eq (7), which splits the breakup cross section into a screened first-Born term integrated up to q0 and an unscreened all-order Glauber term integrated from q0 to infinity. The boundary q0 ~ alpha (m_e mu Z^{1/3})^{1/2} is chosen by an order-of-magnitude estimate after Eq (7), with no matching condition ensuring that the Born term is accurate below q0 and the Glauber term accurate above it, and no interference term between the two. Because the two contributions are computed with different potentials (screened vs unscreened), adding cross sections rather than amplitudes is not equivalent to a single consistent amplitude calculation. If the two contributions overlap or interfere near q0, Eq (7) is not the true cross section, and the spectra in Section 6 inherit this uncertainty. The authors' own admission in Section 6 that the detailed Born-Glauber differences cannot be checked experimentally makes the analytical split the only support for the summed result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers the breakup (ionization) of relativistic Coulomb-bound meson pairs (elementary atoms, EA) in the Coulomb field of a target atom, with application to the DIRAC experiment. The authors present three approximations: a dipole approximation, a first-Born approximation with screened target potential, and an all-order Glauber treatment with an unscreened potential. The central formal result is the closed analytical expression for the ground-state EA breakup amplitude, Eqs. (27) and (34), which accounts for multiple photon exchanges with the target and all Coulomb interactions in the initial and final meson-pair states. To obtain a total cross section, Eq. (7) splits the integral over transverse momentum q into a screened Born term for q<q0 and an unscreened Glauber term for q>q0, with q0 chosen by an order-of-magnitude estimate. The paper compares relative momentum and angular spectra from the three approaches for 1S and 2S states and shows that the simple dipole approximation gives similar shapes, while acknowledging that the detailed differences cannot be checked with existing DIRAC data.","tokens_in":11014,"tokens_out":3137,"duration_ms":34553,"significance":"If the all-order Glauber amplitude in Eqs. (27) and (34) is correct, this is a genuinely useful closed-form result: it would be the first all-order treatment of EA breakup in a Coulomb field that includes all multiple exchanges and all initial/final-state Coulomb interactions, going beyond the first-Born approximation used in earlier DIRAC analyses. The analytical structure of the amplitude, involving a hypergeometric function and contour integrals over the final-state Coulomb parameter ξ, appears to follow from standard eikonal and Coulomb-wave-function methods, and the paper honestly reports the limitations of its comparison to experiment. However, the practical value of the result depends on the q0-split in Eq. (7), which is not derived or validated, and on the ability to produce normalized absolute spectra, which the paper does not demonstrate. The manuscript also provides no numerical cross-checks against independent calculations or against absolute experimental rates, so the central quantitative claim remains unsupported until the split is justified or tested.","major_comments":[{"comment":"The total breakup cross section is written as the sum of a screened first-Born cross section integrated up to q0 and an unscreened all-order Glauber cross section integrated from q0 to infinity, with no interference term and no matching condition. Because the two terms are computed with different potentials (screened versus unscreened), adding the cross sections is not equivalent to a single amplitude-level calculation. The boundary q0 ~ α(m_e μ Z^{1/3})^{1/2} is introduced by an order-of-magnitude estimate after Eq. (7), not derived from a requirement that the Born term is accurate below q0 and the Glauber term accurate above it. This is load-bearing: if the two contributions overlap or interfere near q0, the summed cross section is not the true physical cross section, and the comparison in Section 6 inherits this uncertainty. I request either a derivation or a quantitative demonstration that the result is independent of q0 over a reasonable range, or a full amplitude-level calculation with a screened potential in the Glauber term.","section":"Section 3, Eq. (7)"},{"comment":"The comparison of the dipole, Born, and Glauber spectra is presented only in arbitrary units with no absolute normalization and no numerical cross-check. Since the stated purpose is to compute breakup spectra that can be used in DIRAC analyses, the lack of any validation of the absolute normalization is a significant gap. The authors state themselves that the detailed Born–Glauber differences cannot be verified experimentally; therefore, the analytical split in Eq. (7) is the only support for the summed result, but no q0-stability test is provided. At minimum, the cross sections should be normalized to the total breakup cross section, and the q0 dependence of the integrated result should be shown.","section":"Section 6, Figs. 2–5"},{"comment":"The closed form for the ionization form factor S_p,nlm(q), which is central to the Born cross section in Eq. (16), is obtained after 'simple but cumbersome algebra' omitted between Eq. (52) and Eq. (53). Because this formula is the basis for the first term in Eq. (7) and for Figs. 1–5, the omission is not merely a presentation issue; it prevents the reader from verifying the derivation. I recommend that the intermediate steps be provided in full, either in the text or as supplementary material, or that a clear pointer be given to a derivation in a referenced work.","section":"Appendix, Eq. (53)"},{"comment":"The main Glauber amplitude expression appears to contain notation and factor inconsistencies that need clarification. In Eq. (28), the integrand includes |Γ(iν)|², while Eq. (22) suggests the natural factor is 1/|Γ(1+iν)|²; the relationship between these is not explained. Also, Eq. (27) defines A_fi(q) via a contour integral over t, but the integrand in Eq. (28) depends on κ = p_T(1-t) and on c defined in Eq. (26); the dependence of A_fi on the transverse momentum p_T and on the final-state angle θ is not made explicit. Since these equations constitute the paper's main result, the notation should be made fully consistent and every variable (including q, q, s, κ, and p_T) should be defined with its vector/scalar character.","section":"Section 5, Eqs. (27)–(34)"}],"minor_comments":[{"comment":"There are several typographical errors: 'haronic' should be 'hadronic', and the abstract uses 'h+h−' without proper spacing; these should be corrected in the final version.","section":"Abstract and Introduction"},{"comment":"The normalization constants C and C' are not given explicitly, yet they are needed to compare the dipole spectra with the Born and Glauber spectra even in arbitrary units. Please provide their definitions or state that they are chosen to normalize the plotted curves to unity.","section":"Eq. (5) and Eq. (6)"},{"comment":"The quantity q is described as a two-dimensional transfer momentum, but it also appears in scalar products and as a magnitude in q0 and in the integration measure q dq dφ. Please define q = |q| and consistently use boldface or an explicit magnitude notation to avoid confusion.","section":"Section 3, after Eq. (7)"},{"comment":"The vertical axes are labeled 'a.u.' with no explanation of the normalization; for a quantitative comparison, the spectra should be normalized to the same total cross section or the normalization procedure should be stated.","section":"Figs. 2–5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears squarely within the journal's scope and the formal Glauber result is potentially valuable. However, the central cross-section split in Eq. (7) and the absence of any q0-stability or absolute-normalization check are substantial enough that the paper should not be accepted in its present form. The authors' own admission that the detailed differences cannot be checked experimentally makes the analytical justification of the split essential. I would be willing to review a revised version that addresses the q0 issue and the omitted algebra."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper gives the first closed all-order Glauber ionization amplitude for a ground-state elementary atom in a Coulomb field, and a closed Born form factor for arbitrary n,l,m initial states. Those are real, checkable analytic results. The catch is that the practical cross section they compare with experiment depends on splitting the integral at an order-of-magnitude q0 and adding a screened Born cross section to an unscreened Glauber one with no interference term. The authors don't show the sum is stable against the choice of q0, and they admit the detailed differences can't be observed anyway.\n\nCredit first. Eq (16) is a genuine advance: previous approaches expanded the final-state wave function in spherical harmonics and truncated, whereas this is a finite closed sum over Gegenbauer and Jacobi polynomials, evaluable to arbitrary accuracy. Eqs (27) and (34) are also new; the manipulation of the Coulomb phase and the overlap integral is nontrivial and looks internally consistent. The paper is honest about scope: Section 6 says explicitly that the dipole approach is sufficient for simplified estimates and the Born-Glauber differences cannot be checked experimentally. That doesn't kill the value of the analytic formulas, but it does mean the impact is limited to refining the DIRAC analysis pipeline, not changing conclusions.\n\nSoft spots, in order. Eq (7) is the load-bearing approximation. Two cross sections computed with different potentials (screened vs unscreened) are added with no interference term and no matching condition. The boundary q0 ~ alpha (m_e mu Z^{1/3})^{1/2} is an estimate, not a derived scale. If the Born and Glauber regimes overlap near q0, the sum isn't the true cross section. A q0-stability scan or an argument about why the interference is negligible would close this. The Appendix skips 'simple but cumbersome algebra' leading to Eq (53); for a main result, that's an omission a referee should ask to be filled, even if it is routine. Figures are in arbitrary units and there's no numerical cross-check of the closed form against an independent method, which would have been easy for the Born case. These are separate minor holes.\n\nBottom line: the central amplitude derivation doesn't rest on the q0 split; the split is only used to produce final spectra. So the core math stands. The paper deserves a serious referee, but the referee should ask for a q0-stability check and more detail in the appendix before acceptance.\n\nWho this is for: members of DIRAC and people working on Coulomb breakup of hadronic atoms. Not a broad physics audience. If I were working on that analysis, I'd cite it for the closed form factor, not for the summed cross section.","headline":"Genuinely new closed analytic amplitudes for dimesoatom breakup, but the practical cross section rests on an unproven q0 split; still deserves a serious referee.","tokens_in":11589,"tokens_out":2972,"would_cite":false,"duration_ms":28821,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["36.10.-k","34.80.Dp"],"model":"deepseek-v4-flash","headline":"This paper claims a closed analytical, all-order Glauber expression for the ground-state ionization amplitude of an elementary meson atom in a Coulomb field, giving the breakup spectra used in mesic-atom lifetime measurements.","keywords":["elementary atoms","pionium","pion-kaon atoms","Coulomb breakup","Glauber approximation","Born approximation","multiple photon exchange","DIRAC experiment"],"falsifier":"Evaluate the ground-state breakup cross section by solving the full eikonal impact-parameter integral with a screened Coulomb potential, keeping all photon orders without splitting the q integral, and compare those spectra with the sum of the screened Born (q < q0) and unscreened Glauber (q > q0) terms of Eq (7) for the same target. If they differ by more than the small visible deviations between the Born and Glauber curves in the paper's figures, then the assumed additivity at the boundary q0 is the reason, and the main cross-section claim is falsified.","tokens_in":10597,"feed_emoji":"⚛️","tokens_out":12431,"duration_ms":113642,"temperature":0.7,"pith_summary":"This paper derives closed analytical formulas for the momentum and angular spectra of meson pairs produced when relativistic hydrogen-like meson atoms (pionium or pion-kaon atoms) are ionized by the Coulomb field of a target nucleus. The motivation is the DIRAC experiment, which extracts pion-pion and pion-kaon scattering lengths from the lifetime of these atoms and therefore needs accurate breakup spectra. The paper gives spectra in three approximations: a simple dipole limit, a first-Born treatment with target screening valid for any initial state, and a new all-order Glauber expression for the ground state that keeps every multiple photon exchange with the target as well as all Coulomb interactions inside the initial and final meson pair. The central claim is that Eqs (27) and (34) are the first closed analytical form of that all-order ionization amplitude.","feed_headline":"New closed formula covers all photon exchanges in pionium breakup","feed_subtitle":"An all-order Glauber amplitude gives ground-state breakup spectra the mesic-atom experiment can use directly.","key_machinery":"The central object is the atomic transition form factor S_{p,nlm}(q,p) = ∫ ψ_f^*(r) $e^{{i q·r}}$ ψ_i(r) d³r, which the paper evaluates in closed form as a finite sum over Gegenbauer and Jacobi polynomials using standard hydrogen-like Coulomb wave functions for the bound and continuum states. For the all-order amplitude, the load-bearing mechanism is the eikonal impact-parameter amplitude f(q,s) = (i/2π)∫ d²b [1 – $e^{{iΔχ(b,s)}}$] $e^{{i q·b}}$ with the unscreened Coulomb phase difference Δχ(b,s) = –ν ln[(b² + bs + s²/4)/(b² – bs + s²/4)], ν = Zα/β. Using hypergeometric integral representations, the paper converts this into the closed form of Eq (34), which contains a hypergeometric function F(iν, –iν; 1; 1 – c²q²/f²) together with derivatives with respect to the reduced-mass parameter and the integration variable t. That identity is what carries the claim of accounting for all multiple exchanges.","core_discovery":"The paper's central claim is that the breakup cross section of a relativistic elementary atom in a target Coulomb field separates into two momentum-transfer regimes. At small transfer momenta, target screening matters and a single-photon (first Born) exchange with a screened potential is sufficient; at large transfer momenta, screening can be neglected but all multi-photon exchanges must be included. The cross section is written as the screened first-Born integral up to a boundary q0 ~ α(m_e μ $Z^{{1/3}}$)^{1/2} plus an unscreened all-order Glauber integral from q0 to infinity. The new piece is the closed expression for the ground-state ionization amplitude, Eq (27) with Eq (34), obtained by evaluating the eikonal impact-parameter integral with the full Coulomb phase difference and Coulomb wave functions in the initial and final states. The paper also reduces the Born ionization form factor for arbitrary initial quantum numbers to a finite sum of Gegenbauer and Jacobi polynomials, and the computed spectra show that the dipole, Born, and Glauber approximations give similar shapes, with the relative-momentum peak narrowing as the principal quantum number n increases.","pith_inferences":["Going beyond the paper: replacing the hard q0 split with a smooth matching term derived from a full screened all-order evaluation would make the result more robust and could be checked numerically with the same amplitude machinery.","Going beyond the paper: the q0-split strategy is not specific to mesic atoms and could be tested on other relativistic hydrogen-like systems, such as muonium, where breakup spectra might be measured with different targets.","Going beyond the paper: because the paper's spectra are normalized absolutely, a future measurement of the absolute number of atomic breakup pairs per incident atom as a function of target Z would probe the q0 choice more directly than shape comparisons.","Going beyond the paper: since the final-state Coulomb interactions are fully resummed, the ground-state formula may also serve as a building block for Coulomb de-excitation of excited exotic atoms, where the cascade populates the states whose breakup feeds the measured spectrum."],"forward_implications":["The ground-state breakup spectra can now be computed with all multiple photon exchanges included, giving the mesic-atom experiment an all-order alternative to the first-Born approximation for the key ionization channel.","Because the Born form factor is closed for any initial quantum numbers, the spectra of excited atomic states can be computed to controlled accuracy, not only for low principal quantum numbers.","The similarity of the dipole, Born, and Glauber spectra means the simple dipole formula is adequate for rough estimates, while the exact spectra matter mainly for the lowest-n states that contribute most to the observed pair distribution.","The closed amplitude provides an absolute normalization for ground-state breakup on a target of atomic number Z, which can be compared with the measured number of atomic pairs in the experiment.","The relative-momentum peak lies near the mean momentum of the initial atomic state and narrows as the principal quantum number grows, so measured spectra carry information on the population of atomic states."],"supporting_citations":[{"why":"Anchors the pionium lifetime program that motivates the breakup spectra.","marker":"[1]"},{"why":"Introduces the breakup-pair detection method whose spectra this paper computes.","marker":"[4]"},{"why":"Provides the first-Born cross-section formula and transition form factor that the screened term uses.","marker":"[5]"},{"why":"Gives earlier numerical estimates of breakup spectra that the present closed forms supersede.","marker":"[11]"},{"why":"Source of the dipole-approach derivation reproduced here.","marker":"[12]"},{"why":"Supplies the eikonal Coulomb phase-shift formalism used for multiple exchanges.","marker":"[14]"},{"why":"Adapts the eikonal formalism to relativistic elementary atoms, underpinning the Glauber term.","marker":"[15]"},{"why":"Provides the hydrogen-like bound-state wave functions for the initial atom.","marker":"[16]"},{"why":"Provides the Coulomb continuum wave function of the final meson pair.","marker":"[17]"},{"why":"Contains the algebraic reduction that yields the closed Born form factors.","marker":"[23]"}],"fun_headline_variants":["All-order Glauber formula gives exact breakup spectra for pionium","Screened Born plus all-order Glauber: complete breakup cross section","Closed-form breakup amplitude for pionium with all multiphoton orders","Pionium breakup spectra now analytic for all photon-exchange orders","Exact all-order breakup amplitude for pionium in Coulomb fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the total breakup cross section is exactly the sum of a screened first-Born term for momentum transfers below the boundary q0 ≈ α(m_e μ $Z^{{1/3}}$)^{1/2} and an unscreened all-order term above q0, with no overlap or interference between the two; that boundary is set by an order-of-magnitude estimate rather than derived.","fun_headline_variants_meta":{"raw":{"variants":["All-order Glauber formula gives exact breakup spectra for pionium","Screened Born plus all-order Glauber: complete breakup cross section","Closed-form breakup amplitude for pionium with all multiphoton orders","Pionium breakup spectra now analytic for all photon-exchange orders","Exact all-order breakup amplitude for pionium in Coulomb fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000998,"raw_usage":{"total_tokens":4211,"prompt_tokens":914,"completion_tokens":3297,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":530,"completion_tokens_details":{"reasoning_tokens":3206}},"tokens_in":530,"tokens_out":3297,"duration_ms":23427,"temperature":1.0,"reasoning_tokens":3206,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:17:02.131306+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the ground-state breakup cross section by solving the full eikonal impact-parameter integral with a screened Coulomb potential, keeping all photon orders without splitting the q integral, and compare those spectra with the sum of the screened Born (q < q0) and unscreened Glauber (q > q0) terms of Eq (7) for the same target. If they differ by more than the small visible deviations between the Born and Glauber curves in the paper's figures, then the assumed additivity at the boundary q0 is the reason, and the main cross-section claim is falsified.","supporting_citations":[{"cited_title":"Afanasyev et al","cited_arxiv_id":null,"evidence_quote":"Anchors the pionium lifetime program that motivates the breakup spectra."},{"cited_title":"Nemenov, Yad","cited_arxiv_id":null,"evidence_quote":"Introduces the breakup-pair detection method whose spectra this paper computes."},{"cited_title":"Afanasyev, A.V","cited_arxiv_id":null,"evidence_quote":"Provides the first-Born cross-section formula and transition form factor that the screened term uses."},{"cited_title":"Afanasyev, A","cited_arxiv_id":null,"evidence_quote":"Gives earlier numerical estimates of breakup spectra that the present closed forms supersede."},{"cited_title":"Afanasyev, PhD Thesis, JINR, Dubna (1997)","cited_arxiv_id":null,"evidence_quote":"Source of the dipole-approach derivation reproduced here."},{"cited_title":"Gevorkyan, A.V","cited_arxiv_id":null,"evidence_quote":"Supplies the eikonal Coulomb phase-shift formalism used for multiple exchanges."},{"cited_title":"Afanasyev, S","cited_arxiv_id":null,"evidence_quote":"Adapts the eikonal formalism to relativistic elementary atoms, underpinning the Glauber term."},{"cited_title":"Sommerfeld, Ann","cited_arxiv_id":null,"evidence_quote":"Provides the hydrogen-like bound-state wave functions for the initial atom."},{"cited_title":"Landau, E.M","cited_arxiv_id":null,"evidence_quote":"Provides the Coulomb continuum wave function of the final meson pair."},{"cited_title":"Bakmaev, O","cited_arxiv_id":null,"evidence_quote":"Contains the algebraic reduction that yields the closed Born form factors."}],"review_version":1}