\documentstyle[aipproc,epsfig]{article}
%\documentclass[acus]{JAC2001}
\usepackage{amsmath}
%\usepackage{color}
%\documentstyle[preprint,aps]{revtex}
%\documentstyle[epsfig]{aipproc}
%\def\topfraction{1.0}
%\def\textfraction{0.0}
%\setlength{\titleblockheight}{25mm}
%\renewcommand{\thefootnote}{\fnsymbol{footnote}}     %USE SYMBOLIC FOOTNOTE
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\normalsize\textlineskip
%\pagestyle{empty}
%\setcounter{page}{1}
%
\def\parenbar#1{{\null\!                        % left-hand spacing
   \mathop{\smash#1}\limits%                    % superscript above
   ^{\hbox{\fiverm(--)}}%                       % 5pt type
   \!\null}}                                    % right-hand spacing
%NEW ****************************************
\newcommand{\be}{\begin{equation}}
\newcommand{\ee}{\end{equation}}
\newcommand{\bea}{\begin{eqnarray}}
\newcommand{\eea}{\end{eqnarray}}
\newcommand{\noi}{\noindent}
\newcommand{\oms}{\overline{MS}}


\newcommand{\lsim}{\mathrel{\lower4pt\hbox{$\sim$}}
\hskip-11.5pt\raise1.6pt\hbox{$<$}\;}

\newcommand{\gsim}{\mathrel{\lower4pt\hbox{$\sim$}}
\hskip-11.5pt\raise1.6pt\hbox{$>$}\;}
%========================================*
%       \parenbar puts a bar in small parentheses over a character to
% indicate an optional antiparticle. \nunubar and \ppbar are special
% cases.

\def\parenbar#1{{\null\!                        % left-hand spacing
   \mathop{\smash#1}\limits%                    % superscript above
   ^{\hbox{\tiny(--)}}%                       % 5pt type
   \!\null}}                                    % right-hand spacing

\def\nunubar{\parenbar{\nu}}
\def\ppbar{\parenbar{p}}

% add words to TeX's hyphenation exception list
\hyphenation{author another created financial paper re-commend-ed Post-Script}

%NEW END*********************************************


\def\nunubar{\parenbar{\nu}}
\def\ppbar{\parenbar{p}}
\def\vep{\varepsilon}
\def\sstheff{{\sin^2\theta^{\rm eff}_W}}
\def\ssbtheff{{\sin^2\bar\theta^{\rm eff}_W}}
\def\ssthmzms{{\sin^2\theta_W(m_Z)_{\overline{MS}}}}
\def\ssthmums{{\sin^2\theta_W(\mu)_{\overline{MS}}}}
\def\ssthw{{\sin^2\theta_W}}
\def\ssthzw{{\sin^2\theta^0_W}}
\def\ssbthw{{\sin^2\bar\theta_W}}
\def\mzms{{(m_Z)_{\overline{MS}}}} 
\def\ms{{\overline{MS}}}
\def\ggam{{\gamma\gamma}}
\def\rhoppv{{\rho^\prime_{PV}}}
\def\kappappv{{\kappa^\prime_{PV}}}
\def\rqs{{R(q^2)}}
\def\numue{{\nu_\mu e}} 
\def\bnumue{{\bar \nu_\mu e}}
\def\rhonmenc{{\rho^{(\nu_\mu e)}_{NC}}} 
\def\kappanme{{\kappa^{(\nu_\mu e)}}}
\def\mtms{(m_t)_{\overline{MS}}}
\def\zp{{Z^\prime}}
\def\zpb{{\hbox{$Z^\prime$ bosons}}}
\def\pho{{\phantom{0}}}
\def\newp{``new physics''}
\def\prl{{\it Phys.\ Rev.\ Lett.\ }}
\def\pr{{\it Phys.\ Rev.\ }}
\def\np{{\it Nucl.\ Phys.\ }} 
\def\pl{{\it Phys.\ Lett.\ }}
\def\ub{\underbar}
\newcommand{\lsim}{\mathrel{\lower4pt\hbox{$\sim$}}
\hskip-12.5pt\raise1.6pt\hbox{$<$}\;}
\newcommand{\gsim}{\mathrel{\lower4pt\hbox{$\sim$}}
\hskip-12.5pt\raise1.6pt\hbox{$>$}\;}

