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. Namely, we have included all T2K measurements as of December 2019, the most recent NO$\nu$A antineutrino statistics, and data collected by the Daya Bay and RENO reactor experiments. Concerning the atmospheric and solar sectors, we have also updated our analyses of DeepCore and SNO data, respectively. All in all, these new analyses result in more accurate measurements of $\theta_{13}$, $\theta_{12}$, $\Delta m_{21}^2$ and $|\Delta m_{31}^2|$. The best fit value for the atmospheric angle $\theta_{23}$ lies in the second octant, but first octant solutions remain allowed at $\sim2\sigma$. Regarding CP violation measurements, the preferred value of $\delta$ we obtain is 1.20$\pi$ (1.54$\pi$) for normal (inverted) neutrino mass ordering. These new results should be regarded as extremely robust due to the excellent agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a preference for the normal neutrino mass ordering at the $2.7\sigma$ level. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the significance to $3.1\sigma$ within a conservative approach. A more aggressive data set combination of cosmological observations leads to a stronger preference for normal with respect to inverted mass ordering, at the $3.3\sigma$ level. This cosmological data set provides $2\sigma$ upper limits on the total neutrino mass corresponding to $\sum m_\nu<0.13$ ($0.15$)~eV in the normal (inverted) neutrino mass ordering scenario. These bounds are among the most complete ones in the literature, as they include all currently available neutrino physics inputs.

The ordering of the neutrino masses is a crucial input for a deep understanding of flavor physics, and its determination may provide the key to establish the relationship among the lepton masses and mixings and their analogous properties in the quark sector. The extraction of the neutrino mass ordering is a data-driven field expected to evolve very rapidly in the next decade. In this review, we both analyze the present status and describe the physics of subsequent prospects. Firstly, the different current available tools to measure the neutrino mass ordering are described. Namely, reactor, long-baseline (accelerator and atmospheric) neutrino beams, laboratory searches for beta and neutrinoless double beta decays and observations of the cosmic background radiation and the large scale structure of the universe are carefully reviewed. Secondly, the results from an up-to-date comprehensive global fit are reported: the Bayesian analysis to the 2018 publicly available oscillation and cosmological data sets provides \emph{strong} evidence for the normal neutrino mass ordering versus the inverted scenario, with a significance of 3.5 standard deviations. This preference for the normal neutrino mass ordering is mostly due to neutrino oscillation measurements. Finally, we shall also emphasize the future perspectives for unveiling the neutrino mass ordering. In this regard, apart from describing the expectations from the aforementioned probes, we also focus on those arising from alternative and novel methods, as 21~cm cosmology, core-collapse supernova neutrinos and the direct detection of relic neutrinos.

We present a Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and CMB observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, 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 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 there is a weak-to-moderate preference for logarithmic priors. 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. We also show that the limits on the total neutrino mass $\sum 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.

Neutrino oscillations are a widely observed and well established phenomenon. It is also well known that deviations with respect to flavor conversion probabilities in vacuum arise due to neutrino interactions with matter. In this work, we analyze the impact of new interactions between neutrinos and the dark matter present in the Milky Way on the neutrino oscillation pattern. The dark matter-neutrino interaction is modeled by using an effective coupling proportional to the Fermi constant $G_F$ with no further restrictions on its flavor structure. For the galactic dark matter profile we consider an homogeneous distribution as well as several density profiles, estimating in all cases the size of the interaction required to get an observable effect at different neutrino energies. Our discussion is mainly focused in the PeV neutrino energy range, to be explored in observatories like IceCube and KM3NeT. The obtained results may be interpreted in terms of a light $\mathcal{O}$(sub-eV--keV) or WIMP-like dark matter particle or as a new interaction with a mediator of $\mathcal{O}$(sub-eV--keV) mass.

We present the results of a global neutrino oscillation data analysis within the three-flavour framework. We include latest results from the MINOS long-baseline experiment (including electron neutrino appearance as well as anti-neutrino data), updating all relevant solar (SK II+III), atmospheric (SK I+II+III) and reactor (KamLAND) data. Furthermore, we include a recent re-calculation of the anti-neutrino fluxes emitted from nuclear reactors. These results have important consequences for the analysis of reactor experiments and in particular for the status of the mixing angle $\theta_{13}$. In our recommended default analysis we find from the global fit that the hint for non-zero $\theta_{13}$ remains weak, at 1.8$\sigma$ for both neutrino mass hierarchy schemes. However, we discuss in detail the dependence of these results on assumptions concerning the reactor neutrino analysis.