Sterile neutrinos with a mass in the eV range have been invoked as a possible explanation of a variety of short baseline (SBL) neutrino oscillation anomalies. However, if one considers neutrino oscillations between active and sterile neutrinos, such neutrinos would have been fully thermalised in the early universe, and would be therefore in strong conflict with cosmological bounds. In this study we first update cosmological bounds on the mass and energy density of eV-scale sterile neutrinos. We then perform an updated study of a previously proposed model in which the sterile neutrino couples to a new light pseudoscalar degree of freedom. Consistently with previous analyses, we find that the model provides a good fit to all cosmological data and allows the high value of $H_0$ measured in the local universe to be consistent with measurements of the cosmic microwave background. However, new high $\ell$ polarisation data constrain the sterile neutrino mass to be less than approximately 1~eV in this scenario. Finally, we combine the cosmological bounds on the pseudoscalar model with a Bayesian inference analysis of SBL data and conclude that only a sterile mass in a narrow range around 1~eV remains consistent with both cosmology and SBL data.

In various extensions of the Standard Model of particle physics, and intriguingly even in the three-generation Standard Model without neutrino masses, neutrinos are allowed to have very tiny electric charges. We revisit the emergence of such charges and constraints on them, adding to the repertoire of limits that exist in literature, including new indirect limits on the charges of the muon and tau neutrinos respectively. We also point out the flavor-universal bounds on neutrino charges, motivated by the weak gravity conjecture and based on the observation of neutron stars with very large magnetic fields. The latter set of limits are relatively model independent, relying only on the hypothesis that gravity must be the weakest force and that any charged neutrino must have a coupling to the photon. The weak gravity conjecture gives a lower bound on the possible neutrino charge. Astrophysical observations of Magnetars meanwhile give an upper bound, comparable to recent limits derived from the reactor neutrinos, considering elastic neutrino-electron scattering, and experiments constraining neutrino magnetic moments.

We consider possible effects of neutrino electric charge (millicharge) and charge radius on the neutrino-atom interaction processes such as (i) atomic ionization by neutrino impact and (ii) coherent elastic neutrino-nucleus scattering. The bounds on the neutrino millicharge and charge radius that follow from, respectively, the GEMMA and COHERENT experiments are presented and discussed.

A number of anomalous results in short-baseline oscillation may hint at the existence of one or more light sterile neutrino states in the eV mass range and have triggered a wave of new experimental efforts to search for a definite signature of oscillations between active and sterile neutrino states. The present paper aims to provide a comprehensive review on the status of light sterile neutrino searches in mid-2019: we discuss not only the basic experimental approaches and sensitivities of reactor, source, atmospheric, and accelerator neutrino oscillation experiments but also the complementary bounds arising from direct neutrino mass experiments and cosmological observations. Moreover, we review current results from global oscillation analyses that include the constraints set by running reactor and atmospheric neutrino experiments. They permit to set tighter bounds on the active-sterile oscillation parameters but as yet are not able to provide a definite conclusion on the existence of eV-scale sterile neutrinos.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

The current status of the phenomenology of short-baseline neutrino oscillations induced by light sterile neutrinos in the framework of 3+1 mixing is reviewed.

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 indications in favor of the existence of light sterile neutrinos at the eV scale found in short-baseline neutrino oscillation experiments is reviewed. The future perspectives of short-baseline neutrino oscillation experiments and the connections with beta-decay measurements of the neutrino masses and with neutrinoless double-beta decay experiments are discussed.

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 derive the relation between the amplitudes of short-baseline appearance and disappearance oscillations in 3+$N_{s}$ neutrino mixing schemes which is the origin of the appearance-disappearance tension that is found from the analysis of the existing data in any 3+$N_{s}$ neutrino mixing scheme. We illustrate the power of the relation to reveal the appearance-disappearance tension in the cases of 3+1 and 3+2 mixing using the results of global fits of short-baseline neutrino oscillation data.

