![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-8f9c7fad-d66c-4376-ae65-985a20e872fb.jpg)

Nuclear Physics A

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-0490a55a-f4a9-4c33-aba3-88af75f6335e.jpg)

Nuclear Physics A 00 (2014) 1-4

## Heavy Ion Collision evolution modeling with ECHO-QGP

V. Rolando a,b , G. Inghirami c,d , A. Beraudo , L. Del Zanna e c,d,f , F. Becattini c,d , V. Chandra , d A. De Pace e , M. Nardi e a Dipartimento di Fisica e Scienze della Terra, Universit` di Ferrara, Via Saragat 1, I-44100 Ferrara, Italy a b INFN - Sezione di Ferrara, Via Saragat 1, I-44100 Ferrara, Italy c Dipartimento di Fisica e Astronomia, Universit` di Firenze, Via G. Sansone 1, I-50019 Sesto F .no (Firenze), Italy a d INFN - Sezione di Firenze, Via G. Sansone 1, I-50019 Sesto F.no (Firenze), Italy e INFN - Sezione di Torino, Via P. Giuria 1, I-10125 Torino, Italy f INAF - Osservatorio Astrofisico di Arcetri, L.go E. Fermi 5, I-50125 Firenze, Italy

## Abstract

We present a numerical code modeling the evolution of the medium formed in relativistic heavy ion collisions, ECHO-QGP. The code solves relativistic hydrodynamics in (3 + 1) GLYPH&lt;0&gt; D, with dissipative terms included within the framework of Israel-Stewart theory; it can work both in Minkowskian and in Bjorken coordinates. Initial conditions are provided through an implementation of the Glauber model (both Optical and Monte Carlo), while freezeout and particle generation are based on the Cooper-Frye prescription. The code is validated against several test problems and shows remarkable stability and accuracy with the combination of a conservative (shock-capturing) approach and the high-order methods employed. In particular it beautifully agrees with the semi-analytic solution known as Gubser flow, both in the ideal and in the viscous Israel-Stewart case, up to very large times and without any ad hoc tuning of the algorithm.

Keywords: QGP, Hydrodynamics, Heavy-ion, Heavy-ion collision

## 1. Introduction

During the past years hydrodynamics has been aknowledeged as the most powerful tool to study the evolution of the medium formed in high-energy nuclear collisions. For this reason a variety of calculations and numerical codes have been developed, starting from the simple longitudinal boost-invariant Bjorken (1 + 1) GLYPH&lt;0&gt; D calculation [1], passing through (2 + 1) GLYPH&lt;0&gt; D and (3 + 1) GLYPH&lt;0&gt; D codes solving ideal hydrodynamics [2-4], up to more complex and structured codes simulating a full (3 + 1) GLYPH&lt;0&gt; D viscous evolution for the QGP fluid[5-10].

In order to solve dissipative hydrodynamics, one needs to deal with a set of non-linear partial di GLYPH&lt;11&gt; erential equations, which includes the conservation of the energy momentum tensor ( @GLYPH&lt;22&gt; T GLYPH&lt;22&gt;GLYPH&lt;23&gt; = 0) and of conserved charges like the baryon number ( @GLYPH&lt;22&gt; N GLYPH&lt;22&gt; = 0), together with the time evolution of the bulk viscosity ( GLYPH&lt;5&gt; ) and the shear stress tensor ( GLYPH&lt;25&gt; GLYPH&lt;22&gt;GLYPH&lt;23&gt; ). The system is closed by the choice of a suitable Equation of State (EoS). While the naive relativistic extension of the Navier-Stokes equations is a GLYPH&lt;11&gt; ected by a well known causality problem, a consistent theoretical setup to include dissipative e GLYPH&lt;11&gt; ects is the one formulated by Israel and Stewart [11] .

In this context ECHO-QGP was presented last year [12], equipped with second-order treatment of causal relativistic viscosity e GLYPH&lt;11&gt; ects and unique features such as the possibility of choosing between two di GLYPH&lt;11&gt; erent metric tensors and of solving either purely ideal or viscous hydrodynamics equations. As the other analogous tools, ECHO-QGP is characterized by a modular structure, allowing one to address the modeling of the initial conditions, the solution of the hydrodynamic evolution of the medium - representing the core of the code - and, eventually, the simulation of

the final particle decoupling. ECHO-QGP is highly customizable, allowing the user to choose among a variety of initial conditions (including both optical and Monte Carlo Glauber model) as well as the possibility of defining an energy density (or entropy density) profile from scratch. The same principle applies to the equation of state: ECHOQGP can handle both analytical and tabulated equations of state. The decoupling stage exploits the Cooper-Frye prescription[13], allowing the calculation of a mean spectrum of or the Monte Carlo generation of a discrete set of particles. The hypersurface detection has been recently upgraded with the embedding of a refined algorithm which creates a smooth mesh[14].

