Neutrino Oscillations

Abstract

In the last decades, a very important breakthrough has been brought about in the elementary particle physics by the discovery of the phenomenon of the neutrino oscillations, which has shown neutrino properties beyond the Standard Model. But a full understanding of the various aspects of the neutrino oscillations is far to be achieved. In this paper the theoretical background of the neutrino oscillation phenomenon is described, referring in particular to the paradigmatic models. Then the various techniques and detectors which studied neutrinos from different sources are discussed, starting from the pioneering ones up to the detectors still in operation and to those in preparation. The physics results are finally presented adopting the same research path which has been crossed by this long saga. The problems not yet fixed in this field are discussed, together with the perspectives of their solutions in the near future.

(ProQuest: ... denotes non-US-ASCII text omitted.)

G. Bellini 1 and L. Ludhova 1 and G. Ranucci 1 and F. L. Villante 2

Academic Editor:Elisa Bernardini

1, Dipartimento di Fisica, Università degli Studi and INFN, 20133 Milano, Italy
2, Dipartimento di Scienze Fisiche e Chimiche, Università degli Studi and INFN, 67100 L' Aquila, Italy

Received 13 August 2013; Revised 15 October 2013; Accepted 29 October 2013; 20 January 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP³ .

1. Introduction

Neutrino studies brought us to some of the most relevant breakthroughs in particle physics of last decades. In spite of that, the neutrino properties are still far to be completely understood.

The discovery of the oscillation phenomenon produced quite a revolution in the Standard Model of elementary particles, especially through the direct evidence of a nonzero neutrino mass. The first idea of neutrino oscillations was considered by Pontecorvo in 1957 [1-3], before any experimental indication of this phenomenon. After several-decades-lasting saga of experimental and theoretical research, many questions are still open around the interpretation of this phenomenon and on the correlated aspects, on the oscillation parameters, on the neutrino masses, on the mass hierarchy, on CP violation in the leptonic sector, and on a possible existence of a fourth, sterile neutrino.

The generally accepted MSW model [4-6] to interpret solar neutrino oscillations is presently validated for the oscillation in vacuum and in matter, but not yet in the vacuum-matter transition region. The shape of this transition could be influenced in a relevant way, as suggested by various theories going beyond the Standard Model as, for example, the nonstandard neutrino interactions and a possible existence of a very light sterile neutrino. For this reason, the transition region deserves further and refined experimental studies.

Checks on the neutrino oscillations are under way through several experiments in data-taking phase, while few others are in preparation or even construction. These projects exploit various approaches, for example, neutrino-flavor disappearance and appearance, short and long source-to-detector baselines, and measure neutrinos and/or antineutrinos of various origins, as the solar, atmospheric, accelerator, geo-, and reactor (anti)neutrinos.

Neutrinos interact with matter only through weak interactions and, thus, they can bring to the observer almost undistorted information about their source. For example, by studying solar neutrinos and geoneutrinos, we gather information not only about the character of neutrino itself but also about the Sun's and the Earth's interior.

This paper consists of five sections. In Section 2, the theoretical aspects of neutrino oscillations are reviewed (for a more detailed discussion, see, e.g., [7-9]). The principles and the structures of the detectors employed in neutrino-physics experiments are discussed in Section 3. The most important milestones and the results of neutrino-physics experiments are summarized in Section 4. We briefly discuss the opened problems of neutrino physics and what can be expected in the near future of this exciting research field in Section 5.

2. Neutrino Oscillations

In the Standard Model (SM) neutrinos are neutral, massless fermions. They only interact with other particles via weak interactions, which are described by the charged-current (CC) and neutral current (NC) interaction Lagrangians: [figure omitted; refer to PDF] In the above relation, g is the S U ( 2 _{) L} gauge coupling constant, _{θ W} is the weak angle, and the charged and neutral currents ^{j ρ CC} and ^{j ρ NC} are given by [figure omitted; refer to PDF] where [cursive l] are the charged lepton fields and we have written only the terms containing the neutrino fields _{ν [cursive l]} .

If neutrinos have non-zero masses, the left handed components _{ν α L} of the neutrino fields with definite flavor α (which enter in the CC current definition) can be a superposition of the left handed components _{ν i L} of the neutrino fields with definite masses _{m i} (in this section, we use Greek letters α and β to refer to neutrino flavors and Latin letters i and j to refer to neutrino masses). Assuming that neutrinos are ultrarelativistic, we have [figure omitted; refer to PDF] where U is an unitary matrix. By considering that a field operator creates antiparticles, this implies that a flavor eigenstate | _{ν α} YA; is a superposition of different mass eigenstates | _{ν i} YA; , according to [figure omitted; refer to PDF] For antineutrinos, we obtain correspondingly [figure omitted; refer to PDF] In principle, the number N of massive neutrinos can be larger than three. In this case, however, we must assume that there are sterile neutrinos, that is, light fermions that do not take part in standard weak interactions (1) and (2) and thus are not excluded by LEP results according to which the number of active neutrinos coupled with the ^{W ±} and Z boson is _{N ν} = 2.984 ± 0.008 [10].

In the assumption of 3 massive neutrinos, the neutrino mixing matrix U can be expressed in terms of three mixing angles _{θ 12} , _{θ 23} , and _{θ 13} and one Dirac-type CP phase δ according to [figure omitted; refer to PDF] where _{R i j} ( _{θ i j} ) represents an Euler rotation by _{θ i j} in the i j plane and [figure omitted; refer to PDF] We are considering here the assumption that neutrinos are Dirac particles. In the case of Majorana (or Dirac-Majorana) mass terms, the most general form of the mixing matrix contains two additional phases and it is obtained by U [arrow right] U · _{U M} , where _{U M} = diag ... ( 1 , ^{e i _{[varphi] 1}} , ^{e i _{[varphi] 2}} ) ; see for example, [7, 8]. The Majorana phases _{[varphi] 1} and _{[varphi] 2} , however, have no observable effects on neutrino oscillations [11]. In components, the mixing matrix U is expressed as [figure omitted; refer to PDF]

where _{c i j} = cos ... _{θ i j} and _{s i j} = sin _{θ i j} . We indicate with Δ ^{m i j 2} ...1; ^{m i 2} - ^{m j 2} . As it is usually done, we order the neutrino masses such that Δ ^{m 21 2} > 0 and Δ ^{m 21 2} ...a; | Δ ^{m 31 2} | . With this choice, the ranges of mixing parameters are determined by [figure omitted; refer to PDF] The sign of Δ ^{m 31 2} determines the neutrino mass hierarchy, being Δ ^{m 31 2} > 0 for normal hierarchy (NH) and Δ ^{m 31 2} < 0 for inverted hierarchy (IH).

2.1. Neutrino Evolution Equation

The evolution of a generic neutrino state | ν ( t ) YA; is described by a Schrödinger-like equation: [figure omitted; refer to PDF] where H represents the Hamiltonian operator. The above equation can be expressed in the flavor eigenstate basis { | _{ν α} YA; } . We obtain [figure omitted; refer to PDF] where ^{ν ( f )} ( t ) is the vector describing the flavor content of the neutrino state | ν ( t ) YA; given by [figure omitted; refer to PDF] with _{a α} ( t ) = Y9; _{ν α} |" ν ( t ) YA; , and the matrix _{H f} is given by: [figure omitted; refer to PDF]

In vacuum, the neutrino Hamiltonian _{H vac} is determined in terms of neutrino masses and mixing parameters. We have in fact [figure omitted; refer to PDF] where ^{H vac ( m )} is the representation of the vacuum Hamiltonian in the mass eigenstate basis, given by [figure omitted; refer to PDF] In the last equality, we adopted the ultra-relativistic approximation E [approximate] | p [arrow right] | + ^{m 2} / 2 | p [arrow right] | and we implicitly assumed that the neutrino state | ν ( t ) YA; can be described as a superposition of states with fixed momentum p [arrow right] . This corresponds to the so-called plane-wave approximation which is adequate to describe neutrino evolution when coherence of the different components of the neutrino wave packet is not lost in the detection and/or propagation processes (for a wave-packet description of neutrino oscillations see, e.g., [7]).

The presence of a matter can affect neutrino propagation in a nontrivial way. In fact, as it was first realized by [4], when a neutrino propagates through a medium, its dispersion relation (i.e., its energy-momentum relation) is modified by coherent interactions with background particles. This phenomenon, which in optics is accounted for by introducing a refractive index, can be described by adding an effective potential V in the evolution equation, so that [figure omitted; refer to PDF] In the SM, the effective potential is diagonal in the flavor basis. We thus have [figure omitted; refer to PDF] where we have taken into account that sterile states do not interact with the medium. At low energies, the potentials can be evaluated by taking the average Y9; _{[Hamiltonian (script capital H)] eff} YA; of the effective four-fermion Hamiltonian due to exchange of W and Z bosons over the state describing the background medium. We have [figure omitted; refer to PDF] with [figure omitted; refer to PDF] where _{G F} is the Fermi constant, ^{g V b} and ^{g A b} are the vector and axial vector coupling constants of the various background particles and we have performed a Fierz reshuffling of the fields; see [7, 12]. In the above equations, it is taken into account that normal matter does not contain muons and taus and, consequently, the CC interactions with the medium only affect electron neutrino propagation. For nonrelativistic unpolarized medium, one obtains [figure omitted; refer to PDF] where the CC contribution [figure omitted; refer to PDF] is proportional to difference between the number densities of electrons and positrons. The NC contribution _{A NC} is equal for all active neutrino flavors and it is given by [figure omitted; refer to PDF] where _{s W} ...1; sin _{θ W} while _{n p} and _{n n} ( _{n p ¯} and _{n n ¯} ) are the number densities of protons and neutrons (antiprotons and antineutrons), respectively.

In neutral matter, it necessarily holds ( _{n e} - _{n e ¯} ) = ( _{n p} - _{n p ¯} ) that implies that the first term in the r.h.s. of the above equation vanishes. Moreover, in the absence of sterile neutrinos, the neutral current contribution to the total Hamiltonian is proportional to the identity matrix. As a consequence, it only introduces an overall unobservable phase factor in the evolution of ^{ν ( f )} ( t ) and, thus, can be neglected.

Finally, the evolution equation for antineutrinos is obtained by replacing U [arrow right] ^{U *} in (16) and V [arrow right] - V in (18). We, thus, understand that CP-violating effects are absent in neutrino oscillations, if the mixing matrix is real (i.e., U = ^{U *} ) and neutrinos propagate in vacuum or in a CP-symmetric medium (i.e., V = 0 ).

2.2. Oscillations in Vacuum and in Matter

In formal terms, neutrino oscillations are easily described. Let us assume that a neutrino flavor _{ν α} is created at a time _{t 0} = 0 . In the flavor eigenstate basis, this state is represented by a vector [figure omitted; refer to PDF] with components _{a β} ( 0 ) = _{δ β α} . After a time interval t , the neutrino propagated to a distance x [approximate] t and its flavor content has evolved according to [figure omitted; refer to PDF] where the evolution operator is given by [figure omitted; refer to PDF] and T represents the time ordering of the exponential. In the presence of neutrino mixing and if neutrino masses are not degenerate, the Hamiltonian ^{H ( f )} is not diagonal. Thus, flavor is not conserved and components β ...0; α can appear as a result of the evolution. The probability to detect a neutrino flavor _{ν β} at a distance L from the neutrino production point is given by [figure omitted; refer to PDF] In the following, we discuss the expectations for P ( _{ν α} [arrow right] _{ν β} ) in few relevant cases.

2.3. Vacuum Neutrino Oscillations

In vacuum, the neutrino Hamiltonian H is constant. The evolution operator can be explicitly calculated as [figure omitted; refer to PDF] where ^{S ( m )} is the evolution operator in the mass eigenstate basis given by [figure omitted; refer to PDF] with _{[varphi] i} = - ^{m i 2} x / 2 | p [arrow right] | . The probability to observe the oscillation _{ν α} [arrow right] _{ν β} over a distance L is thus given by [figure omitted; refer to PDF] where _{[varphi] i j} = [ ( ^{m j 2} - ^{m i 2} ) L ] / 2 E and we considered that for a relativistic particle E [approximate] | p [arrow right] | .

The above expression can be recast in few alternative forms that are useful to discuss the property of neutrino oscillations. We obtain, for example, [figure omitted; refer to PDF] which gives the oscillation probability as the sum of a constant and an oscillating term. The oscillating part averages to zero if the phases _{[varphi] i j} vary over ranges Δ _{[varphi] i j} ...b; 1 , as it can be due, for example, to a spread of the neutrino energy E and/or the neutrino baseline L . The constant part represents the "classical" limit that is obtained by neglecting interference among the different components of the neutrino wave packet and by combining probabilities, rather than amplitudes, to derive P ( _{ν α} [arrow right] _{ν β} ) .

Alternatively, we can write [figure omitted; refer to PDF] The first two terms in the r.h.s. of the above equation do not change for U [arrow right] ^{U *} and describe the CP-conserving part of the neutrino oscillation probability. The last part, instead, changes sign introducing a difference between neutrino and antineutrino oscillation probabilities that can be quantified as [figure omitted; refer to PDF] For α = β , this term vanishes showing that CP asymmetry can be measured only in transitions between different neutrino flavors.

If we assume two neutrino mixing, that is, we take only one nonvanishing mixing angle _{θ i j} in (8), the oscillation probability reduces to the well known expression [figure omitted; refer to PDF] where α ...0; β and the involved flavors depend on the mixing angle _{θ i j} (an angle _{θ 12} ...0; 0 induces _{ν e} [arrow right] _{ν μ} oscillations; _{θ 13} ...0; 0 induces _{ν e} [arrow right] _{ν τ} oscillations; _{θ 23} ...0; 0 induces _{ν μ} [arrow right] _{ν τ} oscillations). The survival probability for the case α = β can be simply deduced by considering that, due to unitarity of the mixing matrix, it always holds _{∑ β} ... P ( _{ν α} [arrow right] _{ν β} ) ...1; 1 that, in this specific case, gives [figure omitted; refer to PDF] Equation (36) describes an oscillating function of L . The amplitude of the oscillation is determined by ^{sin 2} ( 2 _{θ i j} ) while the oscillation length is given by [figure omitted; refer to PDF] The oscillation probabilities are unchanged when Δ ^{m i j 2} [arrow right] - Δ ^{m i j 2} or _{θ i j} [arrow right] π / 2 - _{θ i j} showing that two neutrino oscillations in vacuum do not probe the hierarchy of the masses _{m i} and _{m j} (i.e., the states _{ν i} and _{ν j} can be interchanged with no effect on (36) and (37)).

In the three-neutrino case, useful expressions can be derived in the approximation of one-dominant mass scale (i.e., Δ ^{m 21 2} ...a; | Δ ^{m 31 2} | [approximate] | Δ ^{m 32 2} | ) which is motivated by the fact that the mass difference required to explain the atmospheric neutrino anomaly is much larger than that required to solve the solar neutrino problem. In this assumption (note that, due to CPT-invariance, it holds P ( _{ν β} [arrow right] _{ν α} ) = P ( _{ν ¯ α} [arrow right] _{ν ¯ β} ) ) one obtains [figure omitted; refer to PDF] where, following [9], we adopted the notation _{S 23} = ^{sin 2} ( Δ ^{m 32 2} L / 4 E ) and _{S 12} = ^{sin 2} ( Δ ^{m 21 2} L / 4 E ) and we set _{θ 13} = 0 in the coefficients of the _{S 12} terms. The CP-violating part _{P CP} , which enters with opposite sign in the corresponding expressions for antineutrinos, is given by [figure omitted; refer to PDF] where [figure omitted; refer to PDF] showing that CP violation is observed in neutrino oscillations only if all the angles and all the mass differences are nonvanishing. The magnitude of CP violating effects depends on the phase δ , being maximal for δ = π / 2 and δ = 3 π / 2 .

The survival probabilities P ( _{ν α} [arrow right] _{ν α} ) = P ( _{ν ¯ α} [arrow right] _{ν ¯ α} ) are given by [9] [figure omitted; refer to PDF]

We note that, at this level of approximation, there is no sensitivity to neutrino hierarchy since the oscillation probabilities do not depend on the sign of Δ ^{m 31 2} . Moreover, in the limit _{θ 13} [arrow right] 0 , the "atmospheric" mass scale Δ ^{m 32 2} does not produce observable effects on electron neutrino oscillations that can be regarded as two neutrino oscillations, driven by the "solar" mass difference Δ ^{m 21 2} , between _{ν e} and the mixed state _{ν μ τ} = _{c 23 ν μ} - _{s 23 ν τ} . This conclusion also holds in presence of matter.

2.4. Neutrino Oscillations in Matter

The evolution of neutrinos in matter is complicated by the fact that the properties of the medium can change along the neutrino trajectory, thus giving a nonconstant Hamiltonian. The evolution equation reads as [figure omitted; refer to PDF] where, if we neglect sterile neutrinos, the only nonvanishing entry of the matrix ^{V ( f )} ( x ) is [figure omitted; refer to PDF] Here, the "+" sign refers to neutrinos while the "-" sign refers to antineutrinos. In the above equations, we omitted the NC contribution to matter potential that is proportional to the identity matrix. We also assumed that the number density of positrons is negligible. See Section 2.1 for details.

It is convenient to diagonalize the Hamiltonian at each point of the space and discuss the evolution in the basis of the local mass eigenstates defined by the following relation: [figure omitted; refer to PDF] where U ~ ( x ) is the unitary matrix that gives [figure omitted; refer to PDF] with [figure omitted; refer to PDF] In this basis, the evolution equation becomes [figure omitted; refer to PDF] We see that the nondiagonal entries, which may cause the transitions between the local mass eigenstates, are proportional to the derivative of U ~ ( x ) whose magnitude is essentially determined by the rate of change of the electrons number density in the background medium.

