Phenomenology of Neutrino Oscillations
S.M. Bilenky
Joint Institute for Nuclear Research, Dubna, Russia, and
Institut für Theoretische Physik,
Technische Universität München,
D-85748 Garching, Germany
C. Giunti
INFN, Sezione di Torino,
and
Dipartimento di Fisica Teorica,
Università di Torino,
Via P. Giuria 1, I-10125 Torino, Italy, and
School of Physics, Korea Institute for Advanced Study,
Seoul 130-012, Korea
W. Grimus
Institute for Theoretical Physics, University of Vienna,
Boltzmanngasse 5, A-1090 Vienna, Austria
Abstract
|
This review is focused on neutrino mixing and neutrino
oscillations in the light of the recent experimental developments.
After discussing possible types of neutrino mixing for Dirac
and Majorana neutrinos and considering in detail
the phenomenology of neutrino oscillations in vacuum and matter,
we review all existing evidence and indications
in favour of neutrino oscillations that have been obtained in the
atmospheric, solar and LSND experiments. We present the results of
the analyses
of the neutrino oscillation data in the framework of mixing of three
and four massive neutrinos and investigate possibilities to test the
different neutrino mass and mixing schemes obtained in this way.
We also discuss briefly future neutrino oscillation experiments.
|
Journal:
Prog. Part. Nucl. Phys. 43, 1 (1999).
Preprint:
UWThPh-1998-61,
DFTT 69/98,
KIAS-P98045,
SFB 375-310,
TUM-HEP 340/98,
hep-ph/9812360.
Introduction
The strong evidence in favour of oscillations of atmospheric neutrinos
found by the Super-Kamiokande Collaboration
[1,2]
opened a new era in particle physics.
There is no doubt that new experiments are necessary to understand
the nature
of neutrino masses and mixing which are intimately connected with
neutrino oscillations,
but the first decisive step has been done:
massive and mixed neutrinos can now be considered as real physical objects.
The problem of neutrino mass has a long history.
Originally, Pauli considered
the neutrino as a particle with a small but non-zero mass
(smaller than the electron mass) [3]
and
the method for the measurement of the neutrino mass
through the investigation of
the
-spectrum near the end point was proposed
in the first theoretical
papers on
-decay of Fermi [4,5]
and Perrin [6].
The first experiments on the measurement of the neutrino mass,
based on the Fermi-Perrin method,
yielded the upper bound
[7]
which was improved in the fifties to
[8].
Therefore, it became evident that
the neutrino mass (if non-zero at all)
is much smaller than the electron mass.
This was the main reason that in 1957,
after the discovery of parity violation in
-decay,
the authors of the two-component theory of the neutrino
(Landau [9],
Lee and Yang [10],
Salam [11])
assumed that
the neutrino is a massless particle,
the field of which is either
a left-handed field
or a right-handed field
.
In 1958, Goldhaber et al. [12]
measured the helicity of the neutrino.
The result of this experiment was in agreement with the
two-component neutrino theory and
it was established that the neutrino field is
.1The results of the experiment of Goldhaber et al.
could not exclude, however, the
possibility of a small neutrino mass.
In the V
A theory
(Feynman and Gell-Mann [15],
Sudarshan and Marshak [16])
the Hamiltonian of weak interactions
contains the left-handed component of the neutrino field
, and also the left-handed components of all massive fields.
Therefore the possibility for the neutrino to be nevertheless a
massive particle became more natural [17] after the
confirmation of the V
A theory.
In 1957, B. Pontecorvo [18,19]
proposed the idea that
the state of neutrinos
produced in weak interaction processes
is a superposition of states of two
Majorana neutrinos [20] with definite masses
(analogous to the states
and
which are the superposition of
and
,
the states of particles with definite masses and widths).
In this way, B. Pontecorvo arrived at the hypothesis of neutrino oscillations
(analogous to
oscillations).
At that time only one type of neutrino was known.
The possibility of mixing of the two species of neutrinos
and
was considered in
Ref. [21]. All possible types of neutrino oscillations for
this case were investigated by Pontecorvo
in 1967 [22].
Gribov and Pontecorvo proposed in 1969 [23]
the first phenomenological theory of neutrino mixing
and oscillations.
In this theory, the two left-handed neutrino fields
and
are linear combinations of
the left-handed
components of the fields of Majorana neutrinos with definite masses
and
the neutrino mass term contains
only the left-handed fields
and
.
In 1976,
neutrino oscillations were considered in the scheme of mixing of two Dirac
neutrinos based on the analogy between quarks and leptons
[24,25]
and in the same year in the general Dirac-Majorana scheme
[26] (for later works see
Refs. [27,28,29,30]).
The theoretical arguments in favour of non-zero neutrino masses and mixing
are based on the models beyond the Standard Model
(see, for example, Ref. [31]).
In such models
the fields of quarks, charged leptons and neutrinos
are grouped in the same multiplets
and the generation of the masses of quarks and charged leptons
with the Higgs
mechanism as a rule provides also non-zero neutrino masses.
In 1979,
the see-saw mechanism for the generation
of neutrino masses was proposed
[32,33,34].
This mechanism
connects the smallness of neutrino masses with the possible violation of
lepton number conservation at a very large energy scale.
At present the effects of neutrino masses and mixing are investigated
in many different experiments. There are three types of experiments in
which the effects of small neutrino masses (say, of the order of 1 eV
or smaller) and mixing can be revealed (for a review and references
see Ref. [35]):
- Neutrino oscillations experiments.
- Experiments on the search for neutrinoless double
-decay.
- Experiments on the measurement of the electron neutrino
mass
with the precise
investigation of the high energy part of
the
-spectrum of
H.
Three indications in favour of neutrino masses and mixing have been found
so far. These indications were obtained in the following experiments:
- Solar neutrino experiments
(Homestake [36,37,38],
Kamiokande [39,40,41],
GALLEX [42,43],
SAGE [44,45],
Super-Kamiokande [46,47,48]).
- Atmospheric neutrino experiments
(Super-Kamiokande [1,2],
Kamiokande [49],
IMB [50],
Soudan [51], MACRO [52]).
- The accelerator LSND experiment [53,54].
Many other neutrino oscillation experiments with neutrinos
from reactor and accelerators
did not find any evidence for neutrino oscillations. In
the experiments on the search for neutrinoless double
-decay
no indications for non-zero neutrino
masses were found (see Section 6.2).
The present upper bound for the electron neutrino mass
obtained in the Troitsk experiment
[55] is 2.7 eV (see also the Mainz experiment [56]).
The upper limits on the masses of
and
are 170 keV (90% CL) and 18.2 MeV (95% CL),
respectively [35].
Neutrinos play an important role in cosmology and astrophysics
and many bounds on neutrino properties can be derived in this
context. For reviews see, e.g., Refs. [31,57].
In this review we discuss the phenomenological theory of
neutrino mixing (Section 2),
neutrino oscillations in vacuum (Section 3),
neutrino oscillations and transitions in matter (Section 4) and
experimental data and results of analyses of the data (Sections 5 and 6).
We also consider the implications
of the existing experimental results on neutrino oscillations
for experiments in preparation. After the conclusions (Section 7) we
discuss some properties of Majorana neutrinos and fields in Appendix A.
We hope that this review will be useful not only for the physicists
that are working
in the field but also for those who are interested in this exciting field of
physics.
In many cases we present not only results but also
derivations of the results.
For those who start to study the subject we refer to the books
Refs. [58,59,60,31,61]
and the reviews
Refs. [62,63,64,65,66,67,68,69,70,71,72,73].
Neutrino Mixing
All the numerous data on weak interaction processes
are perfectly well described by
the Standard Model
[74,75,76].
The standard weak interactions
are due to the coupling of quarks and leptons
with the gauge
and
vector bosons,
described by the charged-current (CC)
and neutral-current (NC)
interaction Lagrangians
| |
|
 |
(2.1) |
| |
|
 |
(2.2) |
Here
is the SU(2)
gauge coupling constant,
is the weak angle and
the charged and neutral currents
and
are given by the expressions
| |
|
 |
(2.3) |
| |
|
 |
(2.4) |
where the
are the physical charged lepton fields with
masses
and
we have written explicitly only the terms containing the neutrino fields.
The flavour neutrinos
,
are determined by CC weak interactions:
for example, the
is the particle produced in the decay
and so on.
The number of light flavour neutrinos is given by the
invisible width of the
boson [77] in the Standard Model.
It was shown in the famous LEP experiments on the
measurement of the invisible width of the
boson that the number of the
light flavour neutrinos is equal to 3. The most recent experimental
value of the number of neutrino flavours is
[35],
showing that are no other neutrino flavours than the well-known
,
,
.
The CC and NC interactions conserve the
electron
,
muon
and tau
lepton numbers,
which are assigned as shown in Table 2.1.
There are no indications in favour of violation
of the law of conservation of lepton numbers in weak processes
and
very strong bounds on the probabilities of the lepton number violating
processes have been obtained from the experimental data.
The most stringent limits (90% CL) are
(see Ref. [35]):
| |
|
 |
(2.5) |
| |
|
 |
(2.6) |
| |
|
 |
(2.7) |
| |
|
 |
(2.8) |
| |
|
 |
(2.9) |
According to the neutrino mixing hypothesis
[18,19,21],
the conservation of the lepton numbers is only approximate.
It is violated because of
non-zero neutrino masses and neutrino mixing.
In the case of neutrino mixing,
the left-handed flavour neutrino
fields
are superpositions
of
the left-handed components
of the fields of neutrinos with definite masses
:
 |
(2.10) |
where
is a unitary mixing matrix. The number of massive neutrino fields
could be equal or more than 3. If
is larger than 3 there are
sterile neutrinos that do not take part in the
standard weak interactions (2.3) and (2.4).
It is well-established that
quarks take part in CC weak interactions in mixed form with the V
A
current
 |
(2.11) |
where
is the Cabibbo-Kobayashi-Maskawa mixing matrix
[78,79].
The relation (2.10)
is analogous to the mixing in Eq.(2.11).
However, between the mixing of neutrinos and quarks
there can be a fundamental difference.
In fact,
quarks are four-component Dirac particles: quarks and antiquarks have
opposite electric and baryonic charges.
Instead,
neutrinos are electrically neutral particles.
If the total lepton charge
 |
(2.12) |
is conserved,
neutrinos with definite masses
are four-component Dirac particles
like quarks
(in this case a neutrino differs from
an antineutrino by the opposite value of
).
If the total lepton number (2.12)
is not conserved,
massive neutrinos are truly neutral two-component Majorana particles.
These possibilities are realized in different
models and correspond to different neutrino mass terms.
Dirac mass term
A Dirac neutrino mass term can be generated by the Higgs mechanism
with the
standard Higgs doublet which is responsible for the generation of the masses of
quarks and charged leptons.2In this case
the neutrino mass term is given
by3
 |
(2.13) |
and
is a complex
matrix.
The Dirac mass term (2.13) can be diagonalized
taking into account that the complex matrix
can be written as
 |
(2.14) |
where
and
are unitary matrices and
is a positive definite diagonal matrix,
,
with
.
Hence,
the Dirac mass term (2.13)
can be written in the diagonal form
 |
(2.15) |
with
 |
(2.16) |
Therefore,
in the case of the Dirac mass term (2.13)
the three flavour fields
(
)
are linear unitary combinations of
the left-handed components
of three fields
of neutrinos with masses
(
). On the other hand, we also have
 |
(2.17) |
but these right-handed fields do not occur in the standard weak
interaction Lagrangian. Therefore, right-handed singlets are
sterile and are not mixed in this scheme with the active
neutrinos.4
The field
is a Dirac field if not only the mass term
(2.13) but also the total Lagrangian
is invariant under the global U(1) transformation
 |
(2.18) |
where the phase
is the same for all neutrino
and charged lepton fields. Then,
using Noether's theorem,
one can see that
the invariance of the Lagrangian under the transformation
(2.18) implies that
the total lepton charge
is conserved. Therefore,
is the
quantum number that distinguishes a neutrino from an antineutrino.
The Dirac mass term (2.13)
allows processes like
,
.
However,
the contribution of neutrino mixing to the probabilities of such
processes is negligibly small
[81,82,83].
The unitary
mixing matrix
can be written in terms of
3 mixing angles and 6 phases.
However, only one phase is measurable
[79].
This is due to the fact that in the Standard Model the only term in the Lagrangian
where the mixing matrix
enters
is the CC interaction Lagrangian (2.1).
With neutrino mixing the lepton charged-current
is given by
 |
(2.19) |
Because the phases of Dirac fields are arbitrary one can eliminate from
the mixing matrix
five phases
by a redefinition of the
phases of the charged lepton and neutrino fields, with
only one physical phase remaining in
.
The presence of this phase causes the violation of CP invariance
in the lepton sector.
A convenient parameterization of the mixing matrix
is the one proposed in Ref. [84]:
 |
(2.20) |
where
and
and
is the CP-violating phase.
Since in the parameterization (2.20)
of the mixing matrix
the CP-violating phase
is associated with
,
it is clear that CP violation
is negligible in the lepton sector if the mixing angle
is small.
More generally,
it is possible to show that if any of the elements
of the mixing matrix is zero,
the CP-violating phase can be rotated away by a suitable
rephasing of the charged lepton and neutrino
fields.5
Dirac-Majorana mass term
If none of the lepton numbers is conserved and both,
left-handed flavour fields
(
) and
sterile right-handed gauge singlet fields
(
)
enter into the mass term
we have the so-called Dirac-Majorana mass term
 |
(2.21) |
with
Here
,
and
are complex matrices and the indices
run over
values.
Let us stress that
the number of right-handed singlet fields
could be different from the number of neutrino flavours
(see [88,30]).
The charge-conjugate fields are defined by
 |
(2.25) |
where
is the charge conjugation matrix (see Appendix A).
Notice that
and
are right-handed and left-handed components, respectively.
Indeed, with
one
finds
 |
(2.26) |
and the analogous reasoning holds for
.
Furthermore, we have
 |
(2.27) |
The matrices
and
are symmetric. This can be shown with the help of the relation
,
which follows from the fact that
is an antisymmetric
matrix and from the
anticommutation property of fermion fields. An analogous relation holds for
the right-handed fields.
Then, using Eq.(2.27),
one obtains
 |
(2.28) |
Let us introduce the left-handed column vector
 |
(2.29) |
where
is the column vector of right-handed fields.
With the relation
 |
(2.30) |
the Dirac-Majorana mass term
can be written in the form
 |
(2.31) |
with the symmetric
matrix
 |
(2.32) |
The complex symmetric matrix
can be diagonalized with the help of a unitary matrix
[89,90]:
 |
(2.33) |
with
and
.
Using this relation the mass term (2.31)
can be written as
 |
(2.34) |
with
 |
(2.35) |
The fields
satisfy the Majorana condition
 |
(2.36) |
Therefore,
in the general case of a Dirac-Majorana mass term
the fields of particles with definite masses
are Majorana fields.
Majorana particles are particles with spin
and all charges
equal to zero (particle
antiparticle).
Some properties of
Majorana fields are discussed in Appendix A.
We want to emphasize that it is natural that the diagonalization of the mass
term (2.31) leads to fields of Majorana particles with definite
masses: the mass
term (2.31) for the case of a general matrix
is not invariant
under any global phase transformation.
In other words,
in the general case of the Dirac and Majorana mass term
there are no conserved quantum numbers
that allow to
distinguish a particle from its antiparticle.
From Eq.(2.35), for the left-handed components
of the neutrino fields we have
the mixing relations
| |
|
 |
(2.37) |
| |
|
 |
(2.38) |
Thus,
in the case of a Dirac-Majorana mass term the flavour fields
are linear unitary combinations of the left-handed components
of the Majorana fields
of neutrinos with definite masses.
These components are also connected with
the sterile fields
through the relation (2.38).
If the masses
are small,
the mixing relations (2.37) and (2.38) imply that
flavour neutrinos
,
and
can
oscillate into sterile states that are quanta of
the right-handed fields
.
We will discuss such transitions in Section 6.
Let us stress that in order to have all three terms in the
Dirac-Majorana mass term (2.31) not only a Higgs doublet but
also a Higgs triplet
[91,92,93]
and a Higgs singlet are necessary.
Thus,
it can be generated only in
the framework of models beyond the Standard Model.
A typical example is the SO(10) model
(see, for example, [31]). For a discussion of
radiative corrections to the Dirac-Majorana mass term see
Ref. [94].
Up to now we did not make any special assumption about the
Dirac-Majorana mass term (2.31).
Let us now consider the possibility
that CP invariance in the lepton sector holds.
In this case we have
 |
(2.39) |
where
is the operator of CP conjugation,
and
.
For the fields
we take the usual CP transformation
 |
(2.40) |
where
is a diagonal matrix of phase factors.
The mass term
can be written in the form
 |
(2.41) |
Using Eq.(2.40) we obtain
 |
(2.42) |
Hence,
the CP invariance condition (2.39)
is satisfied if
 |
(2.43) |
Up to now
was an arbitrary diagonal matrix of phase factors.
If we choose
,
in this case
CP invariance in the lepton sector implies that
 |
(2.44) |
Taking into account that
is a symmetric matrix,
the condition (2.44) is equivalent to
 |
(2.45) |
The real symmetric matrix
can be diagonalized with
the transformation
 |
(2.46) |
where
is an orthogonal matrix (
)
and
is a diagonal matrix,
.
The eigenvalues
can be positive or negative and the neutrino masses are thus given by
. Therefore, we write the eigenvalues as
 |
(2.47) |
where
is the sign of the
eigenvalue of the matrix
.
Then the relation (2.46) can be rewritten in the form
 |
(2.48) |
with
 |
(2.49) |
and
.
Here
is the diagonal matrix with elements
.
Therefore,
if CP is conserved in the lepton sector,
the neutrino mixing matrix has the simple structure
shown in Eq.(2.49), which implies that
 |
(2.50) |
The CP transformation of the Majorana field
is given by (see Appendix A)
 |
(2.51) |
where
is the CP parity of the Majorana field
which is
(see Ref. [65]).
From Eq.(2.51) we have
 |
(2.52) |
This relation and Eqs.(2.35), (2.40) and (2.50) imply
that the CP parity of the Majorana field
is
[95,96,97]
 |
(2.53) |
Indeed, we have
 |
(2.54) |
Majorana mass term
If only left-handed neutrino flavour fields
(
) enter into the Lagrangian
we can write down the mass term [23,98]
 |
(2.55) |
where
is a symmetric complex matrix (see previous section).
Then we have for the mixing
 |
(2.56) |
where
is the field of Majorana neutrinos with mass
.
In this case
the number of massive Majorana neutrinos is equal
to the number of lepton flavours (three).
Note that the generation
of the Majorana mass term (2.55) requires an
enlargement of the scalar sector of the Minimal Standard Model
[99]:
with a Higgs triplet like in the Majoron models
Ref. [91,92,93] Majorana masses are
obtained at the tree level, whereas
radiative generation is possible at the 1-loop level with a singly
charged Higgs singlet (plus an additional Higgs doublet)
[100,101,102,103,104,31])
or at the 2-loop level with a doubly charged Higgs singlet (plus an
additional singly charged scalar) [105].
Since the Majorana condition (2.36)
does not allow the rephasing
of the neutrino fields,
only three of the six phases in the
mixing matrix
can be absorbed into
the charged lepton fields
in the charged current (2.19).
Therefore,
in the Majorana case the mixing matrix contains three CP-violating phases
[28,29,106] in contrast
to the single CP-violating phase of the Dirac case discussed in
Section 2.1.
However, the additional CP-violating phases in the Majorana case
have no effect
on neutrino oscillations
in vacuum
[28,29,106]
(see Section 3)
as well as in matter
[107].
The one-generation case
Let us consider the Dirac-Majorana mass term in the simplest case of one
generation.
We have
with
 |
(2.58) |
For simplicity we assume CP invariance
in the lepton sector (see Subsection 2.2).
In this case
,
and
are real parameters.
In order to diagonalize the matrix
,
let us write it in the form
 |
(2.59) |
with
 |
(2.60) |
The eigenvalues of the matrix
are
 |
(2.61) |
Furthermore, we have
 |
(2.62) |
where
and
is the orthogonal matrix
 |
(2.63) |
with the mixing angle
given by
 |
(2.64) |
Hence, the matrix
can be written as
 |
(2.65) |
where
and
![\begin{displaymath}
m'_{1,2}
=
\frac{1}{2}
\left[
( m_R + m_L )
\mp \sqrt{ ( m_R - m_L )^2 + 4 \, m_{\mathrm{D}}^2 }
\right]
\,.
\end{displaymath}](img185.gif) |
(2.66) |
The eigenvalues of the matrix
are real but can have positive or negative sign.
Let us write them as
 |
(2.67) |
where
and
is the sign of the
eigenvalue of the matrix
.
As shown in Eq.(2.53),
the CP parity of the Majorana field
with definite mass
is
.
The relation (2.65)
can be written in the form
 |
(2.68) |
with
and
 |
(2.69) |
with
.
Using now the general formulas obtained in
Section 2.2,
we have the mixing relation
 |
(2.70) |
with
 |
(2.71) |
Therefore,
the three parameters
,
,
are related with
the mixing angle
and the neutrino masses
by the relations (2.64) and (2.66), (2.67).
The signs of the eigenvalues of
determine
the CP parities of the massive Majorana fields
.
In the framework of CP conservation, the relations obtained so far are general.
In the following part of this section
we consider some particular cases
with special physical significance.
The see-saw mechanism
Let us consider
the Dirac-Majorana mass term (2.57)
and assume [32,33,34] that
,
,
where
is the mass of a quark or a charged lepton
of the same generation,
and
.
In this case,
from the relations (2.64), (2.66) and (2.67)
we have
These relations imply that approximately
 |
(2.75) |
and the Majorana fields
,
are connected with the fields
and
by
![\begin{displaymath}
\nu_1 = i \, [ \nu_L - (\nu_L)^c ]
\,,
\qquad
\nu_2 = \nu_R + (\nu_R)^c
\,.
\end{displaymath}](img204.gif) |
(2.76) |
The mechanism which we consider here is called see-saw mechanism
[32,33,34].
It is based on the assumption that the conservation of
the total lepton number
is violated by the right-handed Majorana mass term at the scale
that is
much larger than the scale of the electroweak
symmetry breaking.
Several models which implement the see-saw mechanism are possible
(see, for example, [31,66,68]
and references therein)
and the scale
depends on the model.
This scale could be as low as the TeV scale
(for example, in left-right symmetric models
[108,109,31])
or an intermediate scale,
or as high as the grand unification scale
or even the Plank scale
.
The great attractiveness of the see-saw model
lies in the fact that,
through the relation (2.73),
it gives an explanation
of the smallness of neutrino masses with respect to the masses of other
fundamental fermions.6
In the case of three generations
the see-saw mechanism leads to
a spectrum of masses of Majorana particles
with three light neutrino masses
and three very heavy masses
(
)
of the order of the scale of violation of the lepton numbers.
This is realized if the mass matrix
(2.32)
has
,
 |
(2.77) |
and
is such that all its eigenvalues
are much bigger than the elements of
.
In this case
the mass matrix is block-diagonalized
(up to corrections of order
)
by the unitary transformation
 |
(2.78) |
with
 |
(2.79) |
The matrices for the light and heavy masses are given by
[112,113]
 |
(2.80) |
The mass eigenvalues of the light neutrinos
are determined by the specific form
of
and
.
Note that in left-right symmetric models and in SO(10) models the
matrix
in the big matrix
(2.32)
can be important (see, e.g.,
Ref. [114]).
Two simple possibilities are discussed in the literature
(see Refs. [115,116]):
- If
,
where
is the identity matrix,
one obtains the quadratic see-saw,
 |
(2.81) |
and the light neutrino masses are given by
 |
(2.82) |
where
is the mass of a quark or a charged lepton of the
generation.
In this case the neutrino masses
scale as the squares of the
masses
:
 |
(2.83) |
- If
,
where
characterizes the scale of
,
one obtains the linear see-saw
(see, for example, Ref. [117]),
 |
(2.84) |
and the light neutrino masses are given by
 |
(2.85) |
In this case the neutrino masses
scale as the
masses
:
 |
(2.86) |
Let us stress that in any case
the see-saw mechanism implies
the hierarchical relation
 |
(2.87) |
for the three light Majorana neutrino masses.
Effective Lagrangians
In the Standard Model
without right-handed singlet neutrino fields
there are no renormalizable interactions that
give masses to the neutrinos
after the spontaneous breaking of the
symmetry with the Higgs doublet mechanism.
However,
there is a general belief that the Standard Model
is the low-energy manifestation of a more complete theory
[118,119]
(for reviews see Refs. [120,31]).
The effect of this new theory is to induce in the Lagrangian
of the Standard Model
non-renormalizable interactions
which preserve the
symmetry
above the electroweak symmetry breaking scale,
but violate
the conservation of lepton and baryon numbers
(see Ref. [121] and references therein).
These non-renormalizable interactions are operators of dimension
and must be multiplied by coupling constants that have dimension
,
where
is a mass scale characteristic of the new theory.
It is clear that the dominant effects at low energies
are produced by the operators with lowest dimension.
In the Standard Model the lepton number non-conserving operator
with minimum dimension that can generate a neutrino mass is7[122,123,124,125,126,127]
 |
(2.88) |
where
is a
matrix of coupling constants,
are the Pauli matrices,
are the standard leptonic doublets
and
is the standard Higgs doublet:
 |
(2.89) |
The effective operator in Eq.(2.88)
has dimension five
and the coupling constant is proportional to
.
When
is broken by the vacuum expectation value
of
,
the effective interaction (2.88)
generates the Majorana mass term
 |
(2.90) |
where
GeV. In this mass term there is a suppression
factor
which is responsible for the smallness of the
neutrino masses [122,123].
Maximal mixing
The expression (2.64) for the mixing angle
implies that the mixing is maximal,
i.e.,
,
if
.
In this case,
assuming that
,
the Majorana neutrino masses are given by
 |
(2.91) |
The fields
and
are connected with the fields
and
by the relations
 |
(2.92) |
The massive Majorana fields are given by
![\begin{displaymath}
\nu_{1,2}
=
\frac{1}{\sqrt{2}}
\left\{
\left[ \nu_L + (\nu_L)^c \right]
\mp
\left[ \nu_R + (\nu_R)^c \right]
\right\}
\,.
\end{displaymath}](img250.gif) |
(2.93) |
If
and
,
the mass term (2.57) is simply a Dirac mass term.
Applying Eqs.(2.64) and (2.66), we obtain
 |
(2.94) |
Using Eq.(2.71) for the mixing matrix
,
one can see that
the fields
and
are connected to the Majorana fields
and
by the relations
 |
(2.95) |
where
are Majorana fields with the same mass
and with CP parities
and
(see Eq.(2.53)).
Eq.(2.93)
implies for a Dirac field
 |
(2.96) |
Thus we arrive at the well-known result that a Dirac field can always be
represented as an equal mixture of two Majorana fields with the same mass and
opposite CP parities.
One can
see this result also directly:
 |
(2.97) |
where
 |
(2.98) |
are Majorana fields.
Finally, there is
the possibility that
but at least one of the
parameters
is non-zero.
In this case Eqs.(2.95) and (2.96) are approximately valid
and
are
two Majorana neutrinos with opposite CP parities and almost degenerate
masses given by
 |
(2.99) |
The field
is called in this case a pseudo-Dirac neutrino field
[128,62,65,129,130].
Neutrino oscillations in vacuum
The general formalism
If there is neutrino mixing,
the left-handed components of the neutrino fields
(
)
are unitary linear combinations of the left-handed components of the
(Dirac or Majorana)
neutrino fields
(
) with masses
:
 |
(3.1) |
The number
of massive neutrinos
is 3 for the Dirac mass term discussed in
Section 2.1
and for the Majorana mass term discussed in
Section 2.3, in which cases there are only the three
active flavour neutrinos.
The number
of massive neutrinos is more than three in the case of a
Dirac-Majorana mass term discussed in
Section 2.2 with a
mixing of both, active and sterile neutrinos. In general, the
number of light massive neutrinos can be more than three.
We enumerate the neutrino masses in such a way that
 |
(3.2) |
In this section we will consider in detail the phenomenon of neutrino
oscillations in vacuum which is implied by
the mixing relation (3.1).
If all neutrino mass differences are small,
a state of a flavour neutrino
produced in a weak process
(as the
decay, nuclear beta-decays, etc.)
with momentum
is described by the coherent superposition of mass eigenstates
(for a discussion of the quantum mechanical problems
of neutrino oscillations
see Refs. [131,62,132,133,134,135,136,137,61,138,139,140,141,142,143,144,145,146])
 |
(3.3) |
Here
is the state of a neutrino with negative helicity, mass
and energy
 |
(3.4) |
Let us assume that at the production point and at time
the state of a neutrino is described by Eq.(3.3).
According to the Schrödinger equation the mass eigenstates
evolve in time with the phase factors
and at the time
at the detection point
we have
 |
(3.5) |
Neutrinos are detected by observing weak interaction processes.
Expanding the state (3.5)
in the basis of flavour neutrino states
,
we obtain
 |
(3.6) |
where
 |
(3.7) |
is the amplitude of
transitions at the time
at a distance
.
Consequently, the probability of this transition is given by
 |
