Anisotropic flows and the shear viscosity of the QGP within an event-by-event massive parton transport approach

We have developed an event-by-event relativistic kinetic transport approach to study the build up of the anisotropic flows $$v_{n}(p_T)$$vn(pT) for a system at fixed $$\eta /s(T)$$η/s(T). The partonic approach describe the evolution of massless partons which imply $$\epsilon =3p$$ϵ=3p as Equation of State (EoS). We extend previous studies to finite partonic masses tuned to simulate a system that expand with an EoS close to the recent lQCD results. We study the role of EoS and the effect of $$\eta /s(T)$$η/s(T) ratio on the build up of $$v_n(p_T)$$vn(pT) up to $$n=5$$n=5 for two beam energies: RHIC energies at $$\sqrt{s}=200$$s=200 GeV and LHC energies at $$\sqrt{s}=2.76$$s=2.76 TeV. We find that for the two beam energies considered the suppression of the $$v_n(p_T)$$vn(pT) due to the viscosity of the medium have different contributions coming from the cross over or QGP phase. We shows that in ultra-central collisions (0–0.2%) the $$v_n(p_T)$$vn(pT) have a stronger sensitivity to the T dependence of $$\eta /s$$η/s that increases with the order of the harmonic n. Finally, we discuss the results for the integrated flow harmonics $$\langle v_{n} \rangle $$⟨vn⟩ in ultra-central collisions pointing-out how the relative strength of $$\langle v_{n} \rangle $$⟨vn⟩ depend on the colliding energies as well as on the freeze-out dynamics.


Introduction
Ultra-relativistic heavy-ion collisions (uRHICs) experiments conducted at the relativistic heavy-ion collider (RHIC) and at large hadron collider (LHC) are the only tools to access experimentally the properties of hot strongly interacting matter. In the last decades it has been reached a general consensus confirmed by several experimental signatures and supported by several theoretical calculations that the matter created in these collisions consists of a strongly interacting quark-gluon plasma (QGP) rather than hadronic matter. One of the most important measurement was the large elliptic flow v 2 ( p T ) = ( p 2 x − p 2 y )/( p 2 x + p 2 y ) observed [1,2]. The a e-mail: salvatore.plumari@ct.infn.it elliptic flow v 2 is a measure of the anisotropy in momentum space of the emitted particles. It is an observable that encodes information about the Equation of State (EoS) and transport properties of the matter created in these collisions [3][4][5]. The comparison between theoretical calculations and experimental results first within the viscous hydrodynamics framework [6][7][8] and in the recent years also within kinetic transport approach [9][10][11][12][13] have shown that the large value of v 2 is consistent with a matter with a very low shear viscosity to entropy density ratio 4πη/s ∼ 1 − 2. This value is close but larger then the conjectured lower bound for a strongly interacting system, η/s = 1/4π [14].
In the recent years it has been possible to measure the event-by-event angular distribution of emitted particles. These measurements have made possible to extend this analysis to higher order harmonics v n ( p T ) [15][16][17]. The origin of these high order harmonics flows is attributed to the fluctuations in the initial geometry [18][19][20][21][22][23][24]. In the recent years most of the theoretical efforts has been focused on the development of the event-by-event viscous hydrodynamics simulations. The comparison between event-by-event viscous hydrodynamical calculations and the experimental results for v n ( p T ) seems to confirm a finite but not too large average value of η/s with 4π η/s ∼ 1−3 [21,22]. Early works, most hydrodynamic simulations of heavy ion collisions, have assumed a constant η/s but a small value of η/s is not an evidence of the creation of a QGP phase formation. On the other hand, it is known that the η/s of the QGP is expected to have a temperature dependence with a minimum close to the cross over region [25][26][27][28][29]. A phenomenological estimation of its temperature dependence could give further information if the matter created in these collisions undergoes a phase transition [30][31][32].
In the recent years, the CMS, ATLAS, and ALICE collaborations have been able to measure the anisotropic flow coefficients v n in ultra-central heavy ion collisions [17,33,34]. These experiments have shown that the triangular flow n = 3 appears to be comparable or even larger of the elliptic flow n = 2. Such a result poses a problem to hydrodynamical approach that predict a significantly larger v 2 w.r.t. v 3 in ultra-central collisions. This discrepancy is even more intriguing considering that it occurs for ultra-central collisions where hydrodynamics is expected to work at best. In the recent years, several theoretical calculations within hydrodynamic simulations have shown that in ultra-central collisions the correlations that takes place between the initial eccentricities n and the corresponding final anisotropic flows v n are the highest than other centralities. On the other hand, recently within a transport approach it has been pointed-out that at LHC energies and for ultra-central collisions a quite large correlation between n and v n similar to viscous hydrodynamics is present up to n = 5 at variance with RHIC energies and/or other centralities [24]. At the same time a much larger sensitivity to η/s(T ) has been spot. Therefore the study of the anisotropic flows in ultra-central collision offer a powerful tool infer information about the initial geometry of the fireball.
In this paper we study the role of the Equation of State and η/s(T ) on the build up of v n ( p T ) up to the order 5 by using a kinetic transport approach with initial state fluctuations. The paper is organized as follows. In Sect. 2, we introduce the transport approach at fixed shear viscosity to entropy density η/s. In Sect. 3, the implementation of the initial state fluctuations in the transport approach. In Sect. 4, we study the role of the EoS on the anisotropic flows and the time evolution of the anisotropic flows v n . Finally in Sect. 5 we study the effect of the η/s(T ) on the differential v n ( p T ). In this paper we will show results on v n ( p T ) for n = 2, . . . , 5 for two different systems: Au + Au collisions at √ s = 200 GeV and Pb+ Pb collisions at √ s = 2.76 TeV at different centralities.

