{"id":"6abdebb2-c252-481e-8bbe-c5c82aecb0e3","arxiv_id":"1908.05292","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A viscous hydrodynamic study shows that normalized symmetric cumulants of flow harmonics can discriminate a crossover QCD equation of state from a first-order one at 62.4 GeV Au+Au collisions.","lead":"This theory paper uses event-by-event hydrodynamic simulations to test whether correlations between flow harmonics can tell apart two candidate equations of state for the quark-gluon plasma. A generalist might read it because it proposes a measurable 'EoS meter' that could help locate the QCD critical point in the RHIC beam energy scan.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"NSC discrimination between EoSs is demonstrated for one first-order bag-model parameter; robustness to B^(1/4) variation is untested.","rationale":"The reader's weakest_assumption identifies exactly this bag-constant uncertainty, and the full text confirms that the authors explicitly acknowledge the allowed range of 150-300 MeV without testing it. This is the most load-bearing concern because it directly tests the generalizability of the central claim: if the NSC separation disappears or reverses for weaker first-order transitions, the claim that NSC can act as an EoS meter is not established. The proposed concrete test is feasible with the existing code and would either validate the robustness or reveal a hidden dependence on the unconstrained parameter. The additional claim about locating the QCD critical point is also unsupported by the simulations, but it is secondary to the EoS-differentiation result. The reader's CONDITIONAL verdict remains appropriate; the condition should be stated as demonstrating robustness of the NSC separation across the allowed range of first-order EoS parameters and, ideally, against an EoS with a critical point.","tokens_in":14123,"tokens_out":3280,"duration_ms":32299,"concrete_test":"Rerun the 20-30% centrality event-by-event simulations with the same MC-Glauber epsilon_WN initial conditions and eta/s=0 and 1/4pi, but replace the first-order EoS with (i) the same bag model with B^(1/4)=150 MeV and (ii) B^(1/4)=300 MeV, or, preferably, with a modern critical-point EoS (e.g., Refs. [48,49]). If the NSC(3,4) difference relative to the crossover EoS does not retain its sign and at least half its magnitude across these variations, the central claim is conditional on the unconstrained bag parameter.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that NSC can differentiate crossover and first-order EoSs rests on simulations using a single first-order EoS, the bag model with B^(1/4)=230 MeV (EoS discussion after Eq. 5). The authors themselves note that 'the choice of bag parameter used here is not unique, it may vary between B^(1/4)=150-300 MeV [8]'. The latent heat and the length of the mixed phase, which govern the hydrodynamic response, depend strongly on B. The observed differences in NSC(2,3), NSC(2,4), and NSC(3,4) (Figs. 5 and 6) could be an artifact of choosing a particularly strong first-order transition. If the true QCD transition at finite baryon density is weaker, or has a smaller latent heat, the magnitude and even the sign of the NSC separation could change, undermining the 'EoS meter' claim. Furthermore, the paper's extrapolation to locating the QCD critical point assumes that a critical-point EoS would produce a similar NSC signature, but such an EoS is not simulated; only a strongly first-order EoS and a crossover EoS are compared.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an event-by-event (2+1)-dimensional viscous hydrodynamic study at finite net baryon density, comparing two equations of state: a lattice-QCD+HRG crossover EoS and a first-order bag-model EoS with B^{1/4}=230 MeV. The authors compute normalized symmetric cumulants NSC(2,3), NSC(2,4), and NSC(3,4) for charged pions at sqrt(s_NN)=62.4 GeV Au+Au collisions, using MC-Glauber and TRENTo initial conditions, and also scan shear viscosity. They report that NSC values differ between the two EoSs, that the difference persists across initial-condition models and viscosity values, and conclude that NSC(3,4) can serve as an 'EoS meter' and possibly help locate the QCD critical point. The hydrodynamic code is tested against the Riemann problem, Gubser flow, and 200 GeV pion spectra and elliptic flow.","tokens_in":14540,"tokens_out":3455,"duration_ms":36359,"significance":"If the central claim holds, the paper would provide a practically measurable observable (NSC) that is sensitive to the QCD equation of state in heavy-ion collisions, complementing earlier work on machine-learning EoS extraction. The study is a forward sensitivity analysis rather than a fit, so there is no circularity: the NSC separation is an output of hydrodynamics with two distinct input EoSs. The code validation against analytic solutions and 200 GeV data is a genuine strength, as is the test against two different initial-condition models and a viscosity scan. However, the significance is conditional because the first-order EoS is represented by a single bag parameter, the critical-point extrapolation is not simulated, and the 62.4 GeV validation is only claimed, not shown.","major_comments":[{"comment":"The central EoS-discrimination claim rests on a single first-order EoS: the bag model with B^{1/4}=230 MeV, chosen to give Tc=164 MeV. The text immediately notes that B^{1/4} may vary between 150 and 300 MeV. The latent heat and the length of the mixed phase, which are the mechanisms invoked to explain the NSC differences, depend strongly on B. A weaker first-order transition or one with a smaller latent heat could reduce or even reverse the NSC separation shown in Figs. 5 and 6. The authors should explicitly test the robustness of NSC(2,3), NSC(2,4), and NSC(3,4) for at least the extreme values B^{1/4}=150 and 300 MeV, or, if that is beyond the exploratory scope, clearly state that the result is conditional on this specific first-order parameterization.","section":"Section II.A"},{"comment":"The extrapolation to locating the QCD critical point is not supported by the simulations shown. Only a crossover EoS and a strongly first-order EoS are compared; no EoS containing a critical point (e.g., those of Refs. [48,49]) is used. The observed NSC difference may be specific to a long first-order mixed phase, whereas a critical point is expected to produce different dynamics, possibly with enhanced fluctuations and critical slowing down. The statement that NSC(3,4) 'can possibly be used to locate the QCD critical point' should be removed or substantially softened unless a critical-point EoS is simulated and shows a similar or distinguishable NSC signal.","section":"Section III / Conclusion"},{"comment":"The validation of the 62.4 GeV setup is only claimed in the text ('We have checked that the above parameters explains the invariant yield of π− across various centralities') but no figure or quantitative comparison is displayed. Since all NSC results are obtained at 62.4 GeV, and the initial parameters epsilon0 and n0 are not otherwise constrained in the paper, the reader cannot assess the model's accuracy at the energy that matters for the main claim. The authors should show the comparison or reference a previous publication where it appears.","section":"Section II.A / Fig. 3"},{"comment":"The statistical and systematic uncertainties of the NSC differences are not quantified in a way that supports the strength of the conclusion. The bootstrap errors are shown, but the number of events is not given, so the statistical significance of the separation between the two EoSs cannot be judged. In addition, resonance decays are neglected (stated for the 200 GeV validation, with the same code used at 62.4 GeV), and the freeze-out energy density differs slightly between the two EoSs (epsilonF=0.28 vs 0.3 GeV/fm3). An estimate of the possible effect of these choices on the NSC separation is needed before claiming that the observable can reliably 'differentiate' the EoSs.","section":"Section II.B / Figs. 5-6"}],"minor_comments":[{"comment":"The terms 'EoS meter' and 'unique observable' overstate the exploratory nature of the study. The paper demonstrates a sensitivity in a specific model setup; softer wording such as 'potential discriminant' would be more accurate.","section":"Abstract / Introduction"},{"comment":"There are typographical errors in the text, for example 'obtsined' instead of 'obtained' and 'decrribed' instead of 'described'. The notation 'm⁄= n' should be typeset as 'm ≠ n'.","section":"Section II.A"},{"comment":"Fig. 6 shows NSC(v2,v3) and NSC(v3,v4) but not NSC(v2,v4), although the text in Section II.B states that NSC(2,3), NSC(2,4), and NSC(3,4) were computed. Please either include NSC(2,4) in the figure or explicitly state why it is omitted.","section":"Fig. 6"},{"comment":"The definitions of v_n and the NSC are clear, but the pT-integration range used for the final v_n and NSC values is not specified. Since the experimental comparison at 200 GeV uses a particular pT window, stating the integration range for the 62.4 GeV NSC results would improve reproducibility.","section":"Section II.B, Eq. (12)"},{"comment":"The notion 'c(ϵ2,v2)' is introduced and used, but the Pearson correlation is computed with what appears to be a small number of events; the number of hydrodynamic events per centrality/initial-condition combination should be stated in the caption or text.","section":"Section II.B, Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful exploratory study with careful numerical tests, but the advertised conclusions go beyond what is demonstrated. The missing B^{1/4} robustness scan and the unvalidated critical-point extrapolation are the main barriers to acceptance. The scope of the journal may also warrant emphasizing the exploratory nature of the result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a legitimate, useful exploratory study. It does what it says: with matched initial conditions, an event-by-event 2+1D viscous hydrodynamic code at finite baryon density, using two EoSs, it shows NSC(2,3), NSC(2,4), NSC(3,4) differ between crossover and a strong first-order transition, across MC-Glauber and TRENTo initial conditions and for several eta/s values. That specific demonstration is new; symmetric cumulants themselves are established, and EoS sensitivity of flow is known, but the systematic separation in a viscous e-by-e code at finite muB was not in the literature. The code is tested against Riemann and Gubser solutions and reproduces the 200 GeV pion yield and v2, which is real evidence the machinery is sound.\n\nThe main soft spots are proportional. First, the first-order EoS is a single bag model with B^(1/4)=230 MeV, and the authors themselves note the bag constant may be 150-300 MeV. Latent heat and mixed-phase length, which drive the hydrodynamic response, change with B. The paper does not scan B, so we do not know if the NSC separation is a feature of first-order transitions in general or of one particular strong one. That does not kill the central observation; it does limit the 'EoS meter' generalization. Second, the leap from 'NSC(3,4) distinguishes crossover vs first-order' to 'possibly locate the QCD critical point' is an extrapolation. A critical-point EoS is not simulated; it may have a weaker signal or a different sign. The paper should either soften that claim or add a critical-point EoS. Third, the 62.4 GeV validation is asserted but not displayed; Fig. 3 only shows 200 GeV. Minor, but easy to fix. Also, no systematic variation of freeze-out energy density or resonance decay treatment, and the code is not released; those are less important for an exploratory paper but worth stating.\n\nThe citation pattern is fine; the relevant ALICE/Bilandzic and EoS references are there. No circularity: the bag constant is tuned to Tc, not to the NSC.\n\nWho this is for: people working on the RHIC BES and on EoS-sensitive flow observables. It is a useful forward sensitivity study, not a final answer. I would send it to peer review; with a B-scan and softened critical-point language it would be a solid paper. Worth a reading group slot, though I would pair it with a critical-point EoS paper.","headline":"A credible exploratory study showing NSC flow correlators separate crossover from first-order EoS in viscous hydro, with the main caveat being the untested bag-constant range and an overreach in the critical-point claim.","tokens_in":14902,"tokens_out":2003,"would_cite":true,"duration_ms":19761,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Flow correlations can tell a first-order phase transition from a crossover in heavy-ion collisions.","keywords":["normalized symmetric cumulants","QCD phase transition","hydrodynamics","equation of state","critical point","flow harmonics","shear viscosity","beam energy scan"],"falsifier":"Compute NSC(3,4) in the same hydrodynamic setup using an equation of state with a critical point of moderate strength (rather than a full first-order transition) and compare it to the crossover case; if the NSC difference falls below the current statistical errors or changes sign, the claim that NSC(3,4) can locate the QCD critical point would be falsified.","tokens_in":13935,"feed_emoji":"💥","tokens_out":1712,"duration_ms":18546,"temperature":0.7,"pith_summary":"This paper argues that normalized symmetric cumulants of anisotropic flow harmonics, observables that can be measured directly in experiments, are sensitive to the equation of state of the quark-gluon plasma. Using event-by-event viscous hydrodynamics, the authors evolve identical collision conditions with two equations of state: a lattice-QCD crossover and a strong first-order phase transition with a bag-model constant. They find that the cumulants, especially NSC(3,4), differ consistently between the two equations of state across different initial conditions and shear-viscosity values. The authors propose NSC(3,4) as a practical 'EoS meter' that could help locate the QCD critical point in beam-energy-scan data.","feed_headline":"Could a flow correlation expose the quark-gluon transition?","feed_subtitle":"A hydrodynamic study says NSC(3,4) can separate a first-order transition from a crossover, offering a route to the QCD critical point.","key_machinery":"The central object is the normalized symmetric cumulant NSC(m,n), defined as SC(m,n)/(<$v_m^{2}$><$v_n^{2}$>), where SC(m,n) is the four-particle cumulant measuring the correlation between the magnitudes of flow harmonics v_m and v_n. The machinery that carries the argument is an event-by-event 2+1-dimensional viscous hydrodynamic simulation with finite baryon density, run with two equations of state: a lattice-QCD+HRG crossover EoS and a first-order bag-model EoS with B^(1/4)=230 MeV. The flow harmonics are generated from fluctuating initial conditions via MC-Glauber and TRENTo models, and the comparison of NSC values between the two EoSs isolates the effect of the phase transition on the hydrodynamic response.","core_discovery":"The paper's central claim is that the normalized symmetric cumulants NSC(m,n) computed from the correlations between flow harmonics v2, v3, and v4 can distinguish an equation of state with a first-order phase transition from one with a crossover, while all other simulation conditions remain the same. This separation is shown to survive changes in the initial condition model (MC-Glauber wounded nucleons versus TRENTo) and variations in shear viscosity up to eta/s = 1/2. The key quantitative results are that c(epsilon2,v2), the Pearson correlation between initial eccentricity and elliptic flow, drops by roughly 5-10 percent for a first-order transition compared with a crossover, and that NSC(3,4) in particular grows with shear viscosity in the first-order case, making it a robust discriminator. The authors explicitly propose NSC(3,4) as a probe of the QCD equation of state that could be computed from existing experimental data across collision energies to search for a sudden change indicative of the QCD critical point.","pith_inferences":["The paper's comparison uses a very strong first-order transition with zero speed of sound over a long mixed phase; if the true QCD transition at finite baryon density is weaker, the magnitude and sign of the NSC difference could change, so the claimed 'EoS meter' may primarily be sensitive to the latent heat and the length of the mixed phase rather than to the mere presence of a critical point.","A natural testable extension would be to run the same hydrodynamic setup with equations of state that include a critical point of varying strength and location, and check whether NSC(3,4) shows a monotonic response; this would tell whether the observable can localize the critical point or only distinguish first-order from crossover transitions.","Since NSC(3,4) grows with shear viscosity in the first-order case but behaves non-trivially for NSC(2,3), the viscosity dependence itself could be used as an additional discriminating handle in experimental data where eta/s is not precisely known.","The 5-10 percent drop in c(epsilon2,v2) for central collisions suggests that the correlation between initial geometry and final flow is a sensitive, albeit model-dependent, indicator of the phase transition; comparing this quantity across initial-condition models in a systematic way could sharpen the prediction."],"forward_implications":["If NSC(3,4) is a robust EoS discriminator, experimental measurements of this cumulant across the RHIC beam energy scan could reveal a sudden change in magnitude indicative of a first-order transition or critical point.","The consistency of the NSC difference across initial conditions and shear viscosity values implies that the observable is relatively robust to model uncertainties, strengthening its use as an experimental probe.","The finding extends to c(epsilon2,v2), suggesting that initial-geometry-flow correlations, though not directly measurable, can serve as a sensitive theoretical diagnostic for the phase transition.","The increase of the NSC(3,4) difference with shear viscosity suggests that viscous effects amplify, rather than erase, the EoS signature, which is a concrete prediction that can be tested against higher-precision hydrodynamics.","The success of this exploratory study motivates constructing equations of state with a critical point and testing whether NSC(3,4) responds to the critical point's location and strength in the same way it responds to a full first-order transition."],"supporting_citations":[{"why":"Defines the symmetric cumulant SC(m,n) and its normalized version NSC(m,n), establishing the observable the paper uses as an EoS probe.","marker":"[45]"},{"why":"Provides the multiparticle cumulant formalism used to compute SC(m,n) from flow harmonics, including the definition of the cumulant that suppresses non-flow correlations.","marker":"[46]"},{"why":"Supplies the parameterized lattice-QCD+HRG crossover equation of state used for one of the two hydrodynamic simulations.","marker":"[18]"},{"why":"Supplies the first-order phase transition equation of state (the bag-model EoS with B^(1/4)=230 MeV) used as the alternative scenario.","marker":"[21]"},{"why":"Establishes the approximate linear correlation between initial eccentricities and final flow coefficients, motivating the Pearson correlation analysis of c(epsilon2,v2).","marker":"[20]"},{"why":"Provides the TRENTo initial condition model used to test the robustness of the NSC difference across initial conditions.","marker":"[65]"},{"why":"Provides the Glauber model and its Monte Carlo implementation used for the wounded-nucleon initial conditions.","marker":"[64]"},{"why":"Supplies the lattice-QCD based equation of state used for the higher-energy collision test at sqrt(s_NN)=200 GeV.","marker":"[47]"},{"why":"Provides the SHASTA algorithm used to solve the hydrodynamic conservation equations, the numerical backbone of the simulation.","marker":"[50]"},{"why":"The deep-learning study that motivated the goal of connecting the QCD equation of state to experimental observables; the paper positions itself as a complementary approach.","marker":"[17]"}],"fun_headline_variants":["Flow cumulants distinguish QCD phase transition types","NSC(3,4) discriminates first-order vs crossover QCD","Cumulant probe for QCD critical point discovery","Hydro study: cumulants act as EoS meter","Normalized cumulants detect QCD phase transition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire comparison relies on the specific strong first-order bag-model equation of state with a bag constant that yields T_c=164 MeV; if the real QCD transition is weaker or has a smaller latent heat, the size or sign of the NSC difference could change.","fun_headline_variants_meta":{"raw":{"variants":["Flow cumulants distinguish QCD phase transition types","NSC(3,4) discriminates first-order vs crossover QCD","Cumulant probe for QCD critical point discovery","Hydro study: cumulants act as EoS meter","Normalized cumulants detect QCD phase transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1333,"prompt_tokens":958,"completion_tokens":375,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":296}},"tokens_in":574,"tokens_out":375,"duration_ms":3935,"temperature":1.0,"reasoning_tokens":296,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:18:03.609971+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute NSC(3,4) in the same hydrodynamic setup using an equation of state with a critical point of moderate strength (rather than a full first-order transition) and compare it to the crossover case; if the NSC difference falls below the current statistical errors or changes sign, the claim that NSC(3,4) can locate the QCD critical point would be falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the SHASTA algorithm used to solve the hydrodynamic conservation equations, the numerical backbone of the simulation."},{"cited_title":"Normalized symmetric cumulants as a measure of QCD phase transition: a viscous hydrodynamic study","cited_arxiv_id":"1908.05292","evidence_quote":"The deep-learning study that motivated the goal of connecting the QCD equation of state to experimental observables; the paper positions itself as a complementary approach."}],"review_version":1}