The sensitivity of liquid scintillator detectors to CP-violation with atmospheric neutrinos

This paper investigates the potential of multi-kiloton liquid scintillator detectors to measure the CP-violating phase in atmospheric neutrino oscillations by calculating expected event rates and distributions across various sites and detector configurations, ultimately estimating sensitivity through a Poisson likelihood analysis.

The Gist

Imagine the universe is a giant, invisible ocean, and constantly raining down on us from space are tiny, ghostly particles called neutrinos. These particles are so shy they can pass through the entire Earth without bumping into anything. But every once in a while, they change their "costume" (or flavor) as they travel.

Scientists have been trying to solve a massive cosmic mystery: Why does our universe exist at all?

According to the laws of physics, the Big Bang should have created equal amounts of matter (us) and antimatter (the "anti-us"). If that were true, they would have annihilated each other instantly, leaving nothing but empty space. But we are here. Something must have tipped the scales. Scientists suspect a hidden "rule breaker" in the neutrino world, called CP-violation, is the culprit.

This paper is like a blueprint for a new kind of detective game to catch this rule breaker in action. Here is the story of how they plan to do it, using simple analogies.

1. The Detective's Tool: The Liquid Scintillator

The authors propose using a giant tank filled with liquid scintillator. Think of this liquid as a giant, glowing pool of Jell-O.

• When a neutrino (the ghost) finally decides to interact with the Jell-O, it creates a tiny flash of light, like a firefly blinking in a dark room.
• By counting these blinks and measuring how bright they are, scientists can figure out what kind of neutrino it was and where it came from.

2. The Mystery: The "Costume Change"

Neutrinos come in three flavors: Electron, Muon, and Tau. As they travel through the Earth, they oscillate, meaning they switch costumes.

• The Twist: The rules for how neutrinos switch costumes might be slightly different from how antineutrinos (their evil twin counterparts) switch costumes.
• If the rules are different, that difference is CP-violation. Finding this difference is the "Holy Grail" of particle physics because it explains why we exist.

3. The Challenge: The Foggy Window

The problem is that our detectors aren't perfect.

• The Blur: Imagine looking at a beautiful, intricate painting through a foggy window. You can see the general shapes and colors, but the fine details are smeared out.
• In this paper, the "fog" is the detector's resolution. It can't pinpoint the exact energy or direction of the neutrino perfectly. The authors calculated how much this "fog" blurs the signal. They found that while the fine details get lost, the big picture (the overall pattern of the costume changes) is still visible.

4. The Noise: The Background Static

In any experiment, there is "static" or noise.

• The Impostors: Sometimes, a neutrino interacts but doesn't leave a clear "costume" behind. It creates a messy burst of energy that looks like a real signal but isn't. The authors call this background noise (specifically from "neutral current" interactions).
• They calculated how many of these impostors would sneak into their data and figured out how to filter them out, kind of like a bouncer at a club checking IDs to keep the fake ones out.

5. The Superpower: Identifying the Flavor

The most critical part of this detective work is knowing which flavor the neutrino was when it arrived.

• The ID Card: The detector needs to be able to say, "This was an Electron-neutrino," and "That was a Muon-antineutrino."
• The Catch: If the detector is bad at reading these ID cards (low accuracy), the mystery remains unsolved. The paper shows that to get a strong answer (a "3-sigma" or "4-sigma" result, which is scientific speak for "very confident"), the detector needs to be right about 90% to 95% of the time. If it's only right 50% of the time, the signal gets lost in the noise.

6. The Verdict: Where to Build and What to Expect

The authors ran simulations for different locations (like deep underground mines in Canada, Italy, and Japan) and different sizes of detectors (a few thousand tons of liquid).

• Location doesn't matter much: Whether you are in Canada or Japan, the "rain" of neutrinos is similar enough that the location isn't the deciding factor.
• The Result: If we build a detector with a large tank of liquid scintillator and make sure it is really good at identifying neutrino flavors, we can detect this CP-violation with a confidence level of up to 4 sigma (which is very strong evidence, though not quite the "5 sigma" gold standard for a full discovery).

The Bottom Line

This paper is a feasibility study. It says: "Yes, we can catch the universe's biggest secret using a giant tank of glowing liquid, but we have to be very good at telling the difference between the different types of neutrinos."

If we build the right detector and tune it correctly, we might finally answer the question: Why is there something rather than nothing?

Technical Summary

1. Problem Statement

The determination of the CP-violating phase ($\delta_{CP}$) in the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix is a primary goal of next-generation neutrino physics. Current constraints on $\delta_{CP}$ are limited, and measurements often suffer from degeneracies related to the unknown neutrino mass ordering (Normal Ordering vs. Inverted Ordering) and the octant of the mixing angle $\theta_{23}$.

While long-baseline accelerator experiments are the primary tool for measuring $\delta_{CP}$, this paper investigates the potential of atmospheric neutrinos detected by liquid scintillator detectors (typically a few kilotons in mass). These detectors are sensitive to the low-energy end of the atmospheric neutrino flux (< 10 GeV). The core challenge is that detecting CP violation requires distinguishing between neutrinos and antineutrinos and identifying their flavors ($\nu_e, \nu_\mu, \bar{\nu}_e, \bar{\nu}_\mu$) with high precision, while accounting for detector resolution effects, background noise, and systematic uncertainties.