Transport approach at fixed η/s
In this work the evolution of the matter created in these collisions is studied employing the kinetic transport theory. This study was performed using a relativistic transport code developed to perform studies of the dynamics of heavy-ion collisions at both RHIC and LHC energies [9,11,12,[35][36][37][38][39][40].
Recently this transport code has been extended to include the initial state fluctuations in order to study the elliptic flow v 2 and high order harmonics v n ( p T ) with n > 2. For a more detailed discussion about the implementation of the initial state fluctuation see [24]. The development of the partonic transport equation has been done only for massless partons [9,11,12]. This implies that the matter simulated has an equation of state (EoS) − 3 p = 0 which is different to the one evaluated by lQCD [41,42] which clearly exhibit a large trace anomaly (T μ μ = − 3 p = 0) and consequently a sound velocity c 2 S (T ) that is significantly smaller than 1/3 also in the initial stages of a Heavy Ion Collisions reaching about 1/10 at T ≈ T C = 155 MeV . In this paper we extend the solution of the Relativistic Boltzmann Transport (RBT) equation to massive partons which allows to simulate a fluid with an EoS similar to the recent lQCD calculations. Therefore, the evolution of the phase-space distribution function f (x, p, t) is given by solving the following RBT equation: with p μ p μ = m 2 and with f 0 the initial distribution function (that can be of non-equilibrium). While C[ f ] is the Boltzmann-like collision integral. In this paper we have considered only the 2 ↔ 2 processes. The C[ f ] can be written as, where σ (s, Θ) is the differential cross section which is related to the total cross section by σ tot = dΩσ (s, Θ).
The transport theory permits to study the effects of the microscopical processes fixed by matrix element or cross sections on the observables. This is however not the aim of this work where instead we determine an effective cross section σ tot as a tool to fix the η/s of the system. This approach is inspired by the success of the hydrodynamical approach that has shown the key role played by the η/s to describe the experimental data [21,22,43]. Therefore we use the RBT equation as an approach converging to hydrodynamics for small scattering relaxation time τ ∼ 1/σρ (or small η/s) [44][45][46]. On the other hand, it permits to study directly the impact of a T dependence of η/s on observables like the anisotropic flows v 2 ( p T ) and v 3 ( p T ) which is the main focus of this paper. However, this approach it is not based on an ansatz for the viscous corrections for the phase-space distribution function δ f and it is not affected by limitation in the transverse momentum range in order to ensure that δ f / f << 1. Also the RBT approach is naturally valid also at large η/s or for p T >> T and it is possible to include initial non equilibrium effects to the initial distribution function f (x, p) [37,38,47]. Furthermore, as will be discussed below, the kinetic freeze-out can be determined self-consistently by increasing η/s(T ) at low temperature which gives a smooth switching-off of the scattering rates. The current disadvantage of the present approach is that hadronization has not yet been included and it will be the main topics of forthcoming studies. In the transport calculations shown in this paper we use massive particles providing a soften equation of state respect to the massless case with a decreasing speed of sound when the crossover region is approached similar to the one of lQCD calculations [41,42].
To study the expansion dynamics of the fireball with a certain η/s(T ), we determine locally in space and time the The orange area refers to the quasi-particle model predictions for η/s [51]. The three different lines indicate different possible T dependencies used in this paper. Symbols are as in the legend. See the text for more details total cross section σ tot according to the Chapmann-Enskog theory [48,49]. For the results shown in this paper we have considered an isotropic cross section. As shown in [48][49][50] for a massive gas in the Chapmann-Enskog approximation at the first order η is given by with z = m/T while the function f (z) is given by where K n -s are the modified Bessel functions. The entropy density for a massive system is given by s = ρ 4 + zK 1 (z)/K 2 (z) . Notice that in the ultra-relativistic limit z → 0 the f (z) → 1.2 and we recover the massless limit for the η. As shown in [36] the expression for η in Eq. (3) is in quite good agreement at level of 3−5% with the Green-Kubo formula and it describes correctly the η/s of the system in the range of interest of HIC. Notice that in this limit the η/s ratio coincides with the one obtained in the relaxation time approximation. In Fig.1 it is shown a collection of different theoretical calculations about the T dependence of the η/s ratio. In literature there are several indications that η/s should have a particular behaviour with the temperature [25][26][27]30,31,52]. As shown in Fig.1 in general η/s should have a typical behaviour of phase transition with a minimum close to the critical temperature T C [30][31][32]. Estimations of η/s within chiral perturbation theory for a meson gas [25,26], have shown that at temperature lower than the critical temperature the η/s is a decreasing function with the temperature, see down-triangles in Fig.1. Similar results have been obtained from extrapolation of heavy-ion collisions data at intermediate energies, see HIC-IE diamonds in Fig.1. At higher temperature T > T c lQCD calculation on pure gauge have shown that in general η/s becomes an increasing function with the temperature [27][28][29], see up-triangles, circles and squares in Fig.1. However due to the large error bars in the lQCD results for η/s at moment it is not possible to infer a clear temperature dependence in the QGP phase.
In Fig.1 we show the three different T dependencies of η/s studied in this paper. In particular we have considered the following cases: one of constant 4πη/s = 1 during all the evolution of the system shown by dot dashed line in Fig.1 another one of 4πη/s = 1 at higher temperature in the QGP phase and an increasing η/s in the cross over region towards the estimated value for hadronic matter 4πη/s ≈ 6 [26,53] and shown by solid line. The third one is shown in Fig.1 by the dot dashed line. In this case we consider the increase of η/s at higher temperature with a linear temperature dependence and a minimum close to the critical temperature with a temperature dependence similar to that expected from general considerations as suggested by quasi particle models [51,54] or by recent calculations on pure gauge [29]. Notice that in our approach an increase of η/s at lower temperature 0.8T C ≤ T ≤ 1.2T C permits to account for a smooth realistic kinetic freeze-out (f.o.) because at lower T the mean free path λ goes like λ ∝ η s 1 T and the total cross section goes like σ ∝ (η/s) −1 . In the following discussion the term f.o. means that we take into account the increase of η/s at T < 1.2T C .