This observation can be used to introduce the so-called adiabatic approximation that applies with good accuracy to the case of solar neutrino oscillations. Let us indicate with _{L ~ i j} [approximate] 4 π E / | Δ ^{m ~ i j 2} | the length scale over which the components of the neutrino wave packet with masses _{m ~ i} and _{m ~ j} acquire a phase difference Δ _{Φ i j} = 2 π . If we assume that the various _{L ~ i j} are much smaller than the distance over which the medium changes its properties D ...1; ( d ln ... _{n e} ( x ) / d x ^{) - 1} , the second term in the r.h.s. of (52) can be neglected. Thus, the components of the vector ^{ν ( m ~ )} ( x ) remain constant (in magnitude) during the evolution, even if the decomposition of ^{ν ( m ~ )} ( x ) in the flavor basis changes along the neutrino trajectory as a result of the variations of _{n e} . If the length scales _{L ~ i j} are also much smaller than the baseline L over which neutrinos propagate, the oscillation probabilities P ( _{ν α} [arrow right] _{ν β} ) only depends on the properties of U ~ ( x ) at the production point _{x p} and at the detection point _{x d} . They can be, in fact, deduced by incoherently combining probabilities of production and detection, obtaining [figure omitted; refer to PDF]

We now consider the specific case of _{ν e} produced by nuclear reactions occurring at the center of the sun. Let us calculate the electron neutrino survival probability, by first considering a two neutrino scenario in which only _{θ 12} ...0; 0 . The effective mixing angle in matter _{θ ~ 12} can be calculated as [figure omitted; refer to PDF] while the difference between the effective neutrino masses is given by [figure omitted; refer to PDF] with [figure omitted; refer to PDF] Matter effects break the degeneracies Δ ^{m 21 2} [arrow right] - Δ ^{m 21 2} and _{θ 12} [arrow right] π / 2 - _{θ 12} probing the hierarchy in the 1-2 neutrino sector. In particular, when Δ ^{m 21 2} > 0 and _{θ 12} < π / 4 , the system has a resonance . There exists, in fact, a value of the electron number density, defined by the following condition: [figure omitted; refer to PDF] for which the local mixing is maximal (i.e., _{θ ~ 12} = π / 4 ) while the mass difference Δ ^{m ~ 12 2} reaches the minimal value Δ ^{m ~ 12 2} = Δ ^{m 12 2} sin ( 2 _{θ 12} ) . As it was discussed by [5, 6], if the resonance region is sufficiently wide, it is possible to achieve a total conversion of _{ν e} into neutrinos of different flavors. This mechanism is called the MSW effect . Considering that the electron density in the Sun is _{n e} [<, ~] 1 ^{0 26} cm^-3 and the typical solar neutrino energies are E [approximate] 1 MeV, the resonance condition requires Δ ^{m 21 2} cos ... ( 2 _{θ 12} ) [<, ~] 1 ^{0 - 5} eV² .

The evolution equation in the local mass eigenstate basis becomes [figure omitted; refer to PDF] and the adiabaticity condition can be explicitly expressed as [figure omitted; refer to PDF] where the adiabaticity parameter γ is given by the ratio between the differences of diagonal elements and off-diagonal elements of (58): [figure omitted; refer to PDF] If condition (59) is fulfilled, the electron neutrino survival probability can be calculated through (53) obtaining [figure omitted; refer to PDF] where _{θ ~ 12} indicates the mixing angle at neutrino production point and we assumed that neutrinos are detected in vacuum.

In order to understand the specific features of P ( _{ν e} [arrow right] _{ν e} ) , it is useful to define a transition energy ^{E *} , given by: [figure omitted; refer to PDF] where _{n e , [ecedil]9;} is the electron number density at the center of the sun. For E ...a; ^{E *} , matter effects are negligible and (61) reduces to: [figure omitted; refer to PDF] which, in fact, corresponds to vacuum averaged neutrino oscillations. For E ...b; ^{E *} , matter potential becomes dominant so that the "heaviest" mass eigenstates in the center of the Sun coincide with _{ν e} . As a consequence, we obtain cos ... ( 2 _{θ ~ 12} ) = - 1 and [figure omitted; refer to PDF] For the value of _{θ 12} and Δ ^{m 21 2} currently favored by neutrino oscillation analysis (see Section 4), the transition energy ^{E *} is approximately ^{E *} [approximate] 1.2 MeV.

The violations of adiabaticity can be taken into account by introducing the crossing probability _{P C} that represents the probability of a transition between the local mass eigenstates during the neutrino evolution. If _{P c} ...0; 0 , the electron neutrino survival probability becomes [figure omitted; refer to PDF] There are different approaches to calculate _{P C} . For several cases of interest, the following expression holds (see e.g., [9, 13, 14] and references therein): [figure omitted; refer to PDF] where γ ~ is the minimal value of γ ( x ) along the neutrino trajectory (in the presence of a resonance, one can often approximate γ ~ [approximate] _{γ res} where _{γ res} is the value of γ ( x ) at the resonance point, see, e.g., [7, 9]) and the parameter F depends on the adopted electron density profile. In particular, for an exponential density profile _{n e} [proportional, variant] exp ... ( - x ) , which is a good approximation for solar neutrinos, one has F = 1 - ^{tan 2} _{θ 12} .

In the case of three mixed neutrinos, the above picture has to be modified to take into account the possibility that _{θ 23} ...0; 0 and _{θ 13} ...0; 0 . Since matter potentials are equal for muon and tau neutrinos, the rotation R ( _{θ 23} ) in (8) can be reabsorbed in the "mixed" basis { | _{ν e} YA; , _{c 23} | _{ν μ} YA; - _{s 23} | _{ν τ} YA; , _{s 23} | _{ν μ} YA; + _{c 23} | _{ν τ} YA; } . This shows that, when _{θ 13} = 0 , electron neutrinos experience two-neutrino oscillations to a mixed state | _{ν μ τ} YA; = _{c 23} | _{ν μ} YA; - _{s 23} | _{ν τ} YA; and, thus, the electron neutrino survival probability is unchanged. In the presence of _{θ 13} ...0; 0 , we have instead non trivial modifications due to the fact that the state | _{ν e} YA; mixes with the state | _{ν 3} YA; being in fact | Y9; _{ν e} |" _{ν 3} YA; | = _{s 13} . By repeating the previous calculations, one obtains [figure omitted; refer to PDF] where it is assumed that matter effects negligibly modify the _{θ 13} mixing angle (i.e., _{θ 13} [approximate] _{θ ~ 13} ).

We finally remark that the above expression applies to solar neutrinos detected during the day, since these neutrinos do not cross the Earth to reach the detector. Matter effects due to propagation across the Earth can modify (67) by introducing a day-night modulation whose magnitude depends on the specific values of mass and mixing parameters.

3. Neutrino Detectors

The successful series of solar, atmospheric, reactor, and accelerator experiments which led to firmly establish the standard three-flavor neutrino oscillation paradigm involved the realization of sophisticated detectors based on a plurality of techniques. In this paragraph we briefly review their main features, which undoubtedly played a key role in the incredible success of this field.

3.1. Radiochemical Detectors

The emerging hint of the so-called solar neutrino problem at the beginning of the 70s from the first results of the pioneering chlorine experiment (whose final findings are summarized in [24]), carried out by Ray Davies in the Homestake mine, signaled the experimental beginning of the neutrino oscillation saga. The problem, consisting in a sizable discrepancy between the data and the prediction of the Standard Solar Model, persisted for more than 30 years before being explained as a manifestation of the neutrino oscillation phenomenon. A beautiful account of the early stage of this field can be found in the seminal book of Bahcall [25], where all the steps which brought to shape unambiguously the existence of the experimental puzzle are vividly and clearly explained. In the 90s additional evidence of the existence of the solar neutrino problem came from other two radiochemical experiments, GALLEX/GNO [26] at Gran Sasso and SAGE [27, 28] at Baksan. The principle of the radiochemical technique is very simple and elegant: the detection medium is a material which, upon absorption of a neutrino, is converted into a radioactive element whose decay is afterwards revealed and counted. The Homestake experiment used a chlorine solution as a target for inverse β -interactions: [figure omitted; refer to PDF] characterized by a threshold of 0.814 MeV. It is worth reminding that such a technique was proposed independently by two giants of modern physics, Bruno Pontecorvo and Louis Alvarez. The other two experiments, instead, adopted gallium as target, which allows neutrino interaction via [figure omitted; refer to PDF] The threshold of this reaction is 233 keV, low enough to essentially probe the entire solar neutrino spectrum (see Section 4.1 for details) which on the contrary cannot be revealed with the chlorine reaction due to the higher threshold. Due to the similarity of the methodology in both cases of chlorine and gallium, in the following its description is focused to the specific gallium implementation. In GALLEX/GNO the target consisted of 101 tons of a GaCl₃ solution in water and HCl, containing 30.3 tons of natural gallium; this amount corresponds to about 1 ^{0 29} Ga 71 nuclei. The solution was contained in a large tank hosted in Hall A of the underground Gran Sasso Laboratory.

Ge 71 produced by neutrinos is radioactive and decays back by electron capture into Ga 71 . The mean life of a Ge 71 nucleus is about 16 days; thus Ge 71 accumulates in the solution, asymptotically reaching equilibrium when the number of Ge 71 atoms produced by neutrino interactions is just the same as the number of the decaying ones. In this equilibrium condition, about a dozen Ge 71 atoms would be present inside the whole gallium chloride solution. Since the exposure time is in practice limited to four weeks, the actual number of Ge 71 atoms is less than the equilibrium value, but still perfectly predictable. Therefore, the solar neutrino flux above threshold is deduced from the number of Ge 71 produced atoms, using the theoretically calculated cross section. The challenging experimental task is thus to identify the feeble amount of Ge 71 atoms. This is accomplished through a complex procedure which contemplates several steps.

(1) The solution is exposed to solar neutrinos for about four weeks.

(2) The Ge 71 atoms present at the end of the four week period in the solution are in the form of volatile GeCl₄ , which is extracted into water by pumping about 3000 m³ of nitrogen through the solution.

(3) The extracted Ge 71 is converted into gaseous GeH₄ and introduced into miniaturized proportional counters mixed with xenon as counting gas. At the end of the process, a quantity variable between 95 and 98% of the Ge 71 present in the solution at the time of extraction is in the counter; extraction and conversion efficiencies are under constant control using nonradioactive germanium isotopes as carriers.

(4) Decays and interactions in the counter are observed for a period of 6 months, allowing the complete decay of Ge 71 and a good determination of the counter background. The charge pulses produced in the counters by decays are recorded by means of fast transient digitizers.

(5) The data, after application of suitable cuts, are then analyzed with a maximum likelihood algorithm to obtain the most probable number of Ge 71 introduced in the counter, with some final corrections applied to take into account the so-called side reaction, that is interactions in the solution generated by high energy muons from cosmic rays and by natural radioactivity.

The key issue in the overall procedure is the minimization of the possible sources of backgrounds. This is performed through a triple strategy, whose first element is the rigorous application of low-level radioactivity technology in the design and construction of the counters; the second element is the use in the analysis of sophisticated pattern recognition techniques able to perform energy and shape discrimination of the signal and background events; the third and final element is the precise calibration of the counters via an external Gd/Ce X-ray source, to enhance the accuracy of the signal/background discrimination ensured by the pattern recognition method.

Thanks to the effective methodology adopted, the radiochemical experiments were able to provide very important results in the studies of solar neutrinos, demonstrating unambiguously the discrepancy between the measured and predicted solar neutrino fluxes and triggering the subsequent vast theoretical and experimental investigations culminated in the proof of the oscillation effect in the solar neutrino sector.

Such fundamental outputs were achieved despite the incredible challenge of the measurement, which can be well appreciated by considering the smallness of the detected signal. In about two decades of operation, Homestake and SAGE detected 860 and 870 decays, respectively, as reported in [24, 27, 28] (GALLEX/GNO did not publish this number).

In this respect, is worth mentioning another important ingredient of the radiochemical solar neutrino programs, which is the source calibration efforts which were performed to prove unambiguously the validity of the entire neutrino detection concept implied by this technique. In particular, GALLEX and SAGE underwent twice the calibration procedure. GALLEX exploited in both cases a Cr 51 source [34], while SAGE adopted two different isotopes, Cr 51 in the first instance [35] and Ar 37 in the second test [36]. The outcome of the source tests was the definitive validation of the radiochemical approach as an effective method to detect neutrinos. However, the ratio R between the detected and predicted neutrino fluxes is significantly less than 1; taking the four tests together, the global result is R = 0.86 ± 0.05 . This anomaly can be interpreted as a possible indication of _{ν e} disappearance (see, e.g., [37]) within models with additional sterile neutrino states (see Section 4.2.4).

3.2. Cerenkov Detectors

The widespread diffusion of the Cerenkov technique in the field of neutrino physics can be appreciated by considering the many experimental setups based on this method which have been employed to investigate the entire neutrino spectrum, from the lowest to the highest energies.

The Cerenkov radiation is produced in a material with refractive index n by a charged particle if its velocity is greater than the local phase velocity of the light. The charged particle polarizes the atoms along its trajectory, generating time dependent dipoles which in turn generate electromagnetic radiation. If v < c / n , the dipole distribution is symmetric around the particle position, and the sum of all dipoles vanishes. If v > c / n , the distribution is asymmetric and the total time dependent dipole is different from zero and thus radiates.

The resulting light wavefront is conical, characterized by an opening angle whose cosine is equal to 1 / ( β n ) ; the spectrum of the radiation is ultraviolet divergent, being proportional to 1 / ^{λ 2} . The propagation properties of the Cerenkov light are therefore fully equivalent to those of the acoustic Mach cone.

3.2.1. SNO

The SNO experiment [38] is a paradigmatic example of how the Cerenkov light can be used as basis to build a very effective neutrino detector. Since SNO encompasses more experimental features than the other important detector of this kind, the Japanese Super-Kamiokande described in Section 3.2.2, we find convenient, for illustrative purposes, to reverse the historical order (the data taking of Super-Kamiokande actually started before SNO). Located underground, in the Inco mine at Sudbury (Canada), this detector employed heavy water, which acted both as target medium for the neutrinos and as light generating material. The basic idea beyond the choice of heavy water is to perform two independent solar neutrino measurements based on the deuterium target: the first is aimed to detect specifically the electron neutrino component, while the second is sensitive to the all-flavor flux. Thus, the comparison of the two results can permit to unambiguously discern if neutrinos, generated only as electron neutrinos in the core of the Sun, undergo flavor conversion during the path Sun-Earth.

Heavy water makes this possible providing both flavor-specific and flavor-independent neutrino reactions. The first, flavor-specific, reaction is the charged current (CC) reaction [figure omitted; refer to PDF] sensitive only to electron neutrinos. Due to the large energy of the incident neutrinos, the produced electron will be so energetic that it will be ejected at light speed, which is actually faster than the speed of light in water, therefore creating a burst of Cerenkov photons; after traveling throughout the water volume, they are revealed by the spherical array of photomultipliers instrumenting the detector. The amount of light is proportional to the incident neutrino energy, which can be inferred from the number of hits on the PMTs. From the hit pattern, also the angle of propagation of the light can be determined.

The second flavor-independent reaction is the so-called neutral current (NC) reaction [figure omitted; refer to PDF] whose net effect is just to break apart the deuterium nucleus; the liberated neutron is then thermalized in the heavy water as it scatters around. The reaction can eventually be observed due to gamma rays which are emitted when the neutron is finally captured by another nucleus. The gamma rays will scatter electrons, which produce detectable light via the Cerenkov process, in the same manner as discussed before.

The neutral current reaction is equally sensitive to all neutrino types; the detection efficiency depends on the neutron capture efficiency and the resulting gamma cascades. Neutrons can be captured directly on deuterium H 2 ( n , γ ) H 3 , but this is not very efficient. For this reason SNO has employed two separate systems to enhance the detection of NC interactions. In the so-called second SNO phase, Cl 35 has been added to the heavy water in form of 2 tons of NaCl and neutrons were detected through Cl 35 ( n , γ ) Cl 36 interaction. In the third SNO phase, the 36 proportional He 3 counters have been deployed in the core of the detector which enabled the neutron detection based on He 3 ( n , p ) H 3 interaction.

There is also a third reaction occurring in the detector, flavor independent as well, which is the electron scattering (ES): [figure omitted; refer to PDF]

This reaction is not unique to heavy water, being instead the primary mechanism in other light water detectors, like Kamiokande/Super-Kamiokande (see Section 3.2.2). Although this reaction is sensitive to all neutrino flavors, due to the different cross sections involved the electron neutrino dominates by a factor of six. The final state energy is shared between the electron and the neutrino, and thus there is very little spectral information from this reaction. On the other hand, good directional information can be obtained.

The general drawback affecting the Cerenkov technique is that, due to the feeble amount of light produced by the Cerenkov mechanism, the effective neutrino threshold is around 4-5 MeV, thus allowing the detection only of the high energy component of the solar flux, essentially the B 8 neutrinos.

The SNO experiment is now over; its architectural scheme was very simple (see Figure 1), aimed to get the most from the Cerenkov technique: 1000 tons of heavy water were contained in a thick transparent acrylic vessel, surrounded by an external layer of light water shielding from the gammas from the radioactivity in the rock. A spherical array of 10000 8[variant prime][variant prime] phototubes detected the light from both volumes of water. A key issue for the success of the experiment was the long standing effort throughout the construction and the operation phases to reduce the natural radioactivity in the target volume, not only in uranium and thorium, but also in the particular ubiquitous radon gas.

Figure 1: Conceptual architectural scheme of the SNO detector.

[figure omitted; refer to PDF]

As a result of this experimental effort, the multiple, clean, and almost background-free CC, NC, and elastic scattering detection of solar neutrinos provided the unambiguous and model independent proof that neutrinos from the Sun undergo flavor conversion. The specific "smoking gun" indication of the flavor-conversion process was obtained from the comparison of the depleted _{ν e} -only flux of the CC measurement with the all-flavor flux evaluated through the NC reaction. The first publication of this result in 2002 [31] nailed down definitively the explanation of the Solar Neutrino Problem.

3.2.2. Super-Kamiokande

