Abstract: In warm dark matter scenarios structure formation is suppressed on small scales with respect to the cold dark matter case, reducing the number of low-mass halos and the fraction of ionized gas at high redshifts and thus, delaying reionization. This has an impact on the ionization history of the Universe and measurements of the optical depth to reionization, of the evolution of the global fraction of ionized gas and of the thermal history of the intergalactic medium, can be used to set constraints on the mass of the dark matter particle. However, the suppression of the fraction of ionized medium in these scenarios can be partly compensated by varying other parameters, as the ionization efficiency or the minimum mass for which halos can host star-forming galaxies. Here we use different data sets regarding the ionization and thermal histories of the Universe and, taking into account the degeneracies from several astrophysical parameters, we obtain a lower bound on the mass of thermal warm dark matter candidates of m(X) > 1.3 keV, or m(s) > 5.5 keV for the case of sterile neutrinos nonresonantly produced in the early Universe, both at 90% confidence level.

Abstract: We present a novel approach to derive constraints on neutrino masses, as well as on other cosmological parameters, from cosmological data, while taking into account our ignorance of the neutrino mass ordering. We derive constraints from a combination of current as well as future cosmological datasets on the total neutrino mass M-nu and on the mass fractions f(nu),i = m(i)/M-nu (where the index i = 1, 2, 3 indicates the three mass eigenstates) carried by each of the mass eigenstates m(i), after marginalizing over the (unknown) neutrino mass ordering, either normal ordering (NH) or inverted ordering (IH). The bounds on all the cosmological parameters, including those on the total neutrino mass, take therefore into account the uncertainty related to our ignorance of the mass hierarchy that is actually realized in nature. This novel approach is carried out in the framework of Bayesian analysis of a typical hierarchical problem, where the distribution of the parameters of the model depends on further parameters, the hyperparameters. In this context, the choice of the neutrino mass ordering is modeled via the discrete hyperparameter h(type), which we introduce in the usual Markov chain analysis. The preference from cosmological data for either the NH or the IH scenarios is then simply encoded in the posterior distribution of the hyper-parameter itself. Current cosmic microwave background (CMB) measurements assign equal odds to the two hierarchies, and are thus unable to distinguish between them. However, after the addition of baryon acoustic oscillation (BAO) measurements, a weak preference for the normal hierarchical scenario appears, with odds of 4 : 3 from Planck temperature and large-scale polarization in combination with BAO (3 : 2 if small-scale polarization is also included). Concerning next-generation cosmological experiments, forecasts suggest that the combination of upcoming CMB (COrE) and BAO surveys (DESI) may determine the neutrino mass hierarchy at a high statistical significance if the mass is very close to the minimal value allowed by oscillation experiments, as for NH and a fiducial value of M-nu = 0.06 eV there is a 9 : 1 preference of normal versus inverted hierarchy. On the contrary, if the sum of the masses is of the order of 0.1 eV or larger, even future cosmological observations will be inconclusive. The innovative statistical strategy exploited here represents a very simple, efficient and robust tool to study the sensitivity of present and future cosmological data to the neutrino mass hierarchy, and a sound competitor to the standard Bayesian model comparison. The unbiased limit on M-nu we obtain is crucial for ongoing and planned neutrinoless double beta decay searches.

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.