---
title: 'Dimuonium: True Muonium Bound State'
url: https://www.emergentmind.com/topics/dimuonium
type: topic
---

# Dimuonium: True Muonium Bound State

Dimuonium, also called true muonium, is the QED bound state $(\mu^+\mu^-)$, one of the pure QED bound systems of leptons together with positronium and ditauonium. It is a hydrogen-like system with a Bohr radius of about $512$ fm and a ground-state binding energy of $1.4$ keV, and it remains unobserved experimentally despite sustained theoretical attention [2508.14424] [1204.4389]. Contemporary work treats dimuonium simultaneously as a precision-QED system, a muon-sector new-physics probe, and a challenging experimental target whose production may be accessible in fixed-target reactions, low-energy $e^+e^-$ machines, heavy quarkonium decays, photon-photon collisions, and ultraperipheral heavy-ion collisions [2503.01393] [2203.07144].

## 1. Bound-state structure and quantum numbers

Dimuonium is the bound state of a muon and an antimuon held together by electromagnetic interaction. Because its constituents are leptons, it is free of nuclear-structure ambiguities that complicate atoms with composite nuclei, and its dynamics are commonly modeled with Coulombic wavefunctions and NRQED matching [2508.14424] [2203.01131].

For the ground state, the wavefunction at the origin enters essentially every production and annihilation calculation:
\[
|\Psi(0)|^2 = \frac{(\alpha m_\mu)^3}{8\pi}.
\]
This quantity appears explicitly in photon-photon production, collider production, and heavy-quarkonium decay calculations [2508.14424].

The lowest-lying spectroscopic states are the spin-singlet para-dimuonium and spin-triplet ortho-dimuonium. In the notation used across the literature, these are the $^1S_0$ and $^3S_1$ states, respectively. The $^3S_1$ state has $J^{PC}=1^{--}$ and the $^1S_0$ state has $J^{PC}=0^{-+}$ [2503.01393]. Selection rules strongly affect accessible production modes: even-spin, $C$-even states can be formed in $\gamma\gamma$ fusion, while the neutral vector $^3S_1$ state requires three photons in UPC triphoton production or other non-$\gamma\gamma$ mechanisms because of the Landau–Yang constraint [2503.10952] [2604.21838].

A compact summary of the lowest states is given below.

| State | Lifetime | Dominant decay |
|---|---:|---|
| Para-dimuonium $^1S_0$ | $595.4$ fs | $\gamma\gamma$ |
| Ortho-dimuonium $^3S_1$ | $\approx 1.81$ ps | $e^+e^-$ |

The para-state lifetime is also quoted as $\sim 0.6$ ps in production studies, while the triplet lifetime is consistently quoted near $1.81$ ps; the difference reflects the usual precision conventions across papers rather than a physics discrepancy [2503.01393] [2503.10952].

## 2. Decays, spectroscopy, and precision observables

Dimuonium decays are state-dependent and experimentally consequential. Para-dimuonium decays dominantly to $\gamma\gamma$, with branching ratio $\approx 98.5\%$ in the NLO photon-photon production study and $\sim 100\%$ in UPC phenomenology [2508.14424] [2503.10952]. For para-dimuonium in the UPC treatment, the diphoton width is
\[
\Gamma_{(\mu^+\mu^-)_0 \to \gamma\gamma} = \frac{\alpha^5 m_\mu}{2},
\]
numerically yielding approximately $1.105$ meV [2503.10952].

Ortho-dimuonium is longer-lived and phenomenologically more favorable in channels with displaced charged tracks. Its dominant decay is to $e^+e^-$, with branching $\approx 86\%$ in the photon-photon NLO analysis, and the annihilation width is written in the rare-decay study as
\[
\Gamma(\text{o-Dm}\to e^+e^-) = \frac{\alpha^5 m_\mu}{6}.
\]
This is the basis of several discovery strategies centered on low-mass displaced $e^+e^-$ vertices [2508.14424] [1706.03891].

