A Consistent Treatment of Final-State Interactions in NuWro Quasielastic Channel

This paper presents a unified framework within the NuWro Monte Carlo generator that consistently treats final-state interactions in quasielastic lepton-nucleus scattering by combining convolution-based inclusive calculations with an event-level cascade classification, significantly improving agreement with both inclusive electron-scattering data and exclusive MicroBooNE measurements.

The Gist

The Big Picture: The "Neutrino Detective" Problem

Imagine you are a detective trying to solve a crime inside a crowded, chaotic stadium (the atomic nucleus). You throw a ball (a neutrino) into the crowd, hoping it hits a specific person (a proton or neutron) and bounces off cleanly so you can see exactly who it hit and how hard.

This is what scientists do in neutrino experiments. They fire neutrinos at atoms to learn about the fundamental laws of physics. However, the "stadium" (the nucleus) is messy. When the neutrino hits a nucleon, that nucleon doesn't just fly out in a straight line. It often bumps into other people in the crowd, changes direction, loses energy, or even knocks someone else over before it escapes.

These messy collisions are called Final-State Interactions (FSI).

The problem is that current computer programs used to simulate these experiments (like NuWro) were struggling to tell the difference between:

1. The Clean Hit: The nucleon flies out without bumping into anyone (Transparent).
2. The Messy Hit: The nucleon bounces around the crowd before escaping (Non-Transparent).

If the computer gets this wrong, the scientists' calculations for the "crime scene" are off, leading to wrong conclusions about the neutrino itself.

The Solution: A New "Traffic Cop" System

The author, Rwik Dharmapal Banerjee, has updated the NuWro computer program with a smarter way to handle these collisions. Think of it as installing a new set of traffic rules for the particles.

Here is how the new system works, broken down into three simple steps:

1. The "Crystal Ball" Prediction (The Spectral Function)

Before the simulation even starts, the program uses a sophisticated math model called the Spectral Function. Imagine this as a crystal ball that predicts the odds of the nucleon hitting a "clean exit" versus a "messy exit."

• It calculates a number called Nuclear Transparency ($T_A$).
• High Transparency: The nucleon has a high chance of zooming out without hitting anyone.
• Low Transparency: The nucleon is almost guaranteed to get stuck in traffic and crash into others.

2. The "Tagging" System (The Algorithm)

In the old version, the computer would just let the nucleon wander and see what happened. In the new version, the computer tags every single event before it simulates the movement.

• Tag: "Transparent"

  • The Analogy: Imagine a VIP guest at a party who is given a direct, unobstructed path to the exit.
  • The Action: The computer forces the struck nucleon to fly straight out. It is not allowed to hit anyone else. It's a "clean" escape.

• Tag: "Non-Transparent"

  • The Analogy: Imagine a guest who is stuck in a mosh pit. They must bump into people before they can get out.
  • The Action: The computer forces the nucleon to crash into at least one other person inside the nucleus. If it tries to escape without hitting anyone, the computer says, "No, that's not allowed for this tag," and forces it to crash again until it satisfies the rule.

3. The "Double-Check" (Consistency)

The genius of this paper is that it ensures the "Crystal Ball" prediction (the math) matches the "Party Simulation" (the event).

• If the math says 30% of particles should escape cleanly, the simulation ensures exactly 30% of the events are tagged "Transparent."
• This creates a perfect bridge between the big-picture math and the individual particle stories.

The Results: Why It Matters

The author tested this new system in two ways:

1. The "Electron Test" (Inclusive Data):
   They compared the simulation to data from electron scattering (a similar type of experiment).

   • Result: The new model (with the traffic cop) matched the real-world data perfectly. The old model (without the rules) was too "sharp" and didn't account for the messy energy loss. The new model showed the correct "blur" and shift in energy, just like the real data.

2. The "MicroBooNE Test" (Exclusive Data):
   They looked at specific data from the MicroBooNE experiment, which counts exactly how many particles come out of the nucleus.

   • Result: The new model fixed a major error. Previously, the computer predicted too many particles flying straight out. The new model, by forcing the "messy" collisions, reduced the number of straight-line particles and matched the real experiment much better.

The Bottom Line

Think of this paper as upgrading the GPS in a self-driving car.

• Before: The GPS knew the destination but didn't account for traffic jams or roadblocks, so it predicted the car would arrive faster and in a straighter line than reality.
• After: The new GPS (NuWro with FSI) knows exactly when a car will get stuck in traffic and when it will have a clear road. It simulates the bumps and the delays accurately.

By making the computer simulation match the messy reality of the atomic nucleus, scientists can now trust their data more. This is crucial for future experiments trying to understand why the universe is made of matter instead of antimatter, or to detect neutrinos from distant supernovas.

In short: The author taught the computer to stop pretending the atomic nucleus is an empty hallway and start treating it like a crowded room where collisions are inevitable. This makes the predictions much more accurate.

Technical Summary

1. Problem Statement

Accurate modeling of neutrino–nucleus interactions in the few-GeV energy region is critical for current and future neutrino oscillation experiments. The Charged-Current Quasielastic (CCQE) scattering channel is dominant in this regime but suffers from significant theoretical uncertainties due to nuclear effects. Specifically, Final-State Interactions (FSI)—the propagation of the struck nucleon through the nuclear medium—play a crucial role in shaping both inclusive cross-sections and exclusive final-state observables.

