Toponium: Implementation of a toponium model in FeynRules (2504.12634v1)

Abstract: Toponium -- a bound state of the top-antitop pair ($t\bar{t}$) -- emerges as the smallest and simplest hadronic system in QCD, with an ultrashort lifetime ($\tau_t \sim 2.5\times 10^{-25}$~s) and a femtometer-scale Bohr radius ($r_{\text{Bohr}} \sim 7\times 10^{-18}$~m). We present a computational framework extending the Standard Model (SM) with two S-wave toponium states: a spin-singlet $\eta_t$ ($J^{PC}=0^{-+}$) and a spin-triplet $J_t$ ($J^{PC}=1^{--}$). Using nonrelativistic QCD (NRQCD) and a Coulomb potential, we derived couplings to SM particles (gluons, electroweak bosons, Higgs boson, and fermion pairs) and implemented the Lagrangian in FeynRules, generating FeynArts, MadGraph, and WHIZARD models for collider simulations. Key results include dominant decay channels ($\eta_t \to gg/ZH$, $J_t \to W^+W^-/b\bar{b}$) and leading order (LO) cross sections for $pp \to \eta_t(nS) \to {\rm non-}t\bar{t}$ ({66 fb} at 13 TeV). The model avoids double-counting artifacts by excluding direct $t\bar{t}$ couplings, thereby ensuring consistency with perturbative QCD. This work establishes a complete pipeline for precision toponium studies, bridging NRQCD, collider phenomenology, and tests of SM validity at future lepton colliders (e.g., CEPC, FCC-ee, muon colliders) and the LHC. This provides the first publicly available UFO model for toponium, enabling direct integration with MadGraph and WHIZARD for simulations.