# Fluid dynamics with saturated minijet initial conditions

in ultrarelativistic heavy-ion collisions

###### Abstract

Using next-to-leading order perturbative QCD and a conjecture of saturation to suppress the production of low-energy partons, we calculate the initial energy densities and formation times for the dissipative fluid dynamical evolution of the quark-gluon plasma produced in ultrarelativistic heavy-ion collisions. We identify the framework uncertainties and demonstrate the predictive power of the approach by a good global agreement with the measured centrality dependence of charged particle multiplicities, transverse momentum spectra and elliptic flow simultaneously for the Pb+Pb collisions at the LHC and Au+Au at RHIC. In particular, the shear viscosity in the different phases of QCD matter is constrained in this new framework simultaneously by all these data.

###### pacs:

25.75.-q, 25.75.Nq, 25.75.Ld, 12.38.Mh, 12.38.Bx, 24.10.Nz, 24.85.+pThe main goal of ultrarelativistic heavy-ion collisions at the Large Hadron Collider (LHC) and the Relativistic Heavy-Ion Collider (RHIC) is to determine the thermodynamic and kinetic properties of strongly interacting matter. The measured hadronic transverse momentum () spectra at the LHC and RHIC provide convincing evidence for a formation of a strongly collective system and a nearly thermalized quark-gluon plasma (QGP) Heinz:2013th (). In particular, the observed systematics of the Fourier harmonics of the azimuth-angle distributions, are remarkably consistent with a low-viscosity QCD matter whose expansion and cooling are describable with dissipative relativistic fluid dynamics Romatschke:2007mq (); Luzum:2008cw (); Schenke:2010rr (); Gale:2012rq (); Song:2010mg (); Song:2011qa (); Shen:2010uy (); Bozek:2009dw (); Bozek:2012qs (); Niemi:2011ix (); Niemi:2012ry ().

The essential inputs to the fluid dynamics are the initial energy density and flow of the matter created in the collision. However, the final state observables like multiplicities, spectra and , are also strongly affected through the fluid dynamical expansion by the viscosity and the equation of state (EoS). Thus the entire spacetime evolution, including partons in the colliding nuclei, the primary production and thermalization of QCD matter and the subsequent fluid dynamical evolution, becomes highly convoluted. Description of all these dynamics in a coherent way, leading to quantitative predictions and a meaningful determination of the QCD matter properties from the measurements, provides an ultimate challenge in the field. As discussed in this paper, the determination of e.g. the temperature dependence of the shear viscosity-to-entropy ratio calls for a simultaneous theory analysis of all possible bulk (low-) observables at the LHC and RHIC.

Parton saturation is a viable mechanism to control the otherwise unsuppressed production of soft small- quanta in hadronic and nuclear collisions Gribov:1984tu (); Mueller:1985wy (); McLerran:1993ni (); Eskola:1996ce (). In essence saturation means that there exists a semihard scale controlling the particle production in the collision. In the perturbative QCD (pQCD) + saturation framework we consider here, the primary particle production in + collisions is computed in collinear factorization by approaching the saturation at semi-hard scales from the perturbatively controllable high- side Eskola:1999fc (); Paatelainen:2012at (). Perturbative QCD provides an excellent description of hard processes in hadronic and nuclear collisions at interaction scales GeV Eskola:2009uj (). Moreover, this framework allows for a quantification of the particle production uncertainties, and their propagation through the fluid dynamical evolution in nuclear collisions Paatelainen:2012at (). In addition to the internal consistency of the pQCD-based approach, it should be noted that perturbative primary gluon production in heavy-ion collisions is complementary to the Color-Glass Condensate models Gelis:2010nm () which build on soft gluon fields. If these different high-energy QCD approaches produce similarly successful heavy-ion phenomenology, the overall uncertainty in determining the QCD matter properties can be dramatically reduced.

