Theoretical modeling of charged current $\nu_\mu(\bar\nu_\mu)-^{40}Ar$ DIS at DUNE energies

This paper presents a theoretical study of charged current $\nu_\mu(\bar{\nu}_\mu)$-induced deep inelastic scattering on $^{40}\mathrm{Ar}$ at DUNE-relevant energies, utilizing a microscopic framework with nuclear medium effects and NNLO QCD corrections to compute nuclear structure functions and differential cross sections.

The Gist

Imagine you are trying to listen to a specific conversation in a very noisy, crowded room. That's essentially what physicists are doing when they study neutrinos—tiny, ghost-like particles that rarely interact with anything.

This paper is about improving the "listening equipment" for a massive experiment called DUNE (Deep Underground Neutrino Experiment). The goal is to understand the universe's fundamental secrets, like why there is more matter than antimatter. To do this, they shoot a beam of neutrinos at a giant tank of liquid Argon (a noble gas, like the one in lightbulbs, but in liquid form).

Here is a breakdown of what the authors did, using simple analogies:

1. The Problem: The "Crowded Room" Effect

When a neutrino hits a single, lonely atom floating in space, the physics is relatively straightforward. It's like a cue ball hitting a single billiard ball.

However, in DUNE, the neutrino hits an Argon nucleus, which is a tight cluster of 40 protons and neutrons packed together. It's not a cue ball hitting one ball; it's a cue ball hitting a stack of 40 marbles glued together.

• The Issue: The marbles inside the stack are jiggling (Fermi motion), stuck together with glue (binding energy), and constantly bumping into each other (correlations).
• The Result: If you use the "single ball" math to predict what happens in the "stack of 40," your predictions will be wrong. This creates "systematic errors" that could ruin the experiment's ability to measure the universe's secrets.

2. The Solution: A Better "Recipe"

The authors created a new, more sophisticated mathematical model (a "recipe") to predict exactly how neutrinos interact with this Argon stack. They didn't just guess; they built a microscopic simulation that accounts for three main things:

• The "Jiggling" (Spectral Function): They accounted for the fact that the particles inside the Argon aren't sitting still; they are zooming around.
• The "Cloud" (Mesonic Contributions): Imagine the nucleons (protons/neutrons) are surrounded by a fuzzy cloud of virtual particles (pions and rho mesons). Sometimes, the neutrino hits the cloud instead of the core particle. The authors added this "cloud" effect into their math, finding it makes a huge difference, especially at lower energies.
• The "Shadow" (Shadowing/Antishadowing): When you shine a light through a thick forest, the trees block some light (shadowing). In the nucleus, the particles can "shadow" each other from the neutrino's view. The authors included these effects too.

3. The "Filter" Problem (The W-Cut)

In physics, there is a "safe zone" where the math is easy to do (called Deep Inelastic Scattering or DIS). But there is a "messy zone" in the middle where particles act like resonances (vibrating like a guitar string) rather than free particles.

• The Analogy: Imagine trying to sort red and blue marbles. The "messy zone" is where the marbles are purple (a mix of red and blue).
• The Fix: Physicists often put a "filter" on their data, saying, "We will only count the events where the energy is high enough to be definitely 'Deep Inelastic'." This is called a W-cut.
• The Finding: The authors found that applying this filter is like using a very strict sieve. It throws away a massive amount of data, especially for antineutrinos (the anti-matter version of neutrinos).
  • For neutrinos, the filter cuts out about 27% of the data.
  • For antineutrinos, the filter cuts out a whopping 70-94% of the data!
  • Why this matters: If you throw away 90% of your data, your experiment becomes much less precise. The authors are warning that we need better ways to understand the "messy zone" so we don't have to throw away so much data.

4. The Main Takeaways

