Abstract: In models with dark matter made of particles with keV masses, such as a sterile neutrino, small-scale density perturbations are suppressed, delaying the period at which the lowest mass galaxies are formed and therefore shifting the reionization processes to later epochs. In this study, focusing on Warm Dark Matter (WDM) with masses close to its present lower bound, i.e., around the 3. keV region, we derive constraints from galaxy luminosity functions, the ionization history and the Gunn-Peterson effect. We show that even if star formation efficiency in the simulations is adjusted to match the observed UV galaxy luminosity functions in both CDM and WDM models, the full distribution of Gunn-Peterson optical depth retains the strong signature of delayed reionization in the WDM model. However, until the star formation and stellar feedback model used in modern galaxy formation simulations is constrained better, any conclusions on the nature of dark matter derived from reionization observables remain model-dependent.

Abstract: We propose here a novel method which singles out the a priori unavoidable dependence on the underlying cosmological model when extracting parameter constraints, providing robust limits which only depend on the considered dataset. Interestingly, when dealing with several possible cosmologies and interpreting the Bayesian preference in terms of the Gaussian statistical evidence, the preferred model is much less favored than when only two cases are compared. As a working example, we apply our approach to the cosmological neutrino mass bounds, which play a fundamental role not only in establishing the contribution of relic neutrinos to the dark matter of the Universe but also in the planning of future experimental searches of the neutrino character and of the neutrino mass ordering.

Abstract: We present an updated global fit of neutrino oscillation data in the simplest three-neutrino framework. In the present study we include up-to-date analyses from a number of experiments. Concerning the atmospheric and solar sectors, besides the data considered previously, we give updated analyses of IceCube DeepCore and Sudbury Neutrino Observatory data, respectively. We have also included the latest electron antineutrino data collected by the Daya Bay and RENO reactor experiments, and the long-baseline T2K and NO nu A measurements, as reported in the Neutrino 2020 conference. All in all, these new analyses result in more accurate measurements of theta (13), theta (12), Delta m212 and Delta m312. The best fit value for the atmospheric angle theta (23) lies in the second octant, but first octant solutions remain allowed at similar to 2.4 sigma. Regarding CP violation measurements, the preferred value of delta we obtain is 1.08 pi (1.58 pi) for normal (inverted) neutrino mass ordering. The global analysis still prefers normal neutrino mass ordering with 2.5 sigma statistical significance. This preference is milder than the one found in previous global analyses. These new results should be regarded as robust due to the agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a weak/moderate preference for the normal neutrino mass ordering of 2.00 sigma. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the preference to 2.68 sigma within a conservative approach. A more aggressive data set combination of cosmological observations leads to a similar preference for normal with respect to inverted mass ordering, namely 2.70 sigma. This very same cosmological data set provides 2 sigma upper limits on the total neutrino mass corresponding to Sigma m(nu)< 0.12 (0.15) eV in the normal (inverted) neutrino mass ordering scenario. The bounds on the neutrino mixing parameters and masses presented in this up-to-date global fit analysis include all currently available neutrino physics inputs.