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

TL;DR

This paper presents the first lattice QCD calculation of the form factor for the decay $J/ar{psi} o ar{gamma} a$, enabling precise constraints on axion-like particles coupling to charm quarks from experimental data.

Contribution

It introduces a novel lattice QCD computation of the decay form factor for $J/ar{psi} o ar{gamma} a$, improving theoretical accuracy for ALP constraints.

Findings

Form factor determined with less than 2% uncertainty across ALP mass range.

Provides new constraints on ALP masses and couplings from BESIII data.

Enhances the theoretical understanding of $J/ar{psi}$ radiative decays involving ALPs.

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.