![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0024v1-figure-a97e6d3b-cfff-4adf-9b5b-a4c4f7e0ab30.jpg)

Nuclear Physics A

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0024v1-figure-7b48c552-5357-42b4-b0aa-ce6e93326099.jpg)

Nuclear Physics A 00 (2014) 1-4

## Extracting the bulk viscosity of the quark-gluon plasma

Jean-Bernard Rose a , Jean-Francois Paquet , Gabriel S. Denicol , Matthew Luzum a a a,b , Bjoern Schenke c , Sangyong Jeon a , Charles Gale a a Department of Physics, McGill University, 3600 rue University, Montr´ eal, Qu´ ebec H3A 2T8, Canada b Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA c Physics Department, Brookhaven National Lab, Building 510A, Upton, NY, 11973, USA

## Abstract

We investigate the implications of a nonzero bulk viscosity coe GLYPH&lt;14&gt; cient on the azimuthal momentum anisotropy of ultracentral relativistic heavy ion collisions at the Large Hadron Collider. We find that, with IP-Glasma initial conditions, a finite bulk viscosity coe GLYPH&lt;14&gt; cient leads to a better description of the flow harmonics in ultracentral collisions. We then extract optimal values of bulk and shear viscosity coe GLYPH&lt;14&gt; cients that provide the best agreement with flow harmonic coe GLYPH&lt;14&gt; cients data in this centrality class.

## 1. Introduction

Currently, most fluid-dynamical simulations of relativistic heavy ion collisions take into account only dissipative e GLYPH&lt;11&gt; ects originating from shear viscosity. However, QCD is a nonconformal field theory and in principle there is no a priori reason to neglect the e GLYPH&lt;11&gt; ects of bulk viscous pressure. In this contribution, we explore the phenomenological implications of a nonzero bulk viscosity coe GLYPH&lt;14&gt; cient on the azimuthal momentum anisotropy of ultracentral heavy ion collisions at the Large Hadron Collider (LHC). We then extract optimal values of bulk and shear viscosity coe GLYPH&lt;14&gt; cients that are able to describe the flow harmonic coe GLYPH&lt;14&gt; cients in this centrality class. We find that, with IP-Glasma initial conditions, a finite bulk viscosity coe GLYPH&lt;14&gt; cient improves the description of flow harmonics in ultracentral collisions.

## 2. Fluid-dynamical model of the collision

In this work we simulate ultrarelativistic heavy ion collisions in the 0-1% centrality class at LHC energies. The initial energy density profile is obtained using the IP-Glasma initial condition model, with a thermalization time of GLYPH&lt;28&gt; 0 = 0 4 fm [1]. : The time evolution of this system is determined using relativistic dissipative fluid dynamics, with the main equations of motion being the continuity equation for the energy-momentum tensor, T GLYPH&lt;22&gt;GLYPH&lt;23&gt; , @GLYPH&lt;22&gt; T GLYPH&lt;22&gt;GLYPH&lt;23&gt; = 0, where T GLYPH&lt;22&gt;GLYPH&lt;23&gt; = " u u GLYPH&lt;22&gt; GLYPH&lt;23&gt; GLYPH&lt;0&gt; GLYPH&lt;1&gt; GLYPH&lt;22&gt;GLYPH&lt;23&gt; ( P 0 +GLYPH&lt;5&gt; ) + GLYPH&lt;25&gt; GLYPH&lt;22&gt;GLYPH&lt;23&gt; , with " being the energy density, P 0 the thermodynamic pressure, GLYPH&lt;5&gt; the bulk viscous pressure, and GLYPH&lt;25&gt; GLYPH&lt;22&gt;GLYPH&lt;23&gt; the shear-stress tensor. We further introduced the projection operator GLYPH&lt;1&gt; GLYPH&lt;22&gt;GLYPH&lt;23&gt; = g GLYPH&lt;22&gt;GLYPH&lt;23&gt; GLYPH&lt;0&gt; u u GLYPH&lt;22&gt; GLYPH&lt;23&gt; onto the 3-space orthogonal to the fluid velocity. We assume that the baryon number density and di GLYPH&lt;11&gt; usion are zero at all space-time points and our metric convention is g GLYPH&lt;22&gt;GLYPH&lt;23&gt; = diag( + 1 ; GLYPH&lt;0&gt; 1 GLYPH&lt;0&gt; 1 GLYPH&lt;0&gt; 1).

