Parametric Resonance: When Earth's Core Tunes a Neutrino Oscillation

The MSW effect — the matter-enhanced oscillation that resolves the solar neutrino problem — is the textbook example of how the propagation medium affects neutrino mixing. As a neutrino moves through matter, its effective oscillation Hamiltonian picks up a contribution from coherent forward scattering, and the mixing angle in matter differs from the mixing angle in vacuum. At a critical density, the in-matter mixing reaches its maximum value; oscillations near this MSW resonance are dramatically enhanced.

This is a single-medium phenomenon: a homogeneous medium at the right density does the job. The standard textbook treatment of the MSW effect assumes constant or slowly varying density, and the resonance is encountered exactly once during the neutrino’s trip.

But there is a different kind of resonance — parametric resonance — that can occur when a neutrino crosses a medium whose density varies in a particular way: as a step function, or periodically, with a length scale that matches the local oscillation length. The parametric effect builds up coherently across multiple zones in a way that single-zone MSW does not. Its existence has been understood since the 1980s, but its observability in real atmospheric neutrinos turns on the Earth providing a useful step-function density profile in the form of the core-mantle boundary.

This post is about how parametric resonance works, why the Earth interior is one of the cleanest places to look for it, and how the next generation of atmospheric-neutrino experiments will use it to disentangle the mass ordering.

The Earth as a two-zone medium

The Earth is not homogeneous. Its density profile is well-known from seismology and is parameterised in the Preliminary Reference Earth Model (PREM). For neutrino-oscillation purposes the relevant features are:

• The mantle, extending from the surface to 2,890 km depth, with an average density of about 4.5 g/cm³.
• The outer core, from 2,890 km to 5,150 km depth, with an average density of about 11.0 g/cm³.
• The inner core, from 5,150 km to the centre at 6,371 km, with an average density of about 13.0 g/cm³.

The transition from mantle to outer core at 2,890 km is the largest matter-density discontinuity inside the Earth — a factor of more than two increase across a 5-km transition zone. Neutrinos crossing the Earth at zenith angles between about $33°$ and $0°$ (where $0°$ is straight up through Earth’s centre) pass through this boundary twice: once on the way in, once on the way out.

The neutrino’s effective oscillation potential, $V_{\mathrm{CC}}=\sqrt{2}{G}_F{n}_e$, scales linearly with the electron density. Inside the mantle the potential is approximately $1.0\times 10^{-13}$ eV at the relevant electron density; inside the core it jumps to about $2.5\times 10^{-13}$ eV. The neutrino therefore experiences a step-function jump in its effective Hamiltonian at each core-mantle boundary crossing.

How parametric resonance works

In a homogeneous medium, the standard MSW Hamiltonian for two-flavor oscillation is

$H=\frac{1}{4E}\left(\begin{array}{cc}-\Delta m^2\cos 2\theta +2VE& \Delta m^2\sin 2\theta \\ \Delta m^2\sin 2\theta & \Delta m^2\cos 2\theta -2VE\end{array}\right)$

where $V$ is the constant matter potential. In the two-zone case where the neutrino passes through alternating regions of different potentials $V_1$ and $V_2$, the evolution is no longer governed by a single Hamiltonian but by a product of two evolution operators:

$U_{\mathrm{total}}=U_1{U}_2{U}_1$

(for a three-segment path mantle-core-mantle, with $U_1$ being the mantle evolution operator and $U_2$ the core one).

The remarkable feature of this product evolution is that under specific conditions on the segment lengths and the potentials, the total evolution operator can produce near-maximal mixing even when each individual operator $U_1$ and $U_2$ produces only modest mixing. The resonance condition is approximately

${\varphi }_1+{\varphi }_2\approx \pi$

where ${\varphi }_1$ and ${\varphi }_2$ are the oscillation phases accumulated in each zone. When this condition is met, the neutrino emerges with a substantially different flavor content than either single zone would predict.

This is the essence of parametric resonance: a multi-segment medium can produce coherent enhancement that no single segment achieves.

Where it kicks in for Earth-crossing neutrinos

The parametric effect for Earth-crossing atmospheric neutrinos peaks at specific energies and zenith angles. The relevant phenomenon was first analysed in detail by Liu and Smirnov in the 1990s and refined by Petcov and collaborators in the 2000s.

For the atmospheric mass splitting $\Delta m_{31}^2$ and matter potentials at the mantle and outer-core densities, the resonance condition is met for upward-going atmospheric neutrinos at energies of approximately:

$E_{\mathrm{parametric}}\approx 4\text{–}10\mathrm{GeV}$

and zenith angles in the range where the neutrino path passes through the outer core but not deep into the inner core — approximately $\cos {\theta }_z\in (-1,-0.84)$, corresponding to zenith angles between 145° and 180° from the local vertical at the detector.

In this regime, the standard MSW-driven ${\nu }_{\mu }\to {\nu }_e$ appearance probability is approximately 5-15% larger (for one mass ordering) or smaller (for the other) than the homogeneous-mantle calculation would predict. The enhancement is small in absolute terms but is concentrated in a specific energy-and-zenith region of phase space, making it identifiable in event distributions.

The mass-ordering connection

