# Atmospheric Neutrino Tests of Neutrino Oscillation Mechanisms

preprint:   UM-P-98/45 UDHEP-09-98 TMUP-HEL-9812 September 1998 

Recent Super-Kamiokande data on the atmospheric neutrino anomaly are used to test various mechanisms for neutrino oscillations. It is found that the current atmospheric neutrino data alone cannot rule out any particular mechanism. Future long-baseline experiments should play an important role in identifying the underlying neutrino oscillation mechanism.

The atmospheric neutrino anomaly has been confirmed by the Super-Kamiokande experiment. The observed up-down asymmetries of the detected muons indicate that, with the three known neutrino flavours, the most likely solution to this anomaly is $`\nu _\mu \nu _\tau `$ oscillations, although significant mixing with $`\nu _e`$ cannot yet be excluded.

The phenomenon of neutrino flavour oscillations was first proposed as a consequence of nondegenerate neutrino masses . Although this mass mixing of weak eigenstates is the most likely mechanism for the observed atmospheric neutrino oscillations, other possible mechanisms have been proposed and, as stressed in Ref., it is important to let experiments rather than our theoretical prejudice determine which is the correct mechanism. It is in this spirit that we undertake the present study.

The alternative neutrino oscillation mechanisms considered below share a common feature: each requires the existence of an interaction (other than the neutrino mass terms) that can mix neutrino flavours. In this paper we shall focus on such interactions which are mediated by a scalar ($`J=0`$), a vector ($`J=1`$), or a tensor ($`J=2`$) field. Assuming a two-neutrino mixing scheme, the $`\nu _\mu \nu _\tau `$ transition probability for each of these possibilities may be parametrized as 

$$P(\nu _\mu \nu _\tau )=\mathrm{sin}^22\theta _J\mathrm{sin}^2(E_\nu ^{J1}L\delta _J),$$
(1)   
where $`E_\nu `$ is the neutrino’s energy, $`L`$ is the neutrino’s path length (i.e., the distance from where the neutrino is produced to the detector), $`\delta _J`$ is a parameter specific to the neutrino oscillation mechanism, and $`\theta _J`$ is the corresponding mixing angle. The subscript $`J`$ in $`\theta _J`$ and $`\delta _J`$ is simply a label for the different mechanisms, and does not imply that the values of these parameters depend on the value of $`J`$. Note the distinct energy dependence of the oscillation probability for each value of $`J`$. This is the key for determining which neutrino oscillation mechanism is at work.

Proposed mechanisms for the $`J=0`$ case include the mass mixing mechanism for which the parameter $`\theta _0`$ in Eq.(1) represents the mixing angle between the neutrino flavour eigenstates and the mass eigenstates, and the parameter $`\delta _0`$ is given by

$$\delta _0(\mathrm{mass})=\frac{\mathrm{\Delta }m^2}{4}\frac{m_2^2m_1^2}{4}_,$$
(2)   
where $`m_i`$ are the masses of the neutrino mass eigenstates. Another possibility for the $`J=0`$ case is neutrino oscillations induced by nonuniversal couplings of the neutrinos to a massless string dilaton for which $`\delta _0`$ is given by

$$\delta _0(\mathrm{dilaton})=\frac{1}{4}\left[\mathrm{\Delta }m^22\varphi (\alpha _2m_2^2\alpha _1m_1^2)\right].$$
(3)   
Here $`\alpha _i`$ are the coupling strengths of different neutrino gravitational eigenstates to the dilaton field and it has been assumed that the neutrino mass eigenstates and the gravitational eigenstates are the same. The angle $`\theta _0`$ for this case is therefore the mixing angle between the flavour eigenstates and the gravitational eigenstates. Note that the nonuniversal neutrino-dilaton couplings constitute a violation of the principle of equivalence. The parameter $`\varphi `$ in Eq.(3) denotes the local Newtonian gravitational potential. It is assumed to be constant over the neutrino path length since the dominant contribution to $`\varphi `$ appears to come from the great attractor which has been estimated to be: $`|\varphi |3\times 10^5`$. With the constant $`\varphi `$ approximation, the energy dependence in the oscillation probability is the same for this mechanism as it is for the mass mixing mechanism. Neutrino oscillation experiments alone will not be able to distinguish these two possibilities. For the purposes of this paper, we shall refer to them collectively as the scalar mechanism. Note, however, that the mass mixing mechanism requires the neutrino masses to be nondegenerate whereas the dilaton induced oscillations can take place even if the masses are degenerate. We also observe that, if the neutrino-dilaton couplings are universal, i.e., if $`\alpha _1=\alpha _2`$, and if $`\delta _0(\mathrm{mass})0`$, $`\delta _0(\mathrm{dilaton})`$ and $`\delta _0(\mathrm{mass})`$ differ only by a multiplicative constant.