The present work has roots in the so-called EKRT saturation model Eskola:1999fc (), which successfully predicted the multiplicities and spectra in central + collisions at RHIC and LHC Eskola:2001bf (); Eskola:2002wx (); Eskola:2005ue (); Renk:2011gj (), and also the centrality dependence at RHIC Eskola:2000xq () (cf. Fig. 23(a) in Abelev:2008ab ()). Here we use the next-to-leading-order (NLO)-improved pQCD + saturation framework of Paatelainen:2012at () to calculate the initial QGP energy density profiles and formation times, and combine these with viscous fluid dynamics. We analyse the centrality dependence of charged particle multiplicities, spectra and elliptic flow () at the LHC and RHIC in terms of the few physical key-parameters of the framework. We show that a good simultaneous description of all these observables can indeed be obtained without retuning the framework from one collision system (cms-energy, nuclei, centrality) to another. This results in the robust predictive power of the approach, originating from the pQCD calculation of the QGP initial conditions. Most importantly, this predictive power enables us to study and restrict the ratio in the different QCD-matter phases more consistently in a simultaneous multiobservable analysis of the LHC and RHIC data.

Let us then discuss the details of our framework Paatelainen:2012at (). The rigorously calculable part is the minijet production in an + collision, in a rapidity interval and above a scale ,

(1) |

where is the transverse location, the impact parameter, and the standard nuclear thickness function with the Woods-Saxon nuclear density profile. The first -moment of the minijet distribution, , Eskola:1988yh (); Paatelainen:2012at () is in NLO

(2) |

where are the collinearly factorized minijet production cross sections and denote the phase-space differentials for the and cases Paatelainen:2012at (); tuominen:2000 (). We apply the CTEQ6M parton distribution functions (PDFs) Pumplin:2002vw () with the EPS09s impact-parameter dependent nuclear PDFs helenius:2012wd (). The measurement functions and define the hard scattering in terms of the minijet transverse momenta and the cut-off scale , as well as the total minijet produced in :

(3) |

where and is the step function. These functions, analogous to the jet definitions Kunszt:1992tn (), are constructed so that is a well-defined, infrared- and collinear-safe, quantity to compute. The hardness-parameter defines the minimum in the interval . As discussed in Paatelainen:2012at (), any is acceptable for the rigorous NLO computation.

Following the new angle in formulating the minijet saturation Paatelainen:2012at (), the production is expected to cease when the and higher-order partonic processes start to dominate over the conventional processes. For a central collision of identical nuclei of radii this leads to a transversally averaged saturation criterion , with an unknown, -independent, proportionality constant . Generalizing to non-zero impact parameters and localizing in the transverse coordinate plane gives

(4) |

where the l.h.s. is the -dependent NLO pQCD calculation defined in Eq. (1).

For given and , we solve the above equation for , and obtain the total in a mid-rapidity unit at saturation from the r.h.s. as . Once the solution is known, the local energy density is obtained Eskola:2001bf (); Eskola:2005ue () as

(5) |

where the local formation time is .

Fig. 1 shows examples of as a function of , calculated for fixed values of and with and three other fixed impact parameters corresponding to the centrality classes 0-5%, 20-30% and 40-50% in TeV Pb+Pb collisions at the LHC and 200 GeV Au+Au at RHIC. To a very good approximation, the b and s dependence of comes only through . This is due to the weak dependence of the nPDFs near the centres of the nuclei helenius:2012wd (). The approximate power-law scaling behaviour seen at large can then be understood as expained in Eskola:2001rx ().

We identify two main uncertainties in mapping the pQCD + saturation calculation to an initial state for fluid dynamics: (i) The energy density given by Eq. (5) is at a time , i.e. different at each transverse point s, while for fluid dynamics we need the initial condition at a fixed time . (ii) We cannot trust the pQCD calculation down to , but we need to set a minimum scale . Wherever we can use the pQCD calculation, but the other regions, i.e. low density edges, need to be treated separately.

