On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

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$.

The Short BaseLine (SBL) neutrino oscillation anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such a neutrino is incompatible with cosmological data, in particular observations of the Cosmic Microwave Background (CMB) anisotropies. However, this conclusion can change by invoking new physics. One possibility is to introduce a secret interaction in the sterile neutrino sector mediated by a light pseudoscalar. In this pseudoscalar model, CMB data prefer a sterile neutrino mass that is fully compatible with the mass ranges suggested by SBL anomalies. In addition, this model predicts a value of the Hubble parameter which is completely consistent with local measurements.

The theory and phenomenology of light sterile neutrinos at the eV mass scale is reviewed. The reactor, Gallium and LSND anomalies are briefly described and interpreted as indications of the existence of short-baseline oscillations which require the existence of light sterile neutrinos. The global fits of short-baseline oscillation data in 3+1 and 3+2 schemes are discussed, together with the implications for beta-decay and neutrinoless double-beta decay. The cosmological effects of light sterile neutrinos are briefly reviewed and the implications of existing cosmological data are discussed. The review concludes with a summary of future perspectives.

We perform a cosmological analysis in which we allow the primordial power spectrum of scalar perturbations to assume a shape that is different from the usual power-law predicted by the simplest models of cosmological inflation. We parameterize the free primordial power spectrum with a "piecewise cubic Hermite interpolating polynomial" (PCHIP). We consider a 3+1 neutrino mixing model with a sterile neutrino having a mass at the eV scale, which can explain the anomalies observed in short-baseline neutrino oscillation experiments. We find that the freedom of the primordial power spectrum allows to reconcile the cosmological data with a fully thermalized sterile neutrino in the early Universe. Moreover, the cosmological analysis gives us some information on the shape of the primordial power spectrum, which presents a feature around the wavenumber $k=0.002\,\text{Mpc}^{-1}$.

The recent discovery of B-modes in the polarization pattern of the Cosmic Microwave Background by the BICEP2 experiment has important implications for neutrino physics. We revisit cosmological bounds on light sterile neutrinos and show that they are compatible with all current cosmological data provided that the mass is relatively low. Using CMB data, including BICEP-2, we find an upper bound of $m_s < 0.85$ eV ($2\sigma$ Confidence Level). This bound is strengthened to 0.48 eV when HST measurements of $H_0$ are included. However, the inclusion of SZ cluster data from the Planck mission and weak gravitational measurements from the CFHTLenS project favours a non-zero sterile neutrino mass of $0.44^{+0.11}_{-0.16}$ eV. Short baseline neutrino oscillations, on the other hand, indicate a new mass state around 1.2 eV. This mass is highly incompatible with cosmological data if the sterile neutrino is fully thermalised ($\Delta \chi^2>10$). However, if the sterile neutrino only partly thermalises it can be compatible with all current data, both cosmological and terrestrial.

We introduce a cosmological invisible decay of the sterile neutrino with the eV-scale mass indicated by short-baseline neutrino oscillation experiments in order to allow its full thermalization in the early Universe. We show that the fit of the cosmological data is practically as good as the fit obtained with a stable sterile neutrino without mass constraints, which has been recently considered by several authors for the explanation of the observed suppression of small-scale matter density fluctuations and for a solution of the tension between the Planck and BICEP2 measurements of the tensor to scalar ratio of large-scale fluctuations. Moreover, the extra relativistic degree of freedom corresponding to a fully thermalized sterile neutrino is correlated with a larger value of the Hubble constant, which is in agreement with local measurements.

We analyze the most recent cosmological data, including Planck, taking into account the possible existence of a sterile neutrino with a mass at the eV scale indicated by short-baseline neutrino oscillations data in the 3+1 framework. We show that the contribution of local measurements of the Hubble constant induces an increase of the value of the effective number of relativistic degrees of freedom above the Standard Model value, giving an indication in favor of the existence of sterile neutrinos and their contribution to dark radiation. Furthermore, the measurements of the local galaxy cluster mass distribution favor the existence of sterile neutrinos with eV-scale masses, in agreement with short-baseline neutrino oscillations data. In this case there is no tension between cosmological and short-baseline neutrino oscillations data, but the contribution of the sterile neutrino to the effective number of relativistic degrees of freedom is likely to be smaller than one. Considering the Dodelson-Widrow and thermal models for the statistical cosmological distribution of sterile neutrinos, we found that in the Dodelson-Widrow model there is a slightly better compatibility between cosmological and short-baseline neutrino oscillations data and the required suppression of the production of sterile neutrinos in the early Universe is slightly smaller.

We present the results of global analyses of short-baseline neutrino oscillation data in 3+1, 3+2 and 3+1+1 neutrino mixing schemes. We show that the data do not allow us to abandon the simplest 3+1 scheme in favor of the more complex 3+2 and 3+1+1 schemes. We present the allowed region in the 3+1 parameter space, which is located at $\Delta{m}^2_{41}$ between 0.82 and 2.19 $\text{eV}^2$ at $3\sigma$. The case of no oscillations is disfavored by about $6\sigma$, which decreases dramatically to about $2\sigma$ if the LSND data are not considered. Hence, new high-precision experiments are needed to check the LSND signal.

