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: We illustrate how recently improved low-redshift cosmological measurements can tighten constraints on neutrino properties. In particular we examine the impact of the assumed cosmological model on the constraints. We first consider the new HST H-0 = 74.2 +/- 3.6 measurement by Riess et al. (2009) and the sigma(8)(Omega(m)/0.25)(0.41) = 0.832 +/- 0.033 constraint from Rozo et al. (2009) derived from the SDSS maxBCG Cluster Catalog. In a ACDM model and when combined with WMAP5 constraints, these low-redshift measurements constrain Sigma m(v) < 0.4 eV at the 95% confidence level. This bound does not relax when allowing for the running of the spectral index or for primordial tensor perturbations. When adding also Supernovae and BAO constraints, we obtain a 95% upper limit of Sigma m(v) < 0.3eV. We test the sensitivity of the neutrino mass constraint to the assumed expansion history by both allowing a dark energy equation of state parameter w not equal -1 and by studying a model with coupling between dark energy and dark matter, which allows for variation in w, Omega(k), and dark coupling strength xi. When combining CMB, H-0 and the SDSS LRG halo power spectrum from Reid et al. 2009, we find that in this very general model, Sigma m(v) < 0.51 eV with 95% confidence. If we allow the number of relativistic species N-rel to vary in a ACDM model with Sigma m(v) = 0, we find N-rel = 3.76(-0.68)(+0.63)(+1.38 -1.21) for the 68% and 95% confidence intervals. We also report prior-independent constraints, which are in excellent agreement with the Bayesian constraints.

Abstract: Models involving an interaction between the dark matter and the dark energy sectors have been proposed to alleviate the long-standing Hubble constant tension. In this paper, we analyse whether the constraints and potential hints obtained for these interacting models remain unchanged when using simulated Planck data. Interestingly, our simulations indicate that a dangerous fake detection for a non-zero interaction among the dark matter and the dark energy fluids could arise when dealing with current cosmic microwave background (CMB) Planck measurements alone. The very same hypothesis is tested against future CMB observations, finding that only cosmic variance limited polarization experiments, such as PICO or PRISM, could be able to break the existing parameter degeneracies and provide reliable cosmological constraints. This paper underlines the extreme importance of confronting the results arising from data analyses with those obtained with simulations when extracting cosmological limits within exotic cosmological scenarios.

Abstract: One of the targets of future cosmic microwave background (CMB) and baryon acoustic oscillation measurements is to improve the current accuracy in the neutrino sector and reach a much better sensitivity on extra dark radiation in the early Universe. In this paper, we study how these improvements can be translated into constraining power for well-motivated extensions of the standard model of elementary particles that involve axions thermalized before the quantum chromodynamics (QCD) phase transition by scatterings with gluons. Assuming a fiducial Lambda cold dark matter cosmological model, we simulate future data for Stage-IV CMB-like and Dark Energy Spectroscopic Instrument (DESI)-like surveys and analyse a mixed scenario of axion and neutrino hot dark matter. We further account also for the effects of these QCD axions on the light element abundances predicted by big bang nucleosynthesis. The most constraining forecasted limits on the hot relic masses are m(a) less than or similar to 0.92 eV and n-ary sumation m(nu) less than or similar to 0.12 eV at 95 per cent Confidence Level, showing that future cosmic observations can substantially improve the current bounds, supporting multimessenger analyses of axion, neutrino, and primordial light element properties.

Abstract: We consider the phenomenological implications of the violation of the Pauli exclusion principle for neutrinos, focusing on cosmological observables such as the spectrum of Cosmic Microwave Background anisotropies, Baryon Acoustic Oscillations and the primordial abundances of light elements. Neutrinos that behave (at least partly) as bosonic particles have a modified equilibrium distribution function that implies a different influence on the evolution of the Universe that, in the case of massive neutrinos, can not be simply parametrized by a change in the effective number of neutrinos. Our results show that, despite the precision of the available cosmological data, only very weak bounds can be obtained on neutrino statistics, disfavouring a more bosonic behaviour at less than 2 sigma.