As an example of the $`J=2`$ case, consider neutrino oscillations induced by the equivalence principle violation that results from nonuniversal couplings of neutrinos to gravity. In this case $`\delta _2`$ is given by

$$\delta _2(\mathrm{VEP})=|\varphi |\mathrm{\Delta }\gamma ,$$
(4)   
where $`\mathrm{\Delta }\gamma `$ measures the degree of violation of the equivalence principle and, as in Eq.(3), $`\varphi `$ denotes the essentially constant local gravitational potential. With a constant $`\varphi `$, this mechanism is phenomenologically identical to the velocity oscillation mechanism that arises from a breakdown of Lorentz invariance for which $`\delta _2`$ assumes the form

$$\delta _2(\mathrm{velocity})=\frac{\mathrm{\Delta }v}{2}_,$$
(5)   
where $`\mathrm{\Delta }v=v_2v_1`$ is the difference between the velocities of two neutrino velocity eigenstates. We shall refer to these two mechanisms collectively as the tensor mechanism. In the former case, the angle $`\theta _2`$ corresponds to the mixing angle between the flavour eigenstates and the gravitational eigenstates, whereas in the latter case $`\theta _2`$ corresponds to the mixing angle between the flavour eigenstates and the velocity eigenstates. It should be emphasized that the tensor mechanism can lead to neutrino oscillations even if neutrinos are massless (or degenerate). These mechanisms were proposed not so much as a competing mechanism for the mass mixing mechanism, but to point out that neutrino oscillation experiments could be used as high-precision tests of the symmetry principles fundamental to the theories of general and special relativity.

We include the vector mechanism (the $`J=1`$ case in Eq.(1)) in our phenomenological study. Although there are models which can yield energy-independent neutrino oscillations in matter, we are not aware of any such model that can explain the atmospheric neutrino anomaly. The aim of our study is not to test any specific model, but rather to check if an energy-independent oscillation mechanism is compatible with the atmospheric neutrino data.

To test the three classes of neutrino oscillation mechanisms, we compare the measured values of

$$R\frac{(N_\mu /N_e)|_{\mathrm{data}}}{(N_\mu /N_e)|_{\mathrm{MC}}}$$
(6)   
and the up-down asymmetry parameters

$$Y_\alpha ^\eta \frac{(N_\alpha ^\eta /N_\alpha ^{+\eta })|_{\mathrm{data}}}{(N_\alpha ^\eta /N_\alpha ^{+\eta })|_{\mathrm{MC}}}_,\alpha =e,\mu $$
(7)   
with the corresponding predictions of these mechanisms, assuming $`\nu _\mu \nu _\tau `$ transitions. Here $`N_e`$ and $`N_\mu `$ are the number of $`e`$-like and $`\mu `$-like events, respectively. $`N_\alpha ^\eta `$ and $`N_\alpha ^{+\eta }`$ are the number of $`\alpha `$-like events produced in the detector with zenith angle $`\mathrm{cos}\mathrm{\Theta }<\eta `$ and $`\mathrm{cos}\mathrm{\Theta }>\eta `$, respectively. Note that $`\mathrm{\Theta }`$ is defined to be negative for upward going directions and we have chosen $`\eta =0.2`$ in our analysis. The calculational method is identical to that described in  and will not be reproduced here for environmental reasons.