As anticipated before, Super-Kamiokande [39], like its predecessor Kamiokande [40], is conceptually very similar to SNO, the major difference being the use of normal water instead of heavy water. Hence the neutrino detection occurs only via the scattering reaction off the electrons; the afterwards mechanism of Cerenkov light production and detection via an array of PMTs is equal to that already described for SNO.

Another major difference is the quantity of water employed, in total 50 ktons (observed by almost 13000 20[variant prime][variant prime] PMTs), which makes this detector the most massive among the neutrino oscillation experiments built so far. The sufficiently high statistics implied by this huge volume have made a fairly precise reconstruction of the spectrum of the scattered electrons possible, which plays an important role in the subsequent analysis for the interpretation of the data. With its huge mass Super-Kamiokande clearly outperforms the findings of the old Kamiokande (containing only 3000 tons of water), obtained in the data taking period from 1983 to 1994; however Kamiokande maintains a crucial historical role in the fields of neutrino oscillation and of astrophysical neutrinos, in this case with the detection of the neutrinos from the Supernova SN1987A, as witnessed by the 2002 Nobel Prize. In this context, it is appropriate to mention also another historically important Cerenkov experiment, the IMB (Irvine-Michigan-Brookhaven) detector [41], realized with 10000 tons of water, which shared with Kamiokande the success of detecting the SN1987A neutrinos.

The intrinsic high directionality of the scattering reaction, coupled to the directionality of the Cerenkov light, provides this experiment with a powerful tool to fight the background due to trace impurities of natural radioactivity dissolved in water, by associating the reconstructed direction of the Cerenkov photons with the angular position of the Sun. Clearly, this is done on top of the purification procedure of the light water, which as for SNO was focused generally on the whole natural radioactivity, but with special emphasis on the radon, which is the factor limiting the threshold at low energy.

An additional important analysis tool is the typical feature of the Cerenkov light to generate sharp Cerenkov rings in case of muon particles, while electrons make rings with fuzzy edges. Contrary to SNO, Super-Kamiokande is still currently taking data. The long history of this detector started in 1996 and evolved through four phases: the first phase lasted until a major PMT incident in November 2001 and produced a very accurate measure of the B 8 flux via the ES detection reaction. Phase II with reduced number of PMTs, from the end of 2002 to the end of 2005, confirmed with larger error the phase I measurement. After the refurbishment of the detector back to the original number of PMTs, the third phase lasted from the middle of 2006 up to the middle of 2008. Later on, an upgrade of the electronics brought the detector into its fourth, current phase. It is important to highlight the evolution of the energy threshold (total electron energy) in all the phases: 5 MeV in phase I, 7 MeV in phase II, 4.5 MeV in phase III, and 4 MeV for phase IV, thanks to the continuously ongoing effort to reduce the radon content in water.

Undoubtedly Super-Kamiokande played a central role in the long path which led to unveiling of the neutrino oscillation phenomenon, since it has been, and still is, a major player in three of the areas of investigation for neutrino oscillation, that is, those based on solar, atmospheric, and accelerator neutrinos. Actually, it was Super-Kamiokande that in 1998 [16] announced the epochal discovery of neutrino oscillations, which stemmed from the observed anomaly of the number of atmospheric muon neutrino events compared to that of electron neutrino events, and it was Super-Kamiokande that first confirmed the oscillation process with a beam of artificial (accelerator) muon neutrinos in the dedicated K2K experiment [42], which took place from 1999 to 2004. And nowadays this successful story continues with the T2K [43] experiment, another accelerator neutrino experiment which is the successor of K2K.

In the solar neutrino study the results provided by Super-Kamiokande are equally of great importance, as key ingredient of the joint analysis of all the experiments to ascertain the allowed regions of the oscillation parameters [10].

3.3. Scintillation Detectors

Scintillation detectors have a long and established tradition in the area of neutrino physics, starting from the Cowan-Reines Savannah River experiment [44], which performed the first neutrino detection ever. Other pioneer detectors of this kind which deserve to be mentioned for their historical role in the field (but not necessary in the oscillation study) are the Baksan Underground Scintillation Telescope (BUST) [45], which also detected the SN1987A neutrinos, the Liquid Scintillation Detector (LSD) at Mont Blanc [46], the Large Volume Detector (LVD) at Gran Sasso [47] devoted to Supernova search, and the Gosgen [48] and Bugey [49] reactor experiments.

In the following we focus our attention on the more recent implementations of this technique, for the realization of experiments which played a fundamental role in nailing down the neutrino oscillation properties.

3.3.1. Borexino

In the context of the solar neutrino research, the Borexino project was conceived and designed to detect in real time the low energy component of the solar flux, with special emphasis on the neutrinos coming from the Be 7 electron capture in the core of the Sun, exploiting as simple and effective mean to reveal the incoming particles their scattering reaction off the electrons of the target medium.

Specifically, Borexino is a scintillator detector [50] which employs as active detection medium a mixture of pseudocumene (PC, 1,2,4-trimethylbenzene) and PPO (2,5-diphenyloxazole, a fluorescent dye) at a concentration of 1.5 g/L. Because of its intrinsic high luminosity (50 times more than that in the Cerenkov technique) the liquid scintillation technology is extremely suitable for massive calorimetric low energy spectroscopy. The isotropic nature of the scintillation light does not allow inferring the direction of the incoming particles; it is therefore impossible, contrary to what happens in the Cerenkov experiments, to distinguish neutrino scattered electrons from electrons due to natural radioactivity by the association with the direction from the Sun. Thus the key requirement in the technology of Borexino is an extremely low radioactive contamination.

To reach ultralow operating background conditions in the detector, the design of Borexino, as shown in Figure 2, is based on the principle of graded shielding, with the inner scintillating core at the center of a set of concentric shells of increasing radiopurity. The scintillator mass (278 tons) is contained in a 125 μ m thick nylon inner vessel (IV) with a radius of 4.25 m. Within the IV a fiducial mass is software-defined through the estimated events position, obtained from the PMTs timing data via a time-of-flight algorithm.

Figure 2: Sketch of the Borexino experiment, highlighting its major components arranged according to a graded shielding design.

[figure omitted; refer to PDF]

A second nylon outer vessel (OV) with radius 5.50 m surrounds the IV, acting as a barrier against radon and other background contamination originating from outside. The region between the inner and outer vessels contains a passive shield composed of pseudocumene and 5.0 g/L (later reduced to 3.0 g/L) of DMP (dimethyl phthalate), a material that quenches the residual scintillation of PC so that spectroscopic signals arise dominantly from the interior of the IV.

A 6.85 m radius stainless steel sphere (SSS) encloses the central part of the detector and serves also as a support structure for the PMTs. The region between the OV and the SSS is filled with the same inert buffer fluid (PC plus DMP) which is layered between the inner and outer vessels.

Finally, the entire detector is contained in a tank (radius 9 m, height 16.9 m) filled with ultrapure water. The total liquid passive shielding of the central volume from external radiation (such as that originating from the rock) is thus 5.5 m of water equivalent (m.w.e.). The scintillator material in the IV was less dense than the buffer fluid by about 0.1% with the original DMP concentration of 5 g/L; this resulted in a slight upward buoyancy force on the IV, implying the need of thin, low-background ropes made of ultrahigh density polyethylene to hold the nylon vessels in place. This modest buoyancy was further reduced more than a factor 10 by removing via distillation a fraction of the total DMP content in the buffer: the process ended with a final DMP concentration of 3 g/L, still perfectly adequate to suppress the buffer scintillation, while at the same time implying less stress applied to the IV.

The scintillation light is viewed by 2212 8[variant prime][variant prime] PMTs uniformly distributed on the inner surface of the SSS. All but 371 photomultipliers are equipped with aluminum light concentrators designed to increase the collection efficiency of the light from the scintillator, concurrently minimizing the detection of photons not coming from the active scintillating volume. Residual background scintillation and Cerenkov light that escape quenching in the buffer are thus reduced. The PMTs without concentrators can be used to study this background, and help identify muons that cross the buffer and not the inner vessel.

