{"id":"c6022477-42d5-4f97-b60b-d4818189004f","arxiv_id":"2607.14776","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A joint scan over (β2,γ,β4) shows ρ2 tracks quadrupole shape while χ4,22 carries hexadecapole signal at the initial state, but only ρ2 survives full hydrodynamic evolution with the tested statistics.","lead":"This paper maps how changes in nuclear shape alter the flow patterns of particles produced in heavy-ion collisions, varying three deformation parameters at once. It shows which collision observables still carry the shape signal after the hot fireball expands, helping guide future nuclear-structure measurements.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Initial-state scan uses quark-Glauber but final-state iEBE-VISHNU runs use TRENTo, so the claimed loss of beta4 sensitivity in chi4,22 may reflect model mismatch rather than hydrodynamic evolution.","rationale":"The reader's weakest assumption is the proxy assumption that initial-state eccentricities and E/S faithfully represent final flow and mean pT, with the chi4,22 result cited as the main counterexample. I agree that this is an important limitation and that the paper acknowledges it. However, a more immediate and testable confound exists in the paper's setup: the initial-state estimator and the final-state hybrid simulations use different initial condition generators. This is explicit in Sec. II.A (modified quark-Glauber) versus Sec. II.B (TRENTo). The proxy assumption is at least stated and partially self-acknowledged in the Summary, whereas the TRENTo/quark-Glauber mismatch is not discussed. Since the central claim about chi4,22—that a strong initial-state beta4 signal does not survive dynamical evolution—depends on comparing these two different models, this is load-bearing. A single numerical test using TRENTo initial conditions alone (or feeding quark-Glauber profiles into VISHNU) can separate initial-condition-model effects from evolution effects. The initial-state mapping itself is a useful contribution, and the final-state conclusions may well be correct; the concern is about the strength of the evidence. I therefore keep the CONDITIONAL verdict, but the condition should explicitly include a matched-initial-condition check in addition to the statistical and data-release requirements the reader already identified.","tokens_in":14240,"tokens_out":3719,"duration_ms":36729,"concrete_test":"Run the same six Xe configurations through TRENTo alone (no hydro) and compute chi4,22(init) using the same E_n-based definition as Eqs. (17)-(18). Compare its beta4 sensitivity (split between beta4=0.0 and beta4=0.1, and dCor over the 500-point scan or a coarser scan) with the quark-Glauber values in Fig. 3 and with the full iEBE-VISHNU values in Fig. 4. If the TRENTo initial-state split is already negligible, evolution is not the cause; if it is comparable to the quark-Glauber split, then viscous/hadronic effects are genuinely responsible. A complementary check is to feed the same quark-Glauber initial profiles into VISHNU for the six configurations and see whether the beta4 sensitivity survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central comparison between initial-state and final-state sensitivities is not apples-to-apples. All initial-state estimators—Eqs. (16)-(18), the dCor rankings, and the SHAP plots—are computed from the modified quark-Glauber wounded-quark entropy profiles described in Sec. II.A. The final-state iEBE-VISHNU calculations in Sec. II.B instead initialize hydrodynamics with TRENTo, a different initial-condition model with different nucleon/participant weighting and entropy deposition. The key negative result—that the strong beta4 sensitivity of chi4,22 seen at the initial-state level (dCor=0.792, Fig. 3) is 'substantially reduced' after evolution (Fig. 4)—therefore conflates two effects: (1) genuine viscous/hadronic reshaping of the response, and (2) the different response of TRENTo geometry to hexadecapole deformation. Unless the two initial-condition models have identical beta4 sensitivity for these observables, the conclusion that chi4,22's resolving power for beta4 is system-dependent and largely diminished is not established by the presented runs. The same ambiguity affects the qualitative statements about rho224 and rho24. The initial-state mapping itself remains valid, but the paper's practical claim that initial-state sensitivity alone is insufficient for final-state probe power rests on this unmatched comparison.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a multidimensional mapping of how multiparticle heavy-ion observables respond to nuclear deformation, using 129Xe+129Xe collisions at 5.44 TeV as a test case. The authors perform a three-dimensional scan over (β2, γ, β4) with a modified quark-Glauber initial-state model and compute initial-state estimators for ρ2, χ4,22, ρ224, ρ235, and ρ24, characterizing sensitivity with distance correlation and SHAP values. They then run iEBE-VISHNU for six selected deformation configurations and compare final-state responses to the initial-state maps. The main claims are that ρ2 retains its β2/γ sensitivity after full evolution, while the strong initial-state β4 sensitivity of χ4,22 is substantially reduced, and the other higher-order observables cannot be resolved with current statistics.","tokens_in":14523,"tokens_out":4170,"duration_ms":38486,"significance":"If the conclusions are substantiated, this would be a useful contribution to the nuclear-deformation program in heavy-ion collisions: it demonstrates a systematic three-parameter sensitivity analysis rather than isolated one-dimensional scans, and it offers quantitative tools (distance correlation, SHAP) for ranking observable sensitivity. The initial-state scan is based on 500 parameter sets and 4e5 events per set, giving a reasonably dense geometric map. The use of a full iEBE-VISHNU hybrid model for final-state verification is also a strength. However, the central final-state conclusion—that β4 sensitivity of χ4,22 is suppressed after evolution—rests on a comparison between two different initial-condition models, which is a load-bearing issue that must be addressed before the practical recommendations can be accepted.","major_comments":[{"comment":"The initial-state estimators are computed from the modified wounded-quark Glauber model, while the iEBE-VISHNU final-state runs use TRENTo initial conditions. The paper's central negative result—that the β4 sensitivity of χ4,22 (dCor=0.792 at the initial-state level) is 'substantially reduced' after evolution—therefore conflates two effects: genuine viscous/hadronic response and the difference between two initial-condition models. A matched comparison is needed: either recompute the initial-state estimators for the same six configurations using TRENTo, or demonstrate that the quark-Glauber and TRENTo geometries have equivalent β4 sensitivity for these observables. Without this, the statement that initial-state sensitivity is insufficient for final-state probe power is not established.","section":"Sec. II.A–II.B; Fig. 4 and Eqs. (16)–(18)"},{"comment":"The final-state conclusions rely on only six deformation configurations and roughly 10^6 events each, yet the figures show no uncertainty bands or error bars even though the text repeatedly invokes 'statistical uncertainties.' For example, the claim that the β4 difference in χ4,22 is 'indistinguishable within statistical uncertainties' cannot be evaluated without a quantitative uncertainty estimate. Similarly, the null results for ρ224 and ρ24 are presented as 'not resolved with the present statistics,' but the reader has no way to judge the size of those statistics. I recommend adding explicit statistical uncertainties to all final-state figures and to the distance-correlation values, or softening the final-state claims accordingly.","section":"Sec. III.B–III.D and Figs. 2, 4, 6, 9"},{"comment":"The final-state scan covers only two β2 values (0.17, 0.27), three γ values (0°, 30°, 60°), and two β4 values (0, 0.1). This is a reasonable representative set, but the abstract and summary state general conclusions about 'largely preserved' sensitivity and 'substantially reduced' beta-4 sensitivity. The evidence is consistent with those statements, but the sparsity of the grid and the absence of any uncertainty quantification make the strength of the claims disproportionate. I suggest either expanding the final-state configuration set or rewording the conclusions to make the exploratory nature explicit.","section":"Sec. II.B and Table I"}],"minor_comments":[{"comment":"The caption says the upper panel shows γ variation and the lower panel shows β2/β4 variation, but the text and the displayed panel layout appear reversed. Please check and correct the caption/panel order.","section":"Fig. 2 caption"},{"comment":"The numeric entries appear to have lost spaces, e.g., '0.1730°' should probably read '0.17 30°'. Please fix the formatting.","section":"Table I"},{"comment":"There is a typo in 'the then-th eccentricity'—should be 'the n-th eccentricity.'","section":"Sec. II.A"},{"comment":"The dCor values are reported to three decimal places without any statistical uncertainty. Given that they are computed from a finite sample of 500 parameter sets, a bootstrap estimate would help the reader understand the robustness of the ranking.","section":"Sec. II.A and III.A"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the initial-state mapping is a useful contribution. My main concern is that the headline comparison—initial-state sensitivity versus final-state sensitivity—is not apples-to-apples because two different initial-condition models are used. This is fixable in a revision by adding TRENTo-based initial-state estimators for the same configurations, or by clearly limiting the claim to the hybrid model's response and discussing the model-mismatch ambiguity. The absence of uncertainty quantification in the final-state plots is a second issue that should be addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"David — quick take on arXiv:2607.14776. The paper does something genuinely new: instead of the usual one-parameter scans or benchmark configurations, it maps how ρ2, χ4,22, ρ224, ρ235, and ρ24 respond across a joint (β2, γ, β4) space for Xe+Xe, using 500 quark-Glauber configurations at 4e5 events each. The dCor and SHAP rankings are a sensible way to summarize sensitivity, and the headline results — ρ2 dominated by β2/γ, χ4,22 dominated by β4 at the initial-state level — are clearly demonstrated and reproducible in principle. That part is solid and worth having.\n\nThe soft spots are where the paper reaches beyond the initial-state scan. The final-state test uses iEBE-VISHNU with TRENTo, while the initial-state estimators come from a modified quark-Glauber. Those are different entropy deposition models. The claim that χ4,22's β4 sensitivity is 'substantially reduced' after evolution therefore conflates genuine dynamical effects with the different geometry of TRENTo. The paper never acknowledges this mismatch. It may be that the conclusion survives a matched test, but the current runs don't establish it. This is the load-bearing flaw in the final-state part.\n\nAlso, the final-state figures show no uncertainty bands even though the text relies on statistical uncertainties, and only six configurations are used. The higher-order correlators are honestly reported as unresolved, which is fine, but the χ4,22 finding in particular needs either a matched initial-state comparison or a clear statement that the two models are being compared differently.\n\nThe proxy assumption (E_n as a stand-in for V_n) is standard and the paper flags it in the summary, so I don't hold that against them.\n\nBottom line: the initial-state mapping is a real contribution and deserves a serious referee. The final-state comparison needs revision — ideally a matched TRENTo initial-state run or at least a discussion of the model mismatch. I would engage with it.","headline":"A useful 3D sensitivity map of deformation observables, but the initial-to-final comparison is muddied by switching initial-condition models.","tokens_in":15049,"tokens_out":3064,"would_cite":true,"duration_ms":25908,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","25.75.Ld"],"model":"deepseek-v4-flash","headline":"This paper constructs a joint three-dimensional map of how multiparticle flow and momentum correlations respond to quadrupole, triaxial, and hexadecapole nuclear deformation, and finds that only one of the tested probes—ρ2—preserves its ini","keywords":["nuclear deformation","quadrupole deformation","triaxiality","hexadecapole deformation","heavy-ion collisions","flow correlation observables","multiparticle cumulants","initial-state estimators"],"falsifier":"Measure χ4,22 ratios to a spherical reference in 0–5% central 129Xe+129Xe with event statistics high enough to resolve a roughly 10% change: if β4 enhancement of the ratio (≈0.15–0.2 in the initial-state map) is clearly observed, the claim that β4 sensitivity is lost after evolution would be refuted.","tokens_in":14067,"feed_emoji":"⚛️","tokens_out":8337,"duration_ms":66351,"temperature":0.7,"pith_summary":"The paper asks whether multiparticle flow and momentum correlations can jointly constrain three nuclear-shape parameters—quadrupole deformation β2, triaxiality γ, and hexadecapole deformation β4—and maps their sensitivity across the full three-dimensional parameter space using 129Xe+129Xe collisions as a testbed. At the initial-geometry level, the v2–mean-pT correlation ρ2 responds mainly to β2 and γ (distance correlations 0.606 and 0.611 versus only 0.044 for β4), while the nonlinear coefficient χ4,22 responds mainly to β4 (0.792). After full hydrodynamic and hadronic evolution, only ρ2 preserves its β2/γ sensitivity; the β4 signal in χ4,22 is substantially washed out, and the higher-order correlators ρ224, ρ235, and ρ24 cannot be resolved within present statistics. The paper concludes that initial-state sensitivity alone is not proof of final-state probing power, and that joint multidimensional scans rather than one-parameter variations are needed to identify robust deformation observables.","feed_headline":"Only ρ2 keeps its shape signal after full evolution","feed_subtitle":"A 3D scan of quadrupole, triaxial, and hexadecapole deformations shows χ4,22's β4 sensitivity fades in the plasma.","key_machinery":"The central machinery is the initial-state estimator pair: event-by-event eccentricity vectors E_n, computed from the wounded-quark entropy profile, stand in for final flow vectors V_n, and the initial energy-per-entropy ratio E/S stands in for mean transverse momentum [pT]. From these, the paper builds geometric counterparts of ρ2, χ4,22, ρ224, ρ235, and ρ24 (Eqs. 16–19). Sensitivity across the three-dimensional deformation space is quantified with distance correlation (dCor) and SHAP values. Representative deformation points are then pushed through a full viscous-hydrodynamic + hadronic-afterburner hybrid model to test whether each observable's geometric sensitivity survives evolution.","core_discovery":"Using a wounded-quark Monte Carlo Glauber initial condition and an energy-per-entropy proxy for mean transverse momentum, the paper simulates 500 deformation configurations spanning β2∈[0,0.4], γ∈[0°,60°], β4∈[0,0.2] for 0–5% central 129Xe+129Xe collisions. Distance-correlation analysis yields dCor(ρ2,β2)=0.606, dCor(ρ2,γ)=0.611, dCor(ρ2,β4)=0.044, and dCor(χ4,22,β4)=0.792 (vs 0.319, 0.289 for β2,γ). After passing representative configurations through a full viscous-hydrodynamic + hadronic-afterburner hybrid model, the ρ2 sensitivity to β2 and γ survives, the χ4,22 β4 sensitivity is largely lost, and the higher-order correlators are statistically unresolved. This yields a three-dimensional r","pith_inferences":["If the breakdown seen for χ4,22 is generic, then any candidate deformation probe motivated by initial-state geometry should be validated through full evolution before being used to extract nuclear-structure parameters.","The same multidimensional mapping could be applied to octupole deformation β3 or to isobar systems, where correlated deformations may couple in unexpected ways.","A dedicated high-statistics run—either much more than 10^6 events or a larger deformed system—could separate viscous damping from small-system-size effects as the cause of χ4,22's lost β4 sensitivity.","The near-zero dCor of ρ2 with β4 and of ρ235 with β4 suggests those two observables might combine into an approximately orthogonal basis for separating quadrupole/triaxiality from hexadecapole parameters, pending final-state confirmation."],"forward_implications":["ρ2 is a dependable experimental probe of β2 and γ in ultra-central 129Xe+129Xe collisions, because its final-state ratios track the initial-state pattern across the sampled γ and β2 range.","χ4,22 cannot serve as a β4 probe in this system after full evolution; its usefulness is system-dependent, contrasting with results in larger deformed nuclei.","Sensitivity can be conditional: ρ224's γ- and β4-dependence appears mainly at nonzero β2, so single-parameter scans would miss the effect.","dCor and SHAP rankings on geometric estimators offer a scalable screening method for higher-dimensional deformation spaces, including additional parameters beyond (β2,γ,β4).","Resolving the higher-order four-particle correlators requires either higher event statistics in 129Xe+129Xe or studies in larger collision systems."],"fun_headline_variants":["Deformation fingerprints: only ρ2 survives the quark-gluon plasma","Relativistic collisions: β4 signal fades, β2/γ persists in ρ2","Nuclear shapes in heavy-ion data: ρ2 keeps signal, χ4,22 loses it","Three-shape scan shows ρ2 robust, χ4,22 fragile in flow","From nucleus to plasma: only ρ2 retains deformation memory"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that initial-state geometric proxies (eccentricity vectors for flow, energy-per-entropy for mean-pT) faithfully capture how the true final-state flow vector and mean-pT respond to deformation, so that sensitivity rankings derived from those proxies predict which observables will remain deformation-sensitive after viscous evolution and hadronic rescattering.","fun_headline_variants_meta":{"raw":{"variants":["Deformation fingerprints: only ρ2 survives the quark-gluon plasma","Relativistic collisions: β4 signal fades, β2/γ persists in ρ2","Nuclear shapes in heavy-ion data: ρ2 keeps signal, χ4,22 loses it","Three-shape scan shows ρ2 robust, χ4,22 fragile in flow","From nucleus to plasma: only ρ2 retains deformation memory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000582,"raw_usage":{"total_tokens":2662,"prompt_tokens":913,"completion_tokens":1749,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":1645}},"tokens_in":657,"tokens_out":1749,"duration_ms":10421,"temperature":1.0,"reasoning_tokens":1645,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T01:04:21.450978+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure χ4,22 ratios to a spherical reference in 0–5% central 129Xe+129Xe with event statistics high enough to resolve a roughly 10% change: if β4 enhancement of the ratio (≈0.15–0.2 in the initial-state map) is clearly observed, the claim that β4 sensitivity is lost after evolution would be refuted.","supporting_citations":[],"review_version":1}