## What the collective flow excitation function can tell about the quark-gluon plasma

Jussi Auvinen ∗

Frankfurt Institute for Advanced Studies (FIAS), Ruth-Moufang-Strasse 1, D-60438 Frankfurt am Main, Germany

Jan Steinheimer † and Hannah Petersen ‡

Frankfurt Institute for Advanced Studies (FIAS), Ruth-Moufang-Strasse 1, D-60438 Frankfurt am Main, Germany and Institut f¨r u Theoretische Physik, Goethe Universit¨t, a Max-von-Laue-Strasse 1, D-60438 Frankfurt am Main, Germany

## Abstract

Recent STAR data from the RHIC beam energy scan (BES) show that the midrapidity slope dv /dy 1 of the directed flow v 1 of net-protons changes sign twice within the collision energy range 7.7 - 39 GeV. To investigate this phenomenon, we study the collision energy dependence of v 1 utilizing a Boltzmann + hydrodynamics hybrid model. Calculations with dynamically evolved initial and final state show no qualitative difference between an equation of state with a cross-over and one with a first-order phase transition, in contrast to earlier pure fluid predictions. Furthermore, our analysis of the elliptic flow v 2 shows that pre-equilibrium transport dynamics are partially compensating for the diminished elliptic flow production in the hydrodynamical phase at lower energies, which leads to a qualitative agreement with STAR BES results in midcentral collisions. No compensation from transport is found in our model for integrated v 3 , which decreases from ≈ 0 02 at . √ s NN = 27 GeV to ≈ 0 005 at . √ s NN = 7 7 GeV in midcentral collisions. .

PACS numbers: 24.10.Lx,24.10.Nz,25.75.Ld

∗ Electronic address: auvinen@fias.uni-frankfurt.de

† Electronic address: steinheimer@fias.uni-frankfurt.de

‡ Electronic address: petersen@fias.uni-frankfurt.de

## I. INTRODUCTION

In 2010, a beam energy scan program was launched at the Relativistic Heavy Ion Collider (RHIC) to study the features of the QCD phase diagram, in particular the location and the order of the phase transition between the hadronic and QCD matter in the plane of baryochemical potential µ B and temperature T . The main observables are the coefficients v n of the Fourier expansion of the azimuthal angle distribution of final-state particle momenta, which are typically associated with collective flow.

Fluid dynamical calculations have predicted that the slope of the directed flow v 1 of baryons will turn negative and then positive again as a function of energy if there is a first order phase transition between hadronic and QCD matter [1-4]. Qualitatively similar behavior of the midrapidity slope of the net-proton directed flow, dv /dy 1 , has been found by STAR experiment in the RHIC beam energy scan, with a minimum in the energy interval √ s NN = 11 5 . -19 6 GeV [5]. . The elliptic flow v 2 is one of the key observables supporting the formation of a strongly coupled quark-gluon plasma (QGP) at the highest energies of RHIC and the Large Hadron Collider (LHC). However, the measured differential elliptic flow v 2 ( p T ) for charged hadrons remains nearly unchanged from the collision energy range √ s NN = 39 GeV down to 7.7 GeV [6], where the formation of hydrodynamically evolving QCD matter is expected to be considerably diminished compared to top RHIC energies. In addition, the preliminary STAR data suggests that the magnitude of the triangular flow v 3 remains constant at lower collision energies [7].

We study the collision energy dependence of the collective flow in the RHIC BES range with a hybrid approach, where the non-equilibrium phases at the beginning and in the end of a heavy-ion collision event are described by a transport model, while a hydrodynamic description is used for the intermediate hot and dense stage and the phase transition between the QGP and hadronic matter. This provides a consistent framework for investigating both high-energy heavy ion collisions with negligible net-baryon density and notable quark-gluon plasma phase, and the collisions at smaller energies with finite net-baryon density, where less QGP is expected to form. This approach should thus be ideal for beam energy scan studies.

## II. HYBRID MODEL

In this study, the transport + hydrodynamics hybrid model by Petersen et al. [8] is utilized. The transport model describing the initial and final state is the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) string / hadron cascade [9, 10]. The transition time to hydrodynamics is defined as the moment when the two colliding nuclei have passed through each other: t start = 2 R/ √ γ 2 CM -1, where R represents the nuclear radius and γ CM = 1 / √ 1 -v 2 CM is the Lorentz factor.

