Charged-Current Neutrino-Induced Single-Pion Production in the Superscaling Approach and Relativistic Distorted-Wave Impulse Approximation

This paper presents a detailed comparison of the SuSAv2 and RDWIA theoretical models against experimental measurements of charged-current neutrino-induced single-pion production from T2K, MINERvA, and MiniBooNE across a broad energy range, highlighting their respective approaches to modeling nuclear targets and pion production channels.

The Gist

Imagine you are trying to solve a massive cosmic mystery: Where do neutrinos go, and how do they change their identity?

Neutrinos are ghostly, tiny particles that zip through the universe, passing through planets and people without leaving a trace. To catch them, scientists build giant detectors filled with water or oil. But here's the problem: Neutrinos are so shy that they rarely interact. When they do hit something, they smash into the atomic nuclei inside the detector, creating a chaotic explosion of smaller particles.

To understand the neutrino's story, scientists have to reconstruct the crash from the debris. But to do that, they need to know exactly how the nucleus reacts when it gets hit. This is where this paper comes in.

The Big Picture: Two Maps for a Crash

The authors of this paper are like cartographers trying to draw the most accurate map of a nuclear crash. They are comparing two different "theoretical maps" (mathematical models) to see which one better predicts what happens when a neutrino smashes into a nucleus and knocks out a single pion (a type of particle).

The two maps they are comparing are:

1. SuSAv2 (The "SuperScaling" Map): Think of this as a map that looks at the nucleus as a whole crowd of people. It uses a clever trick called "scaling" to predict how the crowd moves based on how a single person moves. It uses data from a model called DCC (Dynamic Coupled-Channels), which is like a very detailed rulebook for how particles interact.
2. RDWIA (The "Distorted Wave" Map): This map looks at the crash more individually. It treats the nucleus like a complex landscape with hills and valleys (a "potential") that distort the path of the particles flying out. It uses a different rulebook called the Hybrid model, which mixes low-energy and high-energy physics rules.

The Experiment: The "Target"

The scientists tested these maps against real-world data from three famous neutrino experiments: MiniBooNE, MINERvA, and T2K.

• Imagine these experiments as different-sized bowling alleys.
• MiniBooNE uses a low-energy "bowling ball" (neutrino).
• MINERvA uses a heavier, faster ball.
• T2K is somewhere in between.
• The "pins" they are trying to knock over are nuclei made mostly of Carbon (like the carbon in your body).

What They Found: The "Ghost" Problem

When the scientists ran their simulations, they found some interesting discrepancies:

1. The Models Disagree on the Details
Just like two weather forecasters might disagree on whether it will rain or snow, the SuSAv2 and RDWIA models often predicted different amounts of pions.

• For some types of pions (like the positive $\pi^+$), the models were close.
• For others (like the neutral $\pi^0$), the models gave very different answers.

2. The "Missing Pieces" (The Underestimation)
In several cases, especially with MiniBooNE and MINERvA, both models predicted fewer pions than the experiments actually saw.

• The Analogy: Imagine you are watching a car crash. You expect to see 5 pieces of glass fly out. Your map predicts 3. But in reality, 5 pieces fly out.
• Why? The paper suggests the models are missing "intra-nuclear cascade effects." Think of this as a game of billiards inside the nucleus. When the neutrino hits a particle, it doesn't just fly straight out. It might bounce off other particles inside the nucleus, changing its path or even turning into a different type of particle (like a $\pi^+$ turning into a $\pi^0$). The current maps don't fully account for this chaotic bouncing.

3. The "Rulebook" Matters More Than the "Landscape"
One of the most surprising findings was that the difference between the two maps (SuSAv2 vs. RDWIA) was actually less important than the difference between the two rulebooks (DCC vs. Hybrid) used inside them.

• Analogy: It doesn't matter if you use a GPS app or a paper map (the framework); if the traffic data (the rulebook) is wrong, your directions will be wrong. The way the fundamental particles interact (the "elementary" physics) is the biggest source of uncertainty.

Why Does This Matter?

You might ask, "Who cares about a few extra pions?"

Neutrinos are the key to understanding the universe's biggest secrets:

