Constraints on axion-like particles using lattice QCD calculations of the rate for $J/\psi \to \gamma a$

Abstract

A key search mode for axion-like particles (ALPs) that couple to charm quarks is $J/\psi \to \gamma a$. Here we calculate the form factor that allows the rate of this process to be determined using lattice QCD for the first time. Our calculations use the relativistic Highly Improved Staggered Quark (HISQ) action for the valence charm quarks on gluon field configurations generated by the MILC collaboration that include $u$, $d$, $s$ and $c$ HISQ quarks in the sea at four values of the lattice spacing and both unphysical and physical sea quark masses. We determine the form factor as a function of ALP mass with an uncertainty of less than 2\% across our full range of ALP masses from zero up to 95\% of the $J/\psi$ mass. This represents a substantial improvement in accuracy of the theoretical picture of this decay compared to the previously used tree-level and $\mathcal{O}(\alpha_s)$ perturbation theory. We use our form factor to determine constraints on ALP masses and couplings to charm quarks and photons in several different scenarios using recent experimental data from BESIII. Our calculation paves the way for further lattice QCD input on new physics constraints from radiative decays.