Precise predictions for the QCD axion contribution to dark radiation with full phase-space evolution

TL;DR

This paper provides precise predictions for the QCD axion's contribution to dark radiation by solving full phase-space Boltzmann equations, improving upon previous approximations and impacting experimental bounds.

Contribution

It introduces a comprehensive method to compute axion-induced dark radiation using full phase-space evolution, accounting for quantum statistics, applicable to various production mechanisms.

Findings

Predictions for ΔN_eff can differ significantly from previous estimates.

Current bounds on axion couplings are relaxed based on precise calculations.

The approach improves accuracy for future CMB measurements.

Abstract

We compute the QCD axion contribution to the energy density of dark radiation, parameterized by $\Delta N_{\rm eff}$, by solving Boltzmann equations for the momentum distribution functions including the effects of quantum statistics for all particles involved in the axion production processes. This approach gives precise prediction for $\Delta N_{\rm eff}$ independently on whether axions are produced via freeze-out or freeze-in. We focus on axions produced via flavor-conserving and flavor-violating interactions with leptons. Our precise predictions for $\Delta N_{\rm eff}$ can differ from those assuming thermal shape for the momentum distribution functions, as commonly done in the literature, by more than the experimental precision of future Cosmic Microwave Background (CMB) observations. Current lower limits on the axion couplings from Planck constraints on $\Delta N_{\rm eff}$ are also affected by our precise computation which, in particular, results in a strongly relaxed bound on flavor-violating axion couplings to tau lepton and muon or electron.