Hydrodynamic modeling of a pureglue initial scenario in highenergy hadron and heavyion collisions
Abstract
Partonic matter produced in the early stage of ultrarelativistic nucleusnucleus collisions is assumed to be composed mainly of gluons, and quarks and antiquarks are produced at later times. The comparable hydrodynamic simulations of heavyion collisions for (2+1)flavor and YangMills equations of state performed by using three different hydrodynamic codes are presented. Assuming slow chemical equilibration of quarks, the spectra and elliptic flows of thermal dileptons and photons are calculated for central Pb+Pb collisions at the LHC energy of TeV. It is shown that a suppression of quarks at early times leads to a significant reduction of the yield of the thermal dileptons, but only to a rather modest suppression of the distribution of direct photons. It is demonstrated that an enhancement of photon and dilepton elliptic flows might serve as a promising signature of the pureglue initial state. Calculations based on Bjorken hydrodynamics suggest that collisions of small systems at intermediate energies available at RHIC or future FAIR facilities may show stronger effects associated with initial pure gluodynamic evolution.
Hydrodynamic modeling of a pureglue initial scenario in highenergy hadron and heavyion collisions
V. Vovchenko, LongGang Pang, H. Niemi, Iu.A. Karpenko, M.I. Gorenstein, L.M. Satarov, I.N. Mishustin, B. Kämpfer, H. Stoecker^{†}^{†}thanks: Speaker.
Frankfurt Institute for Advanced Studies, D60438 Frankfurt, Germany
Institut für Theoretische Physik, Goethe Universität, D60438 Frankfurt, Germany
Department of Physics, Taras Shevchenko National University of Kiev, 03022 Kiev, Ukraine
Bogolyubov Institute for Theoretical Physics, 03680 Kiev, Ukraine
INFN  Sezione di Firenze, I50019 Sesto Fiorentino (Firenze), Italy
National Research Center ”Kurchatov Institute”, 123182 Moscow, Russia
HelmholtzZentrum DresdenRossendorf, D01314 Dresden, Germany
Technische Universität Dresden, Institut für Theoretische Physik, D01062 Dresden, Germany GSI Helmholtzzentrum für Schwerionenforschung GmbH, D64291 Darmstadt, Germany
Email: H.Stoecker@gsi.de
\abstract@cs
1 Introduction
There are many different approaches to describe the initial stage of nucleusnucleus (A+A) collisions. Usually it is assumed that strong nonequilibrium effects take place only during a very short proper time interval , where is the socalled gluon saturation scale [1]. The idea that the gluonic components of colliding nucleons dominate in high energy collisions was originally put forward in Ref. [2], and it was motivated by the fact that the perturbative gluongluon cross sections are larger than the quarkantiquark ones. A twostep equilibration of QGP was proposed in [3, 4, 5] assuming that the gluon thermalization is accomplished already at the early proper time , while the quarkantiquark chemical equilibration proceeds until later times (according to Ref. [6], ). Possible signatures of quark undersaturation in high energy A+A collisions were considered by different authors, see, e.g., Refs. [4, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17].
Recently, the pure glue initial scenario of Pb+Pb collisions at Large Hadron Collider (LHC) energies was discussed in Refs. [18, 19]. Below we describe how a pureglue initial scenario for heavyion collisions can be modeled using the (2+1)–dimensional boostinvariant hydrodynamics. In particular, we demonstrate the difference between hydrodynamic evolution of the full QCD matter and of the pure glue matter by performing simulations within three different hydrodynamic codes. We also introduce the timedependent quarkantiquark fugacity in order to describe the QGP evolution in the absence of the full chemical equilibrium. The new aspects of the present study include constructing the equation of state for chemically nonequilibrated QCD matter via an interpolation of the lattice data, and the solution of hydro equations with a timedependent equation of state.
A significant emphasis is put on electromagnetic probes (thermal photons and dileptons), which may carry an important information about the deconfined phase. This problem has been repeatedly addressed in the literature, see, e.g., [9, 10, 14, 15, 16], however, a definitive conclusion about the role of chemically nonequilibrium evolution is still not obtained. We also analyze the impact of chemical nonequilibrium effects on the momentum anisotropy of photons and dileptons. A more detailed version of this study is given in Ref. [20].
In addition, we simulate the evolution of system created in A+A and p+p collisions at lower energies, by employing the simplified Bjorken hydrodynamics. There we consider only two limiting cases: the fully equilibrated (2+1)flavor evolution and the pure gluodynamic evolution.
2 Formulation of the model
The equations of the longitudinally boostinvariant (2+1)–dimensional ideal hydrodynamics can be written as ()
(2.0) 
where is the energymomentum tensor, is the 4velocity, and are the local restframe energy density and pressure, respectively, and is the metric tensor. with oriented along the beam axis. Below we use the curvilinear lightcone coordinates , where is the proper time and is the spacetime rapidity. In the case of the longitudinal boost–invariant (2+1)–dimensional flow one can represent the fluid’s fourvelocity as [21] , where is the transverse velocity in the symmetry plane and stands for the transverse Lorentz factor. To solve Eq. (2) one needs the equation of state (EoS), namely a relation connecting and . For chemically nonequilibrated matter considered in this paper, , where is the (anti)quark fugacity. In this work we use a simple analytic parametrization for as a function of the proper time (see below Eq. (2)).
Additionally, we calculate the spacetime dependence of the local proper time of a fluid cell element by solving the equation
(2.0) 
where the parameter corresponds to initial longitudinal proper time of the hydrodynamic expansion. Equation (2) is solved simultaneously with Eq. (2). In general, is found to be smaller than the ’global’ time due to the presence of nonzero transverse flow. In the limiting case of the onedimensional longitudinal Bjorken expansion [22], one has and, consequently, .
In order to construct the equation of state we use the lattice QCD calculations for the EoS of the strongly interacting matter in two limiting cases: 1) the chemically equilibrated QCD matter [23, 24], 2) the SU(3) gluodynamics without (anti)quarks [25, 26]. In the following we denote these cases as FQ (Full QCD) and PG (Pure Glue), respectively. The FQ case corresponds to the (2+1)flavour QCD calculations which predict the crossovertype transition at . The PG calculation provides a firstorder deconfinement phase transition at . The temperature dependencies of the pressure and energy density for FQ and PG scenarios are depicted in Fig. 1. Larger values of and in the FQ calculation appear due to the contribution of quarkantiquark degrees of freedom. Note that there is a discontinuity of at in the PG case. We note that very small values of and at in the PG matter originate from large masses of glueballs () which are the constituents of the confined phase [26].
The suppression of the quark and antiquark densities at given temperature is characterized by the (anti)quark fugacity (for details, see Ref. [27]). Generalizing the lattice EoS for the chemically nonequilibrium case with is not a straightforward task. We construct the and as functions of and by using the linear interpolation (LI) between the the PG and FQ cases:
(2.0)  
(2.0) 
After excluding the temperature variable in Eqs. (2)–(2), one gets the relation which is used in hydrodynamic simulations. The limits and correspond to the thermodynamic functions of the PG and FQ matter, respectively. Note that the linear – dependence of and is a characteristic feature of the ideal gas of massless gluons and (anti)quarks studied in Refs. [19, 27].
Using Eqs. (2) and (2) and basic thermodynamic identities, one can calculate the total density of (anti)quarks and the entropy density . The following relations are obtained
(2.0)  
(2.0) 
We assume that at the initial (anti)quark densities vanish in all cells and gluons are in thermal and chemical equilibrium. Similarly to Refs. [19, 27] we postulate that is an explicit function of the local proper time which increases from at to at . The following simple parametrization is used:
(2.0) 
where is a model parameter characterizing the chemical equilibration time of (anti)quarks. There are different estimates for in the literature ranging from [28] to [6]. Note that corresponds to the instantaneous chemical equilibration of quarks and gluons.
3 Hydrodynamic simulations for YangMills and equilibrium QCD equations of state
3.1 CLVisc hydro simulation
As a first step we perform the hydrodynamic simulations of heavyion collisions for two limiting cases: the (2+1)flavor full QCD matter and the quarkless pure YangMills (YM) matter. The former case corresponds to fm/ while the latter implies . Such a comparison is rather instructive and may indicate the possible effects of the pure glue scenario that one can later study in a more realistic approach. For this purpose we use three different hydrodynamic codes. First, we use the CLVisc hydro code [29] and simulate the central Pb+Pb collisions. The initial entropy deposition at fm/ in the transverse plane is proportional to the number of wounded nucleons calculated in the optical Glauber model. The normalization constant is fixed in order to reproduce the maximum initial energy density of 166 GeV/fm, which is reached in the central cell. The resulting spacetime profile of the temperature in the plane is shown in Fig. 2. Here is the transverse coordinate in the reaction plane.
For the equation of state of the full QCD the resulting temperature profile, shown in the right panel of Fig. 2 is rather typical for hydrodynamic simulation in ultrarelativistic collisions: the system smoothly and rather rapidly cools down to the “freezeout” temperatures of MeV within the time interval of fm/. The situation, however, is very different for the pure YM scenario. The system cools down to the critical temperature of MeV initially, then spends a significant amount of time in the mixed phase of the deconfinement phase transition, and finally enters the phase of quickly decaying glueballs at MeV. The resulting lifetime of the system appears to be much longer as compared to the full QCD case: it can be about 23 times larger and is a consequence of undergoing the long deconfinement phase transition during the system evolution.
3.2 SHASTAbased hydro simulation
A similar calculation of the temperature profile performed within the SHASTAbased hydrodynamical solver introduced in Refs. [30, 31] is shown in Fig. 3, where the hydrodynamical calculation using the pure YM EoS is compared to the one using a parametrization [32] of (2+1) flavor QCD EoS. For this calculation the initial energy density profile is taken to be proportional to the density of binary nucleon collisions, but the normalization is fixed such that the final pion multiplicity is approximately the same in both the YM and full QCD calculations. The impact parameter is fm, which roughly corresponds to the % centrality class. Despite some differences in the initial conditions the temperature profiles shown in Fig. 3 are generally consistent with the ones depicted in Fig. 2, which were obtained with the CLVisc hydro code. In Fig. 4a we show the transverse momentum spectra and in Fig. 4b the elliptic flow coefficient of the positively charged pions for both pure YangMills and the (2+1) flavor QCD equations of state. The pion spectra in the pure SU(3) case are calculated by constructing an equation of state of free (noninteracting) glueballs, such that pressure as a function of temperature, as well as its slope is approximately the same around MeV in both pure YangMills and glueball EoS. The number of glueball states is the same as listed in Ref. [33], but the masses of the states are quite freely modified in order to get a good match to the pure YM EoS. The glueball states then decay into pions and kaons.
3.3 vHLLE hydro simulation
It is also instructive to perform a comparable hydrodynamic simulation of heavyion collisions with the YM and the (2+1)flavor QCD EoS within another hydro code, namely the vHLLE hydrodynamic solver [34]. The resulting temperature profiles in the 040% central Pb+Pb collisions are shown in Fig. 5. One can compare these calculations with the corresponding results from the two other codes shown in Figs. 2 and 3. While the initial conditions are slightly different in these three simulations they all show a consistent physical picture: the evolution of the YM matter in heavyion collisions is very different from that of a fully equilibrated QCD matter. The YM matter evolves for a much longer time and spends a significant portion of the spacetime evolution in the region of mixed phase.
4 Numerical results for PbPb collisions at LHC with pure glue initial state
In a more realistic scenario one has to take into account the gradual appearance of (anti)quarks during the system evolution, and also retain the consistence with the available experimental data. We model the Pb+Pb collisions at the LHC with the c. m. energy per nucleon pair TeV. In our calculations we choose fm/ as the initial longitudinal proper time of the hydrodynamic evolution. It is assumed that there is no initial transverse flow, i.e., , and that the initial energy density profile is proportional to the linear combination of the transverse distributions of wounded nucleons and of binary collisions taken from the eventaveraged Monte Carlo Glauber model as implemented in the GLISSANDO code [35]. The coefficient of proportionality in the initial profile is fixed to reproduce the observed hadron spectra within the simulation assuming chemical equilibrium with the full QCD EoS for a corresponding centrality interval (see Ref. [36] for details). We use the same initial energy density profile in the present calculations for the chemical nonequilibrium case. It is also assumed that the initial state is purely gluonic, i.e. that initially the fugacity of (anti)quarks is zero, and that gluons are in chemical equilibrium at .
Equations (2) and (2) are solved using the (2+1)–dimensional version of the vHLLE hydrodynamic code [34]. We consider the 0 –20 % and 20 –40 % central Pb+Pb collisions at LHC.
The contour plot of the quark fugacity in the plane is given in Fig. 6 a. The dashed line shows the isotherm which presumably corresponds to the hadronization hypersurface. One can see that typical lifetimes of the deconfined phase in the considered reaction do not exceed . In Fig. 6 a one observes that deviations from chemical equilibrium () may survive up to the hadronization stage. As discussed in Ref. [27] this may lead to a suppression of (anti)baryontopion ratios observed [37] for the considered reaction. Note that evolves with in Fig. 6 a even in the large regions where there is practically no matter. This behavior follows from applying Eq. (2). In reality, of course, the values of for these very dilute and cold fluid elements are irrelevant and should be ignored.
Figure 6 b shows the contour plot of the temperature in the plane . The solid and dashed curves correspond to and , respectively. One can see that the chemically undersaturated matter is hotter as compared to the equilibrium case (), especially at the earlier times ^{1}^{1}1 Note that in both cases we take the same profile of the energy density at . . This is a consequence of reduced number of degrees of freedom in such a medium. According to Fig. 6 b, typical lifetimes of the mixed phase are rather short, they do not exceed . This is at variance with calculations in the (1+1) dimensional hydrodynamics which predict [38] much larger lifetimes of the mixed phase within the chemically equilibrated bag model. Our nunerical results show that about 25% of the total final entropy is generated during the ideal hydro evolution of the chemically undersaturated QCD matter, in line with results of the semianalytic calculation in Ref. [27].
5 Direct photon and thermal dilepton emission at LHC
The emission of direct ^{2}^{2}2 By direct photons we denote the ’noncocktail’ photons i.e. those which are not produced in decays of , and mesons in the final stage of the reaction. photons from expanding matter created in relativistic A+A collisions has several components [39, 40]: a) ’prompt’ photons from binary collisions of initial nucleons, b) ’thermal’ photons from the hightemperature deconfined phase, c) direct photons from the lowtemperature hadronic phase. The contribution of prompt photons becomes dominant at large transverse momenta. As we will see below, this greatly reduces the sensitivity of photon spectra to chemical nonequilibrium effects. However, the situation with transverse flows of photons is different because of low azimuthal anisotropy of prompt photons. Note that the ALICE experiments [41] reveal large elliptic flows of direct photons, which still can not be explained in the chemically equilibrium scenario [39].
Figure 7 a shows our results for the thermal photon spectrum in the 0–20 % and 20–40 % central Pb+Pb collisions at calculated with the cutoff temperature of MeV. We apply the parametrization QGPUA [20] of the photon production rate in a chemically nonequilibrium QGP. The low spectrum looks very similar in all scenarios, while there is a sizable difference at high . The spectra of direct photons, i.e. those that additionally include the prompt photons, is depicted in Fig. 7 b. The inclusion of prompt photons makes the difference between different scenarios at high rather small.
The direct photon production in Pb+Pb collisions at LHC has been considered in various theoretical models ([43, 44, 39, 45]. As noted in Ref. [42], the present uncertainties in the heavyion photon data at LHC do not allow to conclusively discriminate between the various scenarios.
The photon elliptic flow is calculated by
(5.0) 
The photon spectrum, entering this equation includes both thermal and prompt components. We assume that prompt photons are azimuthally symmetric. Therefore, they contribute only to the denominator of Eq. (5) reducing at large . The results of the calculations of the photon elliptic flow are shown in Fig. 8. In the pureglue scenario the momentum anisotropy is significantly enhanced for the high thermal photons. As seen from Fig. 8 b, the inclusion of prompt photons notably decreases the effect.
We also analyze spectra of thermal dileptons produced in nuclear collisions at the LHC energies. It is assumed that dileptons are produced in the processes. The additional factor is introduced into the dilepton production rate which takes into account the quark suppression in a chemically nonequilibrium QGP. The results of calculating the dilepton mass spectrum in % and % central Pb+Pb collisions at TeV are shown in Fig. 9 for the cutoff temperature . One can see that the initial quark suppression leads to a strong reduction of the dilepton yield at . Note that we do not include contributions of hard (DrellYan) dileptons [40] produced in binary collisions of initial nucleons.
The elliptic flows of thermal dileptons in the same reaction are shown in Fig. 10 for several values of . Similar to direct photons we predict a strong enhancement of the dilepton elliptic flow as compared to the equilibrium scenario ().
6 Pure glue scenario at lower energies
The presented calculations, as well as their comparison to experimental data, show that the pure glue initial scenario does not spoil the existing agreement of hydro simulations with the data at LHC energies. On the other hand, e.g. as seen from Fig. 6b, the effect of initial pure gluodynamic evolution, in particular the firstorder phase transition at MeV, is almost completely washed out in most observables. It is evident that the gluodynamical effects can be more pronounced at lower collision energies where the initial temperature are much closer to the critical temperature MeV of the deconfinement phase transition in the YangMills theory.
In order to investigate this aspect qualitatively, we use the onedimensional Bjorken model and consider just the two limiting cases: the fully equilibrated (2+1)flavor evolution and the pure gluodynamic evolution. The corresponding equations of state are shown in Fig. 1. We also consider both the heavyion () and the smaller systems. Unlike for LHC energy, we adopt fm/ in the present analysis. The initial entropy density at a given collision energy is estimated by using the available data on pion multiplicity. In order to estimate the uncertainties of the obtain results the transverse radius is varied, namely fm for collisions and fm for p+p collisions.
The values of extracted initial temperature for and are shown in Fig. 11. The resulting initial temperature in pure SU(3) case at collision energies GeV is rather close to (slightly above) the phase transition temperature of MeV. This especially concerns the p+p collisions. Thus, these energies look promising for observing the effects of the phase transition in the pure glue scenario.
At smaller collision energies, however, the effect of constituent quarks in colliding nucleons becomes more and more important. In that sense it becomes difficult to consider the matter created in such systems as netbaryon free. On the other hand, the collisions at these energies may be more promising with regards to the pure glue scenario. The PANDA experiment at FAIR [46] is, in principle, capable of running such an experiment in the future, by colliding the proton and (anti)proton beams at GeV. The dependence of the temperature for such a configuration within Bjorken hydrodynamics is depicted in Fig. 12. In the pure glue scenario the matter can enter the mixed phase at MeV at the early stage of evolution and spend a significant amount of time there. In this regard, the future PANDA experiment looks promising in the search of new exotic states of matter.
7 Summary
In summary, the aspects and the calculation results of the hydrodynamic modeling of the pure glue initial scenario for hadron and heavyion collisions has been presented. The calculations performed within three different hydro codes all show a consistent physical picture: the evolution of the pure glue matter in heavyion collisions is very different from that of a fully equilibrated QCD matter. The pure glue matter evolves for a much longer time and spends a significant portion of its spacetime evolution in the region of mixed phase. The calculations performed for the LHC energy show that the pure glue initial scenario does not spoil the existing agreement of hydro with the data, in particular with regards to the direct photon yield. The suppression of dilepton yield and the enhancement of its momentum anisotropy may serve as the promising observables to determine the properties of the initial state in Pb+Pb collisions at LHC. Estimates based on the Bjorken hydrodynamics imply that collisions of smaller systems at smaller collision energies, for instance the collisions at the possible PANDA energy of GeV, are promising for searching new exotic states of matter.
Acknowledgements
This work was supported by HIC for FAIR within the LOEWE program of the State of Hesse. H.St. acknowledges the support through the Judah M. Eisenberg Laureatus Chair at Goethe University. V.V. acknowledges the support from HGSHIRe for FAIR. The work of M.I.G. was supported by the GoalOriented Program of the National Academy of Sciences of Ukraine and the European Organization for Nuclear Research (CERN), Grant CO132016.
References
 [1] L. V. Gribov, E. M. Levin, M. G. Ryskin, Semihard processes in QCD, Phys. Rept. 100 (1983) 1.
 [2] L. van Hove and S. Pokorski, Highenergy hadronhadron collisions and internal hadron structure, Nucl. Phys. B86 (1975) 243.
 [3] S. Raha, Dilepton, diphoton and photon production in preequilibrium, Physica Scripta T32 (1990) 180.
 [4] E. Shuryak, Two stage equilibration in highenergy heavy ion collisions, Phys. Rev. Lett. 68 (1992) 3270.
 [5] J. Alam, S. Raha, and B. Sinha, Successive equilibration in quark  gluon plasma, Phys. Rev. Lett. 73 (1994) 1895.
 [6] Z. Xu and C. Greiner, Thermalization of gluons in ultrarelativistic heavy ion collisions by including threebody interactions in a parton cascade, Phys. Rev. C 71 (2005) 064901 [hepph/04062778].
 [7] T. S. Biró, E. van Doorn, B. Müller, M. H. Thoma, and X.N. Wang, Parton equilibration in relativistic heavy ion collisions, Phys. Rev. C 48 (1993) 1275 [nuclth/9303004].
 [8] M. Strickland, Thermal photons and dileptons from nonequilibrium quark  gluon plasma, Phys. Lett. B331 (1994) 245.
 [9] B. Kämpfer and O. P. Pavlenko, Photon production in an expanding and chemically equilibrating gluon enriched plasma, Z. Phys. C 62 (1994) 491.
 [10] B. Kämpfer, O. P. Pavlenko, A. Peshier, and G. Soff, Dilepton production in a chemically equilibrating, expanding, and hadronizing quark  gluon plasma, Phys. Rev. C 52 (1995) 2704.

[11]
C. T. Traxler and M. H. Thoma, Photon emission from a parton gas at chemical nonequilibrium,
Phys. Rev. C 53 (1996) 1348 [hepph/9601254].  [12] D. M. Elliott, D. H. Rischke, Chemical equilibration of quarks and gluons at RHIC and LHC energies, Nucl. Phys. A671 (2000) 583 [nuclth/9908004].
 [13] D. Dutta, S. V. S. Sastry, A. K. Mohanty, K. Kumar, Hard photon production from unsaturated quarkgluon plasma at twoloop level, Nucl. Phys. A710 (2002) 415 [hepph/0104134].
 [14] F. Gelis, H. Niemi, R. V. Ruuskanen, S. S. Räsänen, Photon production from nonequilibrium QGP in heavy ion collisions, J. Phys. G30 (2004) S1031 [nuclth/0403040].
 [15] F.M. Liu and S.X. Liu, Quarkgluon plasma formation time and direct photons from heavy ion collisions, Phys. Rev. C 89 (2014) 034906 [arXiv:1212.6587].
 [16] A. Monnai, Thermal photon v2 with slow quark chemical equilibration, Phys. Rev. C 90 (2014) 021901 [arXiv:1403.4225].
 [17] P. Moreau, O. Linnyk, W. Cassing, and E.L. Bratkovskaya, Examinations of the early degrees of freedom in ultrarelativistic nucleusnucleus collisions, Phys. Rev. C 93 (2016) 044916 [arXiv:1512.02875].

[18]
H. Stoecker et al., Glueballs amass at RHIC and LHC! The early quarkless firstorder phase transition at T=270 MeV  from pure YangMills glue plasma to Hagedorn glueball states,
J. Phys. G43 (2016) 015105 [arXiv:1509.00160]. 
[19]
H. Stöcker et al., Undersaturation of quarks at early stages of relativistic nuclear collisions:
The hot glue initial scenario and its observable signatures, Astron. Nachr. 336 (2015) 744 [arXiv:1509.07682]. 
[20]
V. Vovchenko, Iu. A. Karpenko, M. I. Gorenstein, L. M. Satarov,
I.N. Mishustin, B. Kämpfer,
and H. Stoecker, Electromagnetic probes of a pureglue initial state in nucleusnucleus collisions
at LHC, Phys. Rev. C94 (2016) 024906 [arXiv:1604.06346].  [21] K. Kajantie, M. Kataja, L. M. McLerran, and P. V. Ruuskanen, Transverseflow effects in dilepton emission, Phys. Rev. D34 (1986) 811.
 [22] J. D. Bjorken, Highly relativistic nucleusnucleus collisions: The central rapidity region, Phys. Rev. D27 (1983) 140.
 [23] S. Borsanyi et al., Full result for the QCD equation of state with 2+1 flavors, Phys. Lett. B730 (2014) 99 [arXiv:1309.5258].
 [24] A. Bazavov [HotQCD Collaboration], The QCD equation of state, Nucl. Phys. A931 (2014) 867.
 [25] G. Boyd et al., Equation of state for the SU(3) gauge theory, Phys. Rev. Lett. 75 (1995) 4169 [heplat/9506025].
 [26] S. Borsanyi et al., Precision SU(3) lattice thermodynamics for a large temperature range, JHEP 1207 (2012) 056 [arXiv:1204.6184].
 [27] V. Vovchenko, M. I. Gorenstein, L. M. Satarov, I. N. Mishustin, L. P. Csernai, I. Kisel, and H. Stöcker, Entropy production in chemically nonequilibrium quarkgluon plasma created in central Pb+Pb collisions at energies available at the CERN Large Hadron Collider, Phys. Rev. C 93 (2016) 014906 [arXiv:1510.01235].
 [28] M. Ruggieri, S. Plumari, F. Scardina, and V. Greco, Quarks production in the quarkgluon plasma created in relativistic heavy ion collisions, Nucl. Phys. A941 (2015) 201 [arXiv:1502.04596].
 [29] L. G. Pang, Y. Hatta, X. N. Wang, and B. W. Xiao, Analytical and numerical Gubser solutions of the secondorder hydrodynamics, Phys. Rev. D91 (2015) 074027 [arXiv:1411.7767].
 [30] E. Molnar, H. Niemi, and D. H. Rischke, Numerical tests of causal relativistic dissipative fluid dynamics, Eur. Phys. J. C65 (2010), 615 [arXiv:0907.2583].
 [31] H. Niemi, G. S. Denicol, P. Huovinen, E. Molnar, and D. H. Rischke, Influence of a temperaturedependent shear viscosity on the azimuthal asymmetries of transverse momentum spectra in ultrarelativistic heavyion collisions, Phys. Rev. C86 (2012) 014909 [arXiv:1203.2452].
 [32] P. Huovinen and P. Petreczky, QCD Equation of State and Hadron Resonance Gas, Nucl. Phys. A837 (2010) 26 [arXiv:0912.2541].
 [33] H. B. Meyer, Glueball regge trajectories, [arXiv:heplat/0508002].
 [34] I. Karpenko, P. Huovinen, and M. Bleicher, 3+1 dimensional viscous hydrodynamic code for relativistic heavy ion collisions, Comput. Phys. Commun. 185 (2014) 3016 [arXiv:1312.4160].
 [35] W. Broniowski, M. Rybczynski, and P. Bozek, GLISSANDO: Glauber initialstate simulation and more.., Comput. Phys. Commun. 180 (2009) 69 [arXiv:0710.5731].
 [36] I. A. Karpenko, Y. M. Sinyukov, and K. Werner, Bulk observables in the hydrokinetic model of AA collisions at RHIC and LHC, Phys. Rev. C 87 (2013) 024914 [arXiv:1204.5351].
 [37] B. Abelev et al. [ALICE Collaboration], Centrality dependence of pi, K, p production in PbPb collisions at c.m. energy 2.76 TeV, Phys. Rev. C 88 (2013) 044910 [arXiv:1303.0737].

[38]
L. M. Satarov, I. N. Mishustin, A. V. Merdeev, and H. Stöcker,
1+1 dimensional hydrodynamics for highenergy heavyion collisions, Phys. Atom. Nucl. 70
(2007) 1773 [hepph/0611099];
Longitudinal fluiddynamics for ultrarelativistic heavyion collisions, Phys. Rev. C 75 (2007) 024903 [hepph/0606074].  [39] J. F. Paquet et al., Production of photons in relativistic heavyion collisions, Phys. Rev. C 93 (2016) 044906 [arXiv:1509.06738].
 [40] O. Linnyk, E. L. Bratkovskaya, W. Cassing, Effective QCD and transport description of dilepton and photon production in heavyion collisions and elementary processes, Prog. Part. Nucl. Phys. 87 (2016) 50 [arXiv:1512.08126].
 [41] D. Lohner [ALICE Collaboration], Measurement of directphoton elliptic flow in PbPb collisions at c.m. energy 2.76 TeV, J. Phys. Conf. Ser. 446 (2013) 012028 [arXiv:1212.3995].
 [42] J. Adam et al. [ALICE Collaboration], Direct photon production in PbPb collisions at c.m. energy 2.76 TeV, Phys. Lett. B754 (2016) 235 [arXiv:1509.07324]

[43]
R. Chatterjee, H. Holopainen, T. Renk, and K. J. Eskola,
Collision centrality and tau0 dependence of the emission of thermal photons from fluctuating initial state in ideal hydrodynamic calculation,
Phys. Rev. C 85 (2012) 064910 [arXiv:1204.2249].  [44] H. van Hees, M. He, and R. Rapp, Pseudocritical enhancement of thermal photons in relativistic heavyion collisions? Nucl. Phys. A933 (2015) 256 [arXiv:1404.2846].
 [45] O. Linnyk, V. Konchakovski, T. Steinert, W. Cassing, and E. L. Bratkovskaya, Hadronic and partonic sources of direct photons in relativistic heavyion collisions, Phys. Rev. C 92 (2015) 054914 [arXiv:1504.05699].
 [46] M. F. M. Lutz et al. [PANDA Collaboration], Physics Performance Report for PANDA: Strong Interaction Studies with Antiprotons, arXiv:0903.3905 [hepex].