2. Methodology

The authors employ a comprehensive simulation and statistical analysis framework to estimate the sensitivity of a typical liquid scintillator detector (using SNOLAB as a reference site) to $\delta_{CP}$.

• Signal Modeling:

  • Flux: Uses the Honda et al. model for unoscillated atmospheric neutrino fluxes.
  • Oscillations: Calculates flavor conversion probabilities ($P_{\beta \to \alpha}$) using nuCraft, incorporating the PREM Earth density profile and PMNS parameters from NuFIT 5.2. The analysis varies $\delta_{CP}$ from $-180^\circ$ to $180^\circ$ and considers both Normal Ordering (NO) and Inverted Ordering (IO), as well as different $\theta_{23}$ octants.
  • Event Rate: Calculates expected charged current (CC) event rates for an exposure of 250 kt-years, integrating over energy ($E_\nu$) and zenith angle ($\theta$).

• Detector Response:

  • Resolution: Models energy resolution as a constant 10% ($\sigma_E = 0.1 E_\nu$) and zenith angle resolution based on extrapolations from JUNO data. These resolutions are applied via a 2D Gaussian smearing, which washes out fine oscillation structures but preserves cumulative deficits/excesses.
  • Backgrounds: Dominated by Neutral Current (NC) interactions of atmospheric neutrinos, which can mimic CC events via high-energy pion tracks. The study assumes a 5% misidentification rate of NC events into signal classes and an 80% signal selection efficiency.

• Flavor Identification:

  • Introduces two key performance metrics: Flavor identification accuracy ($\alpha_f$) and Neutrino/Antineutrino identification accuracy ($\alpha_\pm$).
  • Constructs a response matrix ($E_{\alpha\beta}$) to model how true flavors map to reconstructed flavors, accounting for misidentification.

• Statistical Analysis:

  • Uses a Poissonian likelihood ratio test to compare a true hypothesis ($\delta_{CP}^0$) against test hypotheses ($\delta_{CP}^{test}$).
  • Treats flux normalization ($N_F$) and cross-section uncertainties as nuisance parameters.
  • Cross-section uncertainties are modeled by mixing predictions from GENIE, NuWro, and Nuance generators.
  • Significance is calculated using Wilks' Theorem, defined as the minimum significance of excluding $\delta_{CP} = 0^\circ$ or $180^\circ$.

3. Key Contributions

• Systematic Sensitivity Study: Provides a quantitative assessment of how specific detector capabilities (flavor ID, energy/angle resolution) impact the ability to measure $\delta_{CP}$ using atmospheric neutrinos in liquid scintillators.
• Background and Systematics Integration: Explicitly models the impact of Neutral Current backgrounds and cross-section model variations (GENIE vs. NuWro vs. Nuance) on the final sensitivity.
• Site Independence: Demonstrates that detector location (SNOLAB, Gran Sasso, Kamioka) has a minor impact on sensitivity compared to detector performance metrics.
• Performance Thresholds: Establishes specific accuracy requirements for flavor and particle/antiparticle identification to achieve statistically significant results.

4. Key Results

• Sensitivity Range: For a 250 kt-year exposure, the maximal sensitivity to CP violation (at $\delta_{CP} = \pm 90^\circ$) ranges from $1\sigma$ to $4\sigma$, heavily dependent on flavor identification capabilities.
• Accuracy Requirements:
  • To achieve $>3\sigma$ sensitivity, an overall identification accuracy ($\alpha_f \approx \alpha_\pm$) of approximately 90% is required.
  • To achieve $4\sigma$ sensitivity, accuracies must exceed 95%.
  • If identification accuracy drops significantly, sensitivity can fall below $1\sigma$.
• Impact of Systematics:
  • Cross-section Uncertainty: Switching between different cross-section models (e.g., injecting NuWro data while fitting with GENIE) reduces sensitivity by less than $0.5\sigma$.
  • Flux Normalization: A $\pm 20\%$ uncertainty in flux normalization tightens the requirements for selection performance but does not fundamentally alter the sensitivity curve shape.
• Mass Ordering and Octant:
  • Sensitivity is generally lower for Inverted Ordering (IO) than Normal Ordering (NO).
  • The specific octant of $\theta_{23}$ (upper vs. lower) has a slight effect, with the lower octant (associated with NO in some fits) showing marginally higher sensitivity.
• Detector Location: The choice of site (SNOLAB, LNGS, Kamioka) has a negligible effect on the final sensitivity, as the variation in initial atmospheric flux is outweighed by the oscillation effects and detector performance.

5. Significance

This study validates the potential of multi-kiloton liquid scintillator detectors as competitive tools for probing CP violation using atmospheric neutrinos, provided they achieve high-fidelity flavor and particle/antiparticle identification.

The results indicate that while detector resolution smears out the detailed oscillation patterns, the cumulative effects of CP violation on the energy and zenith angle distributions remain detectable. The paper establishes that the limiting factor is not the detector size or location, but rather the flavor identification performance. If future detectors can achieve $\sim90-95\%$ accuracy in distinguishing $\nu_e$ from $\nu_\mu$ and $\nu$ from $\bar{\nu}$, they can contribute significantly to the global understanding of CP violation, potentially offering an independent cross-check to long-baseline accelerator experiments.