Our results for the $`R`$’s and $`Y`$’s are displayed in Figures 1 and 2 together with the Super-Kamiokande results ,

$`R(\mathrm{sub}\mathrm{GeV})`$ $`=`$ $`0.63\pm 0.03\pm 0.05,`$ (8)   
$`R(\mathrm{multi}\mathrm{GeV})`$ $`=`$ $`0.65\pm 0.05\pm 0.08,`$ (9)   
$`Y_\mu ^{0.2}(\mathrm{sub}\mathrm{GeV})`$ $`=`$ $`0.79\pm 0.05,`$ (10)   
$`Y_\mu ^{0.2}(\mathrm{multi}\mathrm{GeV})`$ $`=`$ $`0.56\pm 0.06.`$ (11)   
The experimental results we use correspond to 535 live days of running. For completeness we mention that the Super-Kamiokande results for the $`e`$-like up-down asymmetries are

$`Y_e^{0.2}(\mathrm{sub}\mathrm{GeV})`$ $`=`$ $`1.14\pm 0.07,`$ (12)   
$`Y_e^{0.2}(\mathrm{multi}\mathrm{GeV})`$ $`=`$ $`0.92\pm 0.12.`$ (13)   
Finally, note that only statistical errors are given for the up-down asymmetries since they should be much larger than possible systematic errors at the moment. We have not shown any figure for $`Y_e`$ \- it is expected to be about 1 because $`\nu _\mu \nu _\tau `$ oscillations have almost no effect on the expected number of electron events.

From Figures 1 and 2 we see that all three mechanisms can fit the data for a range of the parameter $`\delta _J`$. The combined $`\chi ^2`$ fit to $`Y_\mu ,Y_e`$ and $`R`$ yield the allowed regions for the neutrino oscillation parameters shown in Figure 3. The $`\chi ^2`$ function is defined to be

$$\chi _{\mathrm{atm}}^2\chi ^2(R)+\chi ^2(Y),$$
(14)   
where

$$\chi ^2(R)\underset{E}{}\left[\left(\frac{R^{SK}R^{th}}{\delta R^{SK}}\right)^2\right]$$
(15)   
and

$$\chi ^2(Y)\underset{E}{}\left[\left(\frac{Y_\mu ^{SK}Y_\mu ^{th}}{\delta Y_\mu ^{SK}}\right)^2+\left(\frac{Y_e^{SK}Y_e^{th}}{\delta Y_e^{SK}}\right)^2\right]_.$$
(16)   
The sum is over the sub-GeV and multi-GeV data samples. The measured Super-Kamiokande values and errors are denoted by the superscript “SK” and the theoretical predictions for the corresponding quantities are labelled by the superscript “th”. The minimal $`\chi ^2`$ for the scalar, vector, and tensor mechanisms is 5.4, 5.7, and 10.4, respectively, for 4 degrees of freedom. All three mechanisms are consistent with the data, but the tensor mechanism has the worst $`\chi ^2`$.

For check we have performed a second $`\chi ^2`$ analysis using the SK $`\chi ^2`$ function,

$$\chi ^2=\frac{\alpha ^2}{\sigma _\alpha ^2}+\frac{\beta ^2}{\sigma _\beta ^2}+\chi _{\mathrm{sub}\mathrm{GeV}}^2+\chi _{\mathrm{multi}\mathrm{GeV}}^2,$$
(17)   
where

