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: The Phenomenologically Emergent Dark Energy model, a dark energy model with the same number of free parameters as the flat Lambda CDM, has been proposed as a working example of a minimal model which can avoid the current cosmological tensions. A straightforward question is whether or not the inclusion of massive neutrinos and extra relativistic species may spoil such an appealing phenomenological alternative. We present the bounds on M-nu and N-eff and comment on the long standing H-0 and sigma(8) tensions within this cosmological framework with a wealth of cosmological observations. Interestingly, we find, at 95% confidence level, and with the most complete set of cosmological observations, M-nu similar to 0.21(-0.14)(+0.15) eV and N-eff = 3.03 +/- 0.32 i.e. an indication for a non-zero neutrino mass with a significance above 2 sigma. The well known Hubble constant tension is considerably easened, with a significance always below the 2 sigma level. (C) 2020 Elsevier B.V. All rights reserved.