Initial conditions
In the following paper we will consider two systems at different centralities: Au+ Au collisions at √ s N N = 200 GeV produced at RHIC and Pb+ Pb collisions at √ s N N = 2.76 TeV at LHC. In this section we discuss the implementation of the initial state fluctuations in the transport approach, for more details about the implementation of the initial fluctuations in our approach see Ref. [24]. To generate an event by event initial profile we use the Monte-Carlo Glauber model. For a fixed centrality class we determine the average impact parameter by the Monte Carlo Glauber model. We use this to generate event by event N part according to the Monte Carlo Glauber model. We employ the geometrical method to determine the the initial number of participant nucleons. Within this method two nucleons collide each other if the relative distance in the transverse plane is  [8,19,20,43,55]. In our calculations we convert the information of the nucleon distribution into parton density in the transverse plane ρ T (x, y) by the following expression where N part is the number of participants and k is an overall normalization factor fixed in order to reproduce the final hadron multiplicity d N ch /dy. In the results shown in this paper we have fixed the Gaussian width of the fluctuations to σ xy = 0.5 f m. In our calculation we assume a longitudinal boost invariant distribution from y = −2.5 to y = 2.5. The transverse density ρ T (x, y) fixes in each event the initial anisotropy in coordinate space that is quantified in terms of the following coefficients n [19,20]: where r T = x 2 + y 2 and φ = arctan(y/x) is the polar coordinate in the transverse plane. In Fig. 2 it is shown the initial spatial anisotropies n for n = 2, . . . , 5 as a function of the impact parameter. In the inset panel it is shown n as a function of the order of the harmonic for central collision. In particular, as shown in the inset panel, for central collision larger is the order of the harmonics larger is the corresponding initial eccentricities. However this is a general behaviour for the eccentricities for central collisions and similar ordering happens also for different bin selections. For the initialization in momentum space at RHIC (LHC) energies we have considered for partons with transverse momentum p T ≤ p 0 = 2 GeV (3 GeV) a thermalized spectrum in the transverse plane. Assuming the local equilibrium the initial local temperature in the transverse plane T (x, y) is evaluated by using the standard thermodynamical relation for massive Fig. 3 it is shown the initial temperature profile in the transverse plane for a given event and for Au + Au collisions at While for partons with p T > p 0 we have assumed the spectrum of non-quenched minijets according to standard NLO-pQCD calculations with a power law shape [56,57]. The initial transverse momentum of the particles is distributed uniformly in the azimuthal angle. We fix the starting time of the simulation to τ 0 = 0.6 fm/c for RHIC and τ 0 = 0.3 fm/c for LHC as commonly done also in hydrodynamical approaches.
In our simulations we have used N event = 1000 events for each centrality class. This number is enough to get stable results for the spectra, differential elliptic flow and high order flow coefficients v n ( p t ) for n = 2, . . . , 5. Moreover, in the calculation shown in this paper have checked the convergence of our results for v n with the lattice spacing of the calculation grid and the number of the test particles N test . In the results hown in this paper we have used a calculation grid with a transverse area of the cell A T = 0.12 fm 2 and the total number of test particles is fixed to be N test = 2 · 10 6 for each event.