• Nuclear effects are huge: You cannot ignore the fact that the target is a nucleus. It changes the results significantly (sometimes by 20-30%).
• Antineutrinos are trickier: The "messy zone" affects antineutrinos much more than neutrinos. If you don't model this correctly, your measurements of the universe's secrets will be off.
• The "Cloud" matters: The virtual particle clouds around the nucleons add a significant amount of "extra" interaction that previous models might have missed.
• Better Math = Better Science: By using a more complex model that includes all these nuclear effects, the authors provide a more accurate map for the DUNE experiment. This helps the scientists reconstruct the energy of the neutrinos correctly, which is the key to solving the mystery of the neutrino mass hierarchy and CP violation.

In a Nutshell

Think of this paper as the team that built a high-definition 3D map of a crowded city (the Argon nucleus) instead of a flat, 2D sketch. They realized that if you use the flat sketch to navigate the city, you'll get lost. Their new map accounts for the traffic, the buildings, and the shadows, ensuring that the DUNE experiment doesn't get lost while trying to find the secrets of the universe.

Technical Summary

1. Problem Statement

The paper addresses the critical need for precise theoretical modeling of (anti)neutrino-nucleus interaction cross-sections, specifically for Deep Inelastic Scattering (DIS) on Argon-40 ($^{40}\text{Ar}$) targets. This is vital for current and future neutrino oscillation experiments like DUNE (Deep Underground Neutrino Experiment) and the Short-Baseline Neutrino (SBN) program.

• The Challenge: In the few-GeV energy region (relevant for DUNE, $\langle E \rangle \approx 2.5-3$ GeV), the transition between the resonance region (where nucleons are excited) and the DIS region (where partons are probed) is complex. Standard perturbative QCD (pQCD) is not strictly applicable in the "shallow inelastic scattering" (SIS) region ($W \lesssim 2$ GeV, $Q^2 \lesssim 2$ GeV$^2$).
• Uncertainties: Current event generators (e.g., GENIE) rely on phenomenological models with significant limitations in this transition region. Insufficient understanding of nuclear medium effects (Fermi motion, binding energy, correlations, mesonic contributions, shadowing) leads to $\sim 25\%$ uncertainties in cross-sections, which directly impacts the precision of oscillation parameter extraction (e.g., $\delta_{CP}$, mass hierarchy).
• Specific Gap: There is a lack of a comprehensive microscopic framework that simultaneously incorporates higher-order QCD corrections, target mass corrections, and specific nuclear medium effects for Argon targets in the weak sector (charged current).

2. Methodology

The authors employ a microscopic many-body field-theoretical approach to model the charged current DIS process.

• Theoretical Framework:

  • Free Nucleon Baseline: The calculation starts with free nucleon structure functions ($F_{iN}$) derived from Parton Distribution Functions (PDFs). They utilize the MMHT 2014 parameterization, including corrections up to Next-to-Next-to-Leading Order (NNLO) in perturbative QCD and Target Mass Corrections (TMC).
  • Nuclear Medium Effects: The transition from free nucleons to the $^{40}\text{Ar}$ nucleus is handled via a relativistic nucleon spectral function ($S_h$) within the Local Density Approximation (LDA). This accounts for:
    • Fermi Motion: The momentum distribution of nucleons.
    • Binding Energy: The energy required to remove a nucleon.
    • Nucleon Correlations: Short-range correlations between nucleons.
  • Mesonic Contributions: The model includes contributions from virtual meson clouds ($\pi$ and $\rho$) surrounding the nucleons, calculated using a many-body approach.
  • Shadowing/Antishadowing: Corrections for coherent nuclear phenomena at low and intermediate Bjorken-$x$ are incorporated following the Kulagin-Petti formalism (based on Glauber-Gribov theory).
  • Structure Functions: The nuclear structure functions ($F_{iA}$) are constructed by convoluting the free nucleon functions with the spectral function and adding mesonic and shadowing terms.
    • $F_{1A}, F_{2A}$: Include nucleon spectral, mesonic ($\pi, \rho$), and shadowing effects.
    • $F_{3A}$: Includes nucleon spectral and shadowing effects (no mesonic contribution due to valence quark averaging).

• Kinematic Constraints:

  • Calculations are performed for beam energies $E = 4$ GeV and $E = 6$ GeV.
  • The study analyzes differential cross-sections ($d^2\sigma/dxdy$ and $d\sigma/dx$) under different kinematic cuts on the invariant hadronic mass $W$:
    1. No cut (Full kinematic range).
    2. $W \ge 1.7$ GeV.
    3. $W \ge 2$ GeV (Standard "clean" DIS cut).

