{"id":"c6521a26-ca72-485d-a83f-5580bbb3837e","arxiv_id":"2505.10070","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"The paper predicts effective true muonium production cross sections of 1.23 μb in 200 GeV AuAu and 14.2 μb in 5.02 TeV PbPb collisions, corresponding to O(10^4) and O(10^5) atoms per billion central events.","lead":"True muonium, an unobserved atom made of a muon and an antimuon, might be produced in the hot quark-gluon plasma of heavy-ion collisions. The paper's transport calculation predicts tens of thousands of these atoms per billion collisions, which would make RHIC and the LHC discovery machines for this state.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 512-fm-radius true muonium cannot be assumed to survive the QGP: Debye screening and Coulomb breakup invalidate the gain-only Boltzmann equation, and the time-scale argument in the text has the wrong implication.","rationale":"The reader's weakest_assumption identifies the same issue I find most load-bearing. The production amplitudes and the DQPM/MUSIC transport framework are standard, and the numerical implementation appears self-consistent as a formal calculation of vacuum bound-state production in an expanding medium. However, the discovery claim requires those vacuum bound states to survive inside the QGP. The paper's own time-scale comparison, τ_b≫τ_QGP and r_b≫r_QGP, is presented as if it proved decoupling, but it proves the opposite: a state with radius 512 fm and binding energy 1.4 keV is more fragile, not less, in a plasma with screening length about 1 fm. The only medium effect checked is the thermal-mass forbiddenness of annihilation to q qbar; Coulomb ionization by thermal partons is not estimated despite being the dominant destruction mechanism for a weakly bound Coulomb state. Because Eq. (5) is gain-only, even a modest loss rate would reduce the integrated yield by orders of magnitude. The proposed Yukawa bound-state check settles the concern without requiring a full simulation. I therefore keep the reader's REJECT verdict, with no adjustment needed.","tokens_in":8804,"tokens_out":9183,"duration_ms":97568,"concrete_test":"Solve the bound-state problem for the Debye-screened Coulomb potential V(r)=-(α/r)e^{-r/r_D} with r_D=2.11/T at T=320 MeV and T=450 MeV and a_0=512 fm. If the 1s state is absent (or its |ψ(0)|^2 is suppressed relative to vacuum by more than a factor of 10), Eq. (3) cannot be used inside the QGP and the yields collapse; if it binds, add the Coulomb-breakup loss term n_q σ_ion v_rel to Eq. (5) and check whether the integrated yield survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—O(10^4–10^5) true muonium per billion central events—rests on Eq. (5), where all dissociation is dropped as 'safe.' The only medium process checked is annihilation to q qbar, forbidden by thermal quark masses; Coulomb breakup is never estimated. The justification offered, 'because τ_b≫τ_QGP and r_b≫r_QGP, the QGP does not interfere,' has the wrong implication. From the paper's own numbers, r_b=512n fm and τ_b=1.7×10^-21 s, while r_D≈2.11/T≈1 fm and τ_QGP≈10^-23 s. A Yukawa potential with r_D/a_0≈0.002 has no 1s bound state (the critical screening length is ≈0.84 a_0), so the vacuum Coulomb vertex of Eq. (3), obtained from free-photon resummation in Fig. 2, is not the in-medium formation amplitude: the exchanged photon carries Debye mass m_D≈gT and cannot resolve a 512-fm bound state. Even if a state formed, a single thermal parton scattering with momentum transfer above the 0.4 MeV ionization threshold would ionize it; the geometric cross section π(512 fm)^2≈8×10^5 fm^2 implies a dissociation rate far exceeding the QGP expansion rate. Eq. (5) is therefore gain-only in a regime where the loss term dominates, and the quoted cross sections 1.23/14.2 μb are unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8993,"tokens_out":6973,"duration_ms":77972,"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":[{"comment":"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.","section":"Production in heavy ion collisions, Eq. (5)"},{"comment":"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.","section":"Characteristic Scales"}],"minor_comments":[{"comment":"The abstract says 'per billion AA collisions' while the text says 'per billion central collisions'; the centrality condition should be stated consistently.","section":"Abstract"},{"comment":"'Six decades smaller' should read 'six orders of magnitude'.","section":"Introduction"},{"comment":"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.","section":"Fig. 3 caption"},{"comment":"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.","section":"Notation"},{"comment":"Reference [37] is an arXiv preprint; if a published journal version exists, it should be cited instead.","section":"References"}],"recommendation":"reject","confidential_remarks":"The paper's main quantitative conclusion depends on the physically implausible claim that a 512 fm-sized bound state can be produced inside a QGP without being dissociated. Fixing this would require a substantially different calculation—either a proper treatment of in-medium bound-state formation and dissociation, or a two-step scenario in which the state forms after freeze-out and the production amplitude is derived accordingly. This goes beyond a local revision of the present Letter, hence my recommendation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper gives a concrete, falsifiable estimate for true muonium production in heavy-ion collisions. That is genuinely new, and I want to give it credit before the criticism. The authors compute NLO amplitudes for q qbar -> (mu+mu-) + g and q g -> (mu+mu-) + q with DQPM thermal masses, use the Coulomb-resummed vertex for the bound state, feed the production rates into an analytically solved Boltzmann equation, and couple it to MUSIC hydrodynamics. The numbers are clear: effective cross sections of 1.23 micro-barn at RHIC and 14.2 micro-barn at LHC, with O(10^4) and O(10^5) true muonium per billion central events. The step beyond Chen and Zhuang is real - three flavors, temperature-dependent coupling, viscous hydro - and the transport solution is elegant.\n\nThe trouble is the survival assumption, and it is load-bearing, not a minor detail. The authors drop dissociation from the Boltzmann equation because, they say, one can safely ignore it (\"small phase space distribution and relatively long lifetime\"). They only check annihilation to q qbar, which is blocked by thermal quark masses. They never estimate Coulomb breakup by thermal partons. And the numbers they quote actually make the problem worse: for n=1, r_b = 512 fm while the Debye length is about 1 fm, and the binding energy is 1.4 keV. A Yukawa potential with m_D a0 ~ 600 has no 1s bound state, so the vacuum Coulomb vertex used in Eq. (3) is not the in-medium formation amplitude. Even if a state somehow formed, the geometric breakup cross section is about 8 x 10^5 fm^2, which with QGP densities gives a dissociation rate far above the plasma expansion rate. Their time-scale argument, tau_b >> tau_QGP and r_b >> r_QGP, has the wrong implication: a state much larger than the plasma and with a formation time longer than the plasma lifetime is exactly what the plasma will destroy. The statement that the binding energy shift is only 1.5 eV at T = 500 MeV is inconsistent with the same Debye mass.\n\nThe production machinery itself appears sound; the error is isolated and testable. But the central claim - order 10^4-10^5 true muonium per billion events - currently rests on a single untested assumption that fails a quick sanity check. I would not cite these yields, and the discovery claim is not justified. Still, the paper deserves a serious referee: the question matters, the calculation is clean, and an expert referee could either pin down the dissociation rate or send the authors to include a loss term. For a reading group, it is a good cautionary example of how in-medium bound states can bite, but I would not put it high on the list.\n\nRecommendation: send to peer review with the loss term and the in-medium vertex flagged as the things to fix.","headline":"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.","tokens_in":9699,"tokens_out":7702,"would_cite":false,"duration_ms":77688,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","36.10.-k","12.20.-m"],"model":"deepseek-v4-flash","headline":"Heavy-ion collisions could make true muonium detectable for the first time.","keywords":["true muonium","quark-gluon plasma","heavy-ion collisions","Boltzmann transport","viscous hydrodynamics","QED bound states","thermal production","RHIC LHC"],"falsifier":"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.","tokens_in":8437,"feed_emoji":"⚛️","tokens_out":12345,"duration_ms":112614,"temperature":0.7,"pith_summary":"True muonium—the hydrogen-like atom made of a muon and an antimuon—has never been observed. This paper argues that the quark-gluon plasma created in relativistic heavy-ion collisions should produce it in detectable numbers, through quark-antiquark and quark-gluon fusion with an extra gluon in the final state. Using thermal dressed quark masses and a viscous hydrodynamic description of the expanding fireball, the authors obtain effective production cross sections of 1.23 $\\mu$b at RHIC and 14.2 $\\mu$b at the LHC, with yields of roughly $10^4$ and $10^5$ per billion central collisions. That is more than a thousand times higher than an earlier estimate, and it would make heavy-ion collisions a viable discovery channel for the first observation of true muonium. The same signal would also act as a sensitive probe of the early-time temperature of the quark-gluon plasma.","feed_headline":"True muonium production jumps 1,000-fold in heavy-ion collisions","feed_subtitle":"A narrow resonance near 211 MeV may reveal the never-seen muonium atom in RHIC and LHC data.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"establishes the true muonium lifetime and dominant decay to e+e-, used to argue the QGP does not interfere with decay.","marker":"[5]"},{"why":"supplies the n^3S1-to-n^1S0 event ratio used to rescale ultra-peripheral cross sections for comparison.","marker":"[26]"},{"why":"supplies the n^1S0 production cross sections in ultra-peripheral collisions that the QGP result is compared against.","marker":"[27]"},{"why":"provides the earlier QGP production calculation whose approximately 10^3 times smaller results this paper improves on.","marker":"[37]"},{"why":"provides the thermal dressed quark and gluon masses and temperature-dependent coupling used in the amplitudes and thermal distributions.","marker":"[38]"},{"why":"gives the Coulomb-resummed photon-bound-state vertex, equivalent to the wavefunction at the origin.","marker":"[39, 40]"},{"why":"supplies the analytic solution of the relativistic Boltzmann equation used to integrate yields over the expanding medium.","marker":"[42]"},{"why":"provides the viscous hydrodynamic simulation of the quark-gluon plasma expansion.","marker":"[43, 44]"},{"why":"supplies the lattice QCD equation of state and shear viscosity used in the hydrodynamic evolution.","marker":"[45]"}],"fun_headline_variants":["Heavy-ion collisions could finally reveal true muonium","True muonium yield jumps 1000-fold in heavy-ion collisions","Quark-gluon plasma could spawn true muonium in heavy-ion collisions","Heavy-ion collisions may yield true muonium at last","A new path to true muonium: heavy-ion collisions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Heavy-ion collisions could finally reveal true muonium","True muonium yield jumps 1000-fold in heavy-ion collisions","Quark-gluon plasma could spawn true muonium in heavy-ion collisions","Heavy-ion collisions may yield true muonium at last","A new path to true muonium: heavy-ion collisions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000931,"raw_usage":{"total_tokens":3981,"prompt_tokens":936,"completion_tokens":3045,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":552,"completion_tokens_details":{"reasoning_tokens":2959}},"tokens_in":552,"tokens_out":3045,"duration_ms":21756,"temperature":1.0,"reasoning_tokens":2959,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:18:27.260653+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Decay Rate and Hyperfine Structure of the bound mu+ mu- system","cited_arxiv_id":"hep-ph/9706401","evidence_quote":"establishes the true muonium lifetime and dominant decay to e+e-, used to argue the QGP does not interfere with decay."}],"review_version":1}