Clustering of color sources and the shear viscosity of the QGP in heavy ion collisions at RHIC and LHC energies

We present our results on the shear viscosity to entropy ratio (η/s) in the framework of the clustering of the color sources of the matter produced at RHIC and LHC energies. The onset of de-confinement transition is identified by the spanning percolating cluster in 2D percolation. The relativistic kinetic theory relation for η/s is evaluated using the initial temperature (T) and the mean free path (λmfp). The analytic expression for η/s covers a wide temperature range. At T∼150 MeV below the hadron to QGP transition temperature of ∼168 MeV, with increasing temperatures the η/s value drop sharply and reaches a broad minimum η/s∼0.20 at T∼175–185 MeV. Above this temperature η/s grows slowly. The measured values of η/s are 0.204±0.020 and 0.262±0.026 at the initial temperature of 193.6±3 MeV from central Au+Au collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sqrt{s_{NN}}=200~\mathrm{GeV}$\end{document} (RHIC) and 262.2±13 MeV in central Pb+Pb collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sqrt{s_{NN}} = 2.76~\mathrm{TeV}$\end{document} (LHC). These η/s values are 2.5 and 3.3 times the AdS/CFT conjectured lower bound 1/4π but are consistent with theoretical η/s estimates for a strongly coupled QGP.


Introduction
The observation of the large elliptic flow at RHIC in noncentral heavy ion collisions suggest that the matter created is a nearly perfect fluid with a very low shear viscosity [1][2][3][4].
Recently, attention has been focused on the shear viscosity to entropy density ratio η/s as a measure of the fluidity [5][6][7][8]. The observed temperature averaged η/s, based on viscous hydrodynamics analyses of RHIC data, are suggestive of a strongly coupled plasma [9,10]. The effect of the bulk viscosity is expected to be negligible. It has been a e-mail: brijesh@purdue.edu conjectured, based on infinitely coupled super-symmetric Yang-Mills (SYM) gauge theory using the correspondence between Anti de-Sitter(AdS) space and conformal field theory (CFT), that the lower bound for η/s is 1/4π and is the universal minimal viscosity to entropy ratio even for QCD [11]. However, there are theories in which this lower bound can be violated [12]. In this work, we use the color string percolation model (CSPM) [13,14] to obtain η/s as a function of the temperature above and below the hadron to QGP transition. The measured η/s values are for Au + Au collisions at √ s NN = 200 GeV at RHIC and for Pb + Pb collisions at √ s NN = 2.76 TeV at LHC.

Clustering of color sources
Multiparticle production is currently described in terms of color strings stretched between the projectile and the target, which decay into new strings and subsequently hadronize to produce observed hadrons. Color strings may be viewed as small areas in the transverse plane filled with color field created by colliding partons. With growing energy and size of the colliding system, the number of strings grows, and they start to overlap, forming clusters, in the transverse plane very much similar to disks in two dimensional percolation theory. At a certain critical density a macroscopic cluster appears that marks the percolation phase transition. This is the Color String Percolation Model (CSPM) [13,14]. The interaction between strings occurs when they overlap and the general result, due to the SU(3) random summation of charges, is a reduction in multiplicity and an increase in the string tension hence increase in the average transverse momentum squared, p 2 t . We assume that a cluster of n strings that occupies an area of S n behaves as a single color source with a higher color field Q n corresponding to the vectorial sum of the color charges of each individual string Q 1 .
The resulting color field covers the area of the cluster. As Q n = n 1 Q 1 , and the individual string colors may be oriented in an arbitrary manner respective to each other, the average Q 1i Q 1j is zero, and Q 2 n = nQ 2 1 . Knowing the color charge Q n one can obtain the multiplicity μ and the mean transverse momentum squared p 2 t of the particles produced by a cluster of n strings [14] μ n = nS n S 1 μ 0 ; where μ 0 and p 2 t 1 are the mean multiplicity and p 2 t of particles produced from a single string with a transverse area S 1 = πr 2 0 . For strings just touching each other S n = nS 1 , and μ n = nμ 0 , p 2 t n = p 2 t 1 . When strings fully overlap, S n = S 1 and therefore μ n = √ nμ 0 and p 2 t n = √ n p 2 t 1 , so that the multiplicity is maximally suppressed and the p 2 t n is maximally enhanced. This implies a simple relation between the multiplicity and transverse momentum μ n p 2 t n = nμ 0 p 2 t 1 , which means conservation of the total transverse momentum produced.
In the thermodynamic limit, one obtains an analytic expression [13,14] where F (ξ) is the color suppression factor. With F (ξ) → 1 as ξ → 0 and F (ξ) → 0 as ξ → ∞, where ξ = N s S 1 S N is the percolation density parameter. Equation (1) can be written as μ n = nF (ξ )μ 0 and p 2 t n = p 2 t 1 /F (ξ ). The critical cluster which spans S N , appears for ξ c ≥ 1.2 [15]. It is worth noting that CSPM is a saturation model similar to the Color Glass Condensate (CGC), where p 2 t 1 /F (ξ ) plays the same role as the saturation momentum scale Q 2 s in the CGC model [16][17][18].