(3.8) |
This formula has a very simple interpretation.
is the amplitude to find the neutrino mass eigenstate
with energy
in the state of the flavour neutrino
,
the factor
gives the time evolution of the
mass eigenstate and, finally, the term
gives the
amplitude to find the flavour neutrino state
in the mass eigenstate
.
We want to remark that weak interaction processes responsible for
neutrino production and detection involve active
neutrinos. Therefore, strictly speaking,
in the derivation of the Eqs.(3.7) and
(3.8) the indices
and
take only active flavour
indices. However, transitions into sterile neutrinos can be
revealed through neutral current neutrino experiments (disappearance
of active neutrinos). In this sense,
the formulas (3.7) and (3.8) have a meaning
also for transitions between active and sterile states.
Notice
that in order to have a non-negligible
active-sterile transition probability
the sterile fields must have a mixing with the
light neutrino mass eigenfields the number of which must be more than three.
Such a possibility is phenomenologically given by the Dirac-Majorana
mass term (2.21), but it is not realized
in the simple see-saw scheme discussed
in Section 2.5, where the scale of the
right-handed Majorana mass term is large.
However, the see-saw scenario can be modified by
additional assumptions to include light sterile neutrinos
(``singular see-saw'' [147],
``universal see-saw [148,149]).
At this point a remark concerning the unitarity of the mixing matrix
is at order. If some of the mass eigenstates
are so heavy that they are not produced in the standard weak processes
then these mass eigenstates will not occur in the flavour state
(3.3). Let us assume that the
first
mass eigenstates are light (
). Consequently,
only that part of
plays a role in neutrino oscillations
where
. In the following we will always assume that
in the situation described here we can confine ourselves to an
submatrix of
which is unitary to a good
approximation (see, e.g., Ref. [150]). This is realized in
the see-saw mechanism with a sufficiently large right-handed scale.
In the further discussion we will drop the distinction between
and
.
From the relation (3.1) it follows that the state describing a flavour
antineutrino
is given by
 |
(3.9) |
Thus the amplitude of
transitions is given by
 |
(3.10) |
Notice that the amplitude for antineutrino transitions
differs from the corresponding amplitude (3.7)
for neutrinos only by the exchange
.
Using the unitarity relation
 |
(3.11) |
the probability (3.8)
can be written as
![\begin{displaymath}
P_{\nu_\alpha\to\nu_\beta}
=
\left\vert
\delta_{\alpha\beta}...
...lta{m}^2_{k1} L }{ 2 E } \right) - 1 \right]
\right\vert^2
\,,
\end{displaymath}](img302.gif) |
(3.12) |
where
and the ultrarelativistic
approximation (3.4) has been used. Thus the probability
of
transitions depends on the
elements of the mixing matrix, on
independent mass-squared
differences and on the parameter
, whose range is determined
by the experimental setup.
If there is no mixing (
) or/and
for all
,
there are no transitions
(
).
Neutrino transitions can be observed only if
neutrino mixing takes place
and at least one8
satisfies the condition
 |
(3.13) |
In this inequality,
is the neutrino mass-squared in eV
,
is the distance between neutrino source and
detector in m (km)
and
is neutrino energy in MeV (GeV).
Thus, the larger the value of the parameter
,
the smaller are values of
which can be probed in an experiment.
The inequality (3.13) allows to estimate (for large mixing
angles) the sensitivity to
the parameter
of
different types of neutrino oscillation experiments.
These estimates are presented in Table 3.1.
Let us stress that this table gives only a very rough idea of the
sensitivity to
.
For example, in the LSND short-baseline (SBL) experiment
antineutrino energies are between 20 and 60 MeV, the distance is
approximately 30 m and the minimal value of
probed
in this experiment is
eV
.
| Experiment |
(m) |
(MeV) |
(eV ) |
| Reactor SBL |
 |
 |
 |
| Reactor LBL |
 |
 |
 |
| Accelerator SBL |
 |
 |
 |
| Accelerator LBL |
 |
 |
 |
| Atmospheric |
 |
 |
 |
| Solar |
 |
 |
 |
The probability (3.8) and the corresponding one for antineutrinos
are invariant under the phase transformation
 |
(3.14) |
Therefore, it is clear that the probabilities
of
and
transitions do not depend
[28,29,106]
on the Majorana CP-violating phases
discussed at the end of Section 2.3 and it is not
possible to distinguish the Dirac and Majorana cases by the
observation of neutrino oscillations.
Comparing the expressions
(3.7) and (3.10)
for the transition amplitudes
of neutrinos and antineutrinos we see that
.
Therefore,
for the transition probabilities we have
 |
(3.15) |
This relation is a consequence of CPT invariance
inherent in any local field theory [151].
From the equality (3.15)
it follows that the neutrino and antineutrino survival
probabilities are equal:
 |
(3.16) |
On the other hand,
the transition probabilities of neutrinos and antineutrinos in general are
different for
.
They are equal only if there is CP invariance in the lepton sector.
In fact,
in the case of massive Dirac neutrinos the phases of
the neutrino fields and of the charged lepton fields
can be chosen in such a way that
the mixing matrix
is real.
In the case of Majorana neutrinos,
CP invariance implies that
(see Eq.(2.50))
 |
(3.17) |
with
,
where
is the CP parity of the Majorana neutrinos with mass
(see Eq.(2.53)).
Thus,
in the Dirac as well as in the Majorana case,
CP invariance implies that
 |