We fix a minimum saturation scale as GeV. Correspondingly, the maximum formation time in our framework is . Then, we evaluate the energy densities from to using either the Bjorken free streaming (FS) or the Bjorken hydrodynamic scaling solution (BJ). We take these two limits to represent the uncertainty in the early pre-thermalization evolution: In the free streaming case the transverse energy is preserved, while the other limit corresponds to the case where a maximum amount of the transverse energy is reduced by the longitudinal pressure.

To obtain the energy density in the transverse region where , we use an interpolation , where the power with the total inelastic nucleon–nucleon cross-section , and fm. This smoothly connects the FS/BJ-evolved pQCD energy density to the binary profile at the dilute edge.

For the fluid-dynamical evolution, we use the state-of-the art 2+1 D setup previously employed in Ref. Niemi:2011ix (); Niemi:2012ry (); Niemi:2012aj (), assuming longitudinal boost invariance, a zero net-baryon density and thermalization at . The equations of motion are given by the conservation laws for energy and momentum, . The evolution equation of the shear-stress tensor is given by transient relativistic fluid dynamics IS (); Denicol:2012cn (); Molnar:2013lta (),

where the co-moving time derivative is denoted by the dot, is the shear viscosity coefficient, is the shear tensor, is the expansion rate, and the angular brackets denote the symmetrized and traceless projection, orthogonal to the fluid four-velocity . The coefficients of the non-linear terms are taken to be , and , where is the thermodynamic pressure and . For details of the numerical algorithm, see Refs. Niemi:2012ry (); Molnar:2009tx ().

The hadron spectra are calculated with the Cooper-Frye freeze-out procedure Cooper:1974mv () by using Israel’s and Stewart’s 14-moment ansatz for the dissipative correction to the local equilibrium distribution function, , where , with the index indicating different hadron species and the 4-momentum of the corresponding hadron. The freeze-out temperature is here always MeV. After calculating the thermal spectra, we include the contribution from all 2- and 3-particle decays of unstable resonances in the EoS.

We use the lattice QCD and hadron resonance gas (HRG) based EoS -PCE-v1 Huovinen:2009yb () with a chemical freeze-out temperature MeV. Although the rather high leads to an overabundance of protons, it however reproduces the low- region of the -spectra much better than e.g. MeV.

For a rough but realistic (non-constant Csernai:2006zz ()) shear viscosity description, we assume the ratio to decrease linearly as a function of temperature in the hadronic phase, be in a minimum at the matching-temperature 180 MeV of the HRG/QGP phases in the used EoS, and either to increase or stay constant vs. in the QGP phase Niemi:2011ix (); Niemi:2012ry (). Fig 2 shows the which in our framework best reproduce the coefficients simultaneously at RHIC and LHC.

At this point, we have a fixed framework with four correlated unknowns, , to be determined using the LHC and RHIC data on the centrality dependence of the charged particle multiplicities, spectra and . We proceed by scanning the parameters , and . In particular, we vary the minimum value and slopes of , keeping its general shape as in Fig. 2. Both the BJ and FS prethermal evolutions are considered. In practice, for each fixed , the remaining parameter is always tuned such that the multiplicity in the % most central collisions at the LHC is reproduced.

In Fig. 3a we show the computed centrality dependence of the charged hadron multiplicity in Pb+Pb collisions at TeV compared with the ALICE data Aamodt:2010pb (). As demonstrated here, several sets give a good agreement with the measurement. However, the data clearly favours and slightly the FS scenario over the BJ. For comparison, we also show the results obtained with the usual (non-saturation) eBC and eWN Glauber model initial states Kolb:2001qz ().

In Fig. 3b we show the multiplicities for Au+Au collisions at GeV, using the same parameter sets as in panel 3a, and compare with the PHENIX Adler:2004zn () and STAR Abelev:2008ab () data. We note that although the RHIC data would seem to favour a slightly smaller and the BJ case, the overall simultaneous agreement at RHIC and LHC is rather good.