The core problem addressed is the lack of consistency between:

1. Inclusive descriptions: Where FSI is often modeled via convolution with a folding function to modify the cross-section.
2. Exclusive descriptions: Where FSI is simulated via intranuclear cascade models tracking individual particle trajectories.

Previous approaches often treated these levels separately, leading to potential mismatches between the theoretical cross-section predictions and the simulated event-level kinematics.

2. Methodology

The author proposes a unified framework within the NuWro Monte Carlo generator (version 25.11) that establishes a direct mapping between the inclusive FSI formalism and the event-level intranuclear cascade.

A. Theoretical Framework: Spectral Function (SF)

The work utilizes the Spectral Function approach to describe the nuclear ground state, which accounts for shell structure and short-range correlations (superior to simple Fermi gas models).

• Inclusive Formalism: The inclusive cross-section with FSI is calculated using a convolution scheme:
  $$\frac{d\sigma_{FSI}}{d\omega d\Omega} = \int d\omega' f_q(\omega - \omega') \frac{d\sigma_{PWIA}}{d\omega' d\Omega}$$
  Here, $f_q$ is a folding function dependent on Nuclear Transparency ($T_A$) and a finite-width function $F_q(\omega)$ that accounts for broadening.
• Energy Shift: The model incorporates the real part of the optical potential ($U_V$) to shift the energy spectrum of the outgoing nucleon.

B. Event-Level Implementation: The "Tagging" Algorithm

The key innovation is enforcing the inclusive folding function logic directly within the NuWro intranuclear cascade at the event level. Events are classified based on nuclear transparency:

1. Transparent Events:
   • The struck nucleon ($N_1$) is propagated without any rescattering (exits the nucleus undisturbed).
   • If the event is "correlated" (involving a partner nucleon $N_2$), $N_2$ undergoes standard cascade propagation.
2. Non-Transparent Events:
   • The struck nucleon ($N_1$) is forced to undergo at least one intranuclear collision.
   • Algorithmic Logic: The code attempts to propagate $N_1$. If it does not interact, the history is discarded, the mean free path is reduced to increase interaction probability, and the propagation is redrawn until an interaction occurs.
   • Correlated partners ($N_2$) are propagated via the standard cascade.

This ensures that the statistical distribution of "transparent" vs. "non-transparent" events in the simulation matches the theoretical folding function used for inclusive cross-sections.

3. Key Contributions

• Unified Framework: Successfully bridges the gap between inclusive cross-section calculations (convolution-based) and exclusive event generation (cascade-based) within the NuWro generator.
• Consistency: Guarantees that the physics governing inclusive observables (total cross-sections) is identical to the physics governing exclusive final states (particle multiplicities and kinematics).
• Algorithmic Innovation: Introduced a rejection-sampling-like algorithm within the cascade to strictly enforce the "transparent" vs. "non-transparent" classification derived from nuclear transparency parameters.
• Validation: Tested the framework against both electron-scattering data (inclusive) and MicroBooNE neutrino data (exclusive).

4. Results

The implementation was validated against two distinct datasets:

• Inclusive Electron-Scattering (Carbon):

  • Observation: Including FSI led to a visible redistribution of strength, causing a broadening of the cross-section and a shift of the Quasielastic (QE) peak.
  • Outcome: The full calculation (with FSI) showed excellent agreement with experimental data from Barreau et al., whereas the Plane Wave Impulse Approximation (PWIA) without FSI failed to match the data.

• Exclusive Neutrino Data (MicroBooNE):

  • Observable: Double-differential cross-sections for CC1p0$\pi$ (Charged Current, 1 proton, 0 pions) interactions, specifically looking at the transverse kinematic imbalance variable ($\delta p_T$).
  • Observation: The inclusion of FSI resulted in a pronounced reduction of the peak in the $\delta p_T$ distribution.
  • Outcome: The FSI-included predictions showed significantly improved agreement with MicroBooNE measurements compared to the non-FSI baseline.
  • Limitation: A deficit in the cross-section was observed at higher transverse momentum regions, suggesting that the intranuclear cascade strength in NuWro might still require tuning in that specific kinematic region.

5. Significance

This work is significant for the neutrino physics community because:

1. Precision Oscillation Physics: Reducing uncertainties in FSI modeling is essential for minimizing systematic errors in neutrino oscillation experiments (e.g., DUNE, Hyper-K).
2. Model Reliability: By ensuring consistency between inclusive and exclusive descriptions, the framework provides a more robust theoretical tool for interpreting complex neutrino-nucleus data.
3. Generator Improvement: The updated NuWro implementation (v25.11) offers a state-of-the-art tool for simulating QE interactions in Carbon, Oxygen, Argon, and Iron, directly incorporating modern Spectral Function data and a rigorous FSI treatment.

In conclusion, Banerjee's work demonstrates that a consistent, event-level enforcement of inclusive FSI formalisms significantly enhances the predictive power of Monte Carlo generators, bringing theoretical models into closer alignment with experimental reality.