Besides being a powerful shield against external backgrounds ( γ 's and neutrons from the rock), the water tank (WT) is equipped with 208 PMTs and acts as a Cerenkov muon detector. The muon flux, although reduced by a factor 10⁶ by the 3800 m.w.e. depth of the Gran Sasso Laboratory, is still significant (1.1 muon m^-2 h^-1 ) and an additional reduction (by about 10⁴ ) is necessary. Ultralow radioactive contamination is the distinctive feature of Borexino, achieved through a multiple strategy [51] that implied on one hand the careful selection and screening of all the construction materials and components and on the other the purification of the active scintillator to unprecedented purity levels (see Table 1).

Table 1: Radiopurity of the Borexino detector during phase 1 of the experiment.

Name	Source	Typical	Required	Achieved
14 C	Intrinsic PC	~1 0 - 12 g/g	~1 0 - 18 g/g	~2 × 1 0 - 12 g/g
238 U	Dust	1 0 - 5 - 1 0 - 6 g/g	<1 0 - 16 g/g	( 5.0 ± 0.9 ) × 1 0 - 18 g/g
232 Th	--	--	--	( 3.0 ± 1.0 ) × 1 0 - 18 g/g
7 Be	Cosmogenic	~3 × 1 0 - 2 Bq/ton	<1 0 - 6 Bq/ton	Not observed
40 K	Dust, PPO	~2 × 1 0 - 6 g/g (dust)	<1 0 - 18 g/g	Not observed
210 Po	Surface		<7 cpd/ton	May 07 : 70 cpd/ton
	Contamination	--	--	May 09 : 5 cpd/ton
222 Rn	Emanation, rock	10 Bq/L (air, water)	< 10 cpd/100 tons	< 1 cpd/100 tons
	--	100-1000 Bq/kg (rock)	--	--
39 Ar	Air, cosmogenic	17 mBq/m 3 (air)	< 1 cpd/100 tons	...a; 85 Kr
85 Kr	Air, nuclear weapons	~1 Bq/m 3 (air)	< 1 cpd/100 tons	30 ± 5 cpd/100 tons

Clearly, in this respect key factors are the many liquid purification and handling systems designed and installed to ensure the proper manipulation of the scintillator at the incredible degree of cleanliness demanded by the experiment. The exceptional low-background environment achieved in the core of the liquid scintillator allowed the unprecedented and precise sub-MeV measure of the Be 7 component of the solar neutrino flux, which to date is the only direct confirmation of the validity in such a low energy range of the MSW mechanism driving the oscillation of neutrinos produced in the core of the Sun.

Borexino has taken data during the so-called phase 1 (May 2007-July 2010) and started again to collect data (phase 2), after a further campaign of purification of the scintillator, in October 2011. The purification campaign succeeded to further reduce the residual contamination of the scintillator.

3.3.2. Other Scintillation Experiments

Other important scintillator-based experiments which provided milestone results for the understanding of the neutrino oscillation properties are KamLAND [19] and, more recently, Daya Bay [52], RENO [53], and Double Chooz [54]. While in term, of methodology all these experiments are very similar to Borexino, as far as detection criteria, techniques, and architectural scheme are considered, their specific characteristics are the measurement target, constituted by antineutrinos from reactors. KamLAND, in particular, is not close to any specific reactor but rather detects antineutrinos from a number of Japanese power plants located at an average distance of 200 km, thus performing a long-baseline test. Day Bay, RENO, and Double Chooz, instead, are located a 1 km from the reactor, acting thus as medium baseline experiments.

While in general the respective technology resembles closely that of Borexino, there are some variations in the type of liquid scintillator and in the material used for the balloon containing the liquid. The main difference with Borexino stems from the inverse beta reaction which is used to detect antineutrinos (in contrast with the scattering reaction adopted in Borexino to reveal neutrinos): [figure omitted; refer to PDF] After the occurrence of this interaction, the prompt signal is due to a positron which decelerates and then annihilates producing two 511 keV γ rays. The neutron thermalizes and is captured by a free proton, generating a typical 2.2 MeV gamma, the so-called delayed signal. The visible energy _{E vis} of the prompt signal is directly correlated with the kinetic energy of the incident antineutrino _{E ν ¯ e} = _{E vis} + 0.784 [MeV]. The mean time between the positron production and the neutron capture is about 200-260 μ s depending on the scintillator type, and therefore the tight time coincidence between the respective light signals prodices a correlated measurement which ensures a powerful discrimination of a true antineutrino detection with respect to the uncorrelated background events. This kind of signature greatly reduce the requirements for the suppression of the intrinsic radioactivity in the scintillator, marking the major difference between the technology of these reactor experiments and that employed for the solar neutrino detection in Borexino. For example, in KamLAND U 238 has been reduced to ( 1.5 ± 1.8 ) × 1 ^{0 - 19} g/g, Th 232 to ( 1.9 ± 0.2 ) × 1 ^{0 - 17} g/g, K 40 to a limit < 4.5 × 1 ^{0 - 18} g/g, Po 210 to ~2 mBq/m³ , Bi 210 to <1 mBq/m³ , and Kr 85 to ~0.1 mBq/m³ .

Historically, the measurement of KamLAND, together with that of SNO, closed the solar neutrino problem, showing unambiguously that also reactor antineutrinos undergo the oscillation phenomenon, while concurrently determining rather precisely the associated mass squared difference parameter Δ ^{m 21 2} and, jointly with the outputs of all other solar experiments, the mixing angle _{θ 12} .

Daya Bay and RENO (and in future Double Chooz, as well) have the additional characteristics of being equipped with a near detector, so that the far-near arrangement allows determining also the _{θ 13} mixing angle.

As additional remark of this section, it has to be emphasized that also some of the experiments whose outputs are used in the analysis concerning the existence of additional sterile state(s) beyond the established three-neutrino oscillation framework (see Section 4.2.4) are liquid scintillator setups. In particular, the LSND [55] and MiniBooNE [23] detectors, at the center of the current hot debate in this area, share essentially all the distinctive features of the other experiments belonging to the same technical "family." Specifically, LSND was a cylindrical tank containing 167 tons of scintillator viewed by 1220 8 inch PMTs, while MiniBooNE was based on a spherical detector geometry to contain 800 tons of scintillator, though still using a similar number of PMTs, 1280. The peculiarity of both setups was the exploitation of a special scintillator mixture able to produce a comparable amount of Cerenkov and true scintillation light.

The KARMEN experiment [56] was another player in this debate, but on the other side, since it did not detect the same hints of LSND and MiniBooNE. It was a segmented liquid-scintillator detector; the segmentation, technically the more distinctive feature of the set-up, was realized with 1.5 mm thick lucite sheets which ensured the transport of the light to the photomultipliers via total internal reflection. The detector was also instrumented with a veto employing plastic scintillator modules.

Finally, the pioneer MACRO detector [57] at Gran Sasso was, as well, based on segmented liquid-scintillator counters. Actually it comprised three subsystems, being additionally equipped with limited streamer tubes and nuclear track detectors, which altogether provided the experiment with the capability to detect the atmospheric neutrino oscillation phenomenon.

3.4. Further Techniques

To complete the illustration of the techniques adopted for the neutrino oscillation studies, a brief mention is due to other experiments which have shed light on important aspects of the field, while being not ascribable to any of the methodological categories described so far, starting with MINOS [58] and OPERA [59]. We briefly describe also the basic features of the near detectors of the already mentioned T2K experiment [43] (which, we remind, uses Super-Kamiokande as the far detector), as well as of the CHORUS [60], NOMAD [61], and ICARUS [62] experiments.

The MINOS experiment exploits two detectors to register the neutrino interactions: the near detector at Fermilab characterizing the neutrino beam (NuMI, neutrinos at the main injector, beam) is located about 1 km from the primary proton beam target, while the far detector performs similar measurements 735 km downstream. The far detector is located in Soudan, hosted in an inactive iron mine where it is positioned in a cavern excavated on purpose, 705 m underground (2070 meter water equivalent (m.w.e.)), 210 m below sea level.

The rationale of the experiment is to make comparisons between event rates, energies, and topologies at both detectors and to infer from those comparisons the relevant "atmospheric" oscillation parameters. The energy spectra and rates are measured separately for _{ν μ} and _{ν e} charged-current (CC) events, as well as for neutral current (NC) events.

Both the near and far MINOS detectors are steel-scintillator sampling calorimeters, equipped with tracking, energy, and topology measurement possibilities. Such a multiple capability is obtained by alternate planes of plastic scintillator strips and 2.54 cm thick magnetized steel plates.

The 1 cm thick by 4.1 cm wide extruded polystyrene scintillator strips are read out using wavelength-shifting fibers coupled to multianode photomultiplier tubes. Both detectors ensure equal transverse and longitudinal sampling for fiducial beam-induced events.

The far detector comprises 486 octagonal steel planes, with edge to edge dimension of 8 m, interleaved with planes of plastic scintillator strips. The total mass is 5400 tons; the set-up is arranged as two "supermodules" separated by a 1.15 m distance, individually equipped with an independently controlled magnet coil.

The near detector, consisting of 282 planes for a total mass of 980 tons, is located at the extreme of the NuMI beam facility at Fermilab, in a 100 m deep underground cavern under a 225 m.w.e. overburden. It exploits the high neutrino flux at this site to identify a relatively small target fiducial volume for selection of events to be employed for the near/far comparison. The upstream part of the detector, that is, the calorimeter portion, contains the target fiducial volume with all the planes instrumented. The downstream part, the spectrometer section dedicated to the measurement of the momenta of energetic muons, has only one plane every five instrumented with scintillator.

The core of the MINOS detector's active system is thus based on the technique of solid scintillator, whose main features are good energy resolution and hermiticity, excellent transverse segmentation, flexibility in readout, fast timing, simple and robust construction, long-term stability, ease of calibration, reliability, and, last but not least, low maintenance requirements. Furthermore, the whole setup met also safety and practicality of construction requirements.

The performances of both detectors rely on some key parameters which are the steel thickness, the width of scintillator strips, and the degree of readout multiplexing, which were carefully studied and optimized during the design phase.

The MINOS detectors represented a significant increase in size from previous fine grained scintillator sampling calorimeters, and therefore the relevant design and construction efforts ended up with important technical advancements in detector technology of general interest for the field of application of this technique.

This technological effort of the MINOS construction resulted in an impressive scientific success, which brought further evidence to the neutrino oscillation investigation performed by Super-Kamiokande and K2K, sharpening significantly the evaluation of the relevant "atmospheric" oscillation parameters.

The OPERA experiment was designed aiming at the direct observation of _{ν τ} appearance stemming from _{ν μ} [arrow right] _{ν τ} oscillation in a long baseline beam (dubbed CNGS) from CERN to the underground Gran Sasso Laboratory, at a distance of 730 km, where OPERA is located.

The design of OPERA was specifically tailored to identify the τ via the topological observation of its decay, reinforced by the kinematic analysis of the event. This goal is pursued through a hybrid apparatus based on two "pillars": real-time detection techniques ("electronic detectors") and the Emulsion Cloud Chamber (ECC) method. A detector based on the ECC approach is made of passive material plates, used as target, alternated with nuclear emulsion films employed as tracking devices, featuring submicrometric accuracy.

The submicrometric position accuracy, coupled to the adoption of passive material, allows for momentum measurement of charged particles through the detection of multiple Coulomb scattering, as well as for identification and measurement of electromagnetic showers, together with electron/pion separation.

In essence, the main advantage of the ECC technique is the unique property of combining a high accuracy tracker with the capability of performing precise measurements of kinematic variables.

OPERA scaled the ECC technology to an unprecedented size: the basic unit of the experiment is a "brick" realized with 56 plates of lead (1 mm thick) interleaved with nuclear emulsion films, for a total mass of 8.3 kg; 150000 of such target units have been assembled, amounting to an overall mass of 1.25 kton. The bricks are arranged in 62 vertical structures (walls), orthogonal to the beam direction, interleaved with planes of plastic scintillators.

The detector is made of two identical supermodules, each comprising 31 walls and 31 double layers of scintillator planes followed by a magnetic spectrometer.

The electronic detectors accomplish the twofold task to trigger the data acquisition, identify and measure the trajectory of charged particles, and locate the brick where the interaction occurred.

The momentum of muons is measured by the spectrometers, with their trajectories being traced back through the scintillator planes up to the brick where the track originates. In case of no muons observation, the scintillator signals produced by electrons or hadronic showers are used to predict the location of the brick that contains the primary neutrino interaction vertex. The selected brick is then extracted from the target and afterwards the two interface emulsion films attached on the downstream face of the brick are developed. If tracks related to neutrino interaction are observed in these interface films, the films of the brick are developed, too, following the tracks back by fully automated scanning microscopes until the vertex is located.

The analysis of the event topology at the primary vertex leads to the identification of possible τ candidates. Topologies of special interest might include one track that shows a clear "kink" due to the decay-in-flight of the τ (long decays) or an anomalous impact parameter with respect to the primary vertex (short decays) compatible with a decay-in-flight in the first lead plate. Once selected, such topologies are double-checked by a kinematic analysis at the primary and decay vertices.

The modular structure of the target ensures to extract only the bricks actually hit by the neutrinos, therefore achieving an efficient analysis strategy of the interaction, while at the same time minimizing the target mass reduction during the run.

In the overall structure of the OPERA detector each brick wall, containing 2912 bricks and supported by a light stainless steel structure, is followed by a double layer of plastic scintillators (Target Trackers, TT) that provide real-time detection of the outgoing charged particles. The instrumented target is further followed by a magnetic spectrometer, consisting of a large iron magnet instrumented with plastic Resistive Plate Chambers (RPC). The bending of charged particles inside the magnetized iron is measured by six stations of drift tubes (Precision Trackers, PT). Left-right ambiguities in the reconstruction of particle trajectories inside the PT are removed by means of additional RPC, with readout strips rotated by ±45° with respect to the horizontal plane and positioned near the first two PT stations.

What was defined before as a supermodule is actually an instrumented target together with its spectrometer.

Finally, two glass RPC planes mounted in front of the first target allow rejecting charged particles originating from outside the target fiducial region, coming from neutrino interactions in the surrounding materials.

As conclusive remark, OPERA is the first very large scale emulsion experiment: the 150000 ECC bricks include about 110000 m² emulsion films and 105000 m² lead plates; the scanning of the events is performed with more than 30 fully automated microscopes. The success of this impressive machine is witnessed by the unambiguous detection of 3 τ events, so far.

In an arrangement similar to MINOS, T2K [43] employs two near detectors located 280 m from the graphite proton target to measure the properties of the unoscillated neutrino beam.

The INGRID near detector comprises 16 modules, 14 of which are positioned in a cross configuration centered on the beam axis. They are made of iron and scintillator layers, allowing the measure of the neutrino rate and profile in the beam axis direction.

The ND280 off-axis near detector is located off the beam axis in the same direction as SK, being exploited to measure the properties of the un-oscillated off-axis beam. It consists of several subdetectors: the so-called Pi-Zero detector ( P ... D ) is a plastic scintillator-based detector optimized for ^{π 0} detection, followed by a tracking detector made of two fine grained scintillator detector units, in turn sandwiched between three time projection chambers. Both the P ... D and tracker are surrounded by electromagnetic calorimeters, including a module located immediately downstream of the tracker itself. The whole detector is located in a magnet with a 0.2 T magnetic field, serving also as mass for a side muon range detector.

Important predecessors of these efforts were two experiments carried out at CERN in the 90s, CHORUS [60] and NOMAD [61].

The active target of CHORUS was realized with nuclear emulsions (total mass of 770 kg). A scintillating fiber tracker was interleaved, both for timing and for extrapolating the tracks back to the emulsions. The set-up comprised also a hexagonal spectrometer magnet for momentum measurement, a high resolution spaghetti calorimeter for measuring hadronic showers, and a muon spectrometer. The scanning of the emulsions was performed with high-speed CCD microscopes.

NOMAD adopted drift chambers as target and tracking medium. The chambers were 44, located in a 0.4 T magnetic field, for a total fiducial mass of 2.7 tons. They were followed by a transition radiation detector (for e / π separation), by additional electron identification devices and by an electromagnetic lead glass calorimeter. The detector comprised also a hadronic calorimeter, 10 drift chambers for muon identification, and an iron-scintillator calorimeter of about 20 tons.

Finally, looking ahead to the future, it must be mentioned that a very promising technique potentially very useful for neutrino oscillation investigation is that based on liquid argon, developed through a very long research and development effort for the ICARUS detector [62]. Such liquid argon time projection chamber allows calorimetric measurement of particle energy together with three-dimensional track reconstruction from the electrons drifting in the electric field applied to a volume of sufficiently pure liquid argon. The technique, thus, successfully reproduces the extraordinary imaging features of a bubble chamber, but with the advantage of being a full electronic detector, potentially scalable to the huge masses required for the next round of experimental neutrino studies.

4. Experimental Results

The experimental results concerning the neutrino oscillations have been obtained studying neutrinos from several sources: solar and atmospheric neutrinos, reactor antineutrinos, and neutrino and antineutrino accelerator beams. The neutrino experiments make use of a variety of techniques: radiochemical methods, water and heavy water Cerenkov detectors, and liquid and plastic scintillators; in some detectors also streamer chambers and time projection chambers are used, in addition to nuclear emulsions.

The experiments can be classified as disappearance and appearance ones: the first are measuring a reduced flux of neutrinos having the same flavor as that at the source, while the second are looking for neutrinos of different flavor with respect to those emitted by the source.

4.1. Neutrino Sources

The Sun is one important source of neutrinos. Energy in the Sun is, in fact, produced by chains of nuclear reactions whose overall result is the conversion of hydrogen into helium [figure omitted; refer to PDF] Due to lepton number conservation, helium production is accompanied by the production of two-electron neutrinos. The total energy released in reaction (74) is Q = 26.73 MeV and only a small part of it (about 0.6 MeV on average) is carried away by the two neutrinos. The total flux of electron neutrinos arriving on Earth (if they do not oscillate) can be then estimated from the radiative flux K produced by the Sun on the Earth surface, obtaining [figure omitted; refer to PDF] Due to the eccentricity of the Earth's orbit, the solid angle from the Sun to the Earth changes during the year, and thus the solar neutrino flux shows a seasonal variation.

The interpretation of solar neutrino experiments requires a detailed knowledge of the solar neutrino spectrum; see Figure 3. Hydrogen burning in the Sun proceeds through two chains, namely, the p p chain and the CNO bicycle. At the temperature and density characteristic of the solar interior, hydrogen burns with ~99% probability through the p p chain that is predominantly initiated by the p + p [arrow right] d + ^{e +} + _{ν e} reaction. This reaction produces the so-called p p neutrinos which have a continuous spectrum extending up to E = 0.42 MeV and constitutes ≈ 90% of the total neutrino flux. Alternatively, the p p chain can originate with 0.23% probability from the reaction p + ^{e -} + p [arrow right] d + _{ν e} that produces the less abundant monochromatic pep neutrinos with energy E = 1.445 MeV.

Figure 3: The solar neutrino spectrum predicted by the SSM calculation of [15].

[figure omitted; refer to PDF]

The p p chain has three possible different branches ( p p -I, p p -II, and p p -III) whose relative rates depend on the central temperature of the Sun. In the p p -II termination, the electron capture reaction ^{e -} + Be 7 [arrow right] Li 7 + _{ν e} produces the monochromatic Be 7 neutrinos with energy E = 0.863 MeV. This value corresponds to transitions to the Li 7 ground state. With ~10% probability, Li 7 is produced in the first excited states together with a neutrino with energy E = 0.383 MeV. In the p p -III branch, the β -decay B 8 [arrow right] B 8 ^{e *} + ^{e +} + _{ν e} is responsible for the production of the B 8 neutrinos . The flux of B 8 neutrinos is extremely low, being approximately equal to 0.01 % of the total flux, but the spectrum extends up to a maximal energy E [approximate] 15 MeV.

In the CNO cycle, the overall conversion of four protons into helium is achieved with the aid of C, N, and O nuclei present in the Sun. The β -decays N 13 [arrow right] C 13 + ^{e +} + _{ν e} , O 15 [arrow right] N 15 + ^{e +} + _{ν e} , and, to a minor extent, F 17 [arrow right] O 17 + ^{e +} + _{ν e} produce the so-called N 13 , O 15 , and F 17 neutrinos , respectively, all together referred to as CNO neutrinos . These three components of the solar neutrino flux have continuous spectra extending up to E ≈ 1.2 , 1.7, and 1.7 MeV, respectively.

The predictions for each component of the solar neutrino flux are obtained by constructing a Standard Solar Model (SSM) which, according to the definition of [63], is a solution of the stellar structure equations (starting from a chemical homogeneous initial model) that reproduces, within uncertainties, the observed properties of the present Sun, by adopting physical and chemical inputs chosen within their range of uncertainties. In Table 2, we report the neutrino fluxes predicted by two recent SSM calculations that adopt two different assumptions for the admixture of heavy elements in the Sun. Namely, the model labeled GS98 is obtained by using the "old" composition from [29], while the model labeled AGSS09 adopts the "new" admixture of [30]. The reason to consider these two calculations is that, in recent years, a new solar problem, often referred to as solar metallicity puzzle , has emerged. The most recent determinations of the solar photospheric heavy-element abundances (among which [30]) have indicated, in fact, that the solar metallicity is lower by 30 to 40% than previous measurements [29]. However, the internal structure of SSMs calibrated against the newly determined solar surface metallicity does not reproduce the helioseismic constraints; see, for example, [64]. The experimental determination of the solar neutrino fluxes, besides providing crucial information for flavor neutrino oscillations, may help to shed light on the origin of these discrepancies.

Table 2: The predictions of SSMs implementing GS98 [29] and AGSS09 [30] admixtures. See [15] for details.

	AGSS09	GS98
pp	6.03 ( 1 ± 0.006 )	5.98 ( 1 ± 0.006 )
pep	1.47 ( 1 ± 0.012 )	1.44 ( 1 ± 0.012 )
hep	8.31 ( 1 ± 0.30 )	8.04 ( 1 ± 0.30 )
7 Be	4.56 ( 1 ± 0.07 )	5.00 ( 1 ± 0.07 )
8 B	4.59 ( 1 ± 0.14 )	5.58 ( 1 ± 0.14 )
13 N	2.17 ( 1 ± 0.14 )	2.96 ( 1 ± 0.14 )
15 O	1.56 ( 1 ± 0.15 )	2.23 ( 1 ± 0.15 )
17 F	3.40 ( 1 ± 0.17 )	5.52 ( 1 ± 0.17 )

The neutrino fluxes are given in units of 1 ^{0 10} ( p p ), 1 ^{0 9} ( ⁷ Be ), 1 ^{0 8} ( p e p , ¹³ N , ¹⁵ O ), 1 ^{0 6} ( ⁸ B , ¹⁷ F ), and 1 ^{0 3} ( h e p ) c ^{m - 2 s - 1} .

The atmospheric neutrinos are produced by cosmic rays, which collide with the atmosphere at its most external regions. In these collisions triggered mostly by the cosmic protons (plus a 5% of He and some minor contributions of heavier nuclei), pions and, at a much smaller rate, kaons are produced [65-69].

The main sources of the atmospheric neutrinos are the following reactions: [figure omitted; refer to PDF] As a consequence, the produced fluxes are approximately [figure omitted; refer to PDF] Moreover, due to the cosmic ray isotropy and the sphericity of the Earth, the up and down neutrino fluxes (i.e., having the zenith angle θ corresponding to cos ... θ < 0 and cos ... θ > 0 , resp.) are expected to have the same magnitude: [figure omitted; refer to PDF] where _{E ν x} is the energy of neutrino with flavor x .

The atmospheric neutrino flux can be evaluated with an uncertainty < 10% at 1 < E < 10 GeV, while at E < 1 GeV the error is larger. At E < 10 GeV, the relation (78) is valid within 2-3% errors. The accuracy worsens at larger energies due to kaon production. Equation (79) is confirmed at 1% at E < 1 GeV and has an uncertainty <1% at 1 GeV.

Nuclear reactors are a source of electron antineutrinos. The energy spectra of antineutrinos released in the fission of the main isotopes used as the fuel in reactor cores ( U 235 , U 238 , Pu 239 , and Pu 241 ) are shown in Figure 4. The reactor antineutrino flux is different from site to site and strongly depends on the presence of reactors in the neighborhoods. Its evaluation [70-73] has to take into account different reactor characteristics, some of them time dependent, as their thermal power and the power fractions of fuel isotopes. The reactor-detector distance has a strong influence on the shape of the oscillated, electron antineutrino energy spectrum. The mean energy of reactor antineutrinos which can be detected by the inverse beta-decay reaction given in (73) is about 4 MeV.

Figure 4: Energy spectra of antineutrinos released in the fission of the main isotopes ( U 235 (cyan), U 238 (red), Pu 239 (green), and Pu 241 (blue)) used as the fuel in the cores of nuclear power plants.

[figure omitted; refer to PDF]

Supernova explosions represent another possible source of neutrinos and antineutrinos of all flavors. The observation of neutrinos produced by a galactic Supernova could bring important information to comprehend the explosion mechanism and to study neutrino propagation in the dense Supernova environment. Supernova neutrino oscillations have a complex and interesting phenomenology; their potential in neutrino oscillation studies may be affected by the large uncertainties of the astrophysical Supernova models.

Finally, neutrinos and antineutrinos of various energies can be produced by accelerators. At CERN, FNAL, KEK, and Los Alamos Neutron Science Center, neutrino and antineutrino beams are produced for short and long baseline neutrino experiments.

4.2. The Neutrino Oscillation Study

The experimental study of the neutrino oscillations can be divided into several phases: the solar neutrino problem, the first proof of the oscillation phenomenon from atmospheric and solar neutrino experiments, precise measurements of the oscillation parameters Δ ^{m 21 2} and _{θ 12} by studying the nuclear-reactor antineutrinos, extension of the oscillation analysis to the low-energy neutrinos and the vacuum regime, confirmation of the oscillation phenomenon via disappearance and appearance experiments with accelerator beams, measurements of non-zero _{θ 13} , and finally indication of a third Δ ^{m 2} and therefore of a possible sterile neutrino.

4.2.1. Proof of the Neutrino Oscillation Phenomenon

The road towards the first understanding of the neutrino oscillation phenomenon passed several milestones.

(1) The Solar Neutrino Problem: An Apparent Deficit in the Solar Neutrino Flux. Pioneering experiments used the radiochemical techniques (see Section 3.1) applied to the observation of solar neutrinos; they are Homestake [24], GALLEX [26], and SAGE [27, 28]. These experiments are based upon the charged-current interaction of electron-flavor neutrino on a nucleus. Because the solar _{ν e} 's oscillate to different flavors, the experiments, which are sensitive only to electron neutrinos, detect a reduced number of events with respect to the expectations based on the Standard Solar Model (SSM). This lack of signal has been called "Solar Neutrino Problem" and the possible explanations were either a wrong description of the Sun by SSMs or the phenomenon of the "Neutrino Oscillations," a hypothesis introduced by Pontecorvo in 1957 [1-3].

The radiochemical experiments measure the integrated flux from the detection reaction threshold to the upper limit of the solar neutrino energy spectrum. Homestake measurements start from a threshold of ~0.814 MeV and, thus, do not probe the p p neutrino component of the solar neutrino flux; it observes 2.56 ± 0.23 SNU (SNU = solar neutrino unit equals to the neutrino flux producing 1 ^{0 - 36} captures per target atom per second) to be compared with the SSM expectation of 7 . ^{7 - 1.0 + 1.2} SNU [24]. A deficit of the solar neutrino signal was confirmed later by GALLEX, which, with a threshold at ~0.23 MeV, found 83 ± 19 ( stat ) ± 8 ( syst ) SNU to be compared to the expected 127 ± 7 SNU [26]. SAGE is still running and its results agree with the GALLEX's ones. The reduction is higher in the Homestake data (~67%) than in GALLEX (~35%). This difference is partly, but not completely, explained by the dependence of _{ν e} survival probability from the neutrino energy. A recent hypothesis of the existence of a light sterile neutrino [74] could explain the Homestake result.

The solar neutrino problem raised by the radiochemical experiments has been confirmed in 1991 by a real-time experiment based on the water Cerenkov technique, Kamiokande, detecting the ν - e elastic scattering [75]. The ν - e elastic scattering cross section σ is lower for μ , τ flavor neutrino than for the electron-flavor neutrino (for the muon flavor σ ( _{ν μ} - e ) ~ 1 / 7 σ ( _{ν e} - e ) ). Kamiokande finds the solar neutrino flux reduced by 40% with respect to what was expected by the SSM. The measured neutrino energy range includes only the B 8 solar neutrinos, because the threshold in Kamiokande is at ~5.0 MeV of the recoil-electron energy (which corresponds to ~5.2 MeV for the neutrino energy).

(2) The First Experimental Proof of Neutrino Oscillations. The experimental evidence for the existence of the neutrino oscillation phenomenon has been provided by three Cerenkov experiments (see Section 3.2), studying the atmospheric neutrinos with water (Kamiokande and Super-Kamiokande) and the solar neutrinos with heavy water (SNO). Here, we demonstrate the atmospheric neutrino measurements on the Super-Kamiokande results, since they are fully compatible with those of Kamiokande but are based on higher statistics.

Super-Kamiokande observed [16] an important discrepancy in the atmospheric _{ν μ} up and down fluxes, not observed, on the other hand, in the _{ν e} rates. The measured ratio of up and down _{ν μ} fluxes is well different from 1, contrary to what is expected in the absence of neutrino oscillations; see (80).

The results are summarized in Table 3 and in Figure 5. Super-Kamiokande detects the muons produced by _{ν μ} and the electrons produced by _{ν e} : muons and electrons are fast enough to produce a Cerenkov-light cone. The sub-GeV events are fully contained in the detector, while this is not the case for the multi-GeV events.

Table 3: The "up/down" asymmetry for muons and electrons observed in Super-Kamiokande [16]. Here, "up" refers to incident neutrinos within the zenith angle range - 1 < cos ... θ < - 0.2 and "down" within 0.2 < cos ... θ < 1 .

Source	"Up/down" asymmetry
Multi-GeV e -like	1.04 ± 0.03 ± 0.03
Multi-GeV μ -like	0.52 ± 0.05 ± 0.006
Sub-GeV e -like	1.09 ± 0.02 ± 0.03
Sub-GeV μ -like	0.65 ± 0.05 ± 0.001

Figure 5: Zenith angle distributions of muon and electrons for sub-GeV and multi-GeV data from Super-Kamiokande [16]. The hatched regions are the Monte Carlo expectations for no-oscillations. The solid black lines plot the best fits for _{ν μ} [arrow right] _{ν τ} oscillations; in the fit the overall flux normalization if left as a free parameter.

[figure omitted; refer to PDF]

The observed "up/down" asymmetry can be interpreted in terms of _{ν μ} [arrow right] _{ν τ} oscillations in vacuum. The best fit values of the oscillation parameters obtained from these data are 1.9 × 1 ^{0 - 3} e ^{V 2} < | Δ ^{m 23 2} | < 3.0 × 1 ^{0 - 3} eV² and sin 2 _{θ 23} > 0.90 (for Kamiokande data [76], the | Δ ^{m 23 2} | allowed region is ranging between 1.3 × 1 ^{0 - 2} and 2.95 × 1 ^{0 - 3} eV² ).

This oscillation effect can be understood on the basis of the oscillation length in vacuum, corresponding to a neutrino energy of ~1 GeV and to Δ ^{m 23 2} ~ 3 × 1 ^{0 - 3} eV² (see Section 2): [figure omitted; refer to PDF]

The "downgoing" neutrinos are reaching the detector after ~10 km from their production, while the "upgoing" neutrinos travel on average ~6000 km. As a consequence, the distance between production and detection for the downgoing neutrino is too short to observe a relevant flavor change. In Figure 6, the number of events versus L / E ( L is the distance between the neutrino production and the detector and E is the neutrino energy) is shown.

Figure 6: Plot of the number of events versus L / E ( L is the distance between the neutrino production and the detector and E is the neutrino energy) for the Super-Kamiokande data (points with error bars); the histogram is the result of the Monte Carlo simulation for atmospheric neutrino events without oscillations [17].

[figure omitted; refer to PDF]

These Super-Kamiokande results have demonstrated for the first time the existence of an oscillation phenomenon on the atmospheric neutrinos. It is essentially model independent and not influenced by any hypotheses assumed in the cosmic ray simulations. The results obtained by Kamiokande and Super-Kamiokande have been confirmed by MACRO experiment [57] with smaller statistics.

The SNO detector is a heavy water Cerenkov experiment installed in the Sudbury Inco mine (see Section 3.2.1). The use of a deuterium target allowed to study two independent neutrino interactions: charge current (CC, (70)) and neutral current (NC, (71)). In addition, the _{ν x} - e elastic scattering (72) has been detected. The data have been collected during three phases, characterized by different techniques to capture the neutron emitted in the NC reactions. The 5 MeV SNO threshold limits the detectable neutrinos to the B 8 component of the solar neutrino flux. Later SNO has repeated the analysis pushing down the threshold to 3.5 MeV (SNO LETA [77]). The B 8 - ν fluxes measured by SNO are summarized in Table 4.

Table 4: Fluxes of ⁸ B solar neutrinos measured by SNO in the three phases of the data taking.

Data set	Φ CC	Φ ES	Φ NC
× 1 0 6 c m - 2 s - 1	× 1 0 6 c m - 2 s - 1	× 1 0 6 c m - 2 s - 1
Phase 1 (306 live days) [31]	1.7 6 - 0.05 - 0.09 + 0.06 + 0.09	2.3 9 - 0.23 - 0.12 + 0.24 + 0.12	5.0 9 - 0.43 - 0.43 + 0.44 + 0.46
Phase 2 (391 live days) [32]	1.6 8 - 0.06 - 0.09 + 0.06 + 0.08	2.3 5 - 0.22 - 0.15 + 0.22 + 0.15	4.9 4 - 0.21 - 0.34 + 0.21 + 0.38
Phase 3 (385 live days) [32]	1.6 8 - 0.04 - 0.08 + 0.05 + 0.07	1.7 7 - 0.21 - 0.10 + 0.24 + 0.09	5.5 4 - 0.31 - 0.34 + 0.33 + 0.36

The CC interactions are produced only by _{ν e} , while the NC ones are triggered by all-flavor neutrinos. Therefore, it is clear by comparing the results from Table 4 that part of the _{ν e} produced in the nuclear reactions in the Sun's core has been transformed to other flavors. The final estimate of the B 8 neutrino flux from the NC reactions, obtained from a joint analysis of the three phases, is ( 5.25 ± 0.16 ( stat ^{) - 0.13 + 0.11} ( syst ) ) × 1 ^{0 6} cm^-2 s^-1 [38], in a good agreement with the SSM prediction of ( 4.59 ± 0.64 ) × 1 ^{0 6} cm^-2 s^-1 .

The SNO results can be interpreted as direct evidence of matter effects on neutrino oscillations; see Section 2. The second term in the r.h.s. of (56) is, in fact, not negligible due to the high electron density in the solar interior and to the relatively high energy of the neutrinos detected by SNO. Therefore, flavor oscillations are enhanced due to neutrino propagation through the Sun.

Super-Kamiokande also measured the solar neutrinos with a threshold of ~5 MeV. The _{ν e} - e elastic scattering gives a result of 2.32 ± 0.04 ( stat ) ± 0.05 ( syst ) × 1 ^{0 6} cm^-2 s^-1 [78], fully compatible with the SNO measurement.

In a two-neutrino analysis, the allowed regions in the Δ ^{m 21 2} versus ^{tan 2} _{θ 12} plane, obtained by a global fit of the radiochemical plus Cerenkov experiments are shown in Figure 7. The region at top right is called Large Mixing Angle (LMA). This name was originally assigned to distinguish it from another region (not shown in the figure) called Small Mixing Angle (SMA) which has been ruled out by recent solar and KamLAND data. The region at the bottom right is the so-called LOW region.

Figure 7: Allowed regions (68.27%, 95.45%, and 99.73% C.L.) in the oscillation parameter plane obtained fitting the Homestake + GALLEX + Super-Kamiokande + SNO data in the frame of the MSW model. Note that these results are based on solar neutrinos only, without considering KamLAND antineutrino data. From [18].

[figure omitted; refer to PDF]

Finally, other experiments were performed to look for neutrino oscillations with high-energy artificial neutrino beam at short baseline. In particular, the NOMAD [61] and CHORUS [60] experiments at CERN obtained a null result. On the other hand, the LSND experiment (see [55] and Section 4.2.4), which took data at Los Alamos with intermediate energy beam, obtained proof of neutrino oscillation, even if the result is controversial.

4.2.2. Checks and Refinements of the Solar and Atmospheric Neutrino Measurements

The results obtained by SNO and Super-Kamiokande on solar and atmospheric neutrinos have been confirmed by other experiments. We discuss here KamLAND, K2K, OPERA, and MINOS. These experiments, despite using different techniques and different neutrino sources (reactor _{ν ¯ e} , accelerator _{ν μ} beams), have E / L ratios, where E is the neutrino energy and L is the neutrino baseline, which permit to probe the same Δ ^{m 2} region as the "solar" and "atmospheric" data. OPERA is an appearance experiment; KamLAND, K2K, and MINOS are disappearance experiments. KamLAND detects the _{ν ¯ e} 's from the 55 Japanese nuclear reactors; K2K, MINOS, and OPERA study the _{ν μ} beam produced by the KEK, Fermilab, and CERN accelerators, respectively (see Section 3.4). It has to be recalled that it is possible to consider the ν ¯ data in the same framework of the ν results only if the CPT invariance is assumed.

The techniques used by these experiments are very different; see also Section 3. K2K is a long baseline experiment: a 1.3 GeV _{ν μ} beam is sent from KEK to Super-Kamiokande, 350 km apart. The E / L ratio is ~ 5 × ^{1 0 - 3} eV² , very close to the atmospheric Δ ^{m 23 2} range.

KamLAND studies the _{ν ¯ e} conversion by observing the _{ν ¯ e} produced by nuclear reactors with an average baseline L ~ 200 km. The E / L ratio falls just in the range of solar neutrinos.

The goal of the OPERA experiment is direct experimental observation of the _{ν τ} appearance in the _{ν μ} beam via the conversion _{ν μ} [arrow right] _{ν τ} . The E / L for OPERA is on average ~ 2.4 × 1 ^{0 - 2} eV² , partially in the range of atmospheric neutrinos. OPERA is expecting to observe no more than 5 to 8 τ decays.

MINOS is a long baseline experiment with near and far detectors. The measured energy spectrum in the far detector is compared to the predictions obtained on the basis of the near-detector data. In this way many sources of systematic uncertainty cancel out. The E / L is ~ 4 × ^{1 0 - 3} eV² , in the range of atmospheric neutrinos.

The KamLAND experiment had a big impact since it permitted to discriminate the possible solution of the solar neutrino problem. Its results [79], in fact, ruled out the LOW solution which was still allowed by the solar neutrino data only (see Figure 7) and restricted the LMA region. This is demonstrated in Figure 8 from [19]. In the frame of the two-neutrino approach the electron antineutrino survival probability can be written as [figure omitted; refer to PDF] In this approximation, the best fit parameters from the KamLAND data only are Δ ^{m 21 2} = ( 7.5 ^{8 - 0.13 + 0.14} ( stat ) ± 0.15 ( syst ) ) × 1 ^{0 - 5} eV² and ^{tan 2} _{θ 12} = 0.5 ^{6 - 0.07 + 0.10} ( stat ^{) - 0.06 + 0.10} ( syst ) [19]. Combining with solar neutrino data, the best fit parameters are Δ ^{m 21 2} = 7.5 ^{9 - 0.21 + 0.21} × 1 ^{0 - 5} eV² and ^{tan 2} _{θ 12} = 0.4 ^{7 - 0.05 + 0.06} .

Figure 8: Allowed regions [19] for the oscillation parameters Δ ^{m 21 2} and ^{tan 2} _{θ 12} from solar and KamLAND data. The allowed LMA area is the crossing between the solar and the KamLAND allowed regions.

[figure omitted; refer to PDF]

K2K [80, 81] studied both the _{ν μ} disappearance and a possible appearance of _{ν e} : the first to check the oscillation parameters in the atmospheric Δ ^{m 23 2} region and the second to study the _{θ 13} mixing angle. From the study of _{ν μ} disappearance, K2K confirmed the results of Super-Kamiokande with atmospheric neutrinos and obtained fully consistent values of the oscillation parameters 1.5 × 1 ^{0 - 3} < | Δ ^{m 23 2} | < 3.4 × 1 ^{0 - 3} eV² and ^{sin 2} 2 _{θ 23} > 0.92 . In the search of possible conversion of _{ν μ} [arrow right] _{ν e} K2K succeeded to extract only an upper limit for _{θ 13} .

OPERA has completed its data taking because the neutrino beam has been switched off at CERN, where the activity for the upgrading of the LHC energy has started. The data have been collected during 797 beam days and up to date the 2008-2009 data have been already analyzed, while the analysis of the 2010-2012 events is ongoing. OPERA has observed up to now three _{ν τ} candidates [20, 59], one of which is shown in Figure 9. The probability to have observed _{ν τ} appearance from _{ν μ} corresponds to ~3.5 σ C.L.

Figure 9: One of the OPERA _{ν τ} candidates [20]. In the primary vertex the _{ν τ} interacts producing a τ , which decays after 376 μ s into a muon plus neutrinos. In addition, the conversion of a γ , produced in the primary vertex, is visible.

[figure omitted; refer to PDF]

MINOS measures the muon flavor disappearance with a _{ν μ} [58] and a _{ν ¯ μ} [82, 83] enhanced beam. The detected charged-current interactions, _{ν μ} ( _{ν ¯ μ} ) + _{N 1} [arrow right] ^{μ -} ( ^{μ +} ) + _{N 2} , give the opportunities to reject the _{ν ¯ μ} ( _{ν μ} ) background in the _{ν μ} ( _{ν ¯ μ} ) beam, respectively, by analyzing the curvature of the reconstructed muon track in the magnetic calorimeter. The statistics of the neutrino data are more than an order of magnitude higher than the antineutrino data. The no-oscillation hypothesis is disfavored, in the case of _{ν ¯ μ} , at 6.3 σ C.L. By fitting the data in the context of two neutrino oscillations and using independent mass and mixing parameters for the neutrino and antineutrino case, the best fit results are:

(i) for _{ν μ} : | Δ ^{m 23 2} | = 2.3 ^{2 - 0.08 + 0.13} × 1 ^{0 - 3} eV² and ^{sin 2} ( 2 _{θ 23} ) > 0.90 (90% C.L.),

(ii) for _{ν ¯ μ} : | Δ ^{m ¯ 23 2} | = ( 3.3 ^{6 - 0.40 + 0.46} ( stat ) ± 0.06 ( syst ) ) × 1 ^{0 - 3} eV² and ^{sin 2} ( 2 _{θ ¯ 23} ) = 0.8 ^{6 - 0.12 + 0.11} ( stat ) ± 0.01 ( syst ) .

MINOS has analyzed also the neutral current interactions in order to investigate a possible active to sterile neutrino mixing. Because the neutral current interaction cross sections are the same for the three flavors, an observation of a neutral current event depletion between the near and far detectors could be due to the existence of a fourth sterile neutrino. MINOS found that the fraction of _{ν μ} , which may show a transition to a sterile neutrino, is <22% (90% C.L.) [82, 83].

4.2.3. Low-Energy Solar Neutrinos and the Oscillation in Vacuum

The analysis of the solar neutrinos by the Cerenkov experiments has been carried out with an energy threshold at ~5 MeV of detectable energy (for the incident neutrinos this threshold is slightly higher and depends on the reaction). Only SNO tried to push down the threshold to ~3.5 MeV [77], but the obtained results have large uncertainties. Thus, the spectrum analyzed by the Cerenkov technique corresponds to ~0.01% of the total solar spectrum and concerns the matter-enhanced neutrino oscillation.

The reason of such a high-energy threshold is the natural radioactivity, present in the environment and in the materials used to build the detectors. Two main families are present: Th 232 and U 238 ; the highest- Q (2.8 MeV) member is Tl 208 from the Th 232 decay chain. Therefore, in order to safely exclude the natural radioactivity (taking into account also the energy resolution) from the data, a high-energy threshold had to be applied.

The understanding of the solar neutrino oscillations which has been reached thanks to the radiochemical and the Cerenkov experiments can be demonstrated on a plot of the _{ν e} survival probability as a function of energy, shown in Figure 10. The grey band is the prediction of the LMA solution in the framework of the MSW model calculated using the best fit values of the oscillation parameters from a global fit of solar + KamLAND data: the thickness of the band takes into account the uncertainties of the oscillation parameters. The two plateaus, at the low and at the high energy regions, correspond to the oscillation in vacuum and in matter, respectively, as it is explained in Section 2. The intermediate region is called the transition region . The black experimental point in the high energy region is obtained from a proper average of the SNO + Super-Kamiokande data; the other two data points are from the radiochemical experiments. As it can be seen, the LMA solution in the frame of the MSW model is validated with good accuracy only for the matter-enhanced oscillation regime, while checks of increased precision are needed for the vacuum regime and the transition region.

Figure 10: Solar _{ν e} survival probability [21] as a function of neutrino energy. The data points are from the radiochemical and the water Cerenkov experiments. The grey band is the prediction of the LMA solution in the frame of MSW model.

[figure omitted; refer to PDF]

The study of the low-energy neutrinos, say below 2 MeV, needs a strong effort and the development of new techniques to strongly suppress the natural radioactivity background, to purify the active part of the detector and the shielding materials. Only one experiment succeeded to solve these problems: Borexino, installed at the Gran Sasso Underground Laboratory; see Section 3.3.1.

During phase 1 (May 2007-July 2010), Borexino succeeded to measure the Be 7 [84], p e p [22], and B 8 (with the lower threshold down to 3.2 MeV) [85] neutrino fluxes and to obtain an upper limit for the CNO neutrino flux [22]. In Table 5 we summarize the measured rates, while in Table 6 we compare the corresponding fluxes, calculated by using the best fit oscillation parameters, with the predictions obtained by the "low" (AGSS09 [30]) and the "high" metallicity (GS98 [29]) SSMs. The fluxes measured by Borexino and by the Cerenkov experiments are in good agreement with the SSM predictions but are unable to discriminate between high and low metallicity, mainly due to the experimental errors and uncertainties in solar model construction.

Table 5: Solar-neutrino rates as measured by Borexino.

Reaction in the Sun	Rate(counts/day/100 tons)
7 Be	46 ± 1.59 ( stat ) - 1.5 + 1.6 ( syst )
pep	3.13 ± 0.55 ( stat ) ± 0.23 ( syst )
CNO	< 1.4
8 B	0.22 ± 0.04 ( stat ) ± 0.01 ( syst )

Table 6: Comparison between the SSM predictions for the solar neutrino fluxes with high (GS98) and low metallicity (AGSS09) and the experimental results. The CNO flux corresponds to the sum of the ¹³ N, ¹⁵ O, and ¹⁷ F solar neutrino components.

	GS98	AGSS09	Experimental result
pep	1.44 ± 0.017	1.47 ± 0.018	1.6 ± 0.3 (Borexino)
7 Be	5.00 ± 0.35	4.56 ± 0.32	4.87 ± 0.24 (Borexino)
CNO	5.25 ± 0.79	3.76 ± 0.56	< 7.7 (95% C.L.) (Borexino)
8 B	5.58 ± 0.78	4.59 ± 0.64	5.2 ± 0.3 (SNO + SK + Borexino +KamLAND)
			5.25 ± 0.16 ( stat ) - 0.013 + 0.11 ( syst ) (SNO-LETA)

The neutrino fluxes are given in units of 1 ^{0 9} ( ⁷ Be ), 1 ^{0 8} ( p e p , CNO), and 1 ^{0 6} ( ⁸ B ) c ^{m - 2 s - 1} .

The impact of Borexino results on the determination of the solar _{ν e} survival probability is shown in Figure 11. In addition to the new measurements of Be 7 and p e p neutrino fluxes, the constraints on the p p flux have been much improved following the Be 7 - ν flux knowledge. Thus, the plateau corresponding to the vacuum regime is validated and through the Be 7 and the p e p neutrinos (this last even if with large errors) a check of the transition region has started. Finally, the B 8 analysis is extended to lower energies.

Figure 11: Solar _{ν e} survival probability [22] as a function of neutrino energy including all solar (with Borexino) experimental results. The grey band is the prevision of the LMA solution in the frame of MSW model.

[figure omitted; refer to PDF]

A further result of the Borexino phase 1 is the measurement [18] of the day/night asymmetry _{A D N} defined as [figure omitted; refer to PDF] for the rate of the Be 7 neutrinos, _{R N} and _{R D} being the corresponding rates during the day and night. It has been found to be null at ~1% error: _{A D N} = 0.001 ± 0.012 ( stat ) ± 0.007 ( syst ) . In a global fit with solar results only, this result is able to rule out the LOW region (see Figure 7 for the situation before Borexino), at 6.2 σ C.L., and, thus, without assuming the CPT invariance which is instead implicitly assumed when KamLAND antineutrino measurements are taken into account (see Figure 8).

All these results will be improved by Borexino during phase 2 because of the further radiocontaminants reduction and the effort to leave the detector untouched for at least three years. Phase 2 goals are: ( 1 ) the reduction of the uncertainties on the Be 7 and p e p neutrino fluxes; ( 2 ) the direct measurement of the p p neutrino flux; ( 3 ) either an improvement of the upper limit or a direct measurement of the CNO neutrino flux; ( 4 ) a measurement of the solar neutrino flux seasonal variation. The physics impact of these goals concerns the determination of the shape of the vacuum-to-matter transition region of the solar _{ν e} survival probability, which could be influenced by the Nonstandard neutrino Interactions (NSI) [86] or by the existence of an ultralight sterile neutrino [74, 87]. In the same direction goes the effort to measure with reduced errors the B 8 neutrino flux allowing an experimental point in the range 3-5 MeV.

4.2.4. A Third Δ ^{m 2} Range around 1 eV² ?

The possibility of a Δ ^{m 2} with a higher value than the solar and atmospheric ones has been considered in connection with the LSND [55] and MiniBooNE [23] data. Both these experiments are short baseline projects and the short distance between the neutrino source and the detector makes them impossible to observe oscillations driven by "atmospheric" mass difference Δ ^{m 23 2} or by the "solar" mass difference Δ ^{m 21 2} .

LSND has taken data during the periods 1993-1995 and 1996-1998 at Los Alamos Neutron Science Center with a _{ν ¯ μ} beam produced by ^{π +} and ^{μ +} decays, most of which at rest, and a _{ν μ} beam. The energy spectrum of both neutrinos and antineutrinos is broad, the maximum of the spectrum is at ~60 MeV for _{ν μ} and at ~45 MeV for _{ν ¯ μ} . The distance between the beam stop and the detector is 30 m. Therefore E / L is spanning around 1 eV² .

Strong effort has been devoted to reject the electron induced reactions, but in any case the electron background in the beam is very limited. In addition, the energy range is restricted within 20 < E < 200 MeV to study the oscillation _{ν ¯ μ} [arrow right] _{ν ¯ e} and within 60 < E < 200 MeV for _{ν μ} [arrow right] _{ν e} , in order to suppress various background sources. The _{ν ¯ e} 's are identified, as usual, through the inverse beta decay; see (73).

In the runs with _{ν ¯ μ} LSND found an excess of 117.9 ± 22.4 inverse beta-decay interactions. Subtracting from this sample 19.5 ± 3.9 events due to ^{μ -} in the beam and 20.5 ± 4.6 events due to _{ν ¯ μ} + p [arrow right] ^{μ +} + n , a final sample of 87.9 ± 22.4 ± 6 . 0 events remains. The maximum likelihood best fit for Δ ^{m 2} falls in the range 0.2-2.0 eV² [55]. LSND does not find any effect with the _{ν μ} beam.

A check of this result has been carried out by the experiment KARMEN with a detector and an energy similar to those of LSND, but with a distance from the neutrino source of 17.5 m. KARMEN, which is installed at the ISIS facility of the Rutherford Appleton Laboratory, did not find any evidence of _{ν ¯ μ} [arrow right] _{ν ¯ e} oscillation, but its sensitivity is lower than that of LSND. In any case, KARMEN succeeded to rule out a large part (but not all) of the ( Δ ^{m 2} , ^{sin 2} 2 θ ) region [88] allowed by LSND.

More recently a new collaboration, which some LSND members have joined, designed and carried out the experiment MiniBooNE [23] at Fermilab, just to check the LSND results. They use a _{ν μ} beam, peaking at 600 MeV of energy, and a _{ν ¯ μ} beam at 400 MeV, while the detector is located at 541 m from the beam target, thus with an E / L in the same range of LSND. The signature of a possible transition _{ν μ} [arrow right] _{ν e} and _{ν ¯ μ} [arrow right] _{ν ¯ e} is an excess of charged-current quasi-elastic events induced by _{ν e} and _{ν ¯ e} .

MiniBooNE finds 480 events passing the _{ν ¯ e} event selections at the neutrino energy range 200-3000 MeV, to be compared with the background expectation of 399.6 ± 20.0 ( stat ) ± 20.3 ( syst ) . Then the excess is 78.4 ± 28.5 ( 2.8 σ ) [23].

But, contrary to LSND, MiniBooNE observes also _{ν e} excess of 128.8 ± 43.4 in the _{ν μ} beam, in the range 200-475 MeV, while no excess has been observed above 475 MeV. The best fits for the oscillation parameters for the antineutrino mode, quoted by MiniBooNE in various papers, using a two-neutrino approach, vary from some hundreds to some tens of eV² for Δ ^{m 2} and from some tenths to few hundredths for ^{sin 2} 2 θ .

These oscillation parameters can only be allowed by assuming a fourth sterile (anti)neutrino, which does not interact but mixes with the active (anti)neutrinos. The present status of the art on this topic is summarized in Figure 12, where the allowed and excluded regions in the ( Δ ^{m 2} , ^{sin 2} 2 θ ) plane by LSND, MiniBooNE, and KARMEN are shown.

Figure 12: The allowed regions for _{ν ¯ μ} [arrow right] _{ν ¯ e} and _{ν μ} [arrow right] _{ν e} from the MiniBooNE data [23] are shown. In the case of antineutrinos, they are compared with the LSND allowed region and with the KARMEN exclusion plot (to the right of the dashed line).

[figure omitted; refer to PDF]

4.2.5. The Mixing Angle _{θ 13}

In the framework of three-neutrino oscillations, if we neglect for simplicity CP-violating effects (i.e., we set the CP phase δ = 0 ; see Section 2), there are five parameters, two squared mass differences Δ ^{m 21 2} , Δ ^{m 23 2} , and three angles _{θ 12} , _{θ 23} , and _{θ 13} . These parameters have been measured by the experiments described in the previous paragraphs, except the angle _{θ 13} for which, however, an upper limit was obtained. Therefore, in order to measure _{θ 13} , high statistics and very low systematic errors are needed. One important improvement can be reached using in the same experiment two detectors, one close to the neutrino source (near detector ) and another one (far detector ) at a distance of few kilometers (in case of low-energy neutrinos) or some hundreds of kilometers (in case of high-energy neutrinos). In a set-up with both near and far detectors many sources of systematic errors cancel out.

Three experiments, which succeeded to achieve a _{θ 13} measurement, are assembled with near and far detectors: Daya Bay, RENO, and T2K. Daya Bay [89] consists of 6 detectors, exposed at 6 nuclear reactors at 26 different distances ranging from 362 to 1925 m. RENO [53] also detects the reactor _{ν ¯ e} with the two detectors at 294 m and 1393 m from the center of a six-reactor array. In T2K [90], which is a second generation followup to the K2K, a _{ν μ} beam is sent off axis (2-3°, in order to reduce the beam energy below 1 GeV and then to have a proper E / L ) to Super-Kamiokande detector, 295 km away, with a near detector 280 m from the beam target. A further experiment, Double Chooz [54], has taken data with only a far detector (1050 m away) exposed to 2 reactors.

All detectors make use of the liquid scintillator technique, with the only exception of the magnetic tracking system of the T2K near detector. Daya Bay, RENO, and Double Chooz are disappearance experiments, while T2K is looking for appearance of _{ν e} in a _{ν μ} beam. All detectors are installed under some hundreds of meters of water equivalent overburden.

The experiments exposed to reactor antineutrinos determine the _{ν ¯ e} survival probability that in the context of three neutrino oscillations is given by (44). Neglecting the term proportional to _{S 12} , this can be written as [figure omitted; refer to PDF] The _{θ 13} values that are obtained from the experimental data are summarized in Table 7.

Table 7: Results of the _{θ 13} measurements from the three experiments detecting the _{ν - e} from nuclear reactors.

Experiment	N (measured)/ N (expected) events	sin 2 θ 13
Daya Bay	0.944 ± 0.007 ( stat ) ± 0.003 ( syst )	0.089 ± 0.010 ( stat ) ± 0.005 ( syst )
RENO	0.920 ± 0.009 ( stat ) ± 0.014 ( syst )	0.113 ± 0.013 ( stat ) ± 0.019 ( syst )
Double Chooz	--	0.097 ± 0.034 ( stat ) ± 0.034 ( syst )

T2K observed 21 _{ν e} candidates. The oscillation probability P ( _{ν μ} [arrow right] _{ν e} ) can be deduced from (39) obtaining approximately [figure omitted; refer to PDF] From the experimental results, one obtains ^{sin 2} _{θ 13} = 0.104 ± 0.060 ( stat ) ± 0.045 ( syst ) .

4.2.6. A Global Analysis

The most powerful tool to extract information on neutrino parameters is provided by global analysis in which all experimental neutrino data are fitted simultaneously in the context of three-neutrino oscillations. In fact, this approach provided the evidence (and the magnitude) for nonzero _{θ 13} before the direct experimental measurements were performed (see Fogli et al. [91]). We present here the status of global fits by reviewing the results of [33] (for a similar approach see also [92]). In this global analysis the following experimental information is included: the atmospheric neutrino data (Super-Kamiokande [93] phases 1-4), the long baseline accelerator experiments (K2K [42], disappearance and appearance data of MINOS [94], and T2K [90]), the reactor experiments (CHOOZ [95], Palo Verde [96], Double-Chooz [97], Daya Bay [98], RENO [53], and KamLAND [99]), solar radiochemical experiments (Homestake [24], GALLEX [34], and SAGE [27, 28]), and solar real-time experiments (Super-Kamiokande [100], SNO [38], and Borexino [84, 85]). The best fit oscillation parameters and the corresponding allowed regions are shown in Table 8. The two columns correspond to two different assumptions for the reactor neutrinos. In the first, the fit is carried out allowing for a free normalization of the reactor neutrino fluxes and including in the analysis the results of the reactor experiments with very short baseline L ...4; 100 m (RSBL), as Bugey, ROVNO, Krasnoyarsk, ILL, Gösgen, and SRP. The results presented in the second column are, instead, obtained by assuming the recent reactor fluxes calculated by Huber [73] and leaving out RSBL data. These two choices permit to explore the relevance of the reactor neutrino flux normalization in the extraction of neutrino parameters, taking also into account that recent calculations of the reactor fluxes as in Mueller et al. [72] and Huber [73] gives a deficit of about 3% with respect to the previous flux evaluations [70, 71]. The adopted choice slightly affects the _{θ 13} determination, partly due to the tension between the new fluxes and the RSBL data (the reactor anomaly): an increase of statistics of the Daya Bay and Reno will reduce the uncertainty connected with this choice.

Table 8: Neutrino oscillation parameters from a global fit from [33]. Details in text.

Parameter	Unit	(1)	(2)
sin 2 θ 12	--	0.30 2 - 0.012 + 0.013	0.31 1 - 0.013 + 0.013
θ 12	[ ° ]	33.3 6 - 0.78 + 0.81	33.8 7 - 0.80 + 0.82
sin 2 θ 23	--	0.41 3 - 0.025 + 0.037 ... 0.59 4 - 0.022 + 0.021	0.41 6 - 0.029 + 0.036 ... 0.60 0 - 0.026 + 0.019
θ 23	[ ° ]	40 . 0 - 0.1 . 5 + 2.1 ... 50 . 4 - 1.3 + 1.3	40 . 1 - 1.6 + 2.1 ... 50 . 7 - 1.5 + 1.2
sin 2 θ 13	--	0.022 7 - 0.0024 + 0.0023	0.025 5 - 0.0024 + 0.0024
θ 13	[ ° ]	8.6 6 - 0.46 + 0.44	9.2 - 0.45 + 0.41
δ	[ ° ]	30 0 - 138 + 66	29 8 - 145 + 59
Δ m 21 2	[ 1 0 - 5 e V 2 ]	7.5 0 - 0.19 + 0.18	7.5 1 - 0.15 + 0.21
Δ m 31 2 (NH)	[ 1 0 - 3 e V 2 ]	2.47 3 - 0.067 + 0.070	2.48 9 - 0.051 + 0.055
Δ m 32 2 (IH)	[ 1 0 - 3 e V 2 ]	- 2.42 7 - 0.065 + 0.042	- 2.46 8 - 0.065 + 0.073

For ^{sin 2} _{θ 23} , two disconnected 1 σ intervals are shown, the first one corresponds to the absolute minimum while the second one represents a secondary local minimum. The global analysis prefers a nonmaximal value of _{θ 23} ; this result is mostly driven by the MINOS _{ν μ} disappearance results.

We also see that there is a marginal sensitivity to δ provided by the combination of the MINOS disappearance results, the _{ν μ} [arrow right] _{ν e} data from long baseline experiments, and the atmospheric data.

Finally, while the sign of the mass difference Δ ^{m 21 2} is determined from matter effects in solar neutrino oscillation, the sign of the "atmospheric" mass splitting Δ ^{m 31 2} ≈ Δ ^{m 32 2} is not known. Correspondingly, we have two different options, the normal hierarchy (NH) for Δ ^{m 31 2} > 0 and the inverted hierarchy (IH) for Δ ^{m 31 2} < 0 , which provide a fit of very similar quality to that of the available data set.

5. Open Problems and Future Projects

In neutrino physics, there are still many open problems. Concerning neutrino oscillations, the more relevant issues are the determination of neutrino mass hierarchy, the measurement of the CP phase δ , and the precise evaluation of _{θ 23} (see also Section 4.2.6). Of course, there is a much vaster playground to observe considering such exciting possibilities as the existence of sterile neutrino and of non standard neutrino Interactions (NSI), the determination of the origin of neutrino mass (Majorana versus Dirac), the absolute mass of neutrinos, the role of neutrinos in cosmology, and possible connections to dark matter. The phenomenon of neutrino oscillations, to which this review is dedicated, provides still a unique tool to answer several of these questions. This last section is dedicated to some future prospects of physics of neutrino oscillations.

The CP violation was observed so far only in the quark sector and the CP violation in the leptonic sector is still a big unknown. Among the current and near-future experiments, T2K and NO ν A [101] have a limited sensitivity to CP violation via studying the _{ν μ} [arrow right] _{ν e} versus _{ν ¯ μ} [arrow right] _{ν ¯ e} appearance. The NO ν A experiment at NuMI beam is finalizing the construction phase at FNAL. A much improved sensitivity to δ and a strong discovery potential are expected only for the experiments of not so immediate future. T2HK (HK stands for Hyper-Kamiokande) will be a successor of T2K, with upgraded J-PARC beam and with the far detector represented by a next generation underground water Cerenkov detector Hyper-Kamiokande [102] with about 1 Mton of fiducial mass. The Hyper-Kamiokande construction is expected to start in 2016. The DAE δ ALUS [103], a phased neutrino physics program using cyclotron decay-at-rest neutrino sources would have a strong discovery potential when combined with Hyper-Kamiokande as the detector. The long baseline neutrino experiment (LBNE) [104], planning to use the strong neutrino beam from Fermilab travelling 1300 m baseline to 34 kton liquid argon time projection chamber has been recently approved. LAGUNA-LBNO (Large Apparatus for Grand Unification [105] and Neutrino Astrophysics and Long Baseline Neutrino Oscillations [106]) is a European long-baseline project using CERN neutrino beam. The _{ν e} [arrow right] _{ν μ} versus _{ν ¯ e} [arrow right] _{ν ¯ μ} appearance is a project of far future of ν -factories ( ν 's from decays of μ 's from accelerator) and β -beams ( ν 's from β -decays of light nuclei) which would have a decisive sensitivity to δ but face a problem of exceeding costs.

The long-baseline projects mentioned before for the CP violation search in the neutrino sector have also the ability to determine the neutrino mass hierarchy. For example, LBNE, jointly with T2K/NO ν A, can in principle address the mass hierarchy issue with a significance of more than 3 σ by 2030. Hyper-Kamiokande can get a similar discovery reach on a comparable timescale through the combination of atmospheric neutrino data with a shorter baseline measurement. In principle, the LAGUNA-LBNO project can afford a very high sensitivity, more than 5 σ , and with a relatively limited data taking period (few years), thanks to its very long baseline. Several other experiments have been discussed and proposed which may have the capability to test the neutrino mass hierarchy in a time frame shorter than that of the long-baseline projects. Among these alternatives, the most promising approaches seem to be reactor neutrinos (JUNO [107] formerly known as Daya Bay II) and atmospheric neutrinos in ice (PINGU [108] at IceCube) or water (ORCA [109]). Through a significant breakthrough in the technology and in the detector performance, JUNO is targeted to a potential sensitivity of more than 3 σ ( 4 σ ) , depending upon present and future uncertainties on Δ ^{m 32 2} . The experiment has been already approved in China and an international collaboration is being formed for its construction and operation. PINGU, as well as ORCA, could guarantee extremely good statistical sensitivity to the hierarchy, provided the systematic effects are under control and well understood. Actually, the estimates of the sensitivity may vary depending upon the choice of oscillation parameters and hierarchy; in an optimistic configuration, a 4 σ measurement could be obtained after 3 years of data.

In the solar neutrino physics, the most important open issues are the determination of the _{ν e} survival probability P ( _{ν e} [arrow right] _{ν e} ) in the transition region and the measurement of the CNO neutrino flux.

The LMA solution of the solar neutrino problem predicts an upturn in P ( _{ν e} [arrow right] _{ν e} ) at E [approximate] few MeV that corresponds to the transition from matter-enhanced oscillations (high-energy end) to vacuum-averaged oscillations (low-energy end); see Section 2 and the grey band in Figure 11. The observation of this feature would be the final confirmation of the LMA paradigm. However, the experiments that measure solar neutrinos have not observed it yet. The statistics in each individual experiment do not allow firm conclusions, but the effect indicating deviations from the LMA predictions is systematically present in all data. Super-Kamiokande, SNO-LETA, and BOREXINO seem to favor a flat distribution. In particular, Super-Kamiokande data disfavor the expected upturn at 1.3 to 1.9 σ C.L. In addition, the Homestake result is about 1.5 σ below the LMA prediction.

The shape of P ( _{ν e} [arrow right] _{ν e} ) in the transition region could be strongly influenced by the possible existence of the non standard neutrino interactions (NSI) [86] and/or of an ultralight sterile neutrino [74, 87]. As an example, in the presence of an ultralight sterile neutrino that mixes very weakly with active neutrinos, a dip in the P ( _{ν e} [arrow right] _{ν e} ) is expected in the transition region; its precise position is determined by the sterile neutrino mass and its width and depth depend on the mixing angle [74]. The possibility of neutrino NSI with other fermions has been predicted by several extensions of the SM, as for instance the left-right symmetric models and supersymmetric models with R -parity violation. The NSI can be described at low energy by effective four-fermion interactions: [figure omitted; refer to PDF] where _{G F} is the Fermi constant, α and β are the neutrino flavors, f and ^{f [variant prime]} are the electron or the light quarks, _{P C} is the chirality of the operator, C can be L or R , and finally ^{... α , β e , u , d} is a dimensionless number which, coupled with the weak coupling constant, parameterizes the strength of the interaction. Also in this case the shape of P ( _{ν e} [arrow right] _{ν e} ) in the transition region is influenced by the NSI and its study can limit the range of the parameter ... . Even more effective is its study in the ν - e elastic scattering. In order to constrain the shape of the transition region, it is critical to achieve a measurement as precise as possible of p e p , CNO, and B 8 (with a lower energy threshold down to 3 MeV) solar neutrinos. In addition to this, a measurement of CNO neutrinos is important also for other purposes.

The CNO bicycle contributes to ~1% of the energy emitted by the Sun. Despite being subdominant in the Sun, the CNO cycle has, however, a key role in astrophysics, being the prominent source of energy in more massive stars and in advanced evolutionary stages of solar-like stars. The evaluation of CNO efficiency is connected with various interesting problems, like, for example, the determination of globular clusters age from which we extract a lower limit to the age of the universe. At the moment, we still miss direct observational evidence for CNO energy generation in the Sun. The detection of CNO solar neutrinos would clearly provide a direct test of the CNO cycle efficiency. The measurement of the CNO solar neutrino flux can also provide clues to solve the so called "solar composition problem" (see Section 4.1). The flux is, in fact, directly related to the abundance of carbon, nitrogen, and oxygen in the core of the Sun and, so to the admixture of heavy elements. Thus, a determination of the CNO flux can help in solving the solar metallicity puzzle.

In the neutrino physics, one important problem is the possible existence of a fourth sterile neutrino. The LSND [55] and MiniBooNE [23] results (see Section 4.2.4) need independent checks, which could either rule out or confirm a third Δ ^{m 2} around 1 eV² . In addition, another two indications exist, which could favor the hypothesis of the fourth sterile neutrino: the reactor anomaly [72, 73] and the calibration of the radiochemical experiments with artificial sources [110]. The comparison of the antineutrino flux from nuclear reactors with the results of short baseline reactor experiments shows, in fact, a deficit of ~3.5%. A deficit has been also evidenced in the GALLEX [34] and SAGE [35] calibrations campaigns with a Cr 51 source (SAGE also with a Ar 37 source [36]). Each of these deficits is at ~ 3 σ level. A possible interpretation of these anomalies is connected to oscillations into new light sterile neutrino.

The existence of sterile neutrinos appears in many extensions of the Standard Model: they would be simply gauge singlets of the model. The simplest model (3 + 1 scheme) introduces only one sterile neutrino _{ν s} . In this scenario the four flavor eigenstates ( _{ν e} , _{ν μ} , _{ν τ} , _{ν s} ) mix through the matrix elements ( _{U e 4} , _{U μ 4} , _{U τ 4} , _{U s 4} ) with a fourth mass eigenstate _{ν 4} . The Δ ^{m i 4 2} ( i = 1,2 , 3 ) are supposed to be ≈ 1 eV² in order to provide explanation of the observed anomalies and, thus, they are much larger with respect to solar and atmospheric mass splittings.

The neutrino community favors new and decisive experimental tests on this matter [111]. Four experiments will try to address these problems in the near future: Borexino, Borexino-SOX, MicroBooNE, and Icarus-Nessie.

Borexino phase 2 has as a possible goal the experimental reproduction of the survival-probability transition region.

Borexino-SOX [112] is a project, in which the Borexino detector is taking data with a Cr 51 _{ν e} source installed in a small tunnel below it, at ~8 m from the detector center. The E / L is in the same range as that of the LSND and MiniBooNE experiments. It will look for the possible existence of a sterile neutrino, and will give an important check on the NSI, studying the ν - e elastic scattering.

MicroBooNE is an experiment based on a 70 tons fiducial volume liquid argon time projection chamber exposed to the NuMI neutrino beamlines at Fermilab. Its goal is to repeat the same measurements of MiniBooNE with high resolution at low energy, starting below 200 MeV [113].

Finally, Icarus-Nessie is an experiment using the Icarus liquid argon time projection chamber technique coupled to near and far spectrometers exposed to _{ν μ} and _{ν ¯ μ} beams with an E / L ~ 1 eV² , again to check the possible sterile neutrino existence [114].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] B. Pontecorvo, "Mesonium and anti-mesonium," Soviet Physics JETP , vol. 6, pp. 429, 1957.

[2] B. Pontecorvo, "Mesonium and anti-mesonium," Journal of Experimental and Theroretical Physics , vol. 33, pp. 549-551, 1957.

[3] B. Pontecorvo, "Inverse beta processes and nonconservation of lepton charge," Journal of Experimental and Theroretical Physics , vol. 34, pp. 247, 1958.

[4] L. Wolfenstein, "Neutrino oscillations in matter," Physical Review D , vol. 17, no. 9, pp. 2369-2374, 1978.

[5] S. P. Mikheyev, A. Yu. Smirnov, "Resonance amplification of oscillations in matter and spectroscopy of solar neutrinos," Soviet Journal of Nuclear Physics , vol. 42, no. 1, pp. 913-917, 1985.

[6] S. P. Mikheyev, A. Yu. Smirnov, "Resonant amplification of ν oscillations in matter and solar-neutrino spectroscopy," Il Nuovo Cimento C , vol. 9, no. 1, pp. 17-27, 1986.

[7] C. Giunti, C. W. Kim Fundamentals of Neutrino Physics and Astrophysics , Oxford University Press, New York, NY, USA, 2007.

[8] S. M. Bilenky, C. Giunti, W. Grimus, "Phenomenology of neutrino oscillations," Progress in Particle and Nuclear Physics , vol. 43, no. 1, pp. 1-86, 1999.

[9] A. Strumia, F. Vissani, "Neutrino masses, mixings and ...," http://arxiv.org/abs/hep-ph/0606054

[10] J. Beringer, J.-F. Arguin, R. M. Barnett, "Review of particle physics," Physical Review D , vol. 86, no. 1, 2012.

[11] S. M. Bilenk, J. Hosek, S. T. Petcov, "On scillations of neutrinos with dirac and majorana masses," Physics Letters B , vol. 94, no. 4, pp. 495-498, 1980.

[12] M. Blennow, A. Y. Smirnov, "Neutrino propagation in matter," Advances in High Energy Physics , vol. 2013, 2013.

[13] T. K. Kuo, J. Pantaleone, "Nonadiabatic neutrino oscillations in matter," Physical Review D , vol. 39, no. 7, pp. 1930-1939, 1989.

[14] T. K. Kuo, J. Pantaleone, "Neutrino oscillations in matter," Reviews of Modern Physics , vol. 61, no. 4, pp. 937-979, 1989.

[15] A. M. Serenelli, W. C. Haxton, C. Peña-Garay, "Solar models with accretion. I. Application to the solar abundance problem," The Astrophysical Journal , vol. 743, no. 1, pp. 24-28, 2011.

[16] Y. Fukuda, T. Hayakawa, E. Ichihara, "Evidence for oscillation of atmospheric neutrinos," Physical Review Letters , vol. 81, no. 8, pp. 1562-1567, 1998.

[17] Y. Ashie, J. Hosaka, K. Ishihara, "Evidence for an oscillatory signature in atmospheric neutrino oscillation," Physical Review Letters , vol. 93, no. 10, 2004.

[18] G. Bellini, J. Benziger, D. Bick, "Absence of a day-night asymmetry in the ⁷ Be solar neutrino rate in Borexino," Physics Leters B , vol. 707, no. 1, pp. 22-26, 2012.

[19] S. Abe, T. Ebihara, S. Enomoto, "Precision measurement of neutrino oscillation parameters with KamLAND," Physical Review Letters , vol. 100, no. 22, 2008.

[20] G. DeLellis Talk at the meeting of the Gran Sasso Laboratory Scientific Committee, private communication, March 2013

[21] G. Bellini, A. Ianni, G. Ranucci, "Borexino and solar neutrinos," Rivista del Nuovo Cimento , vol. 35, no. 9, pp. 481-537, 2012.

[22] G. Bellini, J. Benziger, D. Bick, "First evidence of pep solar neutrinos by direct detection in Borexino," Physical Review Letters , vol. 108, no. 5, 2012.

[23] A. A. Aguilar-Arevalo, B. C. Brown, L. Bugel, "Improved search for _{[upsilon] - μ} [arrow right] _{[upsilon] ¯ e}     oscillations in the MiniBooNE experiment," Physical Review Letters , vol. 110, no. 16, 2013.

[24] B. T. Clevel, T. Daily, R. Davis Jr., J. R. Distel, K. Lande, C. K. Lee, P. S. Wildenhain, J. Ullman, "Measurement of the solar electron neutrino flux with the homestake chlorine detector," The Astrophysical Journal , vol. 496, no. 1, pp. 505-526, 1998.

[25] J. N. Bahcall Neutrino Astrophysics , Cambridge University Press, Cambridge, UK, 1989.

[26] W. Hampel, J. Handt, G. Heusser, J. Kiko, T. Kirsten, M. Laubenstein, E. Pernicka, W. Rau, M. Wojcik, Y. Zakharov, R. V. Ammon, K. H. Ebert, T. Fritsch, D. Heidt, E. Henrich, L. Stielglitz, F. Weirich, M. Balata, M. Sann, F. X. Hartmann, E. Bellotti, C. Cattadori, O. Cremonesi, N. Ferrari, E. Fiorini, L. Zanotti, M. Altmann, F. V. Feilitzsch, R. Mößbauer, S. Wänninger, G. Berthomieu, E. Schatzman, I. Carmi, I. Dostrovsky, C. Bacci, P. Belli, R. Bernabei, S. D'Angelo, L. Paoluzi, M. Cribier, J. Rich, M. Spiro, C. Tao, D. Vignaud, J. Boger, R. L. Hahn, J. K. Rowley, R. W. Stoenner, J. Weneser, "GALLEX solar neutrino observations: results for GALLEX IV," Physics Letters B , vol. 447, no. 1-2, pp. 127-133, 1999.

[27] J. N. Abdurashitov, V. N. Gavin, S. V. Girin, "Measurement of the solar neutrino capture rate with gallium metal," Physical Review C , vol. 60, no. 5, 1999.

[28] J. N. Abdurashitov, V. N. Gavrin, V. V. Gorbachev, P. P. Gurkina, T. V. Ibragimova, A. V. Kalikhov, N. G. Khairnasov, T. V. Knodel, I. N. Mirmov, A. A. Shikhin, E. P. Veretenkin, V. E. Yants, G. T. Zatsepin, T. J. Bowles, S. R. Elliott, W. A. Teasdale, J. S. Nico, B. T. Cleveland, J. F. Wilkerson, "Measurement of the solar neutrino capture rate with gallium metal. III. Results for the 2002-2007 data-taking period," Physical Review C , vol. 80, no. 1, 2009.

[29] N. Grevesse, A. J. Sauval, "Standard solar composition," Space Science Reviews , vol. 85, no. 1-2, pp. 161-174, 1998.

[30] M. Asplund, N. Grevesse, A. J. Sauval, P. Scot, "The chemical composition of the sun," Annual Reviews of Astronomy and Astrophysics , vol. 47, pp. 481-522, 2009.

[31] Q. R. Amhad, R. C. Allen, T. C. Andersen, "Direct evidence for neutrino flavor transformation from neutral-current interactions in the sudbury neutrino observatory," Physical Review Letters , vol. 89, no. 1, 2002.

[32] B. Aharmim, S. N. Ahmed, J. F. Amsbaugh, "Measurement of the _{ν e} and total ⁸ B solar neutrino fluxes with the sudbury neutrino observatory phase-III data set," Physical Review C , vol. 87, no. 1, 2013.

[33] M. C. Gonzalez-Garcia, M. Maltoni, J. Salvado, T. Schwetz, "Global fit to three neutrino mixing: critical look at present precision," Journal of High Energy Physics , vol. 12, article 123, 2012.

[34] F. Kaether, W. Hampel, G. Heusser, J. Kiko, T. Kirsten, "Reanalysis of the GALLEX solar neutrino flux and source experiments," Physics Letters B , vol. 685, no. 1, pp. 47-54, 2010.

[35] J. N. Abdurashitov, V. N. Gavrin, S. V. Girin, V. V. Gorbachev, T. V. Ibragimova, A. V. Kalikhov, N. G. Khairnasov, T. V. Knodel, V. N. Kornoukhov, I. N. Mirmov, A. A. Shikhin, E. P. Veretenkin, V. M. Vermul, V. E. Yants, G. T. Zatsepin, Y. S. Khomyakov, A. V. Zvonarev, T. J. Bowles, J. S. Nico, W. A. Teasdale, D. L. Wark, M. L. Cherry, V. N. Karaulov, V. L. Levitin, V. I. Maev, P. I. Nazarenko, V. S. Shkol'nik, N. V. Skorikov, B. T. Cleveland, T. Daily, R. Davis Jr., K. Lande, C. K. Lee, P. S. Wildenhain, S. R. Elliott, J. F. Wilkerson, "Measurement of the response of a gallium metal solar neutrino experiment to neutrinos from a ⁵¹ Cr source," Physical Review C , vol. 59, no. 4, pp. 2246-2263, 1999.

[36] J. N. Abdurashitov, V. N. Gavrin, S. V. Girin, "Measurement of the response of a Ga solar neutrino experiment to neutrinos from a ³⁷ Ar source," Physical Review C , vol. 73, no. 4, 2006.

[37] M. A. Acero, C. Giunti, M. Laveder, "Limits on _{ν e} and _{ν e} - disappearance from Gallium and reactor experiments," Physical Review D , vol. 78, no. 7, 2008.

[38] B. Aharmin, S. N. Ahmed, A. E. Anthony, "Combined analysis of all three phases of solar neutrino data from the sudbury neutrino observatory," Physical Review C , vol. 88, no. 2, 2011.

[39] S. Fukuda, Y. Fukuda, T. Hayakawa, "The super-kamiokande detector," Nuclear Instruments and Methods in Physics Research Section A , vol. 501, no. 2-3, pp. 418-462, 2003.

[40] K. Hirata, T. Kajita, T. Kifune, "Observation of ⁸ B solar neutrinos in the kamiokande-II detector," Physical Review Letters , vol. 63, no. 1, pp. 16-19, 1989.

[41] R. Becker-Szendy, C. B. Bratton, D. Casper, "Search for muon neutrino oscillations with the Irvine-Michigan-Brookhaven detector," Physical Review Letters , vol. 69, no. 7, pp. 1010-1013, 1992.

[42] M. H. Ahn, E. Aliu, S. Andringa, "Measurement of neutrino oscillation by the K2K experiment," Physical Review D , vol. 74, no. 7, 2006.

[43] K. Abe, N. Abgrall, Y. Ajima, "Indication of electron neutrino appearance from an accelerator-produced off-axis muon neutrino beam," Physical Review Letters , vol. 107, no. 4, 2011.

[44] F. Reines, C. L. Cowan Jr., "Detection of the free neutrino," Physical Review , vol. 92, no. 3, pp. 830-831, 1953.

[45] E. N. Alekseev, L. N. Alekseeva, V. N. Bakatanov, M. M. Boliev, A. V. Voevodski..., V. A. Doguzhaev, V. N. Zakidyshev, V. A. Kozyarivski..., Y. F. Novosel'tsev, R. N. Novosel'tseva, V. Y. Poddubny..., V. Y. Petkov, O. V. Suvorova, A. K. Tsyabuk, A. E. Chudakov, "The baksan underground scintillation telescope," Physics of Particles and Nuclei , vol. 29, no. 3, pp. 254-256, 1998.

[46] M. Aglietta, G. Badino, G. F. Bologna, C. Castagnoli, W. Fulgione, P. Galeotti, O. Saavedra, G. Trinchero, S. Vernetto, V. L. Dadykin, V. B. Kortchaguin, P. V. Kortchaguin, A. S. Malguin, V. G. Ryassny, O. G. Ryazhskaya, V. P. Talochkin, G. T. Zatsepin, V. F. Yakushev, "Results of the liquid scintillation detector of the mont blanc laboratory," Il Nuovo Cimento C , vol. 9, no. 2, pp. 185-195, 1986.

[47] M. Aglietta, B. Alpat, E. D. Alyea, "The most powerful scintillator supernovae detector: LVD," II Nuovo Cimento A , vol. 105, no. 12, pp. 1793-1804, 1992.

[48] G. Zacek, F. V. Feilitzsch, R. L. Mössbauer, L. Oberauer, V. Zacek, F. Boehm, P. H. Fisher, J. L. Gimlett, A. A. Hahn, H. E. Henrikson, H. Kwon, J. L. Vuilleumier, K. Gabathuler, "Neutrino-oscillation experiments at the Gösgen nuclear power reactor," Physical Review D , vol. 34, no. 9, pp. 2621-2636, 1986.

[49] M. Abbes, B. Achkar, S. Ait-Boubker, "The bugey 3 neutrino detector," Nuclear Instruments and Methods in Physics Research Section A , vol. 374, pp. 164-187, 1996.

[50] G. Alimonti, C. Arpesella, H. Back, "The Borexino detector at the laboratori nazionali delran sasso," Nuclear Instruments and Methods in Physics Research Section A , vol. 600, no. 3, pp. 568-593, 2009.

[51] G. Alimonti, C. Arpesella, M. B. Avanzini, "The liquid handling systems for the Borexino solar neutrino detector," Nuclear Instruments and Methods in Physics Research Section A , vol. 609, no. 1, pp. 58-78, 2009.

[52] F. P. An, J. Z. Bai, A. B. Balantekin, "Observation of electron-antineutrino disappearance at daya bay," Physical Review Letters , vol. 108, no. 17, 2012.

[53] J. K. Ahn, S. Chebotaryov, J. H. Choi, "Observation of reactor electron antineutrinos disappearance in the RENO experiment," Physical Review Letters , vol. 108, no. 19, 2012.

[54] Y. Abe, C. Aberle, J. C. dos Anjos, "First measurement of _{θ 13} from delayed neutron capture on hydrogen in the Double Chooz experiment," Physics Letters B , vol. 723, no. 1-3, pp. 66-70, 2013.

[55] A. Aguilar, L. B. Auerbach, R. L. Burman, "Evidence for neutrino oscillations from the observation of _{ν e} - appearance in a _{ν μ} - beam," Physical Review D , vol. 64, no. 11, 2001.

[56] H. Gemmeke, V. Eberhard, H. Gemmeke, "The high resolution neutrino calorimeter KARMEN," Nuclear Instruments and Methods in Physcis Research A , vol. 289, no. 3, pp. 490-495, 1990.

[57] G. Giacomelli, A. Margiotta, "MACRO results on atmospheric neutrinos," Nuclear Physics B--Proceedings Supplements , vol. 145, no. 1-3, pp. 116-119, 2005.

[58] D. G. Michael, P. Adamson, T. Alexopoulos, "Observation of muon neutrino disappearance with the MINOS detectors in the NuMI neutrino beam," Physical Review Letters , vol. 97, no. 19, 2006.

[59] N. Agafonova, A. Aleksandrov, O. Altinok, "Search for oscillation _{ν μ}   [arrow right]   _{ν τ} with the OPERA experiment in the CNGS beam," New Journal of Physics , vol. 14, 2012.

[60] D. Macina, J. Konijn, R. G. C. Oldeman, "The CHORUS experiment," Nuclear Physics B--Proceedings Supplements , vol. 48, no. 1-3, pp. 183-187, 1996.

[61] J. Altegoer, M. Anfreville, C. Angelini, "The NOMAD experiment at the CERN SPS," Nuclear Instruments and Methods in Physics Research Section A , vol. 404, pp. 96-128, 1998.

[62] C. Rubbia, M. Antonello, P. Aprili, "Underground operation of the ICARUS T600 LAr-TPC: first results," Journal of Instrumentation , vol. 6, 2011.

[63] V. Castellani, S. Degl'Innocenti, G. Fiorentini, M. Lissia, B. Ricci, "Solar neutrinos: beyond standard solar models," Physics Report , vol. 281, no. 5-6, pp. 309-398, 1997.

[64] S. Basu, H. M. Antia, "Helioseismology and solar abundances," Physics Reports , vol. 457, no. 5-6, pp. 217-283, 2008.

[65] M. Honda, T. Kajita, K. Kasahara, S. Midorikawa, "Calculation of the flux of atmospheric neutrinos," Physical Review D , vol. 52, no. 9, pp. 4985-5005, 1995.

[66] M. Honda, K. Kasahara, K. Hidaka, S. Midorikawa, "Atmospheric neutrino fluxes," Physics Letters B , vol. 248, no. 1-2, pp. 193-198, 1990.

[67] G. Barr, T. K. Gaisser, T. Stanev, "Flux of atmospheric neutrinos," Physical Review D , vol. 39, no. 11, pp. 3532-3534, 1989.

[68] V. Agrawal, T. K. Gaisser, P. Lipari, T. Stanev, "Atmospheric neutrino flux above 1 GeV," Physical Review D , vol. 53, no. 3, pp. 1314-1323, 1996.

[69] T. K. Gaisser, T. Stanev, "An improved calculation of the atmospheric neutrino flux in the GeV range," in Proceedings of the 24th International Cosmic Ray Conference, vol. 1, pp. 694-697, Rome, Italy, August 1995.

[70] P. Vogel, G. K. Schenter, F. M. Mann, R. E. Schenter, "Reactor antineutrino spectra and their application to antineutrino-induced reactions. II," Physical Review C , vol. 24, no. 4, pp. 1543-1553, 1981.

[71] P. Huber, T. Schwetz, "Precision spectroscopy with reactor antineutrinos," Physical Review D , vol. 70, no. 5, 2004.

[72] Th. A. Mueller, D. Lhuillier, M. Fallot, "Improved predictions of reactor antineutrino spectra," Physical Review C , vol. 83, no. 5, 2011.

[73] P. Huber, "Determination of antineutrino spectra from nuclear reactors," Physical Review C , vol. 84, no. 2, 2012.

[74] P. C. de Holanda, A. Y. Smirnov, "Solar neutrino spectrum, sterile neutrinos, and additional radiation in the universe," Physical Review D , vol. 83, no. 11, 2011.

[75] K. S. Hirata, K. Inoue, T. Ishida, T. Kajita, K. Kihara, M. Nakahata, K. Nakamura, S. Ohara, N. Sato, Y. Suzuki, Y. Totsuka, Y. Yaginuma, M. Mori, Y. Oyama, A. Suzuki, K. Takahashi, M. Yamada, M. Koshiba, K. Nishijima, T. Suda, T. Tajima, K. Miyano, H. Miyata, H. Takei, Y. Fukuda, E. Kodera, Y. Nagashima, M. Takita, K. Kaneyuki, T. Tanimori, E. W. Beier, L. R. Feldscher, E. D. Frank, W. Frati, S. B. Kim, A. K. Mann, F. M. Newcomer, R. Van Berg, W. Zhang, "Real-time, directional measurement of ⁸ B solar neutrinos in the Kamiokande II detector," Physical Review D , vol. 44, no. 8, pp. 2241-2260, 1991.

[76] K. S. Hirata, T. Kajita, M. Koshiba, "Experimental study of the atmospheric neutrino flux," Physical Letters B , vol. 205, pp. 416-420, 1988.

[77] B. Aharmim, S. N. Ahmed, A. E. Anthony, "Low-energy-threshold analysis of the phase I and phase II data sets of the sudbury neutrino observatory," Physical Review C , vol. 81, no. 5, 2010.

[78] K. Abe, Y. Hayato, T. Iida, "Solar neutrino results in super-kamiokande-III," Physical Review D , vol. 83, no. 5, 2011.

[79] K. Eguchi, S. Enomoto, K. Furuno, "First results from KamLAND: ervidence for reactor arnti-neutrino disappearance," Physical Review Letters , vol. 90, no. 2, 2003.

[80] E. Aliu, S. Andringa, S. Aoki, "Evidence for muon neutrino oscillation in an accelerator-based experiment," Physical Review Letters , vol. 94, no. 8, 2005.

[81] S. Yamamoto, J. Zalipska, E. Aliu, "An improved search for _{ν μ}   [arrow right]   _{[upsilon] e} oscillation in a long-baseline accelerator experiment," Physical Review Letters , vol. 96, no. 18, 2006.

[82] P. Adamson, C. Andreopoulos, D. J. Auty, "First direct observation of muon antineutrino disappearance," Physical Review Letters , vol. 107, no. 2, 2011.

[83] P. Adamson, D. J. Auty, D. S. Ayres, "Improved search for muon-neutrino to electron-neutrino oscillations in MINOS," Physical Review Letters , vol. 107, no. 18, 2011.

[84] G. Bellini, J. Benziger, D. Bick, "Precision measurement of the ⁷ Be solar neutrino interaction rate in Borexino," Physical Review Letters , vol. 107, no. 14, 2011.

[85] G. Bellini, J. Benziger, S. Bonetti, "Measurement of the solar ⁸ B neutrino rate with a liquid scintillator target and 3 MeV energy threshold in the Borexino detector," Physical Review D , vol. 82, no. 3, 2010.

[86] Y. Minakata, C. Peña-Garay, "Solar neutrino observables sensitive to matter effects," Advances in High Energy Physics , vol. 2012, 2012.

[87] P. C. de Holanda, A. Y. Smirnov, "Homestake result, sterile neutrinos, and low energy solar neutrino experiments," Physical Review D , vol. 69, no. 11, 2004.

[88] B. Armbruster, I. M. Blair, B. A. Bodmann, "Upper limits for neutrino oscillations _{ν - μ}   [arrow right]   _{ν - e} from muon decay at rest," Physical Review D , vol. 65, no. 11, 2002.

[89] F. P. An, Q. An, J. Z. Bai, "Improved measurement of electron antineutrino disappearance at Daya Bay," Chinese Physics C , vol. 37, no. 1, 2013.

[90] K. Abe, N. Abgrall, Y. Ajima, "First muon-neutrino disappearance study with an off-axis beam," Physical Review D , vol. 85, no. 3, 2012.

[91] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo, A. M. Rotunno, "Evidence of _{θ 13} > 0 from global neutrino data analysis," Physical Review D , vol. 84, no. 5, 2011.

[92] G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo, A. M. Rotunno, "Global analysis of neutrino masses, mixings and phases: entering the era of leptonic CP violation searches," Physical Review D , vol. 86, no. 1, 2012.

[93] R. Wendell, C. Ishihara, K. Abe, "Atmospheric neutrino oscillation analysis with sub-leading effects in super-kamiokande I, II, and III," Physical Review D , vol. 81, no. 9, 2010.

[94] P. Adamson, D. S. Ayres, C. Backhouse, "Improved measurement of muon antineutrino disappearance in MINOS," Physical Review Letters , vol. 108, no. 19, 2012.

[95] M. Apollonio, A. Baldini, C. Bemporad, E. Caffau, F. Cei, Y. Déclais, H. De Kerret, B. Dieterle, A. Etenko, J. George, G. Giannini, M. Grassi, Y. Kozlov, W. Kropp, D. Kryn, M. Laiman, C. E. Lane, B. Lefièvre, I. Machulin, A. Martemyanov, V. Martemyanov, L. Mikaelyan, D. Nicolò, M. Obolensky, R. Pazzi, G. Pieri, L. Price, S. Riley, R. Reeder, A. Sabelnikov, G. Santin, M. Skorokhvatov, H. Sobel, J. Steele, R. Steinberg, S. Sukhotin, S. Tomshaw, D. Veron, V. Vyrodov, "Limits on neutrino oscillations from the CHOOZ experiment," Physics Letters B , vol. 466, no. 2-4, pp. 415-430, 1999.

[96] A. Piepke, F. Boehma, J. Busenitzb, "Final results from the Palo Verde neutrino oscillation experiment," Progress in Particle and Nuclear Physics , vol. 48, no. 1, pp. 113-121, 2002.

[97] Y. Abe, C. Aberle, J. C. dos Anjos, "Reactor _{ν e} - disappearance in the Double Chooz experiment," Physical Review D , vol. 86, no. 5, 2012.

[98] D. Dwyer, "Improved measurement of electronantineutrino disappearance at Daya Bay," Nuclear Physics B Proceedings , vol. 2, pp. 1-3, 2013.

[99] A. Gando, Y. Gando, K. Ichimura, "Constraint on _{θ 13} from a three-flavor oscillation analysis of reactor antineutrinos at KamLAND," Physical Review D , vol. 83, no. 5, 2011.

[100] J. Hosaka, K. Ishihara, J. Kameda, "Solar neutrino measurements in super-kamiokande-I," Physical Review D , vol. 73, no. 11, 2006.

[101] D. Perevalov, "First data from the NOvA experiment," in Proceedings of the 24th Workshop on Weak Interactions and Neutrinos (WIN'13), Natal, Brasil, 2013.

[102] M. Shiozawa, "The hyper-kamiokande project," Nuclear Physics B--Proceedings Supplements , vol. 237-238, pp. 289-294, 2013.

[103] C. Aberle, A. Adelmann, J. Alonso, "Whitepaper on the DAEdALUS program," http://arxiv.org/abs/1307.2949

[104] C. Adams, T. Akiri, M. Andrews, "Scientific opportunities with the long-baseline neutrino experiment," http://arxiv.org/abs/1307.7335

[105] D. Angus, A. Ariga, D. Autiero, "The LAGUNA design study--towards giant liquid based underground detectors for neutrino physics and astrophysics and proton decay searches," http://arxiv.org/abs/1001.0077

[106] M. Ghosh, P. Ghoshal, S. Goswami, S. K. Raut, "Synergies between neutrino oscillation experiments: an "adequate" configuration for LBNO," http://arxiv.org/abs/1308.5979

[107] Y.-F. Li, J. Cao, Y. Wang, L. Zhan, "Unambiguous determination of the neutrino mass hierarchy using reactor neutrinos," Physical Review D , vol. 88, no. 1, 2013.

[108] M. G. Aartsen, R. Abbasi, Y. Abdou, "PINGU sensitivity to the neutrino mass hierarchy," http://arxiv.org/abs/1306.5846

[109] P. Coyle, "ORCA:oscillation research with cosmics in the Abyss," in Proceedings of the Open Symposium on the European Strategy for Particle Physics (ESPP '12), of European Strategy Group Document under the Chair of the Scientific Secretary of the Strategy Session of CERN Council, Krakow, Poland

[110] C. Giunti, M. Laveder, "Statistical significance of the gallium anomaly," Physical Review C , vol. 83, no. 6, 2011.

[111] K. N. Abazajian, M. A. Acero, S. K. Agarwalla, "Light sterile neutrinos: a white paper," http://arxiv.org/abs/1204.5379

[112] G. Bellini, D. Bick, G. Bonfini, "SOX: short distance neutrino oscillations with Borexino," Journal of High Energy Physics , vol. 2013, article 38, 2013.

[113] H. Chen, J. Farrell, F. Lanni, "A proposal for a new experiment using the Booster and NuMI neutrino beamlines: MicrooBooNE," Fermilab, 2007

[114] M. Antonello, D. Bagliani, B. Baibussinov, "Search for anomalies in the neutrino sector with muon spectrometers and large LArTPC imaging detectors at CERN," in Proceedings of the Open Symposium on the European Strategy for Particle Physics (ESPP '12), pp. 1-10, Krakow, Poland, September 2012.