O.L.G.Peres, A.Yu.Smirnov Nucl.Phys. B680(2004) 479 “Atmospheric neutrinos:LMA oscillations, Ue3 induced interference and CP-violation”より
$\nu_f=\pmatrix{\nu_e\\\nu_\mu\\\nu_\tau}$として、その時間発展は、 で与えられる。ここで、$U$はMNS行列$U_{23}D_{\delta_{CP}}U_{13}U_{12}$で、$D_{\delta_{CP}}=\pmatrix{1&0&0\\0&1&0\\0&0&e^{i\delta_{CP}}}$。$V$が、物質によるポテンシャルの項で$V=\pmatrix{V&0&0\\0&1&0\\0&0&1}$である($V=\sqrt{2}G_FN_e$)。
このHamiltonianを対角化して、固有行列$\tilde\nu=\pmatrix{\tilde\nu_e\\ \tilde\nu_\mu\\ \tilde\nu_\tau}$を求めれば良い。
まず、$\nu_f=U_{23}D_{\delta_{CP}}U_{13}\nu'$で与えられる$\nu'$を考える。すると、$\nu'$に対するHamilitonianは、 となる。この行列にさらに$(\nu_e'-\nu_\tau')$の回転を加えて$H_{13}',H_{31}'$を取り除く。