Earth-Density Effects in Long Baseline Neutrino Experiments

This paper demonstrates that realistic spatial variations in Earth's matter density introduce energy-dependent structures in neutrino oscillation probabilities that cannot be captured by path-averaged approximations, thereby creating degeneracies that obscure CP violation measurements and necessitating the use of spatially resolved density profiles in future long-baseline experiment analyses.

The Gist

Imagine trying to listen to a specific radio station while driving through a city. If the city were perfectly flat and empty, you could predict exactly how the signal would travel. But the Earth isn't empty; it's a layered cake of rock, metal, and magma, with different densities at different depths.

This paper is about how scientists listen to neutrinos—tiny, ghost-like particles that zip through the Earth—and how the "layers" of our planet mess with the signal.

Here is the breakdown of what the authors found, using simple analogies:

1. The Ghosts and the Traffic Jam

Neutrinos come in three "flavors" (like ice cream flavors: electron, muon, and tau). As they travel, they naturally switch flavors, a bit like a chameleon changing colors. Scientists want to measure a specific "twist" in this switching process called CP violation (let's call it the "secret code"). If they get this code wrong, they might think the universe is behaving differently than it actually is.

However, when these neutrinos travel through the Earth, they bump into electrons in the rock. This creates a "traffic jam" (known as the MSW effect) that changes how the neutrinos switch flavors. The denser the rock, the heavier the traffic jam.

2. The "Flat Earth" Mistake

For a long time, scientists studying neutrinos traveling long distances (hundreds or thousands of miles) made a simplifying assumption: they treated the Earth's density as if it were constant.

• The Analogy: Imagine driving from New York to London. To calculate your fuel usage, you assume the road is a perfectly flat, straight line with no hills or valleys.
• The Reality: The Earth is actually a layered cake. The crust is light, the mantle is heavier, and the core is incredibly dense.

The authors asked: Is it okay to pretend the road is flat, or does the "hills and valleys" of the Earth actually change our results?

3. The Short Trip vs. The Long Trip

The team ran simulations to see what happens when neutrinos travel different distances through the Earth.

• The Short Trip (Up to 3,000 miles):
  If the neutrinos travel a shorter distance, the "flat road" assumption works fine. The error in their "secret code" measurement is tiny—less than the width of a hair (less than 0.3 degrees). It's like driving a few miles on a slightly bumpy road; you don't really notice the difference in your fuel calculation.

• The Long Trip (Over 4,000 miles):
  This is where things get messy. As the distance increases, the neutrinos dive deeper into the Earth, hitting the heavy lower mantle and eventually the dense core.

  • The Result: The "flat road" assumption breaks down completely.
  • At 4,300 miles, the error jumps to nearly 18 degrees.
  • At 7,400 miles, the error explodes to 172 degrees.
  • The Analogy: It's like trying to drive across the entire Earth assuming the road is flat. You would end up thinking you are in a completely different country than you actually are. In fact, at the longest distances, the error is so big that scientists might think the "secret code" is the exact opposite of what it really is.

4. Why Does This Happen?

The paper explains that the Earth's layers act like a complex filter. Because the density changes as the neutrino goes deeper, the "traffic jam" changes strength along the way.

• If you pretend the density is constant, you miss these subtle shifts.
• These shifts create a confusion between the "natural" flavor switching and the "Earth-induced" switching.
• The authors found that you can't fix this by just picking an "average" density. It's like trying to average the temperature of a freezer and a furnace; the average doesn't tell you what's actually happening inside either.

5. The Bottom Line

The authors conclude that for future experiments that send neutrinos very far (thousands of miles), we cannot use the simple "constant density" shortcut anymore. It is not a safe simplification; it is a source of fundamental error.

To get the right answer, scientists must use a detailed map of the Earth's layers (called the PREM model), which accounts for every change in density from the crust down to the core. Without this detailed map, our measurements of the universe's fundamental secrets could be completely wrong.

In short: If you want to measure the universe's secrets using neutrinos traveling long distances, you can't pretend the Earth is a uniform block of cheese. You have to respect the layers, or you'll get the recipe wrong.

Technical Summary

Technical Summary: Earth-Density Effects in Long Baseline Neutrino Experiments

Problem Statement
The accurate determination of the CP-violating phase ($\delta_{CP}$) in long-baseline (LBL) neutrino oscillation experiments relies heavily on the precise modeling of Earth matter effects. As neutrinos propagate through the Earth, they undergo the Mikheyev–Smirnov–Wolfenstein (MSW) effect, where coherent forward scattering off electrons creates an effective potential that modifies oscillation probabilities. While current analyses often employ a "constant density" or "path-averaged density" approximation to simplify calculations, this paper investigates whether such simplifications remain valid as experiments enter a precision era. The authors posit that mismodeling the Earth's internal density structure—specifically the transition from the crust and mantle to the core—may introduce systematic biases that mimic or obscure intrinsic CP violation, particularly for baselines exceeding current experimental scales.

Methodology
The authors utilize a full three-flavour neutrino oscillation framework to simulate $\nu_\mu \to \nu_e$ appearance probabilities. The study compares two distinct modeling approaches across baselines ranging from $L = 1000$ km to $L = 12000$ km:

1. Ground Truth: A spatially resolved density profile based on the Preliminary Reference Earth Model (PREM), which accounts for the Earth's stratified layers (crust, upper/lower mantle, outer/inner core).
2. Approximation: A constant density model utilizing the path-length-weighted average density ($\rho_{avg}$) for each specific trajectory.

The simulation employs a matrix exponentiation technique to solve the Schrödinger-like evolution equation in the flavor basis, incorporating the full three-flavour Hamiltonian (vacuum terms derived from the PMNS matrix and mass-squared differences, plus the matter potential). The analysis assumes a normal mass ordering and a true CP phase of $\delta_{CP} = -90^\circ$, a value chosen to maximize leptonic CP asymmetry and align with current experimental hints from T2K and NO$\nu$A.

To quantify the impact, the authors calculate the absolute bias in the reconstructed $\delta_{CP}$ ($|\Delta \delta_{CP}|$) by minimizing a Poisson $\chi^2$ statistic between the "true" PREM-generated event rates and the rates predicted by the constant density model. They also analyze $\Delta \chi^2$ profiles and event rate spectra to examine how density mismodeling distorts the statistical landscape and energy dependence of the oscillation signal.

Key Results
The study reveals a critical threshold in baseline length where the constant density approximation transitions from a negligible simplification to a source of significant systematic error:

• Short-to-Medium Baselines ($L \leq 5000$ km): The constant density approximation introduces a negligible bias of less than $0.3^\circ$ in the reconstructed $\delta_{CP}$. At these distances, neutrinos primarily traverse the upper mantle and crust, where density variations are slow enough that an average value suffices.
• Long Baselines ($L > 5000$ km): Beyond 5000 km, the bias increases sharply. At $L = 7000$ km, the bias reaches $17.8^\circ$. At $L = 12000$ km, where the trajectory penetrates the high-density outer core, the bias explodes to $172.2^\circ$.
• Non-Monotonic Behavior: The bias is not strictly monotonic; it decreases from $17.8^\circ$ at 7000 km to $8.8^\circ$ at 9000 km before rising again. The authors attribute this to oscillation phase degeneracies where the constant density model's $\chi^2$ minimum accidentally coincides with the true value.
• Structural Divergence: At $L = 7000$ km, the constant density model predicts a sharper, more pronounced oscillation peak near 3.0 GeV compared to the broader, suppressed peak predicted by the PREM profile. Conversely, above 4 GeV, the PREM profile yields significantly higher event rates due to the accumulated potential from the denser lower mantle. This divergence exceeds $1\sigma$ statistical uncertainties, indicating the models are statistically distinguishable.
• Fundamental Error: The analysis demonstrates that the bias arises from energy-dependent distortions caused by Earth's internal layering. These distortions create degeneracies between matter-induced effects and intrinsic CP violation that cannot be resolved by marginalizing over a single effective density parameter.

Significance and Claims
The paper concludes that the constant density approximation is not a "conservative simplification" but a source of fundamental systematic error for long-baseline experiments. The authors argue that as experiments push toward higher precision, the Earth's density profile must be treated as a spatially resolved systematic uncertainty rather than a simple numerical constant.

While the authors note that no currently proposed accelerator-based LBL experiment exceeds a baseline of approximately 3000 km (where the approximation remains valid), they emphasize the relevance of their findings for:

1. Future Very Long-Baseline Proposals: Experiments planning baselines beyond 5000 km.
2. Atmospheric Neutrino Analyses: Detectors like IceCube-Upgrade and KM3NeT/ORCA, where neutrinos traverse a continuous distribution of baselines and zenith angles, making accurate density modeling essential for simultaneous determination of $\delta_{CP}$ and mass ordering.
3. Systematic Frameworks: The need to incorporate continuous PREM profiles and propagate Earth model uncertainties (including regional deviations and transition zones) into the analysis frameworks of next-generation experiments to ensure that observed neutrino-antineutrino asymmetries reflect fundamental physics rather than environmental assumptions.