Skip to main content
2 of 2
Added link to a detailed article to augment the terse statement of problem here

The transition probability $P_{\nu_e \rightarrow \nu_{\mu}}$ is indeed decreasing with neutrino production energy, i.e., the survival probability $P_{\nu_e \rightarrow \nu_{e}}=1-P_{\nu_e \rightarrow \nu_{\mu}}$ is increasing. Why? The first sine term ($\sin^2 2 \theta_m $) in the oscillation equation (your first equation) determines the amplitude of the oscillations. The limit of that term as $V \rightarrow \infty$, i.e., as energy or density increases towards infinity, is zero: $$ \begin{align} \sin 2 \theta_m = \frac{\sin 2 \theta}{ \sqrt{ \left(\Delta V / \Delta m^2 - \cos 2 \theta \right)^2 + \sin^2 2 \theta} } \\ \lim_{V\to \infty} \; \left[ \frac{\sin 2 \theta}{ \sqrt{ \left( V / \Delta m^2 - \cos 2 \theta \right)^2 + \sin^2 2 \theta} } \right] = 0 \end{align} $$

As $\sin^2 2 \theta_m$ descends towards zero asymptotically with increasing neutrino production energy (and/or production electron number density) above MSW resonance any oscillations in the second sine term, $\sin^2 \left( 1.27 \Delta m_m^2 \frac{L}{E} \right)$, of the oscillation equation will accordingly be reduced in amplitude, going to zero in the limit as well.

The second sine term is not monotonically decreasing (it might appear to be so if only a few samples are calculated), but is rather oscillating rapidly, hence the smeared oscillations in the graph.

The equations you gave should not be used for solar neutrino analysis except at production energies much lower than the MSW resonance.

You can obtain proper equations in the PDG 2018 Review of Particle physics section 14, a free download at lbl.gov.

If the above terse statement of the problem is insufficient you can view an article demonstrating how to graph using the posted equations (with Python code) and how they malfunction at higher energies (with a link to a proper equation and an article describing that) at this link: MSW flavor calculation