\def\subrightarrow#1{%                          % #1 under arrow
  \setbox0=\hbox{%                              % set a box
    $\displaystyle\mathop{}%                    % no mathop
    \limits_{#1}$}%                             % just limits
  \dimen0=\wd0%                                 % get width
  \advance \dimen0 by .5em%                     % add a bit
  \mathrel{%                                    % space like =
    \mathop{\hbox to \dimen0{\rightarrowfill}}% % arrow to width
       \limits_{#1}}}                           % text below
%%
%%%     EPAC 2002 PARIS - Abstract ID WEPLE070  ,  Session ID ____________
%%%%%
%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\def\mdmatm{\Delta m^2_{32}}
\def\dmatm{$\mdmatm$}
\def\mdmsol{\Delta m^2_{21}}
\def\dmsol{$\mdmsol$}
\def\numunue{$\nu_\mu \rightarrow \nu_e$}
\def\anumunue{$\bar\nu_\mu \rightarrow \bar\nu_e$}
\def\meV{e\mbox{V}}
\def\eV{$\meV$}
\def\MeV{M\eV}
\def\GeV{G\eV}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%
%\setlength{\titleblockheight}{25mm}
%

\begin{document}
\title{Physics  of an Intense Neutrino Beam from BNL to
 a Very Long Baseline Detector} 
\author{Zohreh Parsa\thanks{Supported by US Department of Energy %\\
contract Number~DE-AC02-98CH10886. E-mail: parsa@bnl.gov, 
URL: http://www.neutrinos.bnl.gov, member of BNL Neutrino Working Group }}
\address{       Brookhaven National Laboratory,
        Physics~Dept.,~510~A,
         Upton, NY~11973,~USA. }
\maketitle

\begin{abstract}
 An intense neutrino facility allows probing of the neutrino mixing angles, 
 mass hierarchy, and leptonic CP violation. 
Physics potential, for making precision measurements of all neutrino 
oscillation parameters ($\theta_{ij}$, $\Delta m^2_{ij}$, $\delta$) 
using a wide band $\nu_\mu$ beam, to a (very long baseline) detector  
is presented. Potential of a Neutrino beam from Brookhaven National 
Laboratory to a 2540 km baseline (with 0.5 megaton) detector at 
Homestake Mine in South Dakota,  is (under study by our neutrino 
working group) discussed.   Schemaics of the beam facility for 
the  AGS upgrade to 1 MW with a cycle time of 2.5 and $10^{14}$ 
protons on target at 28 GeV;  and a map with possible detector
 sites are also included.
\end{abstract}
%\vskip-.5cm
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  section 1 Introduction %%%%%%%%%%%%%%%%%%%%%%%%%%%
\section*{Introduction}

Success of the atmospheric and solar neutrino experiments that has 
provided  evidence for non - zero neutrino masses and mixing  
 has increased our interest in neutrino oscillation searches using
 accelerator created neutrinos. Protons from an accelerator
 (e.g. AGS) would hit a target (e.g. Mercury Jet, or graphite), and 
produce bursts of particles e.g., pions, that decay to muons, which then
decay to neutrinos.
To focus the beam a magnetic horn (and/ or solenoid) can be used to keep the 
particles from spreading and to direct  the beam in the detector(s) 
direction. After leaving the horn pions decay into  neutrinos.

Upgraded conventional Neutrino horn beams (Superbeams) are 
being considered (at BNL) for probing of the neutrino masses, mixing angles,
leptonic CP violation, matter effects, new interactions, etc. We discuss in
 the following sections:  Physics \& Extra long baseline experiment; 
AGS Upgrade;  Neutrino Superbeam; Detector; and  Outlook. 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  section 2  physics %%%%%%%%%%%%%%%%%%%%%%%%%%%
%\vskip-.7cm
\section*{ physics}
  
The  Atmospheric Neutrino ``Anomaly'' suggests that
GeV $\nu_\mu$'s (from $p + N \to \pi \to \mu \nu_\mu$)
disappear while traversing the Earth's diameter, indicating
%$\Rightarrow
 $\Delta m^2_{32} = m_{3}^2 - m_{2}^2  = \pm2.0^{+1.0}_{-0.7}\times 10^{-3}
{\r(eV)}^2$ for $\sin^2 2\theta_{23} \simeq 0.85-1.0 $.  The value of  $\Delta m^2_{32}$ has decreased over the years, with recent reductions from \cite{haya}
 $3.0\to 2.5\to 2.0\times10^{-3}$eV$^2$. Fortunately this change is good,
for experiments with very long baselines (  $L\simeq2000$--$4000$)
such as our BNL to HomStake, (WIPP or Henderson) proposal }.
%(Kamiokande, IMB, Soudan-2, MACRO, Super-Kamiokande).

Solar neutrino ($\nu_e\to\nu_e{\rm~and~}\nu_e\to\nu_k$) oscillation experiments
and the Kamland reactor study of $\bar\nu_e$ disappearance prefer \cite{sno}
$\Delta m^2_{21}=m^2_2-m^2_1=7.3\pm1\times10^{-5}{\rm eV}^2 $, and
$\sin^22\theta_{12} \simeq 0.84\pm 0.10 $.

Increased interest in the Neutrino oscillation physics span from
the solar neutrino deficit and some evidence for
$\nu_{\mu} \rightarrow \nu_{e}$, 
oscillations (from the LSND experiment), as well as the exciting atmospheric
neutrino results including measurements of the atmospheric
Muon~-~Neutrino deficit from the SuperK (Superkamiokande) experiment that
has provided convincing evidence for lepton flavor violation.
The experimental results interpeted is based on oscillation of one 
neutrino flavor $\nu_e$,$\nu_\mu$ and $\nu_\tau$, (state  $|\nu_\ell>, \ell=e,\mu,\tau$) into others and are related to the neutrino mass eigenstates $|\nu_i>, i=1,2,3$ (with masses $m_i$) by $U$ a $3\times3$ unitary matrix,  with $
c_{ij} = \cos\theta_{ij}$, and $ s_{ij}=\sin\theta_{ij}$:
\vskip.3cm

%
%\be
%\pmatrix{|\nu_e> \cr |\nu_\mu> \cr |\nu_\tau>} = U \pmatrix{|\nu_1> \cr
%|\nu_2> \cr |\nu_3>}
%\ee

\[
U = \pmatrix{ c_{12}c_{13} & s_{12}c_{13} & s_{13}e^{-i\delta} \cr
-s_{12}c_{23}-c_{12}s_{23}s_{13}e^{i\delta} &
c_{12}c_{23}-s_{12}s_{23}s_{13}e^{i\delta} & s_{23}c_{13} \cr
s_{12}s_{23}-c_{12}c_{23}s_{13}e^{i\delta} &
-c_{12}s_{23}-s_{12}c_{23}s_{13}e^{i\delta} & c_{23}c_{13}} \nonumber
\label{eqone}
\]
\vskip.3cm
%\[
%c_{ij} = \cos\theta_{ij} \quad , \quad s_{ij}=\sin\theta_{ij}
%\]}

%%%%%%%%%%%%%%%%%%%  subsection  Extra Long Baseline %%%%%%%%%%%%%%%%%%%%%%%
%\vskip-.5cm
\subsection*{Extra Long-Baseline Physics} 
Extra-long neutrino flight paths provide the possibility of observing
multiple nodes of the neutrino oscillation (probability) in appearance and 
 disappearance experiments.  Observation of such a pattern will  directly 
demonstrate the oscillatory nature of the flavor changing phenomenon.
For fixed distance L, the oscillation maxima will occur roughly at energies 
of
\begin{eqnarray}
 E_{\nu}(n) = \frac{\Delta{m_{32}^2}L}{2(2n-1)\pi},\nonumber\\
 n= 1,2,3, \ldots\nonumber
\end{eqnarray}

For a given $E_{\nu}$ and $L$,  the oscillation of \numunue{} appearance can be
 described by:% the following equation:  
%\begin{center}
%\begin{alignat}%{2}
%  \label{eq:one}
%\lefteqn{P(\nu_\mu\to\nu_e)  =  4(s^2_2s^2_3c^2_3 +J_{CP}\sin\Delta_{21})
%\sin^2\frac{\Delta_{21}}{2}}\nonumber\\ %\hspace{1em} \nonumber\\
%& +2(s_1s_2s_3c_1c_2c^2_3 \cos\delta -s^2_1s^2_2s^2_3c^2_3) \sin
%\Delta_{31} \sin \Delta_{21}\nonumber \\
%& +4(s^2_1c^2_1c^2_2c^2_3 +s^4_1s^2_2s^2_3c^2_3 
%-2s^3_1s_2s_3c_1c_2c^2_3 \cos\delta \nonumber \\
%& -J_{CP} \sin\Delta_{31})\sin^2\frac{\Delta_{21}}{2} 
% +8(s_1s_2s_3c_1c_2c^2_3 \cos\delta \nonumber \\
%&- s^2_1s^2_2s^2_3c^2_3) \sin^2
%\frac{\Delta_{31}}{2} \sin^2 \frac{\Delta_{21}}{2} + {\rm matter\ effects} \nonumber
%\end{alignat}
%\end{center}
\begin{eqnarray}
  \label{eq:one}
P(\nu_\mu\to\nu_e)  &=&4(s^2_2s^2_3c^2_3 +J_{CP}\sin\Delta_{21})
\sin^2\frac{\Delta_{21}}{2}\nonumber\\ %\hspace{1em} \nonumber\\
&& +2(s_1s_2s_3c_1c_2c^2_3 \cos\delta -s^2_1s^2_2s^2_3c^2_3) \sin
\Delta_{31} \sin \Delta_{21}\nonumber \\
&& +4(s^2_1c^2_1c^2_2c^2_3 +s^4_1s^2_2s^2_3c^2_3 
-2s^3_1s_2s_3c_1c_2c^2_3 \cos\delta \nonumber \\
&& -J_{CP} \sin\Delta_{31})\sin^2\frac{\Delta_{21}}{2} 
 +8(s_1s_2s_3c_1c_2c^2_3 \cos\delta \nonumber \\
&&- s^2_1s^2_2s^2_3c^2_3) \sin^2
\frac{\Delta_{31}}{2} \sin^2 \frac{\Delta_{21}}{2} + {\rm matter\ effects} \nonumber
\end{eqnarray}
Where, $c_{i}\equiv cos\theta_{i}$, $s_{i}\equiv sin\theta_{i}$,   
$J_{CP} \equiv s_1s_2s_3c_1c_2c^2_3\sin\delta$,  
$\Delta_{31}\equiv \Delta m^2_{31}L/2E_{\nu}$, and
$\Delta_{21}\equiv \Delta m^2_{21}L/2E_{\nu}$.

$J_{CP}$ is an invariant that quantifies CP violation in the neutrino
sector; $\Delta_{31}$ is the atmospheric term and $\Delta_{21}$ is the solar term \cite{marc}. 
For $P(\bar\nu_\mu\to\bar\nu_e)$ the above formula holds 
except  $J_{CP}$ terms will have opposite sign and matter effect will change.

The oscillation is primarily due to the first term linear 
in $\sin^2\frac{\Delta_{31}}{2}$, and oscillation probability rises for 
lower energies due to the terms linear in $\sin^2 \frac{\Delta_{21}}{2}$.  

The interference terms involve CP violation and they create an asymmetry 
between neutrinos and anti-neutrinos.  The CP asymmetry grows 
linearly with distance:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   Figure 1
%%%%%%
\begin{figure}%[bh]
\centerline{\epsfig{file=WEPLE070fig1.epsi,width=4.0in }}
\caption{Probability of \numunue{} 
      and \anumunue{}  oscillations at 2540 km
      assuming
 a $45^o$ CP violation phase, including matter effect.
%Here $\theta_{ij}$ notation is used \cite{LOI} rather
% than $\theta_{i}$.
 }
 \label{cp45asym}
\end{figure}
%%%%%%%5

%\bea
%A_{\rm CP}  =  \frac{P(\nu_{\mu} \rightarrow\ \nu_{e})- P(\bar\nu_\mu\to\bar\nu_e)}{P(\nu_{\mu} \rightarrow \nu_{e})+ P(\bar\nu_\mu\to\bar\nu_e)}\\
%  & \simeq & \frac{2s_{1}c_{1}c_{2}sin\delta}{s_{2}s_{3}}
%\frac{\Delta m_{21}^2}{\Delta m_{31}^2}
%\frac{\Delta{m_{31}^2}L}{4E_{\nu}} + O(\Delta_{21}^2) \\
%& & + {\rm matter~effects}. \nonumber\\ }
%\eea
\bea
A_{\rm CP}  &=&  \frac{P(\nu_{\mu} \rightarrow\ \nu_{e})- P(\bar\nu_\mu\to\bar\nu_e)}{P(\nu_{\mu} \rightarrow \nu_{e})+ P(\bar\nu_\mu\to\bar\nu_e)}\nonumber\\
  & \simeq & \frac{2s_{1}c_{1}c_{2}sin\delta}{s_{2}s_{3}}
\frac{\Delta m_{21}^2}{\Delta m_{31}^2}
\frac{\Delta{m_{31}^2}L}{4E_{\nu}} + O(\Delta_{21}^2) \nonumber\\
& & + {\rm matter~effects}. \nonumber\\ }
\eea

\noi or is given by (to  order of $\Delta m^2_{21}$ assuming
$\sin^22\theta_{13}$ is not too small)  

\bea
A_{CP} & \simeq & \frac{\cos\theta_{23} \sin2\theta_{12}
\sin\delta}{\sin\theta_{23} \sin\theta_{13}} \left(\frac{\Delta
m^2_{21}L}{4E_\nu} \right) \nonumber \\
& & + {\rm matter~effects}\nonumber
\eea

\noi In this expression, the asymmetry grows linearly with distance 
and increases as $\theta_{13}$ gets smaller, noted earlier. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Figure 2
\begin{figure}[bh]
\centerline{\epsfig{file=dm20025_sp.ps,width=4.0in }}
\caption{Example of expectected $\nu_\mu$ disappearance spectra,  
without oscillations; (the middle)  with oscillation and (the
  bottom histogram shows) the background contribution to the
  oscillating spectrum. This spectrum has improved with decrease in
 $\Delta m^2_{32}$, with recent reductions from \cite{haya}
 $3.0\to 2.5\to 2.0\times10^{-3}$eV$^2$.  }
\label{nwg-0303081}
\end{figure}
 
 Fig.~\ref{cp45asym} includes the matter effect since matter
 will enhance (suppress) neutrino (anti-neutrino) conversion at
high energies and will  lower (increase) the energy at which the
oscillation maximum occurs, detection of matter enhancement effect  
can be made by measuring the asymmetry between neutrino and anti-neutrino
oscillations (or by measuring the spectrum of electron neutrinos which
 also provide the sign of \dmatm{}).

If both the CP violation and the signal to \numunue{} is  large then effects
 of CP violation can be measured with only the ($\nu_\mu$) 
neutrino beam. It grows linearly with decrease in 
energy or the increase in baseline. For extra-long baseline experiments, 
comparison of the signal strength in the $\pi/2$ node versus 
the $3\pi/2$ (or higher) nodes will provide measurements of CP violation.


%%%%%%%%%%%%%%%%%%%%%%%%%%***  SECTION 3 ******************************
\vskip-.5cm
\section*{ AGS Upgrade and {$\nu$}beam}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Figure 3
\begin{figure}[bh]
\centerline{\epsfig{file=nu-site-goodink.ps,width=4.in }}
\caption{Schematic of the BNL-AGS RHIC facility showing location of
the new beam-line for sending a neutrino beam to Homestake mine in 
South Dakota, and any detector in the Western direction.}
%\vskip-.35cm
\label{AGS-RHIC}
\end{figure}
%%%%%%%%%
The preliminary design of the BNL-AGS upgrades and the new neutrino beam
has been produced by the AGS department \cite{beav} to reach an 
AGS power of e.g. 0.53 MW ($1.2\times{10^{21}}ppp$) in its first phase and
 1.3 MW ($1.2\times {10^{21}}ppp$) in the second phase.% \cite{roser01}. 
 In the first  phase the LINAC will be improved to
inject protons to the booster at 400 \MeV{} (at present it is 200
\MeV{}), and the booster energy increases to 2.5 \GeV{} from 1.8
\GeV{}. The addition of a fixed field accumulator storage ring between
the booster and the AGS main ring will increase the AGS input beam
from the present 4 booster pulses per AGS acceleration to 6 booster
pulses per AGS acceleration and, at the same time, increase the AGS
frequency from 0.6 Hz to 1.0 Hz. The AGS power increase would be from
0.14 to 0.53 MW.    
The new accumulator will be in the same tunnel as the AGS.   
In the second phase of the upgrades the AGS repetition rate will 
be increased to 2.5 Hz to reach a total beam power of 1.3 MW. 

%%%A preliminary design for a beam to Homestake is 
%%%shown in \ref{fig5}{} %and \ref{eleview}.
%%%Due to the constraint at BNL to keep the beam line above the water 
%%%table (which is at a shallow depth of $\sim$ 20 m) on Long Island, the 
%%%beam line is to be constructed on a hill with 11.5 degree slope. 
%And is anticipated to be relatively easy, and
% inexpensive to build such hill  given the flat, and  sandy 
%geology of Long Island.  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Figure 4
%%%\begin{figure}[tbh]
%%%\centerline{\epsfig{file= 3d_view.ps,width=4.in }}
%%%\caption{3 Dimensional Perspective of our BNL proposed Neutrino 
%%%Superbeam}
%%%\label{3d_view}}
%%%\end{figure}
%%%%%%%
%%%%%%%%%%%%%%%%%%%%%  
The proton beam is to be elevated to a target station on top of the hill.
 And the new proposed  fast extracted proton beam line in the U-line tunnel
will 
come off the line feeding RHIC. And will turn  west, a few
hundred meters before the horn-target building. In addition to its 90
degree bend, the extracted proton beam will be bent upward through
13.76 degrees to strike the proton target.  The downward 11.30 degree
angle of the 667.8 ft meson decay region will then be aimed at the
2500 meter level of the Homestake Laboratory. This will require
the construction of a 39 meter hill to support the target-horn building,
so as to  avoid any penetration of the water table.  At its midpoint
(about Lake Michigan) the center of the neutrino beam will be roughly
120 km below the Earth's surface.
(For a shorter baseline e.g., to Lansing NY in approximately the same 
direction as Homestake the hill won't be needed. 
Various combinations of the proton transport and the target station 
for the  extra-long, (short/intermediate) baselines are being considered.) 
%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Figure 3
\begin{figure}[tbh]
\centerline{\epsfig{file=map3names.eps,width=4.0in }}
\caption{Possible extra long neutrino baselines from BNL to  Lead (Homestake) SD ($\sim 2540 Km$, 11.5~degrees dip angle),  to Carlsbad (WIPP) NM ($\sim 2900 km$,  13.0~degrees), and to the Henderson Mine in Colorado.} 
%dip angle).}
\label{map3names}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%  SECTION    %%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\vskip-.7cm