$$\chi _{\mathrm{sub}\mathrm{GeV}}^2=\underset{\alpha =e,\mu }{}\underset{a=1}{\overset{5}{}}\frac{\left(y_a^\alpha n_a^\alpha \right)^2}{n_a^\alpha }$$
(18)   
and similarly for $`\chi _{\mathrm{multi}\mathrm{GeV}}^2`$. Here $`\alpha `$ and $`\beta `$ are the absolute and relative normalizations of the neutrino flux, $`y_a^\alpha `$ are theoretical predictions for each zenith angle bin, $`n_a^\alpha `$ are the experimental data, and we assume the uncertainties are given by $`\sigma _\alpha =\mathrm{}`$ and $`\sigma _\beta =0.12`$ (we assumed the larger value of $`\sigma _\beta `$ given in the third reference of ). It is understood that $`y_a^\mu `$ ($`\mu `$-like events) are multiplied by $`(1+\alpha )(1+\beta /2)`$ while $`y_a^e`$ (e-like events) are multiplied by $`(1+\alpha )(1\beta /2)`$, and $`\chi ^2`$ is optimized with respect to $`\alpha `$ and $`\beta `$. With this alternative $`\chi ^2`$, we also found that all three cases ($`J=0,1,2`$) gave an acceptable fit to the data.

The basic reason that the Super-Kamiokande atmospheric neutrino data cannot distinguish these rather different oscillation probabilities is geometrical. Neutrinos from above typically travel quite short distances ($`50\mathrm{km}`$) whereas neutrinos going up through the Earth travel quite long distances ($`5000\mathrm{km}`$). The atmospheric neutrino data can be explained by assuming that neutrinos from above do not have time to oscillate while neutrinos travelling through the Earth experience averaged oscillations. In fact the data suggest that this occurs for both the sub-GeV and multi-GeV energy range (which is, very roughly, $`0.3\mathrm{GeV}E_\nu 6\mathrm{GeV}`$). Thus, while atmospheric neutrino experiments can sensitively test for oscillations through the zenith angle dependence (or up-down asymmetry) they are really only sensitive to averaged oscillations and are consequently not sensitive to the oscillation mechanism. The explicit energy dependence of the oscillation probability is only important in the region where the oscillations are significant and are not averaged. The proposed long-baseline experiments should be much more sensitive to the oscillation mechanism. The reason is that the neutrino path length is fixed and at a moderate distance (typically about 250 - 600 km). Provided that not all of the oscillations are averaged, and provided that at least some of the neutrinos have time to oscillate, then a comparison of the detected energy distribution with the expected energy distribution should distinguish the oscillation mechanisms.

In summary, we have examined three distinct oscillation mechanisms ($`J=0,1,2`$ in Eq.(1)) and compared these with the current Super-Kamiokande atmospheric neutrino data. We have found that the current data are quite insensitive to the oscillation mechanism responsible for the atmospheric neutrino anomaly. To ultimately determine the underlying neutrino oscillation mechanism, it is important to do a more controlled experiment. This should be possible in the near future with the advent of long-baseline experiments.

Acknowledgement

C.N.L. was supported in part by the U.S. Department of Energy under grant DE-FG02-84ER40163. O. Y. was supported in part by a Grant-in-Aid for Scientific Research of the Ministry of Education, Science and Culture, #09045036, #10140221 and #10640280. R. F. is an Australian Research Fellow.

FIGURE CAPTIONS

Figure 1: $`R`$ as a function of $`\delta _J`$ (see Eq.(1)) for the tensor (bold solid line), vector (solid line), and scalar (dashed line) mechanisms. Maximal mixing has been assumed in each case. The dimensionless variable X is equal to $`\delta _2`$ GeV-km, $`\delta _1`$ km, and $`\frac{\delta _0}{1.27}\mathrm{GeV}^1`$-km, respectively for the tensor, vector, and scalar mechanisms (the choice of X for the mass mixing mechanism corresponds to the familiar $`\mathrm{\Delta }m^2/\mathrm{eV}^2`$). The horizontal dashed lines correspond to Super-Kamiokande’s measured $`R`$ value $`\pm `$ 1$`\sigma `$ error. Figure 1a is for the sub-GeV neutrinos and Figure 1b is for the multi-Gev neutrinos.

Figure 2: Same as Figure 1 except $`Y_\mu `$ (see text) is plotted instead of $`R`$.

Figure 3: The allowed region for the mixing parameter $`\theta _J`$ and the oscillation parameter $`\delta _J`$ obtained with the $`\chi ^2`$ function defined in Eq.(14). Figure 3a, 3b, and 3c are for the scalar, vector, and tensor mechanisms, respectively.