Earth-Density Stratification and Quantum Gravity Corrections in Long-Baseline Neutrino Oscillation Experiments

This paper presents a unified three-flavor analysis demonstrating that Earth's density stratification and Planck-scale quantum-gravity corrections act as correlated systematic uncertainties in long-baseline neutrino oscillations, where their interplay can significantly bias or even reverse the inferred CP-violating phase, particularly at baselines exceeding 5000 km.

The Gist

Imagine trying to listen to a faint radio station while driving through a city. To hear the music clearly, you need to know exactly how the buildings, hills, and tunnels along your route will distort the signal. If you assume the road is perfectly flat and empty, your navigation app will give you the wrong direction, and you'll miss the station entirely.

This paper is about two specific "distortions" that happen when scientists try to listen to neutrinos—tiny, ghost-like particles that zip through the Earth. Scientists use these particles to measure a fundamental property of the universe called CP violation (a kind of "handedness" or asymmetry in how nature works).

Here is the simple breakdown of what the authors found:

1. The "Flat Road" Mistake (Earth's Density)

For a long time, scientists calculated how neutrinos travel through the Earth by assuming the planet is like a giant, uniform ball of rock with the same density everywhere. It's like assuming the entire drive from New York to London is on a flat, empty highway.

• The Reality: The Earth is actually layered like a giant cake. It has a thin crust, a thick mantle, and a super-dense core.
• The Problem: When neutrinos travel short distances (like 3,000 km), this "flat road" assumption works fine. But when they travel very long distances (over 5,000 km), they dive deep into the dense lower mantle and even the core.
• The Result: The authors found that if you keep using the "flat road" math for these long trips, your calculation of the neutrino's "handedness" (CP violation) goes wildly wrong.
  • At 7,000 km, the error is about 18 degrees.
  • At 12,000 km (going almost all the way through the Earth), the error is 172 degrees. This is so bad that it flips the answer completely: if the universe is "left-handed," your math says it's "right-handed."

2. The "Quantum Gravity" Whisper

The second distortion comes from something much more exotic: Quantum Gravity. This is the theory that gravity works differently at the tiniest scales (the Planck scale).

• The Idea: The paper suggests that gravity might slightly tweak the "mass" of the neutrinos as they travel, similar to how a slight breeze might change the pitch of a flute note.
• The Effect: This tweak changes the neutrino's internal rhythm (specifically the difference between two of their mass states). It's a tiny effect, but in the precise world of neutrino physics, even a whisper can be heard.

3. The "Magic Cancellation" (The Big Discovery)

This is the most exciting part of the paper. Usually, when you have two sources of error (the "flat road" mistake and the "quantum gravity" whisper), they just add up to make things worse.

However, the authors discovered a strange degeneracy (a fancy word for a confusing overlap) at a specific distance (around 7,000 km):

• The error from the Earth's layers and the error from quantum gravity can cancel each other out.
• Imagine you are trying to balance a scale. One side has a heavy rock (Earth density error), and the other has a small weight (Quantum gravity error). Usually, the scale tips. But in this specific scenario, the small weight is placed in just the right spot to perfectly counterbalance the heavy rock.
• The Catch: This cancellation only happens if the "secret settings" of the neutrinos (called Majorana phases) are just right. If they are, the two big errors hide each other, making the final result look deceptively accurate. If you ignore one of them, you might think your measurement is perfect when it's actually wrong.

Why This Matters for Future Experiments

The paper focuses on future experiments like DUNE (Deep Underground Neutrino Experiment).

• The Warning: If scientists plan to send neutrinos through the Earth for very long distances (like 7,000 km or more) to get ultra-precise measurements, they cannot use the old "flat Earth" math.
• The Solution: They must use a detailed, 3D map of the Earth's density (called the PREM model) and also account for the tiny quantum gravity effects.
• The Risk: If they don't, they won't just be slightly off; they could completely misinterpret the fundamental laws of the universe, thinking they found a signal that isn't there, or missing one that is.

In short: The Earth is not a uniform ball, and gravity might whisper to neutrinos. If you ignore both, your map of the universe will be wrong. If you ignore just one, you might get tricked by a "magic cancellation" that hides the truth. To see the universe clearly, you need to account for both the layers of the Earth and the whispers of quantum gravity simultaneously.

Technical Summary

Technical Summary: Earth-Density Stratification and Quantum Gravity Corrections in Long-Baseline Neutrino Oscillation Experiments

Problem Statement
Long-baseline (LBL) neutrino oscillation experiments, such as T2K, NOvA, and the upcoming Deep Underground Neutrino Experiment (DUNE), aim to measure the CP-violating phase $\delta_{CP}$ with high precision. As these experiments target sub-degree precision, two distinct sources of systematic uncertainty become critical:

1. Earth Matter-Density Stratification: The standard analysis often employs a constant-density approximation for the Earth's matter potential (MSW effect). However, for baselines exceeding $\sim5000$ km, neutrinos traverse the lower mantle and outer core, where electron densities deviate significantly from path-averaged values. Previous studies treated this in isolation, but its impact on very-long-baseline (VLB) proposals remains under-explored as a fundamental systematic.
2. Planck-Scale Quantum Gravity Corrections: Effective Field Theory (EFT) predicts that non-renormalizable gravitational interactions at the Planck scale ($M_{pl} \sim 1.2 \times 10^{19}$ GeV) introduce corrections to the neutrino sector via the unique dimension-5 Weinberg operator. These corrections shift neutrino mass-squared differences, particularly in degenerate mass spectra.

While these effects have been studied independently, no prior analysis has treated them simultaneously. The central problem addressed is whether these two distinct physical mechanisms—spatially varying matter density and Planck-scale mass shifts—can interfere, creating a "degeneracy regime" where they mimic or cancel each other, thereby complicating the extraction of $\delta_{CP}$.

Methodology
The authors present a unified theoretical framework combining three components into a single evolution equation:

• Unified Hamiltonian: The study constructs a total Hamiltonian $H_f(x)$ that simultaneously incorporates the vacuum oscillation term (via the PMNS matrix), the spatially varying MSW matter potential $V_f(x)$ derived from the Preliminary Reference Earth Model (PREM), and the Planck-scale gravitational mass operator.
• PREM Implementation: Instead of a constant density, the authors utilize a spatially resolved four-shell model of the Earth (crust, upper mantle, lower mantle, outer/inner core) with densities ranging from $\sim2.9$ to $\sim13.0$ g cm$^{-3}$. The neutrino trajectory is discretized into 300–400 steps to numerically integrate the evolution equation $i \frac{d}{dx}|\nu(x)\rangle = H_f(x)|\nu(x)\rangle$ using matrix exponentiation.
• Quantum Gravity Perturbation: The Planck-scale effect is modeled as a dimension-5 operator generating a universal mass shift $\mu \approx 2.5 \times 10^{-6}$ eV. Using first-order perturbation theory, the authors calculate the shifts in mass-squared differences ($\Delta'_{ij}$) and mixing angles for a degenerate mass spectrum ($M \sim 2$ eV), explicitly tracking the dependence on Majorana phases ($\alpha_1, \alpha_2$).
• Statistical Analysis: The authors generate "true" event rates using the full PREM profile and fit them using a constant-density approximation to quantify the bias $|\Delta \delta_{CP}|$. They further perform $\chi^2$ analyses to examine the interplay between the two systematics, defining a "degeneracy fraction" to measure partial cancellation.

Key Contributions

1. Unified Analysis: This is the first work to treat Earth density stratification and Planck-scale corrections as correlated systematic uncertainties within a single three-flavour framework.
2. Identification of a New Degeneracy: The paper identifies a specific regime (around $L \approx 7000$ km) where Planck-scale perturbations to $\Delta'_{21}$ can partially cancel or mimic the biases induced by PREM stratification. This degeneracy depends on the unknown Majorana phases.
3. Quantification of Constant-Density Failure: The study rigorously demonstrates that the constant-density approximation, while adequate for $L \lesssim 5000$ km, introduces catastrophic errors for longer baselines, including a complete sign reversal of the reconstructed CP violation at $L = 12000$ km.

Results

• PREM Bias: The bias in $\delta_{CP}$ reconstruction due to the constant-density approximation is negligible ($<0.3^\circ$) for $L \le 5000$ km. However, it grows sharply to $17.8^\circ$ (Normal Ordering) at $L = 7000$ km and reaches $172.2^\circ$ at $L = 12000$ km. The latter represents a total sign reversal of the CP violation signal. These biases exceed the projected $1\sigma$ precision target of DUNE Phase II ($5\text{--}10^\circ$) by factors of 3.6 and 34.4, respectively.
• Planck-Scale Shifts: For a degenerate mass spectrum ($M \sim 2$ eV), the Planck-scale correction shifts the solar mass-squared difference $\Delta'_{21}$ by $(1.0 \pm 0.5) \times 10^{-5}$ eV$^2$, spanning the full experimental $1\sigma$ band. The atmospheric splitting $\Delta'_{31}$ remains effectively unchanged.
• Degeneracy Analysis: At $L = 7000$ km, the interplay between the two effects is non-trivial. When the Majorana phase $\alpha_1 \approx 90^\circ$, the combined bias is reduced by approximately 30% compared to the quadrature sum of the individual biases (degeneracy fraction $f \approx 0.69$). This indicates that marginalizing over only one systematic leads to an underestimation of the total uncertainty on $\delta_{CP}$.
• Mass Ordering: The results hold for both Normal (NH) and Inverted (IH) mass orderings, though the IH bias at $L=7000$ km is slightly lower ($14.2^\circ$) due to the suppression of antineutrino cross-sections.

Significance and Claims
The paper claims that for next-generation LBL experiments targeting sub-degree precision on $\delta_{CP}$, these systematics are not merely "conservative refinements" but fundamental obstacles.

• Necessity of Unified Treatment: Analyses that marginalize over only one of these effects (e.g., ignoring Planck-scale corrections while accounting for PREM, or vice versa) will misestimate confidence intervals. The authors argue that both must be propagated as correlated systematics in global fits.
• Implications for DUNE: The findings necessitate a replacement of the constant-density matter profiles currently used in DUNE simulation frameworks (like GLoBES) with spatially resolved PREM implementations for any baseline beyond $\sim5000$ km.
• Model Independence: The Planck-scale correction is presented as a universal, model-independent prediction derived from the unique gauge-invariant dimension-5 operator in the Standard Model EFT, making its inclusion in precision oscillation analyses a theoretical imperative rather than a speculative addition.

In conclusion, the authors assert that the interplay between Earth's internal structure and quantum gravity effects creates a new class of degeneracy that must be resolved to achieve the physics goals of future very-long-baseline neutrino experiments.