• Why is there more matter than antimatter? (Why did we exist?)
• What is the hierarchy of neutrino masses?

To answer these, scientists need to measure neutrino oscillations with extreme precision. If their "map" of how neutrinos interact with nuclei is wrong, their measurements of the oscillations will be wrong, and they might draw the wrong conclusions about the universe.

The Conclusion

The paper concludes that while both models are getting better, neither is perfect yet.

• They need to better understand the "bouncing" inside the nucleus (rescattering).
• They need to refine the "rulebooks" for how particles interact, especially the tricky "axial" forces that are hard to measure.
• Future experiments with even bigger detectors (like DUNE and Hyper-K) will provide more data to help these cartographers draw a perfect map.

In short: This paper is a rigorous check-up of our theoretical tools. It tells us, "We are on the right track, but we need to fix the details of the crash reconstruction before we can fully trust our map of the neutrino universe."

Technical Summary

1. Problem Statement

Neutrino oscillation experiments (e.g., T2K, DUNE, Hyper-K) rely on precise reconstruction of neutrino energy to determine oscillation parameters. A major source of uncertainty in these reconstructions is the modeling of neutrino-nucleus interactions, particularly charged-current single-pion production (CC1$\pi$).

• Context: In the energy range of hundreds of MeV to a few GeV, while quasielastic (QE) scattering dominates, inelastic processes like single-pion production (often via $\Delta$ resonance excitation) are crucial.
• The Challenge: Existing theoretical models struggle to consistently describe CC1$\pi$ cross-sections across different experiments (MiniBooNE, MINERvA, T2K) and kinematic regimes. There is a need to reduce nuclear-medium uncertainties and validate the description of the elementary boson-nucleon-pion vertex and subsequent nuclear effects.

2. Methodology

The authors perform a comparative analysis of two distinct theoretical frameworks applied to CC1$\pi$ production on nuclear targets (primarily $^{12}$C and hydrocarbon).

A. Theoretical Frameworks

1. SuSAv2-inelastic (SuperScaling Approach version 2):

   • Nuclear Dynamics: Based on superscaling properties observed in inclusive electron scattering. It expresses the cross-section as a product of a single-nucleon cross-section and a scaling function ($f$) that encapsulates nuclear dynamics.
   • Elementary Vertex: Utilizes the ANL-Osaka Dynamical Coupled-Channels (DCC) model. This allows for a separation of pion production channels ($\pi^+, \pi^-, \pi^0$) and distinguishes between inclusive contributions and specific single-pion channels ($\pi$-DCC).
   • Structure Functions: Uses single-nucleon inelastic structure functions ($W_{1-5}$) derived from the DCC model, accounting for multi-meson final states in the inclusive case but isolating specific SPP channels for comparison.

2. RDWIA (Relativistic Distorted-Wave Impulse Approximation):

   • Nuclear Dynamics: Solves the Dirac equation for the outgoing nucleon using an energy-dependent Relativistic Mean Field (EDRMF) potential. This accounts for Pauli blocking and current conservation naturally. The pion is treated as a plane wave.
   • Elementary Vertex: Employs the Hybrid model (developed by the Ghent group), which combines low-energy models with high-energy Regge-based models to describe the boson-pion-nucleon vertex across a broad energy range.

B. Experimental Data Comparison

The models are tested against flux-averaged differential cross-sections from three major experiments:

• MiniBooNE: Mineral oil ($CH_2$) target; neutrino flux peaking at $\sim 0.5$ GeV.
• MINERvA: Hydrocarbon ($CH$) target; broad flux extending up to 20 GeV (peaking $\sim 3.5$ GeV).
• T2K: $C_8H_8$ target; neutrino flux peaking at $\sim 0.6$ GeV.

The analysis covers various kinematic variables: muon momentum, muon kinetic energy, and scattering angles, for both neutrino and antineutrino beams and different pion charge states.

3. Key Contributions