\section*{Detectors for the very long baseline experiment} 

There is an  interest to convert the Homestake Gold Mine in Lead, 
 South Dakota into a National Underground Science Laboratory (NUSL).
This  will provide unique opportunity for an extra-long baseline 
neutrino oscillation experiments from BNL.   The extra-long  
baseline is  2540 km from the (Brookhaven  National Laboratory) BNL
to Lead, South Dakota.  
The proposed NUSL facility is to  accommodate an array of detectors with 
about  1 Megaton total mass. Most of these will be water Cerenkov 
detectors that can observe neutrino interactions in the desired energy 
range with sufficient energy and time resolution.

Other detector types (e.g. Liquid Argon),
and sites are also being considered,~e.g.,Henderson Mine in Colorado,
~the Waste Isolation Pilot Plant~(WIPP)~located in an ancient salt 
bed at a depth of $\sim 700 m$ near Carlsbad, New Mexico, etc. The 
distance from BNL to WIPP is about 2880 km,. The cosmic ray 
background will be higher at WIPP because the facility is not as deep 
as Homestake (with levels as deep as $\sim 2500 m$). 
\vskip-.8cm
%%%%%%%%%%%%%%%%%%%%%%%%********* SECTION   *********
\vskip-.4cm
\section*{Outlook}
 Four goals of neutrino physics: precise determination of 