Although dimuonium is a pure QED atom in its binding dynamics, rare hadronic channels are not absent. The decay
\[
\text{o-Dm} \rightarrow \pi^0 + \gamma
\]
has branching ratio
\[
\frac{\Gamma(\text{o-Dm}\to\pi^0\gamma)}{\Gamma(\text{o-Dm}\to e^+e^-)} \simeq 0.9 \times 10^{-5},
\]
so approximately one in $10^5$ ortho-dimuonia decays into this channel [1706.03891]. The emitted photon has energy
\[
E \simeq \frac{4m_\mu^2 - m_\pi^2}{4m_\mu} \simeq 63~\mathrm{MeV},
\]
providing a distinctive monochromatic signature unavailable in positronium because of the lower mass threshold [1706.03891].

The same virtual $\pi^0\gamma$ dynamics induces a calculable hyperfine-splitting correction. The real part of the shift is quoted as about $-6$ MHz, corresponding to a $0.3$ part-per-million effect, while the loop also reproduces the leading $\pi^0\gamma$ contribution to the muon anomalous magnetic moment,
\[
\Delta a_\mu(\pi^0\gamma) = 4.8 \times 10^{-11},
\]
as a consistency check [1706.03891]. This establishes that “pure QED” should not be misread as “completely isolated from hadronic vacuum-polarization effects”; rather, hadronic contributions are small, structured, and calculable.

## 3. Production in lepton colliders and heavy quarkonium decays

A major production axis is $\gamma\gamma$ fusion at high-luminosity $e^+e^-$ facilities. The relevant EPA process is
\[
e^+e^- \to e^+e^- + \gamma\gamma \to e^+e^- + (\mu^+\mu^-)\left[^1S_0/^3S_1\right],
\]
with total cross section obtained by convolution with the Weizsäcker–Williams photon flux [2508.14424]. The first full NLO QED calculation for this channel finds that the corrections are negative for both para and ortho states. Quantitatively, the para correction is
\[
\sigma_{\text{NLO}(^1S_0)} \approx 0.9941\, \sigma_{\text{LO}(^1S_0)},
\]
corresponding to $-0.59\%$, while the ortho correction ranges from $-1.56\%$ at BEPCII to $-1.82\%$ at CEPC [2508.14424]. The same study quotes order $89{,}000$ para-dimuonium events/year at Belle II/SuperKEKB and about $35{,}661$ para-dimuonium plus $53$ ortho-dimuonium events/year at STCF, indicating that LO-based discovery projections remain essentially intact after NLO refinement [2508.14424].

Heavy quarkonium decays provide an orthogonal route. In the NRQED/NRQCD analysis of $J/\psi$ and $\Upsilon$ decays, the most favorable dimuonium channel is inclusive ortho-dimuonium production in $J/\psi$ decay:
\[
Br(J/\psi \longrightarrow (\mu^+\mu^-)[^3S_1] +X) = 1.5 \times 10^{-12}.
\]
The corresponding inclusive para-dimuonium branching fraction is
\[
Br(J/\psi \to (\mu^+\mu^-)[^1S_0] +X) = 4.2 \times 10^{-13},
\]
and the radiative para channel is
\[
Br(J/\psi \to (\mu^+\mu^-)[^1S_0] + \gamma) = (3.93 \pm 0.13) \times 10^{-13}.
\]
With the annual production of $3.4 \times 10^{12}$ $J/\psi$ events at STCF, this translates to approximately $5$ ortho-dimuonium events/year and about $1.4$ para-dimuonium events/year under the stated assumptions [2505.02676].

Dedicated threshold machines have also been proposed. The DIMUS concept uses collisions of $408$ MeV electrons and positrons at a $75^\circ$ angle in a compact $23$ m circumference collider, targeting peak luminosity $1.6 \times 10^{32}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$ and $\sim 0.5$ million dimuons per year [2203.07144]. The large crossing angle boosts the produced atom so that the decay point is shifted from the beam-collision region; a related collider design quotes a decay length
\[
l = c\tau_{\mu\mu}\beta\gamma = c\tau_{\mu\mu}\tan\theta,
\]
with $l \approx 2$ mm at $75^\circ$, enabling strong rejection of prompt Bhabha backgrounds [1708.05819]. This suggests that direct threshold production is not only a spectroscopy strategy but also a background-engineering strategy.

## 4. Fixed-target, heavy-ion, and ultraperipheral production channels

