Skip to content

Neutrino Physics

Neutrinos are produced and detected in flavour eigenstates (νe,νμ,ντ)(\nu_e, \nu_\mu, \nu_\tau)But propagate As mass eigenstates (ν1,ν2,ν3)(\nu_1, \nu_2, \nu_3) related by the PMNS mixing matrix UU:

να=iUαiνi|\nu_\alpha\rangle = \sum_i U_{\alpha i}^* |\nu_i\rangle

As a neutrino of flavour α\alpha propagates, the mass eigenstates acquire different phases: exp(imi2L/(2E))\exp(-im_i^2 L/(2E))Leading to oscillations.

Two-flavour oscillation probability:

P(νανβ)=sin2(2θ)sin2(Δm2L4E)P(\nu_\alpha \to \nu_\beta) = \sin^2(2\theta)\sin^2\left(\frac{\Delta m^2 L}{4E}\right)

Where Δm2=m22m12\Delta m^2 = m_2^2 - m_1^2, θ\theta is the mixing angle, LL is the distance, and EE is the Energy.

Evidence: The Solar Neutrino Problem (deficit of νe\nu_e from the Sun, resolved by νeνμ,ντ\nu_e \to \nu_\mu, \nu_\tau oscillations) and atmospheric neutrino oscillations (Super-Kamiokande, 1998).

Neutrino oscillations imply that neutrinos have mass, but the masses are extremely small: mν<0.12\sum m_\nu \lt 0.12 eV (Planck 2018).

In the Standard Model, neutrinos are massless. Their masses require physics beyond the Standard Model, most commonly via the seesaw mechanism:

mνmD2Mm_\nu \sim \frac{m_D^2}{M}

Where mDm_D is a Dirac mass and MmDM \gg m_D is the mass of a heavy right-handed neutrino.

Example 8.1: Atmospheric neutrino oscillation calculation

Atmospheric neutrinos are produced when cosmic rays strike the upper atmosphere, creating Pions that decay: π+μ++νμ\pi^+ \to \mu^+ + \nu_\muFollowed by μ+e++νˉμ+νe\mu^+ \to e^+ + \bar{\nu}_\mu + \nu_e.

Super-Kamiokande (1998) observed that upward-going muon neutrinos (travelling through the Earth, L104L \sim 10^4 km) were significantly depleted relative to downward-going ones (L10L \sim 10 km), while electron neutrinos showed no such deficit.

Using the two-flavour formula with the atmospheric parameters Δm3222.5×103\Delta m^2_{32} \approx 2.5 \times 10^{-3} eV2^2 and sin2(2θ23)1\sin^2(2\theta_{23}) \approx 1 (maximal mixing):

For upward-going νμ\nu_\mu with E=1E = 1 GeV and L=10000L = 10\,000 km:

\frac{\Delta m^2 L}{4E} = \frac{2.5 \times 10^{-3}\;\mathrm{eV}^2 \times 10^4\;\mathrm{km}{4 \times 1\;\mathrm{GeV}}}

Converting to natural units (c=1.973×107\hbar c = 1.973 \times 10^{-7} eV\cdotM): L=107L = 10^7 m, so L/E=107/109=102L/E = 10^7 / 10^9 = 10^{-2} eV1^{-1}.

Δm2L4E=2.5×103×1024=6.25×106  eV2eV1\frac{\Delta m^2 L}{4E} = \frac{2.5 \times 10^{-3} \times 10^{-2}}{4} = 6.25 \times 10^{-6}\;\mathrm{eV}^2\cdot\mathrm{eV}^{-1}

Wait --- we need to be more careful with units. Using the practical formula:

Δm2[eV2]L[km]4E[GeV]=2.5×103×1044×1=254=6.25  rad\frac{\Delta m^2 [\mathrm{eV}^2] \cdot L [\mathrm{km}]}{4E [\mathrm{GeV}]} = \frac{2.5 \times 10^{-3} \times 10^4}{4 \times 1} = \frac{25}{4} = 6.25\;\mathrm{rad}

P(νμνμ)=1sin2(2θ23)sin2(6.25)=11×sin2(6.25)10.0180.98P(\nu_\mu \to \nu_\mu) = 1 - \sin^2(2\theta_{23})\sin^2(6.25) = 1 - 1 \times \sin^2(6.25) \approx 1 - 0.018 \approx 0.98

Hmm, this gives almost no oscillation. Let me reconsider. Actually:

P(νμντ)=sin2(2θ)sin2(Δm2L4E)=sin2(6.25)0.018P(\nu_\mu \to \nu_\tau) = \sin^2(2\theta)\sin^2\left(\frac{\Delta m^2 L}{4E}\right) = \sin^2(6.25) \approx 0.018

This seems small. But at E=0.5E = 0.5 GeV:

Δm2L4E=252=12.5  rad\frac{\Delta m^2 L}{4E} = \frac{25}{2} = 12.5\;\mathrm{rad}

sin2(12.5)sin2(0.35)0.12\sin^2(12.5) \approx \sin^2(0.35) \approx 0.12

