{"id":"2473a631-bf1f-46c0-a548-13e71d25d85a","arxiv_id":"2509.00315","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Applying rejection sampling to 16O configurations constrained by a Fermi density and a two-nucleon distance distribution reproduces the initial-state eccentricities and energy fluctuations of NLEFT and VMC ab-initio models to a few percent.","lead":"A new sampling method (acceptance-rejection) generates 16O nucleon configurations that match the one-body density and the two-nucleon distance correlation from two ab-initio models, then tests them in simulated relativistic collisions. The method reproduces initial-state observables to within a few percent, but key details about how it preserves the one-body density are missing.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"ARM acceptance depends on pair distance, so the one-body density is not preserved automatically; without checking ρ(r) of ARM configurations, the claimed isolation of NN correlations is unproven.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the acceptance-rejection method is not shown to preserve the one-body radial density, and this unverified premise is central to the paper's interpretation that differences are angular. I agree with that assessment. The paper's own description of the method (Section II) makes it clear that acceptance depends on pair distance, and the radial distribution of accepted samples is therefore a biased version of the proposal 3pF distribution unless a special property holds. The absence of any check of ρ_ARM(r) is a concrete, testable flaw. However, the concern does not necessarily invalidate the entire paper: it can be settled by a simple numerical diagnostic, and if the density is distorted, the method can potentially be corrected. Therefore the appropriate verdict remains CONDITIONAL, as the reader originally assigned. No change in verdict is needed; the condition should explicitly require demonstrating that the ARM configurations reproduce the 3pF one-body density (or appropriately reweighting).","tokens_in":13355,"tokens_out":2931,"duration_ms":37708,"concrete_test":"Compute the one-body radial density of the generated ARM configurations from the same sampled nucleon positions used in the TRENTo simulations (before any recentering), for both NLEFT and VMC settings, and compare to the 3pF input by plotting ρ_ARM(r)/ρ_3pF(r) over r ∈ [0,5] fm. If the ratio deviates from 1 by more than ~1–2% (beyond statistical uncertainty), the isolation claim fails. A complementary check is to repeat the correlator analysis after reweighting accepted configurations to exactly enforce ρ(r)=3pF; if the agreement with original configurations changes materially, the radial distortion is responsible.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section II, the authors generate nucleon radial magnitudes from a 3pF density (Eq. 4) and then accept/reject each new nucleon based on the two-body distance Δr to an existing nucleon, regenerating angular coordinates until acceptance. Because the acceptance probability depends on Δr, and Δr depends on the radial magnitude of the candidate nucleon, the marginal distribution of accepted radial magnitudes is p_3pF(r) times an angular-averaged acceptance factor A(r). Unless A(r) is constant, the resulting one-body density differs from 3pF. The paper never computes or displays ρ(r) for the ARM configurations (Fig. 1(a) shows only VMC, NLEFT, and 3pF, not ARM). This matters because later the paper claims that since all configurations start from the same 3pF radial density, differences are purely angular (Section III, around Ref. [13]). If the ARM radial density is biased, the difference between ARM and original configurations contains radial effects, and the central claim that ARM quantitatively decodes NN correlations collapses. This is not a minor technicality: the method's entire interpretive power rests on preserving the one-body density.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":13674,"tokens_out":3945,"duration_ms":53391,"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":[{"comment":"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":"Section II, Eq. (4) and the acceptance step"},{"comment":"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":"Section II, after Eq. (4)"},{"comment":"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":"Section II, ARM algorithm description"},{"comment":"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.","section":"Section II and Fig. 1"}],"minor_comments":[{"comment":"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.","section":"Fig. 1 caption"},{"comment":"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":"Fig. 3 text"},{"comment":"Typos: 'Approch' should be 'Approach'; 'staring from the 3pF density' should be 'starting from'; 'without any constrains' should be 'constraints'.","section":"Section II"},{"comment":"Typos and grammar: 'a lager discrepancy' -> 'a larger discrepancy'; 'form NLEFT' -> 'for NLEFT'; 'more consist' -> 'more consistent'.","section":"Section III"},{"comment":"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.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the main idea is worth pursuing. The density-preservation issue is substantive but appears testable and fixable with a modest additional figure and text. I would not reject, but the current manuscript is not acceptable until this point is resolved. The M = Δr_biggest issue also suggests the sampling details need a careful re-derivation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe punchline: this paper proposes a rejection-sampling method to import two-body correlations from NLEFT and VMC wave functions into initial conditions for 16O+16O, and it does capture trends in several initial-state correlators. But the central claim—that differences between ARM and the full configurations are purely angular—is not supported, because the authors never check the one-body density of their accepted samples.\n\nWhat is new: previous work studied alpha clustering or fitted densities, but using the two-body distance distribution C(Δr) as the only input to a sampler, then testing against TReNTo observables, is a useful step. The NLEFT vs VMC comparison is well chosen; their C(Δr) differ sharply at short distance, and the ARM reproduces the direction of the resulting changes in ε2{2}, var(E), and the covariances. The paper is candid that it does not yet separate short- and long-range correlations.\n\nThe soft spots are real. The acceptance probability depends on Δr, which depends on the candidate nucleon's radial magnitude, so the accepted configurations will not automatically have the 3pF radial profile. The paper never shows ρ(r) for ARM in Fig. 1. Without that, the claim that differences are angular collapses, and the \"decoding\" of NN correlations is unproven. The algorithm is also underspecified: it is not clear how Δr' is drawn from g', how the candidate position is constructed, and how all pair distances are constrained. This matters for reproducibility. The C(Δr) agreement in Fig. 1 is circular—the target is the input—so the only non-circular evidence is the correlator plots, which lack error bars. The choice of M = Δr_biggest is asserted without justification.\n\nNone of these are fatal to the underlying idea. If the one-body density is indeed distorted, the method could still generate configurations with prescribed pair correlations, but the specific claim about angular isolation is wrong. If the authors add a density check, rewrite the method section with enough detail to reproduce, and include uncertainties, the paper would be a solid contribution.\n\nThis is a paper for people working on light-ion initial conditions and ab-initio to heavy-ion links. It deserves peer review, but the referee should insist on the density check.\n\nBest,\n[Your name]","headline":"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.","tokens_in":14122,"tokens_out":4986,"would_cite":false,"duration_ms":52950,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["nucleon-nucleon correlations","acceptance-rejection method","ab-initio nuclear models","initial-state geometry","relativistic heavy-ion collisions","16O+16O","eccentricity fluctuations","NLEFT and VMC"],"falsifier":"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.","tokens_in":13286,"feed_emoji":"⚛️","tokens_out":7810,"duration_ms":78732,"temperature":0.7,"pith_summary":"The paper claims that the two-nucleon distance correlation C(Δr) of a nucleus, combined with a fixed radial density, carries enough nuclear-structure information to reproduce the initial-state observables of relativistic 16O+16O collisions. Using an acceptance-rejection sampling method, the authors reconstruct nucleon configurations that match the eccentricities, energy fluctuations, and correlation ratios of two ab-initio models, NLEFT and VMC, to within a few percent. If true, this reveals that the relevant structural differences between light nuclei lie in angular two-body correlations, not just the radial density. It also provides a simple way to transport ab-initio structure into initial-condition models of heavy-ion collisions.","feed_headline":"Two-body distances reproduce ab-initio oxygen collisions","feed_subtitle":"Pair-distance correlations plus a fixed radial profile capture oxygen's eccentricity and energy fluctuations","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduces C(Δr) and the constrained-sampling approach for nucleon-nucleon correlations that ARM extends.","marker":"[14]"},{"why":"Supplies NLEFT ab-initio configurations of 16O used as target correlations.","marker":"[20]"},{"why":"Supplies the NLEFT method and pinhole algorithm for determining nucleon positions.","marker":"[21]"},{"why":"Supplies VMC configurations built with the Argonne v18 two-nucleon and Urbana X three-nucleon potentials.","marker":"[22]"},{"why":"Provides the approximate functional form C(Δr)=1-g/g' used in Eq. 3.","marker":"[25]"},{"why":"Defines the acceptance-rejection (rejection sampling) statistical method used throughout.","marker":"[29]"},{"why":"Provides the 3pF density parameters for 16O used as the fixed radial profile.","marker":"[30]"},{"why":"Supplies the TRENTo initial-condition model used to compute the correlators.","marker":"[31]"},{"why":"Establishes that relevant differences between configurations are angular rather than radial; the paper's angular conclusion builds on this.","marker":"[13]"}],"fun_headline_variants":["Pair distances alone reproduce oxygen's collision shapes","Two-nucleon distances encode oxygen's eccentricity","Ab-initio oxygen collisions, reduced to pair distances","One correlation captures oxygen's collision observables","Oxygen collisions without ab-initio: just pair distances"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Pair distances alone reproduce oxygen's collision shapes","Two-nucleon distances encode oxygen's eccentricity","Ab-initio oxygen collisions, reduced to pair distances","One correlation captures oxygen's collision observables","Oxygen collisions without ab-initio: just pair distances"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000873,"raw_usage":{"total_tokens":3592,"prompt_tokens":697,"completion_tokens":2895,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":441,"completion_tokens_details":{"reasoning_tokens":2822}},"tokens_in":441,"tokens_out":2895,"duration_ms":25696,"temperature":1.0,"reasoning_tokens":2822,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T13:44:06.690000+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the acceptance-rejection (rejection sampling) statistical method used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 3pF density parameters for 16O used as the fixed radial profile."},{"cited_title":"Heiselberg, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the TRENTo initial-condition model used to compute the correlators."}],"review_version":1}