{"id":"e84cb1d0-feb2-440a-9ca8-8ab44e15f59e","arxiv_id":"1203.2452","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In viscous hydrodynamic simulations of gold and lead collisions, the shear viscosity near the QCD transition temperature (150 to 200 MeV) has the largest influence on elliptic and triangular-to-quartic flow, with the hadronic phase dominating at RHIC and the quark-gluon phase mattering more at LHC.","lead":"A viscous hydrodynamic simulation of heavy-ion collisions shows that the flow of debris is controlled by different temperature windows at different beam energies, and that the region around the quark-gluon transition dominates at all energies. The result sharpens how experiments can convert flow data into a temperature-dependent viscosity of strongly interacting matter.","discovery_kind":"extension","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"Sec. VII's temperature scan, which carries the abstract's headline claim (abstract:s4), never states whether the initial energy density is renormalized to fixed multiplicity, so the reported δvn may mix local viscosity at Ti with global entropy production.","rationale":"The reader's weakest_assumption identifies exactly the right soft spot: the Sec. VII scan is the sole support for abstract:s4, and the text is silent on whether the initial energy density is renormalized per Ti, in contrast to the explicit compensation described in Sec. III C for the four main parameterizations. That asymmetry in bookkeeping is the load-bearing gap, and it is worth flagging because the scan's signal is only a few percent in δv2/v2. I would not, however, move the verdict on it. The paper's own Fig. 19 shows the peak at 150–200 MeV rather than near Tmax, which is the opposite of what a monotonic entropy-production artifact would produce; the qualitative conclusion therefore has internal support independent of the missing statement. The reader's other three points are weaker as load-bearing concerns. The absence of uncertainty quantification is real but is a general practice issue in this literature rather than a defect specific to the argument, and the paper's claims are ordinal (which phase dominates, where the sensitivity peaks) rather than quantitative fits. The per-parameterization normalization of Sec. III C is a tuned knob, but the paper discloses it and shows in Fig. 13 that the grouping of v2(pT) curves survives a Navier-Stokes shear initialization, which is a genuine robustness check. The overlap with Ref. [9] is disclosed in the introduction and in Sec. V, so the novelty is correctly located in the scan and the relaxation-time degeneracy of Fig. 15, not in the RHIC-versus-LHC dominance statement. One correction to the record: the classifier reports 'limitations the authors state: none found', but the text does acknowledge several — the 14-moment δf is called 'merely an ansatz' in Sec. II C, the smooth optical-Glauber initial conditions are flagged in Sec. IV A and Sec. VII, the LHC low-pT v2 undershoot is noted in Sec. IV B, and the Ti = 110 MeV v4 point is explicitly excluded as a freeze-out δf artifact in Sec. VII. Those acknowledgments should be credited rather than re-raised. The remaining unresolved item is the multiplicity matching of the scan, and the test above settles it cheaply. Verdict unchanged: CONDITIONAL, with the Sec. VII bookkeeping as the specific item to confirm.","tokens_in":18672,"tokens_out":4003,"duration_ms":75615,"concrete_test":"Rerun the Sec. VII scan for one energy (RHIC 200 GeV, 20–30% centrality) with the overall normalization Ce of Eq. (12) retuned at each Ti so that the final charged-particle dN/dy at midrapidity matches the constant-η/s = 0.08 baseline to within 1%. Plot δv2/v2 versus Ti exactly as in Fig. 19 and compare peak position and peak height against the published curve. If the peak stays within about 10 MeV and its height within about 20%, the concern does not land and the localization claim stands as stated; if the peak shifts by 20 MeV or more, or flattens appreciably, the claim in abstract:s4 needs to be restated with the multiplicity-matching caveat made explicit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The scan of Sec. VII modifies the baseline η/s = 0.08 by a factor of roughly three over a ±10 MeV window around Ti (Eq. 16 with δT = 10 MeV, Δ = 1.5 MeV), then reads off δvn/vn at each Ti. The problem is that a bump placed at high Ti sits in the early, high-gradient stage where, as the paper itself says in Sec. III C, 'most of the entropy is produced during the early stages of the collision'; a bump placed at low Ti sits in the late, dilute stage. Total entropy production therefore varies systematically and monotonically with Ti across the scan, dragging the fireball size, lifetime, and radial flow with it. Sec. III C explicitly compensates normalization between the four main parameterizations ('We compensate this using different overall normalizations'), but the scan uses a different baseline (constant 0.08 plus a narrow bump) and the text never says whether the same compensation is applied per Ti. If it is not, part of every δvn point is a global size/lifetime effect that grows with Ti, which would bias the extracted 'most sensitive temperature' upward and could in principle manufacture a peak. The effect sizes are small (δv2/v2 of order 2–8% in Fig. 19), so a few-percent systematic is not negligible relative to the signal. There is partial internal counter-evidence: the observed peak sits at 150–200 MeV, far below Tmax = 313/430/504 MeV (Table I), where entropy-production contamination would be largest, so a pure global artifact is unlikely. But the peak's location and width are the quantitative content of abstract:s4, and a 20–30 MeV shift would change the 'narrow band around the transition' reading. This is the one place where the headline claim's method is under-specified.","agreement_with_reader":"agree"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The manuscript studies how a temperature-dependent shear-viscosity-to-entropy-density ratio η/s, different shear relaxation times τπ, and different initial conditions affect the transverse momentum spectra and azimuthal flow asymmetries of charged and identified hadrons in ultrarelativistic heavy-ion collisions. The authors use a (2+1)-dimensional boost-invariant Israel-Stewart hydrodynamic framework with the s95p-PCE-v1 lattice equation of state, Cooper-Frye freeze-out with the 14-moment δf correction, and optical Glauber initial conditions. Four η/s parameterizations (LH-LQ, LH-HQ, HH-LQ, HH-HQ) are compared with PHENIX, STAR, and ALICE data at RHIC 200 GeV, LHC 2.76 TeV, and LHC 5.5 TeV. The main findings are that at RHIC the elliptic flow is dominated by the hadronic viscosity, that only at higher collision energies does the QGP viscosity visibly affect the flow asymmetries, that the shear viscosity near the QCD transition region has the largest impact on collective flow, and that the centrality dependence of v2 is sensitive to the temperature dependence of η/s. A separate temperature scan with a narrow bump in η/s (Sec. VII, Eq. 16) is used to identify the temperature region to which v2 and v4 are most sensitive.","tokens_in":19026,"tokens_out":8966,"duration_ms":166722,"significance":"This is a careful and systematic study that extends the authors' earlier work with a detailed decomposition of the temperature regions probed by v2 and v4 at three collision energies. The finding that flow data constrain η/s mainly in a narrow band around the QCD transition, and that higher harmonics are sensitive to the hadronic viscosity and to the δf correction rather than to the high-temperature viscosity, is a useful and falsifiable prediction for the field. The paper is transparent about its approximations (14-moment δf, smooth initial conditions, optical Glauber, isentropic low-temperature equation of state) and checks several of them explicitly. The main caveat is that the Sec. VII temperature scan does not state whether the initial energy density is renormalized to fixed multiplicity per Ti, which is load-bearing for the quantitative peak location.","major_comments":[{"comment":"The temperature scan does not state whether the initial energy density is renormalized to fixed multiplicity for each Ti. Sec. III C explicitly compensates normalization between the four main parameterizations, but the scan uses a different baseline (constant η/s = 0.08 plus a narrow bump) and the text is silent on this point. Since a bump at high Ti sits in the early high-gradient stage where, as the paper itself says, most entropy is produced, δvn may mix local viscosity at Ti with global entropy production, which would bias the extracted peak location upward. The authors should state whether renormalization was applied, and if not, quantify the systematic by repeating the scan with multiplicity-matched initial conditions. This is load-bearing for the abstract's headline claim [abstract:s4].","section":"Sec. VII, Eq. (16)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The claim that at RHIC the δf corrections clearly dominate v4 rests on the single lowest-temperature point (Ti = 110 MeV) where the bump overlaps the freeze-out surface. At other Ti the difference between Fig. 19 and Fig. 20 is smaller. The claim should be softened or supported by a more systematic decomposition. Moreover, the δf is a 14-moment ansatz, and the conclusion that v4 is dominated by δf is sensitive to the form of δf, which the authors acknowledge is an ansatz. This is load-bearing for the claim that higher harmonics give access to the δf corrections and the hadronic viscosity rather than the high-temperature viscosity.","section":"Sec. VII and Conclusions","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The claim that the centrality dependence of the elliptic flow is sensitive to the temperature dependence of η/s depends on the choice of initial condition. With GLmix, constant η/s = 0.08 describes the most central collisions better, while with BCfit the temperature-dependent hadronic η/s is preferred. The abstract should qualify this dependence on the initial state, or the text should state more explicitly that the conclusion holds within the BCfit initialization.","section":"Abstract [abstract:s5] and Sec. VIII","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The typesetting of the argument of the exponential is garbled; it should read exp((|T − Ti| − δT)/Δ). Please correct.","section":"Eq. (16)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The caption says 'with BC initialization'; this should read 'with BCfit initialization' for consistency with the rest of the paper.","section":"Fig. 12 caption","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The labels 'LHC 2.76 ATeV' and 'LHC 5.5 ATeV' should read 'LHC 2.76 A TeV' and 'LHC 5.5 A TeV' (or 'TeV').","section":"Fig. 16","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The choice λ2 = λ3 = λ4 = 0 is a simplification; the effect of nonzero λ2, λ3, λ4 on the flow is not studied. A sentence acknowledging this would be useful.","section":"Sec. II A","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The scan uses a bump of width ±10 MeV with sharp edges (Δ = 1.5 MeV). The temperature resolution of the scan is limited by the bump width, and the sharp edges may introduce numerical artifacts. The authors should state the temperature resolution and check convergence with respect to Δ.","section":"Sec. VII","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The proton and anti-proton spectra are identical because net-baryon number is not considered. This is acknowledged, but the comparison to PHENIX data in Fig. 5 would benefit from a sentence noting that the data are for p + p̄ and that the model cannot distinguish them.","section":"Sec. IV A","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The hadronic η/s parameterization (Eq. 9) is motivated by Ref. [32], but the authors note that many authors obtain considerably larger values for the shear viscosity of hadronic matter. A brief discussion of how the conclusions would change with a larger hadronic η/s would strengthen the paper.","section":"Sec. III B","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The time evolution of v2 is computed from constant-time hypersurfaces connected to a constant-temperature hypersurface. The paper should state how sensitive the results are to the choice of the connecting surface.","section":"Sec. VI","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The Tmax values are given for the BCfit initialization with the (LH-LQ) parameterization. The text says that for the other initializations Tmax differs by less than 5%. A brief statement of how Tmax varies with the η/s parameterization would be useful.","section":"Table I","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The NS initialization increases entropy production up to 30%, corrected by adjusting the initial energy density. The paper should state how the adjustment is made (e.g., global rescaling of e(τ0, x, y)) and whether it affects the initial eccentricity.","section":"Sec. V","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision"},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"Bottom line: this is a solid, careful viscous-hydro paper, but the headline claim in the abstract is a restatement of the authors' own PRL (Ref. [9]), and they say so. The genuinely new material is the temperature-localization scan of Sec. VII, the τπ-versus-η/s degeneracy in Fig. 15, and the v4 analysis. Those are the parts worth your time; the RHIC-versus-LHC dominance story you already know from their earlier work. The paper is published in PRC and is a standard reference for temperature-dependent η/s, so it is not going anywhere, but read it for the scan, not the abstract. What it does well: the transport coefficients are given in closed form (Eqs. 9, 10, 11), the initial-state parameters are tabulated (Table I), the EoS is a named lattice parameterization, and the freeze-out is standard Cooper-Frye with the 14-moment δf. The authors are candid about their own weak points — the δf ansatz, the smooth optical-Glauber initial conditions, the LHC low-pT v2 undershoot. The τπ degeneracy result (Fig. 15: raising η/s_min from 0.08 to 0.16 can be largely compensated by cτ = 10) is a real caution for anyone doing Bayesian extraction of η/s from v2 alone, and I would cite it for that. Soft spots, in proportion. First, no uncertainty quantification anywhere: every curve is one deterministic run, and the centrality-dependence claims rest on visual comparison. Standard for 2012, but worth flagging. Second, the per-parameterization multiplicity normalization of Sec. III C is a tuned knob that couples to the effect under study. Third, and this is the one that matters for the headline: the Sec. VII scan never states whether the initial energy density is renormalized to fixed multiplicity per Ti, whereas Sec. III C explicitly does so for the four main parameterizations. The stress-test concern lands as an under-specification, not a fatal flaw — the bump is narrow (δT = 10 MeV), and the observed peak at 150–200 MeV sits far below Tmax = 313/430/504 MeV, where entropy-production contamination would be largest, so a pure global artifact is unlikely. But with δv2/v2 of only 2–8%, a few-percent systematic could shift the peak by 20–30 MeV, and that shift is the quantitative content of abstract:s4. A referee should ask for the per-Ti normalization to be stated and the scan shown to be robust to it. Who this is for: hydro practitioners and anyone fitting η/s to flow data. It deserves a serious referee.","headline":"A careful viscous-hydro paper whose real new content — the temperature-localization scan and the τπ degeneracy — is narrower than the abstract implies, and whose Sec. VII scan is under-specified on multiplicity renormalization.","tokens_in":19809,"tokens_out":2059,"duration_ms":39023},"pith_extraction":null,"created_at":"2026-11-13T00:20:32.763390+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"pipeline_version":"pith-pipeline@v1.5.0-archive-selfhost","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":null,"supporting_citations":[],"lane_accuracy":{"status":"not yet measured","lane":{"model":"selfhosted/deepseek-v4.1-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}