# Spectra and Elliptic Flow of (Multi)Strange Hadrons at RHIC and LHC within Viscous Hydrodynamics + Hadron Cascade Hybrid Model

Authors: Xiangrong Zhu
Journal: Advances in High Energy Physics
 (2016)

---

## Abstract

Using the (2
+
1)-dimensional ultrarelativistic viscous hydrodynamics + hadron cascade, VISHNU, hybrid model, we study the p
T-spectra and elliptic flow of Λ, Ξ, and Ω in Au + Au collisions at s
N
N = 200 GeV and in Pb + Pb collisions at s
N
N = 2.76 TeV. Comparing our model results with the data measurements, we find that the VISHNU model gives good descriptions of the measurements of these strange and multistrange hadrons at several centrality classes at RHIC and LHC. Mass ordering of elliptic flow v
2 among π, K, p, Λ, Ξ, and Ω are further investigated and discussed at the two collision systems. We find that, at both RHIC and LHC, v
2 mass ordering among π, K, p, and Ω is fairly reproduced within the VISHNU hybrid model, and more improvements are needed to be implemented for well describing v
2 mass ordering among p, Λ, and Ξ.


---

## Body

## 1. Introduction

Ultrarelativistic heavy-ion collisions at the BNL Relativistic Heavy-Ion Collider (RHIC) and CERN Large Hadron Collider (LHC) are used to produce and study a hot and dense medium consisting of strongly interacting quarks and gluons, namely, Quark-Gluon Plasma (QGP), which is expected to exist in the early stage of the universe, and to understand its properties, such as the equation of state (EoS) and transport coefficients. The hadronic interactions are expected to have less influence on the multistrange hadrons, such asΞ and Ω, due to their much smaller hadronic cross-sections. Therefore, final observables of these multistrange hadrons are more sensitive to the early (partonic) stage of the collision. In the past few decades, different aspects of strange and multistrange hadrons have been investigated theoretically [1–13] and experimentally [14–24].Anisotropic flow, which is considered as an evidence for the QGP formation, typically displays the collective behavior of the final emitted particles. It can be characterized by the coefficients of the Fourier expansion of the final particle azimuthal distribution defined as [25](1)
E
d
3
N
d
3
p
=
1
2
π
d
2
N
p
T
d
p
T
d
y
1
+
2
∑
n
=
1
∞
v
n
cos
⁡
n
φ
-
Ψ
n
,where v
n is nth order anisotropic flow harmonic with its corresponding reaction plane angle Ψ
n and φ is the azimuthal angle of the final emitted particles. Recently, the anisotropic flow and other soft hadron data of all charged and identified hadrons at the RHIC and LHC have been studied by many groups within the framework of hydrodynamics [13, 26–36]. VISHNU is a hybrid model [37] for single-shot simulations of heavy-ion collisions, which connects the (2
+
1)-dimensional viscous hydrodynamics with a hadronic afterburner. Employing the VISHNU hybrid model, the specific QGP shear viscosity values of (
η
/
s
)
Q
G
P are extracted from the elliptic flow measurements of charged hadrons with MC-KLN initial conditions [29]. With the extracted (
η
/
s
)
Q
G
P, the VISHNU provides good descriptions of the soft hadron data of π, K, and p at the RHIC and LHC [30]. Compared with other hadrons, anisotropy flow of (multi)strange particles is mainly produced in the QGP stage and less contaminated by the subsequent hadronic interactions. Meanwhile, the p
T-spectra and elliptic flow for Λ, Ξ, and Ω have been measured in the Au + Au collisions at the RHIC [17–20] and Pb + Pb collisions at the LHC [21–23]. Therefore, it is timely to systematically study these strange and multistrange hadrons at RHIC and LHC via the VISHNU hybrid model.In this paper, we investigate thep
T-spectra and elliptic flow v
2 for (multi)strange hadrons in Au + Au collisions at s
N
N = 200 GeV and in Pb + Pb collisions at s
N
N = 2.76 TeV within the viscous hydrodynamic hybrid model VISHNU. The paper is organized as follows. Section 2 briefly introduces the VISHNU hybrid model and its setup in the calculations. Section 3 compares our VISHNU results in Au + Au collisions and Pb + Pb collisions with the measurements from the STAR at RHIC and ALICE at LHC, respectively, mainly including p
T-spectra and differential elliptic flow for Λ, Ξ, and Ω. In Section 4, the mass ordering of elliptic flow among π, K, p, Λ, Ξ, and Ω is studied and discussed at the RHIC and LHC energies. Finally, we summarize our works and give a brief outlook for the future in Section 5.