Role of the equation of state on the v n
As well known the anisotropic flow coefficients v n ( p T ) are observable sensitive to the bulk properties of the created matter and they encode information about the EoS and transport properties of the medium like for example gives information about η/s ratio. In this section we discuss the role of the EoS on the build up of the v n ( p T ). To extract information about the EoS we evaluate the average energy-momentum tensor.
We evaluate the energy-momentum tensor in a cylinder of radius r T = 1.5 fm and in longitudinal direction |η| < 0.3. We evaluate the energy density by = T 00 and the pressure by assuming that p = 1 3 i T ii . In Fig. 4 it is shown the comparison of p/ ratio obtained dynamically in our simulations with the one of lQCD for both RHIC energies in Au + Au collisions at √ s N N = 200 GeV and at LHC energies in Pb + Pb collisions at √ s N N = 2.76 TeV for central collisions. The results shown a reasonably agreement with lQCD data for both energies. In these calculations we have fixed m = 0.5 GeV that are in agreement with the masses extracted within quasi-particle models [51,59]. As shown, this implies that in our simulations dynamically we describe a systems that approximatively have the equation of state as the one in lQCD calculations. The initial fluctuations in density induce fluctuations in energy density and the corresponding value of p/ can fluctuate event-by-event. In Fig. 4 the blue and green bands show the possible range of trajectories covered in our simulations. The band is due to the fact that locally the relation between pressure and energy density can be modified by non equilibrium effects.
In this paper the anisotropic flows v n ( p T ) with n = 2, . . . , 5 have been calculated using the event plane method as where the momentum space angles Ψ n are given by while the eccentricities are evaluated using Eq. (6) as done in [8,19,20].
As a first step we study the time evolution of the initial eccentricities n and the corresponding integrated anisotropic flows v n . Where the average ... is over all the events.
In the left panel of Fig. 5 we plot the time evolution of the initial anisotropies in coordinate space n for both cases massless (dashed lines) and massive case (solid lines). In the right panel of Fig. 5 instead we show the time evolution for the corresponding final anisotropic flow coefficients v n . These calculation are for LHC energies at mid rapidity and for (20-30)% centrality collisions. Results show that in general for both massless and massive case the v n ( p T ) are built-up within the first few fm/c and for both cases different harmonics have a different formation time with the following ordering t v 2 < t v 3 < t v 4 · · · as shown in Ref. [24], but here we consider the impact of the EoS. Moreover as shown in  Fig. 5 we observe that for massless case the system is more efficient to convert the initial anisotropy in coordinate space into final anisotropy in momentum space where the effect of a realistic EoS is to reduce the v n . This effect increase with the increasing of the order of the harmonics. The larger effect for larger n is related to the fact that higher harmonics have larger formation time and therefore they develop in a region where the p/ or the speed of sound is smaller. In fact for example the second coefficient v 2 is almost completely developed within 3-4 fm/c where the p/ ≈ 0.3 and close to the conformal limit and the effect increases for large n because higher harmonics develop later at lower temperature where p/ < 0.3. In Fig. 6 we compare the final v n ( p T ) for both massless (dashed lines) and massive case (solid lines). These results are for Pb + Pb collisions at √ s N N = 2.76 TeV and for (20-30)% centrality collisions. We observe that for the elliptic flow v 2 a mass ordering at low transverse momentum where larger is the mass smaller is the corresponding elliptic flow which is typical of hydrodynamic expansion that boost strongly massive particles [4,60]. The explanation of the mass ordering is due to the fact that at low p T the elliptic flow is v 2 ( p T ) ∝ p T − β T m T therefore for light particles m T = p 2 T + m 2 ≈ p T and v 2 ( p T ) ∝ p T while for massive particle the coupling of m T with the radial flow gives at the same p T a smaller v 2 ( p T ) and it is translated in a non linear behaviour at low p T . Moreover a similar mass ordering it is shown also for higher harmonics v n ( p T ). At low p T the effect is to give a non linear behaviour for the second and third coefficients. Moreover we observe that the sensitivity to the EoS of the v n ( p T ) increases with the increasing of