A recent fixed-target proposal introduces a new mechanism for forming true muonium in Drell–Yan dimuon production off nuclei in a thin foil target. In this picture, a fraction of the $\mu^+\mu^-$ pairs produced in proton–nucleus Drell–Yan reactions bind via the static Coulomb interaction when produced with small relative momentum [2503.01393]. The discovery signature is an excess of $e^+e^-$ pairs with invariant mass near $2m_\mu \approx 211$ MeV and a displaced vertex downstream of the foil, with decay time consistent with the ortho-dimuonium lifetime [2503.01393]. For $4\times 10^{16}$ POT at $24$ GeV and $40$ GeV with an optimal $^9$Be foil thickness of $0.388$ mm, the surviving $oT\!M$ yield is quoted as $N_{oT\!M} \sim 993$ and $1360$, respectively; after including passage, excitation, and dissociation, the number of $oT\!M$ atoms with $E_{T\!M} > 2$ GeV surviving to vacuum is $\sim 125$ and $\sim 187$ [2503.01393]. Low-$Z$, low-density targets such as Be are identified as optimal.

UPCs furnish both para- and ortho-production channels. In exclusive photon fusion for even-spin resonances, para-dimuonium is treated as an even-spin QED bound state producible in $\gamma\gamma$ fusion. For PbPb at $5.5$ TeV, one study quotes a para-dimuonium cross section of $1.36\,\mu$b and an expected yield of $14{,}000$ events, while FCC PbPb at $39.4$ TeV yields $2.97\,\mu$b and $330{,}000$ events [2503.10952]. A separate LHC study using a direct bound-state projection formalism gives
\[
\sigma_{PM}(1S) = 2.32\, \mu\text{b}, \qquad \sigma_{PM}^{\text{total}} = 2.86\, \mu\text{b}
\]
for Pb+Pb at $\sqrt{s_{NN}} = 5.5$ TeV, and finds that EPA overestimates production at low $P_T$ by $\sim 15\%$ [2307.12387]. These two predictions are model-dependent but numerically compatible at the order-of-magnitude level.

For the vector state, triphoton fusion has emerged as a distinct UPC mechanism:
\[
\gamma + \gamma + \gamma \to (\mu^+\mu^-)[^3S_1].
\]
In Pb+Pb UPCs at $\sqrt{S}=5.02$ TeV, the predicted dimuonium cross section via $3\gamma$ fusion is $2.25\,\mu{\rm b}$ for $x > 10^{-5}$, compared with $1.5\,{\rm nb}$ for the radiative two-photon channel, so triphoton production is dominant by orders of magnitude in that setup [2604.21838]. The same paper quotes, for ortho-dimuonium, LHC values of $\hat{\sigma}\times\mathcal{B}=1.9$--$65\,\mu{\rm b}$ and $(0.61$--$21)\times 10^5$ events/year, depending on $x_{\min}$ and detector-threshold assumptions [2604.21838]. The preferred signature is an extremely low-$p_T$ displaced $e^+e^-$ pair with invariant mass near $211$ MeV [2604.21838].

Dimuonium production has also been studied in the QGP created in relativistic heavy-ion collisions. There the dominant process is identified as
\[
q\bar q\rightarrow (\mu^+\mu^-)g,
\]
with transport evolution governed by a Boltzmann-type equation because dimuonium does not thermalize with the medium [1204.4389]. The overall yield is small, about $1.3\times 10^{-9}$ for central RHIC collisions and three times larger at LHC, but the low-$E_t$ spectrum can be fitted with effective temperatures of $195$ MeV at RHIC and $240$ MeV at LHC, implying that the transverse-energy distribution encodes QGP thermodynamics [1204.4389]. This suggests that, even where total rates are marginal, dimuonium can function as a QED probe of hot QCD matter.

## 5. Matter interactions, foil transport, and computational infrastructure

Passage through matter is central both to detector concepts and to dedicated foil-based spectroscopy proposals. The interaction is modeled as scattering in the external screened Coulomb field of atomic nuclei, typically in the non-relativistic Born approximation [2203.01131] [2405.04460]. In MuMuPy, the core observable is the atomic form factor
\[
F_{n_1 l_1 m_1}^{n_2 l_2 m_2}(\vec{q}) = \int d^3 r \; \varphi_{n_2 l_2 m_2}^*(\vec{r}) \, e^{i\vec{q} \cdot \vec{r}} \, \varphi_{n_1 l_1 m_1}(\vec{r}),
\]
from which transition and breakup cross sections are computed [2203.01131]. The screened target potential is modeled with Thomas-Fermi-Molière screening:
\[
U(q) = 4 \pi Z e \sum_{i=1}^3 \frac{\alpha_i}{q^2 + \beta_i^2}.
\]
MuMuPy implements three independent calculation strategies—Dewangan’s method, Afanasev and Tarasov’s method, and a new hypergeometric-based method—and supports arbitrary hydrogenic quantum numbers, target $Z$ from $1$ to $137$, and states up to $n\approx 13$ [2203.01131].