## 2. Setup of the Calculation

We here give brief descriptions of the inputs and setup ofVISHNU calculations for the soft data at the RHIC and LHC energies. The VISHNU [37] hybrid model consists of two parts, which are the (2
+
1)-dimensional ultrarelativistic viscous hydrodynamics VISH2+1 [39, 40] for the expansion of strongly interacting matter QGP and a microscopic hadronic cascade model (UrQMD) [41, 42] for the hadronic evolution. In the calculations a switching temperature T
s
w of 165 MeV is set for the transition from the macroscopic to microscopic approaches in VISHNU. This switching temperature value is close to the QCD phase transition temperature [43–46]. We input the equation of state (EoS) s95p-PCE [47, 48] for the hydrodynamic evolution above the switching temperature T
s
w. The s95p-PCE, which accounts for the chemical freeze-out at T
c
h
e
m = 165 MeV, was constructed by combing the lattice QCD data at high temperature with a chemically frozen hadron resonance gas at low temperature.Following [29, 30], we input MC-KLN initial conditions [49–51] and start the hydrodynamic simulations at τ
0
=
0.9
f
m
/
c. For improving computational efficiency, we implement single-shot simulations [13, 29, 30, 37, 38, 52] with smooth initial entropy density profiles generated by the MC-KLN model. The smooth initial entropy densities are obtained by averaging over a large number of fluctuating entropy density profiles within a specific centrality class. The initial density profiles are initialized with the reaction plane method, which was once used in [29, 30, 38]. Considering the conversion from total initial entropies to final multiplicity of all charged hadrons, we do the centrality selection through the distribution of total initial entropies that are obtained from the event-by-event fluctuating profiles. Such centrality classification was firstly used by Shen et al. in [53], which is more close to the experimental one defined from the measured multiplicity distributions. The normalization factors for the initial entropy densities in Au + Au collisions and Pb + Pb collisions are, respectively, fixed to reproduce the charged hadron multiplicity density d
N
c
h
/
d
η with 687.4
±
36.6 at the RHIC [54] and 1601
±
60 at the LHC [55] at most central collisions. The λ parameter in the MC-KLN model, which quantifies the gluon saturation scale in the initial gluon distributions [50], is tuned to 0.218 at the RHIC and 0.138 at the LHC for a better description of the centrality dependent multiplicity density for all charged hadrons.In theVISHNU simulations with MC-KLN initial conditions, we set a value of 0.16 for the QGP specific shear viscosity (
η
/
s
)
Q
G
P. Such combined setting in VISHNU calculations once nicely described the elliptic flow of π, K, and p in Au + Au collisions [52] and Pb + Pb collisions [30]. Here, we continue to use it to further study the soft hadron data of strange and multistrange hadrons at both RHIC and LHC. For simplicity of the theoretical calculations, we neglect the bulk viscosity, net-baryon density, and the heat conductivity in the QGP system evolution.


## 3. Spectra and Elliptic Flow

