First Study of the Nuclear Response to Fast Hadrons via Angular Correlations between Pions and Slow Protons in Electron-Nucleus Scattering

We report on the first measurement of angular correlations between high-energy pions and slow protons in electron-nucleus ($eA$) scattering, providing a new probe of how a nucleus responds to a fast-moving quark. The experiment employed the CLAS detector with a 5-GeV electron beam incident on deuterium, carbon, iron, and lead targets. For heavier nuclei, the pion-proton correlation function is more spread-out in azimuth than for lighter ones, and this effect is more pronounced in the $πp$ channel than in earlier $ππ$ studies. The proton-to-pion yield ratio likewise rises with nuclear mass, although the increase appears to saturate for the heaviest targets. These trends are qualitatively reproduced by state-of-the-art $eA$ event generators, including BeAGLE, eHIJING, and GiBUU, indicating that current descriptions of target fragmentation rest on sound theoretical footing. At the same time, the precision of our data exposes model-dependent discrepancies, delineating a clear path for future improvements in the treatment of cold-nuclear matter effects in $eA$ scattering.

💡 Research Summary

This paper presents the first measurement of angular correlations between high‑energy charged pions (π⁺) and slow (“grey”) protons produced in deep‑inelastic electron‑nucleus (eA) scattering. The experiment was carried out with the CLAS detector at Jefferson Lab using a 5 GeV continuous electron beam incident on four targets: liquid deuterium, carbon, iron, and lead. Integrated luminosities of 1.0 × 10⁴ pb⁻¹ (D), 4.3 × 10² pb⁻¹ (C), 1.4 × 10² pb⁻¹ (Fe), and 1.7 × 10¹ pb⁻¹ (Pb) were accumulated.

Event selection followed standard DIS criteria (Q² > 1 GeV², W > 2 GeV, 2.3 GeV < ν < 4.2 GeV) and required a leading pion with fractional energy z₁ > 0.5, momentum p > 2.7 GeV/c, and Cherenkov‑counter confirmation. Protons were required to have p > 350 MeV/c, transverse momentum p_T > 70 MeV/c, and momentum up to ≈2.7 GeV/c. Both particles had to satisfy polar‑angle fiducial cuts and have azimuthal angles reconstructed with ≤4 mrad resolution.

The two‑dimensional correlation function C(Δϕ,ΔY) was defined as the normalized yield of π⁺–proton pairs per leading‑pion event, where Δϕ = ϕ_π − ϕ_p is the azimuthal separation around the virtual‑photon direction and ΔY = Y_π − Y_p is the rapidity difference. Acceptance corrections were derived from a GEANT‑based CLAS simulation (GSIM) tuned with PYTHIA6, and a data‑driven method removed contamination from the aluminum end‑caps of the deuterium cell. The overall normalization constant C₀ was chosen so that the integral of C over Δϕ and ΔY equals unity for deuterium; the same C₀ was used for all targets, making the nuclear‑to‑deuterium ratio R = C_A/C_D independent of single‑proton efficiency.

From C(Δϕ,ΔY) the azimuthal root‑mean‑square width σ was extracted for each target, and a broadening parameter b = ±√(σ_A² − σ_C²) was defined using carbon as the reference. The results show a clear A‑dependence: σ grows from ≈0.42 rad for carbon to ≈0.55 rad for iron and ≈0.61 rad for lead, indicating that the π⁺–proton correlation becomes more diffuse in heavier nuclei. Correspondingly, the broadening b is positive and increases with A (≈+0.13 rad for Fe, +0.19 rad for Pb).

The proton‑to‑pion yield ratio (π/p) also rises with nuclear mass, reflecting enhanced production of low‑energy protons from nuclear fragmentation. However, the increase appears to saturate for lead, suggesting a limit to the number of “grey” protons that can be knocked out in very heavy nuclei.

Three state‑of‑the‑art eA event generators—BeAGLE, eHIJING, and GiBUU—were run through the same analysis chain. All three reproduce the qualitative trends in σ, b, and π/p, but quantitative discrepancies of 5–15 % remain, especially for lead where the simulated azimuthal width is narrower than observed and the π/p ratio is slightly off. These differences point to incomplete modeling of intra‑nuclear color propagation, proton re‑absorption, and low‑energy fragmentation dynamics.

Systematic uncertainties were evaluated by varying acceptance cuts, particle‑identification thresholds, and luminosity inputs, and by comparing the GEANT‑based simulation to data‑driven corrections. The dominant systematic contributions amount to 5–7 % of the measured observables, comparable to the statistical uncertainties.

The paper discusses the broader impact of these findings. Correlating a fast hadron (π⁺) with a slow proton provides a novel probe of the nuclear response to a fast‑moving quark, accessing lower energy scales of hadronization than dipion studies. The results are directly relevant for future Electron‑Ion Collider (EIC) programs, where dedicated “grey‑proton” detectors will be used for centrality determination and for tagging spectator nucleons in e‑d collisions. Moreover, the measured π‑p correlations can serve as benchmarks for neutrino‑nucleus DIS experiments (e.g., MINERvA, NOvA), where similar kinematics and nuclear effects are encountered.

In summary, this work delivers the first quantitative measurement of π⁺–slow‑proton angular correlations in eA DIS, demonstrates a clear nuclear‑mass dependence of the correlation width and the π/p yield ratio, validates current eA Monte‑Carlo generators while exposing their limitations, and establishes a new experimental observable for studying cold‑nuclear‑matter effects, hadronization, and nuclear fragmentation in both electron‑ and neutrino‑induced reactions.