Toponium: The Top Quark Bound State

• Toponium is the quantum bound state of a top quark and antiquark, characterized by a minute Bohr radius (~8×10⁻¹⁸ m) and an ultrashort lifetime (~2.5×10⁻²⁵ s).
• Its binding energy (~3 GeV) and dominant Coulombic interaction in the ultraviolet regime offer a pristine lab for testing perturbative QCD and asymptotic freedom.
• Experimental signatures include resonance enhancements near 343–345 GeV, facilitating high-precision top quark mass measurements and detailed studies of QCD dynamics.

Toponium is the quantum bound state of a top quark ($t$) and its antiquark ($\bar t$), bound by the strong interaction. It constitutes the smallest spatially extended hadronic system in nature, with defining scales characterized by a Bohr radius of $r_B\simeq 8\times 10^{-18}$ m and an ultrashort lifetime $\tau \simeq 2.5\times 10^{-25}$ s. The study of toponium provides a theoretically clean and perturbative laboratory for quantum chromodynamics (QCD) at the highest available energies and shortest distances, offering distinctive insights compared to other quarkonium systems such as charmonium and bottomonium. Unlike all other known hadrons, the toponium system is governed by ultraviolet freedom, probing distances well below the QCD confinement scale and thus directly sensitive to the running of the strong coupling in the asymptotically free regime (Fu et al., 2024).

1. Quantum Mechanical Structure and Fundamental Parameters

The theoretical description of toponium is grounded in the nonrelativistic quantum-mechanical two-body Hamiltonian,

$H = -\frac{\nabla^2}{2\mu} + V(r)$

with reduced mass $\mu = m_t/2$, where $m_t \simeq 172.5$ GeV is the top-quark mass (Fu et al., 2024). The QCD interquark potential is taken to be a sum of a Coulomb term and a linear confinement term,

$V(r) = -\lambda/r + \sigma\,r$

with $\lambda = 0.285 \pm 0.011$ and $\sigma = 0.206$ GeV$\bar t$0. Due to the extremely small $\bar t$1 compared to $\bar t$2, the system is dominantly controlled by the Coulombic interaction, with the linear term yielding only subleading corrections (e.g., $\bar t$3 MeV) to the binding energy.

The leading-order Bohr radius, quantifying the spatial extent of the system, is

$\bar t$4

and the binding energy in the Coulomb approximation is

$\bar t$5

with a corresponding wavefunction at the origin,

$\bar t$6

(Fu et al., 2024).

The three-dimensional ground-state energy, including NNNLO perturbative corrections as well as Higgs and photon exchange, reaches $\bar t$7 GeV. The weak decay of the constituent quarks dominates the total width,

$\bar t$8

with $\bar t$9 GeV, giving $r_B\simeq 8\times 10^{-18}$0 GeV and a resulting lifetime $r_B\simeq 8\times 10^{-18}$1 s, implying that the toponium bound state decays before hadronizing (Fu et al., 2024). In this regime, strong decay channels are suppressed and all hadronic final states arise from the $r_B\simeq 8\times 10^{-18}$2 and $r_B\simeq 8\times 10^{-18}$3 decays themselves.

2. Distinctive Regime: Ultraviolet Freedom and Comparison to Other Quarkonia

Unlike charmonium or bottomonium, which probe QCD in the confining, infrared-dominated regime with $r_B\simeq 8\times 10^{-18}$4, toponium is characterized by $r_B\simeq 8\times 10^{-18}$5, so that $r_B\simeq 8\times 10^{-18}$6 and the entire system remains in the ultraviolet, perturbative regime of QCD—an explicit realization of asymptotic freedom in a bound state (Fu et al., 2024, Thompson, 14 Jul 2025). In particular, toponium functions as a "hydrogen atom" of QCD at the smallest distances, with all hadronic properties calculable in perturbation theory to high accuracy.

Leading contributions to the binding arise from the Coulombic potential with $r_B\simeq 8\times 10^{-18}$7–$r_B\simeq 8\times 10^{-18}$8, tightly matched to the QCD running coupling at the relevant scale, e.g., $r_B\simeq 8\times 10^{-18}$9 for $\tau \simeq 2.5\times 10^{-25}$0 (Fu et al., 2024, Fu et al., 17 Apr 2025).

This regime is unique among hadrons.

• The lifetime $\tau \simeq 2.5\times 10^{-25}$1 is orders of magnitude shorter than any typical hadronic time: $\tau \simeq 2.5\times 10^{-25}$2 s for toponium, compared to $\tau \simeq 2.5\times 10^{-25}$3 s for bottomonium and $\tau \simeq 2.5\times 10^{-25}$4 s for charmonium (Thompson, 14 Jul 2025).
• The binding energy ($\tau \simeq 2.5\times 10^{-25}$5–$\tau \simeq 2.5\times 10^{-25}$6 GeV) is small relative to twice the constituent mass but still larger than the constituent width, enabling a well-defined resonance feature.
• The near-degeneracy of S-wave levels in the absence of confining effects leads to a "Bohr tower" spectrum $\tau \simeq 2.5\times 10^{-25}$7 (Fu et al., 2024, Lopez et al., 4 Aug 2025).