In Figure1, we present the transverse momentum spectra of hadrons Λ, Ξ, and Ω in Au + Au collisions at s
N
N = 200 GeV from the VISHNU hybrid model and compare these results with the STAR measurements. We observe that the VISHNU generally describes the p
T-spectra Ξ but slightly overestimates the production of Ω at all centrality classes. Our VISHNU results of Λ at p
T
<
2
G
e
V
/
c are about 40% lower than Λ from STAR and about 20% lower than Λ
¯ from STAR. This can be understood from the following. In our calculations, the production of Λ is obtained from the original values of strong resonance decays from UrQMD of VISHNU. For the STAR measurements, the Λ spectra are corrected for the feed-down of multistrange baryon weak decays (the feed-down contributions to the Ξ spectra from Ω decays are negligible) [15]. Meanwhile, the STAR Λ spectra are not corrected for the feed-down of σ
0 decays from the channel of Σ
0
→
Λ
+
γ (for Λ
¯, the contribution from Σ
0
¯ is via the channel of Σ
0
¯
→
Λ
¯
+
γ). At LHC, it was found that the contribution from Σ
0 for Λ is about 30% in VISHNU calculations [38]. Furthermore, we notice that STAR measurements of Λ and Ξ
- are slightly larger than their corresponding antiparticles due to nonzero baryon density at this collision energy. In our calculations, however, zero net-baryon density is used, which leads to the same results between these (multi)strange hadrons and their antiparticle partners.Figure 1
Transverse momentum spectra ofΛ(Λ
¯), Ξ
-(Ξ
¯
+), and Ω
-
+
Ω
¯
+ at various centralities in Au + Au collisions at s
N
N = 200 GeV. Experimental data are taken from STAR measurements [15]. Theoretical curves are calculated by the VISHNU hybrid model with the parameters presented in Section 2. From top to bottom the curves correspond to 0–5% (×
1
0
4), 10–20% (×
1
0
3), 20–40% (×
1
0
2), 40–60% (×
1
0
1), and 60–80% (×
1
0
0) centrality, respectively, where the factors in parentheses are the multipliers applied to the spectra for clear separation. The multiplied factor for spectra of Ω at 20–40% is 2
×
1
0
1 instead of 1
×
1
0
1. In VISHNU, the spectra of particles and corresponding antiparticles are the same due to zero net-baryon density used in our calculations. Therefore, only results of particles are shown with solid curves.
(a)
(b)
(c)
(d)
(e)
(f)We also calculate the transverse momentum spectra ofΛ, Ξ, and Ω in Pb + Pb collisions at LHC with our VISHNU hybrid model. The calculations, compared with the measurements from the ALICE Collaboration, are presented in Figure 2. For the ALICE measurements, the production differences between these (multi)strange hadrons and their antiparticles are very small due to very small net-baryon density at the LHC energies. It also shows that the VISHNU generally describes the p
T-spectra of hadrons Λ, Ξ, and Ω at several centrality classes, except for the 60–80% centrality bin. Here, our theoretical calculations of Λ only include strong resonance decays from UrQMD of VISHNU. As a result, they are about 30% lower than the ALICE measurements including contribution from nonweak decays of Σ
0 and Σ
(
1385
) family [21]. The Ω spectra from VISHNU are slightly higher than the experimental data at these centrality classes, as similarly observed in the calculations at the RHIC. Such deviations between theory and experiment are consistent with the model and data differences for the centrality dependent multiplicity shown in [38].Figure 2
Transverse momentum spectra ofΛ
+
Λ
¯, Ξ
-(Ξ
¯
+), and Ω
-(Ω
¯
+) at various centralities in Pb + Pb collisions at s
N
N = 2.76 TeV. Experimental data are taken from LHC [21, 22]. Theoretical curves are calculated with the VISHNU hybrid model with the parameters presented in Section 2. From top to bottom the curves correspond to 0–10% (×
1
0
4), 10–20% (×
1
0
3), 20–40% (×
1
0
2), 40–60% (×
1
0
1), and 60–80% (×
1
0
0) centrality, respectively, where the factors in parentheses are the multipliers applied to the spectra for clear separation. Spectra of Λ start from 0–5% (×
1
0
5) and 5–10% (×
1
0
4), instead of 0–10%. In VISHNU, the spectra of particles and corresponding antiparticles are same due to zero net-baryon density used in our calculations. Therefore, only results of particles are shown with solid curves.
(a)
(b)
(c)
(d)
(e)
(f)From comparisons between our calculations and measurements at the RHIC and LHC, we find that, although theVISHNU cannot fully reproduce the p
T-spectra of these strange and multistrange hadrons in the production amount, it gives nice descriptions of the slopes (distribution shapes) for the spectra of them at various centralities. This can be found from the ratio between model results and data in Figures 1 and 2, which are weakly centrality dependent. Together with the early nice descriptions of the p
T-spectra for π, K, and p [30], it reveals that during the QGP and hadronic evolution the VISHNU hybrid model generates a proper amount of radial flow to push the spectra of various hadrons.Figure3 presents the comparisons of differential elliptic flow of Λ, Ξ, and Ω from the VISHNU model with the STAR measurements in Au + Au collisions at s
N
N = 200 GeV. The theoretical curves are calculated from VISHNU model by using reaction plane initial conditions from the MC-KLN model and (
η
/
s
)
Q
G
P
=
0.16. The experimental data are from the STAR, which are measured with the event plane method [19, 20]. This method covers a fraction of contribution from event-by-event flow fluctuations and nonflow contribution mainly including resonance decays and jets. Compared with STAR measurements, the elliptic flow from the VISHNU generally reproduces the data for Λ at 0–10% and 10–40% and for Ξ and Ω at 0–30% centrality classes. At 40–80% Au + Au collisions, the model gives rough descriptions of the data for Λ at p
T
<
1.0
G
e
V
/
c and for Ξ and Ω at p
T
<
1.5
G
e
V
/
c but underestimates at higher p
T region. Together with the failed descriptions of elliptic flow for charged hadrons and identified lighter hadrons at 40–80% collisions in [28, 29], it reflects that the VISHNU, with the MC-KLN initial conditions, fails to describe the 40–80% semiperipheral collisions. This can be probably interpreted as follows. In 40–80% collisions, the collision lifetimes are shorter, which leave less time to generate the elliptic flow in the fluid dynamic QGP stage, and the highly dissipative effects in hadronic stage cannot compensate for this.Figure 3
Differential elliptic flow of strange hadronsΛ (a) at centralities of 0–10% (squares), 10–40% (diamonds), and 40–80% (circles), and multistrange hadrons Ξ (b) and Ω (c) at centralities of 0–30% (squares) and 30–80% (circles) in Au + Au collisions at s
N
N = 200 GeV. Experimental data (open points) are from STAR [19, 20]. Theoretical results (solid points) are calculated from the VISHNU viscous hydrodynamics hybrid model with the inputs presented in Section 2.
(a)
(b)
(c)In Figure4, we show the differential elliptic flow of Λ, Ξ, and Ω at 10–20%, 30–40%, and 50–60% Pb + Pb collisions, which were first given in [38]. The presented experimental data are measured by the ALICE Collaboration with the scalar product method [23]. The VISHNU theoretical results are calculated with the inputs as presented in Section 2. Figure 4 shows that the VISHNU fairly predicts the elliptic flow data for Λ, Ξ, and Ω at chosen three centrality classes at p
T
<
2
G
e
V
/
c within the statistical error bars. At p
T
>
2
G
e
V
/
c, the descriptions of the elliptic flow for Ξ at 50–60% and for Ω at 30–40% and 50–60% become worse. Together with the worse descriptions of elliptic flow data at high-p
T at the RHIC, as shown in Figure 3, we consider that the viscous corrections probably become too large at high-p
T region, in which the hydrodynamic description in the VISHNU lost its predictive power.Figure 4
Differential elliptic flow of strange hadronsΛ (a), multistrange hadrons Ξ (b) and Ω (c) at 10–20% (squares), 30–40% (diamonds), and 50–60% (circles) centralities in Pb + Pb collisions at s
N
N = 2.76 TeV. Experimental data (open points) are from the ALICE measurements [23]. Theoretical results (solid points) are calculated from the VISHNU hybrid model with the inputs presented in Section 2, which were first shown in [38]
(a)
(b)
(c)


