%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%    CERN to Gran Sasso; An ideal distance for superbeam?
%    by Hisakazu Minakata and Hiroshi Nunokawa
%    Talk presented by H. Minakata at TAUP2001 workshop, 
%    Laboratori Nazionali del Gran Sasso, Italy, September 8-12, 2001.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\documentstyle[twoside,fleqn,espcrc2,psfig]{article}
\newcommand{\lsim}{\raisebox{0.3mm}{\em $\, <$} 
\hspace{-3.3mm} \raisebox{-1.8mm}{\em $\sim \,$}}
\def\VEV#1{\left\langle #1\right\rangle}
\def\lsim{\raise0.3ex\hbox{$\;<$\kern-0.75em\raise-1.1ex
\hbox{$\sim\;$}}}
\def\gsim{\raise0.3ex\hbox{$\;>$\kern-0.75em\raise-1.1ex
\hbox{$\sim\;$}}}
\renewcommand{\topfraction}{1}
\renewcommand{\bottomfraction}{1}
\renewcommand{\textfraction}{0}

\topmargin = -1.0cm

\def\npb#1{Nucl.\ Phys.\ {\bf B\,#1}}
\def\npbps#1{Nucl.\ Phys.\ (Proc. Suppl.) {\bf B\,#1}}
\def\plb#1{Phys.\ Lett.\ {\bf B\,#1}}
\def\prc#1{Phys.\ Rev.\ {\bf C\,#1}}
\def\prd#1{Phys.\ Rev.\ {\bf D\,#1}}
\def\prl#1{Phys.\ Rev.\ Lett. {\bf#1}}
\def\ncim#1{Nuo.\ Cim.\ {\bf #1}}
\def\zpc#1{Z.~Phys.\ {\bf C\,#1}}
\def\epjc#1{Eur.~Phys.~J.\ {\bf C\,#1}}
\def\prep#1{Phys.\ Rep.\ {\bf #1}}
\def\rmp#1{Rev.\ Mod.\ Phys.\ {\bf #1}}
\def\ijmpa#1{Int.\ J.\ Mod.\ Phys.\ {\bf A\,#1}}
\def\progtp#1{Prog.\ Th.\ Phys.\ {\bf #1}}
\def\mpla#1{Mod. Phys. Lett. {\bf A\,#1}}
\def\sjnp#1{Sov. J. Nucl. Phys. {\bf #1}}
\def\JHEP#1{JHEP {\bf #1}}

\textheight = 22.5cm

\begin{document}
%\pagestyle{empty}

\title{CERN to Gran Sasso; An ideal distance for superbeam?
\thanks{Talk presented at 7th International Workshop on Topics in 
Astroparticle and Underground Physics (TAUP2001), 
Laboratori Nazionali del Gran Sasso, Italy, September 8-12, 2001.
}}

\author{Hisakazu Minakata$^a$ and Hiroshi Nunokawa$^{b,c}$ \\
        {\ }\\
$^a$
Department of Physics, Tokyo Metropolitan University, 
1-1 Minami-Osawa, Hachioji \\
Tokyo 192-0397, Japan \\ 
\vglue 0.15cm
$^b$
Instituto de F\'{\i}sica Te\'orica, Universidade Estadual Paulista, 
Rua Pamplona 145 \\ 
01405-900 S\~ao Paulo, SP Brazil \\
\vglue 0.15cm
$^c$
Instituto de F\'{\i}sica Gleb Wataghin, 
Universidade Estadual de Campinas -- UNICAMP\\
P.O. Box 6165, 13083-970 Campinas SP Brazil
} 

\begin{abstract}
We use the CP trajectory diagram as a tool for pictorial 
representation of the genuine CP and the matter effects to 
explore the possibility of an {\it in situ} simultaneous 
measurement of $\delta$ and the sign of $\Delta m^2_{13}$. 
We end up with a low-energy conventional superbeam experiment 
with a megaton-class water Cherenkov detector and baseline length 
of about 700 km. 
%
A picturesque description of the combined ambiguity which may 
arise in simultaneous determination of $\theta_{13}$ and 
the above two quantities is given in terms of CP trajectory diagram. 

\end{abstract}

\maketitle

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Section I %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\newpage
%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction}

Exploring the structure of lepton flavor mixing and the neutrino 
mass pattern is one of the most challenging goals of contemporary 
particle physics. Among other things, the least known is the 
(1-3) sector of the MNS matrix, the angle $\theta_{13}$, the sign of 
$\Delta m^2_{13}$, and the CP violating angle $\delta$. 
(See e.g., ref. \cite {NOW2000mina}.)