Experimental determination of the color suppression factor F (ξ)
The suppression factor is determined by comparing the pp and A + A transverse momentum spectra. To evaluate the initial value of ξ from data for Au + Au collisions, a parameterization of pp events at 200 GeV is used to compute the p t distribution [19][20][21] dN c /dp 2 where a is the normalization factor. p 0 and α are parameters used to fit the data. This parameterization also can be used for nucleus-nucleus collisions to take into account the interactions of the strings [14] dN c /dp 2 The color suppression factor F (ξ) is related to the percolation density parameter ξ .
In pp collisions F (ξ) ∼ 1 at these energies due to the low overlap probability.
In this way the STAR analysis of charged hadrons obtained the preliminary results for the percolation density parameter, ξ at RHIC for several collisions systems as a function of centrality [19]. Figure 1 shows a plot of F (ξ) as a function of charged particle multiplicity per unit transverse area dN c dη /S N for Au + Au collisions at 200 GeV for various centralities for the STAR data [20,21]. The error on F (ξ) is ∼3 %. F (ξ) decreases in going from peripheral to central collisions. The ξ value is obtained using Eq. (5), which increases with the increase in centrality. The fit to the Au + Au points has the functional form The STAR results for Au + Au collisions at √ s NN = 200 GeV can be used to estimate F(ξ ) values for Pb + Pb collisions at different centralities using the fit function given by Eq. (6) for Au + Au. Recently, the ALICE experiment at LHC published the charged-particle multiplicity density data as a function of centrality in Pb + Pb collisions at √ s NN = 2.76 TeV [22]. The ALICE data points are shown in Fig. 1. For central 0-5 % in Pb + Pb collisions ξ = 10.56 as compared to ξ = 2.88 for central Au + Au collisions at 200 GeV. For Au + Au central collisions we have found that the Bjorken energy density ε in the collision is proportional to ξ . To evaluate ε the charged pion multiplicity at mid rapidity and the Schwinger QED2 production time  [19]. The error is smaller than the size of the symbol. The line is fit to the STAR data. The solid blue squares are for Pb + Pb at 2.76 TeV were used [21,23]. Figure 2 shows a plot of energy density as a function of ξ . ε = 0.788ξ for the range 1.2 < ξ < 2.88. The extrapolated value of ε for central Pb + Pb collision at 2.76 TeV is 8.32 GeV/fm 3 as shown in Fig. 2.