\dmatm{}, observation of \numunue{} appearance,
measurement of matter effects, and detection of CP violation are all
 possible with an intense neutrino broad band beam, 
very long distance baseline, and large detector. Both very long O(2500 km)
 and intermediate O(400 km) baseline experiments  can be staged (from
 Brookhaven) as the AGS is upgraded to .5 MW, as much as 2.5 MW or 
higher (4 MW needed for a Neutrino Factory). AGS improvements will  
also allow rare muon and kaon decay  studies, muon EDM measurements, etc. 
Thus providing  additional  windows for discovery.

\vskip-.9cm
\begin{thebibliography}{99}

\bibitem {zp} Z. Parsa, $\nu$~Superbeam for
 BNL-sec Procds.~PAC20001, Chicago, Ill.
; ibid, Procds.  EPAC2002, Paris France (2002). Presentation here are, from [1]-[8] and Refs therein.

\bibitem{bnl} BNL Neutrino Working Group: M. Diwan {\it et al}.,
hep-ex/0211001.%\bibitem{diwan} %M. Diwan {\it et al}.
ibid, Phys.\ Rev.\ {\bf D68}, 012002;  LOI submitted to BNL (April
 2002) and Refs therein.

\bibitem{haya} Y. Hayato, HEP2003 talk, Aachen Germany July 2003. For a
review see V. Barger, D. Marfatia and K. Whisnant, Int.\ J. Mod.\
Phys.\ {\bf E12}, 569 (2003).  %1

