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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v2
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:1101.4638 by the authors.
References
How the literature uses it
48 arXiv papers cite this in 61 sentences
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
-
Simulations with an evolution-averaged $η/s$ based on comparisons with combined experimental data from top-RHIC and LHC energies [ref] extracted a possible range of $1<(η/s)/(1/4π)<5$.
Quark Flavour Dependence of the Shear Viscosity in a Quasiparticle Model · Introduction · 0.91
-
However, in [ref] the initial time $t_i=1$ fm/c used at the LHC was bigger than at RHIC where $t_i=0.6$ fm/c [ref] , based on the requirements of the hydrodynamic simulations to fit the low $p_T$ elliptic flow data.
Holographic light quark jet quenching at RHIC and LHC via the shooting strings · Calculation of $R_AA$ · 0.80
-
Some characteristic values of the relevant parameters for the quark-gluon plasma are $T ∼330 MeV$, $τ_th ∼1 fm/c$, $s ∼85 fm^-3$ [ref] and $η/s ∼0.15 ħ/k_B$ [ref] .
Thermalization, Viscosity and the Averaged Null Energy Condition · The quark-gluon plasma · 0.96
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
-
The larger integrated charged hadron elliptic flow at the LHC implies a larger total momentum anisotropy [ref] and suggests a higher efficiency of the QGP fluid for converting the initial spatial deformation of the collision fireball into anisotropic collective flow.
Elliptic flow in 200 A GeV Au+Au collisions and 2.76 A TeV Pb+Pb collisions: insights from viscous hydrodynamics + hadron cascade hybrid model · Introduction · 0.98
Cites it as background 8 papers
-
Based on these experimental data, lots of studies on properties of hadronic $p_T$ distribution are carried out, which greatly improve people's understandings on property of the created hot nuclear matter and the mechanism of hadron production in relativistic heavy-ion collisions [ref] .
Averaged transverse momentum correlations of hadrons in relativistic heavy-ion collisions · Introduction · 0.81
-
Most hydrodynamic calculations at these energies therefore assume boost invariance and employ a 2+1D hydrodynamic evolution when investigating transverse flow and it's fluctuations [ref] .
Rapidity dependence of mean transverse momentum fluctuation and decorrelation in baryon-dense medium · Introduction · 0.87
-
Furthermore, various hybrid models have recently been developed, where the hydrodynamical evolution is coupled with a hadronic cascade afterburner [ref] .
Systematics of the charged-hadron P_T spectrum and the nuclear suppression factor in heavy-ion collisions from sqrt{s}=200 GeV to sqrt{s} =2.76 TeV · Introduction · 0.95
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.