Universal bound on QCD axions from supernovae

We identify a new production channel for quantum chromodynamics (QCD) axions in supernova environments that contributes to axion emissivity for all models solving the strong <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mi>C</a:mi><a:mi>P</a:mi></a:math> problem. This channel arises at tree-level from a shift-symmetry-breaking operator constructed at next-to-leading order in chiral perturbation theory. In scenarios where model-dependent derivative couplings to nucleons are absent, this sets the strongest model-independent constraint on the axion mass, improving on existing bounds by two orders of magnitude.