From supernovae to neutron stars: a systematic approach to axion production at finite density

We present a systematic study of QCD axion production in environments with finite baryon density and temperature, implying significant changes to axion phenomenology. Within heavy baryon chiral perturbation theory, we derive the effective Lagrangian describing axion interactions with nucleons and mesons up to next-to-leading-order in the chiral expansion. We focus on corrections to the axion-nucleon couplings from higher orders and finite density. These couplings are modified by up to an order of magnitude near nuclear saturation density, significantly impacting axion production in supernovae and neutron stars. Density-dependent corrections enhance the axion luminosity in supernovae by an order of magnitude, strengthening current best bounds by a factor of three. We stress the importance of including all axion production channels up to a given chiral order for a consistent luminosity calculation and classify the missing contributions up to the third chiral order. The modified axion-nucleon couplings also affect neutron star cooling rates via axion emission. A re-evaluation of existing neutron star cooling bounds, constrained to regions where perturbative control is reliable, weakens these bounds by a factor of four. Lastly, our results have implications for terrestrial axion searches that rely on precise knowledge of axion-nucleon couplings.