As long as the centrality dependence of the multiplicity is described, all the scenarios studied here give a very good description of the charged hadron -spectra. More relevant parameters in this case are and which here are kept unchanged from RHIC to LHC. The obtained spectra are shown in Fig. 3c for the LHC and in Fig. 3d for RHIC. The data are from Refs. Abelev:2012hxa () and Ref. Adams:2003kv (); Adler:2003au (), correspondingly.

In Figs. 3e and 3f we show the elliptic flow coefficients at the LHC and RHIC, respectively. The data are from ALICE Aamodt:2010pa () and STAR Bai (). The coefficients depend strongly on the parametrization, and e.g. an ideal fluid description (not shown) does not give a correct . By scanning the as explained above, while keeping of order 1, we observed that a good simultaneous agreement with the measurements is obtained with the cases shown in Fig. 2. We emphasize that at RHIC, where the flow gradients are larger, one should require the agreement in particular in the small- region, where the dissipative corrections to the particle distributions do not grow unphysically large. Note especially that since is considered as a material property, it must not be changed between different collision systems.

To conclude, we computed the energy density profiles and formation times of the produced QGP at the LHC and RHIC in a new NLO-improved pQCD + local saturation framework of considerable predictive power. The subsequent evolution of these initial conditions was described with dissipative fluid dynamics. Identifying the framework uncertainties, a good global agreement with the measured centrality dependence of the low- bulk observables was obtained simultaneously at the LHC and RHIC. In particular, we were able to constrain the parametrization simultaneously by all these data. In the future, we will extend this analysis to include event-by-event fluctuations.

Acknowledgements. This work was financially supported by the Wihuri foundation (RP) and the Academy of Finland, projects 133005 (KJE) and 267842 (KT). We thank I. Helenius, T. Lappi, H. Mäntysaari and H. Paukkunen for useful discussions, and CSC-IT Center Science for supercomputing time.

## References

