{"id":"e9646fda-64d7-4fe4-bba3-e51205bba65e","arxiv_id":"2607.24153","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A 3D entropy deposition model with β≈0.35 and n_BC-dependent rapidity loss, plus ab initio deuteron sampling, reproduces d+Au dNch/dη, spectra, and vn and transfers to p+Au, 3He+Au, and Au+Au.","lead":"A new initial-state entropy recipe with a tunable deposition exponent and collision-number-dependent rapidity loss lets hydrodynamics match charged-particle rapidity shapes in d+Au and related small systems. That matters for upcoming light-ion runs that aim to extract nuclear structure from bulk flow.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The \"universality\" leg of the strongest claim is not parameter-free: β changes 0.35→0.5 for Au+Au, the n_BC-dependent loss is restructured from one-sided to two-sided, and Eq. 18's n_min/n_max normalization is recomputed per system — so the transfer test partially re-fits the very mechanism it claim","rationale":"The reader's weakest_assumption already identifies the circularity burden (β fitted on the data being explained, co-tuned envelope parameters). My scrutiny confirms that is the load-bearing point and sharpens one edge the reader noted but did not fully press: the cross-system transfer — the paper's best defense against circularity — is itself not parameter-free. Appendix A changes both β (0.35→0.5) and the rapidity-loss architecture (one-sided → two-sided n_BC dependence, new baselines), and Eq. 18's min/max normalization is recomputed per system even for p+Au/3He+Au. This does not make the paper wrong; the hydro machinery is standard, the DWF work is careful (if inconsequential for dNch/dη, as the authors themselves show in Figs. 7–8), and the authors explicitly acknowledge the phenomenological, fit-based nature of the framework in Sec. IV. The concern therefore moves neither to REJECT nor to ACCEPT: it is exactly the CONDITIONAL the reader assigned. The proposed cross-fit/hold-out test is cheap (initial-condition generation plus existing CLVisc runs, no new code) and would directly settle whether the transfer constitutes prediction or re-calibration — i.e., whether the universality claim carries evidential weight beyond the fit. If the blind transfers hold, the verdict could reasonably strengthen toward ACCEPT; if they fail, the paper should be reframed as a per-system calibrated module, which is what the CONDITIONAL verdict already anticipates.","tokens_in":28317,"tokens_out":2332,"duration_ms":90436,"concrete_test":"Run a genuine hold-out transfer: freeze ALL of parameter set (d) from d+Au — including the numerical n_min^BC, n_max^BC values in Eq. 18 and the one-sided loss structure — and compute p+Au and 3He+Au dNch/dη blind, changing only the projectile geometry (Woods-Saxon vs single nucleon). Conversely, tune the full parameter set on p+Au+3He+Au (PHENIX) alone and predict the five PHOBOS d+Au centrality curves. If either direction degrades beyond experimental uncertainties (especially the η>2 Au-side shoulder and peripheral bins), the \"universality without re-fit\" claim reduces to per-system calibration and the non-circular support for the d+Au result disappears.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The flagship d+Au agreement is admittedly a fit (β lowered from the motivated 0.5 to 0.35 on the same PHOBOS curves being explained; g_L, g_R, g_P, η_plat, σ_ηgw, σ co-tuned in Table I). The paper's defense against circularity is the transfer to p+Au, 3He+Au, and Au+Au \"without a full re-fit.\" That defense is weaker than stated. (1) For Au+Au (Appendix A), β is changed back to 0.5 AND the rapidity-loss structure is changed: Eq. 18 is applied to both ∆η_L and ∆η_R with new baselines 4.36/4.36, versus d+Au where only the deuteron side carries n_BC dependence with baseline 4.36 and the Au side is fixed at 1.36. Changing both the exponent and the functional architecture of the loss term is a re-fit of the deposition shape, not a transfer. (2) Eq. 18 normalizes ∆η_R by (n_BC − n_min)/(n_max − n_min) with n_min, n_max \"observed across all events\" — i.e., recomputed for each collision system. So even the p+Au/3He+Au application, presented as zero-adjustment, silently re-anchors the projectile-side rapidity loss to each system's own n_BC distribution. For p+Au (one projectile nucleon) this rescaling is doing real work in matching the forward-rapidity slope. (3) The physical reading of β (Sec. II.B, Eqs. 19–21) derives 0.5 from (T_A T_B)^{1/2} energy deposition inherited via free streaming; the fitted small-system value β≈1/3 has no derived counterpart — the summary's interpretation (\"not all energy converted into particles\") is post hoc. The claim that survives scrutiny is narrower: a ~6-parameter ansatz can fit d+Au and, with per-system re-anchoring of its loss normalization and a system-dependent β, also fit p+Au/3He+Au/Au+Au. That is useful phenomenology, but the universality evidence is not the non-circular support the strongest claim asserts. No internal inconsistency — the paper concedes the phenomenological nature in Sec. IV — so this is a calibration-strength issue, not a correctness failure.","agreement_with_reader":"agree"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript addresses a recognized deficiency of (3+1)D hydrodynamic simulations: the failure of factorized (transverse × longitudinal envelope) initial conditions to reproduce charged-particle pseudorapidity distributions in asymmetric d+Au collisions at √s_NN = 200 GeV. The authors (i) sample deuteron configurations from an ab initio Argonne-v18 wavefunction including S–D interference (finding, honestly, that this has negligible effect on dNch/dη), and (ii) introduce a modified 3D entropy deposition ansatz (Eq. 15) with a wounded-nucleon term, a binary-overlap term raised to a power β, and a rapidity-loss shift on the deuteron side that scales linearly with the number of binary collisions n_BC (Eq. 18). With β = 0.35 and the n_BC-dependent loss (set (d), Table I), CLVisc(+SMASH) reproduces PHOBOS dNch/dη in five centrality classes, plus PHENIX identified spectra and v_n. The same framework is then applied to p+Au, ³He+Au (Fig. 13) and Au+Au (Appendix A, with β reset to 0.5 and a two-sided n_BC-dependent loss), and good agreement is reported throughout. The authors claim this demonstrates \"excellent universality\" of the deposition mechanism.","tokens_in":28872,"tokens_out":3538,"duration_ms":197393,"significance":"If the central claims hold at the stated strength, the work is a useful contribution: simultaneous hydrodynamic descriptions of dNch/dη across d+Au centralities have been a persistent failure mode for 3D initial-condition models, and a working, openly specified parametric ansatz with event-by-event CLVisc+SMASH evolution, an explicit centrality-classification cross-check (Figs. 4–6), and coverage of spectra and v_n in four collision systems is of practical value to the community, including for upcoming O+O/Ne+Ne studies. The negative result on the deuteron wavefunction (DWF vs HWF, Figs. 7–8) is also worth publishing. However, the evidence for the headline \"universality\" claim is weaker than stated: the d+Au agreement is an acknowledged multi-parameter fit (β lowered from the motivated 0.5 to 0.35 on the same data being described, with g_L, g_R, g_P, η_plat, σ_ηgw, σ co-tuned), and the Au+Au application changes both β and the functional architecture of the rapidity-loss term. The paper's strength is a well-executed phenomenological fit with partial transfer; it is not a parameter-free or mechanism-validating result, and the abstract and summary should say so.","major_comments":[{"comment":"The claim of 'excellent universality ... without further adjustments' (Sec. III.D) is not supported by the Au+Au application. In Appendix A the authors (i) change β from 0.35 back to 0.5, and (ii) restructure the rapidity-loss term from one-sided (only the deuteron side carries n_BC dependence, with ∆η_L^s fixed at 1.36 for Au) to two-sided with new baselines ∆η_L^s = ∆η_R^s = 4.36. Changing both the exponent and the functional form of the loss term is a re-fit of the deposition shape, not a transfer of a calibrated mechanism. Note also the internal inconsistency this creates: the same Au nucleus is assigned ∆η^s = 1.36 when struck by a deuteron but 4.36 when struck by another Au nucleus. The authors should either (a) demonstrate what the unmodified set (d) predicts for Au+Au and quantify the failure, or (b) remove the universality language from the abstract/Summary and present the Au+Au","section":"Sec. III.D and Appendix A (universality claim)"},{"comment":"The normalization (n_BC − n_min)/(n_max − n_min) with n_min, n_max 'observed across all events' is recomputed separately for each collision system. This means the projectile-side rapidity loss is silently re-anchored to each system's own n_BC distribution even in the p+Au and ³He+Au applications that are presented as zero-adjustment predictions. For p+Au, where the projectile is a single nucleon with a broad n_BC distribution, this rescaling is presumably doing real work in the forward-rapidity slope. The authors should quantify this: show the p+Au/³He+Au dNch/dη obtained with the d+Au-derived ∆η_R^s(n_BC) map applied without re-normalization, so the reader can see how much of the Fig. 13 agreement is genuine transfer versus per-system re-anchoring.","section":"Sec. II.B, Eq. (18)"},{"comment":"The physical motivation for β is internally strained. Eqs. (19)–(21) derive (T_A T_B)^{1/2} scaling for the deposited energy from the energy-flux argument, and the free-streaming paragraph argues the entropy inherits this dependence, motivating β = 0.5. The d+Au fit then requires β = 0.35, for which no derived counterpart exists; the Summary's interpretation ('not all energy from the central fireball is converted into final-state particles') is post hoc and not connected to any mechanism in the text. Since β is the paper's central new ingredient, the authors should either provide a physical argument for a system-dependent β ≈ 1/3 in small systems (e.g., from the transverse-density dependence of energy-to-entropy conversion, which their own free-streaming argument flags as an assumption), or state plainly that β is an empirical exponent and remove the Eqs. (19)–(21) derivation's implied e","section":"Sec. II.B, Eqs. (19)–(21) vs. fitted β = 0.35"},{"comment":"The parameter-tuning procedure is under-documented relative to the weight it carries. Table I shows that sets (a)–(d) vary only β and the loss structure, but the values g_L = g_R = 8.0, g_P = 22.5, η_plat = 1.3, σ_ηgw = 1.3, σ = 2.5, ∆η_L^s = 1.36, ∆η_R^s = 4.36 must themselves have been tuned to the same PHOBOS dNch/dη family, and the stated constraint g_L = g_R 'ensures longitudinal symmetry' is puzzling for an intrinsically asymmetric system. Please (i) describe how these values were obtained and how many effective degrees of freedom the final agreement in Fig. 9(d) represents relative to the five-centrality data; (ii) explain the rationale for g_L = g_R in d+Au; and (iii) justify using the parameter set selected at T_frz = 128 MeV without afterburner (Fig. 9) for the T_frz = 150 MeV CLVisc+SMASH production runs (Fig. 10) without re-checking optimality.","section":"Sec. II.B, Table I and Sec. III.B (parameter tuning)"}],"minor_comments":[{"comment":"The symbol β is used both for the Hulthén wavefunction parameter (Eq. 2, β = 1.18 fm⁻¹) and for the entropy deposition coefficient (Eq. 15). Please rename one of them.","section":"Sec. II.A–II.B (notation)"},{"comment":"In Eq. (18) the symbol ∆η_R^s denotes both the n_BC-dependent function (left-hand side) and the constant baseline (right-hand side). Please distinguish them, e.g., ∆η_R^s(n_BC) = f(n_BC) + ∆η_R^{s,0}.","section":"Eq. (18)"},{"comment":"The fit annotation reads 'R /two.superior= 0.995', presumably a rendering error for R². Also the non-zero intercept (−0.505) is noted but its physical implication for the centrality-classification assumption is not discussed; a sentence would help.","section":"Fig. 5"},{"comment":"The d+Au panel (b) uses PHENIX data with centrality classes 0–5%, ..., 40–60% and the −3.9 < η < −3.1 centrality definition, whereas Figs. 9–10 use PHOBOS classes 0–20%, ..., 80–100% defined via 3.0 < |η| < 5.4. Please state explicitly that the centrality-classification procedure of Sec. II.C was redone for the PHENIX definition, and comment on whether set (d) remains optimal under it.","section":"Fig. 13"},{"comment":"Axis label reads 'Gev' instead of 'GeV'; several spacing artifacts appear in the text ('RESUL TS', 'T RENTo', 'sa mpling', 'CL Visc'). Please proofread.","section":"Fig. 12 caption"},{"comment":"The discussion of the π⁺ underestimate at p_T ≳ 1.6 GeV and of v₂ at high p_T appropriately cites coalescence and subnucleon fluctuations; it would strengthen the paper to state whether these shortcomings are expected to feed back on the fitted β/∆η values if addressed.","section":"Sec. III.C, Figs. 11–12"},{"comment":"TRENTo-2D is used only for Fig. 7, but it is not stated which TRENTo parameters (p, k, σ_w, etc.) were used for that comparison, making the DWF/HWF eccentricity comparison hard to reproduce.","section":"Sec. III.A"}],"recommendation":"major_revision","confidential_remarks":"The technical execution (event-by-event CLVisc+SMASH, centrality cross-check, breadth of observables) is competent and the negative DWF result is honestly reported. The concern is framing: the abstract's \"excellent universality\" overstates what is demonstrated, since the Au+Au leg re-fits both β and the loss architecture, and the Eq. (18) normalization quietly re-anchors the loss per system. This is fixable within the manuscript's scope — either by tempering the claims or by adding the requested no-transfer control calculations — hence major rather than reject. The stress-test note's skepticism on this point does land on the text as written."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real payload here is not the deuteron structure—it is a practical longitudinal entropy formula (β-scaled overlap plus n_BC-dependent rapidity loss) that finally gets d+Au dNch/dη across centralities into decent shape with CLVisc, then carries over to p+Au and 3He+Au with limited retuning.\n\nWhat they do well: the DWF sampling (S–D interference, tensor-force story) is cleanly done and honestly reported as barely moving dS/dη or multiplicity. The hydro pipeline is standard and competent—centrality-from-initial-entropy is cross-checked against final multiplicity ordering, spectra and vn are shown, and they do not hide that bulk/baryon are dropped at 200 GeV. Fig. 9’s parameter scan makes the levers visible: lowering β from 0.5 to 0.35 fixes the periphery; n_BC-dependent Δη helps the forward side in central events. That is useful engineering for people who need asymmetric ICs for light-ion runs.\n\nSoft spots, in proportion. The flagship d+Au rapidity curves are a fit: β, g_L/g_R/g_P, and the envelope widths are co-tuned on the same PHOBOS family they claim to explain. The physical story for β=0.5 from (T_A T_B)^{1/2} is fine motivation; the jump to β≈1/3 is post hoc. The “excellent universality” claim overreaches. For Au+Au they reset β to 0.5 and make the loss two-sided with new baselines; Eq. 18’s n_min/n_max is recomputed per system, so even p+Au/3He+Au silently re-anchor the projectile loss. That is still useful multi-system phenomenology, not a parameter-free transfer. No code or tables of the full IC recipe beyond Table I, so reproducibility is middling.\n\nWho it is for: hydro/IC practitioners and anyone planning O+O, Ne+Ne, Pb+Ne structure-from-flow work who needs a longitudinal handle that does not collapse on asymmetric systems. Math and citations look solid; no internal contradiction—they concede the phenomenological character in the summary.\n\nI would send it to peer review. Ask referees to force a sharper fit-vs-prediction split and to tone the universality language. Worth engaging if you touch small-system ICs; skip if you only care about first-principles deposition.","headline":"Useful calibrated 3D entropy ansatz for asymmetric systems; the deuteron wavefunction is careful but not the lever, and “universality” is weaker than the abstract sells.","tokens_in":29843,"tokens_out":631,"would_cite":true,"duration_ms":23555,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A longitudinal entropy model with deposition coefficient β and collision-number-dependent rapidity loss lets hydrodynamics match charged-particle rapidity distributions in asymmetric d+Au collisions.","keywords":["d+Au collisions","longitudinal entropy deposition","relativistic hydrodynamics","pseudorapidity distributions","binary-collision rapidity loss","small systems","deuteron wave function","anisotropic flow"],"falsifier":"Apply the identical β=0.35 and n_BC-dependent rapidity-loss form, without retuning the longitudinal envelope, to measured dNch/dη in O+O or Ne+Ne at LHC energies; a clear failure across centralities would falsify the claimed universality of the deposition mechanism.","tokens_in":29185,"feed_emoji":"⚛️","tokens_out":1089,"duration_ms":23448,"temperature":0.7,"pith_summary":"Standard relativistic hydrodynamics works well for bulk observables in symmetric heavy-ion collisions but has long failed to reproduce the full charged-particle pseudorapidity distributions in asymmetric systems such as d+Au. This paper argues that the missing piece is mainly in how entropy is laid down along the beam direction in the initial state. The authors sample deuteron nucleon positions from a realistic ab initio wave function and, more decisively, replace the usual factorized longitudinal profile with a three-component entropy density that includes a mid-rapidity interaction term raised to a power β and a rapidity-loss shift that grows with the number of binary collisions on the light-projectile side. With β fixed at 0.35 and that collision-dependent loss, (3+1)D viscous hydrodynamics plus a hadronic afterburner reproduces the measured dNch/dη across five centrality classes at 200 GeV, together with identified-particle spectra and anisotropic flow. The same deposition framework, with only a modest change of β to 0.5 for the larger system, also describes p+Au, 3He+Au, and Au+Au data, suggesting a unified longitudinal initial condition that can be carried over to upcoming light-ion runs.","feed_headline":"Hydrodynamics finally matches d+Au rapidity spectra","feed_subtitle":"A β=0.35 entropy deposition term plus collision-dependent rapidity loss unifies small and large systems at 200 GeV","key_machinery":"The three-component 3D entropy density (Eq. 15): wounded-nucleon Gaussians from each nucleus plus a mid-rapidity plateau term proportional to (sum of left thicknesses × sum of right thicknesses)^β, with beam-directed Gaussians whose centers are shifted by an n_BC-dependent rapidity loss on the light side.","core_discovery":"The paper establishes that charged-particle pseudorapidity distributions in d+Au collisions at 200 GeV are reproduced across centralities once the initial entropy density includes an interaction term scaled by a transverse deposition coefficient β ≈ 0.35 and a rapidity loss on the deuteron side that increases with the number of binary collisions; the same longitudinal deposition form, with β raised to 0.5, also describes p+Au, 3He+Au, and Au+Au without a full re-fit of the longitudinal shape.","pith_inferences":["If β truly tracks system size, a continuous scan from p+A through intermediate systems to A+A should show a smooth rise of the preferred β toward 0.5, offering a diagnostic of when the fireball becomes fully hydrodynamic.","The n_BC-dependent rapidity loss on the light side is effectively a baryon-stopping proxy; the same functional form could be tested against net-proton rapidity distributions once those data are included.","Because the realistic deuteron wave function changes initial eccentricities but barely changes dNch/dη, longitudinal multiplicity is a weak probe of light-nucleus structure, while flow harmonics remain the sharper observable."],"forward_implications":["The same longitudinal entropy prescription can be used as the initial condition for O+O, Ne+Ne, and Pb+Ne collisions at the LHC.","With a reliable longitudinal profile, differences in final-state flow and multiplicity can be attributed more cleanly to the nuclear structure of light projectiles.","Small systems appear to require a smaller entropy deposition coefficient (β≈1/3) than large systems (β=1/2), giving a concrete handle on incomplete energy-to-entropy conversion.","Centrality can be assigned from the initial longitudinal entropy in the forward rapidity window rather than from full hydrodynamic runs, reducing computational cost."],"fun_headline_variants":["β=0.35 entropy term plus n_BC rapidity loss fixes d+Au spectra","Ab initio deuteron sampling and longitudinal deposition unify d+Au data","CLVisc with β-scaled entropy matches d+Au rapidity across centralities","Same deposition form with β=0.5 extends from d+Au to p+Au and Au+Au","n_BC-dependent rapidity loss enables hydro success in asymmetric systems"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The mid-rapidity entropy is assumed to scale as the product of nuclear thicknesses raised to a single adjustable power β that can be lowered from the theoretically motivated value 0.5 down to 0.35 by fitting the same multiplicity data the model is meant to explain.","fun_headline_variants_meta":{"raw":{"variants":["β=0.35 entropy term plus n_BC rapidity loss fixes d+Au spectra","Ab initio deuteron sampling and longitudinal deposition unify d+Au data","CLVisc with β-scaled entropy matches d+Au rapidity across centralities","Same deposition form with β=0.5 extends from d+Au to p+Au and Au+Au","n_BC-dependent rapidity loss enables hydro success in asymmetric systems"]},"model":"grok-4.5","effort":"low","cost_usd":0.003513,"raw_usage":{"total_tokens":1230,"prompt_tokens":858,"num_sources_used":0,"completion_tokens":98,"cost_in_usd_ticks":35128000,"prompt_tokens_details":{"text_tokens":858,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":274,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":858,"tokens_out":98,"duration_ms":6485,"temperature":1.0,"reasoning_tokens":274,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T22:25:55.933951+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Apply the identical β=0.35 and n_BC-dependent rapidity-loss form, without retuning the longitudinal envelope, to measured dNch/dη in O+O or Ne+Ne at LHC energies; a clear failure across centralities would falsify the claimed universality of the deposition mechanism.","supporting_citations":[],"review_version":1}