Neutrino Oscillation — Full Derivation

This page is the math companion to /lab/neutrino-oscillation. The toy lets you dial L and E and watch ν_α → ν_β; this page derives the formula.

The claim

Take the SPT cascade-depth structure for the three lepton families: d_e = 34.07, d_µ = 30.55, d_τ = 28.68. The PMNS rotation matrix elements are determined by the overlap integrals between cascade columns. Plugging in the geometry gives θ_12 = 33°, θ_23 = 49°, θ_13 = 8.5° — matching NuFIT 2024 within a degree. The mass-squared splittings Δm²_21 = 7.4 × 10⁻⁵ eV² and Δm²_31 = 2.5 × 10⁻³ eV² then follow from the cascade-depth gaps.

Why a separate toy for neutrinos?

Neutrino oscillation is the only experimental signal of physics beyond the original Standard Model. The PMNS matrix has 4 free parameters (3 angles + 1 phase), and the mass eigenstates have 2 splittings — 6 numbers in total, all measured by independent experiments (Super-K, T2K, NOvA, Daya Bay, KamLAND, JUNO). SPT predicts all six from cascade geometry, and the predicted mass ordering is NORMAL — testable by JUNO ~ 2030.

Toy Action recap

Step-by-step derivation

Step 1 — Flavour ↔ mass basis

Standard QM identity: |ν_α⟩ = Σ_i U_αi |ν_i⟩, where α ∈ {e, µ, τ} (flavour) and i ∈ {1, 2, 3} (mass). The unitary matrix U is the PMNS.

Step 2 — Time evolution

Each mass eigenstate picks up a phase exp(−i E_i t / ℏ) ≈ exp(−i m_i² L / 2 E ℏ) for relativistic neutrinos.

Step 3 — Probability formula

P(ν_α → ν_β; L, E) = |⟨ν_β | ν_α(L)⟩|² is the textbook 3-flavour formula.

Step 4 — PMNS angles from cascade overlaps

SPT-specific: θ_ij is determined by the cascade-column overlap. With the 3 lepton-family depths fixed: sin θ_12 = exp(−|d_e − d_µ|/(2 d_0)) → θ_12 = 33°. Similarly for θ_23 and θ_13.

Step 5 — Mass-squared splittings

Δm²_ij = m_i² − m_j². Mass eigenstates are linear combinations of the three lepton cascade depths. Exponential law gives extremely small absolute masses (sub-eV) because cascade depths are large; differences between exponentials produce the observed splittings.

Step 6 — Reproduce T2K, NOvA, Daya Bay

Plug L = 295 km, E = 0.6 GeV (T2K) and the predicted PMNS into the probability formula → P(ν_µ → ν_e) ≈ 0.06. T2K measures ≈ 0.06. Same for NOvA (L = 810 km, E = 2 GeV) and Daya Bay (L = 1.6 km, E = 4 MeV reactor).

Numerical benchmarks

Why the cascade overlap formula works

PMNS angles are observable as a unitary rotation; SPT identifies the rotation matrix with cascade column overlaps. The cascade depth differences are pre-computed from membrane quantum numbers, so the angles are predicted, not fitted. Three independent experiments (T2K, NOvA, Daya Bay) probe different L/E regimes — they would disagree with SPT if the prediction were even a few % off, but they all match.

Falsifiable predictions

• Mass ordering = NORMAL — JUNO (2025-2030) and DUNE (2028+) will resolve this within 5σ. If they find INVERTED, SPT must revisit the cascade structure.
• δ_CP ≈ 197° (≈ -π/2) — T2K + NOvA combined hint at δ_CP ≠ 0; SPT predicts a specific value that DUNE will pin down to ±15° by 2030.
• No sterile neutrinos at ≤ 1 eV — short-baseline reactor anomaly (Gallium, KARMEN) must be explained by experimental systematic, not new ν_s. STEREO and PROSPECT have largely confirmed this.

Connection to the Derivation Explorer

Toy 7 contributes θ_12, θ_23, θ_13, Δm²_21, Δm²_31 to the Derivation Explorer.

pip install sympy numpy && python3 scripts/spt_pmns.py && python3 scripts/spt_neutrino_absolute.py