The role of final-state interaction modeling in neutrino energy reconstruction and oscillation measurements

This paper demonstrates that uncertainties in final-state interaction modeling can distort reconstructed neutrino-energy spectra in next-generation long-baseline experiments like DUNE to an extent comparable to or exceeding oscillation parameter variations, thereby creating potential degeneracies that necessitate improved theoretical characterizations and dedicated experimental measurements for robust results.

The Gist

The Big Picture: Trying to Weigh a Ghost

Imagine you are trying to weigh a ghost that flies through a wall. You can't see the ghost itself, but when it hits a pile of bricks (an atom), it knocks some bricks loose. By measuring how fast those loose bricks fly and how much energy they carry, you try to guess how heavy the ghost was when it hit.

In the world of physics, this "ghost" is a neutrino, and the "bricks" are argon atoms inside a giant detector called DUNE. Scientists want to know the exact energy of these neutrinos to understand the fundamental laws of the universe, like why the universe is made of matter instead of antimatter.

The Problem: The "After-Math" Mess

The paper argues that there is a huge problem with how we calculate that weight.

When the neutrino hits the atom, it doesn't just stop. It creates a chaotic explosion of smaller particles (like pions and protons). These particles are like billiard balls bouncing around inside a crowded room (the atomic nucleus) before they finally escape to be seen by the detector.

This bouncing around is called Final-State Interaction (FSI).

• The Analogy: Imagine throwing a tennis ball into a room full of people. The ball hits a person, bounces off a wall, hits another person, and loses some speed before it finally rolls out the door.
• The Issue: If you only measure the ball rolling out the door, you might think it was thrown much slower than it actually was, because you didn't account for all the bumps and bruises it took inside the room.

The Core Discovery: Guessing vs. Reality

The authors of this paper ran a simulation to see what happens if we use different "rules" to guess how those particles bounce around inside the nucleus. They used four different mathematical models (like four different weather forecasters) to predict the chaos.

Here is the shocking result:
Depending on which "rulebook" you use to calculate the bouncing, your estimate of the neutrino's energy changes dramatically.

• The Metaphor: Imagine you are trying to solve a mystery: "Did the suspect run 10 miles or 11 miles?"
  • If you use Rulebook A, you calculate the distance as 10.1 miles.
  • If you use Rulebook B, you calculate it as 10.9 miles.
  • The difference between 10.1 and 10.9 looks like a huge error, but in this experiment, that tiny difference is bigger than the actual mystery they are trying to solve.

Why This Matters: The "Fake" Clues

The scientists are trying to measure very subtle things, like a tiny shift in the neutrino's energy caused by a specific property of the universe (called an "oscillation parameter").

The paper shows that the confusion caused by the bouncing particles (FSI) creates a "fake signal" that looks exactly like the real signal scientists are hunting for.

• The Analogy: Imagine you are listening for a specific bird song in a forest. But, depending on how you interpret the wind rustling the leaves (the FSI), the wind sounds exactly like the bird song. You might think you found the bird, but it was just the wind.

If the scientists don't figure out exactly how the particles bounce (the wind), they might think they have discovered a new law of physics, when they have actually just misunderstood the noise.

The Solution: Better Maps and New Tools

The paper concludes that we cannot just guess anymore. We need:

1. Better Theory: We need to write better "rulebooks" for how particles bounce inside atoms.
2. New Experiments: We need to build special detectors (like the "Near Detector" mentioned in the paper) that sit right next to the neutrino source. These detectors act like a "control group" to watch the particles before they get lost in the main experiment, helping us calibrate our rules.
3. Direct Measurements: We need to shoot particles at argon in a lab and watch exactly what happens, rather than just guessing with math.

Summary

This paper is a warning label for the future of neutrino physics. It says: "Stop! Before you claim you've found the secrets of the universe, make sure you understand how the particles bounce around inside the atom. If you don't, the 'bouncing' will look like a discovery, and you'll be fooled by your own math."

It's a call to action to build better tools and understand the "messy middle" of the experiment so that the final answer is truly accurate.

Technical Summary

1. Problem Statement

Long-baseline (LBL) neutrino oscillation experiments, such as the Deep Underground Neutrino Experiment (DUNE) and Hyper-Kamiokande (HK), aim to precisely measure oscillation parameters (e.g., $\Delta m^2_{32}$, $\delta_{CP}$) and determine the neutrino mass ordering. A critical prerequisite for these measurements is the accurate reconstruction of the incident neutrino energy ($E_{\nu}$) from the visible products of neutrino-nucleus interactions.

The primary challenge addressed in this work is Final-State Interactions (FSI). When a neutrino interacts with a nucleus (e.g., Argon in DUNE), the resulting hadrons (protons, pions, etc.) undergo further scattering, absorption, or charge exchange within the nucleus before exiting. These interactions alter the observable particle multiplicities and kinematics, leading to "invisible" energy loss (primarily carried away by neutrons).

• The Core Issue: Different theoretical models used to simulate FSI yield significantly different predictions for the visible energy spectrum.
• The Risk: If these modeling uncertainties are not properly constrained, the resulting distortions in the reconstructed energy spectrum ($E_{rec}$) could mimic or obscure the effects of genuine neutrino oscillation parameters, creating a degeneracy that limits the precision of next-generation experiments.

2. Methodology

The authors conducted a quantitative case study using the DUNE experimental configuration (1300 km baseline, specific flux profiles) to assess the impact of FSI modeling variations.