We discuss the implications for short-baseline electron neutrino disappearance in the 3+1 mixing scheme of the recent Troitsk bounds on the mixing of a neutrino with mass between 2 and 100 eV. Considering the Troitsk data in combination with the results of short-baseline nu_e and antinu_e disappearance experiments, which include the reactor and Gallium anomalies, we derive a 2 sigma allowed range for the effective neutrino squared-mass difference between 0.85 and 43 eV^2. The upper bound implies that it is likely that oscillations in distance and/or energy can be observed in radioactive source experiments. It is also favorable for the ICARUS@CERN experiment, in which it is likely that oscillations are not washed-out in the near detector. We discuss also the implications for neutrinoless double-beta decay.

We present a complete update of the analysis of electron neutrino and antineutrino disappearance experiments in terms of neutrino oscillations in the framework of 3+1 neutrino mixing, taking into account the Gallium anomaly, the reactor anomaly, solar neutrino data and nu_e-C scattering data. We discuss the implications of a recent 71Ga(3He,3H)71Ge measurement which give information on the neutrino cross section in Gallium experiments. We discuss the solar bound on active-sterile mixing and present our numerical results. We discuss the connection between the results of the fit of neutrino oscillation data and the heavy neutrino mass effects in beta-decay experiments (considering new Mainz data) and neutrinoless double-beta decay experiments (considering the recent EXO results).

There are hints, both from cosmology and from neutrino oscillation experiments, that one or two sterile neutrinos with eV masses are favored. We consider the implications of combining data from short baseline neutrino experiments with cosmological observations. In particular we consider cosmological models where the equation of state of dark energy and the curvature are allowed to vary. We find that the consistency of results from short baseline experiments and cosmology is dependent on the cosmological data sets that are used. When the two are consistent, we find that the inclusion of information from short baseline experiments only mildly affects the preferred values of the cosmological model parameters. A universe with a cosmological constant type dark energy and zero curvature is still fully consistent with the data we use. We also report that a thermalized sterile neutrino in the 10 eV range is fully consistent with a wide range of cosmological data as long as information on the matter power spectrum in not taken into account.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We investigate the possible effects of short-baseline antinu_e disappearance implied by the reactor antineutrino anomaly on the Double-Chooz determination of theta_{13} through the normalization of the initial antineutrino flux with the Bugey-4 measurement. We show that the effects are negligible and the value of theta_{13} obtained by the Double-Chooz collaboration is accurate only if Delta m^2_{41} is larger than about 3 eV^2. For smaller values of Delta m^2_{41} the short-baseline oscillations are not fully averaged at Bugey-4 and the uncertainties due to the reactor antineutrino anomaly can be of the same order of magnitude of the intrinsic Double-Chooz uncertainties.

We present an update of our analysis of short-baseline neutrino oscillation data in the framework of 3+1 neutrino mixing taking into account the recent update of MiniBooNE antineutrino data and the recent results of the MINOS search for nu_mu disappearance into sterile neutrinos (the more complicated 3+2 neutrino mixing is not needed since the CP-violating difference between MiniBooNE neutrino and antineutrino data has diminished). The results of our fits of short-baseline neutrino oscillation data including the MiniBooNE low-energy anomaly (now present both in the neutrino and antineutrino data) leads to a strong tension between appearance and disappearance data. Hence, it seems likely that the low-energy anomaly is not due to nu_mu -> nu_e transitions. Excluding the MiniBooNE low-energy anomaly, appearance and disappearance data are marginally compatible. The global analysis has the best-fit point at Delta m^2_{41} about 5.6 eV^2, which is rather large in comparison with cosmological bounds, but there are three regions within 1 sigma at Delta m^2_{41} about 1.6, 1.2, 0.91 eV^2. We also show that the data on the Gallium neutrino anomaly favor values of Delta m^2_{41} larger than about 1 eV^2.

We present an upgrade of the 3+1 global fit of short-baseline neutrino oscillation data obtained with the addition of KARMEN and LSND nu_e-Carbon scattering data. We discuss the implications for the measurements of the effective neutrino mass in beta-decay and neutrinoless double-beta-decay experiments. We find respective predicted ranges of about 0.1-0.7 eV and 0.01-0.1 eV.

The new data release of OPERA - CNGS experiment, obtained with a shorter spill of protons, confirms the tachyionic behavior expected from the phenomenological model of a Majorana neutrino with a fictitious imaginary mass term acquired during the propagation in the Earth's crust, recently presented by us. We performed numerical simulations of neutrino event detections to compare the properties of these Majorana tachyons with the new OPERA results, finding a good agreement. The possibility of spin-to orbital angular momentum conversion that is expected to give a negative squared mass in a medium, is also briefly discussed.

From the data release of OPERA - CNGS experiment, and publicly announced on 23 September 2011, we cast a phenomenological model based on a Majorana neutrino state carrying a fictitious imaginary mass term, already discussed by Majorana in 1932. This mass term can be induced by the interaction with the matter of the Earth's crust during the 735 Km travel. Within the experimental errors, we prove that the model fits with OPERA, MINOS and supernova SN1987a data. Possible violations to Lorentz invariance due to quantum gravity effects have been considered.