The conservation laws and equation of state P 0 ( " ) are not enough to determine the time evolution of T GLYPH&lt;22&gt;GLYPH&lt;23&gt; . In order to do so, one still needs to provide the time-evolution or constitutive equations satisfied by GLYPH&lt;5&gt; and GLYPH&lt;25&gt; GLYPH&lt;22&gt;GLYPH&lt;23&gt; . For this purpose, we employ relaxation-time equations derived from kinetic theory [2, 3], solved numerically using the music

hydrodynamic simulation [4, 5]. Explicitly, we solve

$$\tau _ { \Pi } \dot { \Pi } + \Pi \ = \ - \zeta \theta - \delta _ { \Pi \Pi } \Pi \theta + \lambda _ { \Pi \pi } \pi ^ { \mu \nu } \sigma _ { \mu \nu }$$

$$\tau _ { \pi } \dot { \pi } ^ { ( \mu \nu ) } + \pi ^ { \mu \nu } \ = \ 2 \eta \sigma ^ { \mu \nu } + 2 \tau _ { \pi } \pi _ { \alpha } ^ { ( \mu } \omega ^ { \nu ) \alpha } - \delta _ { \pi \pi } \pi ^ { \mu \nu } \theta + \varphi _ { \gamma } \pi _ { \alpha } ^ { ( \mu } \pi ^ { \nu ) \alpha } - \tau _ { \pi \pi } \pi _ { \alpha } ^ { ( \mu } \sigma ^ { \nu ) \alpha } + \lambda _ { \pi \Pi } \Pi \sigma ^ { \mu \nu }.$$

The above equations of motion include nonlinear terms that couple bulk viscous pressure to the shear-stress tensor and vice-versa. Recently, they were shown to be in good agreement with solutions of the 0 + 1 Anderson-Witting equation in the massive limit [6]. The transport coe GLYPH&lt;14&gt; cients GLYPH&lt;28&gt; GLYPH&lt;5&gt; , GLYPH&lt;16&gt; , GLYPH&lt;14&gt; GLYPH&lt;5&gt;GLYPH&lt;5&gt; , GLYPH&lt;21&gt; GLYPH&lt;5&gt; GLYPH&lt;25&gt; , GLYPH&lt;28&gt; GLYPH&lt;25&gt; , GLYPH&lt;17&gt; , GLYPH&lt;14&gt; GLYPH&lt;25&gt;GLYPH&lt;25&gt; , ' 7 , GLYPH&lt;28&gt; GLYPH&lt;25&gt;GLYPH&lt;25&gt; , and GLYPH&lt;21&gt; GLYPH&lt;25&gt; GLYPH&lt;5&gt; are complicated functions of temperature. For strongly coupled fluids, such as the quark-gluon plasma created at the LHC, their values are not known. As a first step, we fix the parametric dependence of all transport coe GLYPH&lt;14&gt; cients on the shear viscosity, GLYPH&lt;17&gt; , and bulk viscosity, GLYPH&lt;16&gt; , using formulas derived from the Boltzmann equation near the conformal limit [3].

The shear viscosity is assumed to be proportional to the entropy density, i.e., GLYPH&lt;17&gt; = as , where the parameter a is fixed by fitting the data. In simulations using IP-Glasma initial conditions that include only dissipative e GLYPH&lt;11&gt; ects from shear viscosity, one usually obtains a = 0 20 : GLYPH&lt;0&gt; 0 22 at LHC energies [7, 8]. The bulk viscosity coe : GLYPH&lt;14&gt; cient is parametrized as [3]

$$\zeta = b \times \eta \left ( \frac { 1 } { 3 } - c _ { s } ^ { 2 } \right ) ^ { 2 },$$

where cs is the velocity of sound and b a parameter that will also be adjusted to fit the data. Of course, the actual temperature dependence of each viscosity coe GLYPH&lt;14&gt; cient is not known for QCD matter. Nevertheless, the aforementioned parametrizations should be able to capture the relevant features of dissipation and will be enough for the purposes of this initial study.

The final state hadrons are produced following the usual Cooper-Frye procedure [9], in the same way as outlined in Ref. [8]. The only di GLYPH&lt;11&gt; erence is that we have to specify the nonequilibrium correction to the hadron momentum distribution function due to bulk viscosity, in addition to the usual correction due to shear viscosity . Here, we employ the distribution derived from the Boltzmann equation using the relaxation time approximation, also used by Bozek et al [10],

$$\delta f _ { \Pi } ^ { i } \ = \ - f _ { 0 \text{k} } ^ { i } \left ( 1 \pm f _ { 0 \text{k} } ^ { i } \right ) \left [ \frac { 1 } { u _ { \mu } k _ { i } ^ { \mu } } \frac { m _ { i } ^ { 2 } } { 3 } - \left ( \frac { 1 } { 3 } - c _ { s } ^ { 2 } \right ) u _ { \mu } k _ { i } ^ { \mu } \right ] \frac { \Pi } { C _ { \Pi } },$$

$$C _ { \Pi } \ = \ \frac { 1 } { 3 } \sum _ { i = 1 } ^ { N } m _ { i } ^ { 2 } g _ { i } \int d K _ { i } \, f _ { 0 \text{k} } ^ { i } \left ( 1 \pm f _ { 0 \text{k} } ^ { i } \right ) \left [ \frac { 1 } { u _ { \mu } k _ { i } ^ { \mu } } \frac { m _ { i } ^ { 2 } } { 3 } - \left ( \frac { 1 } { 3 } - c _ { s } ^ { 2 } \right ) u _ { \mu } k _ { i } ^ { \mu } \right ].$$

Above, k GLYPH&lt;22&gt; i , mi , and gi are the 4-momentum, mass, and degeneracy factor of the i -th hadron species, respectively, N is the number of hadrons considered in the calculation, and f i 0 k is the Fermi-Dirac or Bose-Einstein distribution function.

We note that the IP-Glasma initial conditions used in this work include the e GLYPH&lt;11&gt; ect of nucleon-nucleon correlations in the initial state. Such correlations were recently found to be important when simulating central and ultracentral collisions with the IP-Glasma model [8].

## 3. Results

We calculate the integrated elliptic flow coe GLYPH&lt;14&gt; cient vn f 2 g via the cumulant method. In order to compare with CMSdata [12], we use a lower cut o GLYPH&lt;11&gt; of 0.3 GeV when performing integrals in the transverse momentum, pT . All our calculations were done using the parametrization of the lattice QCD equation of state by Huovinen and Petreczky [11]. The calculations with bulk viscosity were all performed in chemical equilibrium and with a freezeout temperature of TFO = 140 MeV. The calculations with only shear viscosity were performed using the same parameters as in [8], with a chemical freezeout temperature of TCF = 150 MeV and TFO = 103 MeV. In all figures, the shaded bands represent the statistical uncertainty.

In Fig. 1 we show the results for a given shear viscosity, GLYPH&lt;17&gt;= s = 0 15, and several values of : b = 15, 45, and 75. We see that increasing the value of bulk viscosity improves the agreement with the flow harmonics measured by CMS. In Fig. 1(b), we compare a simulation which includes only shear viscosity, with GLYPH&lt;17&gt;= s = 0 22, to another with : GLYPH&lt;17&gt;= s = 0 15 :

Figure 1. (Color online) Comparison of the flow harmonics, vn f 2 , calculated using IP-Glasma initial conditions, with with CMS data for the 0-1% g centrality class and for several values of the parameter b.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0024v1-figure-30f09d1c-a417-4b2d-8337-11e5e77e6db0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0024v1-figure-7dbb1f00-36e6-4d3b-838d-e9a941eaf797.jpg)