We have explored in a series of papers \cite {MN} the features of 
interplay between the CP phase and matter effect, and tried to 
develop strategies for measurement of lepton CP violation in 
long-baseline neutrino oscillation experiments. 
See also ref. \cite {cp-matter} for related works.
Along the line of thought, we have introduced recently a 
powerful tool called ``CP trajectory diagram in bi-probability space'' 
which allows us a separate pictorial representation of the genuine 
CP and the matter effects \cite {MNjhep01}. 
%
We pointed out that an ambiguity exists in determination 
of these parameters in a correlated way 
($\delta$ $-$ sign of $\Delta m^2_{13}$). 

We address here the issue of an {\it in situ} simultaneous 
determination of $\delta$ and the sign of $\Delta m^2_{13}$ 
in a single experiment. 
%
In a companion article \cite {nufact01}, which is a 
contribution to Proceedings of NuFACT01, a concise summary 
of the idea of CP trajectory diagram and its use is given and the 
principle of optimizing beam energy is discussed. 

Let us start with a brief remark on the more generic feature 
of the ambiguity problem. 

\section{Clover-leaf ambiguity}

If the value of $\theta_{13}$ is unknown, there exists another 
ambiguity in a correlation ($\delta - \theta_{13}$), 
as pointed out in ref. \cite {B-Castell}.
Together with the ambiguity we have uncovered, there exists the  
combined ambiguity which can be as large as four-fold. 
We now demonstrate that it can be described 
in a simple picturesque way by using the CP trajectory diagram.


In fig. \ref{Fig1} drawn is the CP trajectory diagram for 
four values of the mixing parameters which are given in the 
caption of fig. \ref{Fig1} assuming, for simplicity, 
neutrino beam with the Gaussian type energy distribution
as used in our recent work \cite {MNjhep01}. 
%
The point of fig. \ref{Fig1} is that the four-fold solutions 
are possible for given oscillation probabilities of 
$P(\nu)\equiv P(\nu_{\mu} \rightarrow \nu_{e})$ and 
$P(\bar{\nu}) \equiv P(\bar{\nu}_{\mu} \rightarrow \bar{\nu}_{e})$, 
both at about 1.1 \% in fig. \ref{Fig1}.
One notices from fig. \ref{Fig1} that the name 
"clover-leaf ambiguity" by which the ambiguity is referred at 
TAUP2001 is quite natural and appealing.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[ht]
\vglue -0.2cm 
\centerline{\protect\hbox{
\psfig{file=taup01_fig1.eps,height=6.2cm,width=6.8cm}
}}
\vglue -1.4cm 
\caption{\small
Illustration of the clover-leaf ambiguity in termes of CP 
trajectory diagram; four solutions exist for given values of 
$P(\nu_\mu\to\nu_e)$ and $P(\bar{\nu}_\mu\to\bar{\nu}_e)$ 
at neutrino energy $\langle E \rangle=1.5$ 
GeV and baseline distance $L=295$ km. 
Corresponding values of $\sin^22\theta_{13}$ as well as 
the sign of $\Delta m^2_{13}$ whose absolute value is 
$3 \times 10^{-3}$ eV$^2$ are indicated in the plot.  
%
The remaining mixing parameters are taken as;
$\sin^22\theta_{23} = 1.0$, 
$\Delta m^2_{12} = 5 \times 10^{-5}$ eV$^2$, 
$\sin^22\theta_{12} = 0.8$, and 
$\delta = \pi/2$. 
We take the matter density as $\rho = 2.72$ g/cm$^3$ and 
the electron fraction as $Y_e$ = 0.5. 
}
\label{Fig1}
\vglue -0.5cm 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section {Optimal distance for measuring the sign of $\Delta m^2_{13}$}

Now we turn to the discussion of our original question, i.e., 
how to resolve the ($\delta$ $-$ sign of $\Delta m^2_{13}$) ambiguity. 
Since the sign of $\Delta m^2_{13}$ can be determined by measuring 
interference between the vacuum and the matter effects, it is 
natural to think about neutrino oscillation experiments which 
utilize longer baselines. One can think of these possibilities in 
the context of either  

\noindent
(i) a single detector experiment for {\it in situ} simultaneous 
determination of $\delta$ and the sign of $\Delta m^2_{13}$, or

\noindent
(ii) a two-detector experiment with second supplemental detector 
which is primarily devoted for determination of the sign of 
$\Delta m^2_{13}$.