## 4. Mass Ordering of Elliptic Flow

It is widely accepted that the characteristic mass ordering of differential elliptic flow among various identified hadrons at low-p
T reflects the interplay between radial and elliptic flow, providing more insights into the properties of the QGP fireball. The radial flow creates a depletion in the particle p
T-spectrum at low values, which increases with increasing particle mass. This leads to heavier particles having a smaller v
2 compared to lighter ones at a given value of p
T, giving a mass ordering of the p
T dependent elliptic flow below 1.5
–
2
G
e
V
/
c. Such v
2 mass ordering has been discovered in the experiments at both the RHIC and LHC [19, 23, 56–58], which has also been studied within the framework of hydrodynamics [8, 13, 59–61] and blastwave model [8, 62].Here we investigate the mass ordering of the elliptic flow among various identified hadronsπ, K, p, Λ, Ξ, and Ω in Au + Au collisions and Pb + Pb collisions. For clear presentations, the experimental data and our VISHNU results are plotted in separate panels at 0–80% at the RHIC and 30–40% at the LHC. Together with the calculations of elliptic flow of identified hadrons in [29, 30] and in this paper, we find that the VISHNU generally describes the v
2
(
p
T
) for different identified hadrons at several centrality classes. However, our theoretical results presented in Figure 5, compared with the measurements at RHIC and LHC, show that the VISHNU fairly describes the mass ordering among π, K, p, and Ω. But it is hard to see mass ordering among p, Λ, and Ξ clearly due to their elliptic flow being almost identical at low p
T region. This underprediction for proton leads to an inverse v
2 mass ordering between p and Λ. Effects of hadronic rescattering on elliptic flow are seen to be particle specific, depending on their scattering cross-sections that couple them to the medium [59]. Compared with nonstrange hadrons, the multistrange hadrons are less affected on their differential elliptic flow due to their smaller scattering cross-sections. Therefore, reevaluating the hadronic cross-sections in the UrQMD is helpful to improve the description of elliptic flow of various hadron species. Meanwhile, an initial flow could enhance the radial flow in the hadronic stage, which is also expected to improve the description of mass ordering within the framework of the hybrid model.Figure 5
Differential elliptic flow ofπ, K, p, Λ, Ξ, and Ω at centrality 0–80% in Au + Au collisions at s
N
N = 200 GeV ((a-b), (a) data, (b) VISHNU calculations), and at centrality 30–40% in Pb + Pb collisions at s
N
N = 2.76 TeV ((c-d), (c) data, (d) VISHNU calculations). The measurements of elliptic flow of various hadron species are taken from STAR [15, 20] and ALICE [23].
(a)
(b)
(c)
(d)


