The cumulative emission of axionlike particles (ALPs) from all past core-collapse supernovae (SNe) would lead to a diffuse flux with energies $\mathcal{O}(50)\mathrm{MeV}$. We use this to constrain ALPs featuring couplings to photons and to nucleons. ALPs coupled only to photons are produced in the SN core via the Primakoff process and then converted into gamma rays in the Galactic magnetic field. We set a bound on $g_{a\gamma }\lesssim 5\times {10}^{-10}{\mathrm{GeV}}^{-1}$ for $m_a\lesssim {10}^{-11}\mathrm{eV}$, using recent measurements of the diffuse gamma-ray flux observed by the Fermi-LAT telescope. However, if ALPs couple also with nucleons, their production rate in SN can be considerably enhanced due to the ALPs nucleon-nucleon bremsstrahlung process. Assuming the largest ALP-nucleon coupling phenomenologically allowed, bounds on the diffuse gamma-ray flux lead to a much stronger $g_{a\gamma }\lesssim 6\times {10}^{-13}{\mathrm{GeV}}^{-1}$ for the same mass range. If ALPs are heavier than $\sim \mathrm{keV}$, the decay into photons becomes significant, leading again to a diffuse gamma-ray flux. In the case of only photon coupling, we find, e.g., $g_{a\gamma }\lesssim 5\times {10}^{-11}{\mathrm{GeV}}^{-1}$ for $m_a\sim 5\mathrm{keV}$. Allowing for a (maximal) coupling to nucleons, the limit improves to the level of $g_{a\gamma }\lesssim {10}^{-19}{\mathrm{GeV}}^{-1}$ for $m_a\sim 20\mathrm{MeV}$, which represents the strongest constraint to date.

• Received 1 September 2020
• Accepted 30 October 2020

DOI:https://doi.org/10.1103/PhysRevD.102.123005

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP³.

Published by the American Physical Society