- (1) For a recent review and references, see U. W. Heinz and R. Snellings, Annu. Rev. Nucl. Part. Sci. 63, 123 (2013).
- (2) P. Romatschke and U. Romatschke, Phys. Rev. Lett. 99, 172301 (2007).
- (3) M. Luzum and P. Romatschke, Phys. Rev. C 78, 034915 (2008) [Erratum-ibid. C 79, 039903 (2009)].
- (4) B. Schenke, S. Jeon and C. Gale, Phys. Rev. Lett. 106, 042301 (2011); Phys. Rev. C 85, 024901 (2012); Phys. Lett. B 702, 59 (2011).
- (5) C. Gale et al., Phys. Rev. Lett. 110, 012302 (2013).
- (6) H. Song et al., Phys. Rev. Lett. 106, 192301 (2011) [Erratum-ibid. 109, 139904 (2012)]; Phys. Rev. C 83, 054910 (2011) [Erratum-ibid. C 86, 059903 (2012)].
- (7) H. Song, S. A. Bass and U. Heinz, Phys. Rev. C 83, 054912 (2011) [Erratum-ibid. C 87, 019902 (2013)].
- (8) C. Shen, U. Heinz, P. Huovinen and H. Song, Phys. Rev. C 82, 054904 (2010); Phys. Rev. C 84, 044903 (2011).
- (9) P. Bozek, Phys. Rev. C 81, 034909 (2010); Phys. Rev. C 85, 034901 (2012).
- (10) P. Bozek and I. Wyskiel-Piekarska, Phys. Rev. C 85, 064915 (2012).
- (11) H. Niemi et al., Phys. Rev. C 86, 014909 (2012).
- (12) H. Niemi et al., Phys. Rev. Lett. 106, 212302 (2011).
- (13) L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rept. 100 (1983) 1.
- (14) A. H. Mueller and J. Qiu, Nucl. Phys. B 268, 427 (1986).
- (15) L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994).
- (16) K. J. Eskola and K. Kajantie, Z. Phys. C 75, 515 (1997).
- (17) K. J. Eskola, K. Kajantie, P. V. Ruuskanen and K. Tuominen, Nucl. Phys. B 570, 379 (2000).
- (18) R. Paatelainen, K. J. Eskola, H. Holopainen and K. Tuominen, Phys. Rev. C 87, 044904 (2013).
- (19) K. J. Eskola, H. Paukkunen and C. A. Salgado, JHEP 0904, 065 (2009).
- (20) F. Gelis, E. Iancu, J. Jalilian-Marian and R. Venugopalan, Ann. Rev. Nucl. Part. Sci. 60, 463 (2010).
- (21) K. J. Eskola, P. V. Ruuskanen, S. S. Räsänen and K. Tuominen, Nucl. Phys. A 696, 715 (2001).
- (22) K. J. Eskola, H. Niemi, P. V. Ruuskanen and S. S. Räsänen, Phys. Lett. B 566, 187 (2003).
- (23) K. J. Eskola et al., Phys. Rev. C 72, 044904 (2005).
- (24) T. Renk, H. Holopainen, R. Paatelainen and K. J. Eskola, Phys. Rev. C 84, 014906 (2011).
- (25) K. J. Eskola, K. Kajantie and K. Tuominen, Phys. Lett. B 497, 39 (2001).
- (26) B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 79, 034909 (2009).
- (27) K. J. Eskola, K. Kajantie and J. Lindfors, Nucl. Phys. B 323, 37 (1989).
- (28) K. J. Eskola and K. Tuominen, Phys. Let. B 489, 329 (2000); Phys. Rev. D63, 114006 (2001).
- (29) J. Pumplin et al., JHEP 0207, 012 (2002).
- (30) I. Helenius, K. J. Eskola, H. Honkanen and C. A. Salgado, JHEP 1207, 073 (2012).
- (31) Z. Kunszt and D. E. Soper, Phys. Rev. D 46, 192 (1992).
- (32) K. J. Eskola, K. Kajantie and K. Tuominen, Nucl. Phys. A 700, 509 (2002).
- (33) H. Niemi, G. S. Denicol, H. Holopainen and P. Huovinen, Phys. Rev. C 87, 054901 (2013).
- (34) W. Israel and J. M. Stewart, Proc. R. Soc. A 365, 43 (1979); Ann. Phys. (N.Y.) 118, 341 (1979).
- (35) G. S. Denicol, H. Niemi, E. Molnar and D. H. Rischke, Phys. Rev. D 85, 114047 (2012).
- (36) E. Molnar, H. Niemi, G. S. Denicol and D. H. Rischke, arXiv:1308.0785 [nucl-th].
- (37) E. Molnar, H. Niemi and D. H. Rischke, Eur. Phys. J. C 65, 615 (2010).
- (38) F. Cooper, G. Frye, Phys. Rev. D 10, 186 (1974).
- (39) P. Huovinen, P. Petreczky, Nucl. Phys. A837, 26 (2010).
- (40) L. P. Csernai, J. I. Kapusta, L. D. McLerran, Phys. Rev. Lett. 97, 152303 (2006).
- (41) K. Aamodt et al. [ALICE Collaboration], Phys. Rev. Lett. 105, 252301 (2010).
- (42) P. F. Kolb et al., Nucl. Phys. A 696, 197 (2001).
- (43) S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 71, 034908 (2005) [Erratum-ibid. C 71, 049901 (2005)].
- (44) B. Abelev et al. [ALICE Collaboration], Phys. Lett. B 720, 52 (2013).
- (45) J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 91, 172302 (2003).
- (46) S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 69, 034910 (2004).
- (47) K. Aamodt et al. [ALICE Collaboration], Phys. Rev. Lett. 105, 252302 (2010).
- (48) Y. Bai, Ph.D. Thesis, Nikhef and Utrecht University, The Netherlands (2007); A. Tang [STAR Collaboration], arXiv:0808.2144 [nucl-ex].