Abstract

Using the ()-dimensional ultrarelativistic viscous hydrodynamics + hadron cascade, VISHNU, hybrid model, we study the -spectra and elliptic flow of , , and in Au + Au collisions at = 200 GeV and in Pb + Pb collisions at = 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 among , , , , , and are further investigated and discussed at the two collision systems. We find that, at both RHIC and LHC, mass ordering among , , , and is fairly reproduced within the VISHNU hybrid model, and more improvements are needed to be implemented for well describing mass ordering among , , and .

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 [113] and experimentally [1424].

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]where is th order anisotropic flow harmonic with its corresponding reaction plane angle 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, 2636]. VISHNU is a hybrid model [37] for single-shot simulations of heavy-ion collisions, which connects the ()-dimensional viscous hydrodynamics with a hadronic afterburner. Employing the VISHNU hybrid model, the specific QGP shear viscosity values of are extracted from the elliptic flow measurements of charged hadrons with MC-KLN initial conditions [29]. With the extracted , the VISHNU provides good descriptions of the soft hadron data of , , and 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 -spectra and elliptic flow for , , and have been measured in the Au + Au collisions at the RHIC [1720] and Pb + Pb collisions at the LHC [2123]. 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 the -spectra and elliptic flow for (multi)strange hadrons in Au + Au collisions at = 200 GeV and in Pb + Pb collisions at = 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 -spectra and differential elliptic flow for , , and . In Section 4, the mass ordering of elliptic flow among , , , , , 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 of VISHNU calculations for the soft data at the RHIC and LHC energies. The VISHNU [37] hybrid model consists of two parts, which are the ()-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 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 [4346]. We input the equation of state (EoS) s95p-PCE [47, 48] for the hydrodynamic evolution above the switching temperature . The s95p-PCE, which accounts for the chemical freeze-out at = 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 [4951] and start the hydrodynamic simulations at . 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 with at the RHIC [54] and 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 at the RHIC and at the LHC for a better description of the centrality dependent multiplicity density for all charged hadrons.

In the VISHNU simulations with MC-KLN initial conditions, we set a value of 0.16 for the QGP specific shear viscosity . Such combined setting in VISHNU calculations once nicely described the elliptic flow of , , and 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 Figure 1, we present the transverse momentum spectra of hadrons , , and in Au + Au collisions at = 200 GeV from the VISHNU hybrid model and compare these results with the STAR measurements. We observe that the VISHNU generally describes the -spectra but slightly overestimates the production of at all centrality classes. Our VISHNU results of at 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 decays from the channel of (for , the contribution from is via the channel of ). At LHC, it was found that the contribution from 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.

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 -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 and 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].

From comparisons between our calculations and measurements at the RHIC and LHC, we find that, although the VISHNU cannot fully reproduce the -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 -spectra for , , and [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.

Figure 3 presents the comparisons of differential elliptic flow of , , and from the VISHNU model with the STAR measurements in Au + Au collisions at = 200 GeV. The theoretical curves are calculated from VISHNU model by using reaction plane initial conditions from the MC-KLN model and . 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 and for and at but underestimates at higher 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.

In Figure 4, 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 within the statistical error bars. At , 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- at the RHIC, as shown in Figure 3, we consider that the viscous corrections probably become too large at high- region, in which the hydrodynamic description in the VISHNU lost its predictive power.

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- 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 -spectrum at low values, which increases with increasing particle mass. This leads to heavier particles having a smaller compared to lighter ones at a given value of , giving a mass ordering of the dependent elliptic flow below . Such mass ordering has been discovered in the experiments at both the RHIC and LHC [19, 23, 5658], which has also been studied within the framework of hydrodynamics [8, 13, 5961] and blastwave model [8, 62].

Here we investigate the mass ordering of the elliptic flow among various identified hadrons , , , , , 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 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 , , , and . But it is hard to see mass ordering among , , and clearly due to their elliptic flow being almost identical at low region. This underprediction for proton leads to an inverse mass ordering between 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.

5. Summary and Outlook

In summary, we studied the -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, , and other inputs, VISHNU generally describes the -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 of among hadrons , , , , , 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 among , , , and but fails to reproduce the mass ordering among , , 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.

Competing Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The author gratefully thanks Huichao Song and Hao-jie Xu for fruitful discussions and critical reading of the draft. This work was supported in part by the NSFC and the MOST under Grants no. 11435001 and no. 2015CB856900 and the China Postdoctoral Science Foundation under Grant no. 2015M570878. The author especially acknowledges extensive computing resources provided by Tianhe-1A from the National Supercomputing Center in Tianjin, China.