The next question to ask is; what is the optimal baseline length 
for this purpose? It would be the best situation if we can tune the 
distance in such a way that the matter effect is relatively enhanced 
compared to the genuine CP violating effect. 
To quantify the request we define the asymmetry parameter defined 
by using the ratio 
$R(P) \equiv \langle P(\nu_\mu \to \nu_e) \rangle /
\langle P(\bar{\nu}_\mu \to \bar{\nu}_e) \rangle$ as 
\begin{equation}
A(R) \equiv 
\frac{R(P;\Delta m_{13}^2>0)-R(P;\Delta m_{13}^2<0)}
{R(P;\Delta m_{13}^2>0)+R(P;\Delta m_{13}^2<0)},
\label{asymmetry}
\end{equation}
where the probabilities are averaged over Gaussian
type energy distributions. 
By using the ratio in defining the asymmetry it is 
insensitive to the values of the mixing parameters, 
in particular to $\delta$.

In fig. \ref{Fig2}, the asymmetry $A(R)$ is plotted for 
$\delta = \pi/2$ with the same mixing parameters as in 
fig. \ref{Fig1}, apart from fixing $\sin^22\theta_{13}$ 
to be 0.05. One notices that the asymmetry 
is large at $L = 600-700$ km, and at $1000-1500$ km 
for $E \sim 1$ GeV. It can be explicitly shown that 
the asymmetry is indeed very insensitive to $\delta$ 
\cite {MNjhep01}. Therefore, the baselines $L \sim 700$ km, 
and at $\sim 1000$ km look promising. We take the former 
option and examine the feature of CP-matter interplay by 
using the CP trajectory diagram. 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[ht]
\vglue -0.9cm 
\centerline{\protect\hbox{\hglue -1.0cm 
\psfig{file=taup01_fig2.eps,height=6.7cm,width=7.7cm}
}}
\vglue -1.3cm 
\caption{\small
Asymmetry of the probability ratio defined in Eq. (\ref{asymmetry})
in the text computed for $\delta = \pi/2$.
The mixing parameters are chosen as the same with those of 
fig. \ref{Fig1} and $\sin^22\theta_{13} = 0.05$.
}
\label{Fig2}
\vglue -1.0cm 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\vglue -0.5cm 
\section {Longer baseline option; single vs. two-detector methods}

In fig. \ref{Fig3}, we present (a) CP trajectory diagram 
in bi-probability plane for neutrino beam with the Gaussian type 
energy distribution with averaged energy $\langle E \rangle =1.5$ GeV 
and baseline distance $L=700$ km, and 
(b) CP trajectory diagram on number of events plane with the 
same baseline for appearance channels 
$\nu_{\mu} \rightarrow \nu_{e}$ and 
$\bar{\nu}_{\mu} \rightarrow \bar{\nu}_{e}$.  

One can clearly see in fig. \ref{Fig3}a that the two 
trajectories corresponding to 
$\Delta m^2_{13} >0$ (solid line) and 
$\Delta m^2_{13} <0$ (dashed line) 
are well separated with each other. Therefore, it is 
in principle possible to carry out an {\it in situ}
simultaneous measurement of $\delta$ and the sign of 
$\Delta m^2_{13}$ in such a baseline length. 
It is amusing to note that the length just corresponds to 
either CERN $\rightarrow$ Gran Sasso, or 
Fermilab $\rightarrow$ Soudan mine distances. 


In fig. \ref{Fig3}b, we assume a water Cherenkov detector of 
fiducial volume 0.9 Mton, and 4 MW of proton beam power 
which is planned in the JHF experiment in its phase II \cite {JHF}.
%
We use the narrow band (NB) 3 GeV beam whose neutrino energy peaks 
at $E \sim 1.4$ GeV considered in ref. \cite {JHF}, 
but with intensity multiplied by factor of 3. It is to 
mimic the off-axis (OA) beam which is designed for  $E \sim 1.4$ GeV, 
a bit of higher energy than the one actually prepared for the 
JHF experiment\cite {JHF}. 
Two (six) years of running is assumed for neutrino 
(antineutrino) channel.
We refer ref. \cite {MNjhep01} for a detailed explanation of how the 
computation of number of events is done.  

As you see in fig. \ref{Fig3}b, the numbers of events are sizable, 
some 1000 $-$ 2000, though not gigantic. Resultant 3 $\sigma$ 
contour is small enough so that simultaneous measurement of 
$\delta$ and the sign of $\Delta m^2_{13}$ seems feasible, 
justifying the title of this talk.  

