First inclusive triple-differential measurement of the muon-antineutrino charged-current cross section using the NOvA Near Detector

Using the largest published sample of approximately one million muon antineutrino events from the NOvA Near Detector, this study presents the first triple-differential measurement of the charged-current inclusive cross section, revealing significant energy- and angle-dependent discrepancies between the data and current neutrino event generator predictions.

The Gist

Imagine you are a detective trying to solve a mystery, but instead of looking for a stolen jewel, you are trying to understand the behavior of ghosts.

In the world of physics, these "ghosts" are neutrinos. They are tiny, invisible particles that zip through the universe (and right through your body) without bumping into anything. They are so elusive that to catch them, you need a massive trap and a very specific kind of bait.

This paper is a report from the NOvA Collaboration, a team of scientists who built a giant trap (the NOvA Near Detector) in Illinois to catch these ghosts. Here is what they did, explained simply:

1. The Setup: The Ghost Trap

The scientists used a beam of particles from a giant accelerator called Fermilab. They shot a stream of muon antineutrinos (a specific type of ghost) at a detector filled with oil and plastic.

• The Analogy: Imagine a dark room filled with thousands of tiny, floating dust motes (the oil). You shine a flashlight (the particle beam) into the room. Most of the light passes right through the dust without hitting anything. But occasionally, a photon hits a dust mote, causing a tiny flash of light.
• The Goal: The scientists wanted to count exactly how many times the "ghosts" hit the "dust" and, more importantly, how they hit.

2. The Big Catch: A Million Ghosts

In the past, scientists had only caught a few hundred or thousand of these specific "ghosts" (muon antineutrinos). This paper is a big deal because they caught one million of them.

• The Analogy: Before, it was like trying to guess the weather by looking at a single raindrop. Now, they have a bucket full of rain. With a sample this huge, they can see patterns they never saw before.

3. The New Measurement: A 3D Snapshot

Usually, scientists measure these collisions in just one or two ways (like measuring the speed of the ghost or the angle it came from). This paper is the first time anyone has measured the collision in three dimensions at once.

They looked at:

1. How fast the ghost was moving (Kinetic Energy).
2. The angle it came in at (Scattering Angle).
3. How much energy was left over to create other particles (Available Energy).

• The Analogy: Imagine trying to describe a car crash.
  • Old way: "The car was going 60 mph."
  • New way: "The car was going 60 mph, hit the wall at a 45-degree angle, and the bumper flew off 10 feet to the left."
  • This "3D" view lets them see exactly what happened inside the crash, separating different types of collisions that were previously mixed together.

4. The Mystery: The "Ghost" vs. The "Prediction"

The scientists compared their real-world data (the million ghosts they caught) against what their computer models predicted would happen. These computer models are like "recipe books" that tell physicists how the universe should work.

• The Result: The recipe books were wrong.
  • In some areas (low energy), the models predicted too many hits.
  • In other areas (resonance regions), the models predicted the ghosts would bounce off at sharp angles, but in reality, they were more "slippery" and didn't bounce as much as expected.
  • The Analogy: It's like a weather forecast that says, "It will rain hard at 2 PM," but when you look outside, it's actually drizzling at 1 PM and pouring at 3 PM. The total amount of rain might be right, but the timing and intensity are off.

5. Why Does This Matter?

You might ask, "So what? It's just a computer model."

This is crucial because these "ghosts" are used to study neutrino oscillations—a phenomenon where ghosts change their "flavor" (like a chameleon changing color) as they travel. To measure this change accurately, scientists need to know exactly how the ghosts behave when they hit matter.

• The Analogy: If you are trying to measure how fast a runner is going, but you don't know how much the wind is slowing them down, your speed calculation will be wrong.
• The Impact: Because the current computer models (the "wind" models) are inaccurate, future experiments (like the massive DUNE experiment) might get the wrong answers about the fundamental laws of the universe.

The Conclusion

This paper is a massive step forward. By catching a million ghosts and taking a high-definition 3D photo of their collisions, the NOvA team has shown that our current "rulebooks" for how these particles interact are incomplete.

They have handed the physics community a list of errors in the recipe book. Now, theorists have to rewrite the recipes to match the reality the detectors saw. This will make future experiments much more accurate, helping us understand the deep secrets of the universe, from why the universe is made of matter instead of antimatter, to the very nature of time and space.

Technical Summary

Here is a detailed technical summary of the paper "First inclusive triple-differential measurement of the muon-antineutrino charged-current cross section using the NOvA Near Detector."

1. Problem Statement and Motivation

• Context: Long-baseline neutrino oscillation experiments (like NOvA, T2K, and the upcoming DUNE) rely on precise measurements of neutrino oscillation probabilities. These probabilities are inferred by comparing energy spectra at the source and after propagation.
• The Challenge: A leading source of systematic uncertainty in these measurements is the modeling of neutrino-nucleus interactions. Current models struggle to accurately describe the complex nuclear environment, including quasi-elastic (QE), resonant (Res), deep-inelastic scattering (S/DIS), and multi-nucleon effects (2p2h), as well as final-state interactions (FSI).
• Specific Gap: While muon-neutrino ($\nu_\mu$) cross-section measurements are abundant, muon-antineutrino ($\bar{\nu}_\mu$) measurements are scarce. Existing $\bar{\nu}_\mu$ data are often limited to single-differential measurements (e.g., only as a function of energy) or lack the granularity to distinguish between different interaction processes (QE vs. Res vs. DIS) within the same dataset.
• Goal: To provide a high-statistics, triple-differential measurement of the inclusive $\bar{\nu}_\mu$ charged-current (CC) cross section to constrain interaction models and reduce systematic uncertainties for future oscillation experiments.