## 5. Summary and Outlook

In summary, we studied thep
T-spectra and elliptic flow of strange and multistrange hadrons in Au + Au and Pb + Pb collisions within the VISHNU hybrid model. At both collision systems, we found that, with MC-KLN initial conditions, η
/
s
=
0.16, and other inputs, VISHNU generally describes the p
T-spectra of strange hadron Λ and multistrange hadrons Ξ at some centrality classes but slightly overestimates for Ω at chosen centrality classes. In spite of the normalization issues, the VISHNU well produces the spectra slopes of these three hadrons at chosen centralities. By comparing the elliptic flow of Λ, Ξ, and Ω in Au + Au collisions and Pb + Pb collisions from VISHNU model with the STAR and ALICE measurements, we found that the VISHNU generally describes the elliptic flow except at 40–80% semiperipheral collisions. The failed descriptions at 40–80% collisions is probably due to shorter lifetimes of these collisions, which leave less time to generate the elliptic flow in the fluid dynamic QGP stage.We also compared the mass ordering ofv
2 among hadrons π, K, p, Λ, Ξ, and Ω from VISHNU calculations with the STAR and ALICE measurements. The comparisons showed that the elliptic flow mass ordering among various hadron species is not fully described at both RHIC and LHC. The VISHNU fairly describes the mass ordering of v
2 among π, K, p, and Ω but fails to reproduce the mass ordering among p, Λ, and Ξ due to slight underpredictions of the elliptic flow of protons. The effects from the initial flow and/or improved UrQMD hadronic cross-sections may solve this issue within the framework of VISHNU, which should be investigated in the near future.



---

*Source: 4236492-2016-09-07.xml*