If the detector at baseline $L=700$ km cannot be so massive 
by some reasons, or if the practical site requires much longer 
distance, then it can be viewed as a secondary (farthest) detector, 
assuming that the primary far detector (e.g., Hyper-Kamiokande 
\cite {totsuka} in the case of the JHF experiment) already exists. 
In this case, the secondary detector is primarily for measurement 
of the sign of $\Delta m^2_{13}$ and the requirement of statistics 
can be relaxed.

Finally, it was pointed out in the talk that while the present tunnel 
in the Gran Sasso Laboratory is too small to accommodate a megaton water 
Cherenkov detector, a 180 kton "proto-type" is ready to be 
created upon filling water into the whole tunnel!


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[ht]
\vglue -0.5cm 
\centerline{\protect\hbox{
\psfig{file=taup01_fig3a.eps,height=7.0cm,width=7.6cm}
}}
\vglue -1.0cm 
\centerline{\protect\hbox{
\psfig{file=taup01_fig3b.eps,height=7.0cm,width=7.6cm}
}}
\vglue -1.5cm 
\caption{\small
CP trajectory diagrams in 
(a) bi-probability space, and in 
(b) event number space $N(e^-) - N(e^+)$ 
for baseline distance of $L=700$ km.
%
In (b) a factor of 3 intensified NB 3 GeV beam is used to imitate
OA beam of peak energy 1.4 GeV, and the dotted circles 
correspond to 3 $\sigma$ statistical uncertainty.}
\label{Fig3}
\vglue -0.6cm 
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%




%%%%%%%%%%%%%%%%
\section*{Acknowledgments}
%%%%%%%%%%%%%%%%

This work was supported by the Brazilian funding agency
%Funda\c{c}\~ao de Amparo \`a Pesquisa do Estado de S\~ao Paulo (FAPESP),
FAPESP, 
and by the Grant-in-Aid for Scientific Research in Priority Areas
No. 12047222, Japan Ministry of Education, Culture, Sports, Science
and Technology.



%%%%%%%%%%%%%%%%%%%%%%%%%% Bibliography %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{thebibliography}{99}

\bibitem {NOW2000mina}
H.~Minakata,
Nucl.\ Phys.\ Proc.\ Suppl.\  {\bf 100} (2001) 237.


\bibitem {MN}
H. Minakata and H. Nunokawa, Phys. Rev. {\bf D57} (1998) 4403; 
Phys. Lett. {\bf B413} (1997) 369;  
{\it ibid.} {\bf B495} (2000) 369. 

\bibitem{cp-matter}
J. Arafune and J. Sato, Phys. Rev. {\bf D55} (1997) 1653;
J. Arafune, M. Koike and J. Sato, Phys. Rev. {\bf D56} (1997) 3093;
Erratum {\it ibid.} {\bf D 60} (1999) 119905;
K.~Dick, M.~Freund, M.~Lindner and A.~Romanino,
Nucl.\ Phys.\ B {\bf 562} (1999) 29;
O. Yasuda,   Acta. Phys. Polon. {\bf B 30} (1999) 3089;
M. Koike and J. Sato, Phys. Rev. {\bf D61} (2000) 073012;
Erratum {\it ibid.} {\bf D 62} (2000) 079903;
A.~Cervera, A.~Donini, M.~B.~Gavela, J.~J.~Gomez Cadenas,
P.~Hernandez, O.~Mena and S.~Rigolin,
Nucl.\ Phys. \ {\bf B579} (2000) 17
[{\it Erratum-ibid.}\ {\bf B593} (2000) 731];
S.~J.~Parke and T.~J.~Weiler,
Phys.\ Lett.\ {\bf B501} (2001) 106.


\bibitem {MNjhep01}
H. Minakata and H. Nunokawa, JHEP {\bf 0110} (2001) 001 .

\bibitem {nufact01}
H. Minakata and H. Nunokawa, ; to be published in 
Proceedings of NuFACT01.

\bibitem{B-Castell}
J. Burguet-Castell, M.B. Gavela, J.J. Gomez-Cadenas, P. Hernandez 
and O. Mena,
Nucl. Phys. B {\bf 608} (2001) 301.

\bibitem {JHF}
Y. Itow {\it et al.},
{\it The JHF-Kamioka neutrino project}, .


\bibitem {totsuka}
Y. Totsuka, Talk at Frontiers in Particle Astrophysics and Cosmology; 
EuroConference on Neutrinos in the Universe, Lenggries, Germany, 
September 29 - October 4, 2001.



\end{thebibliography}




\end{document}