Effects of η/s(T ) on the v n ( p T )
Differential anisotropic flow coefficients v n ( p T ) are observables that carry out more information about the fireball created in the heavy ion collisions. Moreover they are observables sensitive to the transport properties of the medium like the η/s ratio. In the following discussion we will study the effect of the η/s on the build up of the elliptic flow v 2 ( p T ) and on the high order harmonics v 3 ( p T ), v 4 ( p T ) and v 5 ( p T ). In particular we study two different systems the one created at RHIC for Au + Au collisions at √ s = 200 GeV and the one at LHC energies for Pb + Pb at √ s = 2.76 TeV. The following calculations are for the massive case which gives a reasonable agreement with the lQCD EoS, see Fig. 4.
In Fig. 7 we show the effect of the η/s(T ) on the v n ( p T ) at mid rapidity and for (20-30)% centrality for both RHIC energies for Au + Au collisions at produced v n ( p T ). These results are qualitatively in agreement with previous results within kinetic transport approach in the conformal limit [24]. Comparing the dot dashed lines with the solid lines in the left panel of Fig. 7, we observe that at RHIC energies the v n ( p T ) are sensitive to the increase of the η/s at lower temperature close to the cross over region. As one should expect by the different formation time, this sensitivity increase with the order of the harmonics. On the other hand at LHC energies, right panel of Fig. 7, the v n ( p T ) are more sensitive to the increase of η/s at high temperature as shown by a larger suppression of the v n ( p T ) for the case where η/s ∝ T in the QGP phase (dashed lines).
In the recent years, it has been possible to access experimentally to the ultra-central collisions. In ultra-central collisions the initial asymmetry in coordinate space measured by n comes from the fluctuations in the initial geometry and there is no coupling with the global overlap shape. Moreover, it has been shown within event-by-event ideal and viscous hydrodynamics calculations [55,61,62] that for these collision centralities the final integrated anisotropic flows coefficient v n are strongly correlated to the initial eccentricities n with a linear correlation coefficient C( n , v n ) ≈ 1. Similar results have been obtained also within event-by-event transport approach [24].
In Fig. 8 Fig. 7 with Fig. 8 one can observe that at low p T both centralities the v n ( p T ) ∝ p n T while at hight p T for ultra-central collisions we observe that the elliptic flow v 2 ( p T ) shows a saturation while v n ( p T ) for n ≥ 3 have an almost linear p T dependence. These results are in qualitative agreement with what has been observed experimentally. Similarly to the results at mid peripheral collisions at RHIC energies the v n ( p T ) are more sensitive to the value of η/s at lower temperature. Moreover, we observe that for central collisions the v n ( p T ) are more sensitive with respect to more peripheral collisions with a suppression of about 30-35% for n > 3 in central collision and of about 20-25% for mid peripheral collisions. On contrary at LHC energies we observe that reduction of v n ( p T ) due to the increase of η/s in the QGP phase (dashed lines) is strongly enhanced for ultra-central collisions. In particular for n ≥ 3 the reduction for central collisions is about 30-35% against a reduction of about 10% for mid peripheral collisions. It is indeed remarkable notice that the 30% effect is determined by a slowly linear rising of η/s with T and with a minimum close T C of η/s ≈ 1/(4π) as the dashed lines in Fig. 1. Moreover in central collisions at LHC energies higher harmonics acquire a larger sensitivity to the value of the viscosity in the QGP phase. Such a sensitivity is not present at RHIC energy and at LHC in non-central collisions. On the other hand this does not appear for the n but only at p T 1 GeV.
The comparison between RHIC and LHC highlights the role played by the EoS and the role of the collision energy 76 TeV corresponding to 0-0.2% centrality. The experimental measurement have shown that the triangular flow n = 3 appears to be comparable or even larger of the elliptic flow n = 2. Such a result poses a problem to hydrodynamical approach that predict a larger v 2 w.r.t. v 3 in ultra-central collisions. In Fig. 9 it is shown v n as a function of the order of the harmonic n for ultra-central collisions and for Au + Au collisions at √ s N N = 200 GeV (left panel) and Pb + Pb collisions at √ s N N = 2.76 TeV (right panel) As shown in the right panel of Fig. 9 by comparing circles with diamonds the effect of a soften Equation of State is to reduce the build-up of the v n . In particular we observe a greater suppression for higher harmonics where for example for the massless case we get that v 4 ≈ v 2 while for the massive case v 4 < v 2 with a suppression of about 15%. On the other hand the role of the collision energy is to reduce the production of the anisotropic flows and moreover it can play a role on the peak for n = 3. The different behaviour between RHIC and LHC is due to the fact that the degree of correlation between the initial n and the final anisotropic flows produced v n is different. In fact the linear correlation coefficient between them C( n , v m ) is a decreasing function with the collision energy and impact parameter as shown for massless partons in [24]. This implies that for ultra-central collisions and at LHC energies the degree of correlation is maximum with C( n , v m ) ≈ 1 for n = 2, . . . , 4. Therefore at higher energies like at LHC energies the v n keep more information about the initial anisotropies in the coordinate space and in particular they tend to keep the initial ordering of the eccentricities where for ultra-central collisions we have 2 < 3 < 4 · · · , see Fig. 2. At LHC it has been observed a v 2 v 3 for ultra-central collisions. Our simulation show that such a feature should disappear at RHIC because of the impact of the increasing η/s(T ) in the crossover region. On the other hand at LHC a v 2 v 3 survives because the freeze-out dynamics implied by η/s(T ) does not affect the v n at high energy.

