Single pion-production and pion propagation in Achilles

We extend the applicability of Achilles (A CHIcagoLand Lepton Event Simulator) by incorporating the single-pion production mechanism in a fully exclusive fashion. The electroweak interaction vertex is modeled by combining the state-of-the-art Dynamical Coupled-Channels approach with realistic hole spectral functions, which account for correlations in both the initial target state and the residual spectator system. Final-state interactions are treated using a semi-classical intranuclear cascade that leverages nuclear configurations to determine the correlated spatial distribution of protons and neutrons. The meson-baryon scattering amplitudes used in the cascade are computed within the Dynamical Coupled-Channels framework, consistent with the electroweak vertex. To model pion absorption, we employ the optical potential approach of Oset and Salcedo. As an alternative approach, we explicitly model the production and propagation of resonances which mediate pion-nucleon scattering and pion absorption. We validate out approach against pion-nucleon and pion-nucleus scattering data, and present comparisons with electron- and neutrino-nucleus measurements from e4$ν$, T2K, MINER$ν$A, and MicroBooNE.

💡 Research Summary

This paper presents a comprehensive extension of the Achilles (A CHIcagoLand Lepton Event Simulator) framework to include fully exclusive single‑pion production and propagation in nuclei. The authors combine the state‑of‑the‑art Dynamical Coupled‑Channels (DCC) model of the ANL‑Osaka collaboration with realistic hole spectral functions (SF) for 12C and 40Ar, thereby providing a quantum‑mechanical description of the electroweak interaction vertex that incorporates both resonant and non‑resonant contributions in a unitary fashion.

The DCC approach solves coupled integral equations for a set of meson‑baryon channels (πN, ηN, KΛ, …) using an energy‑independent Hamiltonian H = H₀ + Γ + v, where Γ generates bare N* and Δ resonances and v encodes non‑resonant meson exchange and Z‑diagram mechanisms. By projecting onto the space of stable particles, the model yields a T‑matrix that is split into a background term t_bg and a resonance term t_res. The electroweak current is built from the same ingredients, guaranteeing consistency between the strong and weak sectors. Vector transition form factors are constrained by electron‑pion production data, while axial form factors follow the partially conserved axial current (PCAC) hypothesis and are tuned to historic ν‑deuterium measurements.

To capture nuclear initial‑state correlations, the authors employ hole spectral functions that give the joint probability of removing a nucleon with momentum k and removal energy E. These SFs are normalized separately for protons and neutrons and provide a realistic distribution of initial nucleon momenta beyond the simple Fermi‑gas picture. The factorization ansatz then allows the hard interaction (described by the DCC amplitudes) to be convoluted with the SF, producing exclusive final states containing a nucleon and a pion (1N1π).

Final‑state interactions (FSI) are modeled with a semi‑classical intranuclear cascade (INC) that propagates all hadrons through a realistic nuclear configuration. Two distinct cascade modes are implemented to explore the sensitivity to pion absorption mechanisms:

1. Optical‑potential mode – The well‑known Oset‑Salcedo optical potential provides a complex self‑energy for the pion; its imaginary part accounts for absorption. In this scheme, Δ and higher resonances are not propagated as explicit degrees of freedom; absorption occurs instantaneously via the optical potential.

2. Resonance‑propagation mode – Inspired by GiBUU and INCL, the Δ (and selected N* states) are treated as dynamical particles that travel through the nuclear medium, decay, and undergo three‑body processes such as π N ↔ Δ, Δ N ↔ NN. This approach captures rescattering before decay, medium‑modified widths, and the long Δ lifetime relative to typical nuclear transit times, offering a more microscopic description of pion absorption and charge‑exchange.

Both cascade modes use meson‑baryon scattering amplitudes computed directly from the DCC model, ensuring that the same underlying physics governs both the production vertex and the subsequent propagation.

The authors validate the extended Achilles in several steps. First, inclusive electron‑nucleus cross sections are reproduced, confirming that the combination of DCC amplitudes and spectral functions respects energy–momentum conservation and yields realistic strength distributions. Next, pion‑nucleus total and differential cross sections are compared to data, showing that the optical‑potential mode reproduces overall absorption rates while the resonance‑propagation mode better captures angular distributions, especially at intermediate angles where rescattering is important.

Finally, the model is confronted with exclusive electron‑ and neutrino‑nucleus measurements from the e4ν, T2K, MINERνA, and MicroBooNE experiments. Comparisons focus on observables sensitive to pion final‑state interactions: pion kinetic‑energy spectra, multiplicities, and lepton–hadron correlation variables. The resonance‑propagation cascade generally yields improved agreement with data, reducing the discrepancy in pion‑absorption fractions and providing a more accurate description of events that would otherwise be mis‑identified as quasielastic.

Key insights from the work include:

• The DCC model provides a unified, unitary treatment of resonant and non‑resonant pion production, preserving interference effects across the entire energy range relevant for accelerator‑based neutrino experiments.
• Realistic hole spectral functions introduce nuclear‑correlation effects that are essential for reproducing the observed lepton kinematics and for correctly modeling the missing‑energy tail.
• Implementing two distinct FSI schemes allows a quantitative assessment of model uncertainties associated with pion absorption; the resonance‑propagation approach, while computationally more demanding, offers a more faithful representation of medium‑modified Δ dynamics.
• The extended Achilles demonstrates superior performance relative to widely used generators (NEUT, GENIE) in reproducing exclusive pion‑production data, suggesting that the approach can substantially reduce systematic uncertainties in oscillation analyses for DUNE, Hyper‑K, and the SBN program.

The paper concludes by outlining future directions: incorporation of multi‑pion and 2p‑2h channels, refinement of axial form factors using modern lattice QCD inputs, and systematic propagation of model uncertainties into oscillation fits. The authors argue that the presented framework constitutes a significant step toward the high‑precision neutrino‑interaction modeling required for next‑generation long‑baseline experiments.