Systematically Testing Singlet Models for $(g-2)_\mu$

We comprehensively study all viable new-physics scenarios that resolve the muon $(g-2)_\mu$ anomaly with only Standard Model singlet particles coupled to muons via renormalizable interactions. Since such models are only viable in the MeV -- TeV mass range and require sizable muon couplings, they predict abundant accelerator production through the same interaction that resolves the anomaly. We find that a combination of fixed-target (NA64$\mu$, $M^3$), $B$-factory (BABAR, Belle II), and collider (LHC, muon collider) searches can cover nearly all viable singlets scenarios, independently of their decay modes. In particular, future muon collider searches offer the only certain test of singlets above the GeV scale, covering all higher masses up to the TeV-scale unitarity limit for these models. Intriguingly, we find that $\mathcal{O}(100~\mathrm{GeV})$ muon colliders may yield better coverage for GeV-scale singlets compared to TeV-scale concepts, which has important implications for the starting center-of-mass energy of a staged muon collider program.