In order to ensure the reliability of the solving algorithm for the evolution stage, in our original publication [12] we showed how ECHO-QGP is able to overcome several numerical tests, displaying a beautiful agreement with some special analytic solutions. In the present contribution the above e GLYPH&lt;11&gt; orts in validationg the code will be briefly reminded and extended to the case of the so called Gubser-flow (both ideal and viscous), meanwhile appeared in the literature [15-17]. The latter represents a very important test, since the code has to reproduce a very non-trivial flow in (2 + 1)D in the presence of non-vanishing viscosity and relaxation time. ECHO-QGP turns out to be able to overcome such a test.

Finally, having at our disposal such a validated tool, we will outline our future programs.

## 2. Testing ECHO-QGP

ECHO-QGP has been widely tested, displaying agreement with analytic and semi-analytic solutions and other publicly available numerical codes.

Among the various test performed, an important role is played by the (2 + 1)-D shock-tube problem. Shock-capturing numerical schemes are designed to handle and evolve discontinuous quantities invariably arising due to the nonlinear nature of the fluid equations. In order to validate these codes, typical tests are the so-called shock-tube problems, performed in the heavy-ion field also by [18]. The test gives a good esteem of the behavior of the fluid in presence of non-vanishing GLYPH&lt;17&gt;= s . In figure 1, we compare the energy density profile for the viscid and unviscid case, where the high accuracy of the results and the absence of numerical spurious oscillations near the shock front in the ideal case can be appreciated.

(a) Shock tube test: velocity profile

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-ed3e24b7-32eb-405e-acc5-216de007d033.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-345693b2-d9d3-4aa4-8de4-7055745cb592.jpg)

(b) Shock tube test: energy density profile

Figure 1: Shock tube test for ECHO-QGP. The test is performed placing an initial diaphragm along the diagonal of a square box (201 points and size 10 fm along x and y) adopting Minkowskian Cartesian coordinates, and letting ECHO-QGP evolve from t = 1 up to t = 4 fm c. = The temperatures in the sides of the diaphragm are chosen to be T L = 0 4 GeV and : T R = 0 2 GeV. It can be appreciated how the viscosity a : GLYPH&lt;11&gt; ects the evolution, and how ECHO-QGP handles shocks in the ideal case without the introduction of dumps.

Another fundamental qualifying test which numerical hydrodynamic codes have recently started addressing is the so called Gubser-flow. The latter is a (2 + 1)D solution for a conformal fluid (with an EoS P = e = 3) characterized by

longitudinal boost invariance and non-trivial azymuthally symmetric radial expansion. It was derived first by Gubser et al. [15, 16] for an ideal fluid and extended only last year to the viscous case, within the Israel-Stewart setup, by Marrochio et al. [17]. The derivation of the solution is based on a rescaling of the metric (Weyl rescaling) and on a change of coordinates, which allows one to exploit at best the symmetries of the system arising from conformal invariance. In particular, one applies to the metric ds 2 = GLYPH&lt;0&gt; d GLYPH&lt;28&gt; 2 + dr 2 + r 2 d GLYPH&lt;30&gt; 2 + GLYPH&lt;28&gt; 2 d GLYPH&lt;17&gt; 2 the rescaling ds 2 ! ds ˆ 2 GLYPH&lt;17&gt; ds 2 =GLYPH&lt;28&gt; 2 .

In addition, one has to perform the coordinate transformation ds ˆ 2 = GLYPH&lt;0&gt; d GLYPH&lt;26&gt; 2 + cosh 2 GLYPH&lt;26&gt; ( d GLYPH&lt;18&gt; 2 + sin 2 GLYPH&lt;18&gt; d GLYPH&lt;30&gt; 2 ) + d GLYPH&lt;17&gt; 2 ; where the new coordinates are given by sinh GLYPH&lt;26&gt; GLYPH&lt;17&gt; GLYPH&lt;0&gt; 1 GLYPH&lt;0&gt; q 2 ( GLYPH&lt;28&gt; 2 GLYPH&lt;0&gt; r 2 ) 2 q GLYPH&lt;28&gt; and tan GLYPH&lt;18&gt; GLYPH&lt;17&gt; 2 qr 1 + q 2 ( GLYPH&lt;28&gt; 2 GLYPH&lt;0&gt; r 2 ) , where q is an arbitrary energy scale. In the new space (with quantities labeled by a hat) the fluid is at rest, ˆ u GLYPH&lt;22&gt; = ( GLYPH&lt;0&gt; 1 0 0 0), and the only equations to solve are: ; ; ;