The (3+1)-D ideal hydrodynamics evolution equations are solved with the SHASTA algorithm [11, 12]. The equation of state (EoS) which is utilized most of this study is from Steinheimer et al. [13]. It is a combination of a chiral hadronic and a constituent quark model and has the important feature of being applicable also at finite net-baryon densities found at lower collision energies.

The transition from hydro to transport (aka 'particlization'), is done when the energy density glyph[epsilon1] reaches the critical value glyph[epsilon1] C = nglyph[epsilon1] 0 , where glyph[epsilon1] 0 = 146 MeV/fm 3 is the nuclear ground state energy density. In this study, the values n = 2 and n = 4 are used. The particle distributions are generated according to the Cooper-Frye formula from the iso-energy density hypersurface, which is constructed using the Cornelius hypersurface finder [14]. The final rescatterings and decays of these particles are then computed in the UrQMD.

## III. RESULTS

To investigate the sensitivity of the directed flow v 1 ( y ) = 〈 p x p T 〉 y i = y on the order of the phase transition, we run the simulations with a first-order phase transition 'Bag model' EoS as an alternative to the above described chiral model EoS which has a cross-over phase transition. The transition from fluid dynamics back to transport happens on an iso-energy density glyph[epsilon1] C = 4 glyph[epsilon1] 0 ≈ 0 6 GeV . / fm 3 hypersurface.

To emulate the earlier fluid calculations, we first utilize a cold nuclear matter initialization, where the colliding nuclei are represented by two distributions of energy and baryon density, which respect boosted Woods-Saxon profiles with a central density of saturated nuclear matter ρ 0 ≈ 0 16 fm . -3 . The starting point of the simulation is just before the two nuclei first make contact; in the early stage of the collision the kinetic energy of the nuclei is

then transformed into large local densities. Figures 1a and 1b show the difference in dv /dy 1 between the two equations of state with a cold nuclear matter initialization and the UrQMD afterburner for Au+Au collisions at impact parameter b = 8 fm. The predicted minimum in dv /dy 1 as a function of √ s NN with a first-order phase transition is clearly observed when using isochronous particlization condition; however, the difference between the two equations of state diminishes greatly when using iso-energy density fluid-to-particles switching condition instead.

Figure 1c shows the result of the full hybrid simulation with the initial non-equilibrium transport phase for the energy dependence of midrapidity slopes of proton and antiproton v 1 . The directed flow was calculated using events with impact parameter b = 4 6 . -9 4 . fm, to approximate the (10 -40)% centrality range of the STAR data. The hybrid model overestimates the experimental data and also the pure UrQMD transport result in the whole examined collision energy range. The two EoS are completely indistinguishable in the hybrid simulations, questioning the validity of v 1 as a signal of the first-order phase transition. Some possible sources for the difference between the model and the experimental data are the momentum transfer between the spectator particles and the fireball, which is not accounted for in this study, and the method used to determine the event plane, as here v 1 was calculated with respect to the reaction plane of the simulation. These uncertainties need further investigation before drawing definite conclusions.

Utilizing the hybrid model with the crossover EoS, the flow coefficients v 2 and v 3 are calculated from the particle momentum distributions using the event plane method [18, 19]. Figure 2a shows the elliptic flow v 2 produced in Au+Au -collisions, integrated over the p T range 0.2 - 2 GeV, compared with the STAR data for the (0-5)%, (20-30)% and (30-40)% centrality classes. In the model these are respectively represented by the impact parameter intervals b = 0 -3 4 fm, . b = 6 7 . -8 2 fm and . b = 8 2 . -9 4 fm. . Critical energy density glyph[epsilon1] C = 2 glyph[epsilon1] 0 is used here, as this value has been found to give a reasonable agreement with the experimental data for particle m T spectra at midrapidity | y | &lt; 0 5 for energies ranging . from E lab = 40 AGeV to √ s NN = 200 GeV [21, 22]. Figure 2b demonstrates the magnitude of v 2 at three different times: just before the hydrodynamical evolution, right after the particlization, and the final result after the hadronic rescatterings have been performed in the UrQMD model.

In the impact parameter range b = 8 2 . -9 4 . fm the rescatterings contribute roughly