\bibitem{sno} SNO Collaboration, Q. Ahmad {\it et al}.,
nucl-ex/0309004; A. Balantekin and H. Yuksel, hep-ph/0309079.  %2

\bibitem{choo} CHOOZ: M. Apollonio {\it et al}., Phys.\ Lett.\
{\bf B466}, 415 (1999), Palo Verde: F. Boehm {\it et al}., Nucl.\
Phys.\ Proc.\ Suppl.\ {\bf91}, 91 (2001)  %3

\bibitem{jarl} C. Jarlskog, Z. Phys.\ {\bf C29}, 491 (1985); Phys.\
Rev.\ {\bf D35}, 1685 (1987).

\bibitem{marc} W. Marciano, ``Extra Long Baseline Neutrino Oscillations
and CP Violation'', hep-ph/0108181, (22~Aug, 2001). %5

\bibitem{beav} D Beavis {\it et al}., BNL Report 52459, AGS proposal 889.  %6

\end{thebibliography}
\enddocument
\end{document}

%\bibitem{3m} \cite{LOI} 

%\bibitem{AGS} \cite{LOI}

%\bibitem{homestk}
%National Underground Science Laboratory at Homestake, Lead, SD,  \\
%http://mocha.phys.washington.edu/NUSL/

%\bibitem{wipp}
%Waste Isolation Pilot Plant, Carlsbad, NM, \\
%http://www.wipp.carlsbad.nm,us/
%
%\bibitem{sk}
%S. Fukuda et al., Phys. Rev. Lett. {\bf 86} 5656, 2001;  
%E.W. Beier, Phys. Lett. {\bf B283}, 446 (1992); 
%T. Kajita and Y. Totsuka, Rev. Mod. Phys. {\bf 73}, 85 (2001).
%
%\bibitem{imb}
%C. McGrew {\it et al.}, Phys.\ Rev.\ D {\bf 59}, 052004 (1999).
%
%\bibitem{sno} 
%Q. R. Ahmad et al., Phys. Rev. Lett. {\bf 87} 071301 (2001). 
%S. Fukuda et al., Phys. Rev. Lett., {\bf 86} 5651 (2001). 
%
%\bibitem{lsnd}
%C. Athanassopoulos et al., Phys. Rev. Lett. {\bf 77} 3082 (1996); 
% C. Athanassopoulos et al., Phys. Rev. Lett. {\bf 81} 1774 (1998) 
%
%\bibitem{boon}
%Booster Neutrino Experiment, Fermi National Laboratory, \\
%http:/www-boone.fnal.gov/
%
%
%\bibitem{k2k}
%S. H. Ahn et al., Phys. Lett. {\bf B 511} 178 (2001).
%
%\bibitem{j2k}
%The JHF-Kamioka neutrino project, Y. Itow et al., 
%arXiv:hep-ex/0106019, June 2001.
%
%\bibitem{kasuga} S. Kasuga {\it et al.}, 
%Phys Lett. {\bf B374}, 238 (1996).
%
%\bibitem{minos}
%Numi MINOS project at Fermi National Accelerator Laboratory, \\
%http:/www-numi.fnal.gov/ 
%
%\bibitem{cngs} 
%CERN Neutrinos to Gran Sasso, \\
%http://proj-cngs.web.cern.ch/proj-cngs/
%
%\bibitem{e889}
%p889 Collaboration, Physics Design Report, BNL No. 52459, April, 1995. \\
%http://minos.phy.bnl.gov/nwg/papers/E889/

