Abstract: We present the cosmological bounds on the thermal axion mass in an extended cosmological scenario in which the primordial power spectrum of scalar perturbations differs from the usual power-law shape predicted by the simplest inflationary models. The power spectrum is instead modeled by means of a “piecewise cubic Hermite interpolating polynomial” (PCHIP). When using cosmic microwave background measurements combined with other cosmological data sets, the thermal axion mass constraints are degraded only slightly. The addition of the measurements of sigma(8) and Omega(m) from the 2013 Planck cluster catalog on galaxy number counts relaxes the bounds on the thermal axion mass, mildly favoring a similar to 1 eV axion mass, regardless of the model adopted for the primordial power spectrum. However, in general, such a preference disappears if the sum of the three active neutrino masses is also considered as a free parameter in our numerical analyses, due to the strong correlation between the masses of these two hot thermal relics.

Abstract: We present here up-to-date neutrino mass limits exploiting the most recent cosmological data sets. By making use of the cosmic microwave background temperature fluctuation and polarization measurements, supernovae Ia luminosity distances, baryon acoustic oscillation observations and determinations of the growth rate parameter, we are able to set the most constraining bound to date, Sigma m(v) < 0.09 eV at 95% C.L. This very tight limit is obtained without the assumption of any prior on the value of the Hubble constant and highly compromises the viability of the inverted mass ordering as the underlying neutrino mass pattern in nature. The results obtained here further strengthen the case for very large multitracer spectroscopic surveys as unique laboratories for cosmological relics, such as neutrinos: that would be the case of the Dark Energy Spectroscopic Instrument survey and of the Euclid mission.

Abstract: We present robust, model-marginalized limits on both the total neutrino mass (E m1,) and abundances (Neff) to minimize the role of parametrizations, priors and models when extracting neutrino properties from cosmology. The cosmological observations we consider are cosmic microwave background temperature fluctuation and polarization measurements, supernovae Ia luminosity distances, baryon acoustic oscillation observations and determinations of the growth rate parameter from the Data Release 16 of the Sloan Digital Sky Survey IV. The degenerate neutrino mass spectrum (which implies the prior sigma m(1), > 0) is weakly or moderately preferred over the normal and inverted hierarchy possibilities, which imply the priors sigma m(1), > 0.06 and sigma m(1), > 0.1 eV respectively. Concerning the underlying cosmological model, the ACDM minimal scenario is almost always strongly preferred over the possible extensions explored here. The most constraining 95% CL bound on the total neutrino mass in the ACDM + sigma m(1), picture is sigma m(1), < 0.087 eV. The parameter N-eff is restricted to 3.08 +/- 0.17 (68% CL) in the ACDM + Neff model. These limits barely change when considering the ACDM + sigma m(1), + Neff scenario. Given the robustness and the strong constraining power of the cosmological measurements employed here, the model -marginalized posteriors obtained considering a large spectra of nonminimal cosmologies are very close to the previous bounds, obtained within the ACDM framework in the degenerate neutrino mass spectrum. Future cosmological measurements may improve the current Bayesian evidence favoring the degenerate neutrino mass spectra, challenging therefore the consistency between cosmological neutrino mass bounds and oscillation neutrino measurements, and potentially suggesting a more complicated cosmological model and/or neutrino sector.

Abstract: Minimal dark energy models, described by the same number of free parameters of the standard cosmological model with cold dark matter plus a cosmological constant to parametrize the dark energy component, constitute very appealing scenarios which may solve long-standing, pending tensions. On the one hand, they alleviate significantly the tension between cosmological observations and the presence of one sterile neutrino motivated by the short-baseline anomalies: we obtain a 95% CL cosmological bound on the mass of a fully thermalized fourth sterile neutrino (N-eff = 4) equal to m(s) < 0.65(1.3) eV within the Phenomenologically Emergent Dark Energy (PEDE) and Vacuum Metamorphosis (VM) scenarios under consideration. Interestingly, these limits are in agreement with the observations at short-baseline experiments, and the PEDE scenario is favored with respect to the Lambda CDM case when the full data combination is considered. On the other hand, the Hubble tension is satisfactorily solved in almost all the minimal dark energy schemes explored here. These phenomenological scenarios may therefore shed light on differences arising from near and far Universe probes, and also on discrepancies between cosmological and laboratory sterile neutrino searches.

Abstract: We present a full Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and Cosmic Microwave Background observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, namely Normal Ordering or Inverted Ordering, with current data. Two possible parametrizations (three neutrino masses versus the lightest neutrino mass plus the two oscillation mass splittings) and priors (linear versus logarithmic) are exhaustively examined. We find that the preference for NO is only driven by neutrino oscillation data. Moreover, the values of the Bayes factor indicate that the evidence for NO is strong only when the scan is performed over the three neutrino masses with logarithmic priors; for every other combination of parameterization and prior, the preference for NO is only weak. As a by-product of our Bayesian analyses, we are able to (a) compare the Bayesian bounds on the neutrino mixing parameters to those obtained by means of frequentist approaches, finding a very good agreement; (b) determine that the lightest neutrino mass plus the two mass splittings parametrization, motivated by the physical observables, is strongly preferred over the three neutrino mass eigenstates scan and (c) find that logarithmic priors guarantee a weakly-to-moderately more efficient sampling of the parameter space. These results establish the optimal strategy to successfully explore the neutrino parameter space, based on the use of the oscillation mass splittings and a logarithmic prior on the lightest neutrino mass, when combining neutrino oscillation data with cosmology and neutrinoless double beta decay. We also show that the limits on the total neutrino mass Sigma m(nu) can change dramatically when moving from one prior to the other. These results have profound implications for future studies on the neutrino mass ordering, as they crucially state the need for self-consistent analyses which explore the best parametrization and priors, without combining results that involve different assumptions.