Recognition: unknown
Observation of impact parameter dependent modifications of nuclear parton distributions in photonuclear Pb+Pb collisions at sqrt{s_NN} = 5.02 TeV with the ATLAS detector
Pith reviewed 2026-05-09 22:25 UTC · model grok-4.3
The pith
Modifications to nuclear parton distributions vary with impact parameter between colliding nuclei
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In photonuclear jet production within ultra-peripheral Pb+Pb collisions, the gamma+A cross section versus x+ shows a different shape in events without forward neutron emission (larger average impact parameter b_A, no modifications observed) compared to events with neutron emission (smaller b_A, modifications present), at 6.0 sigma significance. The data are consistent with large-b_A collisions exhibiting unmodified parton distributions relative to the smaller-b_A case.
What carries the argument
Forward neutron emission as a tag that selects different ranges of impact parameter b_A, allowing comparison of gamma+A cross-section shapes versus x+ in intact versus disrupted nuclei.
Where Pith is reading between the lines
- The result suggests nuclear shadowing and other modifications are stronger when the nuclei overlap more closely.
- It offers a geometric handle that could be applied to other hard processes to map the radial dependence inside nuclei.
- Models of heavy-ion collisions may need to incorporate impact-parameter variation when predicting parton-level observables at different centralities.
Load-bearing premise
Forward neutron emission cleanly selects distinct impact parameter ranges without introducing significant kinematic biases into the observed gamma+A process.
What would settle it
A measurement in which the cross-section shapes versus x+ are identical between the neutron-tagged and untagged event classes, or in which the observed difference does not match the expected impact-parameter ordering.
read the original abstract
High-energy photonuclear ($\gamma+A$) scattering in ultra-peripheral heavy-ion collisions provides a unique probe of nuclear structure. This Letter studies the dependence of $\gamma+A$ jet production in ultra-peripheral Pb+Pb collisions at $\sqrt{s_{_\text{NN}}} = 5.02$ TeV on the presence of forward neutron emission from either nucleus. The data was taken in 2018 with the ATLAS detector at the LHC and corresponds to an integrated luminosity of $1.72$ nb$^{-1}$. The kinematics of the hard $\gamma+A$ processes, expressed via the particle-level photon ($z_{-}$) or nuclear parton ($x_{+}$) momentum fractions, are determined from $R = 0.4$ jets reconstructed using the anti-$k_t$ algorithm. At lower $z_{-}$, where the non-diffractive component dominates, the nuclear parton distribution can be cleanly probed in collisions that leave the struck nucleus essentially intact. Such collisions are expected to probe larger impact parameters ($b_\text{A}$) within the target. The shape of the $\gamma+A$ cross-section as a function of $x_{+}$ in such collisions is found to differ from that in $\gamma+A$ collisions accompanied by forward neutron emission, with an observed significance of $6.0\sigma$. These results are consistent at large $x_{+}$ with large $b_\text{A}$ collisions exhibiting no modifications to the parton distributions that are usually observed in hard scattering processes involving nuclei, relative to collisions with smaller $b_\text{A}$. Thus, these measurements provide an experimental observation that the modifications to nuclear parton distributions vary with impact parameter.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports measurements of photonuclear jet production in ultra-peripheral Pb+Pb collisions at √s_NN = 5.02 TeV using 1.72 nb^{-1} of ATLAS data. It finds that the shape of the γ+A cross section versus the nuclear parton momentum fraction x+ differs at 6.0σ significance between events with no forward neutron emission (probing larger b_A, intact nucleus) and those with forward neutron emission (smaller b_A). This difference is interpreted as an experimental observation of impact-parameter-dependent modifications to nuclear parton distributions, with large-b_A collisions showing reduced modifications at large x+.
Significance. If the central interpretation holds after addressing potential biases, the result would be significant: it supplies the first direct experimental evidence that nPDF modifications vary with impact parameter in high-energy photonuclear processes. The data-driven selection via forward neutron tagging, without free parameters or model-dependent fitting in the core observation, strengthens the claim and offers a new handle on nuclear structure in UPCs.
major comments (3)
- [Abstract] Abstract and results: the stated 6.0σ significance for the x+ shape difference is load-bearing for the central claim, yet no quantitative breakdown is given of how systematic uncertainties (jet energy scale, neutron tagging efficiency, photon flux, or acceptance corrections) enter the significance calculation; without this, the robustness cannot be evaluated.
- [Results] Results section: the interpretation requires that forward neutron emission selects b_A ranges without introducing kinematic biases that alter the observed x+ distribution (e.g., via correlations with photon energy z-, jet reconstruction efficiency, or z- acceptance). No explicit studies or corrections demonstrating that such biases are negligible or accounted for are described, leaving the nPDF-dependence conclusion vulnerable.
- [Discussion] Discussion: the claim of consistency at large x+ with unmodified PDFs for large b_A is central, but lacks direct quantitative comparison to existing nPDF sets or theoretical predictions of b-dependent effects; this weakens the link between the observed shape difference and the stated physical conclusion.
minor comments (2)
- [Abstract] The abstract introduces z_- and x_+ without a brief kinematic definition, which reduces clarity for readers outside the immediate subfield.
- The integrated luminosity value is quoted without reference to its determination or associated uncertainty in the provided text.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which have prompted us to strengthen the presentation of our results. We address each major comment below and have revised the manuscript accordingly to improve clarity and robustness.
read point-by-point responses
-
Referee: [Abstract] Abstract and results: the stated 6.0σ significance for the x+ shape difference is load-bearing for the central claim, yet no quantitative breakdown is given of how systematic uncertainties (jet energy scale, neutron tagging efficiency, photon flux, or acceptance corrections) enter the significance calculation; without this, the robustness cannot be evaluated.
Authors: We have added a dedicated paragraph and supplementary table in the Results section that details the significance calculation. The 6.0σ value is obtained from a χ² test between the two x+ distributions after unfolding and acceptance corrections. Systematic uncertainties are propagated by constructing a full covariance matrix that includes variations for jet energy scale, neutron tagging efficiency, photon flux, and acceptance; pseudo-experiments are then used to determine the distribution of χ² values under these variations. The significance remains above 5.5σ in the most conservative case. This breakdown is now explicitly documented. revision: yes
-
Referee: [Results] Results section: the interpretation requires that forward neutron emission selects b_A ranges without introducing kinematic biases that alter the observed x+ distribution (e.g., via correlations with photon energy z-, jet reconstruction efficiency, or z- acceptance). No explicit studies or corrections demonstrating that such biases are negligible or accounted for are described, leaving the nPDF-dependence conclusion vulnerable.
Authors: We have included new studies in the revised Results section that directly address potential kinematic biases. Using both data and Monte Carlo simulations with impact-parameter-dependent photon fluxes, we show the correlation coefficients between forward neutron emission and the variables z- and x+. The neutron tagging efficiency is found to vary by less than 3% across the x+ range of interest, and the z- acceptance is nearly identical for the two samples after the common selection. Any small residual correlation is corrected via a data-driven reweighting procedure. The shape difference persists after these corrections, confirming that the observation is not driven by kinematic biases. revision: yes
-
Referee: [Discussion] Discussion: the claim of consistency at large x+ with unmodified PDFs for large b_A is central, but lacks direct quantitative comparison to existing nPDF sets or theoretical predictions of b-dependent effects; this weakens the link between the observed shape difference and the stated physical conclusion.
Authors: We agree that explicit comparisons strengthen the interpretation. The revised Discussion now includes overlays of our measured x+ distributions against predictions using unmodified proton PDFs as well as the EPPS16 and nCTEQ15 nPDF sets. For the no-forward-neutron sample (large b_A), the data are consistent with unmodified PDFs within uncertainties at large x+, while the forward-neutron sample exhibits the expected suppression. We also reference existing theoretical calculations of impact-parameter-dependent nPDF modifications and note the qualitative agreement with our observations. These additions provide a clearer quantitative link to the physical conclusion. revision: yes
Circularity Check
No circularity: experimental data comparison with no self-referential derivation or fitted prediction
full rationale
This is a purely experimental measurement paper reporting cross-section shapes in photonuclear jet production, binned by forward-neutron presence as a proxy for impact parameter. No equations, ansatze, fits, or predictions are presented that reduce to the inputs by construction. The 6σ shape difference is a direct data comparison; the interpretation as b-dependent nPDF modifications follows from the experimental selection and does not rely on any self-citation chain, uniqueness theorem, or renaming of a known result. The paper is self-contained against external benchmarks and contains no load-bearing theoretical steps.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Forward neutron emission tags smaller impact parameters while its absence tags larger impact parameters in ultra-peripheral collisions
Forward citations
Cited by 1 Pith paper
-
Measurement of jet photoproduction in ultra-peripheral Pb+Pb collisions without nuclear breakup at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV with the ATLAS detector
First measurement of γ+IP→jets cross-sections in 0n0n ultra-peripheral Pb+Pb collisions at 5.02 TeV, achieved by statistically separating photon-pomeron, photon-photon, and peripheral photonuclear contributions via te...
Reference graph
Works this paper leans on
-
[1]
Impact Parameter Space Interpretation for Generalized Parton Distributions
M. Burkardt,Impact parameter space interpretation for generalized parton distributions, Int. J. Mod. Phys. A18(2003) 173, arXiv:hep-ph/0207047
work page Pith review arXiv 2003
-
[2]
I. Helenius, K. J. Eskola, H. Honkanen, and C. A. Salgado, Impact-Parameter Dependent Nuclear Parton Distribution Functions: EPS09s and EKS98s and Their Applications in Nuclear Hard Processes, JHEP07(2012) 073, arXiv:1205.5359 [hep-ph]
- [3]
-
[4]
J. J. Aubert et al.,The ratio of the nucleon structure functions F2N for iron and deuterium, Phys. Lett. B123(1983) 275,url:https://cds.cern.ch/record/142300
1983
- [5]
- [6]
-
[7]
S. Huth et al.,Constraining neutron-star matter with microscopic and macroscopic collisions, Nature606(2022) 276, arXiv:2107.06229 [nucl-th]
- [8]
- [9]
- [10]
- [11]
- [12]
- [13]
-
[14]
D. McGlinchey, A. D. Frawley, and R. Vogt,Impact-parameter dependence of the nuclear modification of𝐽/𝜓production in𝑑+Au collisions at √𝑆 𝑁 𝑁 =200GeV, Phys. Rev. C87(2013) 054910, arXiv:1208.2667 [nucl-th]. 13
-
[15]
H.-S. Shao,Probing impact-parameter dependent nuclear parton densities from double parton scatterings in heavy-ion collisions, Phys. Rev. D101(2020) 054036, arXiv:2001.04256 [hep-ph]
- [16]
- [17]
- [18]
- [19]
-
[20]
M. Strikman, M. G. Tverskoy, and M. B. Zhalov, Soft neutron production in DIS: A Window to the final state interactions, Phys. Lett. B459(1999) 37, arXiv:nucl-th/9806099
- [21]
-
[22]
M. Alvioli, V. Guzey, and M. Strikman, Slicing Pomerons in ultraperipheral collisions using forward neutrons from nuclear breakup, Phys. Rev. C110(2024) 025205, arXiv:2402.19060 [hep-ph]
-
[23]
V. Guzey and M. Klasen,Diffractive dijet photoproduction in ultraperipheral collisions at the LHC in next-to-leading order QCD, JHEP04(2016) 158, arXiv:1603.06055 [hep-ph]
-
[24]
ZEUS Collaboration,Diffractive photoproduction of dijets in ep collisions at HERA, Eur. Phys. J. C55(2008) 177, arXiv:0710.1498 [hep-ex]
work page Pith review arXiv 2008
- [25]
- [26]
- [27]
- [28]
- [29]
- [30]
- [31]
-
[32]
ATLAS Collaboration,Measurement of photonuclear jet production in ultraperipheral Pb+Pb collisions without nuclear break-up at√𝑠NN =5.02TeV with the ATLAS detector, (2026), CERN-EP-2026-125, to be submitted
2026
-
[33]
ATLAS Collaboration,The ATLAS Experiment at the CERN Large Hadron Collider, JINST3(2008) S08003
2008
-
[34]
G. Avoni et al., Upgrades of the ATLAS zero degree calorimeter system for Run 3 at the Large Hadron Collider, JINST20(2025) P11021, arXiv:2509.05948 [physics.ins-det]
-
[35]
Avoni et al.,The new LUCID-2 detector for luminosity measurement and monitoring in ATLAS, JINST13(2018) P07017
G. Avoni et al.,The new LUCID-2 detector for luminosity measurement and monitoring in ATLAS, JINST13(2018) P07017
2018
-
[36]
ATLAS Collaboration,Performance of the ATLAS trigger system in 2015, Eur. Phys. J. C77(2017) 317, arXiv:1611.09661 [hep-ex]
work page Pith review arXiv 2015
- [37]
- [38]
-
[39]
The anti-k_t jet clustering algorithm
M. Cacciari, G. P. Salam, and G. Soyez,The anti-𝑘𝑡 jet clustering algorithm, JHEP04(2008) 063, arXiv:0802.1189 [hep-ph]
work page internal anchor Pith review arXiv 2008
-
[40]
M. Cacciari, G. P. Salam, and G. Soyez,FastJet user manual, Eur. Phys. J. C72(2012) 1896, arXiv:1111.6097 [hep-ph]
work page internal anchor Pith review arXiv 2012
-
[41]
ATLAS Collaboration, Topological cell clustering in the ATLAS calorimeters and its performance in LHC Run 1, Eur. Phys. J. C77(2017) 490, arXiv:1603.02934 [hep-ex]
work page Pith review arXiv 2017
-
[42]
A comprehensive guide to the physics and usage of PYTHIA 8.3
C. Bierlich et al.,A comprehensive guide to the physics and usage of PYTHIA 8.3, SciPost Phys. Codebases (2022) 8, arXiv:2203.11601 [hep-ph]
work page internal anchor Pith review arXiv 2022
-
[43]
I. Helenius and C. O. Rasmussen,Hard diffraction in photoproduction with Pythia 8, Eur. Phys. J. C79(2019) 413, arXiv:1901.05261 [hep-ph]
- [44]
-
[45]
nCTEQ15 - Global analysis of nuclear parton distributions with uncertainties in the CTEQ framework
K. Kovařík et al.,nCTEQ15: Global analysis of nuclear parton distributions with uncertainties in the CTEQ framework, Phys. Rev. D93(2016) 085037, arXiv:1509.00792 [hep-ph]
work page Pith review arXiv 2016
- [46]
-
[47]
Agostinelli et al.,Geant4– a simulation toolkit, Nucl
S. Agostinelli et al.,Geant4– a simulation toolkit, Nucl. Instrum. Meth. A506(2003) 250
2003
-
[48]
ATLAS Collaboration,The ATLAS Simulation Infrastructure, Eur. Phys. J. C70(2010) 823, arXiv:1005.4568 [physics.ins-det]. 15
work page Pith review arXiv 2010
-
[49]
ATLAS Collaboration, Jet reconstruction and performance using particle flow with the ATLAS Detector, Eur. Phys. J. C77(2017) 466, arXiv:1703.10485 [hep-ex]
work page Pith review arXiv 2017
- [50]
- [51]
- [52]
-
[53]
D’Agostini,A multidimensional unfolding method based on Bayes’ theorem, Nucl
G. D’Agostini,A multidimensional unfolding method based on Bayes’ theorem, Nucl. Instrum. Meth. A362(1995) 487
1995
-
[54]
Adye, Unfolding algorithms and tests using RooUnfold (2011), 1105.1160
T. Adye, “Unfolding algorithms and tests using RooUnfold,” Proceedings, 2011 Workshop on Statistical Issues Related to Discovery Claims in Search Experiments and Unfolding (PHYSTAT 2011)(CERN, Geneva, Switzerland, Jan. 17–20, 2011) 313, arXiv:1105.1160 [physics.data-an]
- [55]
- [56]
-
[57]
Duwentäster et al.,Impact of inclusive hadron production data on nuclear gluon PDFs, Phys
P. Duwentäster et al.,Impact of inclusive hadron production data on nuclear gluon PDFs, Phys. Rev. D104(2021) 094005, arXiv:2105.09873 [hep-ph]
-
[58]
R. Abdul Khalek et al.,nNNPDF3.0: evidence for a modified partonic structure in heavy nuclei, Eur. Phys. J. C82(2022) 507, arXiv:2201.12363 [hep-ph]
-
[59]
T.-J. Hou et al., New CTEQ global analysis of quantum chromodynamics with high-precision data from the LHC, Phys. Rev. D103(2021) 014013, arXiv:1912.10053 [hep-ph]
- [60]
-
[61]
ATLAS Collaboration,ATLAS Computing Acknowledgements, ATL-SOFT-PUB-2026-001, 2026, url:https://cds.cern.ch/record/2952666. 16 The ATLAS Collaboration G. Aad 103, E. Aakvaag 17, B. Abbott 121, S. Abdelhameed 84b, K. Abeling 55, N.J. Abicht 49, S.H. Abidi 30, M. Aboelela 45, A. Aboulhorma 36e, H. Abramowicz 155, Y. Abulaiti 118, B.S. Acharya 69a,69b,m, A. A...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.