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νA measurements, as reported in the Neutrino 2020 conference. All in all, these new analyses result in more accurate measurements of θ₁₃, θ₁₂, $$ \Delta {m}_{21}^2 $$Δ m 21 2 and $$ \left|\Delta {m}_{31}^2\right| $$Δ m 31 2 . The best fit value for the atmospheric angle θ₂₃ lies in the second octant, but first octant solutions remain allowed at ∼ 2.4σ. Regarding CP violation measurements, the preferred value of δ we obtain is 1.08π (1.58π) for normal (inverted) neutrino mass ordering. The global analysis still prefers normal neutrino mass ordering with 2.5σ 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σ. 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σ 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σ. This very same cosmological data set provides 2σ upper limits on the total neutrino mass corresponding to Σm_ν< 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.

Here, we update our previous global fit of neutrino oscillations by including the recent results that have appeared since the Neutrino 2012 conference. These include the measurements of reactor antineutrino disappearance reported by Daya Bay and RENO, together with latest T2K and MINOS data including both disappearance and appearance channels. We also include the revised results from the third solar phase of Super-Kamiokande, SK-III, as well as new solar results from the fourth phase of Super-Kamiokande, SK-IV. We find that the preferred global determination of the atmospheric angle ${\theta }_{23}$ is consistent with maximal mixing. We also determine the impact of the new data upon all the other neutrino oscillation parameters with an emphasis on the increasing sensitivity to the $CP$ phase, thanks to the interplay between accelerator and reactor data. In the Appendix, we present the updated results obtained after the inclusion of new reactor data presented at the Neutrino 2014 conference. We discuss their impact on the global neutrino analysis.

Current neutrino experiments are measuring the neutrino mixing parameters with an unprecedented accuracy. The upcoming generation of neutrino experiments will be sensitive to subdominant neutrino oscillation effects that can in principle give information on the yet-unknown neutrino parameters: the Dirac CP-violating phase in the PMNS mixing matrix, the neutrino mass ordering and the octant of θ₂₃. Determining the exact values of neutrino mass and mixing parameters is crucial to test various neutrino models and flavor symmetries that are designed to predict these neutrino parameters. In the first part of this review, we summarize the current status of the neutrino oscillation parameter determination. We consider the most recent data from all solar neutrino experiments and the atmospheric neutrino data from Super-Kamiokande, IceCube, and ANTARES. We also implement the data from the reactor neutrino experiments KamLAND, Daya Bay, RENO, and Double Chooz as well as the long baseline neutrino data from MINOS, T2K, and NOνA. If in addition to the standard interactions, neutrinos have subdominant yet-unknown Non-Standard Interactions (NSI) with matter fields, extracting the values of these parameters will suffer from new degeneracies and ambiguities. We review such effects and formulate the conditions on the NSI parameters under which the precision measurement of neutrino oscillation parameters can be distorted. Like standard weak interactions, the non-standard interaction can be categorized into two groups: Charged Current (CC) NSI and Neutral Current (NC) NSI. Our focus will be mainly on neutral current NSI because it is possible to build a class of models that give rise to sizeable NC NSI with discernible effects on neutrino oscillation. These models are based on new U(1) gauge symmetry with a gauge boson of mass ≲ 10 MeV. The UV complete model should be of course electroweak invariant which in general implies that along with neutrinos, charged fermions also acquire new interactions on which there are strong bounds. We enumerate the bounds that already exist on the electroweak symmetric models and demonstrate that it is possible to build viable models avoiding all these bounds. In the end, we review methods to test these models and suggest approaches to break the degeneracies in deriving neutrino mass parameters caused by NSI.