%\bibitem{marciano} W. Marciano hep-phy/0108181 (22~Aug, 2001).

%\bibitem{e734}
%L. A. Ahrens et al., Phys. Rev. {\bf D 34}, 75 (1986).
%
%\bibitem{e734d}
%L. A. Ahrens et al., Phys. Rev. {\bf D 41}, 3297 (1990).
%
%\bibitem{arafune} Jiro Arafune, Masafumi Koike 
%and Joe Sato, Phys. Rev {\bf D56}, 3093 (1997).
%
%\bibitem{marciano} William J. Marciano, 
%arXiv: hep-phy/0108181,  22 Aug 2001.
%
%%\bibitem{irina} Irina Mociouiu and
% Robert Schrock, 
%arXiv: hep-ph/0106139v3, 15 Nov. 2001
%
%\bibitem{wolfenstein}
%By L. Wolfenstein (Carnegie Mellon U.). 1978.
%In *West Lafayette 1978, 
%Proceedings, Neutrinos '78*, West Lafayette 1978, 
%C3-C6 and *Washington 1978, 
%Proceedings, Long-distance Neutrino Detection*, 108-112.
%
%%\bibitem{study2}
% S.~Ozaki et al., eds.,
% {\sl Feasibility Study-II of a Muon-Based Neutrino Source}
% (June 14, 2001), %\hfill\break
% http://www.cap.bnl.gov/mumu/studyii/FS2-report.html
%
%