3. Theoretical Approaches and Computational Methodologies

The structure and spectrum of toponium have been computed using a range of techniques:

• Nonrelativistic Schrödinger and Bethe–Salpeter Equations: These provide analytic expressions for $\tau \simeq 2.5\times 10^{-25}$8, $\tau \simeq 2.5\times 10^{-25}$9, and the spatial wavefunctions under the Cornell or pure Coulombic potential (Fu et al., 2024, Llanes-Estrada, 2024, Lopez et al., 4 Aug 2025).
• Quantum Bootstrap Method: Translating the Schrödinger equation into recursion relations for radial moments $H = -\frac{\nabla^2}{2\mu} + V(r)$0 and requiring the positivity of Hankel matrices in these moments yields bounds on the allowed energy eigenvalues, leading to a $H = -\frac{\nabla^2}{2\mu} + V(r)$1 toponium mass $H = -\frac{\nabla^2}{2\mu} + V(r)$2 GeV. Agreement of this prediction with observed $H = -\frac{\nabla^2}{2\mu} + V(r)$3 threshold enhancements at the LHC points to the reliability of this formalism (Lopez et al., 4 Aug 2025).
• QCD Sum Rules: Using operator product expansions up to dimension-eight condensates in two-point current correlators, sum-rule extractions yield spin-singlet and spin-triplet masses $H = -\frac{\nabla^2}{2\mu} + V(r)$4 GeV and $H = -\frac{\nabla^2}{2\mu} + V(r)$5 GeV, respectively. These are in close agreement with experimentally reported excesses and feature negative binding energies $H = -\frac{\nabla^2}{2\mu} + V(r)$6 GeV, signalling genuine bound states below threshold (Najjar et al., 13 Nov 2025).
• Resummed NRQCD Green's Functions: The effect of finite top width and higher-order QCD corrections is incorporated by evaluating the imaginary part of the Green's function of the nonrelativistic QCD Hamiltonian. This approach predicts a broad threshold enhancement rather than a narrow peak, in accordance with the top's large decay width fully smearing out the Bohr tower of levels (Fuks, 6 May 2025, Garzelli et al., 2024).

A summary of representative numerical predictions from different approaches is given below.

Method	$H = -\frac{\nabla^2}{2\mu} + V(r)$7 [GeV]	$H = -\frac{\nabla^2}{2\mu} + V(r)$8 [GeV]	Notes
Potential/NRQCD (Fu et al., 2024, Fu et al., 17 Apr 2025, Bai et al., 17 Jun 2025)	341–343	–3.1 to –2.7	$H = -\frac{\nabla^2}{2\mu} + V(r)$9–$\mu = m_t/2$0
Bootstrap (Lopez et al., 4 Aug 2025)	344.3	–0.69	Cornell, $\mu = m_t/2$1
Sum rules (Najjar et al., 13 Nov 2025)	343.5	–1.6	Up to dim-8

4. Experimental Signatures and Collider Phenomenology

The key observable for toponium is a threshold enhancement in the $\mu = m_t/2$2 invariant mass distribution at $\mu = m_t/2$3–$\mu = m_t/2$4 GeV, as reported by both ATLAS and CMS with measured excesses attaining local significances above $\mu = m_t/2$5 in the dileptonic decay mode (Desai et al., 26 Jan 2026, Fu et al., 2024).

A particularly sensitive observable at lepton colliders is the cross section ratio

$\mu = m_t/2$6

where toponium resonance formation produces a pronounced dip–bump feature via interference with the $\mu = m_t/2$7 and $\mu = m_t/2$8 continuum amplitudes, observable in a threshold energy scan at $\mu = m_t/2$9–$m_t \simeq 172.5$0 GeV (Fu et al., 2024). The measurement of $m_t \simeq 172.5$1 around threshold enables extraction of the 1S top mass with a total uncertainty $m_t \simeq 172.5$2 MeV, a factor of ten better than the current on-shell mass uncertainty (Fu et al., 2024).

At the LHC, toponium appears as a moderate ($m_t \simeq 172.5$3–$m_t \simeq 172.5$4 pb) enhancement in the $m_t \simeq 172.5$5 cross section just below twice the top mass, with no measurable interference between the signal and continuum QCD background at leading order (Djouadi et al., 2024). Exclusive production, primarily in the $m_t \simeq 172.5$6 channel, is extremely rare (few attobarns at the LHC), but potentially observable with forward proton tagging at future high-luminosity facilities (Francener et al., 5 Feb 2025).

Dileptonic final states exhibit distinctive kinematic correlations: small dilepton invariant mass ($m_t \simeq 172.5$7 GeV) and azimuthal separation ($m_t \simeq 172.5$8) sharply distinguish the toponium contribution from the continuum (Fuks et al., 2021).

