## Improved Theory of Neutrino Oscillations in Matter

## Leonard S. Kisslinger

Department of Physics, Carnegie Mellon University, Pittsburgh PA 15213 USA.

November 20, 2014

## Abstract

This is revision of the S-Matrix theory of neutrino oscillations used for many years. We evaluate the transition probability of a µ to e neutrino without an approximation used for many theoretical studies, and find important differences which could improve the extraction of neutrino parameters from experimental data in the future.

.PACS Indices:11.30.Er,14.60.Lm,13.15.+g

Keywords: neutrino, neutrino oscillations, S-Matrix theory

## 1 Introduction

In previous research on neutrino oscillations we used S-Matrix theory with three active neutrinos to evaluate time reversal violation[1], CP violation[2], and recently[3] the transition probability P ( ν µ → ν e ) . One of the main objectives of Ref.[3] was to calculate neutrino oscillations and anti-electron neutrino disappearance to help in the measurement of θ 13 from Daya Bay, Double Chooze, and RENO data. For this it is essential to have an accurate theory of neutrino transition probabilities.

Our previous research using S-Matrix theory with a 3 × 3 mixing matrix for three active neutrinos was based on earlier publications[4, 5, 6, 7]. In carrying out the evaluation of P ( ν µ → ν e ) we used an approximation for an important function, I α ∗ , found in this earlier work, such as Ref[6], whose formulation we used. In the present work we evalute I α ∗ exactly and find corrections to P ( ν µ → ν e ). Our objective is to improve P ( ν µ → ν e ), and thereby the extraction of important parameters from neutrino oscilation data.

## 2 P ( νµ → νe ) Derived Using S-Matrix Theory

Neutrinos are produced as ν f , with f =flavor= e, µ, τ . Neutrinos with definite mass are ν m , m=1,2,3. The flavor neutrinos are related to the neutrinos with definite mass by the 3 × 3 unitary matrix, U ,

$$\nu _ { f } \ = \ U \nu _ { m } \,,$$

where ν , ν f m are 3 × 1 column vectors and U is ( sinθ ij ≡ s ij , etc).

$$U = \left ( \begin{array} { c c c c } c _ { 1 2 } c _ { 1 3 } & s _ { 1 2 } c _ { 1 3 } & s _ { 1 3 } e ^ { - i \delta _ { C P } } \\ - s _ { 1 2 } c _ { 2 3 } - c _ { 1 2 } s _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & c _ { 1 2 } c _ { 2 3 } - s _ { 1 2 } s _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & s _ { 2 3 } c _ { 1 3 } \\ s _ { 1 2 } s _ { 2 3 } - c _ { 1 2 } s _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & - c _ { 1 2 } s _ { 2 3 } - s _ { 1 2 } c _ { 2 3 } s _ { 1 3 } e ^ { i \delta _ { C P } } & c _ { 2 3 } c _ { 1 3 } \end{array} \right ) \\ \mathbb { M } _ { \infty } \, \mathbb { m } _ { \infty } \, \varkappa \, \varepsilon \, \varpi \, \varpi \, = \, \mathbb { M } _ { \infty } \, \varkappa \, \varpi \, = \, \varpi \, \varpi \, = \, \varpi \, \varpi \, \varpi \, = \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \mathbb { M } _ { \infty } \, \cdots = \, \mathbb { M } _ { \infty }$$

We use c 12 = 83 . , s 12 = 56 . , s 23 = c 23 = 7071, . s 13 = 0.19, and δ CP =0.

The transition probability P ( ν µ → ν e ) is obtained from the S-Matrix element S 12 , where i d dt S t, t ( 0 ) = H t S t, t ( ) ( 0 ), with H t ( ) the Hamiltonian:

$$\mathcal { P } ( \nu _ { \mu } \to \nu _ { e } ) \ = \ ( R e [ S _ { 1 2 } ] ) ^ { 2 } + ( I m [ S _ { 1 2 } ] ) ^ { 2 } \.$$

From Ref[3]

with

$$\begin{array} { r c l } \bar { \Delta } & = & \Delta - ( V + \delta ) / 2 \\ \Delta & = & \delta m _ { 1 3 } ^ { 2 } / ( 2 E ) \\ \delta & = & \delta m _ { 1 2 } ^ { 2 } / ( 2 E ) \,, \end{array}$$

where the neutrino mass differences are δm 2 12 = 7 6 . x 10 -5 ( eV ) 2 , δm 2 13 = 2 4 . x 10 -3 ( eV ) 2 , sin 2 θ = s 12 c 12 δ ω , a = s 13 (∆ -s 2 12 δ ), and E is the neutrino energy. Note that t → L , where L is the baseline, for v ν /similarequal c . The neutrino-matter potential V = 1 13 . × 10 -13 eV.

The main quantity of interest in our present work, I α ∗ , is

$$I _ { \alpha ^ { * } } \ = \ \int _ { 0 } ^ { t } d t ^ { \prime } \alpha ^ { * } ( t ^ { \prime } ) e ^ { - i \bar { \Delta } t ^ { \prime } }$$

with α t ( ) = cos ωt ( ) -icos 2 θsin ωt ( ), ω = δ 2 + V 2 -2 δV cos (2 θ 12 ) /

In Ref.[3], as in Ref.[6], one used δ, ω /lessmuch