Determination of the temperature
The connection between the measured ξ and the temperature T (ξ) involves the Schwinger mechanism (SM) for particle production. The Schwinger distribution for massless particles is expressed in terms of p 2 t [24, 25] dn/dp 2 t ∼ e −πp 2 where the average value of the string tension is x 2 . The tension of the macroscopic cluster fluctuates around its mean value because the chromo-electric field is not constant. The origin of the string fluctuation is related to the stochastic picture of the QCD vacuum. Since the average value of the color field strength must vanish, it can not be constant but changes randomly from point to point [26]. Such fluctuations lead to a Gaussian distribution of the string tension for the cluster, which transforms SM into the thermal distribution [26] dn/dp 2 with x 2 = π p 2 t 1 /F (ξ ). The temperature is expressed as [27] Recently, it has been suggested that fast thermalization in heavy ion collisions can occur through the existence of an event horizon caused by a rapid de-acceleration of the colliding nuclei [28]. The thermalization in this case is due to the Hawking-Unruh effect [29,30]. In CSPM the strong color field inside the large cluster produces de-acceleration of the primary qq pair which can be seen as a thermal temperature by means of the Hawking-Unruh effect. The string percolation density parameter ξ which characterizes the percolation clusters measures the initial temperature of the system. Since this cluster covers most of the interaction area, this temperature becomes a global temperature determined by the string density. In this way at ξ c = 1.2 the connectivity percolation transition at T (ξ c ) models the thermal deconfinement transition. We adopt the point of view that the experimentally determined universal chemical freeze-out temperature (T f ) is a good measure of the phase transition temperature, T c [31]. p 2 t 1 is evaluated using Eq. (9) at ξ c = 1.2 with T f = 167.7 ± 2.6 MeV [32]. This gives p 2 t 1 = 207.2 ± 3.3 MeV which is close to 200 MeV used in a previous calculation of the percolation transition temperature [27]. This calibrates the CSPM temperature scale. The dynamics of massless particle production has been studied in QED2 quantum electrodynamics. QED2 can be scaled from electrodynamics to quantum chromodynamics using the ratio of the coupling constants. Here the production time for a boson (gluon) is τ pro = 2.405 mc 2 [25]. This gives τ pro ∼ 1.13 fm for central Au + Au collisions at √ s NN = 200 GeV. The temperature obtained using Eq. (9) was ∼193.6 MeV for Au + Au collisions. For Pb + Pb collisions the temperature is ∼262.2 MeV for 0-5 % centrality, which is expected to be ∼35 % higher than the temperature from Au + Au collisions [21]. A recent summary of the results from Pb + Pb collisions at the LHC has mentioned that the initial temperature increases at least by 30 % as compared to the top RHIC energy [33]. Table 1 gives the CSPM values ξ , T , ε and η/s at T /T c = 0.88, 1, 1.16 and 1.57.
One way to verify the validity of extrapolation from RHIC to LHC energy is to compare the energy density expressed as ε/T 4 with the available lattice QCD results.  Table 1 The measured percolation density parameter ξ , temperature T , energy density ε and η/s for the meson gas [36]; the hadron to QGP transition; Au + Au at 200 GeV and Pb + Pb at 2.76 TeV (estimated). Au + Au is for 0-10 % and Pb + Pb is for 0-5 % central events tice QCD results are from HotQCD Collaboration [34]. It is observed that at LHC energy the CSPM results are in excellent agreement with the lattice QCD results. The lattice and CSPM results are available for T /T c < 2.

