NuFast-Earth: Efficient Atmospheric, Solar, and Supernova Neutrino Propagation Through the Earth

The paper introduces NuFast-Earth, a fast, flexible, and accurate C++ algorithm designed to efficiently compute neutrino oscillation probabilities for atmospheric, nighttime solar, and supernova neutrinos traversing the Earth's varying density profiles, thereby addressing the growing computational demands of next-generation neutrino experiments.

The Gist

Imagine you are trying to predict the path of a ghost as it flies through a giant, multi-layered cake. This ghost is a neutrino, a tiny, almost massless particle that rarely interacts with anything. The cake is Earth, and the layers are different types of frosting and sponge (the core, the mantle, the crust) with varying densities.

Sometimes, these ghosts change their "costume" (flavor) as they fly through the cake. A ghost that started as a "muon" might turn into an "electron" or a "tau" by the time it reaches the other side. This phenomenon is called neutrino oscillation.

However, calculating exactly how they change costumes while flying through a cake that gets denser and denser toward the center is incredibly hard. It's like trying to solve a math problem where the rules of the game keep changing every inch you fly.

This paper, titled "NuFast-Earth," introduces a new, super-fast computer program designed to solve this problem. Here is the breakdown in simple terms:

1. The Problem: The "Slow" Ghost

Scientists are building massive detectors (like DUNE and HyperK) to catch these ghosts. To understand what they see, they need to run millions of simulations.

• The Old Way: Imagine trying to calculate the ghost's path by taking a tiny step, recalculating the rules, taking another step, and recalculating again. If you have to do this for millions of ghosts, millions of times, your computer will overheat and take days to finish.
• The Issue: The Earth isn't a smooth ball; it has sharp boundaries between the core and the mantle. Existing tools were either too slow or couldn't handle these sharp changes accurately.

2. The Solution: The "Smart" Shortcut

The authors (Peter Denton and Stephen Parke) created a new algorithm called NuFast-Earth. Think of it as a GPS for neutrinos that uses a few clever tricks to skip the boring parts of the math.

Here are the three main "magic tricks" they use:

A. The "Middle Man" Trick (The Tilde Basis)

Usually, to calculate the ghost's path, you have to carry around a heavy, complicated backpack of variables (specifically the mixing angles $\theta_{23}$ and the CP phase $\delta$).

• The Trick: The new algorithm says, "Let's take that heavy backpack off for the middle part of the journey." They do the hard math in a simplified "tilde" world where the Earth's layers are just numbers. They only put the backpack back on at the very beginning and the very end of the trip.
• The Result: This makes the math much lighter and faster, especially because scientists often want to test different values for those backpack variables. They can change the backpack without re-doing the whole flight path!

B. The "Mirror" Trick (Symmetry)

When a neutrino flies through the Earth, it goes down to the center and then back up.

• The Trick: The path going down is almost a mirror image of the path going up. The old way calculated both separately. NuFast-Earth realizes, "Hey, I already calculated the bottom half! I can just flip the result to get the top half."
• The Result: This cuts the work in half. It's like realizing you only need to bake one half of a cake and then mirror it to get the whole thing.

C. The "Cheat Sheet" (Caching)