FIG. 1: a) and b) Slope of v 1 of protons around midrapidity | y | &lt; . 0 5 with an equation of state with first-order phase transition (black) and cross-over phase transition (red) for isochronous fluid-toparticles transition hypersurface (a) and iso-energy density transition hypersurface (b) for Au+Au collisions at impact parameter b = 8 fm. c) Midrapidity slope dv /dy 1 of protons (solid symbols) and anti-protons (open symbols) for impact parameter range b = 4 6 . -9 4 . fm, extracted from the hybrid model calculations with a bag model (black) and crossover EoS (red). Compared with standard UrQMD (grey) and experimental data [5, 15, 16] (green). Plots from [17].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0098v1-figure-de49dbf3-c478-4e41-86bd-36da5deb6e3e.jpg)

10% on the final result. The hydrodynamics produce very little elliptic flow at √ s NN ≤ 7 7 GeV; . v 2 below √ s NN = 10 GeV is in practice completely produced by the transport dynamics (resonance formations and decays, string excitations and fragmentation). This initial transport gains importance at lower energies due to the prolonged pre-equilibrium phase. On the other hand, above √ s NN = 19 6 GeV the hydrodynamics are clearly the . dominant source of v 2 .

The simulation results overshoot the experimental data for all collision energies. This suggests that the viscous corrections should be included - indeed, good results have already been achieved using similar hybrid approach with viscous hydro [23]. In the most central collisions below √ s NN = 11 5 GeV, the model deviates from the observed energy dependence; . it is likely that the abrupt change from the pre-equilibrium phase to hydro is not a valid assumption in this regime and a more dynamic procedure should be implemented. For the purposes of this study, however, the most important feature is the good qualitative agreement in the midcentral collisions, as here the flow effects are at their largest.

Based on the above results, it appears that the decrease in the hydrodynamically produced elliptic flow is partially compensated by the increased flow production in the transport phase, and so the observed v 2 has weaker collision energy dependence than one would have naively

FIG. 2: a) Integrated v 2 at midrapidity | η | &lt; 1 0 from hybrid model (circles) in collision energy . range √ s NN = 7 7 . -200 GeV at three impact parameter ranges, compared to the STAR data [6, 20] (stars). b) Integrated v 2 for impact parameter range b = 8 2 . -9 4 at . √ s NN = 5 -200 GeV, at the beginning of hydrodynamical evolution (diamonds), immediately after particlization (squares) and after the full simulation (circles). c) Integrated v 3 at midrapidity | η | &lt; 1 0 for . b = 0 -3 4 fm . and b = 6 7 . -8 2 fm at . √ s NN = 5 -200 GeV. Plots from [21].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0098v1-figure-34ee9771-2508-4fc3-a374-dfa81dc61255.jpg)

expected. To study this phenomenon further, we do the same analysis for the triangular flow v 3 , which originates purely from the event-by-event variations in the initial spatial configuration of the colliding nucleons.

As illustrated by Figure 2c, the p T -integrated v 3 increases from ≈ 0 01 to above 0.015 . with increasing collision energy in the most central collisions. However, in midcentrality b = 6 7 . -8 2 fm there is a rapid rise from . ≈ 0 at √ s NN = 5 GeV to the value of ≈ 0 02 . for √ s NN = 27 GeV, after which the magnitude remains constant. The energy dependence of v 3 in midcentral collisions qualitatively resembles the hydrodynamically produced v 2 in Figure 2b; the viscous medium described by transport smears the anisotropies in the initial energy density profile instead of converting them into momentum anisotropy, and is thus unable to compensate for the lack of ideal fluid described by hydrodynamics. This makes v 3 the clearer signal of the presence of low-viscous medium.

## IV. CONCLUSIONS

Contrary to the earlier fluid calculations, we have found no difference between an equation of state with a first order phase transition and one with a cross-over phase transition for the midrapidity slope of the directed flow v 1 , when utilizing the full hybrid model with an iso-energy density switching criterion between hydrodynamics and final transport. Thus

dv /dy 1 cannot currently be considered as a good signal for the existence of a first-order phase transition. However, there is currently a notable discrepancy between the model results and the experimental data which necessitates further investigation.

We have demonstrated that a hybrid transport + hydrodynamics approach can qualitatively reproduce the experimentally observed behavior of v 2 as a function of collision energy √ s NN . While the v 2 production by hydrodynamics is diminished at lower collision energies, this is partially compensated by the pre-equilibrium transport dynamics. Same does not apply to triangular flow v 3 , which decreases considerably faster, reaching zero in midcentral collisions at √ s NN = 5 GeV. Thus v 3 is the better signal for the formation of quark-gluon plasma in heavy ion collisions.

## V. ACKNOWLEDGEMENTS