3. Key Contributions

• Unified Microscopic Model: Development of a consistent theoretical framework that integrates NNLO pQCD, TMC, and a rigorous treatment of nuclear medium effects (spectral function + mesons + shadowing) specifically for $^{40}\text{Ar}$.
• Quantification of Nuclear Effects: A detailed breakdown of how individual nuclear effects (spectral function vs. mesons vs. shadowing) modify the cross-sections relative to free nucleons.
• Impact of $W$ Cuts: A systematic analysis of how imposing kinematic cuts on $W$ (to exclude resonance regions) affects the DIS cross-sections, highlighting the significant loss of events and the shift in kinematic distributions.
• Comparison with Phenomenology: Direct comparison of their microscopic results with global fits of Nuclear PDFs (nCTEQ15, nCTEQnu), revealing discrepancies that suggest current global fits may not fully capture weak-sector nuclear effects in the transition region.

4. Key Results

• Nuclear Suppression (Spectral Function): The inclusion of the nucleon spectral function (Fermi motion, binding, correlations) leads to a suppression of the differential cross-sections compared to free nucleons.
  • At $E=4$ GeV, suppression is $\sim 10-13\%$ at $y=0.5$ for $0.35 \le x \le 0.75$.
  • This suppression is kinematically dependent and non-uniform across $x$ and $y$.
• Mesonic Enhancement: The inclusion of mesonic clouds ($\pi, \rho$) significantly enhances the cross-section, particularly at low-to-intermediate $x$ ($x \lesssim 0.45$).
  • For $\nu_\mu$, enhancement is $\sim 20\%$ at $y=0.5, x=0.35$.
  • For $\bar{\nu}_\mu$, the enhancement is even more pronounced ($\sim 38-48\%$ in certain kinematic regions), indicating a stronger sensitivity of antineutrinos to mesonic degrees of freedom.
• Shadowing/Antishadowing: These effects were found to be negligible in the specific kinematic range ($Q^2 \ge 1$ GeV$^2$, few GeV energies) studied.
• Effect of $W$ Cuts:
  • Applying a cut of $W \ge 2$ GeV causes a drastic suppression of the cross-section, particularly at higher $x$ and lower $y$.
  • For antineutrinos, the suppression is more severe than for neutrinos (e.g., up to $\sim 94\%$ suppression at $x=0.4, E=4$ GeV).
  • The mesonic contribution becomes negligible for $\nu_\mu$ when $W \ge 2$ GeV is applied, but remains significant for $\bar{\nu}_\mu$.
• Comparison with nCTEQ: The authors' "Total" model yields higher cross-sections than the nCTEQnu parameterization at low-to-intermediate $x$ and lower $y$, but lower at high $x$. This suggests that global fits relying heavily on charged-lepton data may not accurately describe weak interactions in the transition region.

5. Significance

• Precision for DUNE: The results underscore that nuclear medium effects are not negligible in the DUNE energy regime. Ignoring them or using simplified models introduces systematic errors that could compromise the measurement of $\delta_{CP}$ and the neutrino mass hierarchy.
• Event Reconstruction: The study highlights that the "safe" DIS region ($W \ge 2$ GeV) contains significantly fewer events than the full kinematic range, especially for antineutrinos. This implies that oscillation analyses must carefully account for the transition region to avoid double-counting or mis-modeling events.
• Theoretical Guidance: The work provides a benchmark for improving neutrino event generators (like GENIE) by demonstrating the necessity of including microscopic spectral functions and mesonic contributions rather than relying solely on phenomenological corrections.
• Weak vs. Electromagnetic Sector: The findings suggest that nuclear modifications in the weak sector (neutrino scattering) may differ from the electromagnetic sector (charged-lepton scattering), challenging the assumption that nuclear PDFs derived from charged-lepton data are directly transferable to neutrino physics without dedicated weak-sector constraints.