Signals of CPT Violation and Non-Locality
[2mm] in Future Neutrino Oscillation Experiments
[8mm] S. Antusch 111E-mail: firstname.lastname@example.org, E. Fernandez-Martinez 222E-mail: email@example.com,
Max-Planck-Institut für Physik (Werner-Heisenberg-Institut),
Föhringer Ring 6, D-80805 München, Germany
We investigate the sensitivities of future neutrino oscillation experiments for measuring the neutrino mass squared differences and leptonic mixing angles independently with neutrinos and anti-neutrinos. We update the expected sensitivities of Neutrino Factories to the “atmospheric” (anti-)neutrino parameters using an optimized setup. A dedicated -Beam facility, in combination with a SPMIN reactor experiment, could give excellent sensitivities also to the “solar” parameters, for neutrinos and anti-neutrinos respectively. A signal of a different mass matrix for neutrinos and anti-neutrinos would imply CPT violation and non-locality of the underlying particle theory.
With the expected high sensitivities of envisioned neutrino facilities like Neutrino Factories - or -Beams , neutrino physics could enter a new era of precision . One main goal of these experiments is the measurement of the yet unknown low energy parameters of the lepton sector accessible to oscillation physics, i.e. the Dirac CP phase , the mixing angle and the sign of the “atmospheric” mass squared difference (corresponding to whether the mass spectrum is normally ordered or inverted). However, in addition to the determination of these parameters, future precision experiments may also find interesting signals of physics beyond the SM or even signals of violation of fundamental principles such as CPT invariance, Lorentz invariance or locality . One interesting signal of this type would be a difference between the masses and mixing angles measured by neutrino oscillation experiments operating with neutrinos and anti-neutrinos.
The CPT theorem  states that any local Quantum Field Theory (QFT) which is Lorentz invariant and has a Hermitean Hamiltonian must have CPT symmetry. From the derivation of the CPT theorem it also follows that CPT violation implies violation of Lorentz invariance. Bounds on CPT and Lorentz invariance violation from processes involving quarks and charged leptons are quite strong, for example from the system . On the other hand, CPT and Lorentz invariance are tested to a much less precision in the neutrino sector by the current experiments, as we will review below in section 2.1. In addition to the less tight constraints, another motivation to look for signals of CPT and Lorentz invariance violation in the neutrino sector is the potentially different mechanism for generating the small neutrino masses compared to the masses of quarks and charged leptons, which might be especially sensitive to new physics.
One direct consequence of CPT invariance is that neutrinos and anti-neutrinos have the same masses and mixing angles. It has been pointed out in Ref.  that a signal for the violation of this prediction, i.e. different masses and mixing angles for neutrinos and anti-neutrinos, would allow to draw the additional conclusion that in addition to CPT and Lorentz invariance also locality must be violated. More precisely, it has been shown that a difference between the masses for particles and anti-particles implies non-locality in the sense that field (anti-)commutators do not vanish for space-like distances and furthermore that the fields enter terms in the Lagrangian at different space-time points . Non-locality is in general predicted by extensions of the SM towards a unified theory with gravity. Regarding CPT violation, it has been argued in Ref.  that non-local interactions in string theory might generate CPT violation close to present bounds. Specific scenarios of CPT breaking can also be realised in noncommutative geometries . We would like to note, however, that if a signal of CPT violation via different masses of neutrinos and anti-neutrinos would be observed, it could be challenging to distinguish the possible intrinsic CPT violation and non-locality from “fake” signals, caused e.g. by non-standard matter effects. In any case, if the experimental data would point to different masses and/or mixing angles for neutrinos and anti-neutrinos, it would be an intriguing signal of physics beyond the SM.
Various aspects regarding CPT violation in the neutrino sector have been analysed in previous studies. For example, Ref.  has discussed the potential to test CPT by combining solar neutrino data and data from the KamLAND reactor experiment. The sensitivity of experiments with atmospheric neutrinos has been discussed in  and the prospects of the MINOS experiment has been studied in . For the “atmospheric” parameters, high sensitivity to CPT violation could be achieved in Neutrino Factories [15, 16]. Additional tests might be possible with the data from supernova neutrinos  and with neutrinoless double beta decay . One motivation for the study of models of CPT violation for neutrino physics has been the LSND anomaly, where CPT violation in the neutrino sector was proposed as a solution . The case studied by most authors in this context is that of a CPT violating background which leads to CPT and Lorentz invariance violating local effective operators at low energy . Recently, a parameterisation for a three family oscillation analysis in the presence of specific types of CPT violation has been studied in .
In order to search for a signal of CPT violation and non-locality in neutrino oscillation experiments a precise determination of the neutrino mass splittings and/or leptonic mixing angles, independently with neutrinos and anti-neutrinos, is desirable. Recently, significant progress has been made in the design of possible future neutrino oscillation facilities, for example regarding Neutrino Factories . Furthermore, -Beam facilities have been proposed which, in some aspects, could even allow more precise measurements than a Neutrino Factory. In addition, new reactor experiments (for  as well as for precision measurements of -) have been envisioned. In this letter, we therefore investigate the combined sensitivity of such future neutrino oscillation experiments for measuring the “solar” and “atmospheric” mass squared differences and leptonic mixing angles independently with neutrinos and anti-neutrinos. We update the expected sensitivities for the “atmospheric” mass squared difference and mixing angle with Neutrino Factories using an optimized setup and analyse how a dedicated -Beam facility, in combination with a SPMIN reactor experiment -, could improve the sensitivities for the “solar” mass splitting and mixing angle.
2 Potential signals of CPT violation and non-locality
In the following we denote parameters for anti-neutrinos with bars and parameters for neutrinos without bars. Using standard PDG parameterisation of the PMNS matrix , CPT invariance implies and . As discussed above, a violation of these equalities, i.e. for one neutrino mass eigenvalue or for one of the mixing angles, if not confused with a “fake” signal of different new physics, would signal violation of CPT invariance (and Lorentz invariance) as well as non-locality of the underlying particle theory. To analyse the expected sensitivities to such signals of new physics beyond the SM, we study how well neutrino oscillation experiments can determine the neutrino masses and leptonic mixing angles for neutrinos and anti-neutrinos, i.e. using only data from oscillation of neutrinos and anti-neutrinos, respectively. In our analysis we treat the parameters and as independent quantities and assume the standard three family oscillation formulae for neutrinos and anti-neutrinos , thereby testing the consistency of the CPT symmetric description. Bounds on CPT violation are presented as the constraints on differences for the mixing angles and for the mass squared differences (defined as and ) which the considered experiments can impose. Before we turn to the discussion of the expected sensitivities of future neutrino oscillation facilities in this respect, we review the bounds on the difference between neutrino and anti-neutrino parameters from the present data. A graphical summary of the present bounds can be found in Fig. 55 of .
2.1 Present bounds (anti-)neutrino parameters
Regarding the “solar” neutrino parameters , and , , the most relevant data stems from experiments with solar neutrinos - and from the KamLAND experiment  observing anti-neutrinos from nuclear reactors. If CPT invariance is assumed, the complementarity between these two data sets allows for a quite good determination of the “solar” mixing angle as well as of the “solar” mass squared difference. However, separating into neutrino and anti-neutrino parameters shows that the bounds on their difference are comparatively weak. While solar neutrino data allows for a quite precise measurement of , there is very little information on . The reason for the good sensitivity on is that solar matter effects adiabatically convert the produced in the solar core into when leaving the sun, and thus detectors on earth can extract by comparing the measured flux to the theoretically expected flux without oscillations or to the total neutrino flux measured through neutral current interactions. The sensitivity of the solar neutrino experiments to the mass splitting is very poor because the distance between earth and sun is so large that the oscillatory behaviour of the transition probability is averaged out. KamLAND, on the other hand, has a much shorter average baseline of Km. The (almost) vacuum oscillation signal of the reactor and the measurement of the distortion of the anti-neutrino energy spectrum allows a precise measurement of . However, the overall normalization of the neutrino flux is rather difficult due to the fact that anti-neutrinos from many reactors (as well as from other sources) are detected by KamLAND. Consequently the constraints on from KamLAND is comparatively weak. Furthermore, the near vacuum oscillations observed by KamLAND depend on the quantity , which implies that the octant of cannot be determined from the present data on anti-neutrinos. In summary, assuming that lies in the first octant, the present bounds on CPT violation in the “solar” sector are at
Regarding the “atmospheric” parameters , and , , the most important experimental data stems from atmospheric neutrinos observed by SuperKamiokande [34, 35] and from the accelerator experiments sensitive to the “atmospheric” sector, K2K  and MINOS . K2K and MINOS operate with neutrinos and the atmospheric signal is dominated by the neutrino event rates due to the larger cross section. Therefore, the measurements of and in particular of are much less precise than the ones of and (see Fig. 55 of ). The present bounds on CPT violation for the “atmospheric” sector are at [5, 38]:
Finally the bounds on the unknown 1-3 mixing angles and for neutrinos and anti-neutrinos, respectively, are similar. The bound (at ) is obtained from the CHOOZ experiment  and a bound of can be derived from the combined data of experiments on solar, atmospheric and accelerator neutrino oscillations . From the present data, without a measurement of neither nor , the bound on the difference (at ) is given by
We will now discuss how these bounds (or equivalently the sensitivities for discovering a signal) might be improved in future neutrino oscillation facilities.
2.2 Strategies for improving the sensitivities
As has been discussed in previous studies , the bounds on and can be strongly improved at a Neutrino Factory facility -. Excellent sensitivities to are also to be expected since the determination of is one of the main goals of the Neutrino Factory. In the context of the IDS, International Design Study of the Neutrino Factory, a baseline for the Neutrino Factory acceleration complex and detection systems has recently been defined [5, 40]. This baseline setup would store 25 GeV muons, whose (anti-)neutrino fluxes aim at two 50 Kton Magnetized Iron Neutrino Detectors (MIND) located at Km and Km from the source, respectively. The goal luminosity for such a facility is useful muon decays per year per polarity per baseline. The efficiencies, backgrounds and energy resolution of the MIND detector when exposed to such a beam have been studied in . We will present the up-dated bounds with these specifications in section 2.3. In comparison, little work has been done so far with respect to the improvement of the bounds on the “solar” parameters, and we will therefore focus on this issue in the remainder of this section.
As discussed in section 2.1, the present uncertainty on is dominated by the comparatively low precision on , while the error on stems from the weak constraints on . Improving the former is much easier than the latter. A dedicated reactor experiment placed at the minimum of the survival probability (SPMIN) could provide an excellent measurement of -. A precise determination of , however, is much more challenging. The observation of the oscillations driven by the small “solar” splitting requires large values of . The neutrino charged current (CC) interaction, however, is very suppressed at low energies since the neutrons with which they can interact are confined in nuclei and a minimum threshold energy is normally required for the interaction to occur. Moreover, while nuclear reactors provide a convenient source of electron anti-neutrinos, it is challenging to produce low energy electron neutrinos in large amounts on earth.
One possibility to improve the determination of would be an intense beam of GeV from the decay of Ne ions accelerated to at a -Beam facility ,-. However, to optimize the determination of and with neutrinos, we consider also a very long baseline of Km to a Mton class water Cerenkov detector -. Somewhat less sensitivity can be reached with a more conventional 750 Km baseline. We have also studied the possible measurements achievable with the that this facility could also provide from the decays of He. However we found that the performance of the anti-neutrino beam was significantly worse. This asymmetry is caused by the different matter effects for neutrinos and anti-neutrinos. In particular, the electron (anti-)neutrino survival probabilities for are given by (see e.g. ):
with and where the signs apply to the probabilities of anti-neutrinos and neutrinos, respectively. Thus, matter effects reduce the amplitude of the oscillation probability and increase the oscillation frequency. Due to the sign in Eq. (2.7) for anti-neutrinos the effect is enhanced, the oscillation amplitude is smaller and a less precise reconstruction of the parameters follows. This situation is illustrated in Fig. 1. The solid and dashed lines correspond to the and survival probabilities respectively. Consequently, present bounds on and will not be improved by this setup. On the other hand, Eq. (2.7) is no longer symmetric under the octant of . Changing the octant has the same effect as changing the sign of or, equivalently, interchanging the neutrino and anti-neutrino curves in Fig. 1. This means that the from the He decay could distinguish the octant of which would not be possible even at a dedicated SPMIN experiment for a precise measurement of due to the small matter effects.
The -Beam setup with 4000 Km baseline is thus optimized for the measurement of “solar” parameters through disappearance, however it is unrealistic to expect that such a dedicated facility would ever be built. As already mentioned above, we have therefore also considered a setup with a shorter, more conventional baseline of Km. This is the baseline usually chosen for the version of the -Beam, which could provide excellent sensitivities to CP-violation  even outperforming the Neutrino Factory for some regions of the parameter space. A run of the -Beam at this baseline, such as the one we consider here to improve the measurement on , can improve the sensitivity to the mass hierarchy and allow to solve degeneracies when combined with the higher run .
2.3 Expected sensitivities of future facilities
For the experimental setup we consider the “IDS baseline” Neutrino Factory , which provides excellent sensitivities to the “atmospheric” parameters and with both neutrino and anti-neutrino beams. In this facility intense (anti-)neutrino beams would be produced from the decay of GeV (anti-)muons with a luminosity of muon decays per year per muon polarity. This beams illuminate two identical 50 Kton iron calorimeters located at 4000 Km and 7000 Km. We consider 5 years of data taking with each muon polarity. We did not consider the Opera-like detector to observe the silver channel  since it did not improve significantly any bound.
In order to improve the present bounds on the “solar” sector we also considered a -Beam facility producing electron (anti-)neutrinos from the decay of Ne (He) ions accelerated to . We considered 5 years of data taking with luminosities of ion decays per year  for both ions. These neutrino beams would be detected at a Mton class water Cerenkov detector at a 750 or 4000 Km baseline. To describe the detector efficiencies and backgrounds when exposed to these beams we followed Ref. .
Finally a reactor experiment with an exposure of GWKtonyr to a detector located at Km distance was considered to improve the measurement of the “solar” parameters with anti-neutrinos. We considered a 5% systematic error in the normalization of the signal and the energy resolution of the detector was set to . The signal was distributed in 10 energy bins in the range between and GeV.
|quantity||present bound||future (B 4000 Km)||future (B 750 Km)|
For the numerical analysis we independently combined the neutrino and anti-neutrino data from the Neutrino Factory, SPMIN and -Beam setups to derive measurements on the neutrino and anti-neutrino parameters. Since the combination of the three facilities allows to constrain all the oscillation parameters no prior information on any of them was assumed except for the measurement that will not be improved by this setup, namely from solar neutrino oscillations, for which we assumed a uncertainty of 10 %. A 5 % uncertainty in the PREM density profile was also assumed. The Globes 3.0  software was used to perform the numerical analysis. The following (CPT conserving) input values were assumed for the oscillation parameters: , , , , .
Figures 2 and 3 show the 90 %, 95 %, 99 % and contours for the constraints on , , , and with either neutrino data (vertical axes) or anti-neutrino data (horizontal axes) alone. The left (right) panels contain the constraints achievable with the 4000 Km (750 Km) baseline for the -Beam. Comparing it with Fig. 55 of  the dramatic improvement on the constraints on the “atmospheric” parameters with the Neutrino Factory is manifest. The constraints for and shrink by about one order of magnitude, while the improvement in and is about two and three orders of magnitude respectively. The bounds on CPT violation that could be derived from this precision measurements including the 750 Km baseline -Beam would be:
With such a Neutrino Factory setup, the constraints on both and could also be significantly improved by almost three orders of magnitude to
In the “solar” sector the improvements are more modest, especially for the shorter -Beam baseline. For the 4000 Km baseline matter effects allow to measure the octant of , but at 750 Km their strength is not sufficient for this task, as can be seen in the top panels of Fig. 2. With neutrinos an improvement of by an order of magnitude could be accomplished with the -Beam and some improvement in the lower bound of could also be achieved. These improved measurements could test CPT invariance to the level of
for the 4000 Km baseline -Beam or
for the 750 Km baseline -Beam. In the latter case the octant of is not measured. has been assumed to lie in the first octant in order to compare with Eq. (2.2). A summary of present and possible future bounds is presented in Tab. 1.
3 Discussion and Conclusions
In this study we have investigated the sensitivity of future neutrino oscillation experiments for measuring the “solar” and “atmospheric” mass squared differences and leptonic mixing angles independently with neutrinos and anti-neutrinos. If a difference between the parameters for neutrinos and anti-neutrinos would be established, it would imply CPT (and Lorentz invariance) violation as well as non-locality.
To improve the present bounds on this form of CPT violation we have considered three types of possible future neutrino oscillation experiments: an optimized Neutrino Factory, a -Beam (with 750 Km or 4000 Km baseline) pointing to a Mton water Cerenkov detector and a SPMIN reactor experiment at the minimum of the survival probability. The possible improvements with respect to the present bounds are summarised in Tab. 1.
Regarding the “atmospheric” parameters and (vs. ), a dramatic improvement could be accomplished with a Neutrino Factory operating with neutrinos as well as with anti-neutrinos. The sensitivity on could be improved almost three orders of magnitude, the sensitivity on by one order of magnitude and the sensitivity on by more than two orders of magnitude.
To improve the sensitivity with respect to the “solar” parameters, the combination of the -Beam with the SPMIN reactor experiment could also be very successful. While the -Beam experiment could improve the measurements of the parameters with neutrinos (mainly ) and determine the octant of , the SPMIN experiment could further improve the measurement of . While the gain in sensitivity on amounts about 50 %, the sensitivity on could be improved by about one order of magnitude.
In summary, future neutrino oscillation experiments have the potential to measure neutrino mass squared differences and leptonic mixing angles separately with neutrinos and anti-neutrinos to high precision and may detect possible differences between neutrino and anti-neutrino parameters. Such a difference would be a clear signal of new physics beyond the Standard Model of particle physics. If other types of new physics which may lead to “fake” signals, e.g. non-standard matter effects, could be excluded it would signal CPT (and Lorentz invariance) violation and also a non-local nature of the underlying particle theory.
We would like to thank Andrea Donini for helpful discussions and Patrick Huber for kindly providing the description of the IDS baseline for the Neutrino Factory. This work was partially supported by The Cluster of Excellence for Fundamental Physics “Origin and Structure of the Universe” (Garching and Munich).
-  S. Geer, Phys. Rev. D 57 (1998) 6989 [Erratum-ibid. D 59 (1999) 039903] [arXiv:hep-ph/9712290]; A. De Rujula, M. B. Gavela and P. Hernandez, Nucl. Phys. B 547 (1999) 21 [arXiv:hep-ph/9811390].
-  A. Cervera, A. Donini, M. B. Gavela, J. J. Gomez Cadenas, P. Hernandez, O. Mena and S. Rigolin, Nucl. Phys. B 579 (2000) 17 [Erratum-ibid. B 593 (2001) 731] [arXiv:hep-ph/0002108].
-  M. Apollonio et al., arXiv:hep-ph/0210192.
-  P. Zucchelli, Phys. Lett. B 532 (2002) 166.
-  A. Bandyopadhyay et al. [ISS Physics Working Group], arXiv:0710.4947 [hep-ph].
-  See e.g.: V. A. Kostelecky and M. Mewes, Phys. Rev. D 69 (2004) 016005 [arXiv:hep-ph/0309025], and references therein.
-  J. Schwinger, Phys. Rev. 82 (1951) 914; G. Lüders, Kgl. Danske Videnskab. Selskab. Mat.-Fys. Medd. 28, No. 5 (1954); W. Pauli, in Niels Bohr and the Development of Physics, edited by W. Pauli, L. Rosenfeld, and V. Weisskopf (Pergamon, London, 1955), p. 30; J.S. Bell, Proc. R. Soc. London A 231, 479 (1955).
-  Particle Data Group Collaboration, W. M. Yao et al., Review of particle physics, J. Phys. G33 (2006) 1–1232.
-  O. W. Greenberg, Phys. Rev. Lett. 89 (2002) 231602 [arXiv:hep-ph/0201258].
-  V. A. Kostelecky and R. Potting, Nucl. Phys. B 359 (1991) 545; J. R. Ellis, N. E. Mavromatos and D. V. Nanopoulos, Int. J. Mod. Phys. A 11 (1996) 1489 [arXiv:hep-th/9212057]; J. R. Ellis, J. L. Lopez, N. E. Mavromatos and D. V. Nanopoulos, Phys. Rev. D 53 (1996) 3846 [arXiv:hep-ph/9505340].
-  G. Amelino-Camelia and S. Majid, Int. J. Mod. Phys. A 15 (2000) 4301 [arXiv:hep-th/9907110]. I. Mocioiu, M. Pospelov and R. Roiban, Phys. Rev. D 65 (2002) 107702 [arXiv:hep-ph/0108136].
-  J. N. Bahcall, V. Barger and D. Marfatia, Phys. Lett. B 534 (2002) 120 [arXiv:hep-ph/0201211].
-  A. Datta, R. Gandhi, P. Mehta and S. Uma Sankar, Phys. Lett. B 597 (2004) 356 [arXiv:hep-ph/0312027].
-  B. J. Rebel and S. L. Mufson [for the MINOS Collaboration], arXiv:0802.3785 [hep-ph].
-  V. D. Barger, S. Pakvasa, T. J. Weiler and K. Whisnant, Phys. Rev. Lett. 85 (2000) 5055 [arXiv:hep-ph/0005197].
-  S. M. Bilenky, M. Freund, M. Lindner, T. Ohlsson and W. Winter, Phys. Rev. D 65, 073024 (2002) [arXiv:hep-ph/0112226].
-  H. Minakata and S. Uchinami, Phys. Rev. D 72 (2005) 105007 [arXiv:hep-ph/0505133].
-  G. Barenboim, J. F. Beacom, L. Borissov and B. Kayser, Phys. Lett. B 537 (2002) 227 [arXiv:hep-ph/0203261].
-  H. Murayama and T. Yanagida, Phys. Lett. B 520 (2001) 263 [arXiv:hep-ph/0010178].
-  D. Colladay and V. A. Kostelecky, Phys. Rev. D 55 (1997) 6760 [arXiv:hep-ph/9703464]; S. R. Coleman and S. L. Glashow, Phys. Rev. D 59 (1999) 116008 [arXiv:hep-ph/9812418]; for limits on explicit models, see also: I. Mocioiu and M. Pospelov, Phys. Lett. B 534 (2002) 114 [arXiv:hep-ph/0202160].
-  A. Dighe and S. Ray, arXiv:0802.0121 [hep-ph].
-  F. Ardellier et al. [Double Chooz Collaboration], arXiv:hep-ex/0606025.
-  A. Bandyopadhyay, S. Choubey and S. Goswami, Phys. Rev. D 67 (2003) 113011 [arXiv:hep-ph/0302243].
-  H. Minakata, H. Nunokawa, W. J. C. Teves and R. Zukanovich Funchal, Phys. Rev. D 71 (2005) 013005 [arXiv:hep-ph/0407326].
-  A. Bandyopadhyay, S. Choubey, S. Goswami and S. T. Petcov, Phys. Rev. D 72 (2005) 033013 [arXiv:hep-ph/0410283].
-  M. C. Gonzalez-Garcia and M. Maltoni, arXiv:0704.1800 [hep-ph].
-  B. T. Cleveland et al., Astrophys. J. 496 (1998) 505.
-  J. N. Abdurashitov et al. [SAGE Collaboration], Phys. Rev. C 60 (1999) 055801 [arXiv:astro-ph/9907113].
-  W. Hampel et al. [GALLEX Collaboration], Phys. Lett. B 447 (1999) 127.
-  S. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 86 (2001) 5651 [arXiv:hep-ex/0103032].
-  Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87 (2001) 071301 [arXiv:nucl-ex/0106015].
-  S. N. Ahmed et al. [SNO Collaboration], Phys. Rev. Lett. 92 (2004) 181301 [arXiv:nucl-ex/0309004].
-  K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90 (2003) 021802 [arXiv:hep-ex/0212021].
-  Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81 (1998) 1562 [arXiv:hep-ex/9807003].
-  M. Ambrosio et al. [MACRO Collaboration], Phys. Lett. B 517 (2001) 59 [arXiv:hep-ex/0106049].
-  E. Aliu et al. [K2K Collaboration], Phys. Rev. Lett. 94 (2005) 081802 [arXiv:hep-ex/0411038].
-  D. G. Michael et al. [MINOS Collaboration], Phys. Rev. Lett. 97 (2006) 191801 [arXiv:hep-ex/0607088].
-  M. C. Gonzalez-Garcia and M. Maltoni, Phys. Rev. D 70 (2004) 033010 [arXiv:hep-ph/0404085].
-  M. Apollonio et al. [CHOOZ Collaboration], Eur. Phys. J. C 27 (2003) 331 [arXiv:hep-ex/0301017].
-  See http://www.hep.ph.ic.ac.uk/ids/docs/index.html
-  T. Abe et al., arXiv:0712.4129 [physics.ins-det].
-  J. Bouchez, M. Lindroos and M. Mezzetto, AIP Conf. Proc. 721 (2004) 37 [arXiv:hep-ex/0310059]; M. Mezzetto, J. Phys. G 29 (2003) 1771 [arXiv:hep-ex/0302007]; Nucl. Phys. Proc. Suppl. 143 (2005) 309 [arXiv:hep-ex/0410083].
-  A. Donini, E. Fernández-Martínez, P. Migliozzi, S. Rigolin and L. Scotto Lavina, Nucl. Phys. B 710 (2005) 402 [arXiv:hep-ph/0406132], A. Donini, E. Fernández-Martínez and S. Rigolin, Phys. Lett. B 621 (2005) 276 [arXiv:hep-ph/0411402].
-  J. E. Campagne, M. Maltoni, M. Mezzetto and T. Schwetz, JHEP 0704, 003 (2007) [arXiv:hep-ph/0603172].
-  Y. Itow et al. [The T2K Collaboration], arXiv:hep-ex/0106019.
-  A. de Bellefon et al., arXiv:hep-ex/0607026.
-  M. Diwan et al., arXiv:hep-ex/0608023.
-  E. K. Akhmedov, R. Johansson, M. Lindner, T. Ohlsson and T. Schwetz, JHEP 0404 (2004) 078 [arXiv:hep-ph/0402175].
-  J. Burguet-Castell, D. Casper, J. J. Gomez-Cadenas, P. Hernandez and F. Sanchez, Nucl. Phys. B 695 (2004) 217 [arXiv:hep-ph/0312068]; J. Burguet-Castell, D. Casper, E. Couce, J. J. Gomez-Cadenas and P. Hernandez, Nucl. Phys. B 725 (2005) 306 [arXiv:hep-ph/0503021].
-  A. Donini and E. Fernandez-Martinez, Phys. Lett. B 641 (2006) 432 [arXiv:hep-ph/0603261]; E. Fernandez-Martinez, arXiv:hep-ph/0605101.
-  A. Donini, D. Meloni and P. Migliozzi, Nucl. Phys. B 646 (2002) 321 [arXiv:hep-ph/0206034]; J. Phys. G 29 (2003) 1865 [arXiv:hep-ph/0209240]; D. Autiero et al., Eur. Phys. J. C 33 (2004) 243 [arXiv:hep-ph/0305185].
See the talk by M. Lindroos at the second plenary IDS meeting.
-  P. Huber, J. Kopp, M. Lindner, M. Rolinec and W. Winter, Comput. Phys. Commun. 177 (2007) 432 [arXiv:hep-ph/0701187].