This work was supported by GSI and the Hessian initiative for excellence (LOEWE) through the Helmholtz International Center for FAIR (HIC for FAIR). H.P. and J.A. acknowledge funding by the Helmholtz Association and GSI through the Helmholtz Young Investigator Grant No. VH-NG-822. The computational resources were provided by the LOEWE Frankfurt Center for Scientific Computing (LOEWE-CSC).

- [1] D. H. Rischke, Y. P¨ urs¨n, u J. A. Maruhn, H. St¨ ocker and W. Greiner, Heavy Ion Phys. 1 (1995) 309.
- [2] L. P. Csernai and D. R¨hrich, Phys. Lett. B 458 (1999) 454. o
- [3] J. Brachmann et al. , Phys. Rev. C 61 (2000) 024909.
- [4] H. St¨cker, Nucl. Phys. A 750 (2005) 121. o
- [5] L. Adamczyk et al. [STAR Collaboration], arXiv:1401.3043 [nucl-ex].
- [6] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. C 86 (2012) 054908.
- [7] Y. Pandit [STAR Collaboration] Talk at Quark Matter 2012.
- [8] H. Petersen, J. Steinheimer, G. Burau, M. Bleicher and H. St¨cker, Phys. Rev. C 78 (2008) o 044901.
- [9] S. A. Bass, M. Belkacem, M. Bleicher, M. Brandstetter, L. Bravina, C. Ernst, L. Gerland and

- M. Hofmann et al. , Prog. Part. Nucl. Phys. 41 (1998) 255.
- [10] M. Bleicher, E. Zabrodin, C. Spieles, S. A. Bass, C. Ernst, S. Soff, L. Bravina and M. Belkacem et al. , J. Phys. G 25 (1999) 1859.
- [11] D. H. Rischke, S. Bernard and J. A. Maruhn, Nucl. Phys. A 595 (1995) 346.
- [12] D. H. Rischke, Y. Pursun and J. A. Maruhn, Nucl. Phys. A 595 (1995) 383 [Erratum Nucl. Phys. A 596 (1996) 717].
- [13] J. Steinheimer, S. Schramm and H. St¨cker, Phys. Rev. C 84 (2011) 045208. o
- [14] P. Huovinen and H. Petersen, Eur. Phys. J. A 48 (2012) 171.
- [15] C. Alt et al. [NA49 Collaboration], Phys. Rev. C 68 (2003) 034903.
- [16] H. Liu et al. [E895 Collaboration], Phys. Rev. Lett. 84 (2000) 5488.
- [17] J. Steinheimer, J. Auvinen, H. Petersen, M. Bleicher and H. St¨cker, Phys. Rev. C 89 (2014) o 054913.
- [18] A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58 (1998) 1671.
- [19] J. Y. Ollitrault, Preprint nucl-ex/9711003 (1997).
- [20] J. Adams et al. [STAR Collaboration], Phys. Rev. C 72 (2005) 014904.
- [21] J. Auvinen and H. Petersen, Phys. Rev. C 88 (2013) 064908.
- [22] J. Auvinen and H. Petersen, PoS CPOD 2013 (2013) 034.
- [23] I. .A. Karpenko, M. Bleicher, P. Huovinen and H. Petersen, J. Phys. Conf. Ser. 503 (2014) 012040.
Jussi Auvinen Jan Steinheimer Hannah Petersen
What the collective flow excitation function can tell about the quark-gluon plasma
Recent STAR data from the RHIC beam energy scan (BES) show that the
midrapidity slope dv1/dy of the directed flow v1 of net-protons changes sign
twice within the collision energy range 7.7 - 39 GeV. To investigate this
phenomenon, we study the collision energy dependence of v1 utilizing a
Boltzmann + hydrodynamics hybrid model. Calculations with dynamically evolved
initial and final state show no qualitative difference between an equation of
state with a cross-over and one with a first-order phase transition, in
contrast to earlier pure fluid predictions. Furthermore, our analysis of the
elliptic flow v2 shows that pre-equilibrium transport dynamics are partially
compensating for the diminished elliptic flow production in the hydrodynamical
phase at lower energies, which leads to a qualitative agreement with STAR BES
results in midcentral collisions. No compensation from transport is found in
our model for integrated v3, which decreases from approx 0.02 at sqrt(sNN) = 27
GeV to approx 0.005 at sqrt(sNN) = 7.7 GeV in midcentral collisions.