REVIEW 2 major objections 6 minor 50 references
Creating true muonium via charmonium radiative decay
T0 review · 2 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read The paper predicts that the radiative decay J/psi -> gamma + (mu+mu-)_bound produces true para-muonium at a branching fraction of about 7 x 10^-13, and argues that a future super tau-charm facility could detect a few events per year.
desk verdict A clean parameter-free QED prediction for J/psi -> gamma + true muonium at ~7e-13; worth refereeing despite a mislabeled selection rule and a speculative reach section. 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 imaginary part of the nonrelativistic Coulomb Green function for the muon-antimuon pair, evaluated at zero separation. Its bound-state poles give the squared wave function at the origin, |psi_n(0)|^2 = $alpha^{3}$ $m_mu^{3}$/(8 pi $n^{3}$), and its branch cut gives the Sommerfeld-Schwinger-Sakharov enhancement factor for the continuum. The paper's key move is to express the desired decay width as the free-pair width times this Green function's imaginary part, then take the ratio with the J/psi -> mu+mu- width so the charmonium wave function R(0) cancels; the remaining input is the measured J/psi leptonic branching fraction.
What would settle it
Measure the J/psi -> gamma mu+mu- differential rate near threshold at a super tau-charm facility with 3.4 x $10^{12}$ J/psi events, and look for the bound-state peak at E = -$alpha^{2}$ m_mu/4 with the predicted height corresponding to a 7 x $10^{-13}$ branching fraction; alternatively, perform an explicit one-loop calculation of the two 'vanishing' diagrams to test the parity argument.
Extended reading notes
Core claim
On the paper's own terms, the discovery is that J/psi -> gamma + (mu+mu-)_bound is a quantitatively predicted and experimentally meaningful channel for producing true para-muonium. Using a nonrelativistic spin-triplet projector for the J/psi, the final-state-radiation amplitude, and the imaginary part of the Coulomb Green function for the muon pair, the authors obtain R|E<0 ~ 1.18 x $10^{-11}$ for the ratio of the bound-state decay width to the J/psi -> mu+mu- width, which translates into Br(J/psi -> gamma (mu+mu-)) ~ 7.03 x $10^{-13}$ when multiplied by the measured 5.961% leptonic branching fraction. The same machinery gives the above-threshold continuum rate, controlled by the Sommerfeld factor; in a threshold window as narrow as the ground-state binding energy it is comparable to the bound-state contribution (R|E>0 ~ 4.91 x $10^{-12}$), and it grows quickly with the window size. The paper then confronts this rate with detector reality: with 3.4 x $10^{12}$ J/psi per year at a future super tau-charm facility, one expects 2-3 bound-state events per year, and the 1.3 mm lab-frame decay length of the boosted true muonium could be separated with vertex resolution at the 0.4 mm level.
Load-bearing premise
The calculation assumes, without derivation, that parity makes the two Feynman diagrams with the photon emitted from the J/psi itself vanish, so only final-state muon radiation contributes; if that cancellation fails, the predicted branching fraction changes.
Editorial extensions
If this is right
- With 10^10 accumulated J/psi events, BESIII would expect well below one bound-state event, so current data cannot test this rate.
- At the proposed super tau-charm facility, 3.4 x 10^12 J/psi per year translates to 2-3 bound-state events per year, making the search statistically plausible but far from easy.
- The lab-frame decay length of the boosted true muonium is about 1.3 mm, so a vertex resolution near 0.4 mm would separate the signal from prompt QED backgrounds.
- Widening the muon-pair energy window to the MeV scale boosts the continuum contribution toward R ~ 1e-8, meaning thousands of events per year at STCF, though those are not bound-state signals.
- Because the prediction is normalized to the measured J/psi -> mu+mu- rate, the QED part is free of hadronic uncertainties, so a future measurement would test the Coulomb-resummation formalism directly.
Reading between the lines
- The same ratio method should transfer directly to Upsilon decays into gamma plus true tauonium, with the tau decay width kept finite; the paper mentions tauonium only in passing, so a dedicated calculation would be a natural next step.
- The parity argument that eliminates the charmonium-radiation diagrams is asserted without derivation; an explicit one-loop check would either confirm the 7e-13 prediction or reveal an additional amplitude.
- A realistic search might target the threshold-enhanced continuum rather than the bound-state peak, since the paper's Eq. (35) shows the continuum rate grows rapidly with the experimental energy resolution.
- The ratio R cancels the J/psi wave function, so this channel is a rare case where a hadronic decay's QED part is computable almost parameter-free; measuring it would probe Coulomb resummation in an entirely new regime.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper computes the branching fraction for J/ψ decaying into a real photon plus a bound para-muonium state, J/ψ → γ + (μ⁺μ⁻)_{bound}. Using a nonrelativistic spin projector for the charmonium and keeping only final-state-radiation diagrams, the authors derive the differential width near the μ⁺μ⁻ threshold, resum Coulomb ladder exchanges via the Green-function formalism, and integrate over the bound-state poles. The ratio R = Γ(J/ψ→γ(μ⁺μ⁻)_{bound}) / Γ(J/ψ→μ⁺μ⁻) is found to be R|E<0 ≈ 1.18 × 10⁻¹¹, which, multiplied by the measured Br(J/ψ→μ⁺μ⁻) = 5.961%, gives Br(J/ψ→γ(μ⁺μ⁻)_{bound}) ≈ 7.03 × 10⁻¹³. The paper also estimates the above-threshold continuum contribution in a chosen energy window and discusses detection prospects at BESIII and the proposed Super Tau-Charm Facility.
Significance. If correct, the calculation provides a concrete, QED-dominated production mechanism for true muonium with a clean normalization to the measured J/ψ→μ⁺μ⁻ branching fraction: the ratio R is, up to the charm-mass choice, parameter-free, since the nonperturbative wave function R(0) cancels. The phase-space integrals and the Coulomb Green-function treatment are explicit and checkable, and the numerical factor is internally consistent. The C-parity selection rule that eliminates the charmonium-radiation diagrams is physically important and correctly identifies FSR as the only leading contribution, although the paper mislabels it as parity. The result is of direct interest to the true-muonium program and to the planning of future tau-charm facilities. The main weakness is the sensitivity of the numerical prediction to the charm-quark mass, and the experimental feasibility section rests on assumptions about neutral-vertex resolution that are not demonstrated.
major comments (2)
- [Sec. V, Eqs. (36)-(39)] The numerical prediction uses m_c = 1.27 GeV, but the nonrelativistic derivation in Sec. II relies on M_J/ψ ≈ 2m_c, which with M_J/ψ = 3.097 GeV gives m_c ≈ 1.55 GeV. The PDG value 1.27 GeV is an MS-bar running mass, not the pole mass appropriate for the nonrelativistic bound-state expansion used in the spin projector. Since R|E<0 ∝ 1/m_c², replacing m_c by M_J/ψ/2 ≈ 1.55 GeV reduces Eq. (39) from 7.03 × 10⁻¹³ to about 4.7 × 10⁻¹³. The authors should either justify the scale choice explicitly or present the result as a function of m_c with a corresponding uncertainty.
- [Sec. II, paragraph after Fig. 1] The statement that diagrams (a) and (b) vanish "due to the conservation of parity" is not the correct selection rule. Parity does not forbid J/ψ → γγ*; the correct rule is charge-conjugation invariance: C(J/ψ) = -1 while a two-photon state (real plus virtual) has C = +1, so the amplitude J/ψ → γγ* is exactly zero. Because this step is the basis for keeping only the final-state-radiation diagrams, the explanation should be corrected and the asserted explicit verification should either be shown or replaced by a precise reference.
minor comments (6)
- [Sec. IV, after Eq. (32)] The sentence claiming that Eq. (32) was verified "in another approach" is not substantiated by any calculation; please include the cross-check or explicitly state that the derivation is omitted.
- [Sec. V, event yield estimate] The estimate of 2-3 events per year at STCF assumes 100% detection efficiency and does not include the reconstruction efficiency for the γγ decay mode or the associated acceptance; the event yield should be presented as a raw rate with these efficiencies clearly separated.
- [Eq. (2)] The reduction from the trace over the spin projector to the second line of Eq. (2) is not shown; a brief derivation or a specific reference would make the amplitude easier to verify.
- [Sec. IV, Eq. (24)] The symbol Γ is overloaded: it denotes the muon decay width in the Coulomb Green function and also the J/ψ decay width in the ratio R; please use distinct symbols such as Γ_μ.
- [Sec. IV, paragraph on muon lifetime] The statement that the lifetime of true muonium is "on the order of several picoseconds" is imprecise for the 1¹S₀ ground state, whose lifetime is about 0.6 ps; please specify the state and distinguish para from ortho.
- [Reference [11]] Reference [11] appears malformed; the correct citation is D. B. Cassidy and A. P. Mills, Jr., Nature 449, 193 (2007).
Circularity Check
No material circularity: the bound-state branching fraction is a parameter-free QED ratio normalized to measured J/psi->mu+mu-, with only contextual self-citation.
full rationale
The central prediction Br(J/psi -> gamma(mu+mu-)_bound) ~ 7e-13 is obtained as R|E<0 times the measured Br(J/psi -> mu+mu-), where R is computed from QED amplitudes and Coulomb resummation. The J/psi radial wave function R(0) appears in both the numerator and denominator of R and cancels (Eqs. (28)-(32)); no fitted parameter is introduced. The normalization by the measured leptonic branching fraction is a standard external input, not a circular reuse of the quantity being predicted. The only self-citation (Ref. [13]) states that true muonium has not yet been observed experimentally and plays no role in the derivation. The parity/C-parity suppression of the photon-from-charmonium diagrams is a selection rule external to the model, and the m_c dependence is an ordinary parametric uncertainty (R is proportional to 1/m_c^2), not an ansatz fitted to the target rate. The derivation is self-contained QED plus standard NRQCD projector inputs, with the numerical result benchmarked against measured Br(J/psi -> mu+mu-) only as a normalization.
Assumptions & free parameters
free parameters (2)
- m_c (charm quark mass) =
1.27 GeV
- Lambda (energy window for the E>0 continuum) =
alpha^2 m_mu/4 ≈ 1.4 keV, later MeV-scale
assumptions (3)
- domain assumption J/psi is a nonrelativistic ccbar system with M_J/psi ≈ 2 m_c
- ad hoc to paper The diagrams with the photon emitted from the charmonium vanish by parity
- standard math Coulomb rescattering of the muon pair is resummed with the Fadin-Khoze Green function formula
Cite this review
Pith. "Pith review of Creating true muonium via charmonium radiative decay." pith.science (2026). https://pith.science/paper/LR7P6M3F
@misc{pith2026241212592,
author = {Pith},
title = {Pith review of: Creating true muonium via charmonium radiative decay},
year = {2026},
howpublished = {\url{https://pith.science/paper/LR7P6M3F}},
note = {Machine review of arXiv:2412.12592}
}
abstract
True muonium, the bound state of a muon and an antimuon ($\mu^+\mu^-$), has long been theoretically predicted but remains experimentally elusive. We investigate the production of true para-muonium in the radiative decay of $J/\psi$ meson,and analyze the prospects for detecting true muonium in current and future high-energy $e^+e^-$ experiments, particularly focusing on the BESIII experiment and the proposed Super Tau-Charm Facility. Although the events are rare at the super tau-charm facility, the detection of true para-muonium via $J/\psi$ radiative decays could become feasible at its future updates.
Figures
Reference graph
Works this paper leans on
-
[13]
M. A. J. Cassidy DB, Nature 7159, 195 (2007)
work page 2007
-
[1]
In the process, the photon can be radi- ated either from the J/ψ meson or from outgoing muons (µ+ and µ−). However, due to the conservation of par- ity, the first two diagrams, where the photon is radiated directly from the J/ψ meson, yield zero contribution, a result that can be rigorously verified through explicit cal- culations. Consequently, the nonzero...
-
[2]
C. M. G. Lattes, H. Muirhead, G. P. S. Occhialini, and C. F. Powell, Nature 159, 694 (1947)
work page 1947
-
[3]
R. E. Marshak and H. A. Bethe, Phys. Rev. 72, 506 (1947)
work page 1947
-
[4]
C. M. G. Lattes, G. P. S. Occhialini, and C. F. Powell, Nature 160, 486 (1947)
work page 1947
-
[5]
C. M. G. Lattes, G. P. S. Occhialini, and C. F. Powell, Nature 160, 453 (1947)
work page 1947
-
[6]
Deutsch, Phys
M. Deutsch, Phys. Rev. 82, 455 (1951)
1951
-
[7]
Mohoroviˇ ci´ c, Astronomische Nachrichten 253, 93 (1934)
S. Mohoroviˇ ci´ c, Astronomische Nachrichten 253, 93 (1934)
work page 1934
Show all 50 references
-
[8]
L. N. S. Bilen’kii, N. van Hieu and F. Tkebuchava, Sov.J 10, 469 (1969)
1969
-
[9]
V. W. Hughes, D. W. McColm, K. Ziock, and R. Prepost, Phys. Rev. Lett. 5, 63 (1960)
1960
-
[10]
Coombes et al., Phys
R. Coombes et al., Phys. Rev. Lett. 37, 249 (1976)
1976
-
[11]
Hughes and B
V. Hughes and B. B. Maglic, Am. Phys. Soc. 16, 69 (1971)
1971
-
[12]
S. G. Karshenboim, V. G. Ivanov, U. D. Jentschura, and G. Soff, J. Exp. Theor. Phys. 86, 226 (1998)
1998
-
[14]
Francener, V
R. Francener, V. P. Goncalves, B. D. Moreira, and K. A. Santos, Phys. Lett. B 854, 138753 (2024), 2404.11610
2024 arXiv
- [15]
-
[16]
Francener, V
R. Francener, V. P. Goncalves, and B. D. Moreira, Eur. Phys. J. A 58, 35 (2022), 2110.03466
2022 arXiv
-
[17]
L. N. S. Bilenikii, N. van Hieu and F. Tkebuchava, Sov.J.Nucl.Phys. 10, 469 (1969)
1969
-
[18]
P. A. Krachkov and A. I. Milstein, Nucl. Phys. A 971, 71 (2018), 1712.09770
2018 arXiv
-
[19]
Arteaga-Romero, C
N. Arteaga-Romero, C. Carimalo, and V. G. Serbo, Phys. Rev. A 62, 032501 (2000), hep-ph/0001278
2000 arXiv
-
[20]
I. F. Ginzburg, U. D. Jentschura, S. G. Karshenboim, F. Krauss, V. G. Serbo, and G. Soff, Phys. Rev. C 58, 3565 (1998), hep-ph/9805375
1998 arXiv
- [21]
-
[22]
Yu and Y.-D
G.-M. Yu and Y.-D. Li, Chin. Phys. Lett. 30, 011201 (2013)
2013
- [23]
-
[24]
G. Yu, Z. Zhao, Y. Cai, Q. Gao, Q. Hu, and H. Yang (2022), 2209.11439
2022 arXiv
-
[25]
Azevedo, V
C. Azevedo, V. P. Gon¸ calves, and B. D. Moreira, Phys. Rev. C 101, 024914 (2020), 1911.10861
2020 arXiv
-
[26]
Nemenov, Sov.J.Nucl.Phys
L. Nemenov, Sov.J.Nucl.Phys. 15, 582 (1972)
1972
-
[27]
B. M. V.W. Hughes, Bull.Am.Phys.Soc. 16, 65 (1971)
1971
-
[28]
J. W. Moffat, Phys. Rev. Lett. 35, 1605 (1975)
1975
-
[29]
Cid Vidal, P
X. Cid Vidal, P. Ilten, J. Plews, B. Shuve, and Y. Soreq, Phys. Rev. D 100, 053003 (2019), 1904.08458
2019 arXiv
-
[30]
Gargiulo, S
R. Gargiulo, S. Palmisano, E. Di Meco, E. Diociaiuti, I. Sarra, and D. Paesani, J. Phys. G 51, 045004 (2024), 2309.11683
2024 arXiv
-
[31]
S. J. Brodsky and R. F. Lebed, Phys. Rev. Lett. 102, 213401 (2009), 0904.2225
2009 arXiv
- [32]
- [33]
-
[34]
V. S. Fadin, V. A. Khoze, and T. Sjostrand, Z. Phys. C 48, 613 (1990)
1990
- [35]
- [36]
-
[37]
Bogomyagkov, V
A. Bogomyagkov, V. Druzhinin, E. Levichev, A. Milstein, and S. Sinyatkin, EPJ Web Conf. 181, 01032 (2018), 1708.05819
2018 arXiv
- [38]
-
[39]
Ablikim et al
M. Ablikim et al. (BESIII), Phys. Rev. D 109, 052006 (2024), 2111.13881
2024
-
[40]
Sommerfeld, Annalen Phys
A. Sommerfeld, Annalen Phys. 403, 257 (1931)
1931
-
[41]
Petrelli, M
A. Petrelli, M. Cacciari, M. Greco, F. Maltoni, and M. L. Mangano, Nucl. Phys. B 514, 245 (1998), hep- ph/9707223
1998
-
[42]
J. S. Schwinger, PARTICLES, SOURCES AND FIELDS. VOLUME II (1973)
1973
-
[43]
A. D. Sakharov, Zh. Eksp. Teor. Fiz. 18, 631 (1948)
1948
-
[44]
Navas et al
S. Navas et al. (Particle Data Group), Phys. Rev. D 110, 030001 (2024)
2024
-
[45]
V. S. Fadin and V. A. Khoze, JETP Lett. 46, 525 (1987)
1987
-
[46]
Ablikim et al
M. Ablikim et al. (BESIII), Chin. Phys. C 48, 093001 (2024), 2403.06766
2024
- [47]
-
[48]
Ablikim et al
M. Ablikim et al. (BESIII), Nucl. Instrum. Meth. A 614, 345 (2010), 0911.4960
2010 arXiv
-
[49]
Achasov et al., Front
M. Achasov et al., Front. Phys. (Beijing) 19, 14701 (2024), 2303.15790
2024
-
[137]
(37) Using the branching fraction of J/ψ → µ+µ− [44] as follows: Br(J/ψ →µ+µ−) = 5.961%, (38) we can then calculate the branching fraction of J/ψ → γ(µ+µ−)
(36) With these inputs, we get the ratio R for the bound state contributions R|E<0 ≈ 1.18 × 10−11. (37) Using the branching fraction of J/ψ → µ+µ− [44] as follows: Br(J/ψ →µ+µ−) = 5.961%, (38) we can then calculate the branching fraction of J/ψ → γ(µ+µ−). The result reads Br(J...
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.