{"id":"f85aef0b-64e6-4bb8-b2d4-ea8ce36b323a","arxiv_id":"2412.12592","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Predicts Br(J/psi -> gamma + true para-muonium) around 7e-13, implying about 2-3 events per year at the proposed Super Tau-Charm Facility.","lead":"The paper computes how often a J/psi particle decays into a photon plus a bound muon-antimuon atom, called true muonium. It finds the rate is too low for today's experiments but might give a few events per year at a proposed super charm-tau factory.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the 'parity' suppression of the photon-from-charmonium diagrams is actually a rigorous charge-conjugation selection rule, and the m_c ambiguity only shifts the branching fraction by ~50%, not by an order of magnitude.","rationale":"The reader's primary weakest assumption (the vanishing of ISR diagrams) is actually well-founded by charge conjugation, so the central approximation is secure. The secondary fragility (m_c choice) is legitimate but only introduces a ~50% shift, which is within the expected theoretical uncertainty for an order-of-magnitude prediction. The paper's internal algebra is consistent: Eq. (32) follows from Eq. (23) with the appropriate Jacobian factors, and the bound-state spin projection onto para-muonium is enforced by C-parity (ortho-muonium production is forbidden). The experimental feasibility section is clearly conditional on future detector capabilities and does not undermine the physics claim. Therefore, the reader's CONDITIONAL verdict remains appropriate without modification.","tokens_in":10411,"tokens_out":52119,"duration_ms":437160,"concrete_test":"Recompute the ratio R in Eq. (32) with m_c = 1.55 GeV (the J/ψ pole-mass half) instead of 1.27 GeV, and verify that the predicted branching fraction shifts from 7.03×10⁻¹³ to approximately 4.7×10⁻¹³, confirming the result remains at the 10⁻¹³ level. Additionally, check the vanishing of the photon-from-charmonium diagrams by a C-parity argument or by an explicit trace of the two-photon amplitude; if the amplitude is exactly zero, Eq. (39) is unaffected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's weakest assumption — that the diagrams where the photon is emitted from the J/ψ constituent charm quarks vanish — is correct, although the paper attributes it to parity instead of charge conjugation. The initial J/ψ has C = -1; the final state of a real photon (C = -1) plus a virtual photon (C = -1) has total C = +1, so the amplitude for J/ψ → γ γ* is exactly zero. Thus the FSR-only amplitude in Eq. (2) is complete and the central branching fraction of 7×10⁻¹³ is not threatened by this step. The secondary concern about m_c = 1.27 GeV is a real but modest uncertainty: replacing it with M_J/ψ/2 ≈ 1.55 GeV changes the ratio R ∝ 1/m_c² by a factor (1.27/1.55)² ≈ 0.67, giving Br ≈ 4.7×10⁻¹³. This does not alter the order-of-magnitude prediction or the conclusion that STCF could expect O(1) events per year. The experimental feasibility claim is speculative (0.4 mm neutral-vertex resolution is not yet demonstrated), but this is a stated condition, not a flaw in the central physics. I therefore find no load-bearing concern that invalidates the manuscript's main claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":10660,"tokens_out":29918,"duration_ms":269671,"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":[{"comment":"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.","section":"Sec. V, Eqs. (36)-(39)"},{"comment":"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.","section":"Sec. II, paragraph after Fig. 1"}],"minor_comments":[{"comment":"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.","section":"Sec. IV, after Eq. (32)"},{"comment":"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.","section":"Sec. V, event yield estimate"},{"comment":"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.","section":"Eq. (2)"},{"comment":"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 Γ_μ.","section":"Sec. IV, Eq. (24)"},{"comment":"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.","section":"Sec. IV, paragraph on muon lifetime"},{"comment":"Reference [11] appears malformed; the correct citation is D. B. Cassidy and A. P. Mills, Jr., Nature 449, 193 (2007).","section":"Reference [11]"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the central derivation appears sound to me; the two major issues are the unjustified charm-mass input, which changes the headline number by about a third, and the mislabeled symmetry argument for dropping the charmonium-radiation diagrams. Both are fixable within the manuscript's scope. The experimental feasibility claims should be regarded as speculative, as they assume a neutral-vertex resolution that is not currently demonstrated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a clean, internally consistent QED calculation that gives a new prediction Br(J/psi -> gamma + true muonium) ~ 7e-13, normalized to the measured Br(J/psi -> mu+mu-). The bound-state number is parameter-free, and the Coulomb resummation machinery is standard but applied to a new initial state. I think the paper deserves a serious referee.\n\nThe genuine new result is Eq. (39), and the authors do two useful things: they calculate the low-invariant-mass FSR amplitude explicitly, and they cross-check the bound-state width via an independent spin-projector/wavefunction calculation. The phase-space integral and the ratio R check out. The citation to their own earlier work [13] is contextual, not load-bearing.\n\nThe soft spots are mostly presentation and experimental extrapolation, not the physics. The paper says the photon-from-J/psi diagrams vanish by parity. That's the wrong name: it's charge-conjugation invariance (J/psi has C=-1; gamma and the virtual photon each have C=-1, so the J/psi -> gamma gamma* amplitude is C-forbidden). The conclusion is right, so the central branching fraction is not threatened. The m_c choice is a real but modest uncertainty: using M_J/psi/2 shifts the rate down by about a third, not by an order of magnitude. The experimental section is optimistic in the way this kind of paper usually is: 0.4 mm neutral-vertex resolution is not demonstrated, and 2-3 events per year at STCF is a challenge, not a discovery. The authors do state the resolution condition, so it's not hidden.\n\nOne thing I'd want the referee to push on: the Lambda-dependence of the E>0 contribution. The bound-state rate is clean, but the continuum contribution depends strongly on the energy window, and the claim that a few-MeV window could give 1e-8 is extrapolating the small-s1 approximation beyond its validity. That section could use a clearer caveat.\n\nMy verdict: send it to review. It's a solid order-of-magnitude prediction for an unobserved system, the derivation is checkable, and the main numerical result is robust at the factor-2 level. The paper's usefulness to BESIII/STCF planning is real, and the experimental discussion, while speculative, is clearly flagged as such.","headline":"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.","tokens_in":11238,"tokens_out":1891,"would_cite":false,"duration_ms":17419,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["true muonium","J/psi radiative decay","charmonium","Coulomb resummation","Sommerfeld enhancement","QED bound states","super tau-charm facility","branching fraction"],"falsifier":"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.","tokens_in":10172,"feed_emoji":"⚛️","tokens_out":9789,"duration_ms":82283,"temperature":0.7,"pith_summary":"True muonium, the bound state of a muon and an antimuon, is one of the simplest QED atoms and has never been observed. The paper shows that the radiative decay of the charmonium vector meson J/psi into a photon and a para-muonium ground state is a viable production channel, with a predicted branching fraction of about 7 x $10^{-13}$. That rate is too small for current BESIII data, but a proposed super tau-charm facility producing 3.4 x $10^{12}$ J/psi per year would yield 2-3 signal events per year. The calculation is built from final-state radiation amplitudes and a Coulomb Green function resummation, and it cancels the nonperturbative charmonium wave function by normalizing to the measured J/psi -> mu+mu- branching fraction.","feed_headline":"One in 1.4 trillion J/psi decays may create true muonium","feed_subtitle":"The predicted rate is too small for today's data, but a super tau-charm factory could catch a few events per year.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"The analogous e+e- -> (mu+mu-) + gamma process and the method of using a threshold energy window for identifying true muonium; the paper follows this strategy for J/psi decays.","marker":"[29]"},{"why":"Supplies the Coulomb-resummation relation that connects the decay width near threshold to the imaginary part of the Coulomb Green function.","marker":"[34]"},{"why":"Provides the nonrelativistic spin projector for the spin-triplet c cbar state of the J/psi, used to derive the decay amplitude.","marker":"[39]"},{"why":"Gives the explicit formula for the imaginary part of the Coulomb Green function with finite width, which yields both the bound-state poles and the continuum enhancement.","marker":"[43]"},{"why":"Supplies the input values m_mu, m_c, alpha, and the measured Br(J/psi -> mu+mu-) = 5.961% used for the numerical prediction.","marker":"[44]"},{"why":"Documents the ~10^10 accumulated J/psi events at BESIII, setting the current sensitivity limit for the search.","marker":"[45]"},{"why":"Documents the 2.7 x 10^9 psi(3686) events, used to estimate the alternative J/psi sample from psi(3686) -> pi+pi- J/psi.","marker":"[46]"},{"why":"Defines the future super tau-charm facility's expected 3.4 x 10^12 J/psi per year, the basis for the reachable event rate.","marker":"[47]"}],"fun_headline_variants":["Charmonium decay offers a path to true muonium","True muonium could be born from J/psi radiative decay","Rare decay predicts new way to create muonium","Future facility may catch 2-3 true muonium events per year","A bound muon-antimuon pair might be observed via J/psi decay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Charmonium decay offers a path to true muonium","True muonium could be born from J/psi radiative decay","Rare decay predicts new way to create muonium","Future facility may catch 2-3 true muonium events per year","A bound muon-antimuon pair might be observed via J/psi decay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000352,"raw_usage":{"total_tokens":1920,"prompt_tokens":951,"completion_tokens":969,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":880}},"tokens_in":567,"tokens_out":969,"duration_ms":8796,"temperature":1.0,"reasoning_tokens":880,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:58:17.817379+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Coulomb-resummation relation that connects the decay width near threshold to the imaginary part of the Coulomb Green function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the ~10^10 accumulated J/psi events at BESIII, setting the current sensitivity limit for the search."}],"review_version":1}