$$\frac { 1 } { \hat { T } } \frac { d \hat { T } } { d \rho } + \frac { 2 } { 3 } \tanh \rho = \frac { 1 } { 3 } \tanh \rho \bar { \pi } _ { \eta } ^ { \eta }, & \hat { \tau } _ { R } \left [ \frac { d \bar { \pi } _ { \eta } ^ { \eta } } { d \rho } + \frac { 4 } { 3 } \left ( \bar { \pi } _ { \eta } ^ { \eta } \right ) ^ { 2 } \tanh \rho \right ] + \bar { \pi } _ { \eta } ^ { \eta } = \frac { 4 } { 3 } \frac { \hat { \eta } } { \hat { S } \hat { T } } \tanh \rho, & ( 1 )$$

where ¯ GLYPH&lt;25&gt; GLYPH&lt;17&gt; ˆ GLYPH&lt;25&gt;= ˆ ˆ. Ts Physical quantities in Minkowski space can be obtained from the above equations through the mapping

$$T = \hat { T } / \tau, \quad u _ { \mu } = \tau \frac { \partial \hat { x } ^ { \alpha } } { \partial x ^ { \mu } } \hat { u } _ { \alpha }, \quad \pi _ { \mu \nu } = \frac { 1 } { \tau ^ { 2 } } \frac { \partial \hat { x } ^ { \alpha } } { \partial x ^ { \mu } } \frac { \partial \hat { x } ^ { \beta } } { \partial x ^ { \nu } } \hat { \pi } _ { \alpha \beta }$$

and compared to the numerical solution provided by ECHO-QGP, as displayed in Fig. 2.

The results have been obtained without any fine tuning of the parameters or modifications of the reconstruction algorithm. Nonetheless, ECHO-QGP reproduces with very high accuracy (the discrepancy is on average of the order of 0.1%) and up to late times (we verified the agreement up to 10 fm c) the temperature (see Fig. / 2d), all the shearstress tensor components (see Figs. 2a,2b,2c) and obviously the flow (not shown).

## 3. Conclusions

ECHO-QGP guarantees stability and precision, as well as completeness in documentation and ease in the usage. We presented all those ECHO-QGP features in [12] and we remark its suitablity to model heavy-ion collisions showing how well it reproduces the flow and shear stress tensor in the viscous relativistic frame proposed in ref. [15-17]. Moreover, ECHO-QGP has been positively used for the investigation fluctuation propagations[19, 20], and it is currently used to study vorticity e GLYPH&lt;11&gt; ects on the directed flow (in preparation).

## 4. Aknowledgements

We would like to thank G. Denicol and the authors of the reference [17] for interesting discussions and advices about the semi-analytic solution. V.R. would like to thank prof. R. Tripiccione, G. Pagliara and A. Drago for their help in many interesting discussions about the decoupling stage.

This work has been supported by the Italian Ministry of Education and Research, grant n. 2009WA4R8W.

## References

- [1] J. Bjorken, Phys.Rev. D27 (1983) 140-151. doi:10.1103/PhysRevD.27.140 .
- [2] P. F. Kolb, J. Sollfrank, U. W. Heinz, Phys.Rev. C62 (2000) 054909. arXiv:hep-ph/0006129 doi:10.1103/PhysRevC.62.054909 , .
- [3] P. F. Kolb, R. Rapp, Phys.Rev. C67 (2003) 044903. arXiv:hep-ph/0210222 doi:10.1103/PhysRevC.67.044903 , .
- [4] P. F. Kolb, U. W. Heinz, arXiv:nucl-th/0305084 .

.

- [5] M. Luzum, P. Romatschke, Physical Review Letters 103 (26) (2009) 1-4. doi:10.1103/PhysRevLett.103.262302
- [6] R. Baier, P. Romatschke, U. A. Wiedemann, Physical Review C 73 (6) (2006) 18. arXiv:0602249 doi:10.1103/PhysRevC.73.064903 , .
- [7] P. Romatschke, Int.J.Mod.Phys. E19 (2010) 1-53. arXiv:0902.3663 doi:10.1142/S0218301310014613 , .
- [8] S. Ryu, S. Jeon, C. Gale, B. Schenke, C. Young, Nucl.Phys. A904-905 (2013) 389c-392c. doi:10.1016/j.nuclphysa.2013.02.031 .
- [9] C. Gale, S. Jeon, B. Schenke, P. Tribedy, R. Venugopalan, Phys.Rev.Lett. 110 (2013) 012302. doi:10.1103/PhysRevLett.110.012302 .
- [10] B. Schenke, S. Jeon, C. Gale, Phys.Rev. C82 (2010) 014903. arXiv:1004.1408 doi:10.1103/PhysRevC.82.014903 , .
- [11] W. Israel, J. Stewart, Annals Phys. 118 (1979) 341-372. doi:10.1016/0003-4916(79)90130-1 .
- [12] L. Del Zanna, V. Chandra, G. Inghirami, V. Rolando, A. Beraudo, A. De Pace, G. Pagliara, G. Drago, F. Becattini, EPJ C73 (2013) 2524. doi:10.1140/epjc/s10052-013-2524-5 .
- [13] Cooper, Fred and Frye, Graham, Physical Review D 10 (1) (1974) 186-189. doi:{10.1103/PhysRevD.10.186} .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-59ddaf44-7644-498a-9c1c-c43356043003.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-12e75f2f-06f5-40a6-ba42-5880235f944c.jpg)

