We calculate the axion emission rate from reactions involving thermal pions in matter encountered in supernovae and neutron star mergers, identify unique spectral features, and explore their implications for astrophysics and particle physics. We find that it is about 2–5 times larger than nucleon-nucleon bremsstrahlung, which in past studies was considered to be the dominant process. The axion spectrum is also found be much harder. Together, the larger rates and higher axion energies imply a stronger bound on the mass of the QCD axion and better prospects for direct detection in a large underground neutrino detector from a nearby galactic supernova.

We perform general relativistic one-dimensional supernova (SN) simulations to identify observable signatures of enhanced axion emission from the pion-induced reaction ${\pi }^{-}+p\to n+a$ inside a newly born protoneutron star (PNS). We focus on the early evolution after the onset of the supernova explosion to predict the temporal and spectral features of the neutrino and axion emission during the first 10 s. Pions are included as explicit new degrees of freedom in hot and dense matter. Their thermal population and their role in axion production are both determined consistently to include effects due to their interactions with nucleons. For a wide range of ambient conditions encountered inside a PNS, we find that the pion-induced axion production dominates over nucleon-nucleon bremsstrahlung processes. By consistently including the role of pions on the dense matter equation of state and on the energy loss, our simulations predict robust discernible features of neutrino and axion emission from a galactic supernova that can be observed in terrestrial detectors. For axion couplings that are compatible with current bounds, we find a significant suppression with time of the neutrino luminosity during the first 10 s. This suggests that current bounds derived from the neutrino signal from SN 1987A can be improved and that future galactic supernovae may provide significantly more stringent constraints.

Axionlike particles (ALPs) coupled to nucleons might be copiously emitted from a supernova (SN) core. We extend existing bounds on free-streaming ALPs to the case in which these are so strongly interacting with the nuclear matter to be trapped in the SN core. For strongly interacting ALPs, we also extend the bound from the absence of an ALP-induced signal in Kamiokande-II neutrino detector at the time of SN 1987A. We find that combining the different arguments, SNe exclude values of ALP-nucleon coupling $g_{aN}\gtrsim {10}^{-9}$ for ALP masses $m_a\lesssim 1\mathrm{MeV}$. Remarkably, in the case of canonical QCD axion models, the SN bounds exclude all values of $m_a\gtrsim {10}^{-2}\mathrm{eV}$. This result prevents the possibility for current and future cosmological surveys to detect any signatures due to hot dark matter QCD axion mass.