Axion Cooling of Neutron Stars

Armen Sedrakian

Published 2015-12-24Version 1

Cooling simulations of neutron stars and their comparison with the data from thermally emitting X-ray sources puts constraints on the properties of axions, and by extension of any light pseudo-scalar dark matter particles, whose existence has been postulated to solve the strong-CP problem of QCD. We incorporate the axion emission by pair-breaking and formation processes by $S$- and $P$-wave nucleonic condensates in a benchmark code for cooling simulations as well as provide fit formulae for the rates of these processes. Axion cooling of neutron stars has been simulated for 24 models covering the mass range 1 to 1.8 solar masses, featuring non-accreted iron and accreted light element envelopes, and a range of nucleon-axion coupling. The models are based on an equation state predicting conservative physics of superdense nuclear matter that does not allow for onset of fast cooling processes induced by phase transitions to non-nucleonic forms of matter or high proton concentration. The cooling tracks in the temperature vs age plane were confronted with the (time-averaged) measured surface temperature of the central compact object in the Cas A supernova remnant as well as surface temperatures of three nearby middle-aged thermally emitting pulsars.We find that the axion coupling is limited to $f_a/10^{10}\textrm{GeV} > 5\times 10^{-3}$ which translates into upper bound on axion mass $m_a\le 0.1$ eV for nucleon Peccei-Quinn charges $C_n \sim 0.04$.