Abstract
Axions and other light particles appear ubiquitously in physics beyond the Standard Model, with a variety of possible couplings to ordinary matter. Cosmology offers a unique probe of these particles as they can thermalize in the hot environment of the early universe for any such coupling. For sub-MeV particles, their entropy must leave a measurable cosmological signal, usually via the effective number of relativistic particles, Neff. In this paper, we will revisit the cosmological constraints on the couplings of axions and other pseudo-Nambu-Goldstone bosons to Standard Model fermions from thermalization below the electroweak scale, where these couplings are marginal and give contributions to the radiation density of ΔNeff > 0.027. We update the calculation of the production rates to eliminate unnecessary approximations and find that the cosmological bounds on these interactions are complementary to astrophysical constraints, e.g. from supernova SN 1987A. We additionally provide quantitative explanations for these bounds and their relationship.