{"id":"22386897-e0fd-410c-ba43-0f8415c7995a","arxiv_id":"1908.08539","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"First numerical implementation of Hydro+ in an expanding fireball shows critical fluctuations lag behind equilibrium, advect outward, and feed back only mildly on the bulk flow in this simplified model.","lead":"This paper reports the first explicit simulation of Hydro+, a framework that tracks out-of-equilibrium critical fluctuations together with the expanding quark-gluon fluid that carries them. It shows the expected lag of fluctuations, a newly seen outward drift from flow advection, and a small back-reaction on the flow in this simplified model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Back-reaction smallness is only shown for Model-A parameters; the paper's own estimate of the realistic Model-H rate (Gamma0≈0.25 fm^-1) is not tested in Fig. 9.","rationale":"The reader's verdict is CONDITIONAL, and the identified weakest assumption matches ours: the small-back-reaction conclusion is not yet robust because it depends on Model A (z=2), ξmax/ξ0=3, and values of Γ0 that are faster than the effective relaxation appropriate to the Model H dynamics of a real QCD critical point. The paper is transparent about these choices, and the first part of the central claim—that Hydro+ can capture lag, memory, and spatial advection—is convincingly demonstrated within the model. The load-bearing issue is the forward-looking simplification claim, which requires that the few-percent back-reaction persists at slower relaxation and larger ξ enhancement. Our proposed numerical check at Γ0=0.25 fm^-1 and ξmax/ξ0=10 would settle whether that condition holds. We do not see an internal inconsistency or a reason to reject the paper; the conditional verdict stands unchanged.","tokens_in":30427,"tokens_out":6668,"duration_ms":69457,"concrete_test":"Run the identical Hydro+ setup with Γ0=0.25 fm^-1, and separately with ξmax/ξ0=10 (keeping Γ0=1 fm^-1), and compute the maximum of |Δε_B.R./ε| and |Δv_r/v_r| over the (r,τ) grid used in Fig. 9. If either exceeds about 10%, the claim that the back-reaction can be neglected for phenomenological modelling is not robust; if both remain at the few-percent level, the claim is supported in the regime the paper itself flags as more realistic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2's central simplification claim—that the out-of-equilibrium back-reaction on ε and v_r is only a few percent and can be neglected in future modeling—rests on model choices that are not stress-tested in the parameter range the paper itself argues is relevant. The calculation uses Model A dynamics with z=2 (Eq. 2.32) and ξmax/ξ0=3 (Eq. 2.35). At the end of Sec. 3.1 the authors state that for a realistic critical point at nonzero μB, in the Model H universality class (z≈3), the relaxation rate scales as Γ∝(ξ/ξ0)^-3 rather than (ξ/ξ0)^-2, so Hydro+ calculations within this model with Γ0≈0.25 fm^-1 are the appropriate qualitative guide. All back-reaction results in Fig. 9 are presented for Γ0=0.5, 1, 2, and 5 fm^-1; no back-reaction result for Γ0=0.25 fm^-1 (or any Model-H-motivated rate) is shown. Figure 7 itself shows that at Γ0=0.25 the out-of-equilibrium deviations in φ are substantially larger than at Γ0≥0.5, yet the corresponding Δs and pressure feedback are not quantified. No sensitivity to ξmax/ξ0 is given either; since c_V^crit ∝ ξ/ξ0 (Eq. 2.42) and the generalized-pressure corrections grow with ξ, larger enhancement would increase both equilibrium and out-of-equilibrium feedback. The conclusion 'if it persists, this will considerably simplify future modelling' is therefore conditional on an untested corner of parameter space, and the paper's own text identifies that corner as physically relevant for QCD.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first explicit numerical implementation of the Hydro+ formalism in a setting that resembles a heavy-ion collision: a boost-invariant, azimuthally symmetric fireball undergoing radial flow, with a model critical point placed near the μB = 0 axis. The authors fully specify the model ingredients: the equilibrium fluctuation spectrum φ(Q) (Ornstein-Zernike form), the relaxation rate Γ(Q) with Model A dynamics and dynamical exponent z = 2, a temperature-dependent correlation length ξ(T), an equation of state built from a parameterized c_V(T) with a critical contribution, and initial conditions at τ_I = 1 fm. They solve the coupled Hydro+ and viscous hydrodynamic equations and report three main results: (i) the out-of-equilibrium φ(Q) lags behind its equilibrium value, first below it while φbar rises and then above it while φbar falls, including momentum-dependent memory effects; (ii) in the radial profile, a peak in φ(Q) can be advected outward by the radial flow when Γ0 is small; and (iii) the back-reaction of the out-of-equilibrium fluctuations on the energy density and radial velocity is at the level of a few percent for Γ0 = 0.5–2 fm^-1, leading the authors to suggest that future phenomenological modeling may neglect this back-reaction.","tokens_in":31174,"tokens_out":8360,"duration_ms":76755,"significance":"If the results hold, this is a valuable proof-of-principle: it demonstrates that the full Hydro+ system—fluctuation evolution, hydrodynamic evolution, and their mutual feedback—can be solved consistently, and it identifies the momentum range Q ∼ 0.4–0.7 fm^-1 that dominates the feedback integral. The paper is unusually complete in its model specification, and Appendix A gives enough numerical detail for the calculation to be reproduced; the authors also test sensitivity to Γ0 and to box size. The implementation uses the published Hydro+ formalism of Ref. [26] with input parameters chosen a priori, so there is no circularity in the central demonstration. The main significance is the explicit, self-consistent exercise of the formalism and the qualitative guidance it provides. The forward-looking conclusion that back-reaction is negligible is, however, conditional on a parameter region that is not fully explored, and this limits the strength of the paper's main phenomenological suggestion.","major_comments":[{"comment":"The paper's forward-looking conclusion that the back-reaction is small and can be neglected rests on runs with Γ0 = 0.5, 1, 2, and 5 fm^-1 only. Because §3.1 argues that a real QCD critical point at nonzero μB is in the Model H class and should be mimicked in this z = 2 model by a much smaller Γ0 (the Γ0 = 0.25 fm^-1 case shown in Fig. 7 is the strongly out-of-equilibrium case), the parameter point most relevant to the QCD motivation is absent from the back-reaction plots. This is a load-bearing gap: Fig. 7 shows substantially larger deviations of φ(Q) from equilibrium at Γ0 = 0.25 fm^-1, and the feedback integrals in Eqs. (2.11)–(2.15) grow with those deviations. I request either back-reaction results at Γ0 = 0.25 fm^-1 (or at a Model-H-motivated rate) or a clear restriction of the 'considerably simplify future modelling' claim to the tested parameter range.","section":"§3.2, Fig. 9"},{"comment":"The smallness of the back-reaction is also not tested against the assumed correlation-length enhancement, ξmax/ξ0 = 3. Since c_V^crit is proportional to ξ/ξ0 and the generalized-pressure corrections in Eq. (2.15) depend on the critical contribution to the equation of state, a larger ξmax/ξ0 would increase both the equilibrium and the out-of-equilibrium feedback. The manuscript should either provide a scan in ξmax/ξ0 or explicitly state that the percent-level conclusion is contingent on this value.","section":"§2.4, Eq. (2.42)"},{"comment":"The advected peak in φ(Q) at Γ0 = 0.25 fm^-1 is, by the authors' own description, seeded by the initial condition φ = φbar in the outer shell where T is initially near Tc, an initialization they call 'almost certainly unrealistic.' Because this artifact is the basis for the new advection phenomenon highlighted in Fig. 7 and the bright band in Fig. 8, the claim that radial flow transports fluctuations outwards as a generic feature is not yet established. I would like to see a test with a less tuned initialization (for example, φ initialized small everywhere, or with a separate prescription for the shell) or an explicit statement that the advection illustration is conditional on that initial condition.","section":"§2.5 and §3.1, Fig. 7"}],"minor_comments":[{"comment":"The caption says the colored curves come from 'three different values' of Γ0, but the figure and the text present four values: Γ0 = 0.5, 1, 2, and 5 fm^-1.","section":"Fig. 9 caption"},{"comment":"The sentence 'We show φ(Q) at τ = 10.5 fm, when it is at its maximum value, in both panels of Fig. 7' appears to refer to Fig. 5, since Fig. 7 shows τ = 2, 3.5, and 5.5 fm.","section":"§3.1"},{"comment":"There is a grammatical fragment in the sentence beginning 'at a relatively larger, near the edge of the fireball'; this should be rephrased.","section":"§2.5"},{"comment":"The equation has a typesetting artifact, 'sqrt√', which should be cleaned up.","section":"Eq. (2.34)"},{"comment":"The phrase 'the colored \"no B.R.\" energy density curves' is inconsistent with Fig. 9, where the 'no B.R.' case is shown as black dashed curves; the colored curves are the back-reaction cases.","section":"§3.2"},{"comment":"The phrase 'quantitatively capture non-equilibrium phenomena' is stronger than what a single model study with acknowledged simplifications can establish; a phrasing such as 'explicitly compute' or 'self-consistently describe' would better match the caveated statements in the text.","section":"Abstract and §4"}],"recommendation":"major_revision","confidential_remarks":"This is a solid proof-of-principle paper from an authoritative group, and I do not see circularity or any fundamental error in the Hydro+ implementation. My main concern is that the paper's most forward-looking claim—that the back-reaction can be neglected—is not tested in the parameter corner that the authors themselves identify as relevant for a real QCD critical point. If the authors can add back-reaction results at Γ0 = 0.25 fm^-1 (and ideally at a larger ξmax/ξ0), the paper would be substantially strengthened; alternatively, softening the concluding claim would make the manuscript acceptable with minor changes."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is the first explicit Hydro+ simulation in an inhomogeneous, radially expanding fireball, and it works. If you care about critical fluctuations at BES, this is worth a careful read. The central demonstration—lag and memory in phi(Q), feedback into epsilon and v_r, and advection—is credible and well documented. Appendix A gives enough detail to reproduce; the equations are complete; no fitting to a target result, so no circularity worry. The advection claim is appropriately nuanced: the big advected peak in Fig. 7 is seeded by their assumption phi(Q)=phi_eq(Q) initially in the critical shell, which they call unrealistic; they still argue advection will matter generically. That is fair.\n\nThe soft spot is the one the stress-test note hit. The paper's forward-looking simplification claim—that back-reaction is only a few percent and can be neglected—is tested for Gamma_0 = 0.5, 1, 2, 5 fm^-1 with Model A dynamics (z=2) and xi_max/xi_0=3. But the authors themselves note that a realistic QCD critical point at nonzero mu_B is Model H with z~3, with Gamma ∝ (xi/xi_0)^-3, and that within this model the appropriate Gamma_0 is about 0.25 fm^-1. That rate is used in Fig. 7, where it produces much larger out-of-equilibrium deviations in phi, but back-reaction on epsilon and v_r at 0.25 is never shown. And there is no sensitivity scan in xi_max/xi_0 even though all critical effects scale with it. So the 'small back-reaction' conclusion is conditionally true, not demonstrated for the parameter corner the authors identify as QCD-like. They are appropriately hedged in the abstract ('within our model', 'if it persists'), so it is not overclaiming—but the reader should not treat the simplification as established.\n\nMinor point: the numerical regularization at r<0.12 fm is asserted to leave results unaffected but the evidence is not shown; a convergence plot would help. Also, no code is released, though the implementation description is detailed enough that a motivated group could reproduce it.\n\nWho is this for? People building quantitative Hydro+ models for BES observables, and anyone who wants to see how the formalism actually behaves in a finite expanding system. It deserves a serious referee. I would send it out; for publication I would ask for at least a back-reaction estimate at Gamma_0=0.25 or a clear statement that the simplification claim is only for Model A parameters, plus a xi sensitivity point. I would cite it in any Hydro+ context.","headline":"A solid and honest first implementation of Hydro+ with back-reaction in an expanding fireball; the small-back-reaction takeaway is real but demonstrated only in a parameter corner the authors themselves flag as not the QCD-relevant one.","tokens_in":31387,"tokens_out":2368,"would_cite":true,"duration_ms":25301,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","64.60.Ht","25.75.Nq"],"model":"deepseek-v4-flash","headline":"Hydro+ run captures out-of-equilibrium critical fluctuations in QGP","keywords":["QCD critical point","Hydro+","critical fluctuations","out-of-equilibrium dynamics","critical slowing down","heavy-ion collisions","bulk viscosity","beam energy scan"],"falsifier":"Run the same Hydro+ setup with the Model H dynamical exponent $z \\approx 3$ and a correlation-length enhancement $\\xi_{\\max}/\\xi_0$ of ten or more; if the fractional shifts in $\\varepsilon$ and $v_r$ then grow well beyond a few percent, the paper's simplification claim fails.","tokens_in":30216,"feed_emoji":"🔥","tokens_out":8074,"duration_ms":75422,"temperature":0.7,"pith_summary":"This paper reports the first explicit solution of the Hydro+ equations in a simplified but self-consistent model of a heavy-ion collision fireball: a boost-invariant, azimuthally symmetric droplet with a hypothetical critical point placed near zero baryon chemical potential. The authors show that the coupled equations for the hydrodynamic flow and the full spectrum of critical fluctuations can be solved together, and that the resulting dynamics reproduces the expected out-of-equilibrium effects quantitatively: fluctuations grow and decay with a lag induced by critical slowing down, preserving memory of earlier conditions, and are carried outward by radial flow. They further find that the back-reaction of these out-of-equilibrium fluctuations on the energy density and radial flow is small, typically at the one-to-a-few-percent level for relaxation rates $\\Gamma_0$ in the range $0.5$ to $2~\\mathrm{fm}^{-1}$. If that smallness persists in more realistic settings, future phenomenological modeling of the QCD critical-point search could ignore the feedback of fluctuations on the flow and still capture the essential physics.","feed_headline":"Hydro+ run captures out-of-equilibrium critical fluctuations in QGP","feed_subtitle":"First self-consistent simulation shows fluctuations lag, ride radial flow, and feedback at only a few percent.","key_machinery":"The central object is the Wigner transform of the equal-time two-point correlation function of the order-parameter fluctuation, $\\phi_Q(t,x)$, which measures the width of the distribution of the fluctuating field at wave vector $Q$ in the local rest frame. It obeys a relaxation equation $D\\phi_Q = -\\Gamma_Q(\\phi_Q - \\bar{\\phi}_Q)$, with $\\Gamma_Q$ encoding critical slowing down through $\\Gamma_Q \\propto (\\xi/\\xi_0)^{-2}(1+(Q\\xi)^2)$ in Model A dynamics and $\\bar{\\phi}_Q$ the equilibrium Ornstein-Zernike form. The feedback on the flow enters through a generalized entropy $s^{(+)} = s + \\Delta s$, whose dependence on the out-of-equilibrium $\\phi_Q$ modifies the pressure $p^{(+)}$ appearing in the stress-energy tensor; this is what produces the effective bulk viscosity and the modified sound velocity. The model is closed by a temperature-dependent correlation length $\\xi(T)$ peaking at $\\xi_{\\max}/\\xi_0 = 3$, an equation of state built from a $c_V(T)$ with a critical contribution proportional to $\\xi$, and initial conditions with $\\phi_Q$ set to equilibrium at $\\tau_I = 1~\\mathrm{fm}$.","core_discovery":"The central discovery is that Hydro+ can be exercised: the deterministic equations for the equal-time two-point function $\\phi_Q$ of the critical order-parameter fluctuations can be solved alongside second-order viscous hydrodynamics for an expanding, cooling droplet, with each feeding back on the other. In this model the fluctuations rise and fall as the droplet crosses the critical regime, but always lag their equilibrium values because long-wavelength modes relax slowly; at fixed radius the lag is visible as a solid curve that trails the dashed equilibrium curve on both sides of the crossing. The simulation also produces a spatial, memory-dependent pattern of $\\phi_Q$, and shows that radial flow advects a pre-existing peak in the fluctuations outward, a phenomenon not reported before. The quantitative finding is that the resulting modification of the entropy, pressure, energy density and radial velocity is modest: with $\\Gamma_0$ in the range $0.5$--$2~\\mathrm{fm}^{-1}$, the fractional changes in $\\varepsilon$ and $v_r$ are at the few-percent level, adding only about $1.5$--$4.5\\%$ to $\\varepsilon$ relative to the no-critical-point evolution. The authors present this as evidence that the first half of Hydro+, evolving the fluctuations on a fixed hydrodynamic background, may be sufficient for many phenomenological purposes.","pith_inferences":["If the real QCD critical point is in the Model H universality class with $z \\approx 3$, as the authors note, its fluctuations relax more slowly than the $z=2$ law used here, so the lag and back-reaction in a comparable simulation could exceed the few-percent level reported.","The outward advection of a fluctuation peak implies that in event-by-event simulations with lumpy initial conditions, critical fluctuations born in one region could be transported into regions whose thermodynamic history never passed near the critical point, making the correlation between fluctuations and local temperature history a potentially useful diagnostic.","Because the feedback enters through the generalized pressure $p^{(+)}$, the same machinery can be used to isolate the effective bulk viscosity and modified sound speed as functions of $\\xi$; quantifying how those scale with $\\xi_{\\max}/\\xi_0$ would show whether the small back-reaction is a robust property or a consequence of the modest factor-of-three correlation-length enhancement.","A direct comparison of Hydro+ with stochastic simulations of the same model would test how much of the critical dynamics is carried by the deterministic two-point-function truncation and whether higher cumulants, which Hydro+ does not yet evolve, feed back appreciably on the hydrodynamic variables."],"forward_implications":["If the central claim is right, out-of-equilibrium critical fluctuations in a heavy-ion fireball are computable with deterministic Hydro+ equations, not only with stochastic simulations.","The lag and memory effects seen in uniform cooling systems appear in a spatially inhomogeneous, radially expanding fireball, with a new advective transport of fluctuations by the flow.","Back-reaction on $\\varepsilon$ and $v_r$ is small, often below one percent and up to a few percent, so direct fluctuation observables should matter more for the critical-point search than modifications of the bulk flow.","Modeling the QCD critical-point search may be simplified to evolving fluctuations on a hydrodynamic background, omitting the feedback loop, at least within the parameter range studied.","Future 3+1-dimensional Hydro+ with a realistic equation of state and freezeout will be needed before comparison to Beam Energy Scan data, but this paper maps the technical path for doing so."],"supporting_citations":[{"why":"Supplies the Hydro+ formalism itself: the coupled equations for $\\phi_Q$, the generalized entropy $\\Delta s$, and the modified pressure that drive the back-reaction.","marker":"[26]"},{"why":"Establishes the long-standing expectation that critical fluctuations fall out of equilibrium because of critical slowing down, the phenomenon the simulation realizes.","marker":"[23]"},{"why":"Provides the earlier treatment of critical-fluctuation evolution in a cooling uniform plasma whose lag and memory effects the present work reproduces and extends.","marker":"[24]"},{"why":"Supplies the causal second-order viscous hydrodynamic solver and numerical framework that the Hydro+ equations are built upon.","marker":"[41]"},{"why":"Supplies the equilibrium scaling forms for the correlation function, relaxation rate, and specific-heat contribution that parameterize the model.","marker":"[47]"},{"why":"Provides the dynamical universality classification that motivates the $z=2$ Model A relaxation used here.","marker":"[51]"},{"why":"Identifies the QCD critical point with the $z \\approx 3$ dynamical universality class, the contrast against which the paper's $z=2$ choice is calibrated.","marker":"[49]"}],"fun_headline_variants":["Hydro+ run: critical fluctuations lag and advect in QGP fireball","Hydro+ simulation: backreaction from out-of-equilibrium fluctuations is few percent","First Hydro+ run captures memory effects and advection of critical fluctuations","Critical fluctuations lag and ride radial flow in Hydro+ fireball"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative conclusion that the back-reaction is small rests on three modeling choices: a faster relaxation law ($z=2$) than the true QCD critical point likely has, a modest correlation-length growth (only a factor of three), and starting the fluctuations in equilibrium even in the shell initially near the critical temperature; if any of these is more extreme, the lag and back-reaction could be larger.","fun_headline_variants_meta":{"raw":{"variants":["Hydro+ run: critical fluctuations lag and advect in QGP fireball","Hydro+ simulation: backreaction from out-of-equilibrium fluctuations is few percent","First Hydro+ run captures memory effects and advection of critical fluctuations","Critical fluctuations lag and ride radial flow in Hydro+ fireball"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001341,"raw_usage":{"total_tokens":5548,"prompt_tokens":1138,"completion_tokens":4410,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":754,"completion_tokens_details":{"reasoning_tokens":4339}},"tokens_in":754,"tokens_out":4410,"duration_ms":27065,"temperature":1.0,"reasoning_tokens":4339,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:37:30.860694+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same Hydro+ setup with the Model H dynamical exponent $z \\approx 3$ and a correlation-length enhancement $\\xi_{\\max}/\\xi_0$ of ten or more; if the fractional shifts in $\\varepsilon$ and $v_r$ then grow well beyond a few percent, the paper's simplification claim fails.","supporting_citations":[{"cited_title":"Causal Viscous Hydrodynamics for Central Heavy-Ion Collisions","cited_arxiv_id":"nucl-th/0610108","evidence_quote":"Supplies the causal second-order viscous hydrodynamic solver and numerical framework that the Hydro+ equations are built upon."},{"cited_title":"Onuki,Phase Transition Dynamics","cited_arxiv_id":null,"evidence_quote":"Supplies the equilibrium scaling forms for the correlation function, relaxation rate, and specific-heat contribution that parameterize the model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the dynamical universality classification that motivates the $z=2$ Model A relaxation used here."}],"review_version":1}