A dedicated transport study for dimuonium traversing Be, Al, and Pb foils formulates the state populations through the kinetic system
\[
\frac{dN_\alpha}{dl} = -n_Z \sigma(\alpha \to X) N_\alpha + \sum_{\beta \ne \alpha} n_Z \sigma(\beta \to \alpha) N_\beta.
\]
Using $15$ different approximations of the atomic potential, that work finds that the atomic-potential-model-dependent error in the yields of the low lying states is quite small: absolute cross sections can vary by as much as $8\%$, but ratios of cross sections vary by less than $1\%$, so the impact on the final low-lying yields is negligible within the applied Born approximation [2405.04460]. The same analysis reports $2P$ yields up to $12\%$ and $3P$ yields up to $3\%$ of the initial $1S$ population for a representative setup, with good convergence already by $n_{\max}=5$ for low-lying states [2405.04460]. A methodological correction is emphasized: the quantization axis should be chosen along the dimuonium beam, not along the transferred momentum [2405.04460].

Matter effects also enter fixed-target discovery proposals through dissociation. In the Drell–Yan foil mechanism, the breakup cross section is parameterized as
\[
\sigma_{\rm dis} = 1.3 \, Z^2 \times 10^{-23}\, \text{cm}^2,
\]
with survival probability
\[
W(L_{T\!M}) = \exp\left(-\frac{L_{T\!M}}{L_{T\!M}^{\rm eff}}\right).
\]
This makes target optimization a balance between production, transmission, excitation, and breakup [2503.01393].

## 6. Experimental status, theoretical caveats, and outlook

The experimental status remains unchanged: positronium was observed in 1951, whereas the search for dimuonium is still in vain [2508.14424]. What has changed is the breadth of viable production strategies and the maturity of the corresponding theory. Recent calculations argue that observation is tenable at Belle II and especially at STCF in $\gamma\gamma$ production [2508.14424], promising in $J/\psi$ decays at STCF [2505.02676], plausible in dedicated threshold machines such as DIMUS [2203.07144], and potentially accessible in fixed-target Drell–Yan-on-foil searches [2503.01393] and in UPCs at the LHC or FCC [2307.12387] [2604.21838].

Several recurring misconceptions are corrected by the current literature. First, the statement that Coulomb corrections are negligible is only partially true. In high-energy $\mu^+\mu^-$ electroproduction and paradimuonium production by a relativistic electron in an atomic field, Coulomb corrections to cross sections integrated over final-electron momentum are small, but Coulomb corrections to the differential cross sections are large: in the peak region near $p_{2\perp}\sim m_e$, the exact cross section is up to $30\%$ lower than the Born result, whereas for $m_e \ll p_{2\perp} \lesssim m_\mu$ it is about $10\%$ higher [1712.09770]. Second, EPA-based estimates are not uniformly adequate across phase space: the LHC para-dimuonium study finds a $\sim 15\%$ overestimate at low $P_T$ [2307.12387]. Third, large production rates do not automatically imply straightforward reconstruction. In UPC para-dimuonium searches, the decay photons have energies close to $105$ MeV each, and detection is described as extremely challenging because of low energies, small opening angles, and continuum/instrumental backgrounds [2503.10952].

The convergence of theory and experimental design is strongest around displaced-vertex and exclusive low-mass signatures. Ortho-dimuonium searches repeatedly target displaced $e^+e^-$ decays near $m_{e^+e^-}\simeq 211$ MeV [2503.01393] [2604.21838], while para-dimuonium programs emphasize narrow diphoton structure at the same mass scale [2307.12387] [2503.10952]. A plausible implication is that the first observation will be driven less by raw production cross section than by the channel in which kinematic cleanliness, vertexing, and trigger thresholds can be made commensurate with the picosecond-scale lifetimes and sub-GeV visible final states.

Source: https://www.emergentmind.com/topics/dimuonium