We review the theory and phenomenology of neutrino electromagnetic interactions, which give us powerful tools to probe the physics beyond the Standard Model. After a derivation of the general structure of the electromagnetic interactions of Dirac and Majorana neutrinos in the one-photon approximation, we discuss the effects of neutrino electromagnetic interactions in terrestrial experiments and in astrophysical environments. We present the experimental bounds on neutrino electromagnetic properties and we confront them with the predictions of theories beyond the Standard Model.

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 review the results of global analyses of short-baseline neutrino oscillation data in 3+1, 3+2 and 3+1+1 neutrino mixing schemes.

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.

Cosmology and short baseline neutrino oscillation data both hint at the existence of light sterile neutrinos with masses in the 1 eV range. Here we perform a detailed analysis of the sterile neutrino scenario using both cosmological and SBL data. We have additionally considered the possibility that the extra neutrino degrees of freedom are not fully thermalised in the early universe. Even when analyzing only cosmological data we find a preference for the existence of massive sterile neutrinos in both (3+1) and (3+2) scenarios, and with the inclusion of SBL data the evidence is formally at the 3.3sigma level in the case of a (3+1) model. Interestingly, cosmological and SBL data both point to the same mass scale of approximately 1 eV. In the (3+1) framework WMAP9+SPT provide a value of the sterile mass eigenstate m_4 = (1.72 \pm 0.65) eV: this result is strenghtened by adding the prior from SBL posterior to m_4 = (1.27 \pm 0.12) eV (m_4 = (1.23 \pm 0.13) eV when SDSS is also considered in the cosmological analysis). In the (3+2) scheme, two additional, non--fully thermalized, neutrinos are compatible with the whole set of cosmological and SBL data, leading to mass values of m_4 = (0.95 \pm 0.30) eV and m_5 = (1.59 \pm 0.49) eV. The inclusion of Planck data does not change our considerations about the mass scale; concerning the extra neutrino degrees of freedom, invoking a partial thermalisation the 3+1 model is still consistent with the latest data.

Light sterile neutrinos can be excited by oscillations with active neutrinos in the early universe. Their properties can be constrained by their contribution as extra-radiation, parameterized in terms of the effective number of neutrino species N_ eff, and to the universe energy density today \Omega_\nu h^2. Both these parameters have been measured to quite a good precision by the Planck satellite experiment. We use this result to update the bounds on the parameter space of (3+1) sterile neutrino scenarios, with an active-sterile neutrino mass squared splitting in the range (10^{-5} - 10^2 ) eV^2. We consider both normal and inverted mass orderings for the active and sterile states. For the first time we take into account the possibility of two non-vanishing active-sterile mixing angles. We find that the bounds are more stringent than those obtained in laboratory experiments. This leads to a strong tension with the short-baseline hints of light sterile neutrinos. In order to relieve this disagreement, modifications of the standard cosmological scenario, e.g. large primordial neutrino asymmetries, are required.

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

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.

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.

Recent results from short--baseline neutrino oscillation experiments and Cosmic Microwave Background anisotropy measurements suggest the presence of additional sterile neutrinos. In this paper we properly combine these data sets to derive bounds on the sterile neutrino masses in the 3+1 and 3+2 frameworks, finding a potentially good agreement between the two datasets. However, when galaxy clustering is included in the analysis a tension between the oscillation and cosmological data is clearly present.

In this review we discuss the main theoretical aspects and experimental effects of neutrino electromagnetic properties. We start with a general description of the electromagnetic form factors of Dirac and Majorana neutrinos. Then, we discuss the theory and phenomenology of the magnetic and electric dipole moments, summarizing the experimental results and the theoretical predictions. We discuss also the phenomenology of a neutrino charge radius and radiative decay. Finally, we describe the theory of neutrino spin and spin-flavor precession in a transverse magnetic field and we summarize its phenomenological applications.

