Abstract: The extraction of the neutrino mass ordering is one of the major challenges in particle physics and cosmology, not only for its implications for a fundamental theory of mass generation in nature, but also for its decisive role in the scale of future neutrinoless double beta decay experimental searches. It has been recently claimed that current oscillation, beta decay and cosmological limits on the different observables describing the neutrino mass parameter space provide robust decisive Bayesian evidence in favor of the normal ordering of the neutrino mass spectrum [1]. We further investigate these strong claims using a rich and wide phenomenology, with different sampling techniques of the neutrino parameter space. Contrary to the findings of Jimenez et al. [1], no decisive evidence for the normal mass ordering is found. Neutrino mass ordering analyses must rely on priors and parameterizations that are ordering-agnostic: robust results should be regarded as those in which the preference for the normal neutrino mass ordering is driven exclusively by the data, while we find a difference of up to a factor of 33 in the Bayes factors among the different priors and parameterizations exploited here. An ordering-agnostic prior would be represented by the case of parameterizations sampling over the two mass splittings and a mass scale, or those sampling over the individual neutrino masses via normal prior distributions only. In this regard, we show that the current significance in favor of the normal mass ordering should be taken as 2.7 sigma (i.e. moderate evidence), mostly driven by neutrino oscillation data. Let us stress that, while current data favor NO only mildly, we do not exclude the possibility that this may change in the future. Eventually, upcoming oscillation and cosmological data may (or may not) lead to a more significant exclusion of IO.

Abstract: Axions currently provide the most compelling solution to the strong CP problem. These particles may be copiously produced in the early universe, including via thermal processes. Therefore, relic axions constitute a hot dark matter component and their masses are strongly degenerate with those of the three active neutrinos, as they leave identical signatures in the different cosmological observables. In addition, thermal axions, while still relativistic states, also contribute to the relativistic degrees of freedom, parameterized via N-eff. We present the cosmological bounds on the relic axion and neutrino masses, exploiting the full Planck mission data, which include polarization measurements. In the mixed hot dark matter scenario explored here, we find the tightest and more robust constraint to date on the sum of the three active neutrino masses, Sigma m nu < 0.136eV at 95% CL, as it is obtained in the very well-known linear perturbation regime. The Planck Sunyaev-Zeldovich cluster number count data further tightens this bound, providing a 95% CL upper limit of Sigma m nu < 0.126 eV in this very same mixed hot dark matter model, a value which is very close to the expectations in the inverted hierarchical neutrino mass scenario. Using this same combination of data sets we find the most stringent bound to date on the thermal axion mass, m(a) < 0.529 eV at 95% CL.

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 present up-to-date cosmological bounds on the sum of active neutrino masses as well as on extended cosmological scenarios with additional thermal relics, as thermal axions or sterile neutrino species. Our analyses consider all the current available cosmological data in the beginning of year 2014, including the very recent and most precise baryon acoustic oscillation measurements from the Baryon Oscillation Spectroscopic Survey. In the minimal three-active-neutrino scenario, we find Sigma m(nu) < 0.22 eV at 95% C.L. from the combination of cosmic microwave background (CMB), baryon acoustic oscillation, and Hubble Space Telescope measurements of the Hubble constant. A nonzero value for the sum of the three active neutrino masses of similar to 0.3 eV is significantly favored at more than three standard deviations when adding the constraints on s 8 and Om from the Planck cluster catalog on galaxy number counts. This preference for nonzero thermal relic masses disappears almost completely in both the thermal axion and massive sterile neutrino schemes. Extra light species contribute to the effective number of relativistic degrees of freedom, parametrized via N-eff. We found that when the recent detection of B mode polarization from the BICEP2 experiment is considered, an analysis of the combined CMB data in the framework of LCDM + r models gives N-eff = 3.90 +/- 0.42, suggesting the presence of an extra relativistic relic at more than 95% C.L. from CMB-only data.

Abstract: 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 and for Delta(QCD) = 200 MeV, the full data set implies that the axion mass m(a) = 82.2 +/- 1.1 μeV [corresponding to the Peccei-Quinn symmetry being broken at a scale f(a) = (7.54 +/- 0.10) x 10(10) GeV], or m(a) = 76.6 +/- 2.6 μeV [f(a) = (8.08 +/- 0.27) x 10(10) GeV] when we allow for a nonstandard effective number of relativistic species N-eff. We also find a 2 sigma preference for N-eff > 3.046. The limit on the sum of neutrino masses is Sigma m(v) < 0.25 eV at 95% C.L. for N-eff = 3.046, or Sigma m(v) < 0.47 eV when N-eff is a free parameter. Considering extended scenarios where either the dark energy equation-of-state parameter w, the tensor spectral index n(t), or the running of the scalar index dn(s)/d ln k is allowed to vary does not change significantly the axion mass-energy density constraints. However, in the case of the full data set exploited here, there is a preference for a nonzero tensor index or scalar running, driven by the different tensor amplitudes implied by the Planck and BICEP2 observations. We also study the effect on our estimates of theoretical uncertainties, in particular the imprecise knowledge of the QCD scale Delta(QCD), in the calculation of the temperature-dependent axion mass. We find that in the simplest ADM scenario the Planck + WP data set implies that the axion mass m(a) = 63.7 +/- 1.2 μeV for Delta(QCD) = 400 MeV. We also comment on the possibility that axions do not make up for all the dark matter, or that the contribution of string-produced axions has been grossly underestimated; in that case, the values that we find for the mass can conservatively be considered as lower limits. Dark matter axions with mass in the 60-80 μeV (corresponding to an axion-photon coupling G(a gamma gamma) similar to 10(-14) GeV-1) range can, in principle, be detected by looking for axion-to-photon conversion occurring inside a tunable microwave cavity permeated by a high-intensity magnetic field, and operating at a frequency nu similar or equal to 15-20 GHz. This is out of the reach of current experiments like the axion dark matter experiment (limited to a maximum frequency of a few GHzs), but is, on the other hand, within the reach of the upcoming axion dark matter experiment-high frequency experiment that will explore the 4-40 GHz frequency range and then be sensitive to axion masses up to similar to 160 μeV.