and b = 15. Both calculations provide an almost equally good agreement with the flow harmonics measured by CMS, but with rather di GLYPH&lt;11&gt; erent values of shear viscosity. In this case, including bulk viscosity reduced the value of GLYPH&lt;17&gt;= s by 32%, from 0.22 to 0.15.

In Fig. 2(a) we show what happens when the parameter b is fixed as b = 45 and the shear viscosity is varied, taking the values GLYPH&lt;17&gt;= s = 0 12, 0.15, and 0.18. : Note that the parametrization of bulk viscosity used in this work, Eq. (3), is linear in the shear viscosity, so keeping b constant does not imply that the bulk viscosity is kept fixed. In this figure we see that the value of shear viscosity that provides the best overall agreement with the CMS data is GLYPH&lt;17&gt;= s = 0 15, : which was the value already employed in the simulations discussed above.

Figure 2. In the left panel we show the flow harmonics, vn f 2 , calculated using IP-Glasma initial conditions, for g b = 45 and several values of GLYPH&lt;17&gt;= s (left panel). In the right panel we show the e GLYPH&lt;11&gt; ect of GLYPH&lt;14&gt; i f GLYPH&lt;5&gt; on vn f 2 . g

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0024v1-figure-5793be14-4961-454f-a9d2-f1040cc00020.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0024v1-figure-3ff2857a-95ba-4169-8cc6-ddb8f8ef8f50.jpg)

In Fig. 2(b) the results for GLYPH&lt;17&gt;= s = 0 18 and : b = 45 are shown calculated with and without the bulk component of GLYPH&lt;14&gt; f , GLYPH&lt;14&gt; i f GLYPH&lt;5&gt; . In contrast to the shear correction, the bulk correction to the single particle distribution function used in this work lead to an increase in the flow harmonics. Such an increase was systematically stronger for higher harmonics: v 2 f 2 g is increased by 4%, v 3 f 2 g by 8%, v 4 f 2 g by 20%, v 5 f 2 g by 57%, and v 6 f 2 g by 117%. Therefore, in order to properly

describe higher harmonics, it is important to take into account the GLYPH&lt;14&gt; i f GLYPH&lt;5&gt; contribution. However, one should note that di GLYPH&lt;11&gt; erent momentum dependences for this quantity can provide rather di GLYPH&lt;11&gt; erent results.