In this brief review we discuss the generation of Majorana neutrino masses through the see-saw mechanism, the theory of neutrinoless double-beta decay, the implications of neutrino oscillation data for the effective Majorana mass, taking into account the recent Daya Bay measurement of theta_13, and the interpretation of the results of neutrinoless double-beta decay experiments.

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 indications in favor of short-baseline neutrino oscillations, which require the existence of one or more sterile neutrinos, are reviewed. In the framework of 3+1 neutrino mixing, which is the simplest extension of the standard three-neutrino mixing which can partially explain the data, there is a strong tension in the interpretation of the data, mainly due to an incompatibility of the results of appearance and disappearance experiments. In the framework of 3+2 neutrino mixing, CP violation in short-baseline experiments can explain the difference between MiniBooNE neutrino and antineutrino data, but the tension between the data of appearance and disappearance experiments persists because the short-baseline disappearance of electron antineutrinos and muon neutrinos compatible with the LSND and MiniBooNE antineutrino appearance signal has not been observed.

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.

After a brief review of the results of solar, atmospheric and long-baseline neutrino oscillation experiments which led to the current three-neutrino mixing paradigm, we discuss indications of neutrino oscillation experiments in favor of short-baseline oscillations which require the existence of one or more sterile neutrinos. We show that the simplest possibility of existence of one sterile neutrino is not enough to fit all data of short-baseline neutrino oscillation experiments because of two tensions: a tension between neutrino and antineutrino data and a tension between appearance and disappearance data. The tension between neutrino and antineutrino data is eliminated with the addition of a second sterile neutrino which allows CP-violating effects in short-baseline experiments. In this case the tension between appearance and disappearance data is reduced, but cannot be eliminated.

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.

It is shown that neutrino mixing cannot explain the GSI time anomaly, refuting recent claims in this direction. Addendum 1: Remarks on arXiv:0801.1465. Addendum 2: Quantum effects in GSI nuclear decay.

The definition and derivation of flavor neutrino states in the framework of standard Quantum Field Theory is reviewed, clarifying some subtle points. It is shown that a flavor neutrino state that describes a neutrino produced or detected in a charged-current weak interaction process depends on the process under consideration and is appropriate for the description of neutrino oscillations as well as for the calculation of neutrino production or detection rates. Hence, we have a consistent framework for the description of neutrino oscillations and interactions in neutrino oscillation experiments. The standard flavor neutrino states are obtained as approximations which describe neutrinos in experiments that are not sensitive to the dependence of neutrino interactions on the neutrino mass differences. It is also shown that the oscillation probability can be derived either through the usual light-ray time = distance approximation or through an average of the space-time dependent oscillation probability over the unobserved propagation time.

