REVIEW 4 major objections 5 minor 2 cited by
Impacts of isolated nucleon-nucleon correlations in relativistic $^{16}$O+$^{16}$O collisions
T0 review · 4 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read The paper claims that two-nucleon distance correlations, combined with a fixed radial density, quantitatively reproduce the initial-state observables of ab-initio 16O+16O collisions to within a few percent.
desk verdict A promising but underspecified method for injecting ab-initio two-body correlations into initial conditions; the missing one-body density check undermines the central isolation claim. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the acceptance-rejection method (ARM) applied to the two-body distance correlation C(Δr), defined as C(Δr)=1−g(Δr)/g'(Δr), where g is the correlated pair density and g' the uncorrelated pair density. ARM generates candidate nucleon positions from a 3-parameter Fermi radial density with random angles, then accepts or rejects each candidate based on the ratio of the target pair-distance distribution g(Δr) to the proposal distribution g'(Δr) times a constant M. The method's work is to imprint the model's two-body distance correlations onto an otherwise featureless radial profile, so that the differences between configurations are angular in origin.
What would settle it
Compute the one-body radial density ρ(r) of the ARM-accepted configurations and compare it with the input 3pF distribution; if it differs significantly, the claimed isolation of angular correlations fails.
Extended reading notes
Core claim
The central claim is that the acceptance-rejection method (ARM), when seeded with a fixed 3-parameter Fermi radial density and the two-nucleon distance correlation C(Δr) computed from an ab-initio model, produces nucleon configurations whose initial-state observables match those of the original ab-initio configurations to within a few percent. Specifically, the ARM-generated configurations reproduce the eccentricities ε2{2}, ε3{2}, ε4{2}, the energy density variance ⟨δE²⟩, the covariance cov(εn², δE), and the eccentricity fluctuation ratio ε2{4}/ε2{2} for 16O+16O collisions at 200 GeV and 5.36 TeV. The method therefore isolates the angular, two-body part of nuclear structure: differences bet
Load-bearing premise
The load-bearing premise is that the acceptance-rejection sampling leaves the one-body radial density unchanged from the input 3pF profile, a point the paper does not explicitly verify.
Editorial extensions
If this is right
- The initial-state geometry of 16O+16O collisions—eccentricities, energy variance, and flow-correlation ratios—can be captured without full ab-initio wave functions, using only one-body radial density and the two-body distance distribution.
- Structural differences between NLEFT and VMC configurations are mostly angular in origin; the fixed 3pF radial profile common to both is not enough to distinguish them.
- Short-range correlations (Δr ≲ 1 fm), which differ sharply between NLEFT and VMC, are the main source of the deviations in ε2{2} from the 3pF baseline at central collisions.
- Extending the same sampling to three- or four-body distance distributions should recover the remaining few-percent discrepancy between ARM and the original ab-initio configurations.
- The method transfers directly to other light nuclei such as 8Be, 12C, and 20Ne, as the paper states.
Reading between the lines
- If the radial profile is truly preserved, the few-percent residual between ARM and the full ab-initio configurations gives a quantitative upper bound on the information carried by three- and four-body correlations—about 1–5% of the initial-state correlators.
- The same acceptance-rejection pipeline could be inverted: given measured flow correlators in 16O+16O runs at RHIC or the LHC, one could scan over C(Δr) and find the pair-distance distribution the data prefer, turning the method into a model-to-data inversion tool.
- The collision-energy dependence of the OC-to-ARM ratios (slightly larger deviations at 5.36 TeV) hints that the overlap geometry's sensitivity to radial versus angular structure changes with energy, which could be studied systematically by varying the impact-parameter selection or the entropy deposition parameter p.
- Because C(Δr) is derived from a specific Hamiltonian, the method offers a route to compare Hamiltonians directly: two ab-initio models that produce the same C(Δr) would be indistinguishable in all the correlators tested here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces an acceptance-rejection method (ARM) to reconstruct nucleon configurations for 16O from a fixed 3-parameter Fermi radial density and a two-body relative-distance distribution C(Δr) extracted from NLEFT and VMC ab-initio configurations. The resulting configurations are used as initial states in TReNTo for 16O+16O collisions at 200 GeV and 5.36 TeV, and compared to the original configurations through initial-state correlators: εn{2}, var(E), cov(εn^2, δE), and ε2{4}/ε2{2}. The central claim is that ARM quantitatively decodes nucleon-nucleon correlations from ab-initio models, and that differences between configurations can be attributed to angular components because the radial density is fixed to 3pF.
Significance. If the central claim holds, the paper offers a computationally light way to translate two-body correlation information from ab-initio nuclear structure calculations into initial-state observables for relativistic heavy-ion collisions. The non-trivial part is not fitting the correlators—these are genuine predictions once C(Δr) is chosen—and the paper tests a broad set of observables that go beyond the fitted two-body input. The TReNTo simulations are standard and the choice p=0 is clearly stated. The method is simple and applicable to other light nuclei. However, the paper's central interpretation depends on an unverified assumption that the rejection step preserves the one-body radial density; this is a concrete, fixable technical issue rather than a conceptual invalidation. The paper also contains a number of presentation problems that should be corrected.
major comments (4)
- [Section II, Eq. (4) and the acceptance step] The paper never verifies that the acceptance-rejection procedure preserves the 3pF one-body density. The candidate radius r is drawn from 3pF, but acceptance depends on the relative distance Δr to an already placed nucleon, so the marginal density of accepted radii is p_ARM(r) = Z^{-1} p_3pF(r) A(r), where A(r) is the angular average of the acceptance probability over partner positions. Since the target ratio g/g' in Eq. (3) depends on Δr, A(r) is generally not constant. This is load-bearing: Section III and the Conclusion attribute OC vs. ARM differences to angular information, and the claim that 'we start with the same radial density... any differences... are angular' requires this density preservation. Fig. 1(a) shows only the original VMC/NLEFT and the 3pF input, not the ARM radial density. Please compute and display rho(r) for accepted ARM configurations. If the density is biased, t
- [Section II, after Eq. (4)] The choice M = Δr_{biggest} does not appear valid for rejection sampling. The acceptance criterion compares g(Δr') with M g'(Δr'), where g and g' are probability densities over pair distances; the supremum of g/g' is a dimensionless constant. Setting M to a length (fm) is dimensionally inconsistent and no proof is given that this M satisfies the domination inequality. The authors should report the actual value of sup[g(Δr)/g'(Δr)] used in the sampling, with a concrete prescription for how it was computed.
- [Section II, ARM algorithm description] The algorithm description is ambiguous about how pair distances are enforced for all pairs. The text says 'we accept other nucleons with relative distances Δr′ from the proposal distribution', but it is not clear whether each new nucleon is accepted based on its distance to one previously placed nucleon, to all previously placed nucleons, or to the closest one. If only one pair per nucleon is constrained, the full two-body distribution of accepted configurations need not match g(Δr). A concise pseudocode or an explicit conditional acceptance rule is needed for reproducibility and for the density-preservation check requested above.
- [Section II and Fig. 1] The agreement between ARM and the target C(Δr) in Fig. 1(b,c) is enforced by construction because g(Δr), equivalently C(Δr), is the input to the rejection sampling. Presenting this agreement as evidence that ARM 'captures' the correlations is circular. The non-circular evidence is the agreement in the collision correlators of Figs. 3-6, which are not fitted. The paper should explicitly separate the consistency check (Fig. 1) from the predictive test (Figs. 3-6) and avoid wording that treats the former as independent validation.
minor comments (5)
- [Fig. 1 caption] The caption text says 'results from VMC are depicted in panel (a), while the results from NLEFT are shown in panel (b)', but panel (a) is the one-body density, panel (b) is the VMC two-body distribution, and panel (c) is the NLEFT two-body distribution. The caption should be corrected.
- [Fig. 3 text] The sentence 'we present the results for ε3{2} (middle panels) and ε4{2} (left panels)' should probably read 'right panels' for ε4{2}.
- [Section II] Typos: 'Approch' should be 'Approach'; 'staring from the 3pF density' should be 'starting from'; 'without any constrains' should be 'constraints'.
- [Section III] Typos and grammar: 'a lager discrepancy' -> 'a larger discrepancy'; 'form NLEFT' -> 'for NLEFT'; 'more consist' -> 'more consistent'.
- [Throughout] The paper reports ratios and percentage deviations without statistical uncertainties. For a quantitative claim of 'a few percent' agreement, error bars or at least a statement of sampling statistics would be useful.
Circularity Check
The ARM reproduction of C(Δr) is enforced by construction, but the main TRENTo correlator comparisons are independent and not fitted.
-
self definitional
[Section II, acceptance-rejection procedure following Eq. (4); Fig. 1(b)-(c) and surrounding text]
"we accept other nucleons with relative distances ∆r′ as samples from the the target distribution g(∆r′) if U ≤ (g(∆r′)/M · g′(∆r′)) ... As can be seen in Figs.1(c), the configurations obtained by ARM can almost capture the NN correlations of NLEFT."
The acceptance probability is exactly proportional to g(Δr)/g'(Δr), so rejection sampling guarantees the accepted pair-distance distribution is g(Δr) (the input), up to Monte Carlo noise. Therefore the agreement of ARM with the original C(Δr) in Fig.1 is a self-consistency check of the sampler, not an independent confirmation that NN correlations have been decoded. The actual independent content lies in the TRENTo observables (εn{2}, var(E), cov(εn^2,δE), ε2{4}/ε2{2}), which are not fitted to those correlators.
full rationale
The paper's circular element is limited to the Fig. 1 validation of the two-body distance distribution. Because ARM accepts/rejects candidate nucleons using the very g(Δr) (equivalently C(Δr)) taken from NLEFT/VMC, the reproduction of that distribution is guaranteed by construction; presenting it as a successful reproduction is a sampler check, not evidence for the method's interpretative power. The central quantitative claims, however, are the TRENTo initial-state correlators: eccentricities, energy variance, εn–E covariances, and ε2{4}/ε2{2}. These are computed from the generated configurations and compared with the original ab-initio configurations without using those correlators as acceptance criteria, so they provide independent (and non-circular) tests of whether two-body correlations plus a fixed radial profile encode the relevant structure. The paper does not verify that the rejection step preserves the 3pF one-body density, so the assertion that all differences are angular rests on an unproven premise; this is a correctness/validity concern rather than a circularity. Self-citations (e.g., Refs. [11] and [17]) are present but are not load-bearing for the central derivation, and the loaded Ref. [13] is not by the present authors. Overall, one step is circular by construction, but the main derivation has independent content: score 4.
Assumptions & free parameters
free parameters (3)
- 3pF density parameters =
R0 = 2.608 fm, a0 = 0.513 fm, w = -0.051
- TReNTo p parameter =
0
- Rejection constant M =
Δr_biggest
assumptions (5)
- domain assumption NLEFT and VMC configurations faithfully represent the 16O ground state
- domain assumption 3pF charge density approximates point-nucleon density
- ad hoc to paper Rejection sampling on pair distances preserves the one-body radial density
- domain assumption TReNTo with p=0 captures initial-state entropy deposition
- domain assumption Hydrodynamic response relations (Eqs. 5-8) link initial to final observables
Cite this review
Pith. "Pith review of Impacts of isolated nucleon-nucleon correlations in relativistic $^{16}$O+$^{16}$O collisions." pith.science (2026). https://pith.science/paper/3G6KELU2
@misc{pith2026250900315,
author = {Pith},
title = {Pith review of: Impacts of isolated nucleon-nucleon correlations in relativistic $^16$O+$^16$O collisions},
year = {2026},
howpublished = {\url{https://pith.science/paper/3G6KELU2}},
note = {Machine review of arXiv:2509.00315}
}
abstract
Nucleon-nucleon interactions are fundamental to the nuclear forces operating within the nucleus and play a crucial role in shaping the initial conditions of relativistic ion collisions through two-nucleon correlations. In this paper, we introduce an innovative approach to explore these encoded nucleon-nucleon correlations within advanced \textit{ab-initio} models in the context of relativistic $^{16}O$ collisions. Our methodology successfully reproduces the structural properties of the nucleonic configurations generated by these models, as well as the distance correlations between the nucleon pairs, denoted as $C(\Delta r)$. By generating nucleon positions that align with authentic configurations and adhering to the constraints imposed by the probability distribution of relative two-nucleon distances, our goal is to better understand nucleon-nucleon interactions within \textit{ab-initio} frameworks.
Figures
Figures from the paper (3 more)
Forward citations
Cited by 2 Pith papers
-
Nonlinear collective flow reveals the breakdown of quadrupole--hexadecapole scaling in heavy ion collisions
The nonlinear flow coefficient ξ6,222 in simulated U+U collisions separates the four (β2, β4) nuclear topology classes, making the sign of the hexadecapole deformation β4 experimentally accessible.
-
The size of the quark-gluon plasma in ultracentral collisions: impact of initial density fluctuations on the average transverse momentum
Volume variation with multiplicity in ultracentral collisions is small when total entropy scales with nuclear mass number, as shown by relating it to initial density fluctuation profiles.
Reference graph
Works this paper leans on
-
[1]
Finding the α-cluster parameters by controlling the radial one-nucleon density and incorporating nucleon-nucleon correlations at short-distances (in- side a given cluster) and long-distances (between two clusters) [17]
-
[2]
Generating samples of nucleons based on a fixed charge density function ρ(r) and the constraints derived directly from the density function g(∆r) of nucleus, which accounts for the separation of nu- cleon pairs ∆ r. In this paper, we use the second approach to study the various structures of 16O 3 using the acceptance-rejection method (ARM) [29]. This met...
-
[3]
Probing the core of the strong nuclear interaction
A. Schmidt et al. (CLAS), Nature 578, 540 (2020), arXiv:2004.11221 [nucl-ex]
work page Pith review arXiv 2020
- [4]
-
[5]
R. Cruz-Torres, D. Lonardoni, R. Weiss, N. Barnea, D. W. Higinbotham, E. Piasetzky, A. Schmidt, L. B. We- instein, R. B. Wiringa, and O. Hen, Nature Phys. 17, 306 (2021), arXiv:1907.03658 [nucl-th]
arXiv 2021
-
[6]
Local Two- and Three- Nucleon Interactions Within Chiral Effective Field The- ory,
M. Piarulli, J. Bub, and I. Tews, “Local Two- and Three- Nucleon Interactions Within Chiral Effective Field The- ory,” in Handbook of Nuclear Physics, edited by I. Tani- hata, H. Toki, and T. Kajino (2022) pp. 1–33
work page 2022
-
[7]
Q.-R. Lin, Y.-J. Huang, L.-G. Pang, X.-F. Luo, and X.-N. Wang, (2025), arXiv:2503.01128 [hep-ph]
work page Pith review arXiv 2025
-
[8]
and 5.36 TeV [32] center-of-mass frame. TRENTo is an initial condition model of heavy ion collisions, where the entropy density can be expressed as: dS dy ∝ TR = ( T p A + T p B 2 ) 1 p , 4 ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 5.36 TeV 200 GeV ● NLEFT ● VMC 0 10 20 30 40 50 0.97 0.98 0.99 1.00 1.01 1.02 centralit...
Show all 50 references
-
[9]
Piarulli and I
M. Piarulli and I. Tews, Front. in Phys. 7, 245 (2020), arXiv:2002.00032 [nucl-th]
2020 arXiv
-
[10]
Giacalone et al., (2025), arXiv:2507.01454 [nucl-ex]
G. Giacalone et al., (2025), arXiv:2507.01454 [nucl-ex]
2025 arXiv
- [11]
-
[12]
Citron et al., CERN Yellow Rep
Z. Citron et al., CERN Yellow Rep. Monogr. 7, 1159 (2019), arXiv:1812.06772 [hep-ph]
2019 arXiv
-
[13]
Brewer, A
J. Brewer, A. Mazeliauskas, and W. van der Schee, in Opportunities of OO and pO collisions at the LHC(2021) arXiv:2103.01939 [hep-ph]
2021 arXiv
-
[14]
Mehrabpour and A
H. Mehrabpour and A. Saha, (2025), arXiv:2501.14243 [nucl-th]
2025
-
[15]
Giacalone et al., (2024), arXiv:2402.05995 [nucl-th]
G. Giacalone et al., (2024), arXiv:2402.05995 [nucl-th]
2024 arXiv
-
[16]
Zhang, J
C. Zhang, J. Chen, G. Giacalone, S. Huang, J. Jia, and Y.-G. Ma, (2024), arXiv:2404.08385 [nucl-th]
2024 arXiv
-
[17]
Alvioli, H
M. Alvioli, H. J. Drescher, and M. Strikman, Phys. Lett. B 680, 225 (2009), arXiv:0905.2670 [nucl-th]
2009 arXiv
-
[18]
Huang, Z
Y.-J. Huang, Z. Meng, L.-G. Pang, and X.-N. Wang, (2025), arXiv:2504.00790 [nucl-th]
2025 arXiv
-
[19]
Broniowski and E
W. Broniowski and E. Ruiz Arriola, Phys. Rev. Lett.112, 112501 (2014), arXiv:1312.0289 [nucl-th]
2014 arXiv
-
[20]
Mehrabpour, (2025), arXiv:2506.12673 [nucl-th]
H. Mehrabpour, (2025), arXiv:2506.12673 [nucl-th]
2025
- [21]
- [22]
-
[23]
B.-N. Lu, N. Li, S. Elhatisari, D. Lee, E. Epelbaum, and U.-G. Meißner, Phys. Lett. B 797, 134863 (2019), arXiv:1812.10928 [nucl-th]
2019 arXiv
-
[24]
Elhatisari, E
S. Elhatisari, E. Epelbaum, H. Krebs, T. A. L¨ ahde, D. Lee, N. Li, B.-n. Lu, U.-G. Meißner, and G. Rupak, Phys. Rev. Lett. 119, 222505 (2017), arXiv:1702.05177 [nucl-th]
2017 arXiv
-
[25]
Lonardoni, A
D. Lonardoni, A. Lovato, S. C. Pieper, and R. B. Wiringa, Phys. Rev. C 96, 024326 (2017), arXiv:1705.04337 [nucl-th]
2017 arXiv
-
[26]
Alvioli, C
M. Alvioli, C. Ciofi degli Atti, and H. Morita, Phys. Rev. C 72, 054310 (2005), arXiv:nucl-th/0506054
2005 arXiv
-
[27]
V. R. Pandharipande and R. B. Wiringa, Rev. Mod. Phys. 51, 821 (1979)
1979
-
[28]
Rybczy´ nski and W
M. Rybczy´ nski and W. Broniowski, Phys. Rev. C 100, 064912 (2019), arXiv:1910.09489 [hep-ph]
2019 arXiv
-
[29]
S. C. Pieper, R. B. Wiringa, and V. R. Pandharipande, Phys. Rev. C 46, 1741 (1992)
1992
-
[30]
G. Baym, B. Blattel, L. L. Frankfurt, H. Heiselberg, and M. Strikman, Phys. Rev. C 52, 1604 (1995), arXiv:nucl- th/9502038
1995
-
[31]
Heiselberg, Phys
H. Heiselberg, Phys. Rept. 351, 161 (2001), arXiv:nucl- th/0003046. 10
2001
-
[32]
Bayesian Ideas and Data Analysis: An Introduction for Scientists and Statisticians,
R. Christensen, A. Branscum, and T. E Hanson, “Bayesian Ideas and Data Analysis: An Introduction for Scientists and Statisticians,” in Chapman and Hall/CRC Texts in Statistical Science(2010) 1st ed
2010
-
[33]
Angeli and K
I. Angeli and K. Marinova, Atomic Data and Nuclear Data Tables 99, 69 (2013)
2013
-
[34]
J. S. Moreland, J. E. Bernhard, and S. A. Bass, Phys. Rev. C 92, 011901 (2015), arXiv:1412.4708 [nucl-th]
2015 arXiv
-
[35]
Loizides, (2025), arXiv:2507.05853 [nucl-th]
C. Loizides, (2025), arXiv:2507.05853 [nucl-th]
2025
-
[36]
Schenke, P
B. Schenke, P. Tribedy, and R. Venugopalan, Phys. Rev. C 86, 034908 (2012), arXiv:1206.6805 [hep-ph]
2012 arXiv
-
[37]
Schenke and S
B. Schenke and S. Schlichting, Phys. Rev. C 94, 044907 (2016), arXiv:1605.07158 [hep-ph]
2016 arXiv
-
[38]
Giacalone, F
G. Giacalone, F. G. Gardim, J. Noronha-Hostler, and J.-Y. Ollitrault, Phys. Rev. C 103, 024909 (2021), arXiv:2004.01765 [nucl-th]
2021 arXiv
-
[39]
F. G. Gardim, G. Giacalone, M. Luzum, and J.- Y. Ollitrault, Nucl. Phys. A 1005, 121999 (2021), arXiv:2002.07008 [nucl-th]
2021 arXiv
-
[40]
Samanta, S
R. Samanta, S. Bhatta, J. Jia, M. Luzum, and J.-Y. Ollitrault, Phys. Rev. C 109, L051902 (2024), arXiv:2303.15323 [nucl-th]
2024 arXiv
-
[41]
Bozek and W
P. Bozek and W. Broniowski, Phys. Rev. C 85, 044910 (2012), arXiv:1203.1810 [nucl-th]
2012 arXiv
-
[42]
Bo˙ zek and W
P. Bo˙ zek and W. Broniowski, Phys. Rev. C 96, 014904 (2017), arXiv:1701.09105 [nucl-th]
2017 arXiv
- [43]
- [44]
- [45]
-
[46]
Blaizot, W
J.-P. Blaizot, W. Broniowski, and J.-Y. Ollitrault, Phys. Lett. B 738, 166 (2014), arXiv:1405.3572 [nucl-th]
2014 arXiv
-
[47]
C. Gale, S. Jeon, and B. Schenke, Int. J. Mod. Phys. A 28, 1340011 (2013), arXiv:1301.5893 [nucl-th]
2013 arXiv
- [48]
-
[49]
Bozek and H
P. Bozek and H. Mehrabpour, Phys. Rev. C 101, 064902 (2020), arXiv:2002.08832 [nucl-th]
2020 arXiv
-
[50]
J. P. Delaroche, M. Girod, J. Libert, H. Goutte, S. Hi- laire, S. Peru, N. Pillet, and G. F. Bertsch, Phys. Rev. C 81, 014303 (2010), arXiv:0910.2940 [nucl-th]
2010 arXiv
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.