5. Experimental Techniques for Reconstruction and Discrimination

Experimental isolation of the toponium signal in $m_t \simeq 172.5$9 events is challenged by missing energy from neutrinos and combinatorial ambiguity in jet–lepton assignment. Recent results employ Recursive Jigsaw Reconstruction (RJR) to resolve event kinematics, partitioning each event into a hierarchy of rest frames and applying constraints (e.g., equal reconstructed top masses) to fully reconstruct the event topology (Desai et al., 26 Jan 2026, Desai et al., 27 Jan 2026).

Key discriminating observables include:

• The azimuthal angle difference $V(r) = -\lambda/r + \sigma\,r$0 between reconstructed tops.
• The "Chernyak–Zhitnitsky–like" lepton angular variable $V(r) = -\lambda/r + \sigma\,r$1, encoding boosted lepton correlation in rest frames.

Inclusion of these variables in analysis grids yields a local significance up to $V(r) = -\lambda/r + \sigma\,r$2 for the toponium signal, surpassing traditional methods by $V(r) = -\lambda/r + \sigma\,r$3 (Desai et al., 26 Jan 2026, Desai et al., 27 Jan 2026).

Complementary studies in the single-lepton channel exploit differences in lepton $V(r) = -\lambda/r + \sigma\,r$4, missing energy, reconstructed $V(r) = -\lambda/r + \sigma\,r$5, and top recoil in the rest frame to distinguish the signal. With current data (Run 2, $V(r) = -\lambda/r + \sigma\,r$6 fb$V(r) = -\lambda/r + \sigma\,r$7), significance exceeds $V(r) = -\lambda/r + \sigma\,r$8 in the optimal kinematic region (Fuks et al., 3 Sep 2025).

6. Implications for QCD, Mass Measurements, and Fundamental Quantum Tests

Beyond perturbative QCD validation, toponium enables precision determinations of the 1S top mass, reducing the dominant uncertainties to the tens-of-MeV level—orders of magnitude better than alternative techniques (Fu et al., 2024). The bound state also provides a window onto the top Yukawa coupling and the short-distance QCD potential via detailed studies of production rates, resonance shape, and level splittings (Bai et al., 17 Jun 2025).

Toponium formation has further ramifications for quantum foundations. The competition between the ultrashort formation time ($V(r) = -\lambda/r + \sigma\,r$9 s) and the top quark's own lifetime renders the system a real-time testbed for the interplay between quantum superposition and relativistic causality. At future colliders, $\lambda = 0.285 \pm 0.011$0 lineshape measurements can differentiate "instantaneous" superposition (A=0) versus "causal" delayed formation ($\lambda = 0.285 \pm 0.011$1) at $\lambda = 0.285 \pm 0.011$2 significance, offering collider-based probes of wavefunction collapse at the yoctosecond scale (Xiong et al., 8 Jul 2025).

The toponium enhancement must also be discriminated from possible new physics, such as a pseudoscalar Higgs boson near $\lambda = 0.285 \pm 0.011$3. Distinctive features such as the absence of signal–background interference and differing energy-dependence in associated production furnish clear diagnostic criteria (Djouadi et al., 2024).

7. Broader Theoretical and Model-Building Impact

Toponium serves as a sensitive laboratory for probing possible deformations of standard QCD:

• Nonlocal ultraviolet completions: Modified Bethe–Salpeter equations with entire-function regulators produce calculable shifts in resonance mass and width, potentially testable in high-precision threshold scans (Thompson, 14 Jul 2025).
• Holomorphic renormalization group flows: Introduction of a holomorphic IR fixed point in the $\lambda = 0.285 \pm 0.011$4-function yields an $\lambda = 0.285 \pm 0.011$5 upward shift in the threshold cross section, which is unobservable in lighter quarkonia but potentially measurable in toponium (Thompson, 14 Jul 2025).
• Constraints on short-range or contact interactions: The sensitivity of the threshold enhancement to possible new-physics Wilson coefficients places novel bounds inaccessible in other systems (Llanes-Estrada, 2024).

The availability of a complete toponium model in FeynRules, including public UFO modules for MadGraph and WHIZARD, enables direct simulation of both production and decay processes at colliders, facilitating future experimental and phenomenological investigations (Fu et al., 17 Apr 2025).

---

References:

(Fu et al., 2024, Bai et al., 17 Jun 2025, Lopez et al., 4 Aug 2025, Najjar et al., 13 Nov 2025, Llanes-Estrada, 2024, Fuks, 6 May 2025, Djouadi et al., 2024, Thompson, 14 Jul 2025, Francener et al., 5 Feb 2025, Desai et al., 26 Jan 2026, Desai et al., 27 Jan 2026, Fuks et al., 3 Sep 2025, Fu et al., 17 Apr 2025, Garzelli et al., 2024, Xiong et al., 8 Jul 2025, Fuks et al., 2021).