We present the results of fits of short-baseline neutrino oscillation data in 3+1 and 3+2 neutrino mixing schemes. In spite of the presence of a tension in the interpretation of the data, 3+1 neutrino mixing is attractive for its simplicity and for the natural correspondence of one new entity (a sterile neutrino) with a new effect (short-baseline oscillations). The allowed regions in the oscillation parameter space can be tested in near-future experiments. In the framework of 3+2 neutrino mixing there is less tension in the interpretation of the data, at the price of introducing a second sterile neutrino. Moreover, the improvement of the parameter goodness of fit is mainly a statistical effect due to an increase of the number of parameters. The CP violation in short-baseline experiments allowed in 3+2 neutrino mixing can explain the positive antinu_mu -> antinu_e signal and the negative nu_mu -> nu_e measurement in the MiniBooNE experiment. For the CP-violating phase we obtained two minima of the marginal chi^2 close to the two values where CP-violation is maximal.

We analyze the LSND, KARMEN and MiniBooNE data on short-baseline antinu_mu -> antinu_e oscillations and the data on short-baseline antinu_e disappearance obtained in the Bugey-3 and CHOOZ reactor experiments in the framework of 3+1 antineutrino mixing, taking into account the MINOS observation of long-baseline antinu_mu disappearance and the KamLAND observation of very-long-baseline antinu_e disappearance. We show that the fit of the data implies that the short-baseline disappearance of antinu_mu is relatively large. We obtain a prediction of an effective amplitude sin^2 2 theta_{mu mu} >~ 0.1 for short-baseline antinu_mu disappearance generated by 0.2 <~ Delta m^2 <~ 1 eV^2, which could be measured in future experiments.

We analyzed the electron neutrino data of the Gallium radioactive source experiments and the electron antineutrino data of the reactor Bugey and Chooz experiments in terms of neutrino oscillations. We found a hint of a CPT-violating asymmetry of the effective neutrino and antineutrino mixing angles.

We analyze the recent results of the MiniBooNE short-baseline experiment on antinu_mu -> antinu_e oscillations in a minimal model-independent framework of antineutrino mixing in conjunction with the positive LSND signal and the negative KARMEN measurements. We show that the data of the three short-baseline antinu_mu -> antinu_e experiments are compatible. Taking into account also the model-independent constraints due to the limits on short-baseline antinu_e disappearance obtained in reactor antineutrino experiments, we find that the favored region of the effective oscillation parameters lies within 0.002 <~ sin^2 2 theta <~ 0.05 and 0.2 <~ Delta m^2 <~ 2 eV^2.

We consider the interpretation of the MiniBooNE low-energy anomaly and the Gallium radioactive source experiments anomaly in terms of short-baseline electron neutrino disappearance in the framework of 3+1 four-neutrino mixing schemes. The separate fits of MiniBooNE and Gallium data are highly compatible, with close best-fit values of the effective oscillation parameters Delta m^2 and sin^2 2 theta. The combined fit gives Delta m^2 >~ 0.1 eV^2 and 0.11 < sin^2 2 theta < 0.48 at 2 sigma. We consider also the data of the Bugey and Chooz reactor antineutrino oscillation experiments and the limits on the effective electron antineutrino mass in beta-decay obtained in the Mainz and Troitsk Tritium experiments. The fit of the data of these experiments limits the value of sin^2 2 theta below 0.10 at 2 sigma. Considering the tension between the neutrino MiniBooNE and Gallium data and the antineutrino reactor and Tritium data as a statistical fluctuation, we perform a combined fit which gives Delta m^2 \simeq 2 eV and 0.01 < sin^2 2 theta < 0.13 at 2 sigma. Assuming a hierarchy of masses m_1, m_2, m_3 << m_4, the predicted contributions of m_4 to the effective neutrino masses in beta-decay and neutrinoless double-beta-decay are, respectively, between about 0.06 and 0.49 and between about 0.003 and 0.07 eV at 2 sigma. We also consider the possibility of reconciling the tension between the neutrino MiniBooNE and Gallium data and the antineutrino reactor and Tritium data with different mixings in the neutrino and antineutrino sectors. We find a 2.6 sigma indication of a mixing angle asymmetry.

In this review we present the main features of the current status of neutrino physics. After a review of the theory of neutrino mixing and oscillations, we discuss the current status of solar and atmospheric neutrino oscillation experiments. We show that the current data can be nicely accommodated in the framework of three-neutrino mixing. We discuss also the problem of the determination of the absolute neutrino mass scale through Tritium beta-decay experiments and astrophysical observations, and the exploration of the Majorana nature of massive neutrinos through neutrinoless double-beta decay experiments. Finally, future prospects are briefly discussed.

A Neutrino Unbound gem (http://www.to.infn.it/~giunti/NU). Essential information (formulas, figures, tables, references, etc.) on solar neutrinos.