Finally, we show the e GLYPH&lt;11&gt; ect of bulk viscosity on the mean transverse momentum of pions, h p GLYPH&lt;25&gt; T i , in Fig. 3. Our calculations are compared to ALICE data [13] in the 0-5% centrality class and, for this purpose, we used a lower cut-o GLYPH&lt;11&gt; of 0.12 GeV when performing integrals in transverse momentum. The e GLYPH&lt;11&gt; ect of bulk viscosity on h p GLYPH&lt;25&gt; T i is not small and can be used to further restrict the size of the bulk viscosity coe GLYPH&lt;14&gt; cient. Here, we see that the best description of h p GLYPH&lt;25&gt; T i happens when b = 20-40, for the case where GLYPH&lt;17&gt;= s = 0 15. Note that the : h p GLYPH&lt;25&gt; T i has a strong dependence on GLYPH&lt;14&gt; i f GLYPH&lt;5&gt; , so the optimal values we obtained can be modified if another GLYPH&lt;14&gt; i f GLYPH&lt;5&gt; is used.

Figure 3. (Color online) The mean transverse momentum of pions, h p GLYPH&lt;25&gt; T i , for several values of b compared to the ALICE data [13].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0024v1-figure-39be8dd3-cf8e-4a4a-8a6c-d486648d780d.jpg)

## 4. Conclusions

In this contribution we performed a preliminary study on the e GLYPH&lt;11&gt; ects of bulk viscous pressure in ultracentral relativistic heavy ion collisions at the LHC. We found that bulk viscous pressure can have non-negligible e GLYPH&lt;11&gt; ect on the flow harmonic coe GLYPH&lt;14&gt; cients in this centrality class. Furthermore, we showed that, at least for IP-Glasma initial conditions, the inclusion of bulk viscosity leads to a better description of the data. Although describing the anisotropy distribution of ultracentral events remains a challenge, we expect bulk viscosity to be a key ingredient in understanding such observables.

Acknowledgments: The authors thank H. Niemi and J. Jia for fruitful discussions. This work was supported in part by the Natural Sciences and Engineering Research Council of Canada, and by the U. S. DOE Contract No. DEAC02-98CH10886. G.S. Denicol acknowledges support through a Banting Fellowship of the Natural Sciences and Engineering Research Council of Canada, and C. G. acknowledges support from the Hessian Initiative for Excellence (LOEWE) through the Helmholtz International Center for FAIR (HIC for FAIR).

## References

- [1] B. Schenke, P. Tribedy and R. Venugopalan, Phys. Rev. C 86 , 034908 (2012).
- [2] G. S. Denicol, H. Niemi, E. Molnar and D. H. Rischke, Phys. Rev. D 85 , 114047 (2012).
- [3] G. S. Denicol, S. Jeon and C. Gale, accepted for publication in Phys. Rev. C, arXiv:1403.0962 [nucl-th].
- [4] B. Schenke, S. Jeon and C. Gale, Phys. Rev. C 82 , 014903 (2010).
- [5] H. Marrochio, J. Noronha, G. S. Denicol, M. Luzum, S. Jeon and C. Gale, arXiv:1307.6130 [nucl-th].
- [6] G. S. Denicol, W. Florkowski, R. Ryblewski and M. Strickland, arXiv:1407.4767 [hep-ph].
- [7] C. Gale, S. Jeon, B. Schenke, P. Tribedy and R. Venugopalan, Phys. Rev. Lett. 110 , 012302 (2013).
- [8] G. S. Denicol, C. Gale, S. Jeon, J. -F. Paquet and B. Schenke, arXiv:1406.7792 [nucl-th].
- [9] F. Cooper and G. Frye, Phys. Rev. D 10 , 186 (1974).
- [10] P. Bozek, Phys. Rev. C 81 , 034909 (2010).
- [11] P. Huovinen and P. Petreczky, Nucl. Phys. A 837 , 26 (2010).
- [12] S. Chatrchyan et al. [CMS Collaboration], JHEP 1402 , 088 (2014).
- [13] B. Abelev et al. [ALICE Collaboration], Phys. Rev. C 88 , no. 4, 044910 (2013).
J. -B. Rose J. -F. Paquet G. S. Denicol M. Luzum B. Schenke S. Jeon C. Gale
Extracting the bulk viscosity of the quark-gluon plasma
We investigate the implications of a nonzero bulk viscosity coefficient on
the azimuthal momentum anisotropy of ultracentral relativistic heavy ion
collisions at the Large Hadron Collider. We find that, with IP-Glasma initial
conditions, a finite bulk viscosity coefficient leads to a better description
of the flow harmonics in ultracentral collisions. We then extract optimal
values of bulk and shear viscosity coefficients that provide the best agreement
with flow harmonic coefficients data in this centrality class.