In the first part of this report we present the results of a model independent analysis of the existing solar neutrino data. We obtained forbidden regions in the plane of the parameters Delta m^2 and sin^2 2 theta in two cases: A) Without any restrictions on the values of the solar neutrino fluxes from different reactions; B) With some restrictions that take into account the predictions of all the existing solar models. We show that the existing solar neutrino data allow to exclude rather large regions in the plane of the parameters Delta m^2 and sin^2 2 theta (especially in case B). In the second part of this report we present a general method for the analysis of solar neutrino data that can be applied to future solar neutrino experiments (SNO, Super-Kamiokande, Icarus) in which high energy 8B neutrinos will be detected. We show that these experiments will allow: 1) To reveal in a model independent way the presence of sterile neutrinos in the flux of solar neutrinos on the earth and to obtain lower bounds for the probability of transition of nu_e's into sterile states; 2) To obtain directly from the experimental data the initial 8B nu_e flux and the probability of nu_e's to survive (if there are no transitions of nu_e's into sterile states).

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.

Gravitational lensing introduces the possibility of multiple (macroscopic) paths from an astrophysical neutrino source to a detector. Such a multiplicity of paths can allow for quantum mechanical interference to take place that is qualitatively different to neutrino oscillations in flat space. After an illustrative example clarifying some under-appreciated subtleties of the phase calculation, we derive the form of the quantum mechanical phase for a neutrino mass eigenstate propagating non-radially through a Schwarzschild metric. We subsequently determine the form of the interference pattern seen at a detector. We show that the neutrino signal from a supernova could exhibit the interference effects we discuss were it lensed by an object in a suitable mass range. We finally conclude, however, that -- given current neutrino detector technology -- the probability of such lensing occurring for a (neutrino-detectable) supernova is tiny in the immediate future.

After a short review of the current status of three-neutrino mixing, the implications for the values of neutrino masses are discussed. The bounds on the absolute scale of neutrino masses from Tritium beta-decay and cosmological data are reviewed. Finally, we discuss the implications of three-neutrino mixing for neutrinoless double-beta decay.

It is shown that, in spite of the appearances, the standard expression for the oscillation probability of ultrarelativistic neutrinos is Lorentz invariant.

I review the current status of solar, KamLAND, atmospheric and K2K neutrino experiments, implications for three-neutrino mixing and upper bounds for the neutrino masses.

Compelling evidences in favor of neutrino masses and mixing obtained in the last years in Super-Kamiokande, SNO, KamLAND and other neutrino experiments made the physics of massive and mixed neutrinos a frontier field of research in particle physics and astrophysics. There are many open problems in this new field. In this review we consider the problem of the absolute values of neutrino masses, which apparently is the most difficult one from the experimental point of view. We discuss the present limits and the future prospects of beta-decay neutrino mass measurements and neutrinoless double-beta decay. We consider the important problem of the calculation of nuclear matrix elements of neutrinoless double-beta decay and discuss the possibility to check the results of different model calculations of the nuclear matrix elements through their comparison with the experimental data. We discuss the upper bound of the total mass of neutrinos that was obtained recently from the data of the 2dF Galaxy Redshift Survey and other cosmological data and we discuss future prospects of the cosmological measurements of the total mass of neutrinos. We discuss also the possibility to obtain information on neutrino masses from the observation of the ultra high-energy cosmic rays (beyond the GZK cutoff). Finally, we review the main aspects of the physics of core-collapse supernovae, the limits on the absolute values of neutrino masses from the observation of SN1987A neutrinos and the future prospects of supernova neutrino detection.

It is shown that the 99% confidence limits from the analyses of the data of cosmological and neutrino experiments imply a small marginally allowed region in the space of the neutrino oscillation parameters of 3+1 four-neutrino mixing schemes. This region can be confirmed or falsified by experiments in the near future.

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

Using an analogy with the well-known double-slit experiment, we show that the standard phase of neutrino oscillations is correct, refuting recent claims of a factor of two correction. We also improve the wave packet treatment of neutrino oscillations taking into account explicitly the finite coherence time of the detection process.

We present the results of a Bayesian analysis of solar neutrino data in terms of nu_e->nu_{mu,tau} and nu_e->nu_s oscillations, where nu_s is a sterile neutrino. We perform a Rates Analysis of the rates of solar neutrino experiments, including the first SNO CC result, and spectral data of the CHOOZ experiment, and a Global Analysis that takes into account also the Super-Kamiokande day and night electron energy spectra. We show that the Bayesian analysis does not suffer any problem from the inclusion of the numerous bins of the CHOOZ and Super-Kamiokande energy spectra and allows to reach the same conclusions on the favored type of neutrino transitions and on the determination of the most favored values of the oscillation parameters in both the Rates and Global Analysis. Our Bayesian analysis shows that nu_e->nu_s transitions are strongly disfavored with respect to nu_e->nu_{mu,tau} transitions. In the case of nu_e->nu_{mu,tau} oscillations, the LMA region is favored by the data (86% probability), the LOW region has some small chance (13% probability), the VO region is almost excluded (1% probability) and the SMA region is practically excluded (0.01% probability). We calculate also the marginal posterior probability distributions for tan^2 theta and Delta m^2 and we show that the data imply large mixing almost with certainty and large values of Delta m^2 are favored. We present also the results of a standard least-squares analysis and we show that the standard goodness of fit test is not able to rejects pure nu_e->nu_s transitions. The likelihood ratio test allows to reject pure nu_e->nu_s transitions in favor of nu_e->nu_{mu,tau} transitions only in the Global Analysis.

We present the results of a Bayesian analysis of solar neutrino data in terms of nu_e->nu_{mu,tau} oscillations, independent from the Standard Solar Model predictions for the solar neutrino fluxes. We show that such a model independent analysis allows to constraint the values of the neutrino mixing parameters in limited regions around the usual SMA, LMA, LOW and VO regions. Furthermore, there is a strong indication in favor of large neutrino mixing and large values of Delta m^2 (LMA region). We calculate also the allowed ranges of the neutrino fluxes and we show that they are in good agreement with the Standard Solar Model prediction. In particular, the ratio of the 8B flux with its Standard Solar Model prediction is constrained in the interval [0.45,1.42] with 99.73% probability. Finally, we show that the hypothesis of no neutrino oscillations is strongly disfavored in a model independent way with respect to the hypothesis of neutrino oscillations.

Using the model independent method of Villante, Fiorentini, Lisi, Fogli, Palazzo, and the rates measured in the SNO and Super-Kamiokande solar neutrino experiment, we calculate the amount of active nu_mu or nu_tau present in the flux of solar neutrinos on Earth. We show that the probability of nu_e->nu_{mu,tau} transitions is larger than zero at 99.89% CL. We find that the averaged flux of nu_{mu,tau} on Earth is larger than 0.17 times the 8B nu_e flux predicted by the BP2000 Standard Solar Model at 99% CL. We discuss also the consequences of possible nu_e->anti-nu_{mu,tau} or nu_e->anti-nu_e transitions of solar neutrinos. We derive a model-independent lower limit of 0.52 at 99% CL for the ratio of the 8B nu_e flux produced in the Sun and its value in the BP2000 Standard Solar Model.

We calculate with Monte Carlo the goodness of fit and the confidence level of the standard allowed regions for the neutrino oscillation parameters obtained from the fit of the total rates measured in solar neutrino experiments. We show that they are significantly overestimated in the standard method. We also calculate exact allowed regions with correct frequentist coverage. We show that the exact VO, LMA and LOW regions are much larger than the standard ones and merge together giving an allowed band at large mixing angles for all Delta m^2 > 10^{-10} eV^2.

We present a simple but general treatment of neutrino oscillations in the framework of quantum mechanics using plane waves and intuitive wave packet principles when necessary. We attempt to clarify some confusing statements that have recently appeared in the literature.

We present an analysis of the neutrino oscillation solutions of the solar neutrino problem in the framework of four-neutrino mixing where a sterile neutrino is added to the three standard ones. We perform a fit to the full data set corresponding to the 825-day Super-Kamiokande data sample as well as to Chlorine, GALLEX and SAGE and Kamiokande experiments. In our analysis we use all measured total event rates as well as all Super-Kamiokande data on the zenith angle dependence and the recoil electron energy spectrum. We consider both transitions via the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism as well as oscillations in vacuum (just-so) and find the allowed solutions for different values of the additional mixing angles. This framework permits transitions into active or sterile neutrinos controlled by the additional parameter $\cos^2(\vartheta_{23}) \cos^2(\vartheta_{24})$ . We discuss the maximum allowed values of this additional mixing parameter for the different solutions.

This review is focused on neutrino mixing and neutrino oscillations in the light of the recent experimental developments. After discussing possible types of neutrino mixing for Dirac and Majorana neutrinos and considering in detail the phenomenology of neutrino oscillations in vacuum and matter, we review all existing evidence and indications in favour of neutrino oscillations that have been obtained in the atmospheric, solar and LSND experiments. We present the results of the analyses of the neutrino oscillation data in the framework of mixing of three and four massive neutrinos and investigate possibilities to test the different neutrino mass and mixing schemes obtained in this way. We also discuss briefly future neutrino oscillation experiments.

It is shown that the future SNO and Super-Kamiokande experiments, in which high energy $^8\mathrm{B}$ neutrinos will be detected through the observation of CC, NC and $\nu$--$e$ elastic scattering processes, could allow to reveal in a model independent way the presence of sterile neutrinos in the flux of solar neutrinos on the earth. Lower bounds for different averaged values of the probability of transition of solar $\nu_e$'s into sterile states and for the total flux of $^8\mathrm{B}$ neutrinos are derived in terms of measurable quantities. The possibilities to reveal the presence of $\nu_\mu$ and/or $\nu_\tau$ in the solar neutrino flux on the earth are also considered and the case of transitions of solar $\nu_e$'s only into sterile states is discussed. Some numerical results for a simple model with $\nu_e$--$\nu_s$ mixing are given.

Motivated by the indication that both the atmospheric and the solar neutrino puzzles may simultaneously be solved by (vacuum as well as matter-induced resonant) oscillations of two generations of neutrinos with large mixing, we have analyzed the data on the atmospheric and solar neutrinos assuming that all {\it three} neutrinos are maximally mixed. It is shown that the values of $ \Delta m^2 $ obtained from the two-generation analyses are still valid even in the three-generation scheme, i.e. the two puzzles can be solved simultaneously if $ \Delta m_{31}^2 \simeq 10^{-2} \, \mathrm{eV}^2 $ for the atmospheric neutrinos and $ \Delta m_{21}^2 \simeq 10^{-10} \, \mathrm{eV}^2 $ for solar neutrinos in the maximally mixed three-generation scheme.

The results of a solar model independent analysis of the existing solar neutrino data under the assumption of MSW transitions in matter or vacuum oscillations (in the cases of $\nu_{e}$--$\nu_{\mu}(\nu_{\tau})$ and $\nu_{e}$--$\nu_{\mathrm{S}}$ mixing) are presented. The analysis was done in two cases. In the first case no assumptions on the value of the total fluxes of neutrinos from all reactions has been made. In the second case some boundaries for the values of the neutrino fluxes, that take into account the predictions of different standard solar models, are imposed. It is shown that rather large regions of the parameters $ \Delta m^2 $ and $ \sin^2 2\theta $ are excluded by the existing data. It is also shown that the astrophysical solution of the solar neutrino problem can be excluded only in the case in which limitations on the values of the solar neutrino fluxes are imposed.

We discuss here what model independent information about properties of neutrinos and of the sun can be obtained from future solar neutrino experiments (SNO, Super-Kamiokande). It is shown that in the general case of transitions of solar $\nu_e$'s into $\nu_\mu$ and/or $\nu_\tau$ the initial $^8\mathrm{B}$ neutrino flux can be measured by the observation of NC events. From the CC measurements the $\nu_e$ survival probability can be determined as a function of neutrino energy. The general case of transitions of solar $\nu_e$'s into active as well as sterile neutrinos is considered. A number of relations between measurable quantities the test of which will allow to answer the question whether there are sterile neutrinos in the solar neutrino flux on the earth are derived. Transitions of solar $\nu_e$'s into active and sterile states due to neutrino mixing and Dirac magnetic moments or into active left-handed neutrinos and active right-handed antineutrinos due to neutrino mixing and Majorana transition magnetic moments are also considered. It is shown that future solar neutrino experiments will allow to distinguish between the cases of Dirac and Majorana magnetic moments.

The main features of the observables in the SNO and Super-Kamiokande solar neutrino experiments in the case of neutrino spin and/or spin-flavour precession in the magnetic field of the sun are discussed. It is shown without any model dependent assumption that in the case of Majorana transition magnetic moments the NC event rate $N^{\mathrm{NC}}$ does not depend on time and a measurement of $N^{\mathrm{NC}}$ will allow to determine the initial flux of $^8\mathrm{B}$ neutrinos with a theoretical uncertainty of a few \%. In the case of Dirac magnetic moments the NC event rate will depend on time and we obtain a model independent lower bound for the transition probability of initial $\nu_e$'s into right-handed sterile neutrinos.