REVIEW 4 major objections 4 minor 24 references
Dimesoatom breakup in the Coulomb field
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict Genuinely new closed analytic amplitudes for dimesoatom breakup, but the practical cross section rests on an unproven q0 split; still deserves a serious referee. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [Section 3, Eq. (7)] 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 6, Figs. 2–5] 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.
- [Appendix, Eq. (53)] 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 5, Eqs. (27)–(34)] 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.
minor comments (4)
- [Abstract and Introduction] 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.
- [Eq. (5) and Eq. (6)] 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 3, after Eq. (7)] 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.
- [Figs. 2–5] 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.
Circularity Check
No significant circularity: Eqs. (27)/(34) are derived from standard Coulomb wave functions and eikonal theory; the q0 split in Eq. (7) is a heuristic input, not a fitted output.
full rationale
The central analytical result, Eqs. (27) and (34), is derived self-containedly by inserting the ground-state Coulomb wave function (23), the continuum Coulomb wave function (13), and the unscreened Coulomb eikonal phase (19) into the standard Glauber amplitude (8)-(10), followed by explicit contour and hypergeometric-function evaluations. The cited references [14,16,17] supply standard formulas (eikonal phase, Coulomb wave functions) rather than the target result, so the self-citations are not load-bearing in any circular sense. The only free input in the practical cross-section formula is the splitting momentum q0 in Eq. (7), chosen as q0 ~ alpha (m_e mu Z^{1/3})^{1/2} by an order-of-magnitude estimate from the inverse screening radius and EA Bohr momentum. That q0 is a heuristic boundary parameter, not a quantity fitted to the breakup spectra being predicted, and the paper does not fit any parameter to the DIRAC data shown in Fig. 6. The lack of a matching or interference justification around q0 is a genuine correctness risk, because Eq. (7) is an ansatz rather than a derived all-order cross section, but this is not circularity: the formula does not reduce to a fitted parameter renamed as a prediction, and no self-citation is invoked to forbid alternatives or to supply q0. The paper even admits in Section 6 that the detailed Born-Glauber differences cannot be checked experimentally, which further shows the comparison is exploratory rather than a closed loop. Overall, the derivation chain is independent of its own outputs, and any weakness is in the approximation strategy, not in circular reasoning.
Assumptions & free parameters
free parameters (2)
- q0 =
~ alpha (m_e mu Z^{1/3})^{1/2}
- Dipole normalization constants C and C' =
not specified (figures in arbitrary units)
assumptions (5)
- domain assumption Coulomb wave functions with no strong interaction describe the EA bound and continuum states.
- domain assumption Eikonal (Glauber) approximation is valid for the relativistic EA scattering on the target.
- domain assumption Atomic screening can be neglected in the Glauber term at q > q0.
- ad hoc to paper Born and Glauber cross-section terms can be added incoherently with the boundary q0.
- domain assumption Moliere parametrization of the Thomas-Fermi potential describes the target atom.
Cite this review
Pith. "Pith review of Dimesoatom breakup in the Coulomb field." pith.science (2026). https://pith.science/paper/5IOYFEM6
@misc{pith2026190805331,
author = {Pith},
title = {Pith review of: Dimesoatom breakup in the Coulomb field},
year = {2026},
howpublished = {\url{https://pith.science/paper/5IOYFEM6}},
note = {Machine review of arXiv:1908.05331}
}
abstract
Momentum and angular distributions of charged meson pairs $h^+h^-$ ($h=\pi, K$) from elementary atoms (EA) breakup (ionization) in the Coulomb field of a target atom is considered in the Born and Glauber approximations. Exploiting the fact that the atomic screening of the target Coulomb potential is important at small transfer momenta, while multi-photon exchanges are essential at large transfer momenta we express the cross sections of EA breakup as a sum of two terms. In the region of modest transfer momenta the cross section is determined by the single-photon exchange (first Born approximation) accounting for the target atoms screening, whereas at large transfer momenta using the unscreened potential allows to take into account all multi-photon exchanges and obtain the cross section of EA breakup in the close analytical form.
Reference graph
Works this paper leans on
- [1]
-
[2]
DIRAC Collaboration, Phys. Lett. B 704, 24 (2011)
work page 2011
-
[3]
DIRAC Collaboration, Phys. Rev. D 96, 052002 (2017)
work page 2017
-
[4]
L. Nemenov, Yad. Fiz. 41, 980 (1985) [Sov. J. Nucl. Phys. 41, 629 (1985)]
work page 1985
-
[5]
L.G. Afanasyev, A.V. Tarasov, Yad. Fiz. 59, 2212 (1996); Phys. Atom. Nucl. 59, 2130 (1996)
work page 1996
-
[6]
L. Afanasyev et al. , J. Phys. B: At. Mol. Opt. Phys. 37, 4749 (2004)
work page 2004
-
[7]
C. Santamarina, M. Schumann, L.G. Afanasyev, T. Heim, J. Phys. B: At. Mol. Opt. Phys.36, 4273 (2003)
work page 2003
- [8]
Show all 24 references
-
[9]
Heim et al
T.A. Heim et al. , J. Phys. B: At. Mol. Opt. Phys. 34, 3763 (2001)
2001
-
[10]
Schumann et al
M. Schumann et al. , J. Phys. B: At. Mol. Opt. Phys. 35, 2683 (2002)
2002
-
[11]
Afanasyev, A
L. Afanasyev, A. Tarasov, O. Voskresenskaya, Spectrum of pions from π+π−-atom breakup. In Proc. Workshop on Hadronic Atoms HadAtom01, Bern, 2001, p. 7 arXiv:hep- ph/0112293; T.A. Heim et al., Distribution of pions from breakup of pionium , ibidem p. 13
2001
-
[12]
Afanasyev, PhD Thesis, JINR, Dubna (1997)
L. Afanasyev, PhD Thesis, JINR, Dubna (1997)
1997
-
[13]
Amirkhanov, I
I. Amirkhanov, I. Puzynin, A. Tarasov, O. Voskresen- skaya, O. Zeinalova, Phys. Lett. B 452, 155 (1999) arXiv:hep-ph/9810251
1999 arXiv
-
[14]
Gevorkyan, A.V
S.R. Gevorkyan, A.V. Tarasov, O.O. Voskresenskaya, Phys. At. Nuclei 61, 1517 (1998)
1998
-
[15]
Afanasyev, S
L. Afanasyev, S. Gevorkyan, O. Voskresenskaya, Eur. Phys. J. A 53, 78 (2017)
2017
-
[16]
Sommerfeld, Ann
A. Sommerfeld, Ann. Phys. 11, 257 (1931)
1931
-
[17]
Landau, E.M
L.D. Landau, E.M. Lifshitz, Quantum Mechanics (Per- gamon, Oxford, 1977)
1977
-
[18]
Omidvar, Phys
K. Omidvar, Phys. Rev. 188, 140 (1969)
1969
-
[19]
Barut, R
A.O. Barut, R. Wilson, Phys. Rev. A 40, 1340 (1989)
1989
-
[20]
Gradshtein, I.M
I.S. Gradshtein, I.M. Ryzhik, Tables of Integrals, Series and Products. (Nauka Publication, Moscow, 1971)
1971
-
[21]
Abramowitz, I.A
Handbook of Mathematical Functions , edited by M. Abramowitz, I.A. Stegun (Dover, New York, 1972)
1972
-
[22]
Varshalovich, A.H
D.A. Varshalovich, A.H. Moskalev, A.H. Khersonsky, Quantum Theory of Angular Momentum (Nauka Pub- lication, Moscow, 1975); L.C. Biedenharn, L.C. Louck, Angular Momentum in Quantum Physics (Addison-Wesley, Massachusetts, 1981)
1975
-
[23]
Bakmaev, O
S. Bakmaev, O. Voskresenskaya, Part. Nucl. Lett. 3, 33 (2006)
2006
-
[24]
Zhabitsky, DIRAC Note 2007-11; cds.cern.ch/record/ 1369651
M. Zhabitsky, DIRAC Note 2007-11; cds.cern.ch/record/ 1369651
2007
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.