• Systematic Comparison: This work provides the first detailed side-by-side comparison of the SuSAv2-DCC and RDWIA-Hybrid approaches specifically for CC1$\pi$ production across multiple experiments and kinematic regimes.
• Decoupling Nuclear vs. Elementary Effects: The study isolates the impact of the elementary vertex (DCC vs. Hybrid) from the nuclear modeling (SuSA vs. RDWIA). It demonstrates that variations at the nucleon level (between DCC and Hybrid) are generally more significant than those introduced by the nuclear model choice.
• Channel Separation: The use of the DCC model within SuSAv2 allows for a rigorous separation of $\pi^+$, $\pi^0$, and $\pi^-$ channels, facilitating a granular comparison with experimental data.
• Identification of Missing Physics: The authors highlight that current models consistently underestimate data in specific channels (notably $\pi^0$), pointing to the necessity of including intra-nuclear cascade effects (pion rescattering) and refining axial form factors.

4. Key Results

A. Free Nucleon Level (DCC vs. Hybrid)

• $\Delta$ Resonance Region ($W_X < 1.4$ GeV): The DCC model generally predicts higher cross-sections than the Hybrid model for $\pi^+$ and $\pi^-$ channels, while the Hybrid model predicts higher cross-sections for $\pi^0$.
• High Energy: Without the invariant mass cutoff, the Hybrid model generally predicts larger total cross-sections than DCC, except for the $\pi^+$ channel.

B. Nuclear Target Results

• MiniBooNE ($\sim 0.5-2$ GeV):
  • $\pi^0$: Both models significantly underestimate the data (SuSAv2 accounts for only $\sim 50\%$). The discrepancy is attributed to missing intra-nuclear cascade effects (e.g., $\pi^+ \to \pi^0$ conversion) and differences in the elementary vertex.
  • $\pi^+$: Both models agree reasonably well with data at higher muon energies but underestimate data at very forward angles and low-intermediate energies.
• MINERvA (Higher Energy):
  • $\pi^+$: Both models show good agreement with data, though slight overestimation occurs at small scattering angles ($5-10^\circ$).
  • $\pi^0$ (Neutrino): Both models significantly underestimate the data.
  • $\pi^0$ (Antineutrino): The RDWIA-Hybrid model agrees well with data, while SuSAv2-DCC underestimates it. This asymmetry suggests large uncertainties in the axial nucleon-to-resonance transition form factors.
  • $\pi^-$ (Antineutrino): RDWIA-Hybrid agrees well; SuSAv2-DCC underestimates.
• T2K: Results are consistent with MiniBooNE findings; both models track the data closely for $\pi^+$, with RDWIA generally yielding slightly smaller cross-sections than SuSAv2.

C. Nuclear Effects

• The reduction of the cross-section per bound nucleon compared to a free nucleon is observed in both models but is stronger in the SuSAv2 model than in the RDWIA approach, particularly at MINERvA kinematics.

5. Significance and Conclusions

• Model Limitations: Neither model provides a satisfactory description of all datasets. The persistent underestimation of $\pi^0$ production (especially in MiniBooNE and MINERvA neutrino modes) indicates a critical gap in understanding final-state interactions (FSI) and axial contributions.
• Axial Form Factors: The discrepancies between neutrino and antineutrino channels in the $\pi^0$ sector suggest significant uncertainties in the axial form factors for nucleon-to-resonance transitions in the second resonance region.
• Future Directions:
  • Intranuclear Cascades: The authors emphasize the need to implement intra-nuclear cascade effects (rescattering, charge exchange) into event generators (like GENIE or NEUT) using these theoretical frameworks to improve agreement with $\pi^0$ data.
  • Event Generators: The successful implementation of SuSAv2 and RDWIA in the QE and 2p2h regimes should be extended to the full inelastic regime for future oscillation analyses.
  • Data Needs: Future analyses with higher-energy data (DUNE, Hyper-K) and different targets (Argon, Water) are required to further constrain the models and the DCC/Hybrid limits.

In summary, this paper establishes that while the choice of nuclear model (SuSA vs. RDWIA) introduces variations, the choice of the elementary interaction model (DCC vs. Hybrid) is the dominant source of theoretical uncertainty. The work underscores the urgent need to refine axial form factors and incorporate detailed final-state interaction mechanisms to resolve the discrepancies between theory and experimental CC1$\pi$ data.