• Simulation Framework:

  • Generator: GENIE v3.04.00 was used to generate charged-current (CC) neutrino and antineutrino interactions on Argon.
  • Baseline Model: The G18_10a_02_000 configuration was used, employing a local Fermi-gas nuclear ground state, Valencia group models for CCQE and 2p2h, Berger-Seghal for pion production, and GRV98 PDFs for Deep Inelastic Scattering (DIS).
  • FSI Models Compared: Four distinct FSI models available in GENIE were compared to establish an uncertainty envelope:
    1. hA2018: An empirical, single-step transport model (baseline).
    2. hN2018: A full intra-nuclear cascade (INC) model (space-like).
    3. INCL: The Liege Intranuclear Cascade model (time-like, includes quantum effects like Pauli blocking).
    4. G4BC: The GEANT4 Bertini Cascade model (space-like).
  • Oscillation Application: Events were weighted using the Prob3++ software with standard oscillation parameters (Normal Mass Ordering, $\delta_{CP}=0$) to simulate the far detector spectrum.

• Energy Reconstruction Metric:
  The authors defined a proxy for reconstructed neutrino energy ($E_{rec}^{\nu}$) assuming a perfect calorimeter:
  $$E_{rec}^{\nu} = E_{\ell} + \sum T_{p} + \sum E_{\pi, K, \gamma}$$
  Where $E_{\ell}$ is lepton energy, $T_p$ is proton kinetic energy, and $E_{\pi, K, \gamma}$ are total energies of mesons and photons. The relative energy bias was defined as $(E_{rec}^{\nu} - E_{true}^{\nu}) / E_{true}^{\nu}$.

• Comparison Strategy:
  The distortions in the energy spectra caused by switching FSI models were directly compared to spectral shifts induced by varying oscillation parameters ($\Delta m^2_{32}$ and $\delta_{CP}$) at the level of precision projected for next-generation experiments (e.g., $\pm 0.4\%$ for $\Delta m^2_{32}$).

3. Key Contributions

• Quantitative Degeneracy Demonstration: The paper provides the first quantitative demonstration that realistic variations in FSI modeling can alter reconstructed energy spectra by amounts comparable to or larger than the shifts induced by oscillation parameters at projected experimental precisions.
• Model Spread as Uncertainty Envelope: By comparing four distinct FSI models (ranging from empirical to sophisticated cascade models), the authors define a "minimal uncertainty envelope" that highlights the lack of consensus in current theoretical descriptions of hadron transport in nuclei.
• Differentiation of CP Sensitivity: The study distinguishes between the impact of FSI on the discovery of CP violation (rate-based) versus the precision measurement of $\delta_{CP}$ (shape-based).

4. Key Results

• Significant Energy Bias: FSI introduces a substantial bias in energy reconstruction. For the models considered, ~40% of neutrino events and ~50% of antineutrino events exhibit a relative energy bias greater than 10%. The bias distribution is broader for antineutrinos due to their tendency to produce more neutrons in the final state.
• Degeneracy with $\Delta m^2_{32}$:
  • Variations in FSI models produce spectral shifts and smearing near the first oscillation maximum that are indistinguishable from a $\pm 0.4\%$ shift in $\Delta m^2_{32}$.
  • Specifically, the G4BC and INCL models showed deviations exceeding the spectral shifts caused by a 0.4% change in $\Delta m^2_{32}$.
• Impact on $\delta_{CP}$ Measurements:
  • Discovery vs. Precision: If CP violation is maximal ($\delta_{CP} \approx -\pi/2$), the discovery of CPV (driven by rate differences between $\nu_e$ and $\bar{\nu}_e$) remains robust against FSI variations.
  • Precision Limit: However, for precision measurements of $\delta_{CP}$, where the signal relies on the shape of the energy spectrum, FSI modeling uncertainties are critical. When $\delta_{CP} \approx -\pi/2$, variations in $\delta_{CP}$ manifest as shape changes in the spectrum, which closely resemble the distortions introduced by FSI model switching. This creates a potential degeneracy that could bias the extracted value of $\delta_{CP}$.
• Antineutrino Sensitivity: The impact of FSI is generally more pronounced for antineutrinos due to the higher production of neutrons (invisible energy) and the broader distribution of energy bias.

5. Significance and Outlook

• Urgent Need for Characterization: The study underscores that current FSI modeling uncertainties are a limiting factor for the physics goals of DUNE and similar experiments. Without improved characterization, these uncertainties could dominate the systematic error budget.
• Role of Near Detectors: The authors argue that robust constraints require a theory-driven uncertainty parameterization that can be fitted to near detector data. Current parameterizations are insufficient to fully constrain the full range of FSI-induced modifications.
• Future Directions:
  • Experimental: Dedicated measurements of hadron-argon cross-sections (e.g., via LArIAT, ProtoDUNE, and the proposed nuSCOPE experiment) are essential to benchmark FSI models.
  • Theoretical: Development of more sophisticated uncertainty parameterizations is needed to allow near detectors to effectively constrain far detector systematics.
  • Observables: Utilizing alternative observables (e.g., transverse kinematic imbalance) may help disambiguate FSI effects from oscillation signals.

Conclusion: The paper concludes that while FSI variations alone do not preclude the discovery of new physics, they pose a severe threat to the precision of oscillation parameter measurements. Addressing this requires a concerted effort combining dedicated experimental measurements of hadron-nucleus interactions with advanced theoretical modeling to reduce the degeneracy between nuclear effects and fundamental neutrino properties.