2. Methodology

• Data Source:
  • Experiment: NOvA Near Detector (ND), located 1 km from the neutrino production target at Fermilab.
  • Beam: NuMI beam, configured for antineutrino mode.
  • Exposure: $12.5 \times 10^{20}$ protons-on-target (POT) collected between June 2016 and July 2019.
  • Sample Size: Approximately 1.04 million selected $\bar{\nu}_\mu$ CC events (the largest published sample of its kind).
• Detector:
  • The ND is a 300-ton fiducial volume detector made of extruded PVC cells filled with scintillator.
  • It is off-axis, peaking the beam energy around 2 GeV.
  • Event selection requires a contained antimuon candidate with a vertex in the fiducial volume.
• Analysis Variables (Triple-Differential):
  The cross section is measured as a function of three kinematic variables:
  1. $T_\mu$: Kinetic energy of the final-state antimuon.
  2. $\cos \theta_\mu$: Cosine of the scattering angle of the antimuon relative to the beam direction.
  3. $E_{\text{avail}}$: Available energy of the interaction (proxy for visible hadronic energy), defined as the sum of neutral pion energies and kinetic energies of protons and charged pions.
• Unfolding and Corrections:
  • Data were corrected for purity (60–90%), efficiency (20–90%), and detector resolution using an iterative D'Agostini unfolding method (2 iterations).
  • Systematic uncertainties were evaluated by comparing the nominal simulation to shifted simulations for flux, detector response, neutrino-nucleus modeling ($\nu$-A), and neutron propagation.
• Simulation and Models:
  • Primary Generator: GENIE 3.0.6 (NOvA Tune v2), which was tuned using neutrino-mode data to better match the Near Detector.
  • Comparisons: Results were compared against multiple GENIE Comprehensive Model Configurations (CMCs) and other major generators: NuWro, NEUT, and GiBUU.

3. Key Contributions

• First Triple-Differential $\bar{\nu}_\mu$ Measurement: This is the first inclusive measurement of the $\bar{\nu}_\mu$ CC cross section in three dimensions ($T_\mu$, $\cos \theta_\mu$, $E_{\text{avail}}$) with 459 bins.
• Highest Statistics: The dataset contains ~1 million events, providing unprecedented statistical power to isolate specific phase-space regions.
• Process Separation: The use of $E_{\text{avail}}$ allows the authors to cleanly separate regions dominated by different interaction mechanisms:
  • $E_{\text{avail}} < 0.1$ GeV: Dominated by QE and 2p2h.
  • $0.1 < E_{\text{avail}} < 0.3$ GeV: Transition region (QE/2p2h to Res).
  • $0.3 < E_{\text{avail}} < 1.0$ GeV: Dominated by Resonance (Res) production.
  • $E_{\text{avail}} > 1.0$ GeV: Dominated by Shallow/Deep Inelastic Scattering (S/DIS).

4. Key Results

• Comparison with NOvA Tune v2:
  • The tuned GENIE model agrees well with data in the low $E_{\text{avail}}$ region (QE/2p2h dominated).
  • Discrepancy: In the Resonance region ($0.3 < E_{\text{avail}} < 1.0$ GeV), the measured cross sections are consistently lower than the GENIE predictions.
  • Agreement: Good agreement is observed in the high $E_{\text{avail}}$ region (S/DIS dominated).
• Comparison with Other Generators:
  • QE/2p2h Region: Untuned simulations (GENIE G18, AR23, NuWro, NEUT) struggle with both shape and normalization. The SuSAv2 model (used in GENIE G21) shows better agreement with data than the Valencia model, suggesting the Valencia model's treatment of QE and 2p2h needs revision.
  • Resonance Region: All generators using the Berger-Sehgal model show shape discrepancies, particularly in the angular distribution ($\cos \theta_\mu$). Most generators predict a sharper forward peak than observed, indicating a need for additional suppression mechanisms in the models.
  • S/DIS Region: GENIE (using Bodek-Yang + PYTHIA6) agrees well. NuWro and NEUT show discrepancies at forward angles, likely due to differences in Final State Interaction (FSI) tuning. GiBUU reproduces the shape at high angles but has normalization issues at forward angles.
• $\chi^2$ Analysis:
  • The NOvA Tune v2 provides the best overall fit ($\chi^2/\text{ndof} \approx 2.9$).
  • Untuned models generally have higher $\chi^2$ values (e.g., NuWro $\approx 15.3$, GENIE G18 $\approx 9.5$), highlighting significant tensions between current models and the data.

5. Significance

• Model Constraints: The results provide strong constraints on neutrino interaction generators, specifically highlighting deficiencies in the modeling of:
  • 2p2h mechanisms: The Valencia model appears less accurate than SuSAv2 for antineutrinos.
  • Resonance production: Current models overpredict the cross section and fail to capture the correct angular dependence.
• Impact on Oscillation Physics: Since next-generation experiments (like DUNE) will operate in the GeV energy range where these interaction uncertainties are dominant, this measurement is critical for reducing systematic errors in oscillation parameter extraction (e.g., $\delta_{CP}$, mass ordering).
• Antineutrino Specifics: By focusing on antineutrinos, this work addresses a critical gap in knowledge, as antineutrino-nucleus interactions differ significantly from neutrino interactions due to nuclear effects and cross-section magnitudes.

In summary, this paper establishes a new benchmark for antineutrino interaction modeling, demonstrating that while current generators can be tuned to fit specific datasets, significant theoretical improvements are required to universally describe the triple-differential kinematics of $\bar{\nu}_\mu$ interactions across all energy regimes.