√ 2. ∆ to obtain

$$R e [ I _ { \alpha ^ { * } } ] \ & \simeq \ s i n \bar { \Delta } L / \bar { \Delta } \\ I m [ I _ { \alpha ^ { * } } ] \ & \simeq \ ( 1 - c o s \bar { \Delta } L ) / \bar { \Delta } \.$$

On the other hand if one uses Eq(5) to evaluate I α ∗ , one obtains

$$\begin{array} { c c c } \quad \cdot \cdot \\ R e [ I _ { \alpha ^ { * } } ] & = & [ ( \omega - \bar { \Delta } \cos 2 \theta ) \cos \bar { \Delta } L s i n \omega L \\ & & - ( \bar { \Delta } - \omega \cos 2 \theta ) s i n \bar { \Delta } L \cos \omega L ] / ( \omega ^ { 2 } - \bar { \Delta } ^ { 2 } ) \\ I m [ I _ { \alpha ^ { * } } ] & = & [ \bar { \Delta } + \omega \cos 2 \theta - ( \bar { \Delta } + \omega \cos 2 \theta ) \cos \bar { \Delta } L \cos \omega L \\ & & - ( \omega + \bar { \Delta } \cos 2 \theta ) s i n \bar { \Delta } L s i n \omega L ] / ( \omega ^ { 2 } - \bar { \Delta } ^ { 2 } ). \end{array}$$

From this we obtain an improved P ( ν µ → ν e )= (Re[ S 12 ]) 2 +(Im[ S 12 ]) 2 .

$$R e [ S 1 2 ] \ & = \ s _ { 2 3 } a [ \cos ( \bar { \Delta } L ) I m [ I _ { \alpha ^ { * } } ] - s i n ( \bar { \Delta } L ) R e [ I _ { \alpha ^ { * } } ] ] \\ I m [ S 1 2 ] \ & = \ - c _ { 2 3 } s i n 2 \theta s i n \omega L - s _ { 2 3 } a [ \cos ( \bar { \Delta } L ) R e [ I _ { \alpha ^ { * } } ] \\ & \quad + s i n ( \bar { \Delta } L ) I m [ I _ { \alpha ^ { * } } ]$$

Figure 1: P ( ν µ → ν e ) for MINOS(L=735 km), MiniBooNE(L=500m), JHF-Kamioka(L=295 km), and CHOOZ(L=1.03 km) using the improved 3 × 3 mixing matrix. Solid curve for precise I α ∗ ( E,L ) and dashed curve for approximate I α ∗ ( E,L ). s 13 =0.19

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0310v3-figure-6e49f3ec-ef72-442a-a309-a0b58496ba7d.jpg)

## 3 Conclusions

As one can see from the figure, for a large baseline (MINOS) or energy (MiniBooNE) the corrections are small. For some baselines and energies, including JHF-Kamioka and CHOOZ, the corrections to the neutrino transition probability P ( ν µ → ν e ) are large. This will modify the predictions of CP and T violation, and could have a large effect on the extraction of parameters via neutrino oscillations, such as s 13 from P (¯ ν e → ¯ ) ν e by Daya Bay[8], Double Chooz[9], and, RENO[10]. This is an important result in neutrino physics, and the precise expression for I α ∗ we have derived should be used for all future studies of neutrino oscillations.

## Acknowledgements

This work was carried out while LSK was a visitor at Los Alamos National Laboratory, Group P25. The author thanks Dr. William Louis for information about recent and future neutrino oscillation experiments, and Dr. Mikkel B. Johnson for many helpful discussions.

## References

- [1] E.M. Henley, M.B. Johnson, and L.S. Kisslinger, arXiv:1102.5106/ hep-ph; Int. J. Mod. Phys. E20 , 2453 (2011)
- [2] E.M. Henley, M.B. Johnson, and L.S. Kisslinger, arXiv:1105.2741/ hep-ph (2011)
- [3] Leonard S. Kisslinger, Ernest M. Henley, and Mikkel B. Johnson, arXiv:1203.6613/hepph; Int. J. Mod. Phys. E 21 , 1250065 (2012)
- [4] J. Arafune and J. Sato, Phys. Rev. D 55 , 1653 (1997)
- [5] M. Freund, Phys. Rev. D 64 , 053003 (2001)
- [6] E. Akhmedov, P. Huber, M. Linder, and T. Ohlsson, Nucl. Phys B 608 , 394 (2001)
- [7] M. Jacobson and T. Ohlsson, Phys. Rev. D 69 , 013003 (2004)
- [8] F.P. An et.al., Daya Bay Collaboration, arXiv:1210.6327; Chin. Phys. C37 , 01101 (2013)
- [9] M. Kuze (Double Chooz Collaboration), arXiv:1109.0074/hep-ex (2011)
- [10] J., K. Ahn, et.al. (RENO Collaboration), arXiv:1003.1391/hep-ex (2010)
Leonard S. Kisslinger
Improved Theory of Neutrino Oscillations in Matter
This is revision of the S-Matrix theory of neutrino oscillations used for
many years. We evaluate the transition probability of a $\mu$ to $e$ neutrino
without an approximation used for many theoretical studies, and find important
differences which could improve the extraction of neutrino parameters from
experimental data in the future.