I'll try to explain what causes neutrino oscillations.
Corresponding to each flavor of charged lepton is a flavor or weak-interaction state. It's the weak interaction that turns electrons into electron neutrinos and back again, and likewise for muons and tau leptons and their neutrinos.
But each flavor of neutrinos is a mixture of mass states, and each of them oscillates with a different frequency, because each one has a different mass. These oscillations turn a pure flavor into an oscillating mixture of flavors.
I'll try to explain it in the two-flavor case, because the math is easier than for the three-flavor case. I'll imagine only electron and muon neutrinos. They mix as follows:
(e-nu) = cos(θ) * (m1-nu) + sin(θ) * (m2-nu)
(μ-nu) = - sin(θ) * (m1-nu) + cos(θ) * (m2-nu)
The mass states can be found from the flavors as follows:
(m1-nu) = cos(θ) * (e-nu) - sin(θ) * (μ-nu)
(m2-nu) = sin(θ) * (e-nu) + cos(θ) * (μ-nu)
Now let them oscillate over time t, where (m1-nu) has energy E1 and (mu2-nu) has energy E2. The new mass states become
(m1-nu)' = (m1-nu) * exp(-i*E1*t)
(m2-nu)' = (m2-nu) * exp(-i*E2*t)
Then,
(e-nu)' = cos(θ) * (m1-nu)' + sin(θ) * (m2-nu)'
(μ-nu)' = - sin(θ) * (m1-nu)' + cos(θ) * (m2-nu)'
(e-nu)' = cos(θ) * (cos(θ) * (e-nu) - sin(θ) * (μ-nu))*exp(-i*E1*t) + sin(θ) * (sin(θ) * (e-nu) + cos(θ) * (μ-nu))*exp(-i*E2*t)
(μ-nu)' = - sin(θ) * (cos(θ) * (e-nu) - sin(θ) * (μ-nu))*exp(-i*E1*t) + cos(θ) * (sin(θ) * (e-nu) + cos(θ) * (μ-nu))*exp(-i*E2*t)
(e-nu)' = (cos
2(θ)*exp(-i*E1*t) + sin
2(θ)*exp(-i*E2*t)) * (e-nu) + (cos(θ)*sin(θ)*(exp(-i*E2*t) - exp(-i*E1*t)) * (μ-nu)
(μ-nu)' = (sin
2(θ)*exp(-i*E1*t) + cos
2(θ)*exp(-i*E2*t)) * (μ-nu) - (cos(θ)*sin(θ)*(exp(-i*E2*t) - exp(-i*E1*t)) * (e-nu)
Since the observed intensity is proportional to the absolute square of the wavefunction, we must take it.
|e-nu| -> |μ-nu|' = |μ-nu| -> |e-nu|' = T
|e-nu| -> |e-nu|' = |μ-nu| -> |μ-nu|' = 1 - T
where
T = 2*cos
2(θ)*sin
2(θ)*(1 - cos((E2-E1)*t))
Since neutrino masses are much smaller than most reaction energies, we get E2 - E1 = (m2
2 - m1
2)/(2*E) to lowest order.
Neutrino oscillation goes into a lot of technical detail for the three-flavor case.