%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This is a sample input file for your contribution to a
% multi-author book ("proceedings-like) to be published
% by Springer Verlag.
%
% Please use it as a template for your own input, and please
% follow the instructions for the formal editing of your
% manuscript as described in the file "1readme".
%
% Please send the Tex and figure files of your manuscript
% together with any additional style files as well as the
% PS file to the editor of your book.
%
% He or she will collect all contributions for the planned
% book, possibly compile them all in one go and pass the
% complete set of manuscripts on to Springer.
%
% Vers. 07/99
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%RECOMMENDED%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
\documentclass[runningheads,colorBG]{svmult}
%
\usepackage{makeidx}   % allows index generation
\usepackage{graphicx}  % standard LaTeX graphics tool
%                       % for including eps-figure files
\usepackage{subeqnar}  % subnumbers individual equations
%                       % within an array
\usepackage{multicol}  % used for the two-column index
\usepackage{cropmark} % cropmarks for pages without
%                       % pagenumbers - only needed when manuscript
%                       % is printed from paper and not from data
\usepackage{physprbb}  % centered layout of diverse elements, etc.
%\makeindex             % used for the subject index
%                       % please use the style sprmidx.sty with
%                       % your makeindex program
%
%%%upright Greek letters (example below: upright "mu")
\newcommand{\greeksym}[1]{{\usefont{U}{psy}{m}{n}#1}}
\newcommand{\umu}{\mbox{\greeksym{m}}}
\newcommand{\udelta}{\mbox{\greeksym{d}}}
\newcommand{\uDelta}{\mbox{\greeksym{D}}}
\newcommand{\uPi}{\mbox{\greeksym{P}}}
\newcommand{\newc}{\newcommand} 
\newc{\beq}{\begin{equation}} 
\newc{\eeq}{\end{equation}} 
\newc{\barr}{\begin{eqnarray}} 
\newc{\earr}{\end{eqnarray}} 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%OPTIONAL%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%\usepackage{amstex}   % useful for coding complex math
%\mathindent\parindent % needed in case "Amstex" is used
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%AUTHOR_STYLES_AND_DEFINITIONS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%
\usepackage{axodraw} 
%
\def\e6{$E(6)$}
\def\10{$SO(10)$}
\def\21{$SU(2) \otimes U(1) $}
\def\lr{$SU(2)_L \otimes SU(2)_R \otimes U(1)$}
\def\422{$SU(4) \otimes SU(2) \otimes SU(2)$}
\def\321{$SU(3) \otimes SU(2) \otimes U(1)$}
%\def\r#1{\Red1}
%%\def\np#1#2#3{           {\it Nucl. Phys. }{\bf #1} (19#2) #3}
\def\ne{\hbox{$\nu_e$ }}
\def\nm{\hbox{$\nu_\mu$ }}
\def\bne{\hbox{$\bar\nu_e$ }}
\def\bnm{\hbox{$\bar\nu_\mu$ }}
\def\nt{\hbox{$\nu_\tau$ }}
\def\ns{\hbox{$\nu_{s}$ }}
\def\nr{\hbox{$\nu_R$}}
\def\O{\hbox{$\cal O$ }}
\def\L{\hbox{$\cal L$ }}
\def\ed{\end{document}}
\def \nbb {$\beta\beta_{0\nu}$ }
\def\Eb{E_{\rm b}}
\def\Ee{\langle E_{\bar\nu_{\rm e}}\rangle}
\def\Eh{\langle E_{\bar\nu_{\rm h}}\rangle}
%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newcommand{\pnu}[1] {\overset{\smash{\scriptscriptstyle (-)}}{\nu}_{\hskip-3pt #1}}
\newcommand{\eps}{\epsilon_{13}}
\newcommand{\si}{S_{13}}
%%%%%%%%%%%%%%%%%%%%%%%%%%
\newcommand{\snoccfluxshort}{1.76^{+0.06}_{-0.05}\mbox{(stat.)}^{+0.09}_{-0.09}~\mbox{(syst.)}} 
\newcommand{\snoesfluxshort}{2.39^{+0.24}_{-0.23}\mbox{(stat.)}^{+0.12}_{-0.12}~\mbox{(syst.)}} 
\newcommand{\snoncfluxshort}{5.09^{+0.44}_{-0.43}\mbox{(stat.)}^{+0.46}_{-0.43}~\mbox{(syst.)}} 
\newcommand{\snoncfluxunc}{6.42^{+1.57}_{-1.57}\mbox{(stat.)}^{+0.55}_{-0.58}~\mbox{(syst.)}} 
\newcommand{\snomutauflux}{3.41^{+0.45}_{-0.45}\mbox{(stat.)}^{+0.48}_{-0.45}~\mbox{(syst.)}} 
\newcommand{\snoeflux}{1.76^{+0.05}_{-0.05}\mbox{(stat.)}^{+0.09}_{-0.09}~\mbox{(syst.)}} 
\newcommand{\snomutaufluxsk}{3.45^{+0.65}_{-0.62}} 
\newcommand{\snomutaufluxcomb}{3.41^{+0.66}_{-0.64}} 
\newcommand{\nccfit}{1967.7^{+61.9}_{-60.9}} 
\newcommand{\nesfit}{263.6^{+26.4}_{-25.6}} 
\newcommand{\nncfit}{576.5^{+49.5}_{-48.9}}  
\newcommand{\nsigmassno}{5.3} 
\newcommand{\nsigmassk}{5.5} 
\newcommand{\ssmflux}{5.05^{+1.01}_{-0.81}}
\newcommand{\phisk}{2.32\pm0.03\mbox{(stat.)}^{+0.08}_{-0.07}~\mbox{(syst.)}}
%Please reduce your own definitions and macros to an absolute
%minimum since otherwise the editor will find it rather
%strenuous to compile all individual contributions to a
%single book file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{document}
%
%\title*{Neutrinos: theory meets facts}
%\title*{A summary of neutrino theory and experiment}
%\title*{\Red {Neutrinos: A summary of  theory and experiment}}
%\title*{\Red  {Neutrinos, where do we stand? }}
\title*{\Red  {Neutrinos: Summarizing the State-of-the-Art }}
\toctitle{ Neutrinos}
% allows explicit linebreak for the table of content
%
%
\titlerunning{ Neutrinos}
% allows abbreviation of title, if the full title is too long
% to fit in the running head
%
\author{J. W. F. Valle}
%
\authorrunning{J. W. F. Valle}
% if there are more than two authors,
% please abbreviate author list for running head
%
%
\institute{Instituto de F\'{\i}sica Corpuscular -- C.S.I.C., 
Universitat de Val{\`e}ncia \\
Edificio Institutos, Aptdo.\ 22085, E--46071 Val{\`e}ncia, Spain\\
E-mail:valle@ific.uv.es}
\maketitle              % typesets the title of the contribution

\begin{abstract}
    
  I review oscillation solutions to the neutrino anomalies and discuss
  how to account for the required pattern of neutrino masses and
  mixings from first principles. Unification and low-energy bottom-up
  approaches are discussed, the latter open up the possibility of
  testing neutrino mixing at high energy colliders, such as the LHC.
  Large \ne mixing is consistent with Supernova (SN) astrophysics and
  may serve to probe galactic SN parameters at Cherenkov detectors.
%%
  I discuss the robustness of the atmospheric neutrino oscillation
  hypothesis against the presence of Flavor Changing (FC) Non-Standard
  neutrino Interactions (NSI), generally expected in models of
  neutrino mass.  Atmospheric data strongly constrain FC-NSI in the
  \nm-\nt channel, while solar data can be explained by FC-NSI in the
  \ne-\nt channel, or, alternatively, by spin flavor precession.
%
  I illustrate how a neutrino factory offers a unique way to probe for
  FC-NSI and argue that a near-site detector is necessary in order to
  probe for leptonic mixing and CP violation.


\end{abstract}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\date{\today}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{ Solar Neutrinos \cite{sol,ssm,Kachelriess:2001sg,Gonzalez-Garcia:1999aj,Ahmad:2002jz,Ahmad:2002ka}}
\label{sec:solar-neutrinos}

%\section{ Neutrino Anomalies }

For decades, the number of solar neutrinos detected in underground
experiments \cite{sol} has been less than expected from theories of
energy generation in the sun \cite{ssm}, suggesting that either the
understanding of the Sun was incomplete, or that the neutrinos were
changing from one type to another in transit from the core of the Sun
to the detector.
%
At the moment we have a conclusive proof for the latter possibility,
though we still can not pin down the exact neutrino conversion
mechanism implied.

Solar neutrinos have now been detected through the geochemical method
(the \ne + ${}^{37}$Cl $\to {}^{37}$Ar + $e^-$ reaction at the
Homestake experiment and the \ne + ${}^{71}$Ga $\to ^{71}$Ge + $e^-$
reaction at the Gallex, Sage and GNO experiments) and also directly
using Cherenkov techniques at Kamiokande and SuperKamiokande (SK) and
the Sudbury Neutrino Observatory (SNO).  At SK the reaction is $\nu_e
e$ scattering on water, while SNO detects neutrinos using 1000 tonnes
of heavy water (D2O) and is sensitive to the charged current (CC) and
neutral current (NC), as well as to the elastic scattering (ES)
reaction.

At present all experiments observe a deficit of 30 to 60 \% with an
energy dependence mainly due to the lower Chlorine rate.
%
The high statistics of SK after 1496 days of data--taking also
provides very useful information on the recoil electron energy
spectrum and this is presently well described by the flat hypothesis,
placing important restrictions on neutrino parameters.
%
Moreover, SK measures the zenith angle distribution (day/night effect)
which is sensitive to the effect of the Earth matter in the neutrino
propagation, observing no significant effect.
%
The absence of a clear hint of spectral distortion, day-night or
seasonal variation implies that, to a large extent the solar neutrino
problem rests heavily on the rate discrepancy and, from this point of
view, future experiments are most welcome.

The global status of the interpretation of the solar neutrino data in
terms of active neutrino oscillations is illustrated in Fig.
\ref{fig:solchi2}.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=5cm,height=4cm]{chi2solteta.eps}  
\includegraphics[width=5cm,height=4cm]{chi2soldm.eps}
  \caption{The solar $\chi^2$ profiles versus angle and $\Delta m^2$, 
    from \cite{Kachelriess:2001sg}, before
    \cite{Ahmad:2002jz,Ahmad:2002ka}.}
  \label{fig:solchi2}
\end{figure}
%
This figure refers to the 1258 day SK data sample.  As can be seen the
best of the oscillation solutions is the Large Mixing Angle (LMA)
solution, confirming first pre-SNO hints obtained in ref.
\cite{Gonzalez-Garcia:1999aj}. This was derived solely on the basis of
the observed flatness of the recoil electron energy spectrum of SK,
highlighting its importance in shaping up what is now much more
strongly supported.

Solar neutrinos from the decay of $^8$B have been detected at SNO via
the CC reaction on deuterium and by the ES reaction.  Such CC reaction
is sensitive only to \ne's, while the ES reaction also has a small
sensitivity to \nm's and \nt's. SNO has recently observed also NC
neutrino interactions \cite{Ahmad:2002jz}.

Using the standard ${}^{8}$B shape they measure the $\nu_{e}$
component of the ${}^{8}$B solar flux $\phi_{e}$ above 5 MeV as well
as the non-$\nu_{e}$ component $\phi_{\mu\tau}$. As shown in Fig.
\ref{fig:fig3sno} the latter is significantly above zero, providing
strong evidence for solar \ne conversion into another active neutrino
flavor.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=7cm,height=4.5cm]{fig3sno.eps}  
  \caption{SNO neutral current measurement \cite{Ahmad:2002jz}.}
  \label{fig:fig3sno}
\end{figure}
%
This shows the flux of ${}^{8}$B solar neutrinos which are $\mu$ or
$\tau$ flavor versus the \ne flux derived from the three neutrino
reactions in SNO.  The diagonal bands show the total ${}^{8}$B flux as
predicted by the SSM~\cite{ssm} (dashed lines) and that measured with
the NC reaction in SNO (solid band).  The intercepts of these bands
with the axes represent the $\pm 1\sigma$ errors.  The bands intersect
at the fit values for $\phi_{e}$ and $\phi_{\mu\tau}$, indicating that
the combined flux results are consistent with neutrino flavor
conversion with no distortion in the ${}^{8}$B neutrino energy
spectrum.
%
They also find the total flux measured with the NC reaction to be
consistent with solar models.

SNO has also measured day and night solar neutrino energy spectra and
rates \cite{Ahmad:2002ka}.  They give an analysis in terms of
matter-enhanced oscillations of two active flavors showing how it
strongly favors the LMA solution. The update of global fits including
these latest data worsen the status of non-LMA and especially that of
non-active solutions, see ref.  \cite{ssm,solat02}. We eagerly await a
possible confirmation of the LMA hypothesis at KamLAND \cite{KamLAND}.
As we will see later there are also good active non-oscillation
neutrino conversion mechanisms based on flavor-changing neutrino NSI
or spin flavor precession that can accout well for the data.

\section{Atmospheric Neutrinos \cite{atm,atmfluxes,Gonzalez-Garcia:2000sq,Gonzalez-Garcia:2002mu,Fornengo:2000sr,Fornengo:2001pm}}
\label{sec:atmosph-neutr}

Atmospheric neutrinos are produced as the decay products in hadronic
showers from particles produced in collisions of cosmic rays with air
in the upper atmosphere. These have been observed in several
experiments \cite{atm}.
%
Although individual $\nu_\mu$ or $\nu_e$ fluxes are only known to
within $30\%$ accuracy, their ratio is predicted to within $5\%$ over
energies varying from 0.1~GeV to 100~GeV~\cite{atmfluxes}.
%
There has been a long-standing discrepancy between the predicted and
measured $\mu/e$ ratio of the muon ($\nu_\mu + \bar{\nu}_\mu$) over
the electron atmospheric neutrino flux ($\nu_e+\bar{\nu}_e$) both in
water Cherenkov experiments (Kamiokande, Super-Kamiokande and IMB) as
well as in the iron calorimeter Soudan2 experiment. Moreover a strong
\texttt{zenith-angle dependence} has been observed both in the sub-GeV
and multi-GeV energy range, but only for $\mu$--like events, the
zenith-angle distributions for the $e$--like being consistent with
expectation. Such clear deficit of neutrinos coming from below is very
suggestive of \nm oscillations.
%
Zenith-angle distributions have also been recorded for upward-going
muon events in SK and MACRO which are also consistent with the \nm
oscillation hypothesis.

Of all laboratory searches for neutrino oscillation, reactors provide
the highest sensitivity in $\Delta{m}^2_{32}$ and thus are most
relevant to confront with the data from underground experiments.  The
restrictions on $\Delta{m}^2_{32}$ and $\sin^2(2\theta_{13})$ that
follow from the non observation of oscillations at the Chooz reactor
experiment are shown in Fig.~\ref{fig:chooz}.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=5.5cm,height=4cm]{chooz.ps}
  \caption{Region in  $\Delta{m}^2_{32}$ and $\sin^2(2\theta_{13})$  
  excluded by the Chooz reactor, from \cite{Gonzalez-Garcia:2000sq}.}
  \label{fig:chooz}
\end{figure}
%
The curves represent the 90, 95 and 99\% CL excluded region defined
with 2 d.~o.~f. for comparison with the Chooz published results.
%
For large $\Delta{m}^2_{32}$ this gives a stringent limit on
$\sin^2(2\theta_{13})$, but not for low $\Delta{m}^2_{32}$ values.
After combining with the full body of atmospheric neutrino data, one
finds a stringent bound $\sin^2(\theta_{13}) < 0.045$ at $99\%$ CL
\cite{Gonzalez-Garcia:2002mu}.  Neglecting $\theta_{13}$ we display in
Fig.~\ref{fig:atm} the allowed regions of atmospheric oscillation
parameters $\Delta{m}^2_{32}$--$\sin^2(2\theta_{23})$ obtained in a
global analysis.  For previous analyzes, see \cite{Fornengo:2000sr}
and references therein.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=5.5cm,height=4cm]{fig.only_osc.eps}  
  \caption{Atmospheric oscillation parameters from global analysis in 
    \cite{Fornengo:2001pm}.}
    \label{fig:atm}
\end{figure}

\section{Supernova Neutrinos 
  \cite{Kachelriess:2000fe,Minakata:2001cd}.}
\label{sec:supernova-neutrinos}

The large-angle neutrino oscillation between the electron neutrino \ne
favored by current solar neutrino data must be properly taken into
account in the interpretation of the SN neutrino data.
%
Indeed, it will induce a significant deformation of the energy
spectrum of neutrinos coming from a supernova (SN) explosion and on
this basis it has been argued \cite{snperm} to be at odds with the
historical observation of neutrinos from SN1987 \cite{Hirata:1987hu}.

The impact of $\bar\nu_e \leftrightarrow \bar\nu_{\mu,\tau}$ neutrino
oscillations on the observed $\bar\nu_e$ signal of supernova SN 1987A
can be studied by performing a maximum-likelihood analysis. The fit
parameters are the released binding energy $\Eb$ and the average
neutrino energy $\Ee$.  It was found that $\bar\nu_e \leftrightarrow
\bar\nu_{\mu,\tau}$ oscillations with large mixing angles have lower
best-fit values for $\Ee$ than small-mixing angle (SMA) oscillations.
%
Moreover, the inferred value of $\Ee$ is already lower than found in
simulations in the SMA or no-oscillation case. This apparent conflict
has been interpreted as evidence against the large mixing oscillation
solutions to the solar neutrino problem.
%
In order to quantify the degree to which the SN data favor
no-oscillations over the large mixing solutions their likelihood
ratios were used as well as a Kolmogorov-Smirnov test. The result was
that within the range of SN parameters predicted by simulations there
are regions where the LMA solution is either only marginally
disfavored or even favored over the SMA case
\cite{Kachelriess:2000fe}.
%
This implies that the LMA solution is not in conflict with the current
understanding of SN physics. In contrast, vacuum oscillation and LOW
solutions may be excluded at the $4\sigma$ level for most of the SN
parameter ranges found in simulations. The reason for the difference
amongst the large mixing solutions has to do with the role of \ne
regeneration in the Earth matter.

One can go a step further, namely to perform a combined likelihood
analysis of the observed SN1987A neutrino signal and of the solar
neutrino data \cite{Kachelriess:2001sg}.
%
In the calculation of neutrino survival and conversion probabilities a
power-law $1/r^3$ profile typical of supernova envelopes was used and
the level-crossing probability was determined analytically in a form
valid for all mixing angles \cite{Kachelriess:2001bs}.
%
It has been found that, although the spectrum swap argument is
sufficient to rule out vacuum-type solutions for most reasonable
choices of astrophysics parameters, the LOW solution may still be
acceptable.
%
On the other hand it was found that the LMA solution can easily
survive as the best overall solution, although its size is typically
reduced when compared to solar-only fits \cite{Kachelriess:2001sg}.
%
How about a future supernova?

The significant distortion of the energy spectrum of neutrinos coming
from a supernova (SN) explosion expected by active LMA \ne
oscillations can be used to extract some of the properties of SN
driven by gravitational collapse.
%
In ref. \cite{Minakata:2001cd} a simple but powerful method was
proposed to extract the original temperatures of both anti-electron
and muon/tau neutrinos at the neutrinosphere from the observation of
anti-electron neutrinos at a large water Cherenkov detector such as
Super-Kamiokande. Fig. \ref{fig:minak} illustrates the potential of
the method in determining intrinsic SN parameters.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=5cm,height=5cm]{tau-En.eps}
%
\includegraphics[width=5cm,height=5cm]{Eb-En.eps}  
%
  \caption{LMA \ne oscillations allow a determination of 
SN parameters \cite{Minakata:2001cd}.}
  \label{fig:minak}
\end{figure}
%
As can be seen one obtains a rather good determination of the
important astrophysical parameter $\tau \equiv {\langle
  E_{\bar{\nu}_\mu}\rangle}/{\langle E_{\bar{\nu}_e}\rangle }$
characterizing the non-electron flavor SN neutrino spectra. The result
is rather robust in the sense that it is not too sensitive to
astrophysical assumptions such as energy equipartition.
%
Prospects are substantially better for the projected megaton-class SK
successor, HyperKAMIOKANDE.

\section{Basic Neutrino Theory \cite{Weinberg:uk,seesaw,Mohapatra:1979ia,Schechter:1980gr,Schechter:1980gk,Chikashige:1980ui,Schechter:1981cv,Schechter:1981bd,Pal:1981rm,Schechter:1981hw}}
\label{sec:basic-neutr-theory}

The basic gauge theoretical aspects of neutrinos have been developed
in the early eighties
\cite{Weinberg:uk,seesaw,Mohapatra:1979ia,Schechter:1980gr,Schechter:1980gk,Chikashige:1980ui,Schechter:1981cv,Schechter:1981bd,Pal:1981rm,Schechter:1981hw}.
The other basic step was the formulation of the MSW effect itself
\cite{Wolfenstein:1977ue}. With both ingredients one can describe all
the exciting experimental developments of the past few years
\cite{Bilenky:1996rw}.
%
In order to make sense of present neutrino data through neutrino
oscillations, neutrinos should have mass. In analogy with the
Kobayashi-Maskawa (KM) quark mixing matrix \cite{Kobayashi:fv}, let
$K$ denote the lepton mixing matrix.  The simplest form for $K$ in a
gauge theory of the weak interaction is characterized by
\begin{itemize}
\item 
the atmospheric angle $\theta_{23} \equiv \theta_A$
\item 
 the solar
  angle $\theta_{12} \equiv \theta_S$ 
\item 
the reactor angle
  $\theta_{13} \equiv \theta_R$
\item
1 KM-like CP phase $\phi_{13}$
\item
 2 extra (Majorana) phases $\phi_1, \phi_2$ \cite{Schechter:1980gr}
\end{itemize}
The ``Majorana'' phases are physical \cite{Schechter:1980gk}. Although
they do not show up in ordinary ($\Delta L=0$) oscillations, they
affect $\Delta L=2$ processes. Unfortunately, however, the effect is
tiny due to helicity suppression, reflecting the V-A nature of the
weak interaction \cite{Schechter:1980gk}.
%
Also the effect of the ``Dirac'' CP violation due to $\phi_{13}$ is
small since it disappears as $\theta_{13} \to 0$
\cite{Schechter:1979bn}.
%
A new generation of long-baseline neutrino oscillation experiments
using a neutrino beam from the decay of muons in a storage ring is
being considered~\cite{nuFacPapers}.  Studying $\theta_{13}$ and the
resulting leptonic CP violating effects provides a challenge for such
so-called \texttt{neutrino factories}.

The structure of charged and neutral current weak interactions is more
complex in models of neutrino mass containing \21 singlet
leptons~\cite{Schechter:1980gr}: the CC is characterized by a
rectangular (effectively non-unitary) matrix, while the NC is
described by a projective matrix with non-diagonal entries coupling
different mass eigenstate neutrinos with each other. Similarly, if
there are light sterile neutrinos, the nontrivial structure of $K$
must be taken into account.

Current solar and atmospheric data can be understood well with the
usual three massive neutrino flavors. The information on the neutrino
oscillation parameters may be summarized by saying that $\theta_{23}$
is nearly maximal, $\theta_{12}$ is large, but non-maximal, and the
reactor angle $\theta_{13}$ is rather small.  In contrast there is at
the moment no information on the CP phases.
%
It is clear from the above that the atmospheric mass splitting is much
larger than the solar one, $\Delta {m_{atm}^2} \gg \Delta
{m_\odot^2}$. Such hierarchical pattern can be realized in two ways,
illustrated in Fig.~\ref{fig:3nu-patterns}.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=4cm,height=2.6cm]{solatmnorm.ps}
\includegraphics[width=4cm,height=2.6cm]{solatminv.ps}   
  \caption{Normal and inverse-hierarchical neutrino masses. The 
wide and narrow  bands represent the atmospheric  and solar 
mass splittings, respectively.}
  \label{fig:3nu-patterns}
\end{figure}
%
Last, but not least, the oscillation data are sensitive only to
neutrino mass splittings, not to the absolute scale of neutrino mass.
Nor do they have any bearing on the fundamental issue of whether
neutrinos are Dirac or Majorana particles. This brings us to the
significance of the \nbb decay \cite{Morales:1998hu} in deciding the
nature of neutrinos.  The connection between the two is given by the
\texttt{black-box theorem} which states that, in a ``natural'' gauge
theory, whatever the mechanism for inducing \nbb is, it is bound to
also yield a Majorana neutrino mass at some level, and vice-versa, as
illustrated by Fig.  \ref{fig:bbox}.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=5cm,height=4.cm]{blackbox.eps}  
  \caption{The black-box \nbb argument, from \cite{Schechter:1981bd}.}
 \label{fig:bbox} 
\end{figure}
%
Quantifying the implications of the black-box argument is a
model-dependent enterprise, but can be done once a particular model
is specified. It serves at least to highlight the importance of
confirming or refuting the present hint
\cite{Klapdor-Kleingrothaus:2001ke,Aalseth:2002dt} in a more sensitive
experiment such as GENIUS \cite{Klapdor-Kleingrothaus:1999hk}.

Can one predict neutrino properties from first principles?  Such
important enterprise is unfortunately difficult, for at least three
basic reasons.
%
First the underlying \texttt{scale} is unknown, models using anything
between the Planck or String scale, through the \e6 or \10 GUT scale,
or an intermediate scale associated with Peccei-Quinn or Left-Right
symmetry, all the way down to the weak \321 scale itself.
%
Second, the \texttt{mechanism} is unknown, and there are both tree
level, like the so-called seesaw schemes
\cite{seesaw,Schechter:1980gr,Chikashige:1980ui,Schechter:1981cv} as
well as radiative \cite{Zee:1980ai} or hybrid alternatives.
%
Similarly the status of B-L (baryon minus lepton number) symmetry is
unknown, it may either be a broken gauge symmetry, such as in the
left-right or \10 schemes, or simply a global ungauged symmetry, in
which case the corresponding Goldstone boson -- called majoron -- is
physical \cite{Chikashige:1980ui,Schechter:1981cv}.
%
Last, but not least, there is \texttt{no theory of flavour}, although
many Yukawa textures have been suggested which would arise from
suitable family symmetries.

Renormalizable neutrino masses are forbidden by the Standard Model
(SM), given its restricted field content and symmetries.  The simplest
and most basic way to add neutrino masses is to regard the SM only as
the effective low-energy limit of some, unknown, higher energy theory,
which may well contain an L-violating dimension-five term
\cite{Weinberg:uk}
%
\begin{equation}
{\cal L}_5= \frac{\lambda_{ij}}{M_5}(\ell_i \phi)(\ell_j \phi),
\label{5-dim}
\end{equation} 
%
involving the lepton doublet fields $\ell_i$ ($i,j=e,\mu,\tau$), and
the Higgs doublet field $\phi$. 
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=4.5cm,height=3cm]{d-5.ps}  
  \caption{Neutrino masses from dimension-five operator \cite{Weinberg:uk}.}
  \label{fig:d-5}
\end{figure}
%
Here $M_5$ is an arbitrary mass scale, related to the unknown high
energy theory of which the SM is an effective theory and
$\lambda_{ij}$ are dimensionless couplings.
%
After electroweak symmetry breaking, eq. (\ref{5-dim}) leads to
neutrino Majorana masses $m_{\nu}^{ij}$ given by
%
\begin{equation}
m^{\nu}_{ij}=\frac{\lambda_{ij}v^2}{M_5}.
\label{mass}
\end{equation}
%
where $v=174$~GeV is the Higgs vacuum expectation value (vev). Note
that eq.~(\ref{mass}) explicitly violates lepton number symmetry.  

One way to realize this picture is in terms of the so-called seesaw
scenario for neutrino
mass~\cite{seesaw,Mohapatra:1979ia,Schechter:1980gr,Chikashige:1980ui,Schechter:1981cv},
where this operator is induced either through the exchange of heavy
\321 singlet neutrinos or extra scalar bosons related, say, with the
breaking of extended left-right-type gauge symmetries.  The first
challenge for unification seesaw schemes is that the smallness of the
quark KM mixing angles typically suggests that leptonic mixing angles
should also be small.  However this is certainly possible to avoid
owing to the lack of predictivity. Indeed, despite its attractiveness
and the existence of many papers claiming to have ``predictive''
seesaw models of neutrino mass and mixing \cite{seesawnupapers}, such
models at most \texttt{post-dict} the observed angles and splittings
and it is hardly possible to obtain ``smoking-gun'' tests of their
validity.
  
\section{Neutrino Models}
  
Notwithstanding the limitations of theory in predicting neutrino
properties from first principles, I will mention some recent examples
of predictive theoretical schemes for accounting for the neutrino
anomalies.  Given the vast literature \cite{seesawnupapers,giunti}, I
must confine here to a personal selection of recent neutrino models.

\subsection{Neutrino Mass Unification~\cite{Chankowski:2000fp}}
\label{sec:neutr-mass-unif}

The idea of neutrino mass unification is inspired by the paradigm of
unification of fundamental interactions. This ``minimalistic'' ansatz
postulates that the neutrino mass and mixings observed at low energies
take a very simple form at some high energy \texttt{neutrino
  unification} scale $M_X$.  The model adds to the basic Lagrangian
the dimension--five operator in eq.~(\ref{5-dim}) assumed to be
characterized at the scale $M_5 = M_X$ by a single real
dimension--less parameter $\lambda_0$.
%
In this picture neutrino masses arise from a common seed at $M_X$,
while their splittings are induced by renormalization effects, as
illustrated in Fig. \ref{fig:pok}.
\begin{figure}[htbp]
  \centering
\includegraphics[width=6cm,height=3cm]{neunif2.ps}  
  \caption{Artists' view of  neutrino masses arising from a 
    common seed (1 eV) at a large neutrino unification scale
    $M_X \sim 10^{14}$ GeV~\cite{Chankowski:2000fp}.}
  \label{fig:pok}
\end{figure}
%
This model is inconsistent with small solar mixing, thus it
``predicts'' a large, though non-maximal solar angle, according to the
formula $$\Delta m^2_{S}~\approx~2 \Delta m^2_{A}
\left[\cos~2~\theta_{S} + 2~\sin^2 \theta_{R} \right],$$
%
where $\Delta m^2_{S}$ and $\Delta m^2_{A}$ are the solar and
atmospheric splittings while $\theta_{S}$ and $\theta_{R}$ are the
solar and reactor angles, respectively.
%
Neutrinos can lie in the eV range and thus be seen in tritium beta
decays \cite{Osipowicz:2001sq} as well as contribute to the hot dark
matter \cite{Elgaroy:2002bi}.  In contrast to other mechanisms leading
to quasi-degenerate neutrinos \cite{Caldwell:kn,DEG} which may be
unnatural \cite{DEGNAT}, here this is not the case because opposite
neutrino CP signs
\cite{Schechter:1981hw,Wolfenstein:1981rk,Valle:1982yw} suppress the
rate for neutrinoless double beta decay.
%
Taken at face value, an analysis including the most recent NC SNO data
\cite{ssm} suggests that the preferred solution in this case will be
in terms of vacuum oscillations, instead of LMA, which as we saw is
disfavored by SN1987A.  In order to evaluate more carefully our
prediction a three-neutrino global fit of all neutrino data, similar
to that in \cite{Gonzalez-Garcia:2000sq}, would be required.

  
\subsection{Minimalistic Seesaw Plus Gravity~\cite{deGouvea:2001jp}}
\label{sec:minimalism}

The idea here rests on the observation that gravity alone can account
for the solar neutrino data by Planck-mass L-violating effects, while
the atmospheric neutrino anomaly can be due to the existence of a
single right-handed neutrino, instead of
three~\cite{Schechter:1979bn}. One assumes it to be at an intermediate
mass scale between $10^9$ GeV and $10^{14}$ GeV. To first
approximation only one of the three flavor neutrinos acquires mass
\textit{a la seesaw}, identified as the atmospheric scale. The other
two neutrinos remain massless, as illustrated in Fig. \ref{fig:gouv}
(left panel). Their degeneracy is lifted only by tiny gravitational
effects (right panel) due to the dimension--five operator.
\begin{figure}[htbp]
  \centering
\includegraphics[width=4cm,height=3cm]{solatm0.ps}  
%
\includegraphics[width=4cm,height=3cm]{solatmin.ps}
  \caption{Minimalistic neutrino masses. The degeneracy of the 
    massless states on the left is lifted by gravitational effects,
    illustrated in the right panel.}
  \label{fig:gouv}
\end{figure}
%
Such simple scheme solves the current solar and atmospheric neutrino
puzzles with probably the smallest amount of beyond-the-Standard-Model
ingredients. Note that although the neutrino mixing angles are not
exactly predicted, they can be large, which agrees well with the
current experimental situation.  However it requires the vacuum or
just-so solution to the solar neutrino problem, not currently favored.
One test for the model is the search for anomalous seasonal effects at
Borexino \cite{bxino}.

\subsection{Supersymmetry As the Origin of Neutrino Mass 
  \cite{Hirsch:2000ef,Mira:2000gg}}
\label{sec:supersymm-as-orig}

The suggestion that neutrino masses may be genuinely supersymmetric is
rather old \cite{rpnu-spo1,rpnu-exp}. One simply postulates that they
arise as a result of the violation of R parity (RPV),
$$R_{p}=(-1)^{3B+L+2S}$$
where $S, B, L$ denote spin, baryon and total
lepton number, respectively. Standard model particles (including Higgs
scalars) are RP-even, while their superpartners are odd.  R parity
acts as a selection rule, according to which supersymmetric particles
can only be pair-produced, the lightest \texttt{sparticle} being
stable.  These properties have been taken to be the basis of most
searches of supersymmetric particles in the laboratory. However, R
parity is by no means a sacred symmetry and its violation could
substantially affect the expected SUSY signatures while endowing
neutrinos with mass.

There are several ways to realize this as illustrated in Fig.
\ref{fig:rpcase1}.
\begin{figure}[htbp]
  \centering
\includegraphics[width=8cm,height=5cm]{rpcase1.eps}  
  \caption{Ways to break R parity.}
  \label{fig:rpcase1}
\end{figure}
%
The violation of R--parity can arise explicitly \cite{rpnu-exp} as a
as a residual effect from unification that filters down to low
energies, or simply as an electroweak phenomenon generated due to
non-zero \21 singlet sneutrino vevs~\cite{rpnu-spo2}. The latter may
be realized either with gauged or ungauged lepton number, in the
latter case there is a majoron.

Barring the existence of special symmetries, when R--parity breaks
explicitly there can be both bi-linear and tri-linear terms in the
low-energy superpotential
$$W= W_{MSSM} + \mu_\alpha \ell_\alpha H_u + {\rm trilinears}.$$
%
However there may exist symmetries which simplify the otherwise very
complex structure of explicit RPV terms at low-energy, characterized
by too many new independent coupling constants.
%
An example of this possibility arises within the context of theories
with \texttt{anomalous $U(1)_H$ horizontal symmetries}. It has been
shown~\cite{Mira:2000gg} how one can choose the symmetry in such a way
that all 45 trilinear R-parity violating couplings are forbidden by
holomorphy. The symmetry also provides a common origin for the
supersymmetric $\mu$ term and the solution of the solar and
atmospheric neutrino anomalies, potentially explaining also the
required hierarchy between these.

Alternatively, if R parity is violated spontaneously \cite{rpnu-spo2}
necessarily one obtains the bi-linear RPV model (BRPV) as the
effective low-energy theory. Thus the BRPV model is the ``common
denominator'' of the three most attractive routes in Fig.
\ref{fig:rpcase1}, ``closed'' under the renormalization group. It also
constitutes the simplest and most predictive extension of the MSSM
which renders a very systematic investigation of R parity violation
effects not only in neutrino physics but also for supersymmetry
phenomenology at collider experiments.

On general grounds the spectrum of the BRPV model is expected to be
hierarchical, since only one of the three neutrinos picks up mass at
the tree level, by mixing with neutralinos, as indicated in Fig.
\ref{fig:solatrpsusy} (left panel).
\begin{figure}[htbp]
  \centering
\includegraphics[width=4cm,height=3cm]{solatm0.ps}
\includegraphics[width=4cm,height=3cm]{solatmbrpv1.ps}  
  \caption{Neutrino masses from RPV supersymmetry. 
    The degeneracy of the massless states on the left is lifted by
    supersymmetric loop effects, illustrated in the right panel.}
  \label{fig:solatrpsusy}
\end{figure}
%
Their degeneracy is lifted only by calculable loop effects (right
panel) \cite{Hirsch:2000ef} due to supersymmetric particle
exchanges~\cite{Diaz:1998xc}.
%
Both the tree level atmospheric scale as well as the more complex
loop-induced solar mass scale depend quadratically on RPV parameters.
%
On the other hand neutrino mixing angles are given as RPV parameter
ratios.  For example the left panel in Fig. \ref{fig:atmtest} shows
the correlation between $\tan^2_{\theta_A}$ with the
$\Lambda_2/\Lambda_3$ ratio. Similar correlations exist expressing the
solar and reactor angles as suitable RPV ratios \cite{Hirsch:2000ef}.

In summary, the BRPV model provides a powerful dynamical scheme for
neutrino masses which can account for the observed atmospheric and
solar neutrino anomalies in terms large angle neutrino oscillations,
as indicated by current data.


\section{Neutrino Mixing Tests at Colliders 
\cite{Porod:2000hv,Restrepo:2001me}}
\label{sec:neutr-mixing-tests}

Despite the smallness of neutrino masses indicated by the oscillation
interpretation of current neutrino data, R-parity violation can be
observable at high--energy collider experiments such as the
LHC~\cite{Bartl:2001yh}.
%
The first and most direct manifestation of RPV is the decay of the
lightest supersymmetric particle (LSP).  In many supergravity
scenaria, such as the BRPV model, the LSP is typically the lightest
neutralino and its mass could potentially be explored at LHC and other
planned colliders.
%
Since R parity is violated the LSP decays as
$$\chi^0 \to f_1 + \bar{f_2} + \ell$$
where $f_i$ denotes a quark,
lepton or neutrino and $\ell$ denotes either a charged lepton or a
neutrino.
%
\begin{figure}[htbp]
  \centering
%
\includegraphics[width=6cm,height=4.5cm]{lsp.ps}  
  \caption{LSP decay length (in cm) vs its mass (in GeV) in the BRPV model,
    from \cite{Bartl:2001yh}.}
  \label{fig:lspdecaylength}
\end{figure}
%
In Fig. \ref{fig:lspdecaylength} we show the typical LSP decay length
expected in the BRPV model. Clearly the decay length is short enough
to be studied at the LHC. Moreover the decay products are mostly
visible, with the main decay channel being semi-leptonic, into two
b-jets plus missing transverse momentum.

Most importantly, the BRPV model provides an unambiguous test of the
mixing angles involved in the neutrino anomalies at accelerators, by
measuring the decay patterns of the lightest neutralino
\cite{Porod:2000hv}. For example, in the right panel of Fig.
\ref{fig:atmtest} I illustrate how the ratio of semileptonic
neutralino decay branching ratios to muons and tau leptons is
\texttt{predicted} by the measured value of the atmospheric mixing
angle indicated by atmospheric data.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=4cm,height=3cm]{TanAtm2Lam23.eps}
\includegraphics[width=4cm,height=3cm]{BrMuQQpTauQQpTanAtmsq.eps}  
  \caption{Testing the atmospheric angle in LSP decays, 
from \cite{Porod:2000hv}.}
  \label{fig:atmtest}
\end{figure}
%
Similar considerations hold if the LSP in not the lightest neutralino,
but some other supersymmetric particles. For example when the stop
quark is the LSP, similar correlations exist between the stop decay
branching ratios and the solar angle \cite{Restrepo:2001me}.


%
\bigskip
%

\section{Life Beyond Oscillations? }
\label{sec:life-beyond-oscill}

Massive neutrinos typically have Non-Standard Interactions (NSI),
which can be both flavour-changing (FC) and non-universal (NU)
\cite{Valle:1990pk}.

\texttt{Gauge-type NSI} are defined as coming from a nontrivial
structure of CC and NC weak interactions. They arise in seesaw-type
schemes of neutrino mass due to the rectangular (effectively
non-unitary) nature of the charged current lepton mixing matrix and
correspondingly non-trivial neutral current matrix (non-diagonal in
mass eigenstate neutrinos) \cite{Schechter:1980gr}.  In some models
these NSI may be sizable and lead to flavor and CP violation, even
with strictly degenerate massless neutrinos \cite{NSImodels2}.
%
On the other hand FC-NSI may also be induced by the exchange of
spinless bosons, a situation that arises in radiative models of
neutrino mass \cite{Zee:1980ai} as well as low-energy SUSY models with
broken R parity \cite{rpnu-spo1,rpnu-exp,rpnu-spo2}.  They also arise
in some supersymmetric unified models \cite{NSImodels3}, where they
may be calculable by the renormalization group evolution. Finally,
majoron models of neutrino mass lead to pseudoscalar NSI which may
induce invisible neutrino decays \cite{Schechter:1981cv,V}.

In addition to such renormalizable NSI, gauge theories of neutrino
mass may also lead to NSI of higher dimension $\geq 5$, such as
Majorana neutrino transition magnetic moments which lead to the
Spin-Flavour Precession phenomenon~\cite{Schechter:1981hw}.  As we
will see shortly, this may affect solar neutrino propagation in an
important way~\cite{RSFP}, providing two solutions to the solar
neutrino puzzle, one of which is resonant \cite{Miranda:2000bi} while
the other is non-resonant~\cite{Miranda:2001hv}.

In what follows I take an effective model-independent description of
NSI in which they are given as dimension-6 terms of the type
$\varepsilon G_F$, as illustrated in Fig. \ref{fig:nuNSI},
%
\begin{figure}[htbp]
  \centering
\includegraphics[scale=.24]{fc.ps}  
  \caption{Effective description of neutrino non-standard interactions.}
  \label{fig:nuNSI}
\end{figure}
%
where $\varepsilon$ specifies their sub-weak strength. Such NSI may
affect the propagation of neutrinos in a variety of situations.
%
For example they have been considered both in context of the
solar~\cite{Valle:1987gv,Guzzo:1991cp,Bergmann:2000gp,Guzzo:2001mi}
and atmospheric ~\cite{Fornengo:2001pm,Gonzalez-Garcia:1998hj}
neutrino problems, as well as in connection with
astrophysics~\cite{Nunokawa:1996tg,Grasso:1998tt,Fogli:2002xj}.

In the light of such NSI one may ask, how robust are oscillations ?
To make a long story short, the answer is simple, atmospheric
oscillations are robust, solar are not. Let me start with the
atmospheric ones.

\section{Flavor Changing (FC) Neutrino-Matter Interactions}
\label{sec:flavor-changing-fc}

\subsection{Atmospheric \cite{Fornengo:2001pm}}
\label{sec:atmospheric}

Flavor-changing neutrino-matter interactions in the \nm-\nt channel
have been shown to account for the zenith-angle-dependent deficit of
atmospheric neutrinos observed in the SuperKamiokande experiment
contained events \cite{Gonzalez-Garcia:1998hj,Fornengo:1999zp},
without directly invoking neither neutrino mass, nor mixing.

However such FC explanation fails to reconcile the Super-Kamiokande
contained events with Super-Kamiokande and MACRO up-going muons, due
to the lack of energy dependence intrinsic of NSI conversions
\cite{Fogli:2000ak}.
%
The most recent global analysis of the atmospheric neutrino anomaly in
terms of NSI shows that a pure NSI mechanism is now ruled out at 99\%,
while the standard \nm $\to$ \nt oscillations provide a remarkably
good description of the anomaly \cite{Fornengo:2001pm}.
%
The robustness of the atmospheric \nm $\to$ \nt oscillation hypothesis
was shown to provide the most stringent and model-independent limits
on flavour-changing (FC) and non-universal (NU) neutrino interactions,
as illustrated in Fig. \ref{fig:atmnsibds}.
\begin{figure}[htbp]
  \centering
\includegraphics[width=8cm,height=6cm]{fig.chi2_eps.eps}  
  \caption{Atmospheric bounds on non-standard neutrino 
    interactions, from \cite{Fornengo:2001pm}.}
  \label{fig:atmnsibds}
\end{figure}
%
These limits are model independent as they are obtained from pure
neutrino-physics processes. The stability of the neutrino oscillation
solution to the atmospheric neutrino anomaly against the presence of
non-standard neutrino interactions illustrates the robustness of the
near-maximal atmospheric mixing and massive-neutrino hypothesis.


\subsection{Neutrino Factories-Atmospheric \cite{Huber:2001zw}}
\label{sec:nufact}

The proposed neutrino factories aim at probing the lepton mixing angle
$\theta_{13}$ with much better sensitivity \cite{Freund:2001ui} than
possible today by combining reactor and atmospheric data
\cite{Fornengo:2001pm}. Given the LMA solution of the solar neutrino
problem, it follows that there is a chance that indeed neutrino
factories (NUFACT, for short) will also probe the concurrent leptonic
CP violating effects (for example, through the measurement of CP
asymmetries).

The potential of such a generic NUFACT in probing non-standard
neutrino-matter interactions has been determined in ref.
\cite{Huber:2001zw}.  It has been found that the sensitivity to
flavour-changing (FC) NSI can be substantially improved with respect
to present atmospheric neutrino data, especially at energies higher
than approximately 50 GeV, where the effect of the tau mass is small,
as illustrated in Fig. \ref{fig:nufacnsibds}.  For example, a 100 GeV
NUFACT can probe FC neutrino interactions at the level of few
$|\epsilon| < {\rm few} \times 10^{-4}$ at 99 \% C.L.
%
\begin{figure}[htbp]
  \centering
\includegraphics[scale=.65]{sensMTlow.eps}
  \caption{NUFACT sensitivities to non-standard neutrino 
    interactions, from \cite{Huber:2001zw}.}
  \label{fig:nufacnsibds}
\end{figure}

\subsection{Solar \cite{Guzzo:2001mi}}
\label{sec:solar}

Let us now turn to the solar neutrino oscillation hypothesis.  It
provides an excellent fit to the data and a powerful way to determine
neutrino mass and mixing. But, is it solid?  The lack of clear hints
of spectral distortion, day-night or seasonal variation in the solar
neutrino data implies that, despite the increasing weight of such
rate-independent observables, they do not give an unambiguous
smoking-gun signal for any particular solar neutrino conversion
mechanism.  Here I examine two viable alternatives which do not
require neutrinos to be mixed.

%
I will show how one can have hybrid three-neutrino interpretations of
the current neutrino data in which the atmospheric data are accounted
for (mainly) by oscillations, while the solar neutrino data are
understood in terms of some sort of NSI, as illustrated in Fig.
\ref{fig:hyb}.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=3.6cm,height=3.6cm]{solatm00.ps}
  \caption{Hybrid three-neutrino interpretations of current neutrino data.
    The figure illustrates solutions to the solar and atmospheric
    neutrino anomalies in which there is only one oscillation, the
    atmospheric one, while the solar conversions are due to
    non-standard interactions (NSI).}
  \label{fig:hyb}
\end{figure}
%

FC-NSI solutions to the solar neutrino problem have been considered in
detail in ref.~\cite{Guzzo:2001mi}. 
%
Since the NSI solution is energy-independent the spectrum is
undistorted, the global analysis observables are the solar neutrino
rates in all experiments and the day-night measurements.
%
It was found that the NSI description of solar data is slightly better
than that of the oscillation solution and that the allowed NSI regions
are determined mainly by the rate analysis \cite{Guzzo:2001mi}.
%
By using simplified ans\"atze for the NSI interactions it was
explicitly demonstrated that the NSI values indicated by the solar
data analysis are fully acceptable also for the atmospheric data.
%
The required parameters are indicated in Fig. \ref{fig:hybparam}.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=8cm,height=7cm]{bm.fig.quark-U.eps}  
  \caption{Parameters of NSI solution to the solar neutrino anomaly, 
from \cite{Guzzo:2001mi}.}
  \label{fig:hybparam}
\end{figure}

Since within such a three-neutrino scheme the solar and atmospheric
neutrino sectors are connected not only by the neutrino mixing angle
$\theta_{13}$ but also by the flavour-changing parameters accounting
for the solar data, it follows that to some extent NSI and
oscillations may be confused.

\subsection{Neutrino Factories-Solar \cite{Huber:2001de,Huber:2002bi}}
\label{sec:neutrino-factories}

The impact of NSI on the determination of neutrino mixing parameters
at a neutrino factory using the so-called ``golden channels'' \ne to
\nm and \bne to \bnm for the measurement of $\theta_{13}$ has been
given in ref. \cite{Huber:2001de,Huber:2002bi}.  It was shown how a
certain combination of FC interactions in neutrino source and earth
matter can give exactly the same signal as oscillations arising due to
$\theta_{13}$. This implies that information on $\theta_{13}$ can only
be obtained if bounds on NSI are available.  Taking into account the
existing bounds on FC interactions, this leads to a drastic loss in
sensitivity in $\theta_{13}$, of at least two orders of magnitude. A
near detector at a neutrino factory offers the possibility to obtain
stringent bounds on some NSI parameters. Such near-site detector
constitutes an essential ingredient of a neutrino factory and a
necessary step towards the determination of $\theta_{13}$ and
subsequent study of leptonic CP violation.


Additional signatures of theories leading to FC interactions would be
the existence of sizable flavour non-conservation effects, such as
$\mu \to e + \gamma$, $\mu-e$ conversion in nuclei, even in the
massless neutrino limit~\cite{Valle:1990pk}.

\section{Solar Spin-Flavour Neutrino Precession 
\cite{Miranda:2000bi,Miranda:2001hv}}
\label{sec:solar-spin-flavour}

The Spin-Flavour Precession (SFP) of neutrinos \cite{Schechter:1981hw}
can be resonant in matter \cite{RSFP} and this leads to a solution of
the solar neutrino anomaly. The magnetic field profile may be fixed
self-consistently by physical requirements and by demanding it to
solve the solar magneto-hydrodynamics (MHD) equations
\cite{Miranda:2000bi}. The resulting \ne survival probability is
displayed in Fig.  \ref{fig:NMMprob} and compared with the one
expected in the case of neutrino oscillation.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=7cm,height=3.6cm]{NMMprob2.eps}  
  \caption{Best fit oscillation vs spin 
    flavor precession \ne survival probability, from
    \cite{Miranda:2001hv}.}
  \label{fig:NMMprob}
\end{figure}
%
There are two solutions: resonant (RSFP) and non-resonant (NRSFP).
Clearly both provide good global descriptions of the solar neutrino
data. 
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=5.5cm,height=6cm]{chi_4rggsno_b84_nebp00.ps}
\includegraphics[width=5.5cm,height=6cm]{spec_1258_histo.ps}  
  \caption{``Grand Unified'' solar spin-flavor-precession/oscillation 
    plot and resulting recoil electron spectra, including the
    unacceptable SMA spectrum, from \cite{Miranda:2001hv}.}
  \label{fig:sfparam}
\end{figure}
%
In the left panel of Fig. \ref{fig:sfparam} the SFP parameters that
account for the solar neutrino problem are presented (the plot
includes the first SNO CC result, as well as the 1258--day
Super-Kamiokande data). In the MHD approach the resulting scheme has
only three effective parameters: $\Delta m^2$, $\mu B_\perp$ and the
neutrino mixing angle $\theta$.  The rates-only analysis for fixed
$\mu B_\perp$ slightly favors SFP solutions over oscillations. In
addition to the resonant solution (RSFP), there is a non-resonant
solution (NRSFP) in the ``dark-side''\cite{Miranda:2001hv}.  Note that
in the presence of a neutrino transition magnetic moment of $10^{-11}$
Bohr magneton, a magnetic field of 80 KGauss eliminates all large
mixing solutions other than LMA, as they would violate the upper limit
on anti-neutrinos from the Sun \cite{Aglietta:zu}.
%
The right panel in Fig. \ref{fig:sfparam} shows that the predicted
LMA, RSFP and NRSFP recoil energy spectra are in excellent agreement
with the data, in contrast to that of SMA.
%
Finally I note that it has recently been shown that the SFP solution
of the solar neutrino problem can be distinguished from the currently
favored oscillation solutions at Borexino \cite{Akhmedov:2002ti}.

Before closing let me also add that the possibility of fast neutrino
decays~\cite{V} has also been advocated could also play a sub-leading
role in the solar~\cite{Bahcall:1986gq} and/or atmospheric neutrino
problems~\cite{atmdecay}. The data allow for neutrino stability tests,
though not specially stringent.



\section{In a Nutshell}

It is now very hard to dispute the need for physics beyond the
Standard Model in order to explain current solar and atmospheric data.
%
In order to also account for the LSND hint \cite{lsnd} indicating \bnm
to \bne oscillations with a eV-mass-squared difference in the
framework of quantum field theory one needs an additional light
sterile neutrino \ns \cite{giunti,Caldwell:kn,ptv92,pv93}.
%
Currently LSND is neither confirmed nor ruled out by any other
experiment.
%
A unified global analysis of current neutrino oscillation data for
4-neutrino mass schemes has been given in \cite{Maltoni:2001bc} using
data from solar and atmospheric neutrino experiments, as well as
information from short-baseline experiments including LSND. 
%
The summary of current 4-neutrino oscillation parameters is given in
Fig. \ref{fig:4nuoscplot}.
%
\begin{figure}[htbp]
  \centering
\includegraphics[width=10cm,height=6.5cm]{bm.fournu_loc.eps}  
  \caption{``Grand-unified'' neutrino oscillation plot, adapted 
    from ref.  \cite{Maltoni:2001bc}, before
    \cite{Ahmad:2002jz,Ahmad:2002ka}.}
  \label{fig:4nuoscplot}
\vglue -.3cm
\end{figure}
%
One finds that the solar and atmospheric data reject sterile
neutrinos.  For example, recent SNO data rule out conversion to
\texttt{ pure} sterile neutrino states at the $\sim 5\sigma$ level. A
smaller degree of sterile neutrino rejection is found for atmospheric
data \cite{Maltoni:2001bc}.  For this reason, the best 2+2 model model
has now changed from one with the sterile state mainly in the solar
sector \cite{pv93} to one in which it lies mainly in the atmospheric
sector \cite{ptv92}.  However I stress that, in a truly global fit
including also the LSND data \cite{Maltoni:2001bc}, sterile neutrinos
can not be ruled out, since the latter disfavor the Standard Model at
99.9\% CL with respect to the best 4-neutrino model. The best-fit
four-neutrino model is a 3+1 scheme, since in this case one can
decouple \ns from both solar and atmospheric sectors.  However, in
such a global sense, even 2+2 models of the type considered in ref.
\cite{ptv92} can not yet be ruled out, being still allowed at the 90\%
CL \cite{solat02}.  For the moment the sterile neutrino oscillation
hypothesis is therefore alive, though there is a clear tension between
underground and short-baseline data.  We wait anxiously for results
from MiniBooNE to settle the issue~\cite{MiniBooNE} in the near
future.


%
\par

I thank Martin Hirsch for comments on the manuscript.  This work was
supported by the European Commission grants HPRN-CT-2000-00148 and
HPMT-2000-00124, by the ESF \emph{Neutrino Astrophysics Network} and
by Spanish MCyT grant PB98-0693.

\begin{thebibliography}{99}
  
\bibitem{sol} S.~Fukuda {\it et al.}  [SuperKamiokande Collaboration],
%``Solar B-8 and he p neutrino measurements from 1258 days of  Super-Kamiokande data,''
  Phys.\ Rev.\ Lett.\ {\bf 86} (2001) 5651 %%.
%%CITATION = ;%%
 http://www-sk.icrr.u-tokyo.ac.jp/doc/sk/
%
 B. T. Cleveland {\it et al.}, Astrophys. J. {\bf 496}, 505 (1998); R.
 Davis, Prog. Part. Nucl. Phys. {\bf 32}, 13 (1994); 
  %
 SAGE Coll., J. N. Abdurashitov {\it et al.}, Phys. Rev. {\bf C60},
 055801 (1999); http://ewiserver.npl.washington.edu/SAGE/SAGE.html
  %
 GALLEX Coll., W.~Hampel {\it et al.}, Phys. Lett. {\bf B447}, 127
 (1999); K. Lande and V.  Gavrin and E. Belloti (GNO Coll.), Talks
 at \textit{Neutrino 2000}, 15--21 June 2000, Sudbury, Canada,
 (\textit{http://{\-}nu2000.{\-}sno.{\-}laurentian.{\-}ca}).
  %
  Q.~R.~Ahmad {\it et al.}  [SNO Collaboration],
%``Measurement of the charged of current interactions produced by B-8  solar neutrinos at the Sudbury Neutrino Observatory,''
  Phys.\ Rev.\ Lett.\ {\bf 87} (2001) 071301 %%.
%%CITATION = ;%%

\bibitem{ssm}
Bahcall's page http://www.sns.ias.edu/~jnb/ and references therein.

\bibitem{Kachelriess:2001sg}
M.~Kachelriess, A.~Strumia, R.~Tomas and J.~W.~F.~Valle,
%``SN1987A and the status of oscillation solutions to the solar neutrino  problem,''
Phys.\ Rev.\ D {\bf 65}, 073016 (2002) %.
%%CITATION = ;%%

\bibitem{Gonzalez-Garcia:1999aj}
M.~C.~Gonzalez-Garcia, P.~C.~de Holanda, C.~Pena-Garay and J.~W.~F.~Valle,
%``Status of the MSW solutions of the solar neutrino problem,''
Nucl.\ Phys.\ B {\bf 573} (2000) 3 .
%%CITATION = ;%%

\bibitem{Ahmad:2002jz}
Q.~R.~Ahmad {\it et al.}  [SNO Collaboration],
%``Direct evidence for neutrino flavor transformation from neutral-current  interactions in the Sudbury Neutrino Observatory,''
arXiv:.
%%CITATION = ;%%

\bibitem{Ahmad:2002ka}
Q.~R.~Ahmad {\it et al.}  [SNO Collaboration],
%``Measurement of day and night neutrino energy spectra at SNO and  constraints on neutrino mixing parameters,''
arXiv:.
%%CITATION = ;%%

\bibitem{solat02}
M. Maltoni et al, in preparation

\bibitem{KamLAND}
KamLAND page http://www.awa.tohoku.ac.jp/html/KamLAND/

\bibitem{atm}
  Y. Fukuda \textit{et al.},
  Kamiokande Coll., Phys. Lett. B \textbf{335} (1994) 237;
  %
  R. Becker-Szendy \textit{et al.},
  IMB Coll., Nucl. Phys. B (Proc. Suppl.) \textbf{38} (1995) 331;
  %
  W.W.M. Allison \textit{et al.}, Soudan Coll.,
  Phys. Lett. B \textbf{449} (1999) 137.
%
  Super-Kamiokande Coll., Y. Fukuda \textit{et al.}, Phys. Rev. Lett.
  \textbf{81} (1998) 1562; C. McGrew in {\it Neutrino Telescopes
    2001}, Venice, Italy, March 2001; T. Toshito in {\it
    Moriond 2001}, Les Arcs, France, March 2001.
%
  MACRO Coll., M. Ambrosio \textit{et al.}, Phys. Lett. B \textbf{434}
  (1998) 451; M. Spurio {\it et al.}, ; B. Barish, Talk
  given at \textit{Neutrino 2000}, 15--21 June 2000, Sudbury, Canada,
  \verb!http://nu2000.sno.laurentian.ca}!.
%
A.~Surdo, Talk given at \textit{TAUP 2001}, 8--12 September 2001, Gran
Sasso, Italy (\textit{http://www.lngs.infn.it/}).

\bibitem{atmfluxes} 
G.~Battistoni et al
%``Progresses in the validation of the FLUKA atmospheric nu flux  calculations,''
Nucl.\ Phys.\ Proc.\ Suppl.\ {\bf 110} (2002) 336
%.
%%CITATION = ;%%

%%%

\bibitem{Gonzalez-Garcia:2000sq}
M.~C.~Gonzalez-Garcia, M.~Maltoni, C.~Pena-Garay and J.~W.~F.~Valle,
%``Global three-neutrino oscillation analysis of neutrino data,''
Phys.\ Rev.\ D {\bf 63} (2001) 033005 .
%%CITATION = ;%%

\bibitem{Gonzalez-Garcia:2002mu}
M.~C.~Gonzalez-Garcia and M.~Maltoni,
%``Two mass-scale oscillation analysis of atmospheric and reactor data,''
arXiv:.
%%CITATION = ;%%

\bibitem{Fornengo:2000sr}
N.~Fornengo, M.~C.~Gonzalez-Garcia and J.~W.~F.~Valle,
%``Updated global analysis of the atmospheric neutrino data in terms of  neutrino oscillations,''
Nucl.\ Phys.\ B {\bf 580} (2000) 58  and references therein.
%%CITATION = ;%%

\bibitem{Fornengo:2001pm}
N.~Fornengo, M.~Maltoni, R.~T.~Bayo and J.~W.~F.~Valle,
%``Probing neutrino non-standard interactions with atmospheric neutrino  data,''
Phys.\ Rev.\ D {\bf 65} (2002) 013010  and
references therein.
%%CITATION = ;%%

\bibitem{Kachelriess:2000fe}
M.~Kachelriess, R.~Tomas and J.~W.~F.~Valle,
%``Large lepton mixing and supernova 1987A,''
JHEP {\bf 0101} (2001) 030 %%.
%%CITATION = ;%%

\bibitem{Minakata:2001cd}
H.~Minakata, H.~Nunokawa, R.~Tomas and J.~W.~F.~Valle,
%``Probing supernova physics with neutrino oscillations,''
arXiv:.
%%CITATION = ;%%

\bibitem{snperm} A.Y.~Smirnov, D.N.~Spergel and J.N.~Bahcall, 
%``Is large lepton mixing excluded?,'' 
Phys.\ Rev.\ {\bf D49}, 1389 (1994)
%%  .
%%CITATION = ;%%
  B.~Jegerlehner, F.~Neubig and G.~Raffelt, 
%``Neutrino oscillations  and the supernova SN1987A signal,'' 
  Phys. Rev. {\bf D54}, 1194 (1996) %%.
%%CITATION = PHRVA,D54,1194;%%

\bibitem{Hirata:1987hu}
K.~Hirata {\it et al.}  [KAMIOKANDE-II Collaboration],
%``Observation Of A Neutrino Burst From The Supernova Sn1987a,''
Phys.\ Rev.\ Lett.\  {\bf 58} (1987) 1490.
%%CITATION = PRLTA,58,1490;%%
R.~M.~Bionta {\it et al.},
%``Observation Of A Neutrino Burst In Coincidence With Supernova Sn1987a In The Large Magellanic Cloud,''
Phys.\ Rev.\ Lett.\  {\bf 58} (1987) 1494.
%%CITATION = PRLTA,58,1494;%%

\bibitem{Kachelriess:2001bs}
M.~Kachelriess and R.~Tomas,
%``Non-adiabatic level crossing in (non-) resonant neutrino oscillations,''
Phys.\ Rev.\ D {\bf 64} (2001) 073002 %.
%%CITATION = ;%%

\bibitem{Weinberg:uk}
S.~Weinberg, ``Opening Address - Neutrinos '78'' 
{\it Proceedings, West Lafayette 1978, 1-21}; R.~Barbieri, J.~R.~Ellis
and M.~K.~Gaillard,
%``Neutrino Masses And Oscillations In SU(5),''
Phys.\ Lett.\ B {\bf 90}, 249 (1980).
%%CITATION = PHLTA,B90,249;%%

\bibitem{seesaw} M.~Gell-Mann, P.~Ramond, R.~Slansky, in {\sl
Supergravity}, ed. P.~van~Niewenhuizen and D.~Freedman (North
Holland, 1979); T.~Yanagida, in {\sl KEK lectures}, ed.  O.~Sawada and
A.~Sugamoto (KEK, 1979) 

\bibitem{Mohapatra:1979ia}
R.~N.~Mohapatra and G.~Senjanovic,
%``Neutrino Mass And Spontaneous Parity Nonconservation,''
Phys.\ Rev.\ Lett.\  {\bf 44} (1980) 91.
%%CITATION = PRLTA,44,912;%%

\bibitem{Schechter:1980gr}
J.~Schechter and J.~W.~F.~Valle,
%``Neutrino Masses In SU(2) X U(1) Theories,''
Phys.\ Rev.\ D {\bf 22} (1980) 2227.
%%CITATION = PHRVA,D22,2227;%%

\bibitem{Schechter:1980gk}
J.~Schechter and J.~W.~F.~Valle,
%``Neutrino Oscillation Thought Experiment,''
Phys.\ Rev.\ D {\bf 23} (1981) 1666.
%%CITATION = PHRVA,D23,1666;%%

\bibitem{Chikashige:1980ui}
Y.~Chikashige, R.~N.~Mohapatra and R.~D.~Peccei,
%``Are There Real Goldstone Bosons Associated With Broken Lepton Number?,''
Phys.\ Lett.\ B {\bf 98} (1981) 265.
%%CITATION = PHLTA,B98,265;%%

\bibitem{Schechter:1981cv}
J.~Schechter and J.~W.~F.~Valle,
%``Neutrino Decay And Spontaneous Violation Of Lepton Number,''
Phys.\ Rev.\ D {\bf 25} (1982) 774.
%%CITATION = PHRVA,D25,774;%%

\bibitem{Schechter:1981bd}
J.~Schechter and J.~W.~F.~Valle,
%``Neutrinoless Double-Beta Decay In SU(2) X U(1) Theories,''
Phys.\ Rev.\ D {\bf 25} (1982) 2951.\\
%%CITATION = PHRVA,D25,2951;%%
M.~Doi, T.~Kotani and E.~Takasugi,
%``Double Beta Decay And Majorana Neutrino,''
Prog.\ Theor.\ Phys.\ Suppl.\  {\bf 83} (1985) 1.
%%CITATION = PTPSA,83,1;%%

\bibitem{Pal:1981rm}
P.~B.~Pal and L.~Wolfenstein,
%``Radiative Decays Of Massive Neutrinos,''
Phys.\ Rev.\ D {\bf 25} (1982) 766.
%%CITATION = PHRVA,D25,766;%%

\bibitem{Schechter:1981hw}
J.~Schechter and J.~W.~F.~Valle,
%``Majorana Neutrinos And Magnetic Fields,''
Phys.\ Rev.\ D {\bf 24} (1981) 1883 [ibid {\bf 25} (1982) 283]
%%CITATION = PHRVA,D24,1883;%%

\bibitem{Wolfenstein:1977ue}
L.~Wolfenstein,
%``Neutrino Oscillations In Matter,''
Phys.\ Rev.\ D {\bf 17}, 2369 (1978).\\
%%CITATION = PHRVA,D17,2369;%%
S.~P.~Mikheev and A.~Y.~Smirnov,
%``Resonance Enhancement Of Oscillations In Matter And Solar Neutrino  Spectroscopy,''
Sov.\ J.\ Nucl.\ Phys.\  {\bf 42} (1985) 913
%%CITATION = SJNCA,42,913;%%

\bibitem{Bilenky:1996rw}
For recent reviews see S.~M.~Bilenky, C.~Giunti and W.~Grimus,
%``Neutrino mass spectrum from the results of neutrino oscillation  experiments,''
Eur.\ Phys.\ J.\ C {\bf 1} (1998) 247
% .
%%CITATION = ;%%
%
M.~C.~Gonzalez-Garcia and Y.~Nir,
%``Developments in neutrino physics,''
arXiv:.
%%CITATION = ;%%

\bibitem{Kobayashi:fv}
M.~Kobayashi and T.~Maskawa,
%``CP Violation In The Renormalizable Theory Of Weak Interaction,''
Prog.\ Theor.\ Phys.\  {\bf 49} (1973) 652.
%%CITATION = PTPKA,49,652;%%

\bibitem{Schechter:1979bn}
J.~Schechter and J.~W.~F.~Valle,
%``Comment On The Lepton Mixing Matrix,''
Phys.\ Rev.\ D {\bf 21} (1980) 309.
%%CITATION = PHRVA,D21,309;%%

\bibitem{nuFacPapers} C.~Quigg, Nucl.\ Instrum.\ Meth.\ A {\bf 451}
  (2000) 1; C.~Albright {\it et al.}, ; B.~Autin {\it et
    al.}, CERN-SPSC-98-30; S.~Geer, Phys.\ Rev.\ D {\bf 57} (1998)
  6989 [Erratum-ibid.\ D {\bf 59} (1998) 039903] %%.

\bibitem{Morales:1998hu}
For a review on \nbb experiments see A.~Morales,
%``Review on double beta decay experiments and comparison with theory,''
Nucl. Phys.\ Proc.\ Suppl.\ {\bf 77}, 335 (1999) .
%%CITATION = ;%%

\bibitem{Klapdor-Kleingrothaus:2001ke}
H.~V.~Klapdor-Kleingrothaus et al,
%``Evidence for neutrinoless double beta decay,''
                  Mod.\ Phys.\ Lett.\ A {\bf 16}, 2409 (2001)
             %%     .
%%CITATION = ;%%
                  
\bibitem{Aalseth:2002dt} 
For a critical view of the \nbb hint see C.~E.~Aalseth {\it et al.},
%``Comment on 'Evidence for neutrinoless double beta decay',''
arXiv: and F.~Feruglio, A.~Strumia and F.~Vissani,
%``Neutrino oscillations and signals in beta and 0nu 2beta experiments,''
arXiv:.
%%CITATION = ;%%

\bibitem{Klapdor-Kleingrothaus:1999hk}
H. V. Klapdor-Kleingrothaus et al,  [GENIUS Coll.],
%``GENIUS: A supersensitive germanium detector system for rare events,''
arXiv:.
%%CITATION = ;%%

\bibitem{Zee:1980ai}
T.~P.~Cheng and L.~F.~Li,
%``Neutrino Masses, Mixings And Oscillations In SU(2) X U(1) Models Of Electroweak Interactions,''
Phys.\ Rev.\ D {\bf 22} (1980) 2860.
%%CITATION = PHRVA,D22,2860;%%
A.~Zee,
%``A Theory Of Lepton Number Violation, Neutrino Majorana Mass, And Oscillation,''
Phys.\ Lett.\ B {\bf 93} (1980) 389 [Err-ibid.\ B {\bf 95} (1980) 461].
%%CITATION = PHLTA,B93,389;%%
K.~S.~Babu,
%``Model Of 'Calculable' Majorana Neutrino Masses,''
Phys.\ Lett.\ B {\bf 203} (1988) 132.
%%CITATION = PHLTA,B203,132;%%

\bibitem{seesawnupapers} 
For recent seesaw papers see 
http://131.169.91.193/cgi-bin/spiface/find/hep/\\www?rawcmd=t+seesaw
+and+tp+neutrino+and+d+after+1997\&FORMAT=\\WWWBRIEF\&SEQUENCE=

\bibitem{giunti}
See Carlo Giunti's page http://www.to.infn.it/\~{ }giunti/neutrino/

\bibitem{Chankowski:2000fp}
P.~H.~Chankowski et al, 
%``Neutrino unification,''
Phys.\ Rev.\ Lett.\ {\bf 86} (2001) 3488 .
%%CITATION = ;%%

\bibitem{Osipowicz:2001sq}
A.~Osipowicz {\it et al.}  [KATRIN Collaboration],
%``KATRIN: A next generation tritium beta decay experiment with sub-eV  sensitivity for the electron neutrino mass,''
arXiv:.
%%CITATION = ;%%

\bibitem{Elgaroy:2002bi}
O.~Elgaroy {\it et al.},
%``A new limit on the total neutrino mass from the 2dF galaxy redshift  survey,''
arXiv:.
%%CITATION = ;%%

\bibitem{Caldwell:kn}
D.~O.~Caldwell and R.~N.~Mohapatra,
%``Neutrino Mass Explanations Of Solar And Atmospheric Neutrino Deficits  And Hot Dark Matter,''
Phys.\ Rev.\ D {\bf 48} (1993) 259.
%%CITATION = PHRVA,D48,3259;%%

\bibitem{DEG} See e.g., A.  Ioannisian, J.~W.~F. Valle, Phys.\ Lett.\ 
  B {\bf 332} (1994) 93; B. Bamert, C.P.  Burgess, Phys.\ Lett.\ B
  {\bf 329} (1994) 289; D. Caldwell and R. N. Mohapatra, Phys.\ Rev.\ 
  D {\bf 50} (1994) 3477; D. G. Lee and R. N. Mohapatra, Phys.\ Lett.\ 
  B {\bf 329} (1994) 463; A. S. Joshipura, Z.\ Phys.\ C {\bf 64}, 31
  (1994).

\bibitem{DEGNAT}
J.~A.~Casas, J.~R.~Espinosa, A.~Ibarra and I.~Navarro,
%``Nearly degenerate neutrinos, supersymmetry and radiative corrections,''
Nucl.\ Phys.\ B {\bf 569} (2000) 82
J.~Ellis and S.~Lola,
%``Can neutrinos be degenerate in mass?,''
Phys.\ Lett.\ B {\bf 458} (1999) 310
N.~Haba, Y.~Matsui, N.~Okamura and M.~Sugiura,
%``The effect of Majorana phase in degenerate neutrinos,''
Prog.\ Theor.\ Phys.\  {\bf 103} (2000) 145;
E.~Ma,
%``Nearly mass-degenerate Majorana neutrinos: Double beta decay and  neutrino oscillations,''
.

\bibitem{Wolfenstein:1981rk}
L.~Wolfenstein,
%``CP Properties Of Majorana Neutrinos And Double Beta Decay,''
Phys.\ Lett.\ B {\bf 107} (1981) 77.
%%CITATION = PHLTA,B107,77;%%

\bibitem{Valle:1982yw}
J.~W.~F.~Valle,
%``Neutrinoless Double-Beta Decay With Quasi-Dirac Neutrinos,''
Phys.\ Rev.\ D {\bf 27} (1983) 1672.
%%CITATION = PHRVA,D27,1672;%%

\bibitem{deGouvea:2001jp}
A.~de Gouvea and J.~W.~F.~Valle,
%``Minimalistic neutrino mass model,''
Phys.\ Lett.\ B {\bf 501} (2001) 115 .

\bibitem{bxino}
Borexino page, http://almime.mi.infn.it/

\bibitem{Hirsch:2000ef}
M.~Hirsch et al, %M.~A.~Diaz, W.~Porod, J.~C.~Romao and J.~W.~F.~Valle,
%``Neutrino masses and mixings from supersymmetry with bilinear R-parity  violation: A theory for solar and atmospheric neutrino oscillations,''
Phys.\ Rev.\ D {\bf 62} (2000) 113008 .
%%CITATION = ;%%
%M.~Hirsch, W.~Porod, M.~A.~Diaz, J.~C.~Romao and J.~W.~F.~Valle,
%``Neutrino masses and mixings in SUSY with broken R-parity,''
and arXiv:.
%%CITATION = ;%%
J.~C.~Romao et al, %M.~A.~Diaz, M.~Hirsch, W.~Porod and J.~W.~F.~Valle,
%``A supersymmetric solution to the solar and atmospheric neutrino  problems,''
Phys.\ Rev.\ D {\bf 61} (2000) 071703 %%.
%%CITATION = ;%%

\bibitem{Mira:2000gg}
J.~M.~Mira et al,
%``Bilinear R-parity violation and small neutrino masses:  A self-consistent framework,''
Phys.\ Lett.\ B {\bf 492} (2000) 81 .
%%CITATION = ;%%

\bibitem{rpnu-spo1} C.~S.~Aulakh and R.~N.~Mohapatra,
%``Neutrino As The Supersymmetric Partner Of The Majoron,''
Phys.\ Lett.\ B {\bf 119}, 136 (1982).
%%CITATION = PHLTA,B119,136;%%
G.~G.~Ross and J.~W.~F.~Valle,
%``Supersymmetric Models Without R Parity,''
Phys.\ Lett.\ B {\bf 151} (1985) 375.
%%CITATION = PHLTA,B151,375;%%
J.~R.~Ellis, G.~Gelmini, C.~Jarlskog, G.~G.~Ross and J.~W.~F.~Valle,
%``Phenomenology Of Supersymmetry With Broken R Parity,''
Phys.\ Lett.\ B {\bf 150} (1985) 142.
%%CITATION = PHLTA,B150,142;%%
A.~Santamaria and J.~W.~F.~Valle,
%``Spontaneous R Parity Violation In Supersymmetry: A Model For Solar Neutrino Oscillations,''
Phys.\ Lett.\ B {\bf 195}, 423 (1987).
%%CITATION = PHLTA,B195,423;%%
%
%``Supersymmetric Majoron Signatures And Solar Neutrino Oscillations,''
Phys.\ Rev.\ Lett.\  {\bf 60}, 397 (1988).
%%CITATION = PRLTA,60,397;%%
%
%``Solar Neutrino Oscillation Parameters And The Broken R Parity Majoron,''
Phys.\ Rev.\ D {\bf 39}, 1780 (1989).
%%CITATION = PHRVA,D39,1780;%%

\bibitem{rpnu-exp} 
L.~J.~Hall and M.~Suzuki,
%``Explicit R Parity Breaking In Supersymmetric Models,''
Nucl.\ Phys.\ B {\bf 231} (1984) 419.
%%CITATION = NUPHA,B231,419;%%
For more recent references see 
http://131.169.91.193/cgi-bin/spiface/find/hep/www?rawcmd=\\
find+tp+r+parity+and+violation+and+neutrino+and+not+cn+%3EA+and+topcite+50%2B+&FORMAT=WWWBRIEF&SEQUENCE=

\bibitem{rpnu-spo2}
A.~Masiero and J.~W.~F.~Valle,
%``A Model For Spontaneous R Parity Breaking,''
Phys.\ Lett.\ B {\bf 251} (1990) 273.
%%CITATION = PHLTA,B251,273;%%
J.~C.~Romao, C.~A.~Santos and J.~W.~F.~Valle,
%``How To Spontaneously Break R-Parity,''
Phys.\ Lett.\ B {\bf 288} (1992) 311.
%%CITATION = PHLTA,B288,311;%%
%
J.~C.~Romao and J.~W.~F.~Valle,
%``Neutrino masses in supersymmetry with spontaneously broken R parity,''
  Nucl.\ Phys.\ B {\bf 381} (1992) 87

\bibitem{Diaz:1998xc} M.~A.~Diaz, J.~C.~Romao and J.~W.~F.~Valle,
%``Minimal supergravity with R-parity breaking,''
  Nucl.\ Phys.\ B {\bf 524} (1998) 23  and references
  therein.

\bibitem{Porod:2000hv}
W.~Porod, M.~Hirsch, J.~Romao and J.~W.~F.~Valle,
%``Testing neutrino mixing at future collider experiments,''
Phys.\ Rev.\ D {\bf 63} (2001) 115004 %%.
%%CITATION = ;%%

\bibitem{Restrepo:2001me}
D.~Restrepo, W.~Porod and J.~W.~F.~Valle,
%``Broken R-parity, stop decays, and neutrino physics,''
Phys.\ Rev.\ D {\bf 64} (2001) 055011 %%.
%%CITATION = ;%%

\bibitem{Bartl:2001yh}
A.~Bartl et al,
%``Neutralino phenomenology at LEP2 in supersymmetry with bilinear  breaking of R-parity,''
Nucl.\ Phys.\ B {\bf 600} (2001) 39 .

\bibitem{Valle:1990pk} 
For a review see, e.g. J.~W.~F.~Valle,
%``Gauge Theories And The Physics Of Neutrino Mass,''
Prog.\ Part.\ Nucl.\ Phys.\  {\bf 26} (1991) 91.
%%CITATION = PPNPD,26,91;%%

\bibitem{NSImodels2}
R.~N.~Mohapatra and J.~W.~F.~Valle,
%``Neutrino Mass And Baryon-Number Nonconservation In Superstring Models,''
Phys.\ Rev.\ D {\bf 34} (1986) 1642.
%%CITATION = PHRVA,D34,1642;%%
D.~Wyler and L.~Wolfenstein,
%``Massless Neutrinos In Left-Right Symmetric Models,''
Nucl.\ Phys.\ B {\bf 218} (1983) 205.
%%CITATION = NUPHA,B218,205;%%

\bibitem{NSImodels3}
  L.\ J.\ Hall, V.\ A.\ Kostelecky and S.\ Raby,
  Nucl.\ Phys.\ B {\bf 267} (1986) 415
  %%CITATION = NUPHA,B267,415;%%

\bibitem{V} J. W. F. Valle, Phys.\ Lett.\ B {\bf 131} (1983) 87; G.
  Gelmini, J. W. F.  Valle, Phys.\ Lett.\ B {\bf 142} (1984) 181; J.
  W. F. Valle, Phys.\ Lett.\ B {\bf 159} (1985) {49}; A.  Joshipura,
  S. Rindani, Phys. Rev. D {\bf 46} (1992) 3000
  
\bibitem{RSFP}
E.~K.~Akhmedov,
%``Resonant Amplification Of Neutrino Spin Rotation In Matter And The  Solar-Neutrino Problem,''
Phys.\ Lett.\ B {\bf 213} (1988) 64.
%%CITATION = PHLTA,B213,64;%%
C.~S.~Lim and W.~J.~Marciano,
%``Resonant Spin-Flavor Precession Of Solar And Supernova Neutrinos,''
Phys.\ Rev.\ D {\bf 37}, 1368 (1988).
%%CITATION = PHRVA,D37,1368;%%

\bibitem{Miranda:2000bi}
O.~G.~Miranda et al,
%``The simplest resonant spin-flavour solution to the solar neutrino  problem,''
Nucl.\ Phys.\ B {\bf 595} (2001) 360 .
%%CITATION = ;%%

\bibitem{Miranda:2001hv}
O.~G.~Miranda et al,
%``A non-resonant dark-side solution to the solar neutrino problem,''
Phys.\ Lett.\ B {\bf 521} (2001) 299 .
%%CITATION = ;%%

\bibitem{Valle:1987gv}
J.W.~F.~Valle,
%``Resonant Oscillations Of Massless Neutrinos In Matter,''
Phys.\ Lett.\ {\bf 199B} (1987) 432.
 
\bibitem{Guzzo:1991cp}
M.M.~Guzzo et al ,
%``On the matter enhanced transitions of solar neutrinos in the absence of neutrino mixing in vacuum,''
Phys.\ Lett.\ {\bf B271} (1991) 17.

\bibitem{Bergmann:2000gp} S.~Bergmann et al, Phys.\ Rev.\ {\bf D62}
  (2000) 073001 .
  
\bibitem{Guzzo:2001mi}
M.~Guzzo et al, 
%``Status of a hybrid three-neutrino interpretation of neutrino data,''
Nucl.\ Phys.\ B {\bf 629} (2002) 479 .
%%CITATION = ;%%

\bibitem{Gonzalez-Garcia:1998hj} M.C.~Gonzalez-Garcia {\it et al.},
  Phys. Rev. Lett. {\bf 82} (1999) 3202, ; N.~Fornengo,
  M.~C.~Gonzalez-Garcia and J.~W.~F.~Valle,
%``On the interpretation of the atmospheric neutrino data in terms of  flavor changing neutrino interactions,''
  JHEP {\bf 0007} (2000) 006 %%.

\bibitem{Nunokawa:1996tg}
H.~Nunokawa, Y.Z.~Qian, A.~Rossi and J.W.~F.~Valle,
%``Resonant conversion of massless neutrinos in supernovae,''
Phys.\ Rev.\ {\bf D54} (1996) 4356.

\bibitem{Grasso:1998tt}
D.~Grasso, H.~Nunokawa and J.~W.~F.~Valle,
%``Pulsar velocities without neutrino mass,''
Phys.\ Rev.\ Lett.\  {\bf 81} (1998) 2412 %%.
%%CITATION = ;%%

\bibitem{Fogli:2002xj}
G.~L.~Fogli, E.~Lisi, A.~Mirizzi and D.~Montanino,
%``Revisiting nonstandard interaction effects on supernova neutrino flavor  oscillations,''
arXiv:.
%%CITATION = ;%%

\bibitem{Fornengo:1999zp}
N.~Fornengo, M.~C.~Gonzalez-Garcia and J.~W.~F.~Valle,
%``On the interpretation of the atmospheric neutrino data in terms of  flavor changing neutrino interactions,''
JHEP {\bf 0007} (2000) 006  and references therein.
%%CITATION = ;%%

\bibitem{Fogli:2000ak}
G.~L.~Fogli et al,
%``Analysis of oscillations of atmospheric neutrinos,''
Nucl.\ Phys.\ Proc.\ Suppl.\ {\bf 91} (2000) 167 .
%%CITATION = ;%%
  
\bibitem{Huber:2001zw}
P.~Huber and J.~W.~F.~Valle,
%``Non-standard interactions: Atmospheric versus neutrino factory  experiments,''
Phys.\ Lett.\ B {\bf 523} (2001) 151 .
%%CITATION = ;%%

\bibitem{Freund:2001ui}
M.~Freund, P.~Huber and M.~Lindner,
%``Systematic exploration of the neutrino factory parameter space  including errors and correlations,''
Nucl.\ Phys.\ B {\bf 615} (2001) 331 %%.
%%CITATION = ;%%
K.~Dick, M.~Freund, M.~Lindner and A.~Romanino,
%``CP-violation in neutrino oscillations,''
Nucl.\ Phys.\ B {\bf 562} (1999) 29 %%.
%%CITATION = ;%%
A.~Cervera et al, %A.~Donini, M.~B.~Gavela, J.~J.~Gomez Cadenas, P.~Hernandez, O.~Mena and S.~Rigolin,
%``Golden measurements at a neutrino factory,''
Nucl.\ Phys.\ B {\bf 579} (2000) 17 [Erratum-ibid.\ B {\bf 593} (2001)
731] %%.
%%CITATION = ;%%

\bibitem{Huber:2001de}
P.~Huber, T.~Schwetz and J.~W.~F.~Valle,
%``How sensitive is a neutrino factory to the angle Theta(13)?,''
Phys.\ Rev.\ Lett.\  {\bf 88} (2002) 101804 %%.
%%CITATION = ;%%

\bibitem{Huber:2002bi}
P.~Huber, T.~Schwetz and J.~W.~F.~Valle,
%``Confusing non-standard neutrino interactions with oscillations at a  neutrino factory,''
arXiv: Phys. Rev. D
%%CITATION = ;%%

\bibitem{Aglietta:zu}
M.~Aglietta {\it et al.},
%``Upper Limit On The Solar Anti-Neutrino Flux According To Lsd Data,''
JETP Lett.\  {\bf 63} (1996) 791
%[Pisma Zh.\ Eksp.\ Teor.\ Fiz.\  {\bf 63} (1996) 753].
%%CITATION = JTPLA,63,791;%%

\bibitem{Akhmedov:2002ti}
E.~K.~Akhmedov and J.~Pulido,
%``Distinguishing magnetic moment from oscillation solutions of the solar  neutrino problem with Borexino,''
Phys.\ Lett.\ B {\bf 529} (2002) 193 %%.
%%CITATION = ;%%

\bibitem{Bahcall:1986gq} J.~N.~Bahcall, S.~T.~Petcov, S.~Toshev and
  J.~W.~F.~Valle,
%``Tests Of Neutrino Stability,''
  Phys.\ Lett.\ B {\bf 181} (1986) 369; for a recent paper see
  A.~Bandyopadhyay, S.~Choubey and S.~Goswami,
%``MSW mediated neutrino decay and the solar neutrino problem,''
Phys.\ Rev.\ D {\bf 63} (2001) 113019 and references therein.

\bibitem{atmdecay} V. Barger et al, Phys.\ Rev.\ Lett.\ {\bf 82}, 2640
  (1999) and references therein.
  
\bibitem{lsnd}
 A.~Aguilar {\it et al.}  [LSND Collaboration],
%``Evidence for neutrino oscillations from the observation of  anti-nu/e appearance in a anti-nu/mu beam,''
 Phys.\ Rev.\ D {\bf 64} (2001) 112007 .
%%CITATION = ;%%
http://www.neutrino.lanl.gov/LSND/

\bibitem{ptv92}
 J.~T.~Peltoniemi, D.~Tommasini and J.~W.~F.~Valle,
%``Reconciling dark matter and solar neutrinos,''
Phys.\ Lett.\ B {\bf 298} (1993) 383.
%%CITATION = PHLTA,B298,383;%%
A.~Ioannisian and J.~W.~F.~Valle,
%``Light Sterile Neutrino From Extra Dimensions And Four Neutrino Solutions To Neutrino Anomalies,''
Phys.\ Rev.\ D {\bf 63} (2001) 073002 .
%%CITATION = PHRVA,D63,073002;%%
M.~Hirsch and J.~W.~F.~Valle,
%``Reconciling neutrino anomalies in a simple four-neutrino scheme with  R-parity violation,''
Phys.\ Lett.\ B {\bf 495} (2000) 121 %.
%%CITATION = ;%%

\bibitem{pv93}
 J.~T.~Peltoniemi and J.~W.~F.~Valle,
%``Reconciling dark matter, solar and atmospheric neutrinos,''
 Nucl.\ Phys.\ B {\bf 406}, 409 (1993) %%.
%%CITATION = ;%%

\bibitem{Maltoni:2001bc}
M.~Maltoni, T.~Schwetz and J.~W.~F.~Valle,
%``Status of four-neutrino mass schemes: A global and unified approach to  current neutrino oscillation data,''
Phys.\ Rev.\ D {\bf 65} (2002) 093004 .
%%CITATION = ;%%

\bibitem{MiniBooNE} 
MiniBooNE's page http://www-boone.fnal.gov/

\end{thebibliography}
\end{document}