%\bibitem{3m}3M Collaboration: Proposal titled: Megaton 
%Modular Multi-Purpose Neutrino Detector, Nov. 26, 2001.
%
%\bibitem{uno} Physics Potential and Feasibility of UNO,
%UNO collaboration, June 2001. 
%
%\bibitem{lannddp}
% D.B.~Cline, F.~Sergiampietri, J.G.~Learned, K.T.~McDonald,
% {\sl LANNDD, A Massive Liquid Argon Detector for Proton Decay,
%Supernova and
% Solar Neutrino Studies, and a Neutrino Factory Detector}
% (May 24, 2001), astro-ph/0105442 \\
% \bibitem{franco1}
%Also see F.~Sergiampietri,
% {\sl On the Possibility to Extrapolate Liquid Argon Technology to a
%Supermassive
% Detector for a Future Neutrino Factory},
% presented at NuFACT'01 (May 26, 2001),
%
%\bibitem{larhighres}
%F. Arneodo, et al., 
%Nucl. Instrum. Meth. {\bf A461} 324 (2001)  
%
%
%\bibitem{argonprop}
% http://www.hep.princeton.edu/\~mcdonald/nufact/bnl\_loi/argonprop.pdf
% {\sl 1 MW AGS proton driver},
% presented by T.~Roser at Snowmass'01 (June 2001), \hfill\break

%\bibitem{foster} We have consulted
%Bill Foster at Fermilab to design and calculate the cost of the
%accumulator ring proposed here.



%\bibitem{edm} \cite{edm-rsvp}
%\bibitem{edm} %Muon Electric Dipole Moment experiment. \\
%http://www.bnl.gov/edm/

%\bibitem{rsvp} \cite{edm-rsvp}
%\bibitem{rsvp}
%Rare Symmetry Violating Processes, \\
%http://meco.ps.uci.edu/RSVP.html

\end{thebibliography}
\enddocument
\end{document}