Shear viscosity
The relativistic kinetic theory relation for the shear viscosity over entropy density ratio, η/s is given by [9,35] η s T λ mfp 5 (10) where T is the temperature and λ mfp is the mean free path given by n is the number density of an ideal gas of quarks and gluons and σ tr the transport cross section for these constituents. After the cluster is formed it behaves like a free gas of constituents. Equation (10) can be applied to obtain the shear viscosity. In CSPM the number density is given by the effective number of sources per unit volume L is the longitudinal extension of the source, L = 1 fm [27]. The area occupied by the strings is related to ξ through the relation (1 − e −ξ )S N . Thus the effective no. of sources is given by the total area occupied by the strings divided by the effective area of the string S 1 F (ξ).
In general N sources is smaller than the number of single strings. N sources equals the number of strings N s in the limit of ξ = 0. The number density of sources from Eqs. (12) and (13) becomes In CSPM the transport cross section σ tr is the transverse area of the effective string S 1 F (ξ). Thus σ tr is directly pro-portional to F (ξ) and hence to 1 T 2 . The mean free path is given by For a large value of ξ the λ mfp reaches a constant value. η/s is obtained from ξ and the temperature Well below ξ c , as the temperature increases, the string density increases and the area is filled rapidly and λ mfp and η/s decrease sharply. Above ξ c , more than 2/3 of the area are already covered by strings, and therefore the area is not filling as fast and the relatively small decrease of λ mfp is compensated by the rising of temperature, resulting in a smooth increase of η/s. The behavior of η/s is dominated by the fractional area covered by strings. This is not surprising because η/s is the ability to transport momenta at large distances and that has to do with the density of voids in the matter. Figure 4 shows a plot of λ mfp , T and λ mfp × T as a function of ξ . Thus the product T (ξ) × λ mfp will have a minimum in η/s. It has been shown that η/s has a minimum at the critical point for various substances for example helium, nitrogen and water [10]. Thus the measurement of η/s as a function of temperature can indicate the critical point in the QCD phase diagram with T ∼ 175-185 MeV. Figure 5 shows a plot of η/s as a function of T /T c . The estimated value of η/s for Pb + Pb is also shown in Fig. 5 at T /T c = 1.57. The lower bound shown in Fig. 5 is given by AdS/CFT [11]. These results from STAR and ALICE data show that the η/s value is 2.5 and 3.3 times the KSS bound [11].

Results and discussion
The theoretical estimates of η/s has been obtained as a function of T /T c for both the weakly (wQGP) and strongly (sQGP) coupled QCD plasma are shown in Fig. 5 [9]. It is seen that at the RHIC top energy η/s is close to the sQGP. Even at the LHC energy it follows the trend of the sQGP. By extrapolating the η/s CSPM values to higher temperatures it is clear that η/s could approach the weak coupling limit near T /T c ∼ 5.8. The CSPM η/s value for the  [9]. The estimated value for Pb + Pb at 2.76 TeV for 0-5 % centrality is shown as a solid blue square. The red dotted line represents the extrapolation to higher temperatures from the CSPM. The hadron gas value for η/s ∼ 0.7 is shown as solid black circle at T /T c ∼ 0.88 [36] hadron gas is in agreement with the calculated value using measured elastic cross sections for a gas of pions and kaons [37]. η/s has also been obtained in several other calculations for pure glue matter [38], in the semi quark qluon plasma [39] and in quasiparticle description [40]. In pure SU(3) gluodynamics a conservative upper bound for η/s was obtained η/s = 0.134(33) at T = 1.65T c [41]. In the quasiparticle approach also low η/s ∼ 0.2 is obtained for T > 1.05T c and rises very slowly with the increase in temperature [42]. In CSPM also η/s grows with temperature as 0.16T /T c .
The CSPM model calculations have also successfully described the elliptic flow and the nuclear modification factor at RHIC and LHC energies [43]. In addition CSPM has determined the equation of state of the QGP and the bulk thermodynamic value of ε/T 4 and s/T 3 in excellent agreement with Lattice Gauge calculations [21]. This emphasizes the quantitative nature of the CSPM when applied to the data at ∼1 TeV scale.

Summary
In summary the relativistic kinetic theory relation for shear viscosity to entropy density ratio η/s = 1 5 T λ mfp was evaluated as a function of the temperature using the measured transverse momentum spectra and the Color String Percolation Model. The color suppression factor F (ξ) was extracted from the transverse momentum spectrum of charged hadrons. We found η/s = 0.204 ± 0.020 at T /T c = 1.15 (RHIC) and η/s = 0.260 ± 0.020 at T /T c = 1.57 (LHC). In the phase transition region η/s is 2-3 times the conjectured quantum limit for RHIC to LHC energies. The whole picture is consistent with the formation of a fluid with a low shear to viscosity ratio. The percolation framework provides us with a microscopic picture which predicts the early thermalization required for hydrodynamical calculations.
The minimum in η/s can be studied as a function of the beam energy at RHIC that could locate the critical point/crossover in the QCD phase diagram seen in substances like helium, nitrogen and water [6,10]. The accurate determination of η/s is also important for the evaluation of another transport coefficient, the jet quenching parameterq [44,45].