Abstract: We propose here a novel method which singles out the a priori unavoidable dependence on the underlying cosmological model when extracting parameter constraints, providing robust limits which only depend on the considered dataset. Interestingly, when dealing with several possible cosmologies and interpreting the Bayesian preference in terms of the Gaussian statistical evidence, the preferred model is much less favored than when only two cases are compared. As a working example, we apply our approach to the cosmological neutrino mass bounds, which play a fundamental role not only in establishing the contribution of relic neutrinos to the dark matter of the Universe but also in the planning of future experimental searches of the neutrino character and of the neutrino mass ordering.

Abstract: We revisit the joint constraints in the mixed hot dark matter scenario in which both thermally produced QCD axions and relic neutrinos are present. Upon recomputing the cosmological axion abundance via recent advances in the literature, we improve the state-of-the-art analyses and provide updated bounds on axion and neutrino masses. By avoiding approximate methods, such as the instantaneous decoupling approximation, and limitations due to the limited validity of the perturbative approach in QCD that forced to artificially divide the constraints from the axion-pion and the axion-gluon production channels, we find robust and self-consistent limits. We investigate the two most popular axion frameworks: KSVZ and DFSZ. From Big Bang Nucleosynthesis (BBN) light element abundances data we find for the KSVZ axion Delta N-eff < 0.31 and an axion mass bound m(a) < 0.53 eV (i.e., a bound on the axion decay constant f(a) > 1.07 x 10(7) GeV) both at 95% CL. These BBN bounds are improved to Delta N-eff < 0.14 and m(a) < 0.16 eV (f(a) > 3.56 x 10(7) GeV) if a prior on the baryon energy density from Cosmic Microwave Background (CMB) data is assumed. When instead considering cosmological observations from the CMB temperature, polarization and lensing from the Planck satellite combined with large scale structure data we find Delta N-eff < 0.23, m(a) < 0.28 eV (f(a) > 2.02 x 10(7) GeV) and Sigma m(nu) < 0.16 eV at 95% CL. This corresponds approximately to a factor of 5 improvement in the axion mass bound with respect to the existing limits. Very similar results are obtained for the DFSZ axion. We also forecast upcoming observations from future CMB and galaxy surveys, showing that they could reach percent level errors for m(a) similar to 1 eV.

Abstract: We study the impact of assumptions about neutrino properties on the estimation of inflationary parameters from cosmological data, with a specific focus on the allowed contours in the n(s)/r plane, where n(s) is the scalar spectral index and r is the tensor-to-scalar ratio. We study the following neutrino properties: (i) the total neutrino mass M-i = Sigma(i)m(i) (where the index i = 1, 2, 3 runs over the three neutrino mass eigenstates); (ii) the number of relativistic degrees of freedom N-eff at the time of recombination; and (iii) the neutrino hierarchy. Whereas previous literature assumed three degenerate neutrino masses or two massless neutrino species (approximations that clearly do not match neutrino oscillation data), we study the cases of normal and inverted hierarchy. Our basic result is that these three neutrino properties induce < 1 sigma shift of the probability contours in the n(s)/r plane with both current or upcoming data. We find that the choice of neutrino hierarchy (normal, inverted, or degenerate) has a negligible impact. However, the minimal cutoff on the total neutrino mass M-v,M-min = 0 that accompanies previous works using the degenerate hierarchy does introduce biases in the n(s)/r plane and should be replaced by M-v,M-min = 0.059 eV as required by oscillation data. Using current cosmic microwave background (CMB) data from Planck and Bicep/Keck, marginalizing over the total neutrino mass M-v and over r can lead to a shift in the mean value of ns of similar to 0.3 sigma toward lower values. However, once baryon acoustic oscillation measurements are included, the standard contours in the n(s)/r plane are basically reproduced. Larger shifts of the contours in the n(s)/r plane (up to 0.8 sigma) arise for nonstandard values of N-eff. We also provide forecasts for the future CMB experiments Cosmic Origins Explorer (COrE, satellite) and Stage-IV (ground-based) and show that the incomplete knowledge of neutrino properties, taken into account by a marginalization over M-v, could induce a shift of similar to 0.4 sigma toward lower values in the determination of ns (or a similar to 0.8 sigma shift if one marginalizes over N-eff). Comparison to specific inflationary models is shown. Imperfect knowledge of neutrino properties must be taken into account properly, given the desired precision in determining whether or not inflationary models match the future data.