Supplement to "Improved cosmological constraints on axion-lepton interactions"

This repository contains the data and code used to generate the plots presented in "Improved cosmological constraints on axion-lepton interaction". The repository includes the axion distributions $f_a$ derived via the full Boltzmann equation (fBE), the resulting values of $\Delta N_{\mathrm{eff}}$ as functions of $f/|C|$ for all considered production processes, the constraints on $f/|C|$ for various axion masses obtained via MCMC analysis using Planck 2018 and DESI DR 2, and the densities $\omega_a$ of cold axions. The provided Jupyter notebooks enable the reproduction of all figures presented in the main text. The raw distributions can be used in CLASS to compute power spectra and perform MCMC scans with one's preferred sampler.