Pith. Sign in

REVIEW

Hadron spectra and elliptic flow for 200 A GeV Au+Au collisions from viscous hydrodynamics coupled to a Boltzmann cascade

Nuclear Theory · Submitted 24 Jan 2011 · 45 references

Pith's verdict: A solid companion paper that documents the VISHNU validation honestly in the body but oversells its scope in the abstract.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 10, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 8 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID1101.4638 (version 2)
Pith IDZH6IFV6K, Pith's permanent record of this paper and its review
CitationChun Shen, Huichao Song, Steffen A Bass, Tetsufumi Hirano, Ulrich W Heinz. "Hadron spectra and elliptic flow for 200 A GeV Au+Au collisions from viscous hydrodynamics coupled to a Boltzmann cascade." arXiv:1101.4638 (2011). https://pith.science/paper/ZH6IFV6K
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

A single hybrid simulation, with one shear-viscosity number per initial-state model, reproduces the pion and proton momentum spectra and the elliptic flow at every centrality in 200 A GeV Au+Au collisions at RHIC.

The reading

The paper asks whether one model of the hot, dense matter made in a heavy-ion collision can simultaneously get right three things that are measured separately: how many particles come out, how their momenta are distributed, and how strongly their emission is squeezed into an ellipse. The model, called VISHNU, runs the quark-gluon plasma stage as a viscous fluid with a constant ratio of shear viscosity to entropy density, $T_c < T < 2T_c$, then hands the fluid's output at $T_{sw} = 165$ MeV to a microscopic Boltzmann cascade that follows individual hadrons through their final rescattering. Starting from fluctuating initial entropy profiles generated event by event and then averaged, with two rival prescriptions for how the colliding nuclei deposit entropy, the authors find that the same viscosity that was extracted from the centrality dependence of the integrated elliptic flow also reproduces the differential elliptic flow $v_2(p_T)$ and the $p_T$ spectra of pions and protons across the whole centrality range. The payoff, if the argument holds, is that a single number controls the entire soft sector of the collision, and the paper closes by quoting the constraint $1 < 4π(η/s)_{QGP} < 2.5$.

What carries the argument

The load-bearing object is the hybrid VISHNU pipeline and the compensating pair $(η/s, cτ_0)$ that it makes nearly degenerate in the spectra. The QGP phase is a relativistic viscous fluid with constant $η/s$ between $T_c$ and about $2T_c$; at $T_{sw} = 165$ MeV, chosen to match the measured chemical freeze-out temperature, the fluid is converted to particle distributions and evolved with UrQMD. Initial entropy densities come from two Monte Carlo generators, the MC-Glauber model with finite-size nucleons and a two-component soft-plus-hard entropy deposition, and the MC-KLN model; each event is recentered and rotated so its short axis aligns with the impact parameter, binned by $N_{part}$, and superimposed, so that the smooth profile carries the correct average participant eccentricity $ε_{part}$. The normalization is fixed iteratively to the measured charged multiplicity in the most central bin, and the four parameter pairs used are $(0, 0.4$ fm$)$, $(0.08, 0.6$ fm$)$, $(0.16, 0.9$ fm$)$, and $(0.24, 1.2$ fm$)$ — larger $η/s$ brings more radial flow, which must be offset by a later start time to keep the spectra on top of each other.

Extended reading notes

Core claim