If you are running a simulation where you change one tiny thing (like the angle of the ghost's entry), the old code would re-calculate the entire flight path from scratch.

• The Trick: NuFast-Earth keeps a "cheat sheet" (cache) of the parts of the calculation that don't change. If you only change the entry angle, it reuses the massive chunks of math it already did for the Earth's core.
• The Result: Changing certain parameters becomes instant. It's like having a library where you don't have to re-shelve every book every time you want to read a different chapter; you just flip to the page.

3. Why Does This Matter?

• Speed: The new code is orders of magnitude faster than existing tools. It can run simulations in seconds that used to take hours.
• Accuracy: It handles the sharp boundaries of the Earth's core perfectly, which is crucial for detecting subtle effects like the "day-night" effect in solar neutrinos (where neutrinos change slightly more at night because they pass through the Earth).
• Future-Proofing: As new experiments come online, they will need to process data at a speed that only this kind of efficient code can handle.

The Analogy Summary

Imagine you are a tour guide leading a group of tourists (neutrinos) through a complex, multi-story museum (the Earth).

• Old Method: You stop at every single exhibit to explain the history, calculate the distance, and write a report before moving to the next room. It takes forever.
• NuFast-Earth Method: You realize that the first floor and the last floor are identical. You calculate the middle section once, and then you just "mirror" the instructions for the return trip. You also realize that if the tourists just want to change their hats (a minor parameter), you don't need to re-tour the whole museum; you just hand them a new hat and keep walking.

In short: This paper gives scientists a super-fast, highly accurate tool to predict how neutrinos behave as they travel through our planet, unlocking the secrets of the universe much faster than before.

Technical Summary

Here is a detailed technical summary of the paper "NuFast-Earth: Efficient Atmospheric, Solar, and Supernova Neutrino Propagation Through the Earth" by Peter B. Denton and Stephen J. Parke.

1. Problem Statement

As next-generation neutrino experiments (e.g., DUNE, Hyper-Kamiokande, IceCube Upgrade, KM3NeT) come online, the demand for high-precision calculations of neutrino oscillation probabilities has increased dramatically.

• Computational Bottleneck: Extracting oscillation parameters (such as $\theta_{23}$, $\delta_{CP}$, and mass ordering) often requires millions of "throws" (statistical simulations) using techniques like Feldman-Cousins. The calculation of oscillation probabilities through the Earth is frequently the slowest part of this process.
• Complexity: Unlike long-baseline accelerator experiments where the baseline is fixed and matter density is relatively constant, atmospheric, solar (nighttime), and supernova neutrinos traverse the Earth with sharply varying density profiles. This requires handling arbitrary Earth models, varying production heights, and complex trajectories (core-crossing vs. mantle-only).
• Limitations of Previous Methods: Existing algorithms (like the authors' previous NuFast-LBL) were optimized for constant density or long-baseline scenarios but were not efficient enough for the multi-layered, variable-density propagation required for Earth-crossing neutrinos.

2. Methodology

The authors developed NuFast-Earth, a new algorithm and C++ implementation designed to handle arbitrary spherical Earth density models efficiently. The core methodology involves several key innovations:

A. The "Tilde" Basis and Hamiltonian Decomposition

To maximize computational speed, the algorithm separates the oscillation calculation into two parts:

1. The "Tilde" Basis ($\tilde{H}$): The authors rotate the Hamiltonian to remove the $\theta_{23}$ mixing angle and the CP-violating phase $\delta$ from the majority of the calculation.
   • The Hamiltonian in this basis is a real symmetric matrix, meaning its eigenvalues are real and eigenvectors can be chosen to be entirely real.
   • This reduces the computational complexity significantly, as complex arithmetic is minimized.
2. Final Rotation: The $\theta_{23}$ and $\delta$ rotations are applied only at the very beginning (atmosphere entry) and end (detector) of the trajectory. This allows these parameters to be varied without recalculating the entire Earth propagation.

B. Efficient Eigenvalue and Eigenvector Calculation

• Eigenvalues: Instead of using the computationally expensive exact Cardano formula for the cubic characteristic equation, the authors use an excellent analytical approximation (based on $\Delta m^2_{ee}$) followed by Newton-Raphson (NR) iterations.
  • They demonstrate that one NR iteration provides precision sufficient for all next-generation experiments while being $\sim 2\times$ faster than the exact method. Two iterations reach double-precision limits.
• Eigenvectors: Using the "Eigenvectors from Eigenvalues" theorem, the eigenvectors are calculated directly from the eigenvalues and the adjugate of the Hamiltonian, avoiding costly matrix diagonalization routines.

C. Trajectory Segmentation and Symmetry

The neutrino trajectory is divided into four segments:

1. A (Air): Production point to Earth surface (vacuum).
2. S (Surface): Surface to detector depth on the far side.
3. I (In): Detector depth to the trajectory's closest point to Earth's center.
4. O (Out): Closest point back to the detector.

Key Optimization: For a spherically symmetric Earth, the density profile of segment O is the exact reverse of segment I. Since the layer amplitudes are symmetric matrices, the amplitude for segment O is simply the transpose of segment I ($\tilde{A}_O = \tilde{A}_I^T$). This cuts the number of required matrix multiplications by nearly half for core-crossing trajectories.

D. Intelligent Caching

The code tracks which parameters have changed.

• "Fast" Parameters: $\theta_{23}$, $\delta$, production height, and solar density. Changing these requires only minimal recalculation (mostly the final rotation or atmospheric segment).
• "Slow" Parameters: $\theta_{12}$, $\theta_{13}$, $\Delta m^2_{21}$, $\Delta m^2_{31}$, and Earth density model. Changing these requires recalculating the full trajectory.
• Caching: Eigenvalues and eigenvectors are cached for specific energy and layer combinations, allowing rapid scanning over zenith angles without re-solving the characteristic equation.

3. Key Contributions

1. NuFast-Earth Code: A publicly available C++ library that handles atmospheric, solar, and supernova neutrinos propagating through arbitrary spherical Earth models.
2. Speed Optimization: The algorithm achieves speedups of 1 to 2 orders of magnitude compared to existing state-of-the-art codes (e.g., nuSQuIDS, OscProb) by exploiting the real-symmetric nature of the Hamiltonian and the symmetry of Earth trajectories.
3. Flexible Earth Modeling: Supports various models including the Preliminary Reference Earth Model (PREM) with arbitrary layering, constant density models, and integration-based models.
4. Validation: The code has been validated against:
   • Analytic approximations for solar neutrinos.
   • Standard long-baseline results (DUNE-like setups).
   • Comparison with other oscillation codes.

4. Results

• Performance Benchmarks:
  • Calculating eigenvalues/eigenvectors with 1 NR iteration is significantly faster than the exact Cardano method with negligible loss in precision.
  • Scanning over "fast" parameters (like $\theta_{23}$ or $\delta$) is extremely fast, as the bulk of the Earth propagation calculation is reused.
  • Scanning over zenith angles is approximately 2x faster than scanning over energy grids because eigenvalues/eigenvectors depend on energy but not zenith angle (for a fixed Earth model).
• Precision:
  • Using $\sim 32$ total layers in the Earth model yields a mean error of $\sim 10^{-4}$ and a maximum error $< 1\%$, which is sufficient for all current and future experiments.
  • The approximation error from using 1 NR iteration is $< 10^{-5}$, far below experimental sensitivity.
• Physics Insights:
  • The code confirms that detector depth (e.g., 2 km vs. surface) can shift oscillation probabilities by up to 5% for low-energy, horizontal atmospheric neutrinos, emphasizing the need for precise depth modeling.
  • It successfully reproduces the "regeneration effect" for nighttime solar neutrinos and parametric resonances for core-crossing atmospheric neutrinos.

5. Significance

The NuFast-Earth algorithm is a critical tool for the next generation of neutrino physics.

• Enabling Precision: It makes the computationally intensive task of fitting oscillation parameters to massive datasets (millions of events) feasible within reasonable timeframes.
• Versatility: It unifies the treatment of atmospheric, solar, and supernova neutrinos, allowing for consistent analysis of Earth tomography and matter effects across different sources.
• Future-Proofing: By providing a framework that is orders of magnitude faster than current methods, it ensures that experiments like DUNE and Hyper-K can fully utilize their statistical power to resolve the neutrino mass ordering, CP violation, and potentially probe the Earth's interior structure via neutrino tomography.

In summary, the paper presents a highly optimized, mathematically rigorous, and practically implemented solution to the computational bottleneck of neutrino oscillation calculations in variable matter densities, setting a new standard for efficiency in the field.