Conclusions
In this paper we have investigated the build up of the anisotropic flows v n ( p T ) for n = 2, 3, 4 and 5 within an event-by-event transport approach. We have studied the role of the EoS, by using a finite partonic mass, on the time evolution of v n and the initial eccentricities n and on the differential v n ( p T ). We have found that the effect of the mass is to give a non linear behaviour at low p T for the elliptic and triangular flow. In general the mass has the effect to reduce the final v n ( p T ) produced. The system is less efficient in converting the initial anisotropy in coordinate space into final anisotropic flows. Due to the different formation time of the harmonics the efficiency is lower for higher harmonics that start to develop at smaller speed of sound. Moreover, we have also studied the effect the temperature dependence of the η/s ratio on v n ( p T ) for two different beam energies: at RHIC for Au + Au collisions at √ s = 200 GeV and at LHC for Pb + Pb collisions at √ s = 2.76 TeV. We have found that at RHIC the v n ( p T ) are more affected by the value of η/s in the cross over region (T < 1.2T C ) and the sensitivity increases with the order of the harmonics. At LHC energies we get that almost all the v n ( p T ) develop in the QGP phase and they have less contamination of the value of η/s in the cross-over region and we observe a weaker sensitivity to the T dependence of the η/s. These results are in qualitative agreement with the results obtained for the massless case. Such a scenario changes for ultra-central collisions. In particular in these collision an enhancement of the sensitivity to the value of η/s at lower temperature for RHIC is observed while at LHC the v n ( p T ) are more sensitive to the value at high temperature. In general we found an enhancement of the sensitivity of the v n ( p T ) for n = 2, . . . , 5 that can reaches about a 30-35%.