Upper Limit on the QCD Axion Mass from Isolated Neutron Star Cooling

The quantum chromodynamics (QCD) axion may modify the cooling rates of neutron stars (NSs). The axions are produced within the NS cores from nucleon bremsstrahlung and, when the nucleons are in superfluid states, Cooper pair breaking and formation processes. We show that four of the nearby isolated magnificent seven NSs along with PSR J0659 are prime candidates for axion cooling studies because they are coeval, with ages of a few hundred thousand years known from kinematic considerations, and they have well-measured surface luminosities. We compare these data to dedicated NS cooling simulations incorporating axions, profiling over uncertainties related to the equation of state, NS masses, surface compositions, and superfluidity. Our calculations of the axion and neutrino emissivities include high-density suppression factors that also affect SN 1987A and previous NS cooling limits on axions. We find no evidence for axions in the isolated NS data, and within the context of the Kim-Shifman-Vainshtein-Zakharov QCD axion model, we constrain m_{a}≲16  meV at 95% confidence level. An improved understanding of NS cooling and nucleon superfluidity could further improve these limits or lead to the discovery of the axion at weaker couplings.