REVIEW 2 major objections 5 minor 48 references
Searching for True Muonium in Relativistic Heavy Ion Collisions
T0 review · 2 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Heavy-ion collisions could make true muonium detectable for the first time.
desk verdict A careful production calculation whose central numbers rest on an overturned survival assumption: true muonium with a 512-fm Bohr radius will not simply sail through the QGP. 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 object doing the work is the $n^3S_1$ true-muonium bound state, a nonrelativistic QED atom with binding energy $-1.4/n^2$ keV, mass $2m_\mu=211$ MeV, and Bohr radius roughly $512n$ fm. Its production is encoded in a Coulomb-resummed vertex $iZ_n=i(\alpha^2/4\pi)\sqrt{m_\mu^3 m_n/n^3}$, equivalent to evaluating the bound-state wavefunction at zero separation. The two NLO amplitudes (Eq. (4) and its crossing) are integrated with thermal distributions of dressed quarks, antiquarks, and gluons, and the resulting collision terms enter an analytic solution of the relativistic Boltzmann equation (Eq. (7)) over a viscous hydrodynamic background; the sum over principal quantum numbers factorizes into $\zeta_3=\sum_n n^{-3}$. The step function $\Theta(T-T_c)$ restricts production to the deconfined phase.
What would settle it
Compute the thermally averaged Coulomb-dissociation cross section of the $1.4$ keV-bound state by scattering thermal gluons and quarks at $T\sim300$–$500$ MeV and add it as a loss term in Eq. (5); if the surviving yield falls far below the quoted $\mathcal{O}(10^4)$–$\mathcal{O}(10^5)$, the discovery claim fails. A targeted experimental check would be a high-statistics search for a narrow peak near $211$ MeV in the dimuon or dielectron invariant-mass spectrum of central Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.
Extended reading notes
Core claim
The paper's central claim is that the quark-gluon plasma acts as a thermal source of the $n^3S_1$ state of true muonium, with production rates large enough for discovery at RHIC and the LHC. The calculation includes two next-to-leading-order processes, $q\bar q\to(\mu^+\mu^-)g$ and $qg\to(\mu^+\mu^-)q$, with the virtual photon coupled to the bound state through a Coulomb-resummed bubble amplitude, and it follows the produced states through an expanding viscous fluid using an analytic Boltzmann transport solution. The result is an effective cross section of $1.23\,\mu$b for central Au+Au at $\sqrt{s_{NN}}=200$ GeV and $14.2\,\mu$b for Pb+Pb at $\sqrt{s_{NN}}=5.02$ TeV, corresponding to $\mathcal{O}(10^4)$ and $\mathcal{O}(10^5)$ true muonium per billion central collisions. The authors take these values, over $10^3$ times the earlier estimate, as sufficient to allow a discovery of $(\mu^+\mu^-)$ and to make true muonium a new probe of the early quark-gluon plasma.
Load-bearing premise
The calculation assumes true muonium, once formed in the plasma, survives long enough to be counted; the Boltzmann equation drops dissociation and only the annihilation channel is ruled out.
Editorial extensions
If this is right
- A narrow resonance near $211$ MeV should appear in the $\mu^+\mu^-$ (or $e^+e^-$) invariant-mass spectrum of central heavy-ion collisions at RHIC and the LHC.
- The predicted di-lepton yield from true muonium lies between that of $J/\psi$ and $\Upsilon$, within reach of existing heavy-ion detectors.
- Because production is concentrated at early times and the yield changes by a factor of two for a 30 MeV change in temperature, the signal is a sensitive early-QGP thermometer.
- QGP production exceeds ultra-peripheral production by about two orders of magnitude, making heavy-ion collisions the most promising route to first observation.
- The effective cross sections, $1.23\,\mu$b at RHIC and $14.2\,\mu$b at the LHC, are comparable to $\Upsilon$ production.
Reading between the lines
- If Coulomb dissociation of the loosely bound state were included as a loss term, the net yield could drop by orders of magnitude; a quantitative breakup calculation would directly test the paper's central assumption.
- The same thermal-production machinery could be applied to other shallow QED bound states that might form inside a quark-gluon plasma, such as pionium.
- Because the predicted peak sits near $211$ MeV with little background, existing high-statistics dimuon data sets from RHIC and the LHC could be reanalyzed for the signal without new runs.
- The steep temperature dependence suggests true muonium could be combined with direct-photon spectra to separate early-time and late-time temperature measurements of the plasma.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes that true muonium (μ⁺μ⁻) can be produced thermally in the quark-gluon plasma formed in relativistic heavy-ion collisions. The authors compute the NLO scattering amplitudes for q q̄ → (μ⁺μ⁻)g and qg → (μ⁺μ⁻)q using DQPM thermal masses, and couple a gain-only Boltzmann equation to viscous hydrodynamic simulations of the QGP. They report effective production cross sections of 1.23 μb at RHIC and 14.2 μb at the LHC, corresponding to yields of O(10⁴) and O(10⁵) per billion central collisions, and conclude that heavy-ion collisions are a promising discovery channel for true muonium.
Significance. If the yield claim were correct, the paper would establish a new production and detection channel for an as-yet-unobserved QED bound state and would introduce a new early-time thermometer for the QGP. The transport framework, the analytic solution of Eq. (7), and the coupling to realistic hydrodynamic evolution are useful features, and the scattering amplitudes follow standard Feynman rules. However, the central quantitative claim rests on the unsupported neglect of dissociation in Eq. (5). Since true muonium is an extremely fragile object with a 512 fm Bohr radius and 1.4 keV binding energy, the absence of any estimate of Coulomb breakup by thermal partons is a load-bearing gap that undermines the reported yields and cross sections.
major comments (2)
- [Production in heavy ion collisions, Eq. (5)] The gain-only Boltzmann equation is not justified. The text asserts that dissociation can be safely ignored because of the small phase-space distribution and relatively long lifetime, but neither property controls the breakup rate in a hot medium. For a state with Bohr radius r_b ≈ 512 fm, the geometric ionization cross section by a scattering thermal parton is σ ≈ π r_b² ≈ 8×10⁵ fm², and with typical QGP parton densities the resulting dissociation rate is many orders of magnitude larger than the QGP expansion rate. Ruling out only the annihilation channel (μ⁺μ⁻) → q q̄ via thermal quark masses does not address Coulomb breakup. Without a loss term in Eq. (5), the yields computed from Eqs. (7)–(9) are unsupported and could be wrong by orders of magnitude.
- [Characteristic Scales] The time-scale argument has the wrong implication. From the paper's own numbers, τ_b ≈ 1.7×10⁻²¹ s, r_b ≈ 512 fm, τ_QGP ≈ 10⁻²³ s, and r_D ≈ 2.11/T ≈ 1 fm. If τ_b ≫ τ_QGP and r_b ≫ r_QGP, then the pair cannot form a bound state inside the QGP; if it binds only after freeze-out, the in-medium collision term in Eq. (6) with the vacuum Coulomb-resummed vertex of Eq. (3) is not the correct production amplitude. Moreover, a Yukawa potential with r_D/r_b ≈ 0.002 has no 1s bound state, so the vertex of Eq. (3), obtained from free-photon resummation in Fig. 2, is not the in-medium formation amplitude. The quoted 1.5 eV modification of the binding energy is not the relevant criterion: the medium modifies the wave function and the production vertex itself.
minor comments (5)
- [Abstract] The abstract says 'per billion AA collisions' while the text says 'per billion central collisions'; the centrality condition should be stated consistently.
- [Introduction] 'Six decades smaller' should read 'six orders of magnitude'.
- [Fig. 3 caption] The caption assigns T = 500, 300, 200 MeV to the orange, green, and blue curves, but the legend in the figure itself does not show this color–temperature correspondence; please add the temperatures to the legend.
- [Notation] The notation n³S₁ should be explicitly defined as the n-th ³S₁ state of the μ⁺μ⁻ system; the symbols n, m_n, and ζ₃ are used later without a clear definition of the principal quantum number summation range.
- [References] Reference [37] is an arXiv preprint; if a published journal version exists, it should be cited instead.
Circularity Check
No material circularity: the true muonium yield is a genuine model output; the few self-citations are technical tools, not load-bearing inputs.
full rationale
The paper's central output—effective cross sections of 1.23 μb at RHIC and 14.2 μb at LHC, and yields O(10^4) and O(10^5) per billion central collisions—is obtained by integrating NLO QED/QCD scattering amplitudes (Eq. (4)) with DQPM thermal masses and couplings that are calibrated to lattice entropy from Ref. [38], and by transporting those rates through a MUSIC hydrodynamic medium whose initial conditions are fixed by measured charged-particle multiplicities. These inputs are external to the present work and do not encode the target yield. The only places where prior work by the same authors enters are the analytic Boltzmann solution in Eq. (7), cited to Refs. [37,42], and the comparison baseline Ref. [37]; these are mathematical or reference tools, not assumptions that define the predicted yield. The decision to drop dissociation in Eq. (5), while physically questionable for a 512-fm-radius state in a Debye-screened plasma, is a modeling approximation rather than a circular reduction: no dissociation parameter is fitted and then relabeled as a prediction, and omitting a loss term cannot make the gain-only result equal to any input by construction. No equation reduces to another by definition, no fitted parameter is renamed as a prediction, and no uniqueness theorem is imported from the authors' own prior work. The self-citations are ancillary to the main derivation, so the score is 1 rather than 0 only to reflect that the transport solution is reused from the authors' earlier framework without independent re-derivation; it remains standard and checkable.
Assumptions & free parameters
free parameters (6)
- lambda in g^2(T) =
2.42
- T* in g^2(T) =
1.19 T_c
- low-T power in g^2(T) =
3.1
- strange quark mass offset =
0.045 GeV
- shear viscosity to entropy ratio =
0.08
- initial proper time and temperature profile =
tau0=0.6 fm/c, Tmax=320 MeV (RHIC); tau0=0.4 fm/c, Tmax=450 MeV (LHC)
assumptions (5)
- domain assumption The QGP does not modify the muon-antimuon bound-state wavefunction, so the vacuum Coulomb-resummed vertex iZ_n is valid.
- domain assumption Thermal quark masses from DQPM are large enough that the leading-order q qbar to muon pair process is kinematically forbidden, leaving only NLO channels.
- ad hoc to paper Dissociation of true muonium inside the QGP is negligible, so the Boltzmann equation has only gain terms.
- domain assumption Quarks, antiquarks, and gluons in the QGP follow thermal distributions with DQPM dressed masses, and the emitted gluon has only transverse polarizations.
- domain assumption The QGP spacetime evolution is described by MUSIC with eta/s=0.08, zero bulk viscosity, and the Bernhard lattice EoS, with initial conditions tuned to final charged-hadron multiplicity.
Cite this review
Pith. "Pith review of Searching for True Muonium in Relativistic Heavy Ion Collisions." pith.science (2026). https://pith.science/paper/PSX2346G
@misc{pith2026250510070,
author = {Pith},
title = {Pith review of: Searching for True Muonium in Relativistic Heavy Ion Collisions},
year = {2026},
howpublished = {\url{https://pith.science/paper/PSX2346G}},
note = {Machine review of arXiv:2505.10070}
}
abstract
We investigate the production of the as-yet-undetected true muonium within the quark-gluon plasma formed in relativistic heavy-ion collisions, employing a relativistic Boltzmann transport framework coupled to viscous hydrodynamic simulations. The obtained effective cross sections for central collisions are 1.23~$\mu b$ in AuAu collisions with $\sqrt{s_{\rm NN}}=200$~GeV and 14.2~$\mu b$ in PbPb collisions with $\sqrt{s_{\rm NN}}=5.02$~TeV, resulting in a yield of $\mathcal{O}(10^4)$ and $\mathcal{O}(10^5)$ true muonium per billion $AA$ collisions at RHIC and the LHC, respectively. This establishes heavy-ion collisions as a promising process for detecting true muonium.
Figures
Reference graph
Works this paper leans on
-
[1]
S. M. Bilenky, V. K. Nguyen, L. L. Nemenov and F. G. Tkebuchava, SLAC-TRANS-0097
-
[2]
S. M. Bilenky, V. H. Nguyen, L. L. Nemenov and F. G. Tkebuchava, Yad. Fiz.10, 812-814 (1969)
work page 1969
-
[3]
J. E. Malenfant, UMI-87-11690
-
[4]
U. D. Jentschura, G. Soff, V. G. Ivanov and S. G. Karshenboim, Phys. Rev. A56, 4483 (1997) [arXiv:physics/9706026 [physics]]
work page Pith review arXiv 1997
-
[5]
U. D. Jentschura, V. G. Ivanov, G. Soff and S. G. Karshenboim, Phys. Lett. B424, 397-404 (1998) [arXiv:hep-ph/9706401 [hep-ph]]
work page Pith review arXiv 1998
-
[6]
Y. Ji and H. Lamm, Phys. Rev. A94, no.3, 032507 (2016) [arXiv:1607.07059 [physics.atom-ph]]
work page Pith review arXiv 2016
-
[7]
Muonic hydrogen and MeV forces
D. Tucker-Smith and I. Yavin, Phys. Rev. D83, 101702 (2011) [arXiv:1011.4922 [hep-ph]]
work page Pith review arXiv 2011
-
[8]
Applying Bayesian Inference to Galileon Solutions of the Muon Problem
H. Lamm, Phys. Rev. D94, no.11, 115007 (2016) [arXiv:1609.07520 [hep-ph]]
work page Pith review arXiv 2016
Show all 48 references
-
[9]
Crane, D
T. Crane, D. Casperson, P. Crane, P. Egan, V. W. Hughes, R. Stambaugh, P. A. Thompson and G. Zu Putlitz, Phys. Rev. Lett.27, 474-476 (1971)
1971
-
[10]
J. W. Moffat, Phys. Rev. Lett.35, 1605 (1975)
1975
-
[11]
Holvik and H
E. Holvik and H. A. Olsen, Phys. Rev. D35, 2124 (1987)
1987
-
[12]
Arteaga-Romero, C
N. Arteaga-Romero, C. Carimalo and V. G. Serbo, Phys. Rev. A62, 032501 (2000) [arXiv:hep-ph/0001278 [hep- ph]]
2000 arXiv
-
[13]
S. J. Brodsky and R. F. Lebed, Phys. Rev. Lett.102, 213401 (2009) [arXiv:0904.2225 [hep-ph]]
2009 arXiv
-
[14]
P. J. Fox, S. R. Jindariani and V. D. Shiltsev, JINST18, no.08, T08007 (2023) [arXiv:2203.07144 [hep-ex]]
2023 arXiv
-
[15]
Gargiulo, S
R. Gargiulo, S. Palmisano, E. Di Meco, E. Diociaiuti, I. Sarra and D. Paesani, J. Phys. G51, no.4, 045004 (2024) [arXiv:2309.11683 [hep-ph]]
2024 arXiv
- [16]
-
[17]
Banburski and P
A. Banburski and P. Schuster, Phys. Rev. D86, 093007 (2012) [arXiv:1206.3961 [hep-ph]]
2012 arXiv
-
[18]
Gargiulo, S
R. Gargiulo, S. Palmisano and E. Di Meco, Phys. Rev. D110, no.9, 092015 (2024) [arXiv:2409.11342 [hep-ex]]
2024 arXiv
-
[19]
L. L. Nemenov, Yad. Fiz.15, 1047-1050 (1972)
1972
-
[20]
M. I. Vysotsky, Yad. Fiz.29, 845-846 (1979)
1979
-
[21]
Malenfant, Phys
J. Malenfant, Phys. Rev. D36, 863-877 (1987)
1987
-
[22]
G. A. Kozlov, Sov. J. Nucl. Phys.48, 167-171 (1988) JINR-E2-87-557
1988
-
[23]
Ji and H
Y. Ji and H. Lamm, Phys. Rev. D98, no.5, 053008 (2018) [arXiv:1706.04986 [hep-ph]]
2018 arXiv
-
[24]
Ji and H
Y. Ji and H. Lamm, Phys. Rev. D99, no.3, 033008 (2019) [arXiv:1810.00233 [hep-ph]]
2019 arXiv
-
[25]
Cid Vidal, P
X. Cid Vidal, P. Ilten, J. Plews, B. Shuve and Y. Soreq, Phys. Rev. D100, no.5, 053003 (2019) [arXiv:1904.08458 [hep-ph]]
2019 arXiv
-
[26]
I. F. Ginzburg, U. D. Jentschura, S. G. Karshenboim, F. Krauss, V. G. Serbo and G. Soff, Phys. Rev. C58, 3565-3573 (1998) [arXiv:hep-ph/9805375 [hep-ph]]
1998 arXiv
-
[27]
Azevedo, V
C. Azevedo, V. P. Gon¸ calves and B. D. Moreira, Phys. Rev. C101, no.2, 024914 (2020) [arXiv:1911.10861 [hep- ph]]
2020 arXiv
-
[28]
G. M. Yu and Y. D. Li, Chin. Phys. Lett.30, 011201 (2013)
2013
-
[29]
J. P. Dai and S. Zhao, Phys. Rev. D109, no.5, 054022 (2024) [arXiv:2401.04681 [hep-ph]]. 6
2024 arXiv
-
[30]
The STAR Collaboration, Nature632, 1026–1031 (2024)
2024
-
[31]
The STAR Collaboration, Nature473, 353-356 (2011)
2011
-
[32]
Adams etal
The STAR Collaboration, J. Adams etal. Nucl. Phys. A 757, 102-183 (2005)
2005
-
[33]
The STAR Collaboration, Nature548, 62-65 (2017)
2017
-
[34]
Stankus, Annu
P. Stankus, Annu. Rev. Nucl. Part. Sci.55, 517 (2005)
2005
-
[35]
The NA60 Collaboration, Eur. Phys. J. C59, 607–623 (2009)
2009
-
[36]
The STAR Collaboration, [arXiv:2402.01998 [nul-ex]]
- [37]
-
[38]
Berrehrah, E
H. Berrehrah, E. Bratkovskaya, W. Cassing, P. B. Gossi- aux, J. Aichelin and M. Bleicher, Phys. Rev. C89, no.5, 054901 (2014) [arXiv:1308.5148 [hep-ph]]
2014 arXiv
-
[39]
Kong and F
X. Kong and F. Ravndal, Phys. Lett. B450, 320-324 (1999) [erratum: Phys. Lett. B458, 565-565 (1999)] [arXiv:nucl-th/9811076 [nucl-th]]
1999 arXiv
-
[40]
Kong and F
X. Kong and F. Ravndal, Nucl. Phys. A665, 137-163 (2000) [arXiv:hep-ph/9903523 [hep-ph]]
2000 arXiv
-
[41]
Braaten, E
E. Braaten, E. Johnson and H. Zhang, JHEP02, 150 (2018) [arXiv:1708.07155 [hep-ph]]
2018 arXiv
-
[42]
J. Zhao, K. Zhou, S. Chen and P. Zhuang, Prog. Part. Nucl. Phys.114, 103801 (2020) [arXiv:2005.08277 [nucl- th]]
2020 arXiv
-
[43]
Schenke, S
B. Schenke, S. Jeon and C. Gale, Phys. Rev. Lett.106, 042301 (2011) [arXiv:1009.3244 [hep-ph]]
2011 arXiv
-
[44]
Schenke, S
B. Schenke, S. Jeon and C. Gale, Phys. Rev. C82, 014903 (2010) [arXiv:1004.1408 [hep-ph]]
2010 arXiv
-
[45]
J. E. Bernhard, J. S. Moreland, S. A. Bass, J. Liu and U. Heinz, Phys. Rev. C94, no.2, 024907 (2016) [arXiv:1605.03954 [nucl-th]]
2016 arXiv
-
[46]
N. J. Abdulameeret al.[PHENIX], Phys. Rev. C109, no.4, 044912 (2024) [arXiv:2203.17187 [nucl-ex]]
2024
-
[47]
C. Shen, U. W. Heinz, J. F. Paquet and C. Gale, Phys. Rev. C89, no.4, 044910 (2014) [arXiv:1308.2440 [nucl- th]]
2014 arXiv
-
[48]
Rapp and H
R. Rapp and H. van Hees, Phys. Lett. B753, 586-590 (2016) [arXiv:1411.4612 [hep-ph]]
2016 arXiv
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.