(3.18) |
The results of neutrino oscillation experiments are usually analysed
under the simplest assumption of oscillations between two neutrino
types. In this case, for the transition
probability (3.12) we get
![\begin{displaymath}
P_{\nu_\alpha\to\nu_\beta}
=
\left\vert
\delta_{\alpha\beta}...
...{ \Delta{m}^2 L }{ 2 E } \right) - 1 \right]
\right\vert^2
\,,
\end{displaymath}](img335.gif) |
(3.19) |
where
and
,
are
,
,
or
,
,
etc
Thus in the simplest case of
transitions between two neutrino types the probability is determined
only by the elements of
which connect flavour neutrinos with
(or
). It is obvious that phases drop out in the expression
(3.19). This is an illustration of Eq.(3.14)
and of the fact that no information about CP violation can be obtained
in the case of transitions between only two neutrino types.
If we put
then we have
and the transition and survival
probabilities are given by the following standard expressions
valid for both, neutrinos and antineutrinos:
| |
|
![$\displaystyle P_{\stackrel{\makebox[0pt][l]
{$\hskip-3pt\scriptscriptstyle(-)$}...
...ft( 1 - \cos \frac{ \Delta{m}^2 L }{ 2 E } \right)
\qquad
(\alpha\neq\beta)
\,,$](img343.gif) |
(3.20) |
| |
|
![$\displaystyle P_{\stackrel{\makebox[0pt][l]
{$\hskip-3pt\scriptscriptstyle(-)$}...
...\scriptscriptstyle(-)$}}{\nu_{\beta}}}
\,.\setlength{\arraycolsep}{\templength}$](img344.gif) |
(3.21) |
The transition probability (3.20)
can be written in the form
![\begin{displaymath}
P_{\stackrel{\makebox[0pt][l]
{$\hskip-3pt\scriptscriptstyle...
...rac{ \Delta{m}^2 L }{ E } \right)
\qquad
(\alpha\neq\beta)
\,.
\end{displaymath}](img345.gif) |
(3.22) |
where
is the source - detector distance expressed in m (km),
is the neutrino energy in MeV (GeV) and
is the neutrino mass-squared difference in eV
.
Thus, the transition probability is a periodic function of
. This
phenomenon is called neutrino oscillations.
The amplitude of oscillations is
and the oscillation length is given by
 |
(3.23) |
The condition (3.13) can be rewritten in the form
 |
(3.24) |
Therefore,
neutrino oscillations can be observed if the oscillation length
is not much larger than the source - detector distance
.
The oscillatory behaviour of the transition probability (3.20)
with
is shown
in Fig. 3.1,
where we have plotted it as a function of
.
The grey line represents the transition probability (3.20),
whereas, in order to demonstrate the effect of energy averaging,
the black line represents the transition probability (3.20)
averaged over a Gaussian energy distribution with mean value
and
standard deviation
.
The averaged probability is
the measurable quantity in neutrino oscillation experiments.
One can see that the averaging over the energy spectrum practically
reduces the probability to the constant
for
.
The expressions (3.20) and (3.21)
are usually employed in analyses of the data of
neutrino oscillation experiments.
In many SBL experiments with neutrinos from reactors
and accelerators, no indication in favour of neutrino oscillations was found.
The data of these experiments give an upper bound
for the transition probability which implies
an excluded region in the space of the parameters
and
.
A typical exclusion plot is presented in Fig. 3.2 [152].
This plot shows the exclusion curves in the
channel
obtained in the
CDHS [153],
FNAL E531 [154],
CHARM II [155],
CCFR [156],
CHORUS [157]
and
NOMAD [158]
experiments.
The excluded region lies on the right of the curves.
The two most stringent exclusion curves
in Fig. 3.2 have been obtained in the
CHORUS [157]
and
NOMAD [158]
experiments, which are operating at CERN using the
neutrino beam from the SPS (with an average energy of about 30 GeV).
800 kg of emulsions are used in the CHORUS experiment as target.
The production and decay of
's
in the emulsion is searched for.
In the NOMAD experiment a magnetic detector is used and
the production of
's is identified with kinematical criteria.
Figure 3.3
shows the exclusion curves obtained in the
CHOOZ [159],
Gösgen [160],
Krasnoyarsk [161]
and
Bugey [162]
reactor
experiments.
The region allowed by the results of the Kamiokande
atmospheric neutrino experiment [49]
(see Section 5.1)
is also depicted in this figure.
The Gösgen, Krasnoyarsk and Bugey experiments
are SBL reactor experiments,
whereas
the recent CHOOZ experiment is the first LBL reactor experiment.
In this experiment the detector
(5 tons of liquid scintillator loaded with Gd)
is at the distance of about 1 km from the CHOOZ power station, which has
two water reactors with a total thermal power of 8.5 GW.
The antineutrinos are detected through the observation of
the reaction
 |
(3.25) |
(the photons from annihilation of the positron and
the delayed photons
from the capture of neutron by Gd are detected).
No indications in favour of neutrino oscillations were found in the CHOOZ
experiment. The ratio
of the numbers of measured antineutrino
events and of expected antineutrino
events in the case of absence of neutrino oscillations
in the CHOOZ detector is [159]
 |
(3.26) |
Since the average value of
in the CHOOZ experiment is
approximately 300
(
,
),
in this
experiment it is possible to probe the value of
the relevant
neutrino mass-squared difference down to
.
With the help of the expression (3.22),
it is possible to understand
qualitatively the general features of exclusion curves.
In the region of large
such that
the oscillation length is much smaller than the
source - detector distance
,
the cosine in the expression (3.22)
oscillates very rapidly as a function of the neutrino energy
.
Since in practice all neutrino beams have an energy spectrum
and
the neutrino sources and detectors are extended in space,
only the average transition probability
![\begin{displaymath}
\langle
P_{\stackrel{\makebox[0pt][l]
{$\hskip-3pt\scriptscr...
...(-)$}}{\nu_{\beta}}}
\rangle
=
\frac{1}{2} \, \sin^22\vartheta
\end{displaymath}](img366.gif) |
(3.27) |
can be determined in the region of large
.
The average probability is independent from
and, therefore, from
an experimental upper bound
on
one obtains the vertical-line part of the exclusion curve.
At
 |
(3.28) |
where
is the average energy and
is the average distance,
the parameter
has the minimal value
![\begin{displaymath}
(\sin^22\vartheta)_{\mathrm{min}} =
\langle
P_{\stackrel{\m...
...]
{$\hskip-3pt\scriptscriptstyle(-)$}}{\nu_{\beta}}}
\rangle_0
\end{displaymath}](img372.gif) |
(3.29) |
on the boundary curve.
Note that we are using here and in the rest of the section the same
units as in Eq.(3.22).
Typically, the upper bound
is much less than one. Then, in the region where
is
large
the expression (3.22) for the transition probability
can be approximated by
![\begin{displaymath}
P_{\stackrel{\makebox[0pt][l]
{$\hskip-3pt\scriptscriptstyle...
...2\vartheta
\left( 1.27\, \frac{ \Delta{m}^2 L }{ E } \right)^2
\end{displaymath}](img373.gif) |
(3.30) |
and in this region the boundary curve in the exclusion plot is given by
 |
(3.31) |
Therefore, this part of the exclusion curve
is a straight line in the
-
plot
as can be seen from Fig. 3.2.
From Eq.(3.31) it follows that the minimal
value of the parameter
that can be probed by
an experiment
is
 |
(3.32) |
It corresponds to
.
Neutrino oscillations and transitions in matter
It has been pointed out by Wolfenstein [163] and by
Mikheyev and Smirnov [164] that the neutrino
oscillation pattern in vacuum can get significantly modified by the
passage of neutrinos through matter because of the effect of coherent
forward scattering. This effect can be described by an effective
Hamiltonian. Starting with neutrino oscillations in vacuum,
one can easily check that the transition
probability (3.12) can be obtained by considering the evolution of
the state vector
of the neutrino types with the
``Schrödinger equation''
 |
(4.1) |
where
is the effective Hamiltonian for
neutrino oscillations in vacuum
with
being the mixing matrix,
the neutrino energy and
the diagonal neutrino mass matrix. The presence of matter will give a
correction to Eq.(4.1). Note that the variable in
Eq.(4.1) is not time but space and the components of
are
given by the amplitudes denoted by
for the neutrino
types
. The
corresponding effective Hamiltonian for antineutrinos
is obtained from
by the exchange
. After the seminal work of Wolfenstein
about neutrino oscillations in matter (see also Ref. [165])
and the discovery of the
importance of
,
the analogue of
in matter,
for the solar neutrino problem
[164] many rederivations of the effective matter Hamiltonian
[166,167] have been presented using a coupled
system of Dirac equations [168,169], field theory [170]
or studying coherent forward scattering in more detail [171].
We want to sketch a derivation using the Dirac equation and
following Ref. [169]. The starting point in most
derivations is the expectation value of the currents for
isotropic non-relativistic
matter given by [163,164,172]
 |
(4.2) |
where
is the number density of the particles represented
by the field
. Note that the
term does not
contribute for non-aligned spins [173].
Eq.(4.2) is a good approximation even for electrons in the
core of the sun with a temperature
K
[174,175]
and thus an average velocity of the electrons of
around 10% of the velocity of light and corrections of order
.
Starting from the weak interaction Lagrangians
(2.1) and (2.2) one gets for low-energy neutrino
interactions of flavour
with the background matter
 |
(4.3) |
where
is the Fermi coupling constant,
the weak
mixing angle,
the eigenvalue of the field
of the third
component of the weak isospin (
in the Standard
Model) and
is the charge of
. In the matter Lagrangian (4.3)
the charged current interaction is represented by the
Kronecker symbol
saying that for neutrinos
of flavour
the charged current only contributes when
background matter containing charged leptons of the same flavour
is present. Concentrating now on realistic matter with
electrons, protons and neutrons which is electrically neutral,
i.e.,
, and taking into account that we
have
and
,
for electrons, protons and
neutrons, respectively,9we get an effective Hamiltonian
 |
(4.4) |
acting on the vector
of the flavour neutrino wave functions,
where we have defined
,
and
.
is
the non-diagonal mass matrix. The last
term in Eq.(4.4) is called the matter potential term.
Several remarks concerning the Hamiltonian (4.4), which is
the point of departure for deriving the effective matter
Hamiltonian in Refs. [168,169], are at order:
has spinor and flavour indices, it leads thus to a
system of Dirac equations coupled via the mass term.
- The neutral current contributions of electrons and protons
cancel for realistic matter
because of opposite
quantum numbers and electric
charges (see discussion above before Eq.(4.4)).
- The Hamiltonian (4.4) is valid for Dirac neutrinos. For
antineutrinos,
is replaced by
and
the matter potential term of
has the opposite sign and the projector
instead of
. This is a
consequence of proceeding along the same lines as before but
using the charge-conjugate fields to describe antineutrinos
(
).
- In the Majorana case, in order to obtain the Hamiltonian
the antineutrino matter potential has
to be added to
and the neutrino field vector
is
subject to the Majorana condition
.
This condition is compatible with the time evolution governed by
[169].
- Sterile neutrinos can easily be incorporated into
and
by simply adding the entry 0 to the matter potential term. Note that
the probability of active - active transitions does not depend on the
neutron density
, which is, however, important for active -
sterile transitions (see Table 4.1).
The essence of deriving an effective Hamiltonian for neutrino
oscillations in matter is to get rid of the spinor indices in
and to obtain an equation involving only the indices
for the different
neutrino types as in the
vacuum case (4.1). To this end,
let us now assume that the neutrino propagation proceeds along
the
axis, that the neutrino momentum
corresponding to propagation in vacuum is much larger than the
matter potentials and that the neutrinos are ultrarelativistic.
Thus we consider a one-dimensional problem
from now on. Defining a wave function
via
 |
(4.5) |
this wave function changes little over distances of the order of
the de Broglie wave length of the neutrino. With respect to
we can decompose the Hamiltonian (4.4) into
which are the
parts commuting or anticommuting with
.
consists solely of the mass term and
contains the rest of
. With
being a large parameter we can perform a Foldy - Wouthuysen
transformation where we truncate the series at
leading to
the Hamiltonian [176,169]
such that all terms in
commute with
and, therefore, all matrices in
with spinor indices
can be diagonalized at the same time in order to separate
positive and negative energy states. Denoting this unitary
diagonalization matrix by
we can achieve
,
and
.
Note that there are no
transitions between left and right-handed states through
. For neutrinos with magnetic moments
(electric dipole moments) the method described above would also
yield the appropriate left-right transitions [169].
Going back to the Foldy - Wouthuysen transform of the wave
function
instead of
removes the
from the
Hamiltonian (4.6). The final step in deriving the effective
matter Hamiltonian consists of considering stationary states and
splitting off the plane wave part by
 |
(4.7) |
Then the wave function
, taken for positive energies
and left-handed neutrinos, fulfills the first order differential
equation (4.1) with
replaced by [163]
 |
(4.8) |
correct up to order
. This effective Hamiltonian is valid
for left-handed Dirac neutrinos or Majorana neutrinos, whereas
in the case of right-handed (Dirac) antineutrinos or
right-handed Majorana neutrinos we have
 |
(4.9) |
because for (Dirac) antineutrinos one has to replace
by
in the Hamiltonian (4.6) and for Majorana neutrinos
holds (see also points 3 and 4 after
Eq.(4.4)). With
 |
(4.10) |
only the left-handed mixing matrix
enters and connection
between Eqs.(4.8) and (4.1) is made.
Thus in the final results (4.8) and (4.9)
there is no distinction between the
Dirac and Majorana case for ultrarelativistic neutrinos.
This is an illustration of the fact that with neutrino
oscillation experiments the Dirac or Majorana nature cannot be
distinguished [98,177].
The Hamiltonians (4.8) and (4.9) have been used to
investigate neutrino oscillations in the sun, in the earth and
in supernovae. In the following we will only be concerned
with the first two subjects. Application limits of the neutrino
evolution equations in matter have been discussed in
Ref. [178]. Elastic and inelastic neutrino scattering
introduces quantum damping into the evolution equations which
is proportional to the neutrino interaction rate [179]. In
the sun and the earth this effect is negligible, in particular, for
low energies, whereas in the early universe it is of crucial
importance [179]. Density fluctuations have been found to
influence considerably neutrino propagation in the sun
[180], however, more realistic considerations with
helioseismic waves as density fluctuations show no observable effect
on the solar neutrino problem discussed in terms of neutrino
oscillations [181]. Solar neutrinos are also influenced by
their passage through the earth [182].
Let us now concentrate on left-handed neutrinos and specify the effective
Hamiltonian (4.8) to two neutrino types. Thus for
two-neutrino oscillations in matter with the definitions
 |
(4.11) |
and
 |
(4.12) |
the differential equation
has to be studied [163,164,178,183,172,184],
where
,
are the amplitudes for the neutrino types
,
(
), respectively.
The evolution equation for antineutrinos is obtained from
Eq.(4.13) by
. In the two-neutrino case
there is one
and the mixing matrix
 |
(4.14) |
is real without loss of generality.
The list of all possible matter densities which determine
and
occur in the
different oscillation channels is given in Table 4.1.
Evidently,
oscillations proceed
as in vacuum because
and the term proportional to the unit
matrix in Eq.(4.13) has no effect on
transitions.
As mentioned before, in the rest of this section we will have in mind that
Eq.(4.13) is applied to neutrino propagation in the sun
and the earth.
Let us first define the eigenfunctions of the effective
Hamiltonian as
 |
(4.15) |
where the matter angle
and the energy eigenvalues
are functions of
. In
the limit of vanishing matter density
is identical with
the vacuum angle
(4.14). The eigenvalues
and the matter angle are given by
and
 |
(4.17) |
respectively.
The general solution of the differential equation
(4.13) can be represented as
 |
(4.18) |
with the coefficients
fulfilling
 |
(4.19) |
and
.
The adiabatic solution is defined as an approximate
solution where
(
).10Postponing
the discussion of the question under which
condition adiabaticity is fulfilled, we consider temporarily the
case of an arbitrary number
of neutrino flavours or types and
define the mixing matrix in matter
via
 |
(4.20) |
where the eigenvectors
generalize Eq.(4.15). One
readily obtains the generalization of the vacuum oscillation
amplitude Eq.(3.12) for neutrino oscillations in matter if the
evolution of the neutrino wave function is adiabatic:
 |
(4.21) |
In this formula neutrino production and detection happen at
and
, respectively. The adiabatic phases
,
which are defined by
 |
(4.22) |
are necessary for the correct evolution of
in the
adiabatic limit. Usually, they can be absorbed into the eigenvectors
[185], except for special matter densities, where these phases
acquire a topological meaning, together with the presence of CP
violation [186]. Note that for real vectors
one has
and therefore
.11
This is so in the two-dimensional case Eq.(4.15).
Evaluating Eq.(4.21) for two neutrino types and assuming that an
averaging over neutrino energies takes place such that
 |
(4.23) |
it is easy to show that the averaged two-neutrino survival probability can be
written in the adiabatic case as
 |
(4.24) |
Let us now derive a condition for adiabaticity in the two-neutrino case.
The evolution of the neutrino state in matter is adiabatic if
the right-hand side of Eq.(4.19) can be neglected.
A formal solution of Eq.(4.19) is given by
 |
(4.25) |
where
denotes the matrix on the right-hand side of
Eq.(4.19). Defining the variable
 |
(4.26) |
and the adiabaticity parameter (as a function of
)
[187,188,189,184]
 |
(4.27) |
one can determine an upper bound to the typical integral which occurs in all
the terms of the sum (4.25) on the right end. In an interval
where
is monotonous there is a value
such that [190]
 |
(4.28) |
according to the second mean value theorem of integral calculus
(
). This consideration allows to put the bound
![\begin{displaymath}
\left\vert \int_{y_0}^{y_1} dy\, \dot{\vartheta}_m e^{2i\alp...
...gamma_\mathrm{min} \equiv \min_{y \in [y_0,y_1]} \gamma(y) \,.
\end{displaymath}](img509.gif) |
(4.29) |
Having found this bound, in the rest of the integrals in the terms in
Eq.(4.25) one can simply take the absolute values of the
integrands and one is left with integrations of the type
. With
the refinement that the interval
possibly has to be divided into
several parts labelled by
such that in each part interval
is monotonous we can define
 |
(4.30) |
and obtain the exact bound
 |
(4.31) |
where the symbol
denotes the
-th element of the vector on the
right-hand side of this inequality and
has a
contribution
for every interval where
has a definite sign. (
is the
difference between the angle
in the initial and the
final point of such an interval.)
Thus with the exact bound (4.31) we have found the following
sufficient condition:
 |
(4.32) |
It is fulfilled if on the scale of the oscillation length in matter
given by
the matter angle
changes
very little, i.e.,
.
Eq.(4.31) allows to get an upper bound on the crossing or
jumping probability from
to
. Assuming the
initial conditions
and
we get
 |
(4.33) |
which illustrates once more the importance of
for
adiabaticity.
In the context of the solar neutrino problem the discovery that
a resonance in the passage of
through
the sun is possible [164] gave a major boost to the
investigation of the propagation of neutrinos in matter. The
possibility of a resonance is most easily understood in the
adiabatic approximation by looking at Eqs.(4.17) and
(4.24). For different neutrino masses we can always label them
in such a way that
. If on the way from the creation
point
in matter of
to a
point
in vacuum the neutrino passes
through a point
where the resonance condition
 |
(4.34) |
is fulfilled then
has to be between
and
(for
the quantity
is positive!) whereas
according to Eq.(4.17) the matter angle
is found between
and
.12Consequently, the product of cosines in the
survival
probability (4.24) is negative and the probability to
find a
after the passage through the sun is less than
1/2. The interesting phenomenon is that this can happen even for
relatively small mixing angles, provided the matter potential at
the production point is large enough such that there is a point
where the condition (4.34) is fulfilled.
Indeed, for
close to
, one can even have
close to
and thus
(4.24).
It has been shown that the resonance is also effective in a
certain area in the
-
plane
where the evolution of the neutrino state is non-adiabatic. For
antineutrinos a resonance is possible if
.
If there is a resonance, then for reasonable matter densities the
adiabaticity parameter (4.27) is smallest at the resonance
and thus adiabaticity is most likely violated there. Therefore,
considering
at the resonance,
from Eq.(4.17) one easily computes
. A
suitable measure of adiabaticity in the case of a resonance
is thus given by [187,188]
 |
(4.35) |
In the solar interior the electron density is maximal in
the center with
where
and
[174].
This leads to
 |
(4.36) |
Considering the resonance condition (4.34) one can read
off from this equation for which
a resonance is possible.
Whereas the density variation of
and
in the sun is smooth,
in the earth due to its layered structure (lithosphere - mantle -
core) the densities can be approximated as step functions
[191,192]. Furthermore, approximately
is valid in all layers. In the lithosphere, which is relevant for LBL neutrino
oscillation experiments, the average mass density
is around 3 g/cm
. It increases to around 13 g/cm
in the core.
Thus we can conveniently write
 |
(4.37) |
The analogue of the adiabatic formula (4.21) in terms of
probabilities is given by
 |
(4.38) |
with
where
is the diagonal matrix of energy eigenvalues present in
Eq.(4.20) (see Ref. [183]).
is a unitary
matrix which is
diagonal in the adiabatic limit corresponding to
(see Eq.(4.22)).
Confining ourselves now to two neutrino types
and having in mind neutrino production in matter and detection in
vacuum we use the notation
and
. With the crucial assumption that
averaging over neutrino energies and the neutrino production region
all terms other than probabilities can be dropped we get
 |
(4.40) |
where
is the crossing probability from
to
and in the second part of Eq.(4.40) unitarity has been
used. Such averaging procedures have been shown to be
effective in the context of solar neutrinos, e.g., in
Refs. [193,194]. Inserting the relations (4.40) into
Eq.(4.38) specialized to the survival probability of
one obtains after some algebra with trigonometric functions the
Parke formula
 |
(4.41) |
which is the
generalization of Eq.(4.24) [189,183].
This formula was first derived in the context of an exact solution of
Eq.(4.13) for matter densities linear in
[188,195] (see also Refs. [196,187]).
The probability
can be estimated with the Landau - Zener method
[197,196]. Following the derivation of Landau, the idea
is to make an analytic
continuation of the matter densities and therefore of the effective
matter Hamiltonian (4.8) into the complex plane by
where
is a complex variable. Then also the
neutrino state
(4.18) is analytic. Considering the
energy eigenvalues
(4.16) as functions of
we find two branching points
and
(
) defined by the equation
 |
(4.42) |
We consider
and
for definiteness and follow its evolution along a
curve
which starts at a value
in the vicinity of the real
axis and goes along a negatively oriented loop around
such that we end
at the same point
but now located on the second sheet.13Along this path
changes into
[197] and
into
.14 The crucial observation
is that, though
and
(
) are two-valued
analytic objects, the analytic continuation
of the solution
(4.18) of the differential equation (4.13) is
single-valued [198]. From this fact it can easily be shown
with Eq.(4.18) that [198]
 |
(4.43) |
where
is the analytic continuation of
in the
vicinity of the real axis and then along
. The Landau - Zener
crossing probability is derived from the following
consideration. Starting on the real axis at
,
where the evolution of
is adiabatic and where the boundary
conditions
,
hold, and going along a path in
the complex plane which passes above
and returning to the real
axis at a point
on the second sheet
this path can be chosen such
that the evolution in the complex variable
is adiabatic because we
have passed the resonance in safe distance and we have thus
. On the other hand, we can go from
on the first sheet to the
point
on the second sheet via
and
with Eq.(4.43) we thus obtain
 |
(4.44) |
neglecting deviations of
from 1.
From this formula we can read off that integrations along the real axis do
not contribute to
. Because of analyticity, the path
can be deformed such that it goes along the real axis from
to
from where it leads to
, circles around the branching point with infinitely small radius
and then goes back to
and
on the second
sheet. In this way, we obtain from Eq.(4.44)
the final result for the Landau - Zener crossing probability
![\begin{displaymath}
\ln P_c = -\frac{1}{E}\, \mbox{Im} \int_{x_\mathrm{res}}^{z_...
...cos 2\vartheta)^2 + (\Delta m^2 \sin 2\vartheta)^2 ]^{1/2}
\,.
\end{displaymath}](img594.gif) |
(4.45) |
To get this formula we have used the explicit form of the energies of
the adiabatic states (4.16).
With the variable
transformation
and
it is easy to
evaluate Eq.(4.45) for a linear density
[187,188,189,172] and is also possible for an exponential
density [199]. These calculations can be summarized by
[189,183]
 |
(4.46) |
with
for the linear and
for the
exponential case. For an exponentially varying matter density
the dependence of
on
is particularly simple because
with
which is independent of the vacuum mixing angle.
For a more general discussion of the crossing probability (4.45)
see Ref. [200].
Exact solutions of the differential equation (4.13) for
neutrino oscillations in matter exist not only for the linear case
[196,187,188,195] in terms of Weber functions but also for the
exponentially varying matter density in terms of Whittaker functions
[201,193]. This case is of particular importance
for solar neutrinos because it approximates the real density variation
in the sun. Further exact solutions are known for
varying with
[202] and
[189].
In Refs. [189,183]
a list of the factors
for all these cases is given, calculated with
the Landau - Zener formula (4.45).
The survival probability (4.41) with
in the Landau -
Zener approximation does not reproduce well the exact survival
probability in the extremely non-adiabatic region [195,193]
(see also Refs. [189,183,172]). In these references the
example of an extremely dense medium (
and thus
) with a sharp boundary to the vacuum is given. In this case
the neutrino does not oscillate in the medium and therefore
stems purely
from vacuum oscillations. This has to be compared with
,
because of the jump in density one gets
, inserted into
Eq.(4.41). Obviously, the resulting expression
does not agree
with the previous one.
A remedy of
this deficiency was found in the framework of the exact solution for
an exponentially varying matter density leading to the following
modification [193] of the Landau - Zener crossing probability
(4.45):
 |
(4.47) |
Numerical calculations for solar neutrinos with an exponential
electron density using this formula agree
very well with numerical integrations of the differential equation
(4.13), typically the agreement is in the percent range for
relevant mixing parameters [203]. It has been conjectured
in Ref. [189] that the form (4.47) of
holds for a wide
class of matter density profiles. In any case, in the limit
one gets
which,
when inserted into formula (4.41),
correctly describes the survival probability in the above example of
an extremely non-adiabatic evolution.
Concluding this section, we consider again an arbitrary number of neutrino
flavours or types and envisage the interesting case
where one of the eigenfunctions
(4.20) of
labelled by the index
has an adiabatic evolution whereas the other part of
has an
arbitrary evolution. We assume
exact adiabaticity for
which amounts to
for
. This
allows to write
 |
(4.48) |
where we have defined
That
is smaller than 1
follows from the Cauchy - Schwarz inequality and from the
above assumption for
because after dropping the row and
the column labelled
by
in
the remaining matrix is again unitary.
Eq.(4.49) and similar cases are useful if neutrino masses
differing by orders of magnitude occur.
Indications of neutrino oscillations
Atmospheric neutrino experiments
In 1912 it was discovered by V.F. Hess [204] in a manned balloon flight
that the intensity of the ionizing radiation in the atmosphere
as a function of the altitude did not conform with the idea
that this ionization was caused by radioactive elements in the surface
of the earth but rather
pointed to an extraterrestrial origin. In the following decades,
before the advent of accelerator physics, this
radiation, which was called first ``ultraradiation'' and later
baptized ``cosmic rays'' by R.A. Millikan, proved to be one of the most
fruitful means for doing particle physics experiments. At the end of
the first half of the 20th century such experiments had lead to the
discovery of the positron, the pion and the muon
and also the first particles with strangeness were found with
cosmic rays [205]. Eventually,
in the beginning of the fifties, proton beams from accelerators
replaced the cosmic proton flux as an experimental tool.
However, after many years
where accelerator physics was dominating in particle physics, at
the end of the 20th century cosmic rays play again a major role
through atmospheric neutrinos which allow to use the whole
globe as a neutrino physics laboratory and to probe neutrino
mass-squared differences down to a few
eV
. In this way
convincing evidence for the existence of neutrino oscillations and thus
for non-zero neutrino masses has been obtained [2].
In a simplified picture, the production of atmospheric neutrinos
[206,207,208,209,210,211,212,213,214,215,216,217,218]
proceeds in three steps [219].
In the first step the primary cosmic rays [219,220]
hit the nuclei in the
atmosphere, thereby producing charged pions and kaons, either
directly or via intermediate particles. In the second
step, the decay of these particles gives rise to part of the
atmospheric
and
neutrino fluxes:
 |
(5.1) |
In the third step, the
and
fluxes and further
and
fluxes are produced by
 |
(5.2) |
There is also a contribution to the neutrino fluxes from the decays
 |
(5.3) |
with
and the charge conjugate processes
[219,221,215] which does not conform with the simple 3-step
picture. At low energies the most important process in
Eq.(5.1) is the pion decay. The contribution of
decays (5.3) to the neutrino fluxes is small
[222,223]. Let us now describe
some useful details of the production of the atmospheric neutrino flux
and its properties.
Galactic cosmic rays enter the solar system as an isotropic flux of
particles. There is convincing evidence that the bulk of the radiation with
energies less than
GeV comes from our galaxy. Though there is
no conclusive proof of the origin of cosmic rays yet, plausible
mechanisms range from material ejected by supernovae to interstellar
medium accelerated in supernova shock waves [224].
Solar cosmic rays are emitted irregularly by major solar flares on the sun.
Apart from electrons, primary cosmic rays consist of protons and bare
nuclei (mostly He). Roughly speaking,
the chemical composition of the galactic cosmic ray particles
is given approximately by 90 % H, 9 % He and less than 1 % heavier
nuclei. However, the chemical composition varies with energy. At
energies of around 100 MeV per nucleon the particle number ratio H/He
is less than 5, it increases to 10 at 1 GeV and is around 30 at 100
GeV [225,226].
The absolute flux of cosmic protons is not very large: it is of the order of
for energies of a few GeV and above the atmosphere. As a function of
the energy the proton flux above a few GeV is well described
by the power law
until
GeV [225].
Above an energy of 100 GeV per nucleon the cosmic ray
fluxes are less precisely known [225,226].
The solar cosmic rays (solar wind)
are important in so far as they weaken and modulate the flux of galactic
cosmic ray particles with the solar activity. The stronger the solar wind
the more difficult it is for the low energy galactic cosmic rays to
enter the solar sphere of influence. This modulation is noticeable for
kinetic energies of around 10 GeV per nucleon or less [224,215].
This variation of the primary cosmic ray flux with the approximate
11-year cycle of solar activity can be
parameterized as a function of the neutron monitor at Mt. Washington
[215]. It induces a corresponding modulation of the
low energy atmospheric neutrino flux but its influence is small. The effect
of the solar wind on cosmic ray particles
depends on the rigidity of the nuclei which is defined by
momentum/charge. It is strongest for
cosmic ray particles with small rigidity. For
neutrino energies
GeV it is practically negligible, for
energies around 1 GeV only appreciable at high geomagnetic latitudes
because the geomagnetic cut-off admits cosmic ray particles with
smaller rigidity there.
The geomagnetic field prevents primary cosmic particles with low
rigidity from entering the atmosphere
[227,215,228,218]. In the atmospheric
neutrino flux calculations it is assumed that the momenta of
the neutrinos have the same direction as the primary cosmic ray
particles responsible for their production [215,217]. This is
a good approximation for
MeV as shown in
Ref. [229]. In a first
approximation, for a given direction, there is a cut-off in
rigidity for the primary cosmic ray particles below which such a
particle cannot reach the top of the atmosphere in a vertical
altitude of around 20 km. The value of the geomagnetic cut-off
can be obtained by a computer simulation with the help of the
``backtracking technique''. In this technique, to establish if a
particle with a given charge and momentum coming from a
certain direction can reach a point on top of the atmosphere, one
integrates the equations of motion for a particle with opposite
charge and reflected momentum starting from its final position.
If the backtracked particle reaches infinity it is assumed that
the trajectory is allowed. As infinity one can take a distance
of around 10 times the radius of the earth where the geomagnetic
field has decreased to the level of the interstellar magnetic
field with a strength
Tesla
[215]. More refined calculations replace the rigidity
cut-off by a probability distribution [227,228,218]. In
the rigidity cut-off approximation one can give for every
experimental location a contour map depicting curves with
constant cut-off in a plot with the coordinate axes representing
azimuth and zenith angles of the primary cosmic ray particle
direction. For Kamioka such a plot is found in Ref. [215].
At the geomagnetic latitude of
the
geomagnetic cut-off disappears, therefore experiments situated at
high geomagnetic latitudes have a higher flux of low energy
atmospheric neutrinos. The magnetic field of the earth is not
exactly of dipole form but also has a considerable multipole
contribution which affects the local rigidity cut-off maps and
thus the atmospheric neutrino flux in particular for down-going
neutrinos [215].
After entering the atmosphere the cosmic ray particles collide
with the nuclei in the air.
Though the fractions of other particles than protons in the cosmic
rays are small, their importance for the atmospheric neutrino
flux gets
enhanced because the predecessors of the neutrinos, pions and kaons,
are created through the hadronic interactions of the cosmic rays with
the air nuclei. These processes depend rather on the number of
nucleons than on the number of nuclei. The ratio of the
charged kaon to the charged pion flux depends sensitively on the
description of hadronic interactions. In this context the quantities
 |
(5.4) |
where
and
comes from the
power law of the primary cosmic protons,
is the fraction of the proton
momentum carried by the meson and
is the distribution of the
charged mesons produced by collisions of protons with nuclei in the
atmosphere, are important indicators for the fraction of the muon neutrino flux
originating from the
meson flux.
The ratio
ranges from 0.10 to
0.15 in different calculations [217]. As a consequence, at
neutrino energies below 100 GeV the pions dominate as neutrino sources
[221,216].
Some of the characteristic properties of the atmospheric
neutrino fluxes are simple consequences of the production
mechanisms (5.1) and (5.2), whereas others follow from
a Monte Carlo or analytic calculation [221] of the neutrino fluxes.
As mentioned above, the flux calculations are
one-dimensional and energy losses in the air, in particular, for the
muons, have to be taken into account [219].
- a..
- Approximately, the
fluxes have a power law
energy spectrum
for
GeV
whereas the
fluxes decrease like
[215].
- b..
- Fixing a neutrino energy
, the largest contribution to
the flux of this energy originates from cosmic ray particles
with energy
[215,216].
- c..
- Denoting neutrino fluxes by
then it follows
immediately from Eqs.(5.1) and (5.2) that
 |
(5.5) |
However, it can easily be estimated that, at neutrino energies
larger than 1 GeV, muons from the reactions (5.1)
start to reach the surface of the earth before they decay.
Therefore the neutrino fluxes from muon decay decrease, in particular,
the
flux,
and the ratio (5.5) begins to rise [219,215].
- d..
- For the same reason muons cease to be the dominant source of
above
GeV and the
reactions (5.3) take over [219,221].
- e..
- Pions cease to be the dominant source of the
flux for neutrino
energies above 100 GeV [221]. The reason is given by the different
masses and lifetimes of the parent mesons
in the decays
.
For relativistic parent mesons and for a given energy
,
the lower bound on the energy of the parent meson in such a
decay is given by
.
Thus, neutrinos with energy
from the decay of the heavier
parent meson need in
average a lower parent meson energy where the corresponding meson flux
is higher. This together with the shorter kaon lifetime leads to
the prevailing of kaons as the source of muon neutrinos at
sufficiently high neutrino energies [230,221,215,216].
- f..
- The ratio
is smaller than 1
because among the primary cosmic ray nucleons the protons by far
dominate which
leads to an enhancement of positively charged pions over
negative ones in the hadronic interactions and,
consequently,15
[221,215].
- g..
- On the other hand, for
GeV, where the muons
have time to decay in the atmosphere, one has
because the muons supply the
antineutrinos (neutrinos) to the neutrinos (antineutrinos) from pions
and kaons. However, for neutrino energies above a few GeV, when part
of the muons reaches the surface of the earth before they decay, this
ratio becomes smaller than one because the
are more numerous
than the
(see above).
- h..
- For
GeV the geomagnetic effects become negligible.
The largest error in the atmospheric neutrino flux calculations arises
from an overall uncertainty of primary cosmic ray measurements of
the order of
[215,217]. Therefore, it is
mandatory to use quantities which are ratios such that this
uncertainty cancels or, in the case of a fit to the data, to include
the overall normalization of the primary cosmic ray flux in the
set of parameters to be determined by the fit. Among the flux calculations of
the different groups the major source of differences comes from
the treatment of pion production in collisions of protons with light nuclei
[217,231]. Thus, the calculated absolute fluxes have
uncertainties of the order of
and even larger
uncertainties at
neutrino energies below 1 GeV, whereas flux ratios
differ in general only by a few percent (see the comparison
between different flux calculations in Refs. [217,231]).
The direction of the atmospheric neutrino flux at the location
of an experiment is described by the zenith angle
and
the azimuth angle
. Neutrinos going vertically downward
have
whereas those coming vertically upward through the earth
have
. The azimuth angle indicates the direction
of the flux in the horizontal plane. There are two natural causes for an
angular dependence [218] of the atmospheric neutrino flux:
- a..
- The development of cosmic ray showers in the atmosphere
depends on the density of the air. At zenith angles
, where along the line of flight the increase
in density is less steep, pion and kaon decay is enhanced
compared to vertical directions. It is easily checked that the
situation is symmetric with respect to
. Therefore, the generation of pions and
kaons in the atmosphere caused by primary cosmic rays induces a
dependence of the atmospheric neutrino flux on
.
In other words, this mechanism does not introduce an up-down
asymmetry in the flux.
- b..
- The geomagnetic field, which acts on the primary cosmic
ray flux, is the cause for a
and
dependence of
the atmospheric neutrino flux. Because of the positive charge of
the primary cosmic rays the neutrino flux is highest (lowest) for
directions coming from the west (east). The nature of the
geomagnetic effects also generates an up-down asymmetry for
low energy neutrinos [215,218].
A third cause for a dependence on the zenith angle is possibly
given by neutrino oscillations which arises because varying
from 0 to
the neutrino path length varies from
around 10 km [232] to around 13000 km. Clearly, this dependence
is not up-down symmetric. If it is disentangled from the
up-down asymmetry caused by the geomagnetic effects, e.g., by using
the high energy component of the neutrino flux, it provides us
with valuable information on neutrino masses and mixing.
Early efforts for detecting atmospheric neutrinos (see Ref. [233]
for a summary) concentrated
on neutrino-induced upward muons, i.e., muons with zenith angles
, or horizontal muons,
i.e., muons with a zenith angle around
, using the process
 |
(5.6) |
where
is a nucleon in the rock surrounding the detector
located underground [234,235]. This method was proposed
to distinguish muons generated by atmospheric neutrinos from
muons originating in the meson decays (5.1) in the
atmosphere. There were two experiments, the Kolar Gold Field
experiment in India [236] and an experiment in South
Africa [237], reporting the first evidence for
atmospheric muon neutrinos. Since these experiments could not
distinguish up and down directions they used horizontal muons
crossing the detector
which lay, due to the great depth of the experimental sites, in
the large zenith angle intervals
and
for the Indian and the South
African experiment, respectively. These early experiments were
accompanied already by atmospheric neutrino flux
calculations [238,239,230,240].
Recent detectors [241] are divided into two classes: water Cherenkov
detectors where the neutrino target is a large volume of water
surveyed by a huge array of photomultiplier tubes sitting on the
surface of the volume (the Kamiokande [242,243,49],
Super-Kamiokande [2] and
IMB [244,50,245] collaborations)
and iron plate calorimeters where
neutrino-induced charged particles ionize the gas between the
plates and the particle paths are reconstructed electronically
(the Fréjus [246], NUSEX [247] and Soudan-2
[51] collaborations).
In contrast to the early detectors
the recent detectors are sensitive to the direction of tracks and can thus
distinguish between up and down through-going tracks. However,
they cannot measure the charge of the leptons
and thus
cannot distinguish between
and
.
Therefore, at low energies they are approximately sensitive the flux
combination
(
)
because at low energies quasi-elastic scattering is predominant
and the ratio of quasi-elastic antineutrino to neutrino
cross sections is approximately 1/3 [233,248,217]. In
atmospheric neutrino physics it is important to distinguish
-like and
-like events which is accomplished by
distinguishing between showers or diffuse Cherenkov rings
(
) and
tracks or sharp Cherenkov rings (
and charged
pions, kaons, protons etc.). The separation between
-like and
-like events is very good, e.g., for Super-Kamiokande its
efficiency is estimated to be 98% or better [2].
For the deep underground detectors two different event classes are
defined [241]. Events in which the neutrino interacts with the material
inside the detector and where all particles from the neutrino interaction
deposit their energies inside the detector are called
contained events. The second class refers to events where muon
neutrinos interact with the material surrounding the detector via
charged current interactions such that the high energy muons enter the
detector [249]. In this way one distinguishes
through-going muons and
stopping muons. Recent measurements of the upward muon
flux were performed by the Baksan [250], Kamiokande
[251,252], IMB [50], Fréjus [246], MACRO
[52] and Super-Kamiokande [1] collaborations.
Super-Kamiokande - and before also Kamiokande
- has an inner detector surrounded by an outer detector. This allows
to further subdivided the contained events into fully contained
events (FC) with all energy of an event deposited in the inner
detector and partially contained events (PC) which have exiting
tracks detected also in the outer detector.
In Kamiokande and Super-Kamiokande, FC
events are separated into those having a visible energy
GeV, the sub-GeV events, and those with
GeV, the multi-GeV events.
Among the contained events in Cherenkov detectors the single-ring
events are well understood, they are predominantly produced by
quasi-elastic scattering of electron and muon neutrinos. In
Kamiokande and Super-Kamiokande for the analyses of FC events only
single-ring events are used with the additional criteria
MeV and
MeV for electron and muon
momenta, respectively, in the case of sub-GeV events.
To quote the numbers of Super-Kamiokande
[2], a Monte Carlo simulation has shown that 88% (96%) of
the sub-GeV
-like (
-like) events are charged current
interactions whereas for multi-GeV events the number is 84%
(99%). The remainder is given by neutral current events.
The PC events were estimated to be 98%
-induced
events for single
and multi-ring configurations and, therefore, all PC events are used for
the analyses.
For the analyses of atmospheric neutrino data it is very important
to have a good understanding of the neutrino interactions in the
detector. The deep-inelastic scattering (DIS) formulas are only valid
for sufficiently high momentum transfer
from the leptons
to the hadrons. For neutrinos with a few GeV interacting with
the detector material, the lowest multiplicity exclusive channels
represent an important fraction of the cross section. Therefore
the following decomposition of the neutrino cross section
[253] has been proposed:
 |
(5.7) |
where QEL indicates the quasi-elastic
and
cross sections and
single pion production. In the case of the latter it is
assumed that
dominates. To avoid double counting in
, the
maximal mass of the pion - nucleon system in the final state is
taken to be
GeV for single pion production whereas the
kinematical region for DIS is bounded by
,
where
is the mass of the hadronic final state.
The first observable to be measured in recent atmospheric neutrino
experiments was the ratio of
-like to
-like events
denoted by
. As discussed earlier, in
flux ratios the large uncertainty in the overall normalization
of the primary cosmic ray flux cancels and there is also
some cancellation of errors in the theoretical calculation.
However, the above ratio is only a limited reflection of the
corresponding ratio of atmospheric neutrino fluxes because of
detector efficiencies and event selection criteria. Thus for the
expected ratio one has to fold the theoretical flux calculations
with the cross sections of the neutrino interactions in the
detector and the detection efficiencies and apply the event
selection criteria. Quoting the result of the Monte Carlo
calculation of Super-Kamiokande as an example, this
collaboration obtains
and
for the sub-GeV ratio of FC events and the ratio considering multi-GeV
FC and PC events, respectively [2], using the neutrino flux
calculations of Ref. [215]. Note that these numbers
significantly deviate from the naive expectation 2. Therefore,
the actual physically relevant quantity is given by the double ratio
 |
(5.8) |
The first indication that this ratio is smaller than 1 was
reported more than ten years ago [254]. In the meantime
the most impressive measurements of
are represented by the
Kamiokande and Super-Kamiokande results:16
 |
(5.9) |
The PC events have been added to the multi-GeV data. The
results
of
the IMB Collaboration [50,255] and
of Soudan-2 [256] are in agreement with those of the
experiments in Kamioka. However, the early iron calorimeter
experiments found values of
compatible with 1, namely
(Fréjus Coll. [246]) and
(NUSEX Coll. [247]).
Apart from the latter two experiments, all others hint at a reduction
of
compared to the expectation
.
Kamiokande and, in particular,
Super-Kamiokande have enough statistics to study the zenith angle
dependence of the measured
and
fluxes.
To this end,
the
interval
is divided into five bins of
length 0.4. Kamiokande [49] has observed a zenith angle
variation of
for the FC multi-GeV + PC events with indications that the zenith
angle variation rather comes from
-like events than from the
-like events. All this is amply confirmed by Super-Kamiokande with
much more statistics and with a significant zenith angle variation in
the
-like events for both, sub-GeV and
multi-GeV.
It is important to study the zenith angle variation independent
of the double ratio
in order to disentangle the zenith angle
dependencies of the electron and muon neutrino fluxes. To this end,
in addition to the oscillation parameters,
also the normalization of primary cosmic ray flux has to be fitted.
Super-Kamiokande has performed a statistical analysis of
the data under the assumption of
oscillations. With
-like and
-like events in five
bins and seven
momentum bins there are altogether 70 data points and three
quantities to be determined: the mixing angle
, the
neutrino mass-squared difference
and the overall neutrino
flux normalization. The best fit gives
and
eV
with
for 67 DOF.
The regions in the
-
plane allowed at
68%, 90% and 99%
confidence level in the case of
oscillations
are shown in Fig. 5.1
[2].
At 90% CL the mass-squared
difference lies in the interval
.
For the simulation of neutrino oscillations the profiles for the
neutrino production heights of Ref. [232] were used.
An analogous procedure with the hypothesis of
oscillations, taking
into account matter effects in the earth, gives a
poor fit with
for 67 DOF.
It is interesting to note that,
as shown in Fig. 5.1,
a fit for Kamiokande with
oscillations gives
with 90% CL [49] for the multi-GeV data
whereas the sub-GeV data show no indication for a zenith angle
dependence of the number of events.
However, Kamiokande has a lower statistics
than Super-Kamiokande.
In Refs. [257,258,259,260,261]
analyses of all available experiments were performed using the
three possible oscillation hypotheses
with similar
results as described above and showing in addition that also the
hypothesis of oscillations into sterile neutrinos gives a satisfactory fit.
For analyses with three neutrinos and including all neutrino
oscillation experiments see [262,263,264,265,266,261].
In Ref. [245] (IMB) no zenith angle variation of
was seen though
with rather small statistics.
In the Super-Kamiokande experiment a significant up-down
asymmetry for the
-like events was found.
The measured value of the asymmetry
as a function of momentum
for
-like and
-like events is shown in Fig. 5.2
[2].
Here
is the number of upward-going events with zenith
angles in the range
and
the number
of downward-going events with
.
The value of the asymmetry for
FC and PC multi-GeV
-like events is
[2]
 |
(5.10) |
The
geomagnetic effect which is one of the causes of an up-down
asymmetry contributes less than
for multi-GeV events
[2]. For the PC events with a mean neutrino energy of
15 GeV [2] this effect is even less important. Note
that for multi-GeV events
the average angle between the incoming neutrino direction
and the charged lepton seen in the single-ring events is of the
order of
or less. Since
corresponds to
, the asymmetry
of the
-like events is in effect an up-down asymmetry of
the
flux. The
result (5.10) constitutes the best indication in favour
of neutrino oscillations found so far with
being 6 standard
deviations away from 0. This is corroborated by
the plots of
as a function of the lepton momentum
in the case of
the single-ring events
shown in Fig. 5.2
[2]. For
-like events,
is a flat function
consistent with no
oscillations. The asymmetry for
-like events analogous to
Eq.(5.10) is compatible with zero:
.
For the
-like events,
starts with
zero and decreases, thus indicating that at small energies both
up and down-going neutrinos oscillate with averaged
probabilities 1/2 whereas at muon energies in the GeV range the
probability of
survival
for the down-going neutrinos approaches 1. This behaviour is
in agreement with the
dependence of the oscillation
probabilities.
It is interesting to note that the experimental values of the
asymmetries
can be used to discriminate between
different oscillation scenarios
[267,268,269,270].
The importance of the results of the atmospheric neutrino
oscillation experiments requires further scrutiny to test
the interpretation in terms of neutrino oscillations. It has
been proposed for Super-Kamiokande to use ratios of charged
current events (CC) to neutral current (NC) events [271,272] in
the spirit of the SNO experiment [273]
in the context of solar neutrinos.
The basic idea is that in Super-Kamiokande NC events could be
seen through
 |
(5.11) |
whereas CC reactions are likely to produce a single charged pion via
 |
(5.12) |
To produce a test
with the reactions (5.11) and (5.12) it is necessary
to discuss their experimental signatures and the contaminations
of these signature with other processes.
Neutral pions are selected by taking events with two diffuse rings from the two
decay photons with an invariant mass between 90 and 180 MeV and
momenta smaller 400 MeV in order to separate the two rings. Such events
were already observed in Kamiokande [274].
Because of the selection criteria for
's,
single pion production is important only for neutrino energies
of order 1 GeV, therefore it involves the sub-GeV events and no
tau lepton will produced in the CC reaction. Defining
via this procedure then the two-ring
ratios [272]
 |
(5.13) |
act as measures of the CC to NC event ratios because in the
ideal situation one can make the identification
,
and
since electrons produce
diffuse and muons and charged pions sharp rings. This picture is
blurred [272,275] by detector efficiencies (e.g.,
it is 0.77 for the
identification of neutral pions via two diffuse rings),
misidentifications (like misidentifying a
whose two
photons cannot be resolved with an electron) and contaminations
(most notably, NC events with
in the final state
contribute to
, electronic CC events with a
whose photons
cannot be resolved add to
and muonic CC events with
the same
configuration add to
). Taking all this into
account [272,275], one has nevertheless reason the expect that
with increasing statistics in Super-Kamiokande the two-ring events can
be used to obtain information on distinguishing between oscillation of
the muon neutrino into tau or sterile neutrinos and to discriminate
between different regions of the parameter space of neutrino mixing
[272,275,276]. Furthermore, it was suggested to use
the asymmetry
for up and down-going
's coming from the NC reaction
(5.11) as an observable to distinguish muon neutrino
oscillations into active neutrinos from those into sterile
neutrinos [277].
One has to assume in this case that
up and down
's originate from up and down neutrinos,
respectively, with high probability. This is not so obvious, however,
since for the identification of the
the events should rather have low
energy as explained above.
The further tests discussed here concern the stopping and through-going
muon events where special efforts have been made to calculate the
fluxes [249,278,253]. Whereas the FC events have neutrino
energies of around 1 GeV the stopping muon events correspond to a mean
neutrino energy of 10 GeV and the through-going muons to 100 GeV
[249,231]. Thus we are now discussing a different range of energy
compared to the discussion above (with the exception of the PC sample).
Obviously, also the zenith angle distribution of upward stopping or
through-going muons can be used to test the neutrino oscillation
hypothesis [279]. Among other experiments (see above)
upward through-going muons have been studied by Kamiokande [252]
and upward muons by MACRO [52]. Kamiokande
has 372 such events above an energy threshold of 1.6 GeV. Fitting the
data to the
oscillation hypothesis yields a best fit with
eV
agreeing
rather well with the Super-Kamiokande result. The analysis of the
MACRO Coll. based on 479 events gives a similar result for
,
however, the zenith angle distribution does not fit very well with the
oscillation hypothesis into tau neutrinos. An attempt has been made to
explain the zenith angle distribution of the MACRO experiment with
oscillations where matter
effects in the earth play a crucial role [280].
In the earth, the density profile can approximately be represented
by constant densities in the mantle and the core, respectively
(see end of Section 4.3). Such a profile can lead to an enhancement of
neutrino transitions due to
the matter effect17in the earth if atmospheric
neutrinos cross the core18 such that the phase picked up by
a neutrino wave function traversing the mantle for the first time and
the phase acquired by traversing the core are each approximately equal
to
. Such an effect for atmospheric neutrinos was recently
considered in detail in Refs. [280,282,283,284,285].
It has also been proposed to exploit this
effect to discriminate between
(no matter effects) and
transitions of atmospheric neutrinos
[280] and to explain an excess of
-like events [286]
possibly seen in the Super-Kamiokande experiment [2].
It was suggested in Ref. [244] to use the observable
 |
(5.14) |
where
and
are the numbers of
stopping and though-going muons, respectively, as an indicator for
neutrino oscillations because this ratio is reduced for neutrino
oscillations with respect to the no-oscillation hypothesis
[253,287]. Also the ratios of ``horizontal'' to ``vertical''
muon events for stopping and through-going muons defined by
 |
(5.15) |
with
 |
(5.16) |
are useful observables as shown in Ref. [287]. In this work a
study was made for Super-Kamiokande, taking
as
the boundary between ``horizontal'' and ``vertical'', showing that the
prospects are good for confirming atmospheric neutrino oscillations found with
the FC single-ring and PC events by using stopping and through-going
muons and the observables
(5.14),
and
(5.15). However, for a precise determination of
these variables are not suitable.
The explanation of the atmospheric neutrino anomaly in terms of
neutrino oscillations
can be checked with long-baseline neutrino oscillation experiments.
The first long-baseline reactor experiment CHOOZ [159]
(see Section 3.2 and Fig. 3.3)
has already excluded atmospheric
oscillations with a large mixing angle for
.
Two other long-baseline reactor experiments are under construction:
Palo Verde [288,289]
and
Kam-Land [290,291].
The Palo Verde experiment has a setup and a sensitivity similar to
the CHOOZ experiment, whereas the Kam-Land experiment,
which is the result of the conversion of the old Kamiokande detector
to a liquid scintillator detector,
will detect
's produced by Japanese reactors
150-200 km away
and will be sensitive to
and a large mixing angle.
The Borexino experiment
(see the end of the Section 5.2 and
[292,291])
will allow to perform a similar measurement.
Accelerator long-baseline experiments will study the oscillation
channels
.
The K2K [271]
experiment,
with a baseline of about 235 km from KEK to Super-Kamiokande
and a neutrino beam with 1.4 GeV average energy,
will be sensitive to
disappearance and
transitions
with
.
A near 1 kton water-Cherenkov detector
will be placed at a distance of about 1 km from the beam dump and
will allow to measure the initial flux and energy spectrum of
's.
This experiment is under construction and
is planned to begin taking data in the year 1999.
Also the MINOS [293] experiment is under construction.
This experiment will have a near detector at Fermilab
and a baseline of about 730 km from Fermilab
to the Soudan mine,
where the far detector will be placed.
The neutrino beam will be produced by protons from the
new Main Injector at Fermilab
and will have an average energy of about 10 GeV.
The far detector is an 8 kton sampling calorimeter made of
magnetized iron and scintillators.
This experiment will be sensitive to
disappearance
and
,
,
transitions,
with the possibility to distinguish the different channels,
for
(the possibility to extend the sensitivity to
lowering the neutrino energy is under study).
In particular
transitions can be revealed through the measurement of a deficit in the
NC/CC ratio.
The MINOS experiment is scheduled to start data-taking around the year 2003.
The ICARUS experiment [294] in Gran Sasso,
constituted of a 0.6 kton liquid argon detector is
scheduled to start in the year 2000.
In the future three new modules with a total mass of 2.4 kton will be installed.
This detector will be sensitive to atmospheric and solar neutrinos
and will allow to reveal long-baseline
oscillations using a neutrino beam
produced at CERN about 730 km away.
Since the average energy of the neutrino beam is rather high,
about 25 GeV,
in order to allow the detection of
through the CC production of
a
,
this experiment will be sensitive to
.
Four other detectors for future LBL CERN-Gran Sasso experiments,
OPERA [295],
NOE [296],
AQUA-RICH [297]
and
NICE [298],
have been proposed and
are under consideration
(see Ref. [299]),
together with the
possibility of a new atmospheric neutrino detector
consisting of a large-mass and high-density
tracking calorimeter
[300].
In Ref. [301] a comparison is made between the possibilities using
atmospheric neutrinos and LBL neutrino experiments for the
determination of the oscillation parameters.
The feasibility to distinguish between atmospheric
and
oscillations using a combination of the results of
future atmospheric, LBL and SBL experiments is discussed in
Ref. [302].
Solar neutrino experiments
The earliest indication in favour of neutrino oscillations
was obtained about 30 years ago in the radiochemical
solar neutrino experiment by R. Davies
et al. [36]. The flux of solar electron
neutrinos measured in this experiment was significantly less than the
predicted one. This phenomenon was called solar neutrino problem.
The existence of this problem was confirmed in
all five solar neutrino experiments
(Homestake [36,37,38],
Kamiokande [39,40,41],
GALLEX [42,43],
SAGE [44,45]
and
Super-Kamiokande [46,47,48])
which measure a flux of electron neutrinos
significantly smaller than the one
predicted by the Standard Solar Model (SSM)
[174,60,303,175,304,305,306,307,308,309,310].
The solar neutrino problem
(see, for example,
[60,311,312,313,314,315,316])
arose in the Homestake experiment by the low counting rate showing
that the flux of
the high-energy
neutrinos
(
)
and of
the medium-energy
neutrinos
(
)
is suppressed by a factor of about 1/3
with respect to the SSM prediction.
In 1988 the solar neutrino problem was confirmed by the
results of the real-time
water-Cherenkov Kamiokande experiment [39]
which measured a flux of
neutrinos of
about half of the SSM flux.
The measurements of the Kamiokande experiment
proved that the observed neutrinos arrive
at the detector from the direction of the sun.
In 1992 the radiochemical GALLEX [42]
and SAGE experiments [45]
succeeded in measuring the neutrino flux
with a low energy threshold
,
which allowed to detect low-energy
neutrinos
produced by the fundamental reaction of the
cycle.
Also these experiments measured a neutrino flux of about half
of the one predicted by the SSM.
Finally,
the Super-Kamiokande experiment
has recently confirmed [46,47,48]
with high statistics
the suppression of the
neutrino flux with respect to the SSM one
by a factor of about 1/2.
The energy of the sun is produced in
the reactions of the thermonuclear
and CNO cycles
shown in Figs. 5.3 and 5.4 (see, e.g.,
Ref. [60]).
The overall result of both cycles is the transition
 |
(5.17) |
where
is the energy release.19Hence,
the production of energy in the sun is accompanied
by the emission of electron neutrinos.
The main part of the solar energy
is radiated through photons
and
a small part (about 2%)
is emitted through neutrinos.
The sources of solar neutrinos are
listed in Table 5.1.
The
,
,
Be,
B and
reactions belong to the
cycle
(see Fig. 5.3),
whereas the
N,
O and
F
reactions
belong to the CNO cycle
(see Fig. 5.4),
which produces only about 2% of the solar energy.
The average and maximum neutrino energies
listed in Table 5.1
are taken from Refs. [317,318].
The neutrino fluxes and predictions
for the neutrino capture rates in the chlorine Homestake experiment
and in the gallium GALLEX and SAGE experiments
given by the
Bahcall-Pinsonneault 1998 (BP98) [304] SSM
are listed in Table 5.2.
Source  |
Reaction |
| Average Neutrino |
Energy
(MeV) |
|
| Maximum Neutrino |
| Energy (MeV) |
|
 |
 |
 |
| |
 |
 |
 |
| |
Be |
 |
|
| |
B |
 |
 |
| |
 |
 |
 |
| |
N |
 |
 |
| |
O |
 |
 |
| |
F |
 |
 |
| |
Source  |
Flux  |
(
) |
|
 |
(
) |
|
 |
| (SNU) |
|
 |
(
) |
|
 |
| (SNU) |
|
 |
 |
- |
- |
 |
| |
 |
 |
 |
 |
 |
| |
Be |
 |
 |
 |
 |
| |
B |
 |
 |
 |
 |
| |
 |
 |
 |
 |
 |
| |
N |
 |
 |
 |
 |
| |
O |
 |
 |
 |
 |
| |
F |
 |
 |
 |
 |
| |
| Total |
|
|
 |
|
| |
As it is seen from the Tables 5.1 and
5.2,
the major part of solar neutrinos are low energy neutrinos
coming from the
reaction.
Monoenergetic neutrinos with
intermediate energy
are produced in the capture of electrons by
Be
and in the
reaction.
High energy neutrinos
are produced in the decay of
B
(the flux of
neutrinos is so small that its contribution
to the event rates of solar neutrino experiments
is negligible).
The flux of
B neutrinos is much smaller
than the fluxes
of
,
Be and
neutrinos.
However,
as we will see later, these neutrinos
give the major contribution to the event rates of experiments
with a high energy detection threshold.
The CNO
,
,
reactions are sources of intermediate energy neutrinos
with a spectrum that extends up to about
.
Their contribution to the event rates of solar neutrino experiment
is small but not negligible.
The neutrino flux
coming from each source
as a function of the neutrino energy
can be written as
 |
(5.18) |
where
is the total flux
and
is the energy spectrum
(
).
The energy spectrum
for each source
is known with negligible uncertainties
[320,317]
because it is determined by the weak interactions
and it is practically independent from solar physics.
On the other hand,
the total flux
of each source
must be calculated with a solar model
and the resulting uncertainties
represent one of the main problem for the interpretation
of the experimental results.
However,
there are some relations
that allow to extract model-independent information
on the neutrino fluxes from the experimental data.
The main one is the luminosity constraint,
which is based on the assumption that the sun is in a stable state
(the energy is produced in the central region
of the sun and for its electromagnetic part
it takes more than
years to reach the surface,
whereas neutrinos escape the sun in about two seconds).
Let us consider a solar neutrino with energy
.
The luminous energy released together with this neutrino
is
.
Multiplying this quantity with the total flux of neutrinos
and
integrating over the neutrino energy one obtains the luminosity
constraint20
 |
(5.19) |
where
is the average energy of the neutrinos from the source
(see Table 5.1),
[35]
is the luminosity of the sun,
[35]
is the sun - earth distance.
The luminosity constraint can be rewritten in the compact form
 |
(5.20) |
where
and
is the solar constant.
Let us emphasize that the luminosity relation
is valid under the assumption that solar
's
on their way to the earth do not transform into other states.
Neglecting the small fraction of energy carried away by neutrinos,
the luminosity constraint gives the approximate value of the
total solar neutrino flux
:
 |
(5.21) |
Since the
nuclei necessary for the
formation of
and
are created by the
or
reactions,
there is another model-independent constraint for
the solar neutrino fluxes of the
cycle
(see [321]):
 |
(5.22) |
Let us now consider the experimental data.
The results of five solar neutrino experiments
are available at present
and are listed in Table 5.3.
| Experiment |
Result |
Theory |
| |
| Homestake [38] |
 |
( ) |
|
 |
 |
| GALLEX [322] |
 |
( ) |
|
 |
 |
| SAGE [323] |
 |
( ) |
|
 |
 |
| Kamiokande [41] |
 |
( ) |
|
 |
 |
| Super-Kamiokande [48] |
 |
(
) |
|
 |
 |
Homestake [36,37,38],
GALLEX [42,43]
and
SAGE [44,45]
are radiochemical experiments.
In the pioneering chlorine Homestake experiment
of R. Davis et. al.,
which started in 1967,
the detector is a tank with a volume of
liters filled with C
Cl
.
Radioactive atoms of
Ar
are produced by solar electron neutrinos through the reaction
[324,325]
 |
(5.23) |
which has an energy threshold
.
The radioactive
Ar atoms that
are created during the time of exposition of each run
(about two months) are extracted from the detector
by purging it with
He
and counted in small proportional counters
which detect the Auger electron produced in the
electron-capture of the
Ar nuclei.
About 0.5 atoms of
Ar are produced every day by solar neutrinos
and about 16 atoms are extracted in each run
(this number is smaller than 30 because of the
Ar lifetime of about 35 days
and because of the extraction efficiency of about 90%).
Since the energy threshold is above the end-point of the
neutrino
spectrum
and the cross section of the detecting process (5.23)
grows with the neutrino energy,
the main contribution
to the counting rate in the
Homestake experiment comes from
B and
Be neutrinos.
According to the BP98 SSM [304]
the event rate in the Homestake experiment should
be21
,
of which
5.9 SNU come from
B neutrinos,
1.15 SNU are produced by
Be neutrinos
and
0.7 SNU are due to
and CNO neutrinos
(see Table 5.2).
The Homestake event rate presented in
Table 5.3
is the event rate averaged over 108 runs [38].
In the radiochemical gallium experiments GALLEX and SAGE
electron neutrinos from the sun are detected through the observation
of radioactive
Ge that is produced in the process
 |
(5.24) |
The GALLEX detector
is a tank containing 30.3 tons of
Ga
(100 tons of a water solution of gallium chloride),
whereas the SAGE experiment
uses about 57 tons of
Ga in metallic form.
The event rate measured in the GALLEX experiment is
events per day.
(Since the
Ge lifetime is about 16.5 days, only around
7 atoms of
Ge are extracted
in each 3-weeks run.)
Since the threshold of the process (5.24) is
,
neutrinos from all sources are detected
in gallium experiments.
According to the SSM,
the contributions to the total predicted event rate
from
,
Be and
B neutrinos are [304]
54%, 27% and 10%,
respectively.
As can be seen from Table 5.3,
both gallium detectors measure an event rate
that is about one half of the SSM prediction.
The weighted average of the GALLEX and SAGE results
yields the event rate
 |
(5.25) |
Taking into account the experimental and theoretical errors,
this rate differs from the SSM prediction by about seven standard deviations!
Both the GALLEX and SAGE detectors have been calibrated using an intense
neutrino source.
The ratio of observed and expected events is
for GALLEX [326]
and
for SAGE [45].
These results demonstrate the absence of unexpected systematic
errors at the 10% level
in both experiments.
In addition,
the GALLEX Collaboration calibrated the detector
by introducing a known number of
radioactive
atoms
in the target solution
[326,327].
The atoms of
resulting from
decay have been extracted in the usual way,
and the As tests prove, at the 1% level, the reliability of the technique
[328] (the number of
atoms
produced in these tests is of the order of
,
whereas the number produced in the
tests is of the order of ten per day).
As emphasized by the GALLEX Collaboration [326,327],
these results rule out the presence of unexpected radiochemical effects
that could explain the deficit of solar
's
measured by the gallium experiments.
It is necessary to emphasize that
the gallium experiments are not only very important for the assessment of
the solar neutrino problem but also for the theory of
thermonuclear energy production in the sun:
They have provided the
first observation of low-energy solar neutrinos
produced in the
reaction
that is the basic reaction of the
cycle and, therefore,
the first direct experimental confirmation of the theory of
thermonuclear origin of solar energy production.
In the
Kamiokande [39,40,41]
and
Super-Kamiokande [46,47,48]
experiments
water-Cherenkov detectors are used for the detection of
solar neutrinos
through the observation of the Cherenkov light emitted by
the recoil electrons in the elastic-scattering process
 |
(5.26) |
Since the direction of the recoil electron is peaked in the direction
of the incoming neutrino,
water-Cherenkov detectors measure the direction of the neutrino flux
and the results of
Kamiokande
and
Super-Kamiokande
have proven that there is a flux of
high-energy neutrinos coming from the sun
with an intensity about half of that predicted by the SSM.
The energy threshold of water-Cherenkov detectors is given by the threshold
for the detection of the recoil electron in the reaction (5.26).
It is higher than in other solar neutrino experiments
because of the large background at low energies:
in the Kamiokande experiment and
in the Super-Kamiokande experiment.
Therefore,
only
neutrinos can be detected in water-Cherenkov experiments.
If nothing happens to electron neutrinos during their trip from
the core of the sun to the earth
(i.e., no neutrino oscillations or other transitions),
the results of the Kamiokande and Super-Kamiokande experiments provide
a measurement of the total flux
of
neutrinos.
Therefore, the results of the Kamiokande and Super-Kamiokande experiments
are usually presented in terms of the measured
flux of
neutrinos. In the Super-Kamiokande experiment
it was found
 |
(5.27) |
which is about one half of the flux predicted by the SSM
(see Table 5.3).22
From the experimental results and
the most updated theoretical predictions [304]
listed in Table 5.3
one can see that the observed event rates in all solar neutrino experiments
are significantly smaller than the predicted rates
(see also [329]
where the comparison with other theoretical predictions is discussed).
This discrepancy constitutes the solar neutrino problem.
The SSM is robust and it has been recently
tested in a convincing way by comparing its predicted
value for the sound speed in the interior of the sun
with precise helioseismological measurements
(see [330,331,304,315,332,333]).
However,
it is clear that a model-independent proof of the existence
of the solar neutrino problem
would be more convincing.
Such a model-independent approach has been discussed in several papers
[334,308,335,336,337,338,339,340,341,321,312,342,329,343]
using the model-independent luminosity constraint (5.20)
and the fact that the energy spectra
Eq.(5.18)
of the various neutrino sources are practically independent
from solar physics
[320,317].
From the luminosity constraint it is possible to obtain a model-independent
lower bound on the gallium event rate.
Indeed,
since
neutrinos have the smallest average energy
and
the
-
cross section
increases with the neutrino energy
,
we have
 |
(5.28) |
where
is the average cross section of neutrinos from the source
(see Table 5.2).
Since the luminosity constraint (5.20) implies that
 |
(5.29) |
we obtain
 |
(5.30) |
where we used the value of
given in Table 5.2
and
.
The lower bound (5.30)
is just compatible with the
combined result (5.25) of the gallium experiments.
This means that the results of the GALLEX and SAGE experiments
can be explained if practically only
neutrinos are emitted by the sun.
This possibility is incompatible with any solar model
constrained by the helioseismological data
[304,315].
More stringent model-independent conclusions
on the existence of the solar neutrino problem
can be obtained by comparing the results of different solar
neutrino experiments.
If the survival probability of solar
's is equal to one,
the result of the Super-Kamiokande experiment gives
the value (5.27) for the flux of
neutrinos,
whose contributions to the event rates of the chlorine and
gallium experiments are
| |
|
 |
(5.31) |
| |
|
 |
(5.32) |
Subtracting the contribution (5.31) of
's from the event rate measured in the Homestake experiment
(see Table 5.3),
one obtains
 |
(5.33) |
Since the quantity
represents the contribution to the chlorine event rate due to
,
,
and CNO neutrinos,
it must be positive (or zero)
but the result in Eq.(5.33)
is only marginally compatible with a positive value
(the probability is less than 26%).
Furthermore,
the result Eq.(5.33)
shows that the fluxes of intermediate energy neutrinos
(
,
and CNO neutrinos)
are strongly suppressed with respect to the SSM prediction
(see Table 5.2).
This is incompatible with any solar model
constrained by the helioseismological data
[304,315]
and with the fact that,
since both the reactions that produce
and
neutrinos
originate from
nuclei
(see Fig. 5.3),
it is very difficult
to have a suppression of the
neutrino flux
(with respect to the SSM value)
that is stronger than the suppression of the
neutrino flux.23
Let us consider now the gallium experiments.
Subtracting the contribution of
neutrinos from the luminosity constraint (5.20),
we have
 |
(5.34) |
Hence,
following the same reasoning as in Eq.(5.28)
we obtain
 |
(5.35) |
and the contribution of
,
,
and CNO neutrinos
to the gallium event rate cannot be bigger than
 |
(5.36) |
This result is compatible with a positive value
with a probability smaller than 6%.
Again,
the result (5.36) implies that
the fluxes of intermediate energy neutrinos
(
,
and CNO neutrinos)
are strongly suppressed with respect to the SSM prediction,
in contradiction with any solar model
constrained by the helioseismological data
and with the moderate suppression of the
neutrino flux.
Further model-independent methods for proving
the existence of the solar neutrino problem
on the basis of the data of solar neutrino experiments
are nicely discussed in Ref. [312].
In the following
we will assume that there is a solar neutrino problem
that is caused by neutrino oscillations.
The solar neutrino data have been analysed in many papers
under the assumption of two-neutrino mixing
(see [342,344,329]
and references therein)
and in a few papers
under the assumption of three-neutrino mixing
(see [345,346,347,348]
and references therein).
Here we will briefly review the results
of the most updated two-generation analysis
[329]
of all the solar
neutrino data,
including those obtained in the first 504 days of operation
of the Super-Kamiokande experiment
(see Table 5.3).
The deficit of solar
's can be explained
in terms of two-generation neutrino mixing
either through vacuum oscillations
or through MSW resonant transitions in matter [164,163].
Furthermore,
the two mixed neutrinos can be the electron neutrino
and another active (
or
) neutrino
or
the electron neutrino and a sterile neutrino.
There are two differences between transitions of
solar
's into active and sterile neutrinos:
- The probability of MSW transitions
is different in the two cases,
because sterile neutrinos do not have the weak-interaction
potential due to elastic forward scattering
that is present for active neutrinos.
- The active neutrinos
(
's or
's)
that are produced by the oscillation mechanism
(either in vacuum or in matter)
contribute to the event rates measured in the Kamiokande and Super-Kamiokande
water-Cherenkov experiments
(notice, however,
that the cross section of
-
and
-
scattering
is about six times smaller than the
-
cross section).
The formulas for the survival probability of solar
's
in the case of vacuum oscillations and MSW transitions
are given in Eq.(3.21) (with
)
and
in Eq.(4.41) together with Eq.(4.47) [193],
respectively.
These formulas depend on the two mixing parameters
and
.
The allowed regions in the
-
plane,
obtained in Ref. [329]
from the fit
of the measured event rates listed in
Table 5.3
by using the BP98 SSM [304] and the analytic formulas in
Ref. [193] for the MSW case,
are shown in Figs. 5.5-5.6.
Figure 5.5
[329]
shows the three allowed regions
in the case of
MSW
or
transitions.
They are the small mixing angle (SMA-active) region at
 |
(5.37) |
and the large mixing angle (LMA) region at
 |
(5.38) |
The best fit, with a confidence level of 19%,
is obtained in the SMA region,
whereas the LMA region has a confidence level of 4%.
In Ref. [329] there is also
the so-called low mass (LOW) region at
and
,
however, the LOW region is only marginally acceptable,
with a confidence level of 0.7%.
The allowed regions
in the
-
plane
in the case of
or
vacuum oscillations are shown in Fig. 5.5
[329].
These regions extend over large ranges of
and
around the best fit values
 |
(5.39) |
The confidence level of the fit is 3.8%.
There is no allowed region in the case of
vacuum oscillations
(the fit has a confidence level of 0.05%)
[329,349].
As shown in Fig. 5.7
[329],
only the small mixing angle (SMA-sterile) region at
 |
(5.40) |
is allowed in the case of
MSW
transitions,
with a confidence level of 19%
(the large mixing angle and low mass solutions
have a confidence level of 0.001%
and 0.003%,
respectively) [329,350].
The authors of Ref. [329]
performed also global fits of the solar neutrino data
including the energy spectrum and zenith-angle
distribution
of the recoil electrons
measured in the Super-Kamiokande experiment
[48].
As noted in Ref. [329],
since these data are still preliminary,
the results of this analysis are less robust than those
obtained fitting only the global rates.
However,
it is interesting to note that
in the case of
MSW
or
transitions
only the SMA-active region remains allowed,
as shown in Fig. 5.8
[329],
with a confidence level of 7%.
During the night, solar neutrinos
pass through the earth and the matter effect can
cause a regeneration of
's
(see Ref. [316] and references therein) and, therefore, a
zenith angle dependence of the solar neutrino flux.
The size of this effect depends on the values of
and
and is sizable only for large values of the mixing angle.
The preliminary value of the day-night asymmetry
of solar neutrino events measured in the Super-Kamiokande experiment is
[48]
 |
(5.41) |
compatible with zero.
Hence this experimental result
tends to exclude solutions of the solar neutrino problem
with a large mixing angle.
As was shown in Refs. [283,351] the
step-like profile of the matter density of the earth could lead to an
enhancement of neutrino transitions not only for atmospheric neutrinos
(see Section 5.1.3) but also for solar neutrinos.
In conclusion of this section,
we would like to emphasize that the results of solar neutrino experiments
provide a rather strong indication on favour of neutrino mixing
and several experiments are under construction
[273,352,322]
or
in project
(see [353,354] and references therein).
In particular,
the measurement of the electron neutrino spectrum
at SNO may allow to obtain model-independent information
on the neutrino oscillation probability and on the flux
of
neutrinos
[355,356,357,358,359],
the measurement of the flux of
neutrinos
on the earth in the Borexino experiment
and
the results of Super-Kamiokande, SNO, Borexino and GNO
will allow to distinguish the different
possible solutions of the solar neutrino problem
(see
Refs. [350,360,361,362,363] and references therein).
The Sudbury Neutrino Observatory (SNO) [273,364]
is located 6800 feet below ground in the Creighton mine,
near Sudbury in Ontario (Canada). In the SNO experiment
a Cherenkov detector
with 1 kton of heavy water (D
O) contained in an
acrylic vessel of 12 m diameter will be used.
The Cherenkov light is detected with a geodesic array
of
photomultiplier tubes
surrounding the heavy water vessel.
The heavy water detector is immersed in normal water in order to reduce the
background.
Solar neutrinos will be observed in real-time through the
charged-current (CC), neutral-current (NC) and
elastic scattering (ES) reactions
Since the energy threshold for the observation of the recoil
electrons in the CC and ES
processes is about 5 MeV
and the neutrino energy threshold for the NC reaction is 2.2 MeV,
only
B neutrinos can be observed.
The event rates predicted by the SSM are around
23 per day for the CC reaction,
7 per day for the NC reaction and
3 per day for the ES reaction.
From the hit pattern of the photomultipliers, the direction and
energy of the neutrino
in the CC reaction can be reconstructed,
allowing a direct measurement of the spectrum of the
solar electron neutrino flux on the earth.
The observation of a distortion of this spectrum
with respect to the one calculated without neutrino oscillations
will represent a model-independent
proof of the occurrence of neutrino oscillations.
The NC reaction will allow to detect all active neutrinos
with the same cross section
(because of the
-
-
universality of weak interactions),
whereas the cross sections of
and
scattering
in the ES reaction is about six times smaller than that of
.
The measurement of the neutral current reaction will
allow to determine the total flux of active neutrinos from the sun
[355,359],
which coincides with the produced flux if there are only
active-active neutrino oscillations,
which can be revealed trough a measurement of
a deficit in the NC/CC event rate.
On the other hand,
if the NC/CC event rate will not show any deficit and
a distortion of the electron neutrino spectrum will be observed,
it will mean that there are transitions of solar
's into sterile states
[356]
(of course, complicated scenarios with transitions
solar
's into active and sterile states are not excluded
and must be taken into account).
The SNO experiment
is expected to release the first results in the year 1999.
The Borexino experiment [352,365]
is designed to detect low-energy
solar neutrinos
in real-time
through the observation of the elastic scattering
process (5.44)
with an energy threshold for the recoil
electron of 250 keV.
It will use 300 tons of liquid scintillator
in an unsegmented detector with 2000 photomultiplier tubes.
The event rate predicted by the SSM for a fiducial volume of about 100 tons
is about 50 events per day,
mostly generated by the 0.86 MeV monoenergetic line of
Be solar neutrinos.
Since this line gives a characteristic spectral signature in the ES process,
the flux on the earth of
Be solar neutrinos will be determined
and it will be possible to check if
it is suppressed with respect to the one predicted by the SSM
as suggested by the results of current experiments.
The neutrino oscillation solutions predict a measurable flux:
events per day for the SMA MSW solution,
events per day for the LMA MSW solution
and
events per day for the vacuum oscillation solution
(see [360]).
The Borexino experiment is
under construction in the Laboratori Nazionali del Gran Sasso in Italy
and is scheduled to start data taking around the year 2001.
The Liquid Scintillator Neutrino Detector (LSND) experiment
[366,53,367,368,369,54,370,371]
at the Los Alamos Meson Physics Facility (LAMPF) is a SBL neutrino
oscillation experiment with a baseline of about 30 meters.
It utilizes the LAMPF proton beam with a kinetic energy of 800 MeV which
produces pions by hitting a 30 cm long water target
located about 1 m upstream from a copper beam stop.
Most of the
's
are stopped in the target
and decay into muons which come to
rest and decay in the target as well.
The decay at rest of the positively charged muons allows
to investigate
oscillations
[366,53,367,370,371].
A small fraction of the positively charged pions (
) decays in
flight in the 1 m long space between the target
and the beam stop
and is used for the study of
oscillations
[54,369,370,371]
(only
of the
's
decay in flight and produce a small contamination of
's in the
beam).
The detector is an approximately cylindrical tank
8.3 m long by 5.7 m in diameter,
filled with a scintillator medium
consisting mainly of mineral oil (CH
).
In the following we review
the results obtained by the LSND Collaboration
from the analyses of the data from
decay at rest
and
decay in flight.
[366,53,367,370,371].
In this case one has the decay chain
 |
(5.45) |
Note that in these decays no
appears,
which makes the
study of
transitions possible.
In
decay at rest the upper bound on the
energy is
MeV (
is the muon mass)
and the energy spectrum of the
flux is very well known,
being determined by the kinematics of muon decay.
The electron antineutrinos are detected by
 |
(5.46) |
The
signature is a coincidence between
and a delayed
from the capture
.
There is a background of
from
the decay chain symmetrical to Eq.(5.45)
starting with
. This
flux
is suppressed relative to the
flux from the
chain by a factor of
for three
reasons:
production in the beam stop is a factor of 8
smaller than
production, only 5 % of the
decay
in flight (all stopped
are absorbed) and of the
resulting
only 12 % decay before they are captured by
nuclei. In the latter case they give rise to
but not to
.
The LSND apparatus does not distinguish between
and
.
Therefore, the reaction
gives a signal like a
but can only occur for
MeV.
LSND has seen an excess of
events with energies between 20
and 60 MeV of
corresponding to a
transition
probability of
[53].
In Refs. [370,371] preliminary results can be
found for the data taken from 1993 to 1997.
The value for the transition probability
agrees with the previous one,
except for a reduced statistical error of
instead of
.
The plot of the favoured region
in the
-
plane obtained by the LSND collaboration
is shown in Fig. 5.9. In the positron energy range between
36 and 60 MeV the background is reduced because in this
energy range there are no electrons from
induced by electron neutrinos from
decay at rest [53]. LSND has
observed 22 events with such positron energies compared to a background
of only
events [53].
[54,369,370,371].
The
decay in flight generates a
beam with which
oscillations are studied. The electron
neutrinos are detected by electrons from the inclusive reaction
 |
(5.47) |
where the electrons have an energy
between 60 and 200 MeV.
In this case the signal is a single electron. The
electron neutrinos from
decay in flight are suppressed by the longer muon lifetime and
the kinematics of the three-body decay, whereas those from
by the small branching ratio of
. The lower bound of 60 MeV on
lies
well above
MeV, the endpoint of the electron
spectrum from muon decay at rest.
Such muons can be induced by beam-related
events. Also
from
decay at rest (5.45) generate
electrons below this energy. LSND has seen 40 electron events
compared to a background of
. This corresponds to
a
transition probability of
[54],
consistent with the one
determined by muon decay at rest.
The favoured region
in the
-
plane obtained by the LSND collaboration
is shown in Fig. 5.10 as the area between the two solid lines
[54].
The regions in Fig. 5.10 enclosed by dashed lines
are the LSND favoured regions obtained with the 1996
muon decay at rest analysis [53,367].
One can see that the two results are perfectly compatible.
Recently, the LSND collaboration has reported [376]
new evidence in favour of
oscillations using the
beam produced by
decay in flight
and detecting the
's through the observation of the reaction
 |
(5.48) |
where
is the ground state of
,
which decays with a lifetime of 15 ms into
.
A temporal (
)
and
spatial (
)
correlation between the observations of the
and
allows an almost background-free measurement.
The LSND collaboration reported [376]
the observation of 5 events,
with an expected background of 0.5 events.
These numbers are in agreement with the results
of the other two analyses of the LSND data
summarized above.
The region of neutrino mixing parameters indicated by the results
of the LSND experiments is presently explored
also by the KARMEN experiment
[377,378,379,375,380,381,382].
The KARMEN experiment is located at
the ISIS spallation neutron facility of the Rutherford Laboratories (UK).
At ISIS a pulsed proton beam is used:
two 100 ns wide proton pulses separated by 330 ns
are produced every 20 ms.
The proton pulses are directed on a target where they produce
positive pions
(the negative pions are absorbed in the source before they decay)
which decay at rest according to the chain (5.45)
producing an equal number of
,
and
.
The KARMEN Collaboration searches for
at a distance of about 18 m. The time structure of the neutrino beam
is important for the identification of the neutrino induced reactions
in the KARMEN detector
and for the effective suppression of the cosmic ray background.
The KARMEN experiment started in 1990 and ran until 1995 as KARMEN 1
[378].
In 1996 the experiment was upgraded to KARMEN 2,
eliminating the main cosmic ray induced background component
in the search for
oscillations
[379].
So far, the KARMEN 2 experiment
measured no events
[375,380,381,382],
also no background events, of which
events were expected.
The resulting exclusion curve obtained with a Bayesian analysis
[375]
is shown in Fig. 5.9.
One can see that the null result of the KARMEN 2 experiment
excludes part of LSND-allowed region in favour of
oscillations, but at present
the sensitivity of the KARMEN experiment is not enough to
check all the LSND-allowed region.24
The KARMEN experiment is expected to reach a sensitivity
that will allow to cover the LSND-allowed region in one or two years
[379].
If the KARMEN experiment will find an oscillation signal
compatible with the one observed in the LSND experiment,
the LSND evidence in favour of
oscillations will be confirmed.
On the other hand,
if the KARMEN experiment will not find any neutrino oscillation signal,
then the KARMEN exclusion curve will lie on the left of the
LSND allowed region in the
-
plane,
but not far away from it [379].
Hence,
it is important to test the LSND result with
a new experiment with much higher sensitivity.
Four such experiments have been proposed and are under study:
BooNE [387] at Fermilab,
I-216 [388] at CERN,
ORLaND [389] at Oak Ridge
and
NESS at the European Spallation Source [390].
Among these proposals BooNE is approved and will start in 2001.
We would like to finish this section by emphasizing the importance
of the results obtained in the LSND experiment for
neutrino physics.
As we will discuss in the next section,
the LSND evidence in favour of SBL
oscillations
taken together with the indications in favour of
solar and atmospheric neutrino oscillations
require the existence of at least one sterile neutrino,
which represents a manifestation of new physics beyond the Standard Model.
Hence it is important that the
LSND indication is checked as soon as possible by another experiment
with high sensitivity.
Analysis of neutrino oscillation data
As we have discussed in the previous section,
at present there are three
indications in favour of non-zero neutrino masses and mixing.
The analysis of the experimental data
(based on the simplest assumption of two neutrino mixing)
indicate that there must be at least three
different
,
each one relevant for the oscillations of neutrinos in one type of
experiments:
(MSW)
or
(Vac. Osc.)
is relevant for oscillations of solar neutrinos,
is relevant for oscillations of atmospheric neutrinos,
and
is relevant
for SBL neutrino oscillations according to the LSND experiment.
Three different
require the mixing of at least four massive
neutrinos.
Here we consider two possible scenarios:
- Only the indications in favour of neutrino mixing obtained in
atmospheric and solar neutrino experiments are confirmed by future
experiments.
In this case the mixing of
three massive neutrinos can accommodate all data.
- All the existing indications in favour of neutrino mixing,
including that obtained in the LSND experiment,
are confirmed by future experiments.
In this case,
the minimal mixing scheme that allows to accommodate all data
is one with mixing of four neutrinos.25.
Mixing of three massive neutrinos
Let us consider the mixing of three massive neutrinos,
i.e.,
in Eqs.(3.3) and (3.9),
 |
(6.1) |
and assume the neutrino mass hierarchy26
 |
(6.2) |
that is indicated by the see-saw mechanism
(see Section 2.5)
and by the models beyond the Standard Model
(see Section 2.6).
In this case,
the only available possibility for the explanation
of the solar and atmospheric neutrino
anomalies through neutrino oscillations is that
is relevant for the transitions of
solar neutrinos and
is
relevant for the transitions of atmospheric neutrinos:
| |
|
 |
(6.3) |
| |
|
 |
(6.4) |
Let us consider the oscillations of atmospheric and LBL neutrinos
and assume that
 |
(6.5) |
The distance
travelled by atmospheric neutrinos
ranges from
for downward-going neutrinos
to
for upward-going neutrinos
and the energy
relevant for the oscillations of
atmospheric neutrino experiments
ranges from
to
.
Therefore,
if the solar neutrino problem
is due to MSW transitions,
the condition (6.5)
is not satisfied by atmospheric neutrinos with energies below
.
In this case,
the following discussion applies only to high-energy
atmospheric neutrinos.
On the other hand,
the inequality (6.5) is valid in any case for
reactor LBL experiments with a source - detector distance
and
for accelerator LBL experiments
with a source - detector distance
and neutrino energy
.
If the inequality (6.5) is satisfied,
from Eq.(3.12), we have
for the probability
of
transitions
![\begin{displaymath}
P_{\nu_\alpha\to\nu_\beta}^{(\mathrm{atm,LBL})}
=
\left\vert...
...lta{m}^2_{31} L }{ 2 E } \right) - 1 \right]
\right\vert^2
\,.
\end{displaymath}](img1027.gif) |
(6.6) |
Therefore,
a mass hierarchy of three neutrinos implies that
the transition
probabilities in atmospheric and LBL experiments
depend only on the largest mass squared difference
and on the
elements
that connect flavour neutrinos with
the heaviest neutrino
.
From Eq.(6.6) we obtain
[402,403,404,405,192]
| |
|
 |
(6.7) |
| |
|
 |
(6.8) |
with the oscillation amplitudes
| |
|
 |
(6.9) |
| |
|
 |
(6.10) |
The expressions (6.7) and (6.8)
have the same dependence on the
quantity
as the standard expressions
(3.20) and (3.21)
that describe neutrino oscillations in the case of two flavours,
with
replaced by
.
This is due to the fact that only one mass-squared difference is
relevant for the oscillations of atmospheric and LBL neutrinos in the case of
the mass hierarchy (6.2).
The expressions (6.9) and (6.10) describe, however,
all possible oscillations between three flavour neutrinos:
,
and
.
The unitarity of the mixing matrix implies that
and
LBL and atmospheric neutrino oscillations in the scheme under consideration are
described by only three parameters:
,
and
.
Let us also stress that the transition probabilities (6.7)
are independent on the phases of the mixing matrix
.
Thus, in the case of a hierarchy of neutrino masses
the relation (3.18) is automatically satisfied and effects of CP
violation in the lepton sector cannot be revealed
in neutrino oscillation experiments.
As it is seen from Eq.(6.8),
the
survival probability
in reactor LBL experiments
depends on the parameters
and
.
This fact allows to obtain information on the value of
in LBL reactor experiments
and, in particular,
from the results of the first of them,
CHOOZ [159],
in which no indications in favour of neutrino oscillations were found.
The exclusion curve in the plane of the two-neutrino mixing parameters
and
that
was obtained in the CHOOZ experiment
is shown in Fig. 3.3.
Comparing Eqs.(3.21) and (6.8),
one can see that
in the three-neutrino scheme under consideration in this section
 |
(6.11) |
Therefore,
we can obtain information on the parameter
from the exclusion plot obtained in the CHOOZ experiment.
From Eq.(6.11) we have
 |
(6.12) |
and the exclusion curves shown in Fig. 3.3
imply the constraint
 |
(6.13) |
with a value of
that depends on
.
From Eqs.(6.12) and (6.13)
it follows that
must satisfy one of the two inequalities
 |
(6.14) |
with
 |
(6.15) |
From the CHOOZ exclusion curve in Fig. 3.3 (90% CL) we find
for
eV
and therefore
 |
(6.16) |
Solar neutrinos
Large values of
as those implied by the second inequality
in Eq.(6.16) are excluded by the data of solar
neutrino experiments.
In order to prove this statement,
let us derive the expression
for the averaged survival probability of solar
electron neutrinos in the case of vacuum oscillations.
From Eq.(3.12) with
,
taking into account that for solar neutrinos
,
the survival probability averaged over the unobservable
fast oscillations due to
is given by
 |
(6.17) |
This expression can be rewritten in the form
 |
(6.18) |
where
is the two-generation survival probability
 |
(6.19) |
(see Eq.(3.21))
with
and the mixing angle
determined by
 |
(6.20) |
The
relation (6.18) is also valid if matter effects are taken into
account (see Eq.(4.48) and Ref. [406]).
In this case,
is the two-neutrino
survival probability in matter given in Eq.(4.41),
with
replaced by
and
by
.
Considering now the constraints (6.16) on
,
if
is large (
),
from Eq.(6.18)
one can see that
.
Such a large value of
is incompatible with the deficit of solar
's observed in all solar neutrino experiments.
Therefore,
from the results of the CHOOZ experiment we reach the conclusion that
the element
is small:
 |
(6.21) |
This limit implies that
 |
(6.22) |
Hence,
the oscillations of solar neutrinos are practically independent
of the value of
and are decoupled from the oscillations of atmospheric
and LBL neutrinos,
which depend only on the elements
of the mixing matrix
(see Eqs.(6.9) and (6.10)).
This means that
the two-generation analyses of the solar neutrino data
are appropriate also
in the three-neutrino scheme
and they give information on the values of
,
and
[407].
According to a recent analysis of the solar neutrino data
[344],
which include preliminary data from Super-Kamiokande
[46,408],
the ranges of the mixing parameters allowed at 90% CL
for the small and large mixing angle MSW solutions
and for the vacuum oscillation solution are,
respectively,
Therefore,
taking into account that
is small
and assuming that
(this choice is necessary only for the MSW solutions in order to fulfill
the resonance condition (4.34)),
we have
Notice that these ranges are statistically rather stable.
For example,
the range of
allowed at 99% CL
in the case of the large mixing angle MSW solution is
[344],
which implies
,
(compare with Eq.(6.27)).
Atmospheric neutrinos
The atmospheric neutrino data of the Super-Kamiokande experiment
give an additional indication that
is small.
Indeed, the up-down asymmetry of electron
events
that could be
generated by
or
(with
)
transitions is compatible with zero (see Section 5.1.3) [2]:
 |
(6.29) |
In the three-neutrino scheme under consideration,
the Super-Kamiokande atmospheric neutrino data can be explained by
oscillations [2,409,261].
The analysis of the data under the assumption of two-generation mixing
yields the following constraints:
 |
(6.30) |
from which it follows with Eqs.(6.7) and (6.9)
that
and
must be close to 0.5 and therefore
small.
Coming back to Eq.(6.21) which implies that
we see that
the CHOOZ upper bound for
is not very small.
It is suggestive, however, to assume that
 |
(6.31) |
i.e. that the element
is so small
that it is unimportant for neutrino oscillations
[410,407].
Such an assumption has been recently adopted by several authors
[411].
Let us infer some consequences of the assumption (6.31).
First of all,
it is easy to see that if the inequality (6.31)
is satisfied,
matter effects are not essential for
transitions of atmospheric and LBL neutrinos.
Indeed,
the evolution equation for neutrinos in matter in the flavour
representation is (see Section 4)
 |
(6.32) |
with
 |
(6.33) |
and
 |
(6.34) |
Remember that for anti-neutrinos the matter potential term
must be replaced by
.
The inequality (6.5)
implies that the phase generated by
can be neglected for atmospheric neutrinos
and
can be approximated by
 |
(6.35) |
In this case
(taking into account that
the phases of the matrix elements
can be included in the charged lepton fields)
we have
 |
(6.36) |
If the condition (6.31) is satisfied,
the elements of the matrix
with indices
or/and
can be neglected.
This means that the equation for
the wave function
is decoupled from those of
,
:
| |
|
 |
(6.37) |
| |
|
 |
(6.38) |
These equations imply that
and
is given by the two-generation expression for neutrino oscillations in vacuum,
 |
(6.39) |
(see Eq.(3.20)),
with
and
the mixing angle
determined by
 |
(6.40) |
Therefore,
if the condition (6.31) is satisfied,
there are no
oscillations
in atmospheric and LBL neutrino experiments
and matter effects are not important for
transitions.
This means that
the two-generation analyses of
the atmospheric neutrino data in terms of
are appropriate
in the three-neutrino scheme under consideration and
yield information on the values of the parameters
and
.
According to a recent analysis [409]
of the atmospheric neutrino data,
the ranges of
and
for
oscillations
allowed at 90% CL
by the Super-Kamiokande multi-GeV data
and by the data of atmospheric neutrino oscillation experiments
are, respectively,
Taking into account the comments after Eq.(6.5),
the ranges (6.42) of the two-generation mixing parameters
can be used in order to constrain the values of
and
only in the case of the vacuum oscillation solution
of the solar neutrino problem.
If the solar neutrino problem is due to MSW transitions,
we must use the less restrictive ranges (6.41).
Thus, assuming that
,
we have
As in the case of the ranges (6.26)-(6.28),
also the ranges (6.43)-(6.44)
are statistically rather stable.
For example,
the range of
allowed at 99% CL
by all atmospheric neutrino data is
[409],
which implies
,
(compare with Eq.(6.44)).
The mixing matrix in the case
We have seen that,
if we assume that
is very small,
the results of the two-generation
analyses of the data of solar and atmospheric neutrino experiments allow to
determine the moduli of the elements of the first line and of the third
column of the mixing matrix
.
The moduli of the other elements can be
determined from the unitarity of
.
The simplest way to do it is
to start from the parameterization (2.20)
of the
mixing matrix
.
A very small
implies that
.
Since in the parameterization (2.20) the CP-violating phase
is associated with
,
it follows that CP violation
is negligible in the lepton sector (see also footnote 5)
and we have
 |
(6.45) |
Using the information on
and
given by the two-generation analyses of the results of solar
and atmospheric neutrino experiments, we obtain
for the moduli of the elements of the mixing matrix [407]:
In the case of the small mixing angle MSW solution
of the solar neutrino problem
could be of the same order of magnitude as
.
The best fit of Super-Kamiokande data corresponds to maximal mixing,
and the vacuum oscillation solution of the solar neutrino problem
includes the maximal mixing case
.
The combination of these two possibilities gives the
bi-maximal mixing matrix
 |
(6.49) |
that has been assumed recently by several authors [411].
Notice that bi-maximal mixing
is not compatible with the large mixing angle MSW solution
of the solar neutrino problem
[412].
It is interesting to note that,
because of the large mixing of
and
with
,
the transitions of solar
's
in
's and
's
are of comparable magnitude.
However,
this phenomenon
and the values of the entries
in the
-
sector of the mixing matrix
cannot be checked with solar neutrino experiments
because the low-energy
's and
's
coming from the sun can be detected only with neutral-current
interactions,
which are flavour-blind.
Moreover,
it will be very difficult to check
the values of
,
,
and
in laboratory experiments
because of the smallness of
.
In the derivation of Eqs.(6.46)-(6.48)
we have assumed that
and
.
The other possibilities,
and
,
are equivalent, respectively,
to an exchange of the first and second columns
and
to an exchange of the second and third rows
in the matrices (6.46)-(6.48).
Unfortunately,
these alternatives are hard to distinguish experimentally
because of the above mentioned difficulty
to measure directly the values of
,
,
and
.
Only the choice
,
which is necessary for the MSW solutions
of the solar neutrino problem,
could be confirmed by the results of the new generation
of solar neutrino experiments
(Super-Kamiokande,
SNO,
ICARUS,
Borexino,
GNO
and others
[353])
if they will allow to exclude the vacuum oscillation solution.
Accelerator long-baseline experiments
Future results
from reactor long-baseline neutrino oscillation experiments
(CHOOZ [413],
Palo Verde [288],
Kam-Land [290])
could allow to improve the upper bound (6.21)
on
.
In this section we discuss how
an improvement of this upper bound
could be obtained by future
accelerator long-baseline neutrino oscillation experiments
that are sensitive to
transitions
(K2K [271],
MINOS [293],
ICARUS [294]
and others [296,297,295,298];
for a review see Ref. [414]).
If matter effects are not important,
in the scheme under consideration
the parameter
measured in
long-baseline experiments is given by
(see Eq.(6.9))
 |
(6.50) |
If accelerator long-baseline neutrino oscillation experiments
will not observe
transitions
and will place an upper bound
,
it will be possible to obtain the limit
 |
(6.51) |
where
is the minimum value of
allowed by the solution of the atmospheric neutrino anomaly
and
by the observation of
long-baseline transitions.
For example,
taking
(see Eq.(6.44))
we have
.
If a value of
,
corresponding to the sensitivity of the ICARUS experiment
for one year of running [294],
will be reached,
it will be possible to
put the upper bound
.
The observation of
transitions in long-baseline experiments
will allow to establish a lower bound for
because
the parameter
is given in the scheme under consideration by
(see Eq.(6.9))
 |
(6.52) |
From the unitarity relation
it follows that
an experimental lower bound
allows to constrain
in the range
 |
(6.53) |
If
is found to be close to one,
as suggested by the solution
of the atmospheric neutrino problem
(see Eqs.(6.41) and (6.42)),
the lower bound
is close to
.
If matter effects are important,
the extraction of an upper bound for
from the data of
accelerator long-baseline experiments
is more complicated.
In this case the probability
of
oscillations
is given by [192]
 |
(6.54) |
where
is the neutrino energy and
is the distance of propagation.
This probability depends on the neutrino energy
not only through the explicit
in the denominator of the phase,
but also through the energy dependence of
.
For long-baseline neutrino beams
propagating in the lithosphere of the earth
the charged-current effective potential
is practically constant:
and
(see Eq.(4.37)).
If long-baseline experiments will not observe
transitions
(or will find that they have an extremely small probability)
for neutrino energies such that
,
it will mean that
is small
and a fit of the experimental data with the formula (6.54)
will yield a stringent upper limit for
(taking into account
the lower limit
obtained from the solution of the atmospheric neutrino anomaly
and
from the observation of
long-baseline transitions).
On the other hand,
the non-observation of
transitions
for neutrino energies such that
does not provide any information on
because in this case the transition probability (6.54)
is suppressed by the matter effect.
Hence,
in order to check the hypothesis
,
as well as to have some possibility to observe
transitions if this hypothesis is wrong,
it is necessary that a
substantial part of the energy spectrum of the neutrino beam
lies below
 |
(6.55) |
This requirement will be satisfied in the
accelerator long-baseline experiments under preparation
(K2K [271],
MINOS [293],
ICARUS [294]
and others [296,297,295,298];
see [414])
if
is not much smaller than
.
Mixing of four massive neutrinos
Let us consider now the case of
mixing of four neutrinos [415,416,417],
i.e.,
in Eqs.(3.3) and (3.9),
 |
(6.56) |
which allows to accommodate all the three existing indications
in favour of neutrino mixing.
Here
is a sterile
neutrino state (see Section 2.2).
We do not make any assumptions
about the mass spectrum of four neutrinos and will consider all possibilities.
The six possible four-neutrino mass spectra
with three different scales of neutrino mass-squared differences
(atmospheric, solar and LSND)
are shown in Fig. 6.1.
In all these mass spectra there are two groups
of close masses separated by the ``LSND gap'' of the order of 1 eV.
The small mass-squared
differences correspond to
(the smallest one,
in schemes I and B,
in schemes II and IV,
in schemes III and A)
and
(
in schemes I and II,
in schemes III and IV,
in scheme A,
in scheme B)
and the largest mass squared difference
is relevant for the oscillations observed in the LSND experiment.
Let us consider SBL oscillations in the case of neutrino mass spectra
presented in Fig. 6.1.
Taking into account that
 |
(6.57) |
from the general formula (3.12)
we obtain that the probabilities of
transitions
and of
survival
in SBL experiments are given by
| |
|
 |
(6.58) |
| |
|
 |
(6.59) |
The oscillation amplitudes
and
depend on the elements of
the mixing matrix and on the type of neutrino mass spectrum:
| |
|
 |
(6.60) |
| |
|
 |
(6.61) |
where the index
runs over the indices
of the first or (because of the unitarity of the mixing matrix)
of the second group
of neutrino masses (see Fig. 6.1).
The expressions (6.58) and (6.59)
have the same form as the corresponding two-neutrino formulas
(3.20) and (3.21)
and are a direct
generalization of the expressions
(6.7) and (6.8)
presented in Section 6.1
for the case of three massive neutrinos and
a neutrino mass hierarchy.
From Eq.(6.60) it is obvious that the amplitude of
oscillations
is equal to
and
 |
(6.62) |
Thus,
CP violation in the lepton sector cannot be revealed in SBL
neutrino oscillation experiments if only one mass-squared difference is
relevant for SBL oscillations.
The expression (6.61)
allows to obtain information on the quantities
and
from the results of SBL
and
disappearance experiments,
whose exclusion curves imply the constraint
 |
(6.63) |
The upper bound
depends on the value of
.
Equations (6.59) and (6.63) imply that each
for
must satisfy one of the two inequalities
 |
(6.64) |
with
 |
(6.65) |
The most stringent limits on
and
are given by the exclusion curves of the
Bugey [162] reactor
experiment
and
of the CDHS [153] and CCFR [418]
accelerator
experiments.
The values of
and
corresponding to the 90% CL
exclusion curves of these experiments
are shown in Fig. 6.2.
One can see that
| |
|
 |
(6.66) |
| |
|
 |
(6.67) |
Hence,
from Eq.(6.64)
it follows that the quantities
and
are either small or large
(close to one).
We will show now that from the four possibilities for the two
quantities
and
(small-small, small-large, large-small, large-large)
for each neutrino mass spectrum in Fig. 6.1
only one possibility is compatible with the results of
solar and atmospheric neutrino experiments
[419,416,420].
Let us consider, for example,
the hierarchical mass spectrum
 |
(6.68) |
which corresponds to type I in Fig. 6.1.
In this case,
the inequalities (6.64) become
 |
(6.69) |
The survival probability of solar
's
averaged over the fast, unobservable oscillations due to
and
is given by [421]
 |
(6.70) |
where
is the
survival probability
due to the mixing of
with
and
.
This probability has the two-generation form
(3.21) and depends on the neutrino energy,
and
 |
(6.71) |
Though Eqs.(6.70) and (6.71) were derived here for
vacuum oscillations, the considerations in Section 4.4 show that they
are also valid including matter effects.
If
,
from Eqs.(6.66) and (6.70)
one can see that
practically does not depend on the neutrino energy
and
for all solar neutrino energies.
This is not compatible with the solar neutrino
data [329,349,422].
Therefore,
we reach the conclusion that from the two possibilities
(6.69) for
only
 |
(6.72) |
is compatible with the results of solar neutrino experiments.
Let us now consider the
survival probability of atmospheric
's.
From the general formula (3.12)
we obtain
 |
(6.73) |
where we have taken into account that for atmospheric
neutrinos27
.
The probability averaged over the fast unobservable oscillations
due to
is given by
 |
(6.74) |
Using the unitarity relation
 |
(6.75) |
the survival probability (6.74)
can be written as [416,420] (see also Section 4.4)
 |
(6.76) |
where
is the survival probability of atmospheric
's
due to the mixing of
with
and
,
.
This probability has the two-generation form
 |
(6.77) |
(see Eq.(3.21)),
with
 |
(6.78) |
The zenith angle dependence of
-like events
in atmospheric neutrino experiments
is due to
, which is multiplied by
in the expression (6.76) for
.
If
is large,
as given by the second inequality in Eq.(6.69),
,
we have
for
(see Eq.(6.67)).
A zenith-angle variation of
-like events smaller than
is incompatible with the up-down asymmetry
of multi-GeV
-like events observed by the Super-Kamiokande
experiment (see Section 5.1).
Therefore,
from the two inequalities (6.69) for
only the first one is compatible
with the results of the atmospheric neutrino experiments:
 |
(6.79) |
The same conclusion can be reached for all the neutrino mass spectra
in Fig. 6.1
in which one mass is separated by the other three
by the LSND gap,
i.e.,
for the spectra I, II, III, IV.
Hence,
in the framework of these spectra the results of SBL disappearance experiments
and of solar and atmospheric neutrino experiments
imply that
 |
(6.80) |
with
in the schemes I and II
and
in the schemes III and IV.
These inequalities imply that
transitions in SBL experiments
are strongly suppressed.
Indeed,
from Eqs.(6.60) and (6.80)
we have
 |
(6.81) |
Thus,
the amplitude of SBL
oscillations is quadratic in the small quantities
and
.
The upper bound (6.81) obtained from the values of
and
presented in Fig. 6.2
is shown in Fig. 6.3
by the solid line passing through the black circles.
The solid, dotted, dashed, dash-dotted and dash-dot-dotted
curves in Fig. 6.3
are the 90% CL
exclusion curves obtained in the
Bugey [162],
BNL E734 [372],
BNL E776 [373],
CCFR [374]
and
KARMEN (Bayesian analysis) [375]
experiments.
The shadowed regions are allowed by the results of the LSND
[370,371]
experiments at 90% CL.
From this figure one can see that the result of the LSND experiment,
taken together with the results of SBL disappearance experiments
and solar and atmospheric neutrino experiments,
disfavours the schemes I, II, III and IV of
Fig. 6.1.28
Let us consider now the spectra A and B in Fig. 6.1.
In this case the inequalities (6.64) apply
with
the index
running over
or
(the two choices are equivalent because of the unitarity of the mixing matrix).
We consider explicitly only the spectrum A,
but all the following results are valid also for scheme B
and can be obtained with the index exchange
.
In the case of spectrum A,
the survival probability of solar
's
is given by [421]
 |
(6.82) |
where
is the survival probability of solar
's
due to the mixing of
with
and
.
This probability has the two-generation form
(3.21) and depends on the neutrino energy,
and
 |
(6.83) |
If
,
i.e., the second inequality (6.64)
for
and
is satisfied,
from Eq.(6.66) we have
,
which implies that
is constant
and
.
Since these constraints are incompatible with the solar neutrino data,
we reach the conclusion that only the first inequality (6.64)
for
and
is compatible with the results of solar neutrino experiments,
i.e.,
.
Now we consider the survival probability of atmospheric
's.
In scheme A it is given by [421]
 |
(6.84) |
where the probability
has the two-generation form (3.21)
with
and
 |
(6.85) |
If
,
which corresponds to the second of the two inequalities (6.64)
for
and
,
the upper bound (6.67) for
implies that
is practically constant
and incompatible with the zenith-angle dependence of
multi-GeV
-like events
observed in the Super-Kamiokande experiment [2].
Hence,
only the first inequality (6.64)
for
and
is compatible with the results of
the Super-Kamiokande atmospheric neutrino experiments
and we have
.
Summarizing,
in the framework of scheme A
the results of SBL disappearance experiments
and of solar and atmospheric neutrino experiments
imply that
 |
(6.86) |
The corresponding bounds in scheme B
are obtained with the exchange
:
 |
(6.87) |
The constraints (6.86) and (6.87)
imply the following upper bound
for the amplitude of
transitions in SBL experiments:
 |
(6.88) |
(here the index
runs over
or over
).
This bound is linear in the small quantities
and
and turns out to be compatible with the results of the LSND experiment.
Actually,
this bound is worse than the unitarity bound
 |
(6.89) |
which causes the exclusion of the large
part of the LSND allowed region
(see Fig. 6.3, where
is represented by the Bugey curve).
We want to note that the bound
ceases to exist shortly below
eV
(see Fig. 6.2). On the other hand,
it follows from the LSND and Bugey experiments that
eV
(see Fig. 6.3). Thus, for
the small interval
our arguments leading to schemes A and B have to be
considered with caution and it would be desirable to have
experimental information on
in this interval.
Let us now discuss some consequences of
the schemes A and B,
taking into account the constraints (6.86) and (6.87).
First of all,
we will show that these bounds
lead to constraints on the neutrino mass measured
in tritium
-decay experiments
and on the effective Majorana mass that determines
the probability of neutrinoless double-
decay.
The electron spectrum in the decay
is given by
 |
(6.90) |
where
is the energy release,
and
are the electron momentum and energy,
,
is the Fermi function
and
is a constant.
Taking into account that
,
the first inequality in Eq.(6.86) implies that
in scheme A the effective neutrino mass
determined from the measurement of the high-energy part of the
-spectrum of
is
practically equal to the ``LSND mass''
:
 |
(6.91) |
If massive neutrinos are Majorana particles,
the matrix element of neutrinoless double-
decay
is proportional to the effective Majorana mass
[106,95] (see Appendix A)
 |
(6.92) |
In the case of scheme A,
the first inequality in Eq.(6.86) implies that
 |
(6.93) |
where
is the difference of the phases of
and
.
Since
is unknown,
from Eq.(6.93) we obtain the constraints
 |
(6.94) |
where the upper and lower bounds
correspond to
or
and
,
respectively.29Taking into account Eqs.(6.83) and (6.86),
we have
 |
(6.95) |
and Eq.(6.94)
can be written as
 |
(6.96) |
If
,
which corresponds to the MSW solution
with a small mixing angle,
we have
,
independent of the value of
and of the conservation of CP.
In the case of a large value of the parameter
,
which correspond to the MSW solution
with a large mixing angle
or to the vacuum oscillation solution,
in the future
it will be possible to obtain information
about
the violation of
CP in the lepton sector
if both the
-decay
and
neutrinoless double-
decay
experiments
will obtain positive results.
Indeed,
from the measurements of
,
and
with the help of Eq.(6.96)
it will be possible to determine the value of
(if
is large).
Let us emphasize that
in the case of CP conservation
the relative CP parities of
and
can be determined:
the CP parities are equal
if
and opposite if
.
Considering now scheme B,
the first inequality in Eq.(6.87) implies that
both
and
are strongly suppressed:
| |
|
 |
(6.97) |
| |
|
 |
(6.98) |
Therefore, only
if scheme A is realized in nature,
there is some possibility to measure the effect of the LSND mass
in future tritium decay experiments
and neutrinoless double-
decay experiments.
We discuss now the implications of the schemes A and B for
LBL neutrino oscillation experiments [423].
The smallness of
in scheme A
and of
in scheme B
implies that the electron neutrino has a
small mixing with the neutrinos whose mass-squared difference is
responsible for the oscillations of atmospheric neutrinos
(i.e.,
,
in scheme A and
,
in scheme B).
Hence,
as shown in the following part of this section,
the transition probability of
electron neutrinos and antineutrinos
into other states
in atmospheric and LBL experiments
is suppressed.
Taking into account that in scheme A we have
and
in LBL experiments,
we obtain from the general formula (3.12)
the following expression for the measurable
probability of
transitions
averaged over the fast oscillations due to
:
 |
(6.99) |
The corresponding probability in scheme B is given by the exchange
and the probability of
transitions
is obtained changing
.
Let us consider reactor LBL experiments
which measure the survival probability
.
We will show that the negative results of SBL reactor experiments
and
the results of solar neutrino experiments
imply strong constraints on the transition probability
of LBL reactor
's into other states.
Indeed,
from Eq.(6.99) and the corresponding one in scheme B,
for the survival probability of LBL reactor
's
we have the lower bounds
| |
|
 |
(6.100) |
| |
|
 |
(6.101) |
in schemes A and B,
respectively.
From the inequalities (6.86) and (6.87)
it follows that in both schemes
 |
(6.102) |
Since
is small
(see Fig. 6.2 and Eq.(6.66)),
the lower bound (6.102)
implies that the survival probability
is close to one.
For the transition probability of
's into any other state in LBL experiments,
 |
(6.103) |
we have the upper bound
 |
(6.104) |
The value of this bound,
which depends on
through the dependence of
on
(see Fig. 6.2),
is shown in Fig 6.4 (the solid line),
where it is compared with the upper bound for
obtained in the CHOOZ experiment [159]
(dash-dotted line)
and with the final sensitivity of the CHOOZ experiment
(dash-dot-dotted line).
The shadowed region in Fig 6.4
corresponds to the range of
allowed at 90% CL
by the results of the LSND and other experiments.
One can see from Fig 6.4
that if
lies in the LSND range,
the upper bound for
is much less than the sensitivity of the CHOOZ experiment
and it will be very difficult to measure
the disappearance of reactor
in future LBL experiments.
Let us consider now
transitions, which will be investigated in the near future
by the
K2K [271],
MINOS [293],
ICARUS [294]
and/or other [296,297,295,298]
experiments
(see [414]).
With the help of the
Cauchy-Schwarz inequality,
from Eq.(6.99) and Eq.(6.60)
we obtain
 |
(6.105) |
where the index
runs over
in the scheme A
and over
in the scheme B.
From the constraints (6.86)and (6.87),
in both schemes A and B we have
 |
(6.106) |
where
is the upper bound for the
oscillation amplitude
measured in SBL experiments.
Since the CPT relation (3.15)
and the conservation of probability, relation (6.103),
imply that
 |
(6.107) |
another upper limit for
can be obtained from the inequality (6.104):
 |
(6.108) |
Combining the upper bounds (6.106) and (6.108)
and taking into account that they are valid also for antineutrinos
[423],
we have
![\begin{displaymath}
P^{(\mathrm{LBL})}_{\stackrel{\makebox[0pt][l]
{$\hskip-3pt\...
...A_{\mu;e}^0
\, , \,
a_e^0 \left( 2 - a_e^0 \right)
\right)
\,.
\end{displaymath}](img1354.gif) |
(6.109) |
The curve corresponding
to this limit
obtained from the 90% CL
exclusion plots of the Bugey
[162]
experiment for
and
of the
BNL E734
[372],
BNL E776
[373],
CCFR
[374]
and
KARMEN (Bayesian analysis)
[375]
experiments
for
is represented by the dashed line
in Fig. 6.5.
For comparison,
the expected sensitivities
of the LBL accelerator neutrino experiments
K2K [271],
MINOS [293]
and
ICARUS [294] are also indicated
(the three straight vertical lines).
The dash-dotted and dash-dot-dotted
lines in Fig. 6.5 represent,
respectively,
the upper bound for
in the K2K and MINOS, ICARUS experiments
corrected for the matter effects [423] due to the
propagation of the neutrino beam in the crust of the earth.
These effects are different for the K2K and MINOS, ICARUS experiments
because of the different value of
in this experiments:
in K2K
and
in MINOS and ICARUS.
(It is obvious that the bound for other LBL experiments
with a neutrino beam from CERN to Gran Sasso
is the same as that for the ICARUS experiment.)
The region between the two horizontal solid lines in Fig. 6.5
corresponds to the range of
allowed at 90% CL
by the results of the LSND and other experiments.
One can see that if
is in the LSND region,
the sensitivity of the K2K experiment is not sufficient
for the observation of
transitions,
whereas experiments with higher sensitivities
like MINOS and ICARUS could reveal these transitions.
Finally,
another upper bound on
can be obtained from Eq.(6.106) noting that
 |
(6.110) |
with
in scheme A and
in scheme B.
Taking into account the inequalities (6.86) and (6.87),
we obtain
![\begin{displaymath}
P^{(\mathrm{LBL})}_{\stackrel{\makebox[0pt][l]
{$\hskip-3pt\...
...yle(-)$}}{\nu_{e}}}
\leq
a_e^0 + a_\mu^0 - 2 a_e^0 a_\mu^0
\,.
\end{displaymath}](img1361.gif) |
(6.111) |
This bound is stable against matter effects [423]
and its value is represented by the dotted curve in Fig. 6.5.
For
this bound is about half of that given by Eq.(6.108).
However,
since
is only small
in the same range of
where
is small,
numerically the bound
(6.111)
turns out to be worse than
the bound (6.109)
(dashed line in Fig. 6.5)
and
the corresponding matter-corrected bound
(dash-dotted and dash-dot-dotted lines in Fig. 6.5).
Summarizing the results presented in this section,
we have shown that only the two four-neutrino schemes A and B
are compatible with the results of all neutrino oscillation experiments.
If scheme A is realized in nature,
the effect of non-zero neutrino masses could be observed in the near future in
tritium
-decay experiments
and neutrinoless double-
decay experiments,
whereas these effects are strongly suppressed in scheme B.
Furthermore,
we have shown that
the results of SBL
disappearance experiments
and solar neutrino experiments imply that
in both schemes A and B the electron neutrino has a
small mixing with the neutrinos whose mass-squared difference is
responsible for the oscillations of atmospheric neutrinos
(i.e.,
,
in scheme A and
,
in scheme B).
As a consequence,
the transition probability of
electron neutrinos and antineutrinos
into other states
in atmospheric and LBL experiments
is strongly suppressed.
Furthermore,
it has been shown in Refs. [417,424]
that in the schemes A and B, if
the standard Big-Bang Nucleosynthesis constraint
on the number of light neutrinos
(see [425] and references therein) is less than 4,
then one has a stringent limit on the mixing of the sterile neutrino
with the two massive neutrinos
that are responsible for the oscillations
of atmospheric neutrinos.
In this case,
is mainly mixed with the
two massive neutrinos that contribute to solar neutrino oscillations
(
and
in scheme A
and
and
in scheme B)
and the unitarity of the mixing matrix implies that
is mainly mixed with the
two massive neutrinos that contribute to the oscillations
of atmospheric neutrinos.
Hence,
the two schemes have the form shown in Fig. 6.6
and have the following testable implications
for solar, atmospheric, long-baseline and
short-baseline neutrino oscillation experiments:
- The solar neutrino problem
is due to
oscillations.
This prediction will be checked by future solar neutrino
experiments
that can measure the ratio of
neutral-current and charged-current events
[356].
- The atmospheric neutrino anomaly is due to
oscillations.
This prediction will be investigated by LBL experiments.
Furthermore,
the absence of
atmospheric neutrino oscillations may be checked in the future
in the Super-Kamiokande atmospheric neutrino experiment
[271,272,275,276].
-
and
oscillations
are strongly suppressed in SBL experiments
[424].
If these prediction will be falsified by future experiments
it could mean that
some of the indications
given by the results of neutrino oscillations experiments
are wrong and neither of the two four neutrino schemes A and B
is realized in nature,
or that Big-Bang Nucleosynthesis occurs with
a non-standard mechanism
[426].
On the other hand,
if the standard Big-Bang Nucleosynthesis constraint
on the number of light neutrinos is less than 4 and one of
the two four neutrino schemes depicted in Fig. 6.6
is realized in nature,
at the leading order
the
neutrino mixing matrix has an extremely simple structure
in which the
and
sectors are decoupled:
in scheme A
 |
(6.112) |
and in scheme B
 |
(6.113) |
where
and
are,
respectively,
the two-generation mixing angles relevant in
solar and atmospheric neutrino oscillations.
Therefore,
the oscillations of solar
and atmospheric neutrinos are independent
and
the two-generation analyses of solar and atmospheric neutrino oscillations
yield correct information on the mixing of four neutrinos.
The
transitions observed in the LSND experiment represent only a slight
distortion of this simple picture.
Conclusions
In this review, we have discussed the phenomenological aspects
of neutrino oscillations. The phenomenon of neutrino oscillations
is possible if neutrinos are massive and mixed particles.
The problem of neutrino masses and mixing is the central problem of
today's neutrino physics and astrophysics.
More than 40 different experiments all over the
world are dedicated to the investigation of this problem
and several new
experiments are in preparation.
The investigation of the properties of neutrinos is considered as
one of the most important direction for the search of new physics
(see Ref. [31]).
Massive neutrinos are among the plausible candidates
for dark matter particles and
the number of massive light neutrinos plays
crucial role in Big-Bang Nucleosynthesis [425].
The Super-Kamiokande collaboration obtained recently a rather convincing
evidence in favour of transitions
of atmospheric
into other neutrino states (
or
).
Other indications
in favour of neutrino mixing were obtained in
all solar neutrino experiments and in the accelerator LSND experiment.
If all these indications will be confirmed by future experiments,
it will mean that the number of massive light neutrinos is larger
than the number of flavour neutrinos,
i.e. that sterile neutrinos exist.
This would be an evidence that
neutrino masses and mixing are phenomena
due to physics beyond the Standard Model with right-handed
neutrino singlets.
If the indications in favour of a ``large''
obtained in the LSND experiment will not be confirmed by future
experiments,
the mixing of three massive neutrinos with a mass
hierarchy is a possible and plausible scenario.
In this case,
the Dirac or Majorana nature of massive neutrinos
can be determined only increasing
the sensitivity of the experiments on the search for
neutrinoless double-
decay by almost two orders of magnitude.
It has been shown by an
analysis of the existing data [427,428]
that the effective Majorana
mass in the case of three massive Majorana neutrinos and a mass hierarchy is
smaller than
,
whereas the present upper bound is about
[429,430]
and the sensitivity of the next generation of experiments
will be of the order of
[431].
Future neutrino oscillation experiments will make possible to investigate in
detail the regions of
and of the elements of the neutrino mixing matrix
allowed by the existing results of neutrino oscillation experiments.
Future solar neutrino experiments
(SNO,
ICARUS,
Borexino,
GNO
and others
[353])
will provide model-independent information on
the transition of high energy
neutrinos into other states
and
investigate in detail the fluxes of medium energy
,
and CNO
neutrinos.
Future accelerator LBL experiments
(K2K [271],
MINOS [293],
ICARUS [294]
and others [296,297,295,298,414])
will explore
the region of
indicated by the atmospheric neutrino anomaly.
Besides
survival, also
and
transitions will be studied.
A very important task
for the neutrino oscillation experiments of the next generation
is the check of the indication in favour of neutrino oscillations
obtained in the LSND experiment.
The KARMEN experiment will reach the LSND sensitivity
in one or two years [379].
The proposed
SBL experiments
BooNE [387] at Fermilab,
I-216 [388] at CERN,
ORLaND [389] at Oak Ridge
and
NESS at the European Spallation Source [390]
plan to achieve a sensitivity about two
orders of magnitude better than the one of the LSND experiment.
It is now widely believed that the recent measurement
of an up-down asymmetry of atmospheric
-like events
in the Super-Kamiokande detector represents a strong evidence in
favour of neutrino oscillations
generated by new physics beyond the Standard Model
(see, for example, Ref. [432]).
We think that in order to understand the
origin of neutrino masses and mixing
it is necessary to investigate in much more detail
all the phenomena generated by these properties:
neutrino oscillations, neutrinoless double-
decay,
distortion of the high-energy part of the electron spectrum of
tritium
decay,
neutrino magnetic moments, etc.
There is no doubt that the next
generation of neutrino experiments will lead to a great
progress in the understanding the properties of neutrinos,
the most puzzling among the known particles.
After this review was finished and presented to the editor (December 15,
1998), new Super-Kamiokande data were reported at WIN99 (Capetown,
January 24 - 30, 1999)
by Y. Suzuki and at the VIIIth International Workshop on Neutrino
Telescopes (Venice, February 23 - 26, 1999) by K. Scholberg and
K. Inoue. Also Soudan-2 published an updated result.
Below we will shortly present the new development.
For the double ratio
(5.8) the new values based on 736 days
for FC events and 685 days for PC events are
For the up-down asymmetry (5.10) of
-like events in the
multi-GeV region the new value
is more than 7 standard deviations away from zero. Thus, the new data
confirm the Super-Kamiokande evidence in favour of oscillations of
atmospheric neutrinos. The updated dependence of the up-down asymmetry
as
a function of the
lepton momentum for
-like and
-like events is given in Fig. A.
The new allowed region in the plane of the oscillation parameters
under the assumption of
oscillations is shown in Fig. B.
The region allowed by Kamiokande and the allowed region found from the
analysis of previous Super-Kamiokande data
(see Fig. 5.1)
are also depicted.
As it is seen from Fig. B,
the new Super-Kamiokande data favour higher values of
than the previous ones.
The best-fit values of the parameters are
This fit corresponds to
for 67 DOF
(the previous best-fit value of
was
).
The new allowed range for
at 90% CL
is
.
The Super-Kamiokande Collaboration also measured
the zenith angle dependence of the number of up-going muon events
(for a discussion see Section 5.1.2).
Also these data are compatible with the hypothesis of
oscillations.
The best-fit values of the oscillation parameters
obtained from the data on the angular dependence
of through-going muon events are
and
(corresponding to
for 8 DOF).
The atmospheric neutrino FC and PC data can also be described if one
assumes that
transitions take place.
In this case the best-fit
values of the oscillation parameters are given by
and
eV
with
for 67 DOF.
The Super-Kamiokande Collaboration reported also the observation
of 231.8 NC
events (see Eq.(5.11)),
after subtraction of background,
with two
-like tracks with an invariant mass
in the range from 90 to 180 MeV.
As a preliminary result it was found that
In the case of
transitions
this ratio should be equal to one,
whereas in the case of
transitions
this ratio should be about 0.75.
The large systematic error in the present data does not allow to reach
a definite conclusion on the transition channel.
The updated value for the double ratio
of the Soudan-2 experiment is
(W.W.M. Allison et al.,
hep-ex/9901024).
In the region of recoil electron energies between 6.5 and 20 MeV,
during 708 days, the Super-Kamiokande Collaboration has observed
events (about 13.5 events per day). For the
value of the
B neutrino flux it was found
cm
s
. The ratio of this flux to the flux predicted
by the SSM [304] is equal to
.
Now Super-Kamiokande has also data with an energy
of the recoil electrons in the range from 5.5 to 6.5 MeV. This
lowering of the electron energy threshold is the
most important difference of the new data with respect to the
previous ones.
The best fit to the recoil energy
spectrum is given by vacuum oscillations with the best-fit parameters
and
eV
. This is mainly because of the observed excess of events with
respect to the SSM prediction in the high energy region of the
recoil electron spectrum.
The contribution to the recoil spectrum in the region
MeV comes from
neutrinos, which according to the SSM is small
(see Table 5.2). If the flux of
neutrinos is considered as a
free parameter, then from a fit of the Super-Kamiokande data under the
assumption of absence of neutrino oscillations it was found
. The contribution of
neutrinos in the
Super-Kamiokande recoil energy spectrum was analysed in detail in the
papers of J. Bahcall and P. Krastev, Phys. Lett. B 436, 243
(1998) and G. Fiorentini et al., Phys. Lett. B 444,
387 (1998).
Super-Kamiokande also presented a new measurement of the day-night
asymmetry with the result
.
S.M.B. acknowledges
the support of the ``Sonderforschungsbereich 375-95 für
Astro-Teilchenphysik der Deutschen Forschungsgemeinschaft''.
C.G. would like to express his gratitude to the Korea Institute
for Advanced Study (KIAS)
for the kind hospitality during part of this work.
Furthermore, we thank V.A. Naumov for reading a preliminary version of the
manuscript and providing several useful comments. We are particularly
grateful to K. Scholberg for providing us with the latest
Super-Kamiokande figures.
In this Appendix we summarize the main properties of Majorana
particles associated with Majorana fields. These
particles have spin 1/2 and all charges equal to zero.
A free Majorana field
satisfies the Dirac equation
 |
(A..1) |
together with the Majorana condition
 |
(A..2) |
where
is the charge-conjugate field and
is
the charge-conjugation matrix which fulfills the relations
 |
(A..3) |
From Eqs.(A.1) and (A.2)
it follows that the Fourier expansion of the Majorana field is
![\begin{displaymath}
\nu(x)
=
\int \frac{ \mathrm{d}^3p }{ \sqrt{(2\pi)^3 2 p_0} ...
... \, e^{-ipx}
+
a_r^\dagger(p) \, v^r(p) \, e^{ipx}
\right]
\,,
\end{displaymath}](img1404.gif) |
(A..4) |
where the spinors
and
are related by
 |
(A..5) |
The operators
and
annihilate and create
Majorana particles with momentum
and helicity
(particle
antiparticle for a Majorana field), respectively.
Let us consider the quantities
 |
(A..6) |
where
and
are two Majorana fields and
 |
(A..7) |
Taking into account the minus sign that
appears under interchange of two fermion operators
and the Majorana condition (A.2),
we obtain
 |
(A..8) |
From Eq.(A.3) we derive
 |
(A..9) |
With the help of this equation and Eq.(A.3) it follows that
 |
(A..10) |
with
 |
(A..11) |
Thus, for a single Majorana field
we obtain
 |
(A..12) |
This means that a Majorana particle has neither electric charge
nor magnetic moments or electric dipole moments.
For any fermion field there is a relation between its
left-handed component and the right-handed component of its
charge-conjugate field (see Section 2.2):
 |
(A..13) |
In the case of a Dirac field,
is different from
and, therefore, the left and right-handed components are
independent, whereas from Eq.(A.13) it follows for a Majorana
field
that
 |
(A..14) |
and thus the left and right-handed components of a Majorana field
are not independent
(in other words, a Majorana field is a two-component field).
Though Dirac and Majorana neutrinos differ in the number of
degrees of freedom, it is nevertheless difficult to
distinguish between the two in practice. The reason is that
in the standard weak interaction Lagrangians (2.1) and (2.2)
only the left-handed components
(
) of the flavour neutrino fields
are present. These are connected with the left-handed components
of the neutrino mass eigenfields by the relation
(see Section 2). The
mass eigenfields could be of Dirac or Majorana type. In the
first case we have
, where the
are independent degrees of freedom, whereas in the
second case we have
. The only
difference between the Dirac and Majorana case with respect to
the standard weak interactions resides in the
mass term. In the
limit of vanishing neutrino masses the right-handed Dirac fields
decouple from the standard weak interaction Lagrangians and the
distinction between Dirac and Majorana gets lost. Then the
physical fields are simply given by the
. In
summary, we have arrived at the following conclusion [433]:
With the
standard weak interactions one cannot distinguish between the
Dirac and Majorana nature of the neutrinos in the limit of
vanishing neutrino masses. In other words, for the
standard weak interactions, Dirac neutrinos with negative
helicity and Dirac antineutrinos with positive helicity are the
same as Majorana neutrinos with positive and negative
helicities, respectively, in this limit.
In neutrinoless double-
decay, which can occur only with
Majorana neutrinos (see Section 6), one needs to calculate the
propagator
 |
(A..15) |
With the Majorana mass eigenfields
satisfying the relation
this propagator is obtained
by the following steps:
Finally,
we want to prove that the CP parity of a Majorana field is
.
The CP transformation of a Majorana field
is given by
 |
(A..17) |
where the phase factor
is the CP parity of the field
,
is the unitary operator of CP conjugation,
and
.
The relation (A.17) can be written as
 |
(A..18) |
Using the Majorana condition (A.2)
one obtains
 |
(A..19) |
From the comparison of the expressions (A.17) and (A.19)
one finds that
,
which implies that
 |
(A..20) |
and
.
As it can be seen from Eq.(A.18) written as
 |
(A..21) |
in general the product of the CP parities of particles and antiparticles is
equal to
.
The relation (A.20) results from the fact that in the Majorana case
there is no difference between particles and antiparticles.
-
- 1
-
T. Kajita,
Super-Kamiokande Coll., Talk presented at Neutrino '98 [434],
1998.
- 2
-
Y. Fukuda et al.,
Super-Kamiokande Coll., Phys. Rev. Lett. 81, 1562 (1998).
- 3
-
W. Pauli,
in Neutrino Physics, p. 1, edited by K. Winter, Cambridge University
Press, 1991.
- 4
-
E. Fermi,
Ricerca Scient. 2, 12 (1933).
- 5
-
E. Fermi,
Z. Physik. 88, 161 (1934).
- 6
-
F. Perrin,
Comptes Rendues 197, 1625 (1933).
- 7
-
O. Kofoed-Hansen, Phys. Rev. 71, 451 (1947); G.C. Hanna and B.
Pontecorvo, Phys. Rev. 75, 451 (1949).
- 8
-
L.L. Langer and R.J.D. Moffat,
Phys. Rev. 88, 689 (1952).
- 9
-
L. Landau,
Nucl. Phys. 3, 127 (1957).
- 10
-
T.D. Lee and C.N. Yang,
Phys. Rev. 105, 1671 (1957).
- 11
-
A. Salam,
Il Nuovo Cim. 5, 299 (1957).
- 12
-
M. Goldhaber, L. Grodzins and A.W. Sunyar,
Phys. Rev. 109, 1015 (1958).
- 13
-
V.L. Telegdi,
History of Original Ideas and Basic Discoveries in Particle Physics,
edited by H.B. Newman and T. Ypsilantis, p. 421, Plenum Press, New York,
published in cooperation with NATO Scientific Affairs Division, 1996.
- 14
-
L. Grénacs et al.,
Am. J. Phys. 50, 931 (1982).
- 15
-
R.P. Feynman and M. Gell-Mann,
Phys. Rev. 109, 193 (1958).
- 16
-
E.C.G. Sudarshan and R. Marshak,
Phys. Rev. 109, 1860 (1958).
- 17
-
L.D. Landau, unpublished (1957), private communication of B. Pontecorvo to S.M.
Bilenky.
- 18
-
B. Pontecorvo,
J. Exptl. Theoret. Phys. 33, 549 (1957)
[Sov. Phys. JETP 6, 429 (1958)].
- 19
-
B. Pontecorvo,
J. Exptl. Theoret. Phys. 34, 247 (1958)
[Sov. Phys. JETP 7, 172 (1958)].
- 20
-
E. Majorana,
Il Nuovo Cim. 14, 170 (1937).
- 21
-
Z. Maki, M. Nakagawa and S. Sakata,
Prog. Theor. Phys. 28, 870 (1962).
- 22
-
B. Pontecorvo,
Zh. Eksp. Teor. Fiz. 53, 1717 (1967)
[Sov. Phys. JETP 26, 984 (1968)].
- 23
-
V. Gribov and B. Pontecorvo,
Phys. Lett. B 28, 493 (1969).
- 24
-
S.M. Bilenky and B. Pontecorvo,
Phys. Lett. B 61, 248 (1976).
- 25
-
H. Fritzsch and P. Minkowski,
Phys. Lett. B 62, 72 (1976).
- 26
-
S.M. Bilenky and B. Pontecorvo,
Lett. Nuovo Cim. 17, 569 (1976).
- 27
-
V. Barger, P. Langacker, J.P. Leveille and S. Pakvasa,
Phys. Rev. Lett. 45, 692 (1980).
- 28
-
S.M. Bilenky, J. Hošek and S.T. Petcov,
Phys. Lett. B 94, 495 (1980).
- 29
-
I.Yu. Kobzarev et al.,
Yad. Fiz. 32, 1590 (1980)
[Sov. J. Nucl. Phys. 32, 823 (1980)].
- 30
-
J. Schechter and J.W.F. Valle,
Phys. Rev. D 22, 2227 (1980).
- 31
-
R.N. Mohapatra and P.B. Pal,
Massive Neutrinos in Physics and Astrophysics, World Scientific
Lecture Notes in Physics, Vol. 41, World Scientific, Singapore, 1991.
- 32
-
M. Gell-Mann, P. Ramond and R. Slansky,
in Supergravity, p. 315, edited by F. van Nieuwenhuizen and D.
Freedman, North Holland, Amsterdam, 1979.
- 33
-
T. Yanagida,
Proc. of the Workshop on Unified Theory and the Baryon Number of the
Universe, KEK, Japan, 1979.
- 34
-
R.N. Mohapatra and G. Senjanovic,
Phys. Rev. Lett. 44, 912 (1980).
- 35
-
C. Caso et al.,
Particle Data Group, Eur. Phys. J. C 3, 1 (1998).
- 36
-
R. Davis, Harmer and K.C. Hoffman,
Phys. Rev. Lett. 21, 1205 (1968).
- 37
-
R. Davis,
Prog. Part. Nucl. Phys. 32, 13 (1994).
- 38
-
B.T. Cleveland et al.,
Astrophys. J. 496, 505 (1998).
- 39
-
K.S. Hirata et al.,
Kamiokande Coll., Phys. Rev. D 38, 448 (1988).
- 40
-
K.S. Hirata et al.,
Kamiokande Coll., Phys. Rev. D 44, 2241 (1991).
- 41
-
K.S. Hirata et al.,
Kamiokande Coll., Phys. Rev. Lett. 77, 1683 (1996).
- 42
-
P. Anselmann et al.,
GALLEX Coll., Phys. Lett. B 285, 376 (1992).
- 43
-
W. Hampel et al.,
GALLEX Coll., Phys. Lett. B 388, 384 (1996).
- 44
-
A.I. Abazov et al.,
SAGE Coll., Phys. Rev. Lett. 67, 3332 (1991).
- 45
-
D.N. Abdurashitov et al.,
SAGE Coll., Phys. Rev. Lett. 77, 4708 (1996).
- 46
-
K. Inoue,
Super-Kamiokande Coll., Talk presented at the
International
Workshop on Topics in Astroparticle and Underground Physics, Gran
Sasso, Italy, September, 1997
(http://www-sk.icrr.u-tokyo.ac.jp/doc/sk/pub/pub_sk.html).
- 47
-
Y. Fukuda et al.,
Super-Kamiokande Coll., Phys. Rev. Lett. 81, 1158 (1998).
- 48
-
Y. Suzuki,
Super-Kamiokande Coll., Talk presented at Neutrino '98 [434],
1998.
- 49
-
Y. Fukuda et al.,
Kamiokande Coll., Phys. Lett. B 335, 237 (1994).
- 50
-
R. Becker-Szendy et al.,
IMB Coll., Nucl. Phys. B (Proc. Suppl.) 38, 331 (1995).
- 51
-
W.W.M. Allison et al.,
Soudan-2 Coll., Phys. Lett. B 391, 491 (1997).
- 52
-
M. Ambrosio et al.,
MACRO Coll., Phys. Lett. B 434, 451 (1998).
- 53
-
C. Athanassopoulos et al.,
LSND Coll., Phys. Rev. Lett. 77, 3082 (1996).
- 54
-
C. Athanassopoulos et al.,
LSND Coll., Phys. Rev. Lett. 81, 1774 (1998).
- 55
-
V. Lobashov, Talk presented at Neutrino '98
[434].
- 56
-
C. Weinheimer,
Talk presented at Neutrino '98 [434], 1998.
- 57
-
E.W. Kolb and M.S. Turner,
The Early Universe, Addison-Wesley, 1990.
- 58
-
F. Boehm and P. Vogel,
Physics of Massive Neutrinos, Cambridge University Press, Cambridge,
1987.
- 59
-
B. Kayser, F. Gibrat-Debu and F. Perrier,
The Physics of Massive Neutrinos, World Scientific Lecture Notes in
Physics, Vol.25, World Scientific, Singapore, 1989.
- 60
-
J.N. Bahcall,
Neutrino Physics and Astrophysics, Cambridge University Press, 1989.
- 61
-
C.W. Kim and A. Pevsner,
Neutrinos in Physics and Astrophysics, Contemporary Concepts in
Physics, Vol. 8, Harwood Academic Press, Chur, Switzerland, 1993.
- 62
-
S.M. Bilenky and B. Pontecorvo,
Phys. Rep. 41, 225 (1978).
- 63
-
P.H. Frampton and P. Vogel,
Phys. Rep. 82, 339 (1982).
- 64
-
J.D. Vergados,
Phys. Rep. 133, 1 (1986).
- 65
-
S.M. Bilenky and S.T. Petcov,
Rev. Mod. Phys. 59, 671 (1987).
- 66
-
J.W.F. Valle,
Prog. Part. Nucl. Phys. 26, 91 (1991).
- 67
-
L. Oberauer and F. von Feilitzsch,
Rept. Prog. Phys. 55, 1093 (1992).
- 68
-
G. Gelmini and E. Roulet,
Rept. Prog. Phys. 58, 1207 (1995).
- 69
-
J. Brunner,
Fortsch. Phys. 45, 343 (1997).
- 70
-
K. Zuber,
Phys. Rep. 305, 295 (1998).
- 71
-
J.M. Conrad,
Talk presented at the 29
International Conference on
High-Energy Physics (ICHEP 98), Vancouver, Canada, 23-29 July, 1998
(hep-ex/9811009).
- 72
-
J.W.F. Valle,
hep-ph/9809234 (1998).
- 73
-
G. Raffelt,
hep-ph/9902271 (1999).
- 74
-
S.L. Glashow,
Nucl. Phys. 22, 597 (1961).
- 75
-
S. Weinberg,
Phys. Rev. Lett. 19, 1264 (1967).
- 76
-
A. Salam,
Proc. of the 8
Nobel Symposium on Elementary Particle
Theory, Relativistic Groups and Analyticity, edited by N. Svartholm, 1969.
- 77
-
R. Bertlmann and H. Pietschmann, Phys. Rev. D 15, 683 (1977); J. Ellis
et. al., Ann. Rev. Nucl. Part. Sci. 32, 443 (1982).
- 78
-
N. Cabibbo,
Phys. Rev. Lett. 10, 531 (1963).
- 79
-
M. Kobayashi and T. Maskawa,
Prog. Theor. Phys. 49, 652 (1973).
- 80
-
J. Schechter and J.W.F. Valle,
Phys. Rev. D 24, 1883 (1981)
[Erratum: D 25, 283 (1982)].
- 81
-
S.T. Petcov,
Yad. Fiz. 25, 641 (1977)
[Sov. J. Nucl. Phys. 25, 340, 698(E) (1977)].
- 82
-
T.P. Cheng and L.F. Li,
Phys. Rev. D 16, 1425 (1977).
- 83
-
B.W. Lee and R.E. Shrock,
Phys. Rev. D 16, 1444 (1977).
- 84
-
L.L. Chau and W.Y. Keung,
Phys. Rev. Lett. 53, 1802 (1984).
- 85
-
C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985); Z. Phys. C 29,
491 (1985); C. Jarlskog and R. Stora, Phys. Lett. B 208, 268 (1988);
C. Jarlskog, Proc. of CP Violation, edited by C. Jarlskog, Advanced
Series in High Energy Physics Vol. 3, p. 3, World Scientific, Singapore,
1989.
- 86
-
D.D. Wu, Phys. Rev. D 33, 860 (1986); I. Dunietz, O.W. Greenberg and
D.D. Wu, Phys. Rev. Lett. 55, 2935 (1985).
- 87
-
J.D. Bjorken and I. Dunietz, Phys. Rev. D 36, 2109 (1987); I. Dunietz,
Ann. Phys. 184, 350 (1988).
- 88
-
J. Schechter and J.W.F. Valle,
Phys. Rev. D 21, 309 (1980).
- 89
-
I. Schur,
Am. J. Math. 67, 472 (1945).
- 90
-
B. Zumino,
J. Math. Phys. 3, 1055 (1962).
- 91
-
G. Gelmini and M. Roncadelli,
Phys. Lett. B 99, 411 (1981).
- 92
-
H. Georgi et al.,
Nucl. Phys. B 193, 297 (1983).
- 93
-
J.E. Kim,
Phys. Rep. 150, 1 (1987).
- 94
-
W. Grimus and H. Neufeld,
Nucl. Phys. B 325, 18 (1989).
- 95
-
L. Wolfenstein,
Phys. Lett. B 107, 77 (1981).
- 96
-
S.M. Bilenky, N.P. Nedelcheva and S.T. Petcov,
Nucl. Phys. B 247, 61 (1984).
- 97
-
B. Kayser,
Phys. Rev. D 30, 1023 (1984).
- 98
-
S.M. Bilenky and S.T. Petcov,
Rev. Mod. Phys. 59, 671 (1987).
- 99
-
W. Konetschny and W. Kummer,
Phys. Lett. B 70, 433 (1977).
- 100
-
A. Zee,
Phys. Lett. B 93, 387 (1980).
- 101
-
A. Zee,
Phys. Lett. B 161, 141 (1985).
- 102
-
A. Zee,
Nucl. Phys. B 264, 99 (1986).
- 103
-
M.B. Voloshin,
Sov. J. Nucl. Phys. 48, 512 (1988).
- 104
-
K.S. Babu and R.N. Mohapatra,
Phys. Rev. D 43, 2278 (1991).
- 105
-
K.S. Babu,
Phys. Lett. B 203, 132 (1988).
- 106
-
M. Doi et al.,
Phys. Lett. B 102, 323 (1981).
- 107
-
P. Langacker et al.,
Nucl. Phys. B 282, 589 (1987).
- 108
-
J.C. Pati and A. Salam, Phys. Rev. D 10, 275 (1975); R.N. Mohapatra
and J.C. Pati, Phys. Rev. D 11, 566 and 2558 (1975).
- 109
-
R.N. Mohapatra and G. Senjanovic,
Phys. Rev. D 23, 165 (1981).
- 110
-
G. Ecker, W. Grimus and M. Gronau,
Nucl. Phys. B 279, 429 (1987).
- 111
-
G.C. Branco, W. Grimus and L. Lavoura,
Nucl. Phys. B 312, 492 (1989).
- 112
-
K. Kanaya,
Prog. Theor. Phys. 64, 2278 (1980).
- 113
-
J. Schechter and J.W.F. Valle,
Phys. Rev. D 25, 774 (1982).
- 114
-
D. Chang, R.N. Mohapatra and M.K. Parida, Phys. Rev. Lett. 52, 1072
(1984); D. Chang and R.N Mohapatra, Phys. Rev. D 32, 1248 (1985);
D.G. Lee and R.N. Mohapatra, Phys. Lett. B 329, 463 (1994).
- 115
-
S.A. Bludman, D.C. Kennedy and P.G. Langacker,
Phys. Rev. D 45, 1810 (1992).
- 116
-
S.A. Bludman, D.C. Kennedy and P.G. Langacker,
Nucl. Phys. B 374, 373 (1992).
- 117
-
E. Witten,
Phys. Lett. B 91, 81 (1980).
- 118
-
J.C. Pati and A. Salam,
Phys. Rev. Lett. 31, 661 (1973).
- 119
-
H. Georgi and S.L. Glashow,
Phys. Rev. Lett. 32, 438 (1974).
- 120
-
P. Langacker,
Phys. Rep. 72, 185 (1981).
- 121
-
F. Wilczek and A. Zee,
Phys. Rev. Lett. 43, 1571 (1979).
- 122
-
S. Weinberg,
Phys. Rev. Lett. 43, 1566 (1979).
- 123
-
S. Weinberg,
Phys. Rev. D 22, 1694 (1980).
- 124
-
H.A. Weldon and A. Zee,
Nucl. Phys. B 173, 269 (1980).
- 125
-
S. Weinberg,
Int. J. Mod. Phys. A 2, 301 (1987).
- 126
-
E.Kh. Akhmedov, Z.G. Berezhiani and G. Senjanovic,
Phys. Rev. Lett. 69, 3013 (1992).
- 127
-
E.Kh. Akhmedov, Z.G. Berezhiani, G. Senjanovic and Z. Tao,
Phys. Rev. D 47, 3245 (1993).
- 128
-
L. Wolfenstein,
Nucl. Phys. B 186, 147 (1981).
- 129
-
M. Kobayashi, C.S. Lim and M.M. Nojiri,
Phys. Rev. Lett. 67, 1685 (1991).
- 130
-
C. Giunti, C.W. Kim and U.W. Lee,
Phys. Rev. D 46, 3034 (1992).
- 131
-
S. Nussinov,
Phys. Lett. B 63, 201 (1976).
- 132
-
B. Kayser,
Phys. Rev. D 24, 110 (1981).
- 133
-
C. Giunti, C.W. Kim and U.W. Lee,
Phys. Rev. D 44, 3635 (1991).
- 134
-
C. Giunti, C.W. Kim and U.W. Lee,
Phys. Rev. D 45, 2414 (1992).
- 135
-
C. Giunti, C.W. Kim and U.W. Lee,
Phys. Lett. B 274, 87 (1992).
- 136
-
C. Giunti, C.W. Kim, J.A. Lee and U.W. Lee,
Phys. Rev. D 48, 4310 (1993).
- 137
-
J. Rich,
Phys. Rev. D 48, 4318 (1993).
- 138
-
K. Kiers, S. Nussinov and N. Weiss,
Phys. Rev. D 53, 537 (1996).
- 139
-
W. Grimus and P. Stockinger,
Phys. Rev. D 54, 3414 (1996).
- 140
-
J.E. Campagne,
Phys. Lett. B 400, 135 (1997).
- 141
-
Y. Grossman and H.J. Lipkin,
Phys. Rev. D 55, 2760 (1997).
- 142
-
K. Kiers and N. Weiss,
Phys. Rev. D 57, 3091 (1998).
- 143
-
C. Giunti, C.W. Kim and U.W. Lee,
Phys. Lett. B 421, 237 (1997).
- 144
-
C. Giunti and C.W. Kim,
Phys. Rev. D 58, 017301 (1998).
- 145
-
W. Grimus, P. Stockinger and S. Mohanty,
Phys. Rev. D 59, 013011 (1998).
- 146
-
L. Stodolsky,
Phys. Rev. D 58, 036006 (1998).
- 147
-
E.J. Chun, C.W. Kim and U.W. Lee,
Phys. Rev. D 58, 093003 (1998).
- 148
-
Y. Koide,
Phys. Rev. D 57, 5836 (1998).
- 149
-
Y. Koide and H. Fusaoka,
Phys. Rev. D 59, 053004 (1998).
- 150
-
S.M. Bilenky and C. Giunti,
Phys. Lett. B 300, 137 (1993).
- 151
-
N. Cabibbo,
Phys. Lett. B 72, 333 (1978).
- 152
-
D. Orestano,
NOMAD Coll., Talk presented at the CAPP 98 Workshop, CERN, June, 1998
[435].
- 153
-
F. Dydak et al.,
CDHS Coll., Phys. Lett. B 134, 281 (1984).
- 154
-
N. Ushida et al.,
Phys. Rev. Lett. 57, 2897 (1986).
- 155
-
P. Vilain et al.,
CHARM II Coll., Z. Phys. C 64, 539 (1994).
- 156
-
K.S. McFarland et al.,
CCFR Coll., Phys. Rev. Lett. 75, 3993 (1995).
- 157
-
CHORUS Coll., Talk presented at the 29
International
Conference on High-Energy Physics (ICHEP 98), Vancouver, Canada, 23-29 July
1998 (preprint hep-ex/9807024).
- 158
-
J. Altegoer et al.,
NOMAD Coll., Phys. Lett. B 431, 219 (1998).
- 159
-
M. Apollonio et al.,
CHOOZ Coll., Phys. Lett. B 420, 397 (1998).
- 160
-
G. Zacek et al.,
Phys. Rev. D 34, 2621 (1986).
- 161
-
G.S. Vidyakin et al.,
JETP Lett. 59, 390 (1994).
- 162
-
B. Achkar et al.,
Bugey Coll., Nucl. Phys. B 434, 503 (1995).
- 163
-
L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); ibid.
20, 2634 (1979).
- 164
-
S.P. Mikheyev and A.Yu. Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J.
Nucl. Phys. 42, 913 (1985)]; Il Nuovo Cim. C 9, 17 (1986);
Zh. Eksp. Teor. Fiz. 91, 7 (1986) [Sov. Phys. JETP 64, 4
(1986)].
- 165
-
V. Barger, K. Whisnant, S. Pakvasa and R.J.N. Phillips,
Phys. Rev. D 22, 2718 (1980).
- 166
-
S.P. Rosen and J.M. Gelb,
Phys. Rev. D 34, 969 (1986).
- 167
-
H.A. Bethe,
Phys. Rev. Lett. 56, 1305 (1986).
- 168
-
A. Halprin, Phys. Rev. D 34, 3462 (1986); P.D. Mannheim, Phys. Rev. D
37, 1935 (1988).
- 169
-
W. Grimus and T. Scharnagl,
Mod. Phys. Lett. A 8, 1943 (1993).
- 170
-
P.B. Pal and T.N. Pham,
Phys. Rev. D 40, 259 (1989).
- 171
-
J. Liu,
Phys. Rev. D 45, 1428 (1992).
- 172
-
P.B. Pal,
Int. J. Mod. Phys. A 22, 5387 (1992).
- 173
-
C. Giunti, C.W. Kim and W.P. Lam,
Phys. Rev. D 43, 164 (1991).
- 174
-
J.N. Bahcall and R. Ulrich,
Rev. Mod. Phys. 60, 297 (1988).
- 175
-
J.N. Bahcall and M.H. Pinsonneault,
Rev. Mod. Phys. 67, 781 (1995).
- 176
-
L.L. Foldy and S.A. Wouthuysen,
Phys. Rev. 78, 29 (1950).
- 177
-
P.G. Langacker et al.,
Nucl. Phys. B 282, 589 (1987).
- 178
-
S.P. Mikheyev and A.Yu. Smirnov,
Prog. Part. Nucl. Phys. 23, 41 (1989).
- 179
-
L. Stodolsky,
Phys. Rev. D 36, 2273 (1987).
- 180
-
H. Nunokawa, A. Rossi, V.B. Semikoz and J.W.F. Valle, Nucl. Phys. B
472, 495 (1996); A.B. Balantekin, J.M. Fetter and F.N. Loreti, Phys.
Rev. D 54, 3941 (1996); C.P. Burgess and D. Michaud, Ann. Phys.
256, 1 (1997).
- 181
-
P. Bamert, C.P. Burgess and D. Michaud,
Nucl. Phys. B 513, 319 (1998).
- 182
-
A.J. Baltz and J. Weneser,
Phys. Rev. D 37, 3364 (1988).
- 183
-
T.K. Kuo and J. Pantaleone,
Rev. Mod. Phys. 61, 937 (1989).
- 184
-
W.C. Haxton,
Ann. Rev. Astron. Astrophys. 33, 459 (1995).
- 185
-
L.I. Schiff,
Quantum Mechanics, McGaw-Hill, Singapore, 1987.
- 186
-
V. A. Naumov, Zh. Eksp. Teor. Fiz. 101, 3 (1992) [Sov. Phys. JETP.
74, 1 (1992)]; Int. J. Mod. Phys. D 1, 379 (1992); Phys.
Lett. B323, 351 (1994).
- 187
-
S.J. Parke,
Phys. Rev. Lett. 57, 1275 (1986).
- 188
-
W.C. Haxton, Phys. Rev. Lett. 57, 1271 (1986); Phys. Rev. D
35, 2352 (1987).
- 189
-
T.K. Kuo and J. Pantaleone,
Phys. Rev. D 39, 1930 (1989).
- 190
-
M. Born and V. Fock,
Z. Phys. 51, 165 (1928).
- 191
-
E. Lisi and D. Montanino,
Phys. Rev. D 56, 1792 (1997).
- 192
-
C. Giunti, C.W. Kim and M. Monteno,
Nucl. Phys. B 521, 3 (1998).
- 193
-
S.T. Petcov, Phys. Lett. B 200, 373 (1988); Phys. Lett. B
214, 139 (1988).
- 194
-
S.T. Petcov and J. Rich,
Phys. Lett. B 224, 426 (1989).
- 195
-
S.T. Petcov,
Phys. Lett. B 191, 299 (1987).
- 196
-
C. Zener,
Proc. Roy. Soc. Lon. A 137, 696 (1932).
- 197
-
L. Landau,
Phys. Z. Sov. 2, 46 (1932).
- 198
-
A. Joye and Ch.-E. Pfister,
preprint math-ph/9807031 (1998).
- 199
-
P. Pizzochero,
Phys. Rev. D 36, 2293 (1987).
- 200
-
A.B. Balantekin and J.F. Beacom, Phys. Rev. D 54, 6323 (1996); A.B.
Balantekin, J.F. Beacom and J.B. Fetter, Phys. Lett. B 427, 317
(1998).
- 201
-
T. Kaneko, Prog. Theor. Phys. 78, 532 (1987); S. Toshev, Phys. Lett. B
196, 170 (1987); M. Ito, T. Kaneko and M. Nakagawa, Prog. Theor.
Phys. 79, 13 (1988).
- 202
-
D. Nötzold,
Phys. Rev. D 36, 1625 (1987).
- 203
-
P.I. Krastev and S.T. Petcov,
Phys. Lett. B 207, 64 (1988)
[Erratum: B 214, 661 (1988)].
- 204
-
V.F. Hess,
Phys. Zeitschr. 13, 1084 (1913).
- 205
-
A. Pais,
Inward bound, Oxford University Press, Oxford, 1994.
- 206
-
L.V. Volkova,
Jad. Fiz. 31, 784 (1980)
[Sov. J. Nucl. Phys. 31, 1510 (1980)].
- 207
-
K. Mitsui, Y. Minorikawa and H. Komori,
Il Nuovo Cim. C 9, 995 (1986).
- 208
-
E.V. Bugaev and V.A. Naumov,
Izv. Akad. Nauk SSSR, Ser. Fiz. 50, 2239 (1986)
[Bull. Acad. of Sci. of the USSR, Phys. Ser. 50, 156 (1986)].
- 209
-
E.V. Bugaev and V.A. Naumov,
Yad. Fiz. 45, 1380 (1987)
[Sov. J. Nucl. Phys. 45, 857 (1987)].
- 210
-
A.V. Butkevich, L.G. Dedenko and I.M. Zeleznykh,
Yad. Fiz. 50, 142 (1989)
[Sov. J. Nucl. Phys. 50, 90 (1989)].
- 211
-
G. Barr, T.K. Gaisser and T. Stanev,
Phys. Rev. D 39, 3532 (1989).
- 212
-
H. Lee and Y.S. Koh,
Il Nuovo Cim. 105 B, 883 (1989).
- 213
-
E.V. Bugaev and V.A. Naumov,
Phys. Lett. B 232, 391 (1989).
- 214
-
M. Honda, K. Kasahara, K. Hidaka and S. Midorikawa,
Phys. Lett. B 248, 193 (1990).
- 215
-
M. Honda, T. Kajita, K. Kasahara and S. Midorikawa,
Phys. Rev. D 52, 4985 (1995).
- 216
-
V. Agrawal, T.K. Gaisser, P. Lipari and T. Stanev,
Phys. Rev. D 53, 1314 (1996).
- 217
-
T.K. Gaisser, M. Honda, K. Kasahara, S. Midorikawa, V. Naumov and T. Stanev,
Phys. Rev. D 54, 5578 (1996).
- 218
-
P. Lipari, T. Stanev and T.K. Gaisser,
Phys. Rev. D 58, 073003 (1998).
- 219
-
T.K. Gaisser,
Cosmic Rays and Particle Physics, Cambridge University Press,
Cambridge, 1990.
- 220
-
V.A. Dogiel, V.S. Berezinsky, S.V. Bulanov and V.S. Ptuskin,
Astrophysics of Cosmic Rays, edited by V.L. Ginzburg, North-Holland,
Amsterdam, 1990.
- 221
-
P. Lipari,
Astropart. Phys. 1, 195 (1993).
- 222
-
A. Misaki, V.A. Naumov, T.S. Sinegovskaya and S.I. Sinegovsky, Proc. of the
25th International Cosmic Ray Conference, Durban, South Africa,
edited by M. S. Potgieter et al., p. 129, 1997.
- 223
-
V.A. Naumov, T.S. Sinegovskaya and S.I. Sinegovsky,
Il Nuovo Cim. A 111, 129 (1998).
- 224
-
J.A. Simpson,
Ann. Rev. Nucl. Part. Sci. 33, 323 (1983).
- 225
-
W.R. Webber and J.A. Lezniak,
Astrophys. Space Sci. 30, 361 (1974).
- 226
-
E.-S. Seo et al.,
Proc. of the 22nd International Cosmic Ray Conference, Dublin,
Ireland, edited by M. Cawley et. al., Vol. 2, p. 627, Reprint Ltd.,
Dublin, Ireland, 1991.
- 227
-
E.V. Bugaev and V.A. Naumov, Proc. of the Second International Symposium
Underground Physics'87, Baksan Valley, August 1987, edited by G. V.
Domogatsky et al., Nauka, Moscow, p. 255, 1988.
- 228
-
P. Lipari and T. Stanev,
Proc. of the 24th International Cosmic Ray Conference, Rome, Vol. 1,
p. 694, 1995.
- 229
-
H. Lee and S.A. Bludman,
Phys. Rev. D 37, 122 (1988).
- 230
-
J.L. Osborne, S.S. Said and A.E. Wolfendale,
Proc. Phys. Soc. 86, 93 (1965).
- 231
-
T.K. Gaisser,
Proc. of Neutrino '96, Helsinki, June 1996, edited by K. Enqvist
et al., p. 211, World Scientific, Singapore, 1997.
- 232
-
T.K. Gaisser and T. Stanev,
Phys. Rev. D 57, 1977 (1998).
- 233
-
T.K. Gaisser,
Phil. Trans. Roy. Soc. Lon. A 346, 75 (1994).
- 234
-
M.A. Markov,
Proc. of 1960 Annual International Conference on High Energy Physics,
Rochester, edited by E.C.G. Sudarshan et al., 1960.
- 235
-
K. Greisen,
Ann. Rev. Nucl. Sci. 10, 63 (1960).
- 236
-
C.V. Achar et al.,
Phys. Lett. 18, 196 and 19, 78 (1965).
- 237
-
F. Reines et al.,
Phys. Rev. Lett. 15, 429 (1965).
- 238
-
M.A. Markov and I.M. Zhelenykh,
Nucl. Phys. 27, 385 (1961).
- 239
-
G.T. Zatsepin and V.A. Kuzmin,
Sov. J. Nucl. Phys. 14, 1294 (1962).
- 240
-
R. Cowsik, Y. Pal and S.N. Tandon,
Proc. Ind. Acad. Sci. 63, 217 (1966).
- 241
-
E.W. Beier and E.D. Frank,
Phil. Trans. Roy. Soc. Lon. A 346, 63 (1994).
- 242
-
K.S. Hirata et al.,
Kamiokande Coll., Phys. Lett. B 205, 416 (1988).
- 243
-
K.S. Hirata et al.,
Kamiokande Coll., Phys. Lett. B 280, 146 (1992).
- 244
-
R. Becker-Szendy et al.,
IMB Coll., Phys. Rev. Lett. 69, 1010 (1992).
- 245
-
R. Clark et al.,
IMB Coll., Phys. Rev. Lett. 79, 345 (1997).
- 246
-
K. Daum et al.,
Fréjus Coll., Z. Phys. C 66, 417 (1995).
- 247
-
M. Aglietta et al.,
NUSEX Coll., Europhys. Lett. 8, 611 (1989).
- 248
-
T.K. Gaisser, F. Halzen and T. Stanev,
Phys. Rep. 258, 174 (1995).
- 249
-
W. Frati et al.,
Phys. Rev. D 48, 1140 (1993).
- 250
-
M.M. Boliev et al., Proc. of the 24th International Cosmic Ray
Conference, ICRC '95, Rome, Vol. 1, p. 686, Arti Grafiche Editoriali,
Urbino, Italy, 1995; ibidem, p. 722; S.P. Mikheyev, Proc. of the
8
Rencontres de Blois on Neutrinos, Dark Matter, and
the Universe, Blois, France, 1996; M.M. Boliev et al., Proc. of the
International Workshop on Aspects of Dark Matter in Astrophysics and
Particle Physics, DARK '96, Heidelberg, Germany, 1996.
- 251
-
Y. Totsuka et al.,
Kamiokande Coll., Nucl. Phys. B 31 (Proc. Suppl.), 428 (1993).
- 252
-
S. Hatakeyama et al.,
Kamiokande Coll., Phys. Rev. Lett. 81, 2016 (1998).
- 253
-
P. Lipari, M. Lusignoli and F. Sartogo,
Phys. Rev. Lett. 74, 4384 (1995).
- 254
-
T.J. Haines et al.,
IMB Coll., Phys. Rev. Lett. 57, 1986 (1986).
- 255
-
E. Kearns,
Nucl. Phys. B (Proc. Suppl.) 70, 315 (1999)
(hep-ex/9803007).
- 256
-
T. Kafka,
Nucl. Phys. B (Proc. Suppl.) 70, 340 (1999)
(hep-ex/9712281).
- 257
-
M.C. Gonzalez-Garcia et al.,
Phys. Rev. D 58, 033004 (1998).
- 258
-
M.C. Gonzalez-Garcia et al.,
preprint FTUV/98-56, IFIC/98-57, hep-ph/9807305 (1998).
- 259
-
R. Foot, R.R. Volkas and O. Yasuda,
Phys. Rev. D 58, 013006 (1998).
- 260
-
O. Yasuda,
Phys. Rev. D 58, 091301 (1998).
- 261
-
G.L. Fogli, E. Lisi, A. Marrone and G. Scioscia,
Phys. Rev. D 59, 033001 (1999).
- 262
-
G.L. Fogli, E. Lisi and D. Montanino,
Phys. Rev. D 49, 3626 (1994).
- 263
-
G.L. Fogli, E. Lisi and D. Montanino,
Astropart. Phys. 4 177 (1995).
- 264
-
G.L. Fogli and E. Lisi,
Phys. Rev. D 52, 2775 (1995).
- 265
-
G.L. Fogli, E. Lisi, D. Montanino and G. Scioscia,
Phys. Rev. D 55, 4385 (1997).
- 266
-
G.L. Fogli, E. Lisi, A. Marrone and D. Montanino,
Phys. Lett. B 425, 341 (1998).
- 267
-
J.N. Flanagan, J.G. Learned and S. Pakvasa,
Phys. Rev. D 57, R2649 (1998).
- 268
-
J. Bunn, R. Foot and R. Volkas,
Phys. Lett. B 413, 109 (1997).
- 269
-
R. Foot, R. Volkas and O. Yasuda,
Phys. Rev. D 58, 013006 (1998).
- 270
-
R. Foot, C.N. Leung and O. Yasuda,
Phys. Lett. B 443, 185 (1998).
- 271
-
Y. Suzuki,
Proc. of Neutrino 96, Helsinki, June 1996, edited by K. Enqvist
et. al., p. 237, World Scientific, Singapore, 1997.
- 272
-
F. Vissani and A.Yu. Smirnov,
Phys. Lett. B 432, 376 (1998).
- 273
-
A. McDonald,
SNO Coll., Talk presented at Neutrino '98 [434], 1998.
- 274
-
Y. Fukuda et al.,
Kamiokande Coll., Phys. Lett. B 388, 397 (1996).
- 275
-
K.R.S. Balaji, G. Rajasekaran and S. Uma Sankar,
preprint hep-ph/9806256 (1998).
- 276
-
F. Vissani,
Talk presented at the Ringberg Euroconference New Trends in Neutrino
Physics, Tegernsee, Germany, May 1998, preprint DESY-98-135, 1998.
- 277
-
J.G. Learned, S. Pakvasa and J.L. Stone,
Phys. Lett. B 435, 131 (1998).
- 278
-
E. Akhmedov, P. Lipari and M. Lusignoli,
Phys. Lett. B 300, 128 (1993).
- 279
-
G.L. Fogli, E. Lisi and A. Marrone,
Phys. Rev. D 57, 5893 (1998).
- 280
-
Q.Y. Liu and A.Yu. Smirnov,
Nucl. Phys. B 524, 505 (1998).
- 281
-
V.K. Ermilova, V.A. Tsarev and V.A. Chechin, Kr. Soob., Fiz. [Short Notices of
the Lebedev Institute] 5, 26 (1986); E.Kh. Akhmedov, Yad. Fiz.
47, 475 (1988) [Sov. J. Nucl. Phys. 47, 301 (1988)].
- 282
-
Q.Y. Liu, S.P. Mikheyev and A.Yu. Smirnov,
hep-ph/9803415 (1998).
- 283
-
S.T. Petcov,
Phys. Lett. B 434, 321 (1998).
- 284
-
E.Kh. Akhmedov,
Nucl. Phys. B 538, 25 (1999).
- 285
-
M. Chizhov, M. Maris and S.T. Petcov,
hep-ph/9810501 (1998).
- 286
-
E.Kh. Akhmedov, A. Dighe, P. Lipari and A.Yu. Smirnov,
hep-ph9808270 (1998).
- 287
-
P. Lipari and M. Lusignoli,
Phys. Rev. D 57, 3842 (1998).
- 288
-
F. Boehm et al., The Palo Verde experiment, 1996
(http://www.cco.caltech.edu/~songhoon/Palo-Verde/Palo-Verde.html).
- 289
-
F. Boehm et al.,
Palo Verde Coll., Nucl. Phys. B (Proc. Suppl.) 70, 191 (1999).
- 290
-
F. Suekane,
preprint TOHOKU-HEP-97-02 (1997).
- 291
-
A. Suzuki,
Talk presented at NOW'98 [376], 1998.
- 292
-
S. Schönert et al.,
Nucl. Phys. B (Proc. Suppl.) 70, 195 (1999).
- 293
-
D. Ayres et al.,
MINOS Coll., report NUMI-L-63 (1995).
- 294
-
P. Cennini et al.,
ICARUS Coll., report LNGS-94/99-I (1994).
- 295
-
H. Shibuya et al.,
OPERA Coll., report CERN-SPSC-97-24 (1997).
- 296
-
M. Ambrosio et al., NOE Coll., Nucl. Instr. Meth. A 363, 604
(1995); G.C. Barbarino et al., NOE Coll., report INFN/AE-96/11
(1996).
- 297
-
T. Ypsilantis et al.,
AQUA-RICH Coll., report CERN-LAA/96-13 (1996).
- 298
-
G. Baldini et al.,
NICE Coll., report LNGS-LOI 98/13 (1998).
- 299
-
P. Picchi and F. Pietropaolo,
Talk presented at Neutrino '98 [434] (preprint hep-ph/9812222)
(1998).
- 300
-
M. Aglietta et al.,
preprint CERN-SPSC-98-28 (1998).
- 301
-
G. Battistoni and P. Lipari,
to appear in the Proc. of the 1998 Vulcano Workshop on Frontier Objects
in Astrophysics and Particle Physics, 1998
(hep-ph/9807475).
- 302
-
A. Geiser,
Contribution to NOW'98 [376], 1998
(hep-ph/9810493).
- 303
-
J.N. Bahcall and M.H. Pinsonneault,
Rev. Mod. Phys. 64, 885 (1992).
- 304
-
J.N. Bahcall, S. Basu and M.H. Pinsonneault,
Phys. Lett. B 433, 1 (1998).
- 305
-
S. Turck-Chièze, S. Cahen, M. Cassé and C. Doom,
Astrophys. J. 335, 415 (1988).
- 306
-
S. Turck-Chièze and I. Lopes,
Astrophys. J. 408, 347 (1993).
- 307
-
S. Turck-Chièze et al.,
Phys. Rep. 230, 57 (1993).
- 308
-
V. Castellani, S. Degl'Innocenti and G. Fiorentini,
Astron. Astrophys. 271, 601 (1993).
- 309
-
V. Castellani et al.,
Phys. Lett. B 324, 425 (1994).
- 310
-
F. Ciacio, S. Degl'Innocenti and B. Ricci et al.,
Astron. Astrophys. Suppl. Ser. 123, 449 (1997).
- 311
-
M. Koshiba,
Phys. Rep. 220, 229 (1992).
- 312
-
V. Castellani et al.,
Phys. Rep. 281, 309 (1997).
- 313
-
V.S. Berezinsky,
Invited lecture at the
International Cosmic Ray
Conference, Durban, 28 July - 8 August, 1997
(astro-ph/9710126).
- 314
-
J.N. Bahcall,
Lectures presented at the XXV SLAC Summer Institute on Particle Physics,
Physics of Leptons, August 4-15, 1997
(hep-ph/9711358).
- 315
-
J.N. Bahcall,
Talk presented at Neutrino '98 [434], 1998
(hep-ph/9808162).
- 316
-
A.Yu. Smirnov,
Talk presented at Neutrino '98 [434], 1998
(hep-ph/9809481).
- 317
-
J.N. Bahcall et al.,
Phys. Rev. C 54, 411 (1996).
- 318
-
J.N. Bahcall,
Phys. Rev. C 56, 3391 (1997).
- 319
-
J.N. Bahcall,
Phys. Rev. D 49, 3923 (1994).
- 320
-
J.N. Bahcall,
Phys. Rev. D 44, 1644 (1991).
- 321
-
J.N. Bahcall and P.I. Krastev,
Phys. Rev. D 53, 4211 (1996).
- 322
-
T. Kirsten,
GNO Coll., Talk presented at Neutrino '98 [434], 1998.
- 323
-
V.N. Gavrin,
SAGE Coll., Talk presented at Neutrino '98 [434], 1998.
- 324
-
B. Pontecorvo,
Chalk River Report PD-205 (1946).
- 325
-
L.W. Alvarez,
University of California Radiation Laboratory Report UCRL-328 (1949).
- 326
-
W. Hampel et al.,
GALLEX Coll., Phys. Lett. B 420, 114 (1998).
- 327
-
W. Hampel et al.,
GALLEX Coll., preprint DAPNIA-SPP-98-07 (1998).
- 328
-
M. Cribier et al.,
GALLEX Coll., Nucl. Phys. B (Proc. Suppl.) 70, 284 (1999).
- 329
-
J.N. Bahcall, P.I. Krastev and A.Yu. Smirnov,
Phys. Rev. D 58, 096016 (1998).
- 330
-
J.N. Bahcall, M.H. Pinsonneault, S.Basu and J. Christensen-Dalsgaard,
Phys. Rev. Lett. 78, 171 (1997).
- 331
-
B. Ricci et al.,
Phys. Lett. B 407, 155 (1997).
- 332
-
D. Gough,
Talk presented at Neutrino '98 [434], 1998.
- 333
-
A.S. Brun, S. Turck-Chièze and P Morel,
preprint astro-ph/9806272 (1998).
- 334
-
M. Spiro and D. Vignaud,
Phys. Lett. B 242, 279 (1990).
- 335
-
V. Castellani, S. Degl'Innocenti and G. Fiorentini,
Phys. Lett. B 303, 68 (1993).
- 336
-
N. Hata, S.A. Bludman and P. Langacker et al.,
Phys. Rev. D 49, 3622 (1994).
- 337
-
V.S. Berezinsky,
Comm. Nucl. Part. Phys. 21, 249 (1994).
- 338
-
J.N. Bahcall,
Phys. Lett. B 338, 276 (1994).
- 339
-
S. Degl'Innocenti, G. Fiorentini and M. Lissia,
Nucl. Phys. B (Proc. Suppl.) 43, 66 (1995).
- 340
-
S. Parke,
Phys. Rev. Lett. 74, 839 (1995).
- 341
-
K.H. Heeger and R.G.H. Robertson,
Phys. Rev. Lett. 77, 3270 (1996).
- 342
-
N. Hata and P. Langacker,
Phys. Rev. D 56, 6107 (1997).
- 343
-
H. Minakata and H. Nunokawa,
Phys. Rev. D 59, 073004 (1999).
- 344
-
G.L. Fogli, E. Lisi and D. Montanino,
Astropart. Phys. 9, 119 (1998).
- 345
-
M. Narayan, M.V.N. Murthy, G. Rajasekaran and S. Uma Sankar,
Phys. Rev. D 53, 2809 (1996).
- 346
-
S. Goswami, K. Kar and A. Raychaudhuri,
Int. J. Mod. Phys. A 12, 781 (1997).
- 347
-
G.L. Fogli, E. Lisi and D. Montanino,
Phys. Rev. D 54, 2048 (1996).
- 348
-
P. Osland and G. Vigdel,
Phys. Lett. B 438, 129 (1998).
- 349
-
P.I. Krastev and S.T. Petcov,
Phys. Rev. D 53, 1665 (1996).
- 350
-
P.I. Krastev, Q.Y. Liu and S.T. Petcov,
Phys. Rev. D 54, 7057 (1996).
- 351
-
M. Maris and S.T. Petcov, Phys. Rev. D 56, 7444 (1997); Phys. Rev. D
58, 113008 (1998).
- 352
-
L. Oberauer,
Borexino Coll., Talk presented at Neutrino '98 [434], 1998.
- 353
-
Y. Suzuki (SuperKamiokande), Proc. of the Fourth International Solar
Neutrino Conference, Heidelberg, Germany, 8-11 April 1997, edited by W.
Hampel, p. 163, Max-Planck-Institut für Kernphysik, 1997; R. Meijer Drees
(SNO), ibid., p. 210; F. von Feilitzsch (Borexino), ibid.,
p. 192; E. Bellotti (GNO), ibid., p. 173; K. Lande (Homestake Iodine),
ibid., p. 228; C. Tao (HELLAZ), ibid., p. 238; A.V. Kopylov
(Lithium), ibid., p. 263; Yu.G. Zdesenko (Xenon), ibid.,
p. 283; P. Cennini et al. (ICARUS), LNGS-94/99-I, May 1994; T.J.
Bowels (GaAs), Proc. of Neutrino 96, Helsinki, Finland, 13-19 June
1996, edited by K. Enqvist et al., p. 83, World Scientific,
Singapore, 1997.
- 354
-
R.E. Lanou,
Talk presented at Neutrino '98 [434], 1998
(hep-ex/9808033).
- 355
-
S.M. Bilenky and C. Giunti,
Phys. Lett. B 311, 179 (1993).
- 356
-
S.M. Bilenky and C. Giunti, Phys. Lett. B 320, 323 (1994); Astrop.
Phys. 2, 353 (1994); Z. Phys. C 68, 495 (1995).
- 357
-
W.K. Kwong and S.P. Rosen,
Phys. Rev. D 51, 6159 (1995).
- 358
-
W.K. Kwong and S.P. Rosen,
Phys. Rev. D 54, 2043 (1996).
- 359
-
F.L. Villante, G. Fiorentini and E. Lisi,
Phys. Rev. D 59, 013006 (1999).
- 360
-
J.N. Bahcall and P.I. Krastev,
Phys. Rev. C 55, 929 (1997).
- 361
-
J.N. Bahcall and P.I. Krastev,
Phys. Rev. C 56, 2839 (1997).
- 362
-
G.L. Fogli, E. Lisi and D. Montanino,
Phys. Lett. B 434, 333 (1998).
- 363
-
B. Faid, G.L. Fogli and E. Lisi D. Montanino,
Astropart. Phys. 10, 93 (1999).
- 364
-
SNO WWW page: http://www.sno.queensu.ca.
- 365
-
Borexino WWW pages: http://almime.mi.infn.it/;
http://pupgg.princeton.edu/~borexino/welcome.html.
- 366
-
C. Athanassopoulos et al.,
LSND Coll., Phys. Rev. Lett. 75, 2650 (1995).
- 367
-
C. Athanassopoulos et al.,
LSND Coll., Phys. Rev. C 54, 2685 (1996).
- 368
-
C. Athanassopoulos et al.,
LSND Coll., Nucl. Instrum. Methods Phys. Res. A 388, 149 (1997).
- 369
-
C. Athanassopoulos et al.,
LSND Coll., preprint nucl-ex/9706006 (1997).
- 370
-
D.H. White,
Talk presented at Neutrino '98 [434], 1998.
- 371
-
LSND WWW page: http://www.neutrino.lanl.gov/LSND/.
- 372
-
L.A. Ahrens et al.,
BNL E734 Coll., Phys. Rev. D 36, 702 (1987).
- 373
-
L. Borodovsky et al.,
BNL E776 Coll., Phys. Rev. Lett. 68, 274 (1992).
- 374
-
A. Romosan et al.,
CCFR Coll., Phys. Rev. Lett. 78, 2912 (1997).
- 375
-
B. Zeitnitz,
KARMEN Coll., Talk presented at Neutrino '98 [434], 1998.
- 376
-
Neutrino Oscillation Workshop NOW'98, Amsterdam, 7-9 September 1998; WWW page:
http://www.nikhef.nl/pub/conferences/now98/.
- 377
-
G. Drexlin et al.,
KARMEN Coll., Nucl. Instrum. Methods Phys. Res. A 289, 490 (1990).
- 378
-
G. Drexlin et al.,
KARMEN Coll., Prog. Part. Nucl. Phys. 32, 375 (1994).
- 379
-
B. Armbruster,
KARMEN Coll., Talk presented at the XXXIII
Rencontres
de Moriond: Electroweak Interactions and Unified Theories, Les Arcs 1800
(France), March 14-21, 1998
[382].
- 380
-
K. Eitel and B. Zeitnitz,
KARMEN Coll., Contribution to Neutrino '98 [434], 1998
(hep-ex/9809007).
- 381
-
B. Armbruster,
KARMEN Coll., Contribution to NOW'98 [376], 1998
(see also [382]).
- 382
-
KARMEN WWW page:
http://www-ik1.fzk.de/www/karmen/karmen_e.html.
- 383
-
G.J. Feldman and R.D. Cousins,
Phys. Rev. D 57, 3873 (1998).
- 384
-
R.D. Cousins,
Am. J. Phys. 63, 398 (1995).
- 385
-
C. Giunti,
Phys. Rev. D 59, 053001 (1998).
- 386
-
C. Giunti,
in the Summary of the Phenomenology Working Group at NOW'98 [376], 1998.
- 387
-
Booster Neutrino Experiment (BooNE), http://nu1.lampf.lanl.gov/BooNE.
- 388
-
I-216
proposal at CERN,
http://chorus01.cern.ch/~pzucchel/loi/.
- 389
-
Oak Ridge Large Neutrino Detector, http://www.phys.subr.edu/orland/.
- 390
-
NESS: Neutrinos at the European Spallation Source,
http://www.isis.rl.ac.uk/ess/neut%5Fess.htm.
- 391
-
C.Y. Cardall and G.M. Fuller,
Phys. Rev. D 53, 4421 (1996).
- 392
-
K. Kang, J.E. Kim and P. Ko,
Z. Phys. C 72, 671 (1996).
- 393
-
A. Acker and S. Pakvasa,
Phys. Lett. B 397, 209 (1997).
- 394
-
T. Teshima, T. Sakai and O. Inagaki,
preprint hep-ph/9801276 (1998).
- 395
-
R.P. Thun and S. McKee,
Phys. Lett. B 439, 123 (1998).
- 396
-
G. Barenboim and F. Scheck,
Phys. Lett. B 440, 332 (1998).
- 397
-
G. Conforto, M. Barone and C. Grimani,
Phys. Lett. B 447, 122 (1998).
- 398
-
G.L. Fogli, E. Lisi, D. Montanino and G. Scioscia,
Phys. Rev. D 56, 4365 (1997).
- 399
-
G.L. Fogli,
in the Summary of the Phenomenology Working Group at NOW'98 [376], 1998.
- 400
-
R. Foot, Mod. Phys. Lett. A 9, 169 (1994); R. Foot and R.R. Volkas,
Phys. Rev. D 52, 6595 (1995); Z. Berezhiani and R.N. Mohapatra,
ibid 52, 6607 (1995); J. Bunn, R. Foot and R.R. Volkas, Phys.
Lett. B 413, 109 (1997); B. Brahmachari and R.N. Mohapatra, Phys.
Lett. B 437, 100 (1998); J.P. Bowes and R.R. Volkas, J. Phys. G
24, 1249 (1998); M. Collie and R. Foot, Phys. Lett. B 432,
134 (1998); R.R. Volkas and Y.Y.Y. Wong, Phys. Rev. D 58, 113001
(1998); D. Suematsu, Prog. Theor. Phys. 99, 483 (1998);
hep-ph/9805305; W. Krolikowski, Acta Phys. Polon. B 29, 629 (1998);
ibid. 29, 755 (1998); hep-ph/9803323; hep-ph/9808207;
hep-ph/9808307; Y. Koide and H. Fusaoka, Phys. Rev. D 59, 053004
(1999); A. Geiser, Phys. Lett. B 444, 358 (1999).
- 401
-
A. Geiser,
in the Summary of the Phenomenology Working Group at NOW'98 [376], 1998.
- 402
-
A. De Rujula et al.,
Nucl. Phys. B 168, 54 (1980).
- 403
-
V. Barger and K. Whisnant,
Phys. Lett. B 209, 365 (1988).
- 404
-
S.M. Bilenky, M. Fabbrichesi and S.T. Petcov,
Phys. Lett. B 276, 223 (1992).
- 405
-
S.M. Bilenky, C. Giunti and C.W. Kim,
Astrop. Phys. 4, 241 (1996).
- 406
-
X. Shi and D.N. Schramm,
Phys. Lett. B 283, 305 (1992).
- 407
-
S.M. Bilenky and C. Giunti,
Phys. Lett. B 444, 379 (1998).
- 408
-
R. Svoboda,
Super-Kamiokande Coll., Talk presented at the Conference on Solar
Neutrinos: News About SNUs, Santa Barbara, California, 2-6 December, 1997
(http://www.itp.ucsb.edu/online/snu/).
- 409
-
M.C. Gonzalez-Garcia et al.,
Phys. Rev. D 58, 033004 (1998).
- 410
-
F. Vissani,
preprint hep-ph/9708483 (1997).
- 411
-
V. Barger, S. Pakvasa, T.J. Weiler and K. Whisnant, Phys. Lett. B 437,
107 (1998); A.J. Baltz, A.S. Goldhaber and M. Goldhaber, Phys. Rev. Lett.
81, 5730 (1998); Y. Nomura and T. Yanagida, Phys. Rev. D
59, 017303 (1999); G. Altarelli and F. Feruglio, Phys. Lett. B
439, 112 (1998); M. Tanimoto, Phys. Rev. D 59, 017304
(1999); E. Ma, Phys. Lett. B 442, 238 (1998); N. Haba, Phys. Rev. D
59, 035011 (1999); H. Fritzsch and Z.Z. Xing, Phys. Lett. B
440, 313 (1998); M. Jezabek and Y. Sumino, Phys. Lett. B
440, 327 (1998); H. Georgi and S.L. Glashow, hep-ph/9808293; S.
Davidson and S.F. King, Phys. Lett. B 445, 191 (1998); R.N.
Mohapatra and S. Nussinov, Phys. Lett. B 441, 299 (1998);
hep-ph/9809415; G. Altarelli and F. Feruglio, JHEP 9811, 021 (1998); C.
Jarlskog, M. Matsuda, S. Skadhauge and M. Tanimoto, hep-ph/9812282.
- 412
-
C. Giunti,
Phys. Rev. D 59, 077301 (1999).
- 413
-
M. Apollonio et al., CHOOZ Coll., Phys. Lett. B 420, 397; C.
Bemporad, Talk presented at the Ringberg Euroconference New Trends in
Neutrino Physics, Ringberg Castle, Tegernsee, Germany, May, 1998.
- 414
-
K. Zuber,
Invited talk at the International Workshop on Simulations and Analysis
Methods for Large Neutrino Detectors, Zeuthen, Germany, 6-9 Jul, 1998
(hep-ex/9810022).
- 415
-
J.T. Peltoniemi, D. Tommasini and J.W.F. Valle, Phys. Lett. B 298, 383
(1993); J.T. Peltoniemi and J.W.F. Valle, Nucl. Phys. B 406, 409
(1993); D.O. Caldwell and R.N. Mohapatra, Phys. Rev. D 48, 3259
(1993); E. Ma and P. Roy, ibid 52, R4780 (1995); S. Goswami,
ibid 55, 2931 (1997); A.Yu. Smirnov and M. Tanimoto,
ibid 55, 1665 (1997); N. Gaur et al., ibid
58, 071301 (1998); E.J. Chun, C.W. Kim and U.W. Lee, ibid
58, 093003 (1998); V. Barger et al., ibid
58, 093016 (1998); K. Benakli and A.Yu. Smirnov, Phys. Rev. Lett.
79, 4314 (1997); E.J. Chun et al., Phys. Lett. B
357, 608 (1995); V. Barger, T.J. Weiler and K. Whisnant, ibid
427, 97 (1998); S.C. Gibbons et al., ibid
430, 296 (1998); J.J. Gomez-Cadenas and M.C. Gonzalez-Garcia, Z.
Phys. C 71, 443 (1996); Q.Y. Liu and A.Yu. Smirnov, Nucl. Phys.
524, 505 (1998); S. Mohanty, D.P. Roy and U. Sarkar, Phys. Lett. B
445, 185 (1998).
- 416
-
S.M. Bilenky, C. Giunti and W. Grimus,
Eur. Phys. J. C 1, 247 (1998).
- 417
-
N. Okada and O. Yasuda,
Int. J. Mod. Phys. A 12, 3669 (1997).
- 418
-
I.E. Stockdale et al.,
CCFR Coll., Phys. Rev. Lett. 52, 1384 (1984).
- 419
-
S.M. Bilenky, C. Giunti and W. Grimus,
Proc. of Neutrino '96, Helsinki, June 1996, edited by K. Enqvist
et al., p. 174, World Scientific, Singapore, 1997
(hep-ph/9609343).
- 420
-
S.M. Bilenky, C. Giunti and W. Grimus,
Proc. of the 5
Workshop on High Energy Physics Phenomenology,
IUCAA, Pune, India, January, 12-26, 1998
(hep-ph/9805368).
- 421
-
S.M. Bilenky, C. Giunti, C.W. Kim and S.T. Petcov,
Phys. Rev. D 54, 4432 (1996).
- 422
-
G. Conforto, A. Marchionni, F. Martelli and F. Vetrano,
Phys. Lett. B 427, 314 (1998).
- 423
-
S.M. Bilenky, C. Giunti and W. Grimus,
Phys. Rev. D 57, 1920 (1998).
- 424
-
S.M. Bilenky, C. Giunti, W. Grimus and T. Schwetz,
preprint hep-ph/9804421 (1998)
[to be published in Astropart. Phys.].
- 425
-
S. Sarkar, Rept. Prog. Phys. 59, 1493 (1996); D.N. Schramm and M.S.
Turner, Rev. Mod. Phys. 70, 303 (1998); K.A. Olive, preprint
astro-ph/9707212; Nucl. Phys. (Proc. Suppl.) 70, 521 (1999); A.D.
Dolgov, preprint astro-ph/9807134.
- 426
-
R. Foot, M.J. Thomson and R.R. Volkas, Phys. Rev. D 53, 5349 (1996);
X. Shi, ibid. 54, 2753 (1996); R. Foot and R.R. Volkas,
ibid. 55, 5147 (1997); R. Foot, preprint hep-ph/9809315.
- 427
-
S.M. Bilenky, C. Giunti, C.W. Kim and M. Monteno,
Phys. Rev. D 57, 6981 (1998).
- 428
-
S.M. Bilenky, C. Giunti and W. Grimus,
Proc. of Neutrino '98, Takayama, Japan, June, 1998
(hep-ph/9809368).
- 429
-
M. Günther et al., Phys. Rev. D 55, 54 (1997); Phys. Lett.
B 407, 219 (1997).
- 430
-
G. Pantis, F. Simkovic, J.D. Vergados and A. Faessler, Phys. Rev. C
53, 695 (1996); F. Simkovic, J. Schwieger, G. Pantis and A.
Faessler, Found. Phys. 27, 1275 (1997).
- 431
-
Talks presented by A. Morales, H. Ejiri, F. Piquemal, H.V.
Klapdor-Kleingrothaus and O. Cremonesi at Neutrino '98 [434].
- 432
-
F. Wilczek,
Talk presented at Neutrino '98 [434], 1998
(hep-ph/9809509).
- 433
-
C. Ryan and S. Okubo,
Il Nuovo Cim. Suppl. 2, 234 (1964).
- 434
-
International Conference on Neutrino Physics and Astrophysics Neutrino
'98, Takayama, Japan, June 1998; WWW page:
http://www-sk.icrr.u-tokyo.ac.jp/nu98.
- 435
-
NOMAD WWW page: http://nomadinfo.cern.ch/.
This document was generated using the
LaTeX2HTML translator Version 2002 (1.62)
Copyright © 1993, 1994, 1995, 1996,
Nikos Drakos,
Computer Based Learning Unit, University of Leeds.
Copyright © 1997, 1998, 1999,
Ross Moore,
Mathematics Department, Macquarie University, Sydney.
The command line arguments were:
latex2html -html_version 4.0 -local_icons -white -split 0 -image_type gif -nonavigation 9812360.tex
The translation was initiated by Carlo Giunti on 2003-11-11
Footnotes
- ....1
- Notice that in the Goldhaber et al.
experiment the helicity of the electron neutrino was measured. The
helicity of the muon neutrino was measured in several experiments (for
the references see the review of V.L. Telegdi [13]). The
best accuracy in the measurement of the muon neutrino was achieved in
the experiment by Grénacs et al. [14].
- ... leptons.2
- One needs
to assume the existence of right-handed SU(2)
singlet fields
.
If these fields do not exist,
it is not possible to construct
a Dirac mass term for the neutrinos,
i.e., neutrinos are massless particles.
This model is sometimes called Minimal Standard Model
[74,75,76].
- ...
by3
- Note that in this review we use four-component spinors. For
the equivalent alternative of two-component spinors see, e.g.,
Refs. [30,80].
- ...
neutrinos.4
- In the Dirac case it is possible to conceive
also the existence of left-handed sterile singlet fields mixed with the
active flavour neutrino fields, in contrast to the right-handed
singlets which cannot be mixed with the active neutrinos because
of lepton number conservation.
- ...
fields.5
- This can also be seen
by noticing that the Jarlskog rephasing-invariant parameter
[85,86,87]
vanishes if one of the elements of the mixing matrix is zero.
- ... fermions.6
- Note that the see-saw mechanism can also be
realized for Dirac neutrinos [110,111].
- ... is7
- This
operator can be written also as
- ... one8
- We use the symbol
for a generic mass-squared difference.
Hence
can be any of the
's.
- ... respectively,9
- For protons and neutrons
this is ensured by CVC (conserved vector current)
for the
part and conservation of the electromagnetic current.
- ...).10
- We use
the dot above a symbol as abbreviation for
.
- ....11
- The scalar product in Eq.(4.22) is
imaginary for arbitrary vectors
normalized to one
, therefore, for real vectors it must be zero.
- ....12
- Without loss of generality we consider
.
- ... sheet.13
- It
is convenient to imagine cuts from
to
and
to
in the complex plane. Crossing these
cuts, one moves from one sheet to the other.
- ....14
- In general this change is accompanied by a
factor of geometrical origin. It has no effect in our case for the
jumping probability [198].
- ...
consequently,15
- In general, we denote particle fluxes by putting
the particle symbol within parentheses.
- ... results:16
- For recent results on
atmospheric neutrinos see Note added.
- ... effect17
- Note that for the case of a periodical
matter density the effect of enhancement was considered
in Refs. [281].
- ... core18
- Neutrinos with
- ... release.19
- Here
is the proton mass,
is the electron mass
and
is the mass of the
-nucleus,
where
is the neutron mass
and
is the binding energy of the
-nucleus.
- ...
constraint20
- The corrections to this relation
due to the generation of gravitational energy
[174]
and due to the fact that the abundance of
nuclei is out of equilibrium in the outer regions of the solar core
[321]
are estimated to be less than 1%
and thus negligible at the present level of accuracy.
- ...
be21
- One Solar Neutrino Unit (SNU) is defined as
.
- ... experiments).22
- For
recent Super-Kamiokande results on solar neutrinos see Note added.
- ... flux.23
- Non-standard solar models
with low central temperature yield a suppression of the
neutrino flux that is stronger than that
of the
neutrino flux
(see [312] and references therein).
- ... region.24
- The KARMEN 2 result
has also been analysed using the Unified Approach [383],
which is a frequentist method with proper coverage
(see Ref. [384]).
The resulting exclusion curve seems to forbid almost all of the
LSND-allowed region [380,382],
but this result depends on the non-observation of the expected
background events.
The analysis of the KARMEN 2 null result with an alternative
frequentist method with proper coverage [385]
leads to an exclusion curve compatible
with the Bayesian exclusion curve shown in Fig. 5.9
[385,386].
- ... neutrinos.25
- From time to time
special schemes
with three neutrinos have been proposed
in order to accommodate all neutrino oscillation data
[391,392,393,394,395,396,397]
However,
none of these schemes can fit all data in a satisfactory way
[398,399].
Schemes with more than four neutrinos have been considered in the papers in
Refs. [130,400]
(see the review Ref. [401]).
- ... hierarchy26
- Actually,
the phenomenology of neutrino oscillations discussed in this section
depends only on the hierarchy
of mass-squared differences.
Hence,
the discussion of neutrino oscillations presented in this section
applies, for example, also to the case of three almost degenerate
neutrinos with masses in the eV range
that give a contribution to the dark matter in the universe.
Furthermore,
with a cyclic permutation of the indices 1, 2, 3,
the discussion of neutrino oscillations presented in this section
can be applied to all cases in which
,
as, for example,
the case of an inverted mass hierarchy
.
- ...
neutrinos27
- As explained after Eq.(6.5),
in the case of the MSW solutions of the solar neutrino problem
this inequality and the following equations are
valid only for multi-GeV neutrinos.
- ...4spectra.28
- The presence of
a small allowed region in Fig. 6.3 for
is due to the fact that
for
there are no constraints on
from the results of
disappearance experiments.
- ...
respectively.29
- If
CP is conserved, we have
or
if the CP parities of
and
are equal
and
if they are opposite.
Indeed, if CP is conserved from Eq.(2.50)
we get
where
and the CP parity of
is
(see Eq.(2.53))
[96,97].
Writing
,
we obtain
or
if
and
if
.
Hence,
we have
or
if
and
if
.
Carlo Giunti
2003-11-11