The central claim is that the hybrid code VISHNU, which couples a viscous hydrodynamic description of the quark-gluon plasma to a Boltzmann cascade for the late hadronic stage, gives a uniformly good description of charged and identified hadron spectra and elliptic flow in 200 A GeV Au+Au collisions at RHIC. Concretely, with $(η/s)_{QGP} = 0.08$ for MC-Glauber initial conditions and $(η/s)_{QGP} = 0.16$ for MC-KLN initial conditions — values slightly below but consistent with the $0.10$ and $0.20$ extracted in the companion Letter from the fluctuation- and non-flow-corrected integrated $v_2$ — the calculated $p_T$ spectra of pions and protons and the differential $v_2(p_T)$ of pions, protons, and all charged hadrons agree with STAR and PHENIX data over essentially all centralities. The paper stresses that the spectra are nearly insensitive to $η/s$ once the initial entropy normalization and the starting time $τ_0$ are readjusted together, so the discrimination comes from $v_2/ε$, which degrades visibly when $η/s$ is raised by $1/4π$. The authors also state plainly where the description fails: the two most peripheral centrality bins, where the calculated $v_2$ falls below the data and they say they do not know how to recover the measured values, and the two most central bins, where STAR and PHENIX disagree with each other by about 30% and the PHENIX data would demand a larger effective $η/s$.

Load-bearing premise

The whole extraction rests on the initial fireball eccentricity computed from MC-Glauber or MC-KLN, two models that differ by up to 25% depending on centrality, and the paper itself says this uncertainty completely dominates the error on the extracted viscosity; the appendix shows the preferred value moving from 0.08 toward 0.16 when one unmeasured input of the Glauber model, the transverse smearing area, is changed.

Editorial extensions

If this is right

  • The soft hadron sector at RHIC — multiplicity, $p_T$ spectra of pions and protons, and both integrated and differential elliptic flow — can be organized by a single viscosity parameter once the initial-state model is fixed.
  • The extracted range $1 < 4π(η/s)_{QGP} < 2.5$ becomes a quantitative target for calculations of transport coefficients in hot QCD.
  • Because the spectra are nearly $η/s$-independent after retuning, the paper argues that viscosity should be extracted from the $p_T$-integrated momentum anisotropy rather than from the shape of $v_2(p_T)$.
  • The failure at peripheral centralities and the near-central STAR/PHENIX disagreement are localized, identifiable discrepancies rather than diffuse disagreement with the whole data set.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper's own appendix shows that shrinking the MC-Glauber transverse smearing area from 42 mb to 4.2 mb pushes the preferred $η/s$ from 0.08 toward 0.16 and spoils the slope of $v_2/ε$ versus $(1/S)dN_{ch}/dy$; an independent, non-model determination of the initial eccentricity would therefore be worth more to this program than any improvement in the flow data themselves.
  • If the $v_2/ε$ versus $(1/S)dN_{ch}/dy$ curve is truly universal at fixed collision energy and depends only on $η/s$, then repeating the same comparison at LHC energies tests whether that universality survives a factor-of-two jump in the multiplicity density, which is a sharper statement than simply checking whether $η/s$ runs.
  • The peripheral-collision failure, which the authors note also appears in an ideal-hydrodynamics cascade using a different afterburner, suggests the missing physics sits in the late hadronic stage or in the treatment of the dilute edge of the fireball rather than in the QGP viscosity.
  • The 30% STAR/PHENIX split in the most central bin, which the authors attribute to a possible 2% centrality-definition shift, is a concrete bookkeeping question: resolving it would decide whether the near-central data really call for a larger effective $η/s$ than the mid-central data.
Original abstract

It is shown that the recently developed hybrid code VISHNU, which couples a relativistic viscous fluid dynamical description of the quark-gluon plasma (QGP) with a microscopic Boltzmann cascade for the late hadronic rescattering stage, yields an excellent description of charged and identified hadron spectra and elliptic flow measured in 200 A GeV Au+Au collisions at the Relativistic Heavy-Ion Collider (RHIC). Using initial conditions that incorporate event-by-event fluctuations in the initial shape and orientation of the collision fireball and values eta/s for the specific shear viscosity of the quark-gluon plasma that were recently extracted from the measured centrality dependence of the eccentricity-scaled, p_T-integrated charged hadron elliptic flow, we obtain universally good agreement between theory and experiment for the p_T-spectra and differential elliptic flow v_2(p_T) for both pions and protons at all collision centralities.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 2,550 input tokens, $0.0001 reported.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.06 · human participants 0.02 · controlled comparison 0.06 · randomized assignment 0.01 · states quantitative result 0.35 · compares with prior results 0.20 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