And at the first oscillation maximum, L/E=2π/(Δm2)=2π/(2.5×103)2513L/E = 2\pi/(\Delta m^2) = 2\pi/(2.5 \times 10^{-3}) \approx 2513 km/GeV. For E=1E = 1 GeV, Losc2513L_{\mathrm{osc} \approx 2513} km, which is comparable to the Earth”s diameter (12700\sim 12\,700 km). The observed deficit is an average over many oscillations and energies, Giving roughly P1/2\langle P\rangle \approx 1/2 for maximal mixing, consistent with the Super-Kamiokande observation of approximately half the expected upward-going νμ\nu_\mu flux.

ParameterValue (best fit)ExperimentRole
Δm212\Delta m^2_{21}7.5×1057.5 \times 10^{-5} eV2^2Solar (SNO, Borexino)Drives solar νeνμ,τ\nu_e \to \nu_{\mu,\tau}
Δm322\Delta m^2_{32}2.5×1032.5 \times 10^{-3} eV2^2Atmospheric (Super-K)Drives νμντ\nu_\mu \to \nu_\tau oscillations
sin2(2θ12)\sin^2(2\theta_{12})0.86Solar (SNO)Solar mixing angle
sin2(2θ23)\sin^2(2\theta_{23})1.0 (maximal)Atmospheric (Super-K)Atmospheric mixing angle
sin2(2θ13)\sin^2(2\theta_{13})0.092Reactor (Daya Bay, RENO, Double Chooz)Non-zero, enables CP violation
  • Confusing flavour and mass eigenstates. Neutrinos are produced and detected as flavour eigenstates but propagate as mass eigenstates. Fix: The PMNS matrix UU relates the two bases; oscillations arise from phase differences between mass components.
  • Assuming all oscillations average to zero. While fast oscillations average over energy and baseline, the survival probability for solar νe\nu_e is 0.55\approx 0.55, not 0.50.5, due to the MSW matter effect in the Sun. Fix: Matter effects modify the effective mixing angle at high densities.
  • Forgetting the practical unit conversion. When using the oscillation formula Δm2L/(4E)\Delta m^2 L / (4E) with LL in km and EE in GeV, the result is in radians directly: Δm2[eV2]L[km]/(4E[GeV])\Delta m^2 [\mathrm{eV}^2] L [\mathrm{km}] / (4E [\mathrm{GeV}]). Fix: Use L/EL/E in km/GeV for quick estimates.
  • Thinking neutrinos are massless in the Standard Model. While the SM predicts massless neutrinos, oscillations prove they have mass. Fix: The seesaw mechanism extends the SM with heavy right-handed neutrinos.
  • Solar neutrino spectroscopy: Precise measurement of solar neutrino fluxes (pp, 7^7Be, 8^8B) tests solar models and constrains the MSW effect transition between vacuum and matter-dominated oscillations.
  • Reactor neutrino monitoring: Antineutrino detectors at nuclear reactors (Daya Bay, Double Chooz) measure θ13\theta_{13} and can monitor reactor power and fuel composition for non-proliferation.
  • Supernova neutrinos: Core-collapse supernovae release 99%\sim 99\% of their gravitational binding energy as neutrinos. Detecting these (SN 1987A, 20\sim 20 events) tests models of stellar death and neutron star formation.
  • Neutrino telescopes: IceCube and KM3NeT detect high-energy astrophysical neutrinos from blazars, gamma-ray bursts, and possibly dark matter annihilation.
  • Cosmology: The sum of neutrino masses mν<0.12\sum m_\nu < 0.12 eV (Planck + BAO) affects structure formation; future surveys (Euclid, DESI) will tighten constraints and potentially determine the mass hierarchy.
Neutrino sourceTypical energyBaseline LLOscillation probedKey experiment
Solar0.1-10 MeV1.5×1081.5 \times 10^8 kmΔm212\Delta m^2_{21}SNO, Borexino, Super-K
Atmospheric0.1-100 GeV10-104^4 kmΔm322\Delta m^2_{32}Super-Kamiokande
Reactor1-10 MeV0.1-100 kmθ13\theta_{13}, Δm212\Delta m^2_{21}Daya Bay, RENO, KamLAND
Accelerator0.1-10 GeV100-1000 kmδCP\delta_{\rm CP}T2K, NOν\nuA, DUNE
  • Mass hierarchy: Is the ordering m1<m2<m3m_1 < m_2 < m_3 (normal) or m3<m1<m2m_3 < m_1 < m_2 (inverted)? Future experiments like JUNO and DUNE aim to resolve this via matter effects in oscillation probabilities.
  • CP violation in the lepton sector: The phase δCP\delta_{\rm CP} in the PMNS matrix determines whether neutrinos and antineutrinos oscillate differently. A non-zero δCP\delta_{\rm CP} could help explain the matter-antimatter asymmetry via leptogenesis.
  • Dirac vs. Majorana nature: Are neutrinos their own antiparticles? Neutrinoless double-beta decay (0νββ0\nu\beta\beta) experiments search for this; a positive signal would prove the Majorana nature and fix the absolute mass scale.
  • Absolute mass scale: Oscillations only measure mass-squared differences. KATRIN measures the electron neutrino mass via tritium beta decay, currently constraining mνe<0.8m_{\nu_e} < 0.8 eV.