{"id":"ea02d606-787e-402f-be88-e86aa642f9c2","arxiv_id":"2411.19709","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A new observable, W, built from ratios of flow cumulants and transverse energies in same-nucleus collisions at RHIC and LHC, is proposed and validated as a measure of how close a collision system is to hydrodynamic behavior.","lead":"Physicists propose a new measurement using oxygen-oxygen collisions at two different energies to tell whether particles in small collision systems move together because of a quark-gluon fluid or just a few final bounces. The new observable, W, could settle a long-standing debate about collective flow in small proton and oxygen collisions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Geometry cancellation is the load-bearing premise: W requires eccentricity cumulants to match between RHIC and LHC OO, currently verified only in trento and shown to fail for PbPb/AuAu; test c{4}/c{2}^2 in OO data first.","rationale":"I agree with the reader's identification of the weakest assumption. The geometry cancellation in Eq. (15) is not a calibration detail; it is the step that makes W measurable at all. If the eccentricity cumulants do not cancel between RHIC and LHC OO, Eq. (19) contains an unknown additive geometry term of the same order as the signal. The paper validates the premise for OO only with trento and explicitly documents its breakdown for PbPb/AuAu, which is the strongest internal evidence that the premise is fragile. The nonconformal and calibration-curve issues are real but can in principle be improved with better theory; the geometry cancellation is a prerequisite for any interpretation of W. I therefore keep the reader's conditional verdict: the proposal is promising and internally consistent, but the central claim should not be treated as established until either OO data or an independent initial-state model verifies the same-nucleus geometry cancellation.","tokens_in":21439,"tokens_out":6703,"duration_ms":67602,"concrete_test":"In the upcoming OO datasets, before quoting W, compute the flow-cumulant ratio R = c_ep{4}/c_ep{2}^2 for matched centrality classes at RHIC (200 GeV) and LHC (7 TeV). By Eq. (14), this ratio equals the initial-state cumulant ratio c_epsilon{4}/c_epsilon{2}^2 and is independent of the response; if R_RHIC and R_LHC differ by more than the combined experimental uncertainty, Eq. (15) is violated and geometry contamination in W is non-negligible. A model-level cross-check would be to recompute c_epsilon{2}|RHIC / c_epsilon{2}|LHC with IP-Glasma initial conditions rather than trento; if the deviation exceeds about 5%, the OO geometry premise is model-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that W isolates the dynamical response rests on Eq. (15), which cancels the initial-state geometry by assuming c_epsilon{2}|RHIC / c_epsilon{2}|LHC is close to 1 for the same nucleus at two energies. The paper validates this premise for OO only within trento (Sec. II E, Fig. 2), while the same paper shows the premise fails for PbPb vs AuAu at roughly the 10% level, spoiling W (Sec. II E, Fig. 7). A few-percent geometry mismatch is already consequential: from Eq. (19), W = 4 Delta log(kappa) / Delta log(dE_perp/deta) + 4 Delta log(c_epsilon{2}) / Delta log(dE_perp/deta). Using the Table II opacity ratio of about 0.72 between RHIC and LHC OO, Delta log(dE_perp/deta) is approximately -1.3, so a 5% eccentricity-cumulant mismatch shifts W by about 0.15, a large fraction of the intended 0-to-1 range. The authors themselves note in Sec. II E that other factors, such as nucleon size or nuclear density profiles in IP-Glasma analyses, could produce more significant geometry discrepancies. Until the same-nucleus geometry equality is tested independently, W cannot be claimed as a calibrated measure of hydrodynamization.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an observable W, built from ratios of elliptic-flow cumulants and transverse-energy yields in collisions of the same ion species at two center-of-mass energies. The authors argue that W approximates the logarithmic derivative d log kappa / d log gamma_hat of the flow response curve, with limits W=1 in the free-streaming limit and W=0 in the ideal-hydrodynamic limit, and they calibrate W using Padé fits to conformal RTA kinetic-theory simulations. The construction is validated in event-by-event simulations for OO, PbPb, and AuAu collisions, tested against a nonconformal hydrodynamic setup, and applied to existing PbPb data at 2.76 and 5.02 TeV as a first exploratory extraction.","tokens_in":21890,"tokens_out":6616,"duration_ms":59931,"significance":"If W can be measured reliably for OO collisions at RHIC and LHC, it would provide a genuinely new experimental discriminator between hydrodynamic and few-rescattering explanations of small-system collectivity. The paper is strong in that it derives the observable from a clear physical decomposition of geometry and response, validates the response-factorization assumption on event-by-event kinetic-theory simulations across a wide opacity range, and makes the plot data publicly available. It also honestly identifies the main failure mode of the construction in the PbPb/AuAu comparison. The main significance is therefore conditional: the proposal is attractive and well-motivated, but its utility hinges on the same-nucleus geometry-cancellation premise, which is currently checked only within a single initial-state model.","major_comments":[{"comment":"The central premise that initial-state geometry cancels between the RHIC and LHC ensembles is validated only within the trento model for OO. The paper itself shows in Fig. 7 that the analogous cancellation fails for PbPb vs AuAu at the ~10% level, which spoils W for those systems. Since W is linear in the geometry mismatch through the additive term (2/k) Delta log c_epsilon{2k}/Delta log(E_perp), a few-percent eccentricity-cumulant mismatch can shift W by roughly 0.1, a significant fraction of the intended 0-to-1 range. The authors mention in Sec. II E that IP-Glasma-based analyses could introduce larger geometry differences, but they do not provide an independent, model-insensitive test of the same-nucleus equality. I would ask the authors to make the test explicit: before W is used, the experimental OO cumulant ratio c{4}/c{2}^2 should be compared between RHIC and LHC, since Eq. (14) shows this ratio is geometry-dominated and insensitive to the response mechanism.","section":"Sec. II D–E, Eq. (19), Fig. 7"},{"comment":"The calibration curves d log kappa/d log gamma_hat and f_work(gamma_hat) are Padé fits to the same simulation ensembles from which W is then computed. The agreement in Figs. 5 and 6 is therefore to a significant extent a self-consistency check rather than an independent validation. The comparison with the kappa curve from the previous study in Fig. 1 is only a partial external check, because that curve is scaled by an arbitrary factor 0.93. Consequently, the opacity estimates in Table II do not include systematic uncertainties from the choice of fit ansatz or from the particular centrality samples used for calibration. I request an out-of-sample validation, for example fitting the Padé forms to half of the centrality classes and testing on the other half, or quoting a systematic band that accounts for the fit-ansatz freedom.","section":"Sec. II E, Eqs. (8)–(9), (A6), Figs. 5–6"},{"comment":"The nonconformal test shows that a constant rescaling of the conformal response curve fails at RHIC energies: in Fig. 10 the ratio kappa_LHC/kappa_RHIC acquires an additional centrality dependence already in ideal hydrodynamics, and the nonconformal results deviate from the scaled conformal curve in central and peripheral classes. The authors rescue interpretability by restricting to mid-central collisions or by rescaling with the ideal-hydrodynamic ratio, and the resulting calibration band in Fig. 11 is quite broad. Since the paper's central proposal includes OO collisions at RHIC, this means the current calibration is not yet quantitative for the RHIC end of the program. The abstract and conclusion should either state this limitation explicitly or be accompanied by a more systematic nonconformal calibration, for instance with different equations of state and switching-time scans that go beyond the single vHLLE implementation.","section":"Sec. III, Figs. 9–11"}],"minor_comments":[{"comment":"The meaning of k and the factor 2/k in the definition of W should be spelled out: a reader can easily confuse c{2k} with the 2k-th moment rather than the 2k-th order cumulant, and the k-dependence of the geometry-mismatch term is otherwise obscure.","section":"Eq. (19)"},{"comment":"There is a typo in the caption: 'unadultered' should be 'unadulterated'.","section":"Fig. 10 caption"},{"comment":"The conversion from particle-number-weighted v2 to energy-weighted flow via Ref. [20] is central to the PbPb extraction, but the paper only states the rescaling factors (1.33 and 1.34) without giving the formula. Please quote or reference the relevant equation of that work.","section":"Sec. II F"},{"comment":"The approximate inverse calibration gamma_hat = 2.5 (1-W)^0.78 / W is presented without derivation, validity range, or uncertainty estimate. It should be stated where this approximation deviates from the exact numerical inversion of the calibration curve.","section":"Eq. (24)"},{"comment":"The columns labelled 'actual gamma_hat' and 'mean gamma_hat estimates' are confusing: it should be clarified that the estimates come from plugging the simulation W values at 4 pi eta/s = 1.5 into the inverse calibration, and whether the quoted asymmetric errors include only propagated statistical errors or also calibration-systematic effects.","section":"Table II"}],"recommendation":"major_revision","confidential_remarks":"This is a well-constructed phenomenological proposal with honest reporting of its limitations. My main reservation is that the same-nucleus geometry cancellation is the load-bearing assumption for the flagship OO application, and it is currently verified only within trento, with the PbPb/AuAu case demonstrating a concrete failure mode. The revision should therefore either tighten the calibration and validation strategy or explicitly reframe the paper as proposing a conditional observable whose applicability requires an empirical geometry test. I see this as fixable within the paper's scope, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper proposes an observable W that compares flow cumulants and transverse energies for the same nucleus at two beam energies (O-O at RHIC vs LHC) to place a system on a scale from free streaming (W=1) to ideal hydrodynamics (W=0). The construction is new as far as I know, and the central derivation holds together. It deserves a serious referee.\n\nThe two-energy log-difference construction is a real advance over the earlier cumulant-ratio cancellation of Giacalone et al. The paper validates W in event-by-event conformal RTA kinetic theory across multiple systems and viscosities, and the collapse of the flow response onto a universal curve is honestly presented. They also probe nonconformal effects with a hydrodynamic code and find that at LHC the response is roughly a constant rescaling of the conformal curve, while at RHIC there is extra centrality dependence. They even extract W from PbPb data and get compatibility with the calibration curve. The raw data are public, which makes the tests checkable.\n\nThe main soft spot is the load-bearing premise that the initial geometry is nearly identical between the two energy ensembles. That is plausible for the same nucleus but is validated only with trento. The paper itself shows the PbPb vs AuAu comparison fails at the 10% level, and the stress-test arithmetic is right: a 5% eccentricity-cumulant mismatch shifts W by about 0.15, which is large relative to the intended 0-to-1 range. So the geometry cancellation should be verified before W is sold as a calibrated hydrodynamization measure. A direct check is to compare c{4}/c{2}^2 between RHIC and LHC O-O in data or in an independent initial-state model.\n\nA second caveat is that the calibration curves for kappa(γ-hat) and f_work(γ-hat) are Pade fits to the same simulation ensembles from which W is computed. The agreement with the derivative of those fits is a self-consistency check, not an independent validation. The PbPb data extraction is a useful external anchor, but the opacity scale is tied to one initial state model. Third, the nonconformal test is a single implementation, and the RHIC deviations mean the calibration needs centrality-dependent adjustment; the paper acknowledges this, but it does limit the current quantitative reach.\n\nNet: the idea is promising, the derivation is careful, and the limitations are stated plainly rather than hidden. The paper is for heavy-ion theorists and experimentalists planning the OO run. With a revision that either validates the geometry cancellation independently or narrows the claim accordingly, it would be a solid publication. I would send it to peer review rather than desk reject.","headline":"A genuinely new two-energy ratio observable for hydrodynamization, but the geometry-cancellation premise needs an experimental check before the calibration can be trusted.","tokens_in":22331,"tokens_out":2608,"would_cite":true,"duration_ms":24596,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single flow ratio W, calibrated from 0 (ideal fluid) to 1 (free streaming), measures how hydrodynamic a collision is, and kinetic-theory simulations confirm it.","keywords":["collective flow","hydrodynamization","small collision systems","oxygen-oxygen collisions","opacity","flow response","anisotropic flow cumulants","kinetic theory"],"falsifier":"Measure $W$ in oxygen\\u2013oxygen collisions at RHIC (200 GeV) and LHC (7 TeV) using the proposed cumulant and transverse-energy ratios; if the extracted values do not follow the predicted calibration curve $d\\log\\kappa/d\\log\\hat{\\gamma}$ as a function of the mean opacity, or fall outside the $[0,1]$ interval, the universal-response interpretation fails.","tokens_in":21239,"feed_emoji":"⚛️","tokens_out":11285,"duration_ms":84388,"temperature":0.7,"pith_summary":"Small collision systems such as proton\\u2013nucleus and proton\\u2013proton collisions show collective flow, but it has been unclear whether that flow comes from a hydrodynamic quark\\u2013gluon plasma or from a few final-state rescatterings. This paper proposes an experimental observable, $W$, that settles the question by comparing collisions of the same nucleus, oxygen-16, at two different energies: RHIC and LHC. The comparison cancels the uncertain initial geometry almost completely, leaving only the dynamical response of the system to its own shape. $W$ is calibrated so that 0 means ideal hydrodynamic behavior and 1 means completely noninteracting free streaming, with a monotonic scale in between. If the construction works in data, it will directly show how close small systems are to hydrodynamics and which model families are viable.","feed_headline":"A single flow ratio reveals when quark-gluon plasma turns fluid","feed_subtitle":"Comparing oxygen-16 collisions at two energies cancels geometry, exposing the flow response.","key_machinery":"The load-bearing object is the universal flow response curve $\\kappa(\\hat{\\gamma})$, the ratio of final elliptic flow to initial ellipticity as a function of the opacity $\\hat{\\gamma} = \\frac{1}{5\\eta/s}\\left(\\frac{R}{\\pi}\\frac{dE_0^\\perp}{d\\eta}\\right)^{1/4}$. This curve is linear at small opacity, saturates to a constant at large opacity, and is extracted from kinetic-theory data with a Pad\\'e fit. The proposed observable is $W = \\frac{2}{k}\\frac{\\log\\left(c_{\\varepsilon_p}\\{2k\\}|_A / c_{\\varepsilon_p}\\{2k\\}|_B\\right)}{\\log\\left(\\langle dE_\\perp/d\\eta\\rangle_A / \\langle dE_\\perp/d\\eta\\rangle_B\\right)}$, which approximates $d\\log\\kappa/d\\log\\hat{\\gamma}$ evaluated between the two energy ensembles $A$ and $B$. Two cancellations carry the argument: cumulant ratios such as $c_{\\varepsilon_p}\\{4\\}/c_{\\varepsilon_p}\\{2\\}^2$ eliminate the response and expose geometry, while ratios between RHIC and LHC ensembles of the same nucleus eliminate geometry and expose the response. The transverse-energy ratio in the denominator stands in for the opacity ratio under the assumption that the two ensembles have the same specific shear viscosity and radius.","core_discovery":"Starting from the observation that the elliptic flow response coefficient $\\kappa = \\varepsilon_p/\\epsilon_2$ collapses onto a single universal curve when plotted against a dimensionless opacity $\\hat{\\gamma}$, the paper derives a cumulant ratio whose logarithm is the logarithmic slope of that curve. In event-by-event simulations in conformal relaxation-time kinetic theory, this ratio is shown to follow the predicted calibration curve for oxygen\\u2013oxygen collisions at RHIC and LHC energies, across a range of shear viscosities and centrality classes. The quantity $W$ therefore provides a calibrated, experimentally accessible measure of the degree of hydrodynamization: $W=0$ in the ideal hydrodynamic limit, $W=1$ in the noninteracting limit, and monotonic values in between, with the paper suggesting $W \\lesssim 0.5$ as the hydrodynamic threshold. The same procedure applied to published lead\\u2013lead data at 2.76 and 5.02 TeV gives values compatible with the calibration curve, placing lead\\u2013lead collisions in the hydrodynamic regime.","pith_inferences":["The same two-energy ratio logic should apply to higher harmonics ($v_3$, $v_4$) and to identified-particle flow, which would test whether the universal response curve is truly harmonic-independent and extend the calibration beyond elliptic flow.","A practical consequence the paper leaves implicit is that the limiting experimental systematic for $W$ will be the precision of $dE_\\perp/d\\eta$ at both energies, since the denominator enters through a logarithm; future runs should prioritize transverse-energy measurements in OO collisions.","If the LHC OO run yields $W \\lesssim 0.5$ in mid-central classes, the small-system collectivity debate would shift decisively toward hydrodynamic explanations, and escape-based models would have to reproduce the same $W$ to remain viable.","The failure of the geometry cancellation for PbPb versus AuAu suggests a built-in cross-check: if the cumulant ratio $c\\{4\\}/c\\{2\\}^2$ in OO collisions at RHIC and LHC differs beyond a few percent, the model-independent version of $W$ should be mistrusted."],"forward_implications":["In oxygen\\u2013oxygen collisions at RHIC and LHC, a measured $W$ yields the mean opacity $\\langle\\hat{\\gamma}\\rangle$ of each centrality class, effectively a measurement of the local interaction rate that controls hydrodynamization.","Cumulant ratios such as $c_{\\varepsilon_p}\\{4\\}/c_{\\varepsilon_p}\\{2\\}^2$ become clean probes of initial-state geometry and nuclear structure, independent of how the flowing matter responds.","The calibration curve maps $W$ to a hydrodynamization threshold: the paper\\u2019s kinetic-theory results suggest that systems with $W \\lesssim 0.5$ (opacity $\\hat{\\gamma} \\gtrsim 3$) can be treated as hydrodynamic.","Existing lead\\u2013lead data already produce $W$ values compatible with the calibration curve, indicating that large systems at LHC sit in the hydrodynamic regime; the same test in OO collisions will settle the small-system question.","Nonconformal effects from a realistic equation of state do not destroy the ordering of systems by interaction rate, so $W$ remains interpretable once an adjusted calibration curve is used."],"supporting_citations":[{"why":"Supplies the universal opacity-dependent flow response curve $\\kappa(\\hat{\\gamma})$ and its linear and ideal-hydro limits, which the observable is built to measure.","marker":"[5]"},{"why":"Supplies the work function $f_{\\rm work}(\\hat{\\gamma})$ describing transverse-energy loss, used for the improved calibration curve.","marker":"[6]"},{"why":"Companion paper providing the simulation setup, centrality classes, and rationale behind the two-energy comparison.","marker":"[7]"},{"why":"Provides the quantum Monte Carlo oxygen-16 nucleon configurations used for event-by-event initial states.","marker":"[9]"},{"why":"Establishes that flow response cancels in ratios of cumulants, the step that lets geometry be read directly.","marker":"[14]"},{"why":"Earlier demonstration that $v_2\\{6\\}/v_2\\{4\\}$ and $\\epsilon_2\\{6\\}/\\epsilon_2\\{4\\}$ coincide, motivating the response-cancellation approximation.","marker":"[15]"},{"why":"Supplies lead\\u2013lead $v_2$ measurements at 2.76 and 5.02 TeV used for the first experimental extraction of $W$.","marker":"[16]"},{"why":"Supplies lead\\u2013lead transverse-energy measurements at the two energies, entering the denominator of $W$.","marker":"[17]"},{"why":"Provides the conversion scheme from particle-number-weighted to energy-weighted flow harmonics used in the experimental comparison.","marker":"[20]"},{"why":"Hydrodynamic code used to run the conformal and nonconformal equation-of-state tests of $W$.","marker":"[21]"}],"fun_headline_variants":["Flow ratio distinguishes hydrodynamic flow from rescatterings","Oxygen collisions pin down origin of collective flow","Universal flow response curve marks fluid onset","New flow ratio calibrates quark-gluon plasma fluidity","Threshold W=0.5 signals hydrodynamic behavior"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The construction assumes that the initial-state geometry of oxygen\\u2013oxygen collisions is essentially the same at RHIC and LHC, so the ratio of flow cumulants cancels the eccentricities and leaves only the response ratio; the paper verifies this within its initial-state model for OO, but shows the same cancellation fails for PbPb versus AuAu at roughly the ten-percent level.","fun_headline_variants_meta":{"raw":{"variants":["Flow ratio distinguishes hydrodynamic flow from rescatterings","Oxygen collisions pin down origin of collective flow","Universal flow response curve marks fluid onset","New flow ratio calibrates quark-gluon plasma fluidity","Threshold W=0.5 signals hydrodynamic behavior"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000285,"raw_usage":{"total_tokens":1622,"prompt_tokens":833,"completion_tokens":789,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":449,"completion_tokens_details":{"reasoning_tokens":717}},"tokens_in":449,"tokens_out":789,"duration_ms":7297,"temperature":1.0,"reasoning_tokens":717,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:55:16.209588+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $W$ in oxygen\\u2013oxygen collisions at RHIC (200 GeV) and LHC (7 TeV) using the proposed cumulant and transverse-energy ratios; if the extracted values do not follow the predicted calibration curve $d\\log\\kappa/d\\log\\hat{\\gamma}$ as a function of the mean opacity, or fall outside the $[0,1]$ interval, the universal-response interpretation fails.","supporting_citations":[],"review_version":1}