Axion cold dark matter: status after Planck and BICEP2
  We investigate the axion dark matter scenario (ADM), in which axions account for all of the dark matter in the Universe, in light of the most recent cosmological data. In particular, we use the Planck temperature data, complemented by WMAP E-polarization measurements, as well as the recent BICEP2 observations of B-modes. Baryon Acoustic Oscillation data, including those from the Baryon Oscillation Spectroscopic Survey, are also considered in the numerical analyses. We find that, in the minimal ADM scenario, the full dataset implies that the axion mass m_a = 82.2 pm 1.1 {\mu}eV (corresponding to the Peccei-Quinn symmetry being broken at a scale f_a = (7.54 pm 0.10)*10^10 GeV), or m_a = 76.6 pm 2.6 {\mu}eV (f_a = (8.08 pm 0.27)*10^10 GeV) when we allow for a non- standard effective number of relativistic species Neff . We also find a 2{\sigma} preference for Neff > 3.046. The limit on the sum of neutrino masses is Sum m_{\nu} < 0.25eV at 95% CL for Neff = 3.046, or Sum m_{\nu} < 0.47 eV when Neff is a free parameter. Considering extended scenarios where either the dark energy equation-of-state parameter w, the tensor spectral index nt or the running of the scalar index dns/dlnk are allowed to vary does not change significantly the axion mass-energy density constraints. However, in the case of the full dataset exploited here, there is a preference for a non-zero tensor index or scalar running, driven by the different tensor amplitudes implied by the Planck and BICEP2 observations. 
E. Di Valentino, E. Giusarma, M. Lattanzi, A. Melchiorri, O. Mena
astro-ph.CO hep-ph