⚛️ phenomenology

NeutrinoOsc3Flavor: CP Phase Dependence in Three-Flavor Neutrino Oscillations: A Numerical Study in Vacuum and Matter

This paper introduces NeutrinoOsc3Flavor, a lightweight, NumPy-based Python framework that provides a transparent and pedagogical numerical solution for three-flavor neutrino oscillations in vacuum and matter, featuring exact CP phase dependence, cross-platform compatibility, and an independent analytical validation of matter effects using the Cardano method.

Original authors: Baktiar Wasir Farooq, Bipin Singh Koranga, Ansh Prasad, Imran Khan

Published 2026-02-04

📖 4 min read🧠 Deep dive

CC0 1.0

Original authors: Baktiar Wasir Farooq, Bipin Singh Koranga, Ansh Prasad, Imran Khan

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how a group of three friends (let's call them the "Flavor Trio") changes their personalities as they travel together on a long road trip. Sometimes they travel through empty space (Vacuum), and sometimes they travel through a crowded, busy city (Matter).

The paper you provided introduces a new, lightweight computer program called NeutrinoOsc3Flavor. Here is what it does, explained simply:

1. The Big Problem: "Black Box" Tools

In the world of particle physics, scientists use huge, complex computer programs to simulate how these "friends" (neutrinos) change. These big programs are like high-tech, sealed black boxes. They give you the answer, but you can't see how they got there. The math is hidden deep inside layers of code, making it hard for other scientists to double-check the work or understand the specific steps.

2. The Solution: A "Glass House" Calculator

The authors built NeutrinoOsc3Flavor to be the opposite of a black box. Think of it as a glass house.

Transparency: Every single step of the math is visible. You can see exactly how the program builds the equations, solves them, and gets the final result.
Simplicity: It doesn't need a massive, complicated installation. It runs on almost any computer (Windows, Mac, Linux) using just two common tools (NumPy and Matplotlib). It's like a simple, portable calculator rather than a supercomputer.
Purpose: It isn't meant to replace the big, heavy-duty tools used for massive experiments. Instead, it's a teaching tool and a "truth detector." It's designed so scientists can say, "Let's check the big tools with this simple one to make sure they aren't making mistakes."

3. How It Works: The "Magic" of Changing

The program studies how neutrinos "oscillate" (change flavors).

The Travelers: Neutrinos come in three types: Electron, Muon, and Tau.
The Journey: As they travel, they can morph from one type to another.
The "Secret Ingredient" (CP Phase): The program pays special attention to a hidden variable called the CP Phase. Imagine this as a "twist" or a "mood" in the universe. If you change this mood, the way the friends interact changes. The program tests many different "moods" (0°, 45°, 90°, etc.) to see how sensitive the results are to this twist.

4. Two Ways to Solve the Puzzle

The program is unique because it solves the same problem in two different ways to make sure the answer is right:

The Computer Way (Numerical): It uses standard math tricks (like a calculator) to find the answer. This is like using a GPS to find a route.
The Human Way (Analytical/Cardano Method): It also uses an old-school, exact mathematical formula (called the Cardano method) to solve the same puzzle by hand (but on the computer). This is like a mathematician solving the route on a piece of paper.

The Result: The paper shows that the "GPS" answer and the "Paper" answer match perfectly. This proves the computer isn't glitching and that the math is solid.

5. What They Found

In Empty Space (Vacuum): When the neutrinos travel through nothing, changing the "mood" (CP Phase) changes the probability of them switching flavors, but it doesn't change how likely they are to stay the same.
In a Crowd (Matter): When they travel through Earth's matter (like the crust), the crowd affects them differently. The "mood" (CP Phase) now changes both the switching probability and the staying-the-same probability. The program successfully maps out these changes.

Summary

NeutrinoOsc3Flavor is a clear, simple, and transparent computer program that lets scientists watch the math of neutrino changes happen in real-time. It uses a "double-check" system (computer math vs. old-school formulas) to ensure accuracy. It's not a replacement for the heavy-duty tools used in big experiments, but it is the perfect reference guide to make sure those big tools are telling the truth, especially when studying the mysterious "mood" (CP violation) of the universe.

Technical Summary: NeutrinoOsc3Flavor

Problem Statement
Modern neutrino oscillation studies, particularly those involving matter effects and CP-violating phase dependence, increasingly rely on large-scale, highly optimized software frameworks designed for experimental simulation. While indispensable for data analysis, these tools often function as "black boxes" where underlying numerical procedures, intermediate steps, and specific physics assumptions are embedded within complex software stacks. This opacity makes independent verification, equation-level inspection, and theoretical cross-checks difficult. There is a distinct need for transparent, minimal, and verification-oriented reference implementations that allow researchers to trace every step of the oscillation calculation explicitly.

Methodology
The paper presents NeutrinoOsc3Flavor, a lightweight, fully transparent computational framework written in pure Python. It is designed to solve the Schrödinger evolution equation for exact three-flavor neutrino oscillations in both vacuum and constant-density matter.

Formalism: The code operates within the PMNS formalism in the flavor basis. It explicitly constructs the effective Hamiltonian ( $H_f$ ) for vacuum and matter cases. For matter, it incorporates the charged-current potential ( $V_{CC}$ ) arising from electron density, while neglecting neutral-current effects (which cancel out in oscillation probabilities) and neutrino absorption.
Numerical Approach: The framework numerically diagonalizes the effective Hamiltonian using the NumPy library to obtain eigenvalues and eigenvectors, thereby constructing the evolution operator $S(L) = e^{-iH L}$ . This allows for the direct calculation of oscillation and survival probabilities.
Analytical Validation: A unique feature of the implementation is an independent analytical determination of the matter-modified Hamiltonian eigenvalues using the Cardano method. The code solves the characteristic polynomial of the $3\times3$ Hamiltonian matrix to derive eigenvalues analytically, providing a benchmark against the numerical diagonalization results.
Dependencies and Portability: The software relies solely on NumPy for linear algebra and Matplotlib for visualization. It requires no external C++/FORTRAN compilers, specialized libraries, or platform-specific installations, ensuring native execution on Linux, Windows, and macOS.
CP-Phase Sensitivity: The framework treats the CP-violating phase ( $\delta_{CP}$ ) as a primary variable, systematically varying it (e.g., $0^\circ, 45^\circ, 90^\circ, 135^\circ, 180^\circ, 217^\circ$ ) to study its impact on oscillation probabilities in both vacuum and matter.

Key Contributions

Transparent Reference Implementation: The code provides equation-level accessibility to the entire oscillation calculation, from Hamiltonian construction to the final probability evaluation, distinguishing it from opaque, large-scale toolkits.
Internal Consistency Check: By implementing the Cardano method for analytical eigenvalue determination alongside numerical diagonalization, the framework offers a robust internal validation mechanism. The authors demonstrate excellent agreement between the analytical and numerical eigenvalues.
CP-Phase Diagnostic Tool: The framework is specifically structured to use CP-phase dependence as a sensitive diagnostic for numerical stability and physical correctness, rather than treating it as a fixed parameter.
Pedagogical and Benchmarking Utility: The code is designed for theoretical cross-checks, educational use, and benchmarking more complex simulation software, rather than replacing full experimental simulation frameworks.

Results

Eigenvalue Agreement: The study confirms that the analytically derived eigenvalues (via Cardano's method) match the numerically computed eigenvalues (via NumPy diagonalization) to high precision (e.g., differences on the order of $10^{-16}$ to $10^{-13}$ ), validating the numerical solver's stability.
Vacuum vs. Matter Behavior:
- In vacuum, the oscillation probabilities vary significantly with the CP phase, while survival probabilities remain largely unaffected by the phase variation in the specific configurations shown.
- In matter, both oscillation and survival probabilities exhibit sensitivity to the CP phase, with matter effects inducing slight changes in survival probabilities that are not present in the vacuum case.
Validation: The results for muon neutrino survival probabilities in matter were compared against the NuOscProbExact toolkit, showing consistent behavior.

Significance
The authors position NeutrinoOsc3Flavor not as a replacement for comprehensive experimental simulation packages, but as a specialized, verification-oriented tool. Its significance lies in its analytical traceability and minimalist design, which lower the barrier to independent verification. By prioritizing the "how" of the calculation over raw throughput, it serves as a reliable foundation for:

Theoretical cross-checks of complex oscillation codes.
Educational applications where understanding the explicit steps of the Schrödinger evolution is paramount.
Benchmarking numerical stability, particularly regarding CP-violating phase dependence in three-flavor systems.

The framework's ability to run without specialized dependencies ensures reproducibility across different computing environments, making it a valuable complement to existing large-scale tools in neutrino physics.