The reason parametric resonance matters for current neutrino physics is its mass-ordering sensitivity. The standard matter-potential MSW resonance in the mantle alone enhances neutrino oscillation for one ordering and antineutrino oscillation for the other; the parametric effect amplifies this asymmetry by a factor of 1.5-2 in the relevant energy-zenith region.

For atmospheric muon neutrinos, the observable signature is:

• For normal ordering ($\Delta m_{31}^2>0$): enhanced ${\nu }_{\mu }\to {\nu }_e$ appearance at 4-10 GeV through the core, with no corresponding ${\bar{\nu }}_{\mu }$ enhancement.
• For inverted ordering ($\Delta m_{31}^2<0$): the reverse — enhanced ${\bar{\nu }}_{\mu }\to {\bar{\nu }}_e$ at the same energies, no enhancement for neutrinos.

A detector that can distinguish atmospheric ${\nu }_{\mu }$ from ${\bar{\nu }}_{\mu }$ events as a function of zenith angle and energy can in principle determine the mass ordering by measuring this pattern. The challenge is that conventional atmospheric-neutrino detectors (Super-K, IceCube) cannot distinguish neutrinos from antineutrinos event-by-event — they see the same Cherenkov ring or cascade in either case.

Two strategies are used in practice:

Statistical separation through the kinematic distribution of events: neutrinos and antineutrinos have slightly different angular and energy distributions in their interactions, allowing a statistical disentangling on the population level even without event-by-event identification.

Magnetic separation in detectors with embedded magnetic fields (such as the never-built INO experiment or hypothetical magnetised LArTPCs) provides direct event-by-event identification but is technologically demanding.

KM3NeT-ORCA and IceCube Upgrade analyses

The two operating experiments that exploit parametric resonance for mass-ordering measurement are KM3NeT-ORCA in the Mediterranean and the IceCube Upgrade at the South Pole, plus the original IceCube DeepCore subarray.

ORCA’s instrumented volume of about 7 megaton is densely packed for GeV-scale event reconstruction, with sensors spaced approximately 9 metres vertically and 20 metres horizontally — much finer than IceCube’s 17-metre by 125-metre spacing. The dense spacing allows event reconstruction down to a few GeV, accessing the parametric-resonance energy range from below.

ORCA’s predicted mass-ordering sensitivity, including the parametric-resonance enhancement, reaches 3-sigma within 5 years of full operation. The IceCube Upgrade’s dense in-fill provides comparable performance at the South Pole.

Super-Kamiokande’s recent atmospheric-oscillation analyses also include the parametric effect explicitly, with the predicted oscillation probabilities computed using the full PREM density profile rather than the historical constant-density approximation. The parametric contribution shifts the central value of the extracted ${\theta }_{23}$ by a small but non-negligible amount, particularly in the upper octant.

The PREM profile and its uncertainties

The PREM density model used for these calculations is based on seismic-wave-propagation data and has been validated by decades of geophysical research. Its quoted uncertainties on the mantle and core densities are around 1-2%, which translates to a similar-level uncertainty on the matter potential and consequently on the parametric-resonance phenomenon.

The 1-2% PREM uncertainty is much smaller than the current oscillation-parameter precision, so it does not significantly limit the mass-ordering analyses. But for future precision atmospheric oscillation measurements at the per-cent level, the Earth-interior density profile becomes a relevant systematic, and dedicated seismic and neutrino-tomography measurements may eventually be needed to refine it further.

Why this is a useful pedagogical example

Parametric resonance is one of the cleaner pedagogical demonstrations that neutrino-oscillation phenomenology contains more than just the textbook MSW effect. The two-zone structure of the Earth — a feature of basic geology — produces a quantum-mechanical phenomenon that is observable in actual data and that interacts with the mass-ordering measurement in a controlled way.

The fact that essentially every modern atmospheric-neutrino analysis includes parametric-resonance effects, computed from a realistic Earth-density model, illustrates how oscillation physics has matured from the simple textbook formulas of the 1990s into a precise predictive framework that takes proper account of the geometric and material structure of the medium the neutrinos traverse.

Summary

Parametric resonance is a phenomenon that arises in multi-zone media where the period or step-scale of the medium matches the local oscillation length. For atmospheric neutrinos crossing the Earth at zenith angles below 145 degrees, the three-zone mantle-core-mantle structure produces parametric resonance at energies of 4-10 GeV, enhancing ${\nu }_{\mu }\to {\nu }_e$ appearance probabilities by 5-15% relative to a homogeneous-medium calculation. The enhancement appears with opposite sign for the two mass orderings, providing a mass-ordering signature that atmospheric neutrino experiments — KM3NeT-ORCA, the IceCube Upgrade, Super-Kamiokande’s continued atmospheric analyses — exploit through statistical or event-by-event neutrino-antineutrino discrimination. The PREM Earth-density model provides the input for these calculations at the 1-2% level, sufficient for current oscillation precision. Parametric resonance is one of the cleaner examples of how the textbook MSW formula extends into more complex multi-medium geometries, and one of the reasons that core-crossing atmospheric neutrinos provide a complementary probe of the mass ordering to the reactor and accelerator experiments discussed elsewhere on this blog.