It is shown that the recently developed hybrid code VISHNU, which couples a relativistic viscous fluid dynamical description of the quark-gluon plasma (QGP) with a microscopic Boltzmann cascade for the late hadronic rescattering stage, yields an excellent description of charged and identified hadron spectra and elliptic flow measured in 200 A GeV Au+Au collisions at the Relativistic Heavy-Ion Collider (RHIC). Using initial conditions that incorporate event-by-event fluctuations in the initial shape and orientation of the collision fireball and values eta/s for the specific shear viscosity of the quark-gluon plasma that were recently extracted from the measured centrality dependence of the eccentricity-scaled, p_T-integrated charged hadron elliptic flow, we obtain universally good agreement between theory and experiment for the p_T-spectra and differential elliptic flow v_2(p_T) for both pions and protons at all collision centralities.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 It is shown that the recently developed hybrid code VISHNU, which couples a relativistic viscous fluid dynamical description of the quark-gluon plasma (QGP) with a microscopic Boltzmann cascade for the late hadronic rescattering stage, yields an excellent description of charged and identified hadron spectra and elliptic flow measured in 200 A GeV Au+Au collisions at the Relativistic Heavy-Ion Collider (RHIC). result 0.84 0.89 0.87 asserted 1.00
s1 Using initial conditions that incorporate event-by-event fluctuations in the initial shape and orientation of the collision fireball and values eta/s for the specific shear viscosity of the quark-gluon plasma that were recently extracted from the measured centrality dependence of the eccentricity-scaled, p_T-integrated charged hadron elliptic flow, we obtain universally good agreement between theory and experiment for the p_T-spectra and differential elliptic flow v_2(p_T) for both pions and protons at all collision centralities. result 0.88 0.93 0.29 asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.99 · evidence from simulation or computation 0.99.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1: (Color online)
Figure 2. FIG. 2: (Color online) Eccentricity-scaled elliptic flow
Figure 3. FIG. 3: (Color online) Same as Fig. 2, but with
Figure 4. FIG. 4: Comparison of the STAR [19] and PHENIX
Figure 5. FIG. 5: (Color online) Same as Fig. 2, but for identified pions
Figure 6. FIG. 6: (Color online) Integrated charged hadron elliptic fl
Figure 7. shows the eccentricity (panel (a)) and trans￾verse area (panel (b)) of the participant-plane averaged initial entropy density distribution as a function of the number Npart of participant (“wounded”) nucleons in Au+Au collisions, for two values of the smearing area σs. We see that a smaller smearing area (more pointlike particle production) leads to larger initial fireball eccen￾tricities εpart and smaller transverse areas S. While the effect of varying σs on S is simply an offset, for εpart 0 0.1 0.2 0.3 0.4 0.5 0.6 εpart -- σs = 42 mb σs = 4.2 mb 0 100 200 300 400 Npart 0 5 10 15…
Figure 8. FIG. 8: (Color online) (a) Centrality dependence of
Figure 9. FIG. 9: (Color online) Comparison of the universal theoret

Figures reproduced from arXiv:1101.4638 by the authors.

References
How the literature uses it

48 arXiv papers cite this in 61 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 40 sentences without a confident label