(b)

Figure 2: Comparison between the semi-analytic solution, eq. 1 and an ECHO-QGP (2 + 1) GLYPH&lt;0&gt; D simulation within the Bjorken coordinate metrics. The test is performed letting the fluid evolution be initialized with a symmetric energy density profile together, created ad hoc. The simulation has been performed with a grid of 0.025 fm in space and 0.001 fm in time. The shear viscosity to entropy density ratio is set to GLYPH&lt;17&gt;= s = 0 2, while the shear relaxation time is : GLYPH&lt;28&gt; R = 5 GLYPH&lt;17&gt;= ( e + P ). The energy scale is fixed at q = 1 fm GLYPH&lt;0&gt; 1 while the equation of state must be the one of a conformal fluid, we choose e = 3 P .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-e6784e0c-24ad-4922-a51e-bb22871ad880.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0186v2-figure-7619b4e0-990b-4f32-b284-037b7d9354ba.jpg)

The excellent agreement between the two can be appreciated for every component of the shear stress tensor (we show here just GLYPH&lt;25&gt; xx , GLYPH&lt;25&gt; xy and GLYPH&lt;25&gt; GLYPH&lt;17&gt;GLYPH&lt;17&gt; , but all the components present the same agreement), and for the thermodynamic variables (again we just show the temperature). For the sake of clearness, we show here just two close time steps ( GLYPH&lt;28&gt; = 2 and GLYPH&lt;28&gt; = 4), but the quality of the result is maintained up to much higher times. We remark then ECHO-QGP reproduces the semi-analytic solution without any need of fine-tuning.

- [14] P. Huovinen, H. Petersen, EPJ A48 (2012) 171. arXiv:1206.3371 doi:10.1140/epja/i2012-12171-9 , .
- [15] S. S. Gubser, A. Yarom, Nucl.Phys. B846 (2011) 469-511. arXiv:1012.1314 doi:10.1016/j.nuclphysb.2011.01.012 , .
- [16] S. S. Gubser, Phys.Rev. D82 (2010) 085027. arXiv:1006.0006 doi:10.1103/PhysRevD.82.085027 , .
- [17] H. Marrochio, J. Noronha, G. S. Denicol, M. Luzum, S. Jeon, et al., arXiv:1307.6130 .
- [18] E. Molnar, H. Niemi, D. Rischke, EPJ C65 (2010) 615-635. arXiv:0907.2583 doi:10.1140/epjc/s10052-009-1194-9 , .
- [19] S. Floerchinger, U. A. Wiedemann, A. Beraudo, L. Del Zanna, G. Inghirami, V. Rolando, arXiv:1312.5482 .
- [20] S. Floerchinger, U. A. Wiedemann, arXiv:1405.4393 .
Valentina Rolando Gabriele Inghirami Andrea Beraudo Luca Del Zanna Francesco Becattini Vinod Chandra Arturo De Pace Marzia Nardi
Heavy Ion Collision evolution modeling with ECHO-QGP
We present a numerical code modeling the evolution of the medium formed in
relativistic heavy ion collisions, ECHO-QGP. The code solves relativistic
hydrodynamics in $(3+1)-$D, with dissipative terms included within the
framework of Israel-Stewart theory; it can work both in Minkowskian and in
Bjorken coordinates. Initial conditions are provided through an implementation
of the Glauber model (both Optical and Monte Carlo), while freezeout and
particle generation are based on the Cooper-Frye prescription. The code is
validated against several test problems and shows remarkable stability and
accuracy with the combination of a conservative (shock-capturing) approach and
the high-order methods employed.
  In particular it beautifully agrees with the semi-analytic solution known as
Gubser flow, both in the ideal and in the viscous Israel-Stewart case, up to
very large times and without any ad hoc tuning of the algorithm.