Uses its method or tool 4 papers
  • In this work, we opt for a mature, well-tested set of event-by-event models [ref] with an established track record of describing diverse RHIC and LHC data [ref] .

    Quantifying properties of hot and dense QCD matter through systematic model-to-data comparison · Model · 0.87

  • In order to get the quenched leading hadron spectra at high $p_T$, this parton-level distribution is convoluted with LO pQCD pp spectra (for gluon and ligh quark) [ref] or FONLL pp spectra (for charm and bottom) [ref] , folded over the MC-Glauber A+A initial hard scattering probability distribution [ref] , and hadronized with parton fragmentation functions (LO KKP [ref] for gluon and light quark[There are systematic studies that confront current NLO parton fragmentation functions with inclusive charged-particle spectra at hadron colliders and compare the varied consistencies between NLO FFs and data, cf. e.g. [ref] .

    Bridging Soft-Hard Transport Properties of Quark-Gluon Plasmas with CUJET3.0 · Perturbative QCD jet quenching: DGLV opacity expansion · 0.87

  • We neglect rescatterings in the hadronic phase [ref] , but implement strong decays of all hadronic resonances into stable hadrons using an adapted version of the decay code of [ref] .

    Hydrodynamic predictions for 5.44 TeV Xe+Xe collisions · Hydrodynamic evolution · 0.99

Uses its data or values 3 papers
Compares with it 4 papers
  • In Fig. [ref], the $v_2$ measurements for two centrality intervals, the 10–20$%$ in the left column and the 40–50$%$ interval in the right column, are compared to hydrodynamic calculations coupled to a hadronic cascade model (VISHNU) [ref] .

    Elliptic flow of identified hadrons in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV · Comparison with hydrodynamic calculations · 0.94

  • For example, the extracted specific QGP viscosity $1/4π<(η/s)_QGP< 2.5×1/4π$ from the VISHNU calculations with MC-Glauber and MC-KLN initial conditions [ref] and the implemented $η/s = 0.095$ (with the same bulk viscosity parametrization ) in the MUSIC simulations with the IP-Glasma initialization [ref] are both consistent with this quantitative extracted results from the DUKE-OSU collaborations.

    Collective flow and hydrodynamics in large and small systems at the LHC · Quantitative extractions of the QGP shear and bulk viscosity with massive data evaluations · 0.96

  • Figure [ref] shows the measured pion and proton elliptic flow measured by STAR compared to VISHNU [ref] model calculations.

    Elliptic Flow: A Brief Review · The Perfect Liquid · 0.93

Cites it as support 1 paper
Cites it as background 8 papers

What it cites · 24 sentences without a confident label

Builds on 1 paper
  • In this Appendix we add a few aspects that, due to space limitations, were left out from the discussion in Ref. [ref] of the almost universal dependence of the eccentricity-scaled elliptic flow $v_2/$ on the charged hadron multiplicity density $(1/S)$ on which our extraction of the QGP shear viscosity $(η/s)_QGP$ from RHIC data was based.

    200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid · Section · 0.89

Uses the method or tool of 1 paper
  • In a recent article [ref] we extracted the shear viscosity to entropy density ratio $(η/s)_QGP$ of the quark-gluon plasma (QGP) created in heavy-ion collisions at RHIC by comparing experimental data for the eccentricity-scaled elliptic flow $v_2/$ with calculations performed with [ref] , a hybrid model that describes the QGP stage of the expansion of the collision fireball macroscopically with viscous hydrodynamics (in which $(η/s)_QGP$ enters as an input parameter) but switches to a microscopic description in the late hadronic phase where we solve the Boltzmann equation with [ref] .

    200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid · Introduction · 0.85

Compares with 1 paper
  • Figure [ref] demonstrates excellent agreement between and the experimental data, over the entire range of centralities except for the two most peripheral bins, if we use $(η/s)_QGP0.08$ for MC-Glauber (left column) and $(η/s)_QGP0.16$ for MC-KLN initial conditions (right column).[Fig. [ref] also shows that in both cases the agreement is destroyed when increasing $(η/s)_QGP$ by $1/4π0.08$.] These values are a little smaller than, but consistent with the corresponding values $(η/s)_QGP ≃ 0.1$ for MC-Glauber and $(η/s)_QGP ≃ 0.2$ for MC-KLN that were extracted in [ref] from the $p_T$-integrated, non-flow and fluctuation corrected charged hadron $v_2$.

    200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid · Results: comparison of spectra and $v_2$ to data · 0.96

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.