Hot Dense Matter : Deconfinement and Clustering of Color Sources in Nuclear Collisions

Within the first few microseconds from after the Big Bang, the hot dense matter was in the form of the Quark Gluon Plasm (QGP) consisting of free quarks and gluons. By colliding heavy nuclei at RHIC and LHC at a velocity close to the speed of light, we were able to create the primordial matter and observe the matter after expansion and cooling. In this report we present the thermodynamics and transport coefficients obtained in the framework of clustering of color sources in both hadron-hadron and nucleus-nucleus collisions at RHIC and LHC energies. Multiparticle production at high energies can be described in terms of color strings stretched between the projectile and target. At high string density single strings overlap and form color sources. This addition belongs to the non-perturbative domain of Quantum Chromo Dynamics (QGP) and manifests its most fundamental features. The Schwinger QED2 mechanism produces color neutral qq̄ pairs when color source strings break. Subsequent hardonization produces the observed hadrons. With growing energy and atomic number of the colliding nuclei the density of strings grows and more color sources form clusters in the transverse plane. At a certain critical density a macroscopic cluster appears, which marks the percolation phase transition. This is the Color String Percolation Model (CSPM). The critical density is identified as the deconfinement transition and happens at the hadronization temperature. The stochastic thermalization in pp and A-A is a consequence of the quantum tunneling through the event horizon introduced by the confining color fields, the Hawking-Unruh effect. The percolation approach within CSPM is successfully used to describe the crossover phase transition in the soft collision region. The same phenomenology when applied to both hadron-hadron and nucleus-nucleus collisions emphasizes the importance of color string density, creating a macroscopic cluster which identifies the connectivity required for a finite droplet of the QGP.


Introduction
What is the behavior of matter when we increase its density?The observed densities of our world expand over many orders of magnitude, from ∼10 −6 nucleons/cm 3 on average in our universe to ∼10 38 nucleons/cm 3 inside a nucleus and 10 39 nucleons/cm 3 in a neutron star.The study of the high density limit, specifically the confinement transition from hadrons to quarks and gluons can be regarded as the place where a high energy collision (two body) probes the short range limit and therefore meets the thermodynamics (many body) of this short distance dynamics.
In the case of Quantum-Chromo-Dynamics (QCD) we have confinement, which means the color field is confined to a small space region, forming hadrons which are colorless bound states of colored quarks.In the confining phase the quarks are dressed with gluons, acquiring an effective constituent mass up to 300 MeV (1/3 of the proton mass).For light quarks the constituent mass is much larger than the mass appearing in the QCD Lagrangian which is close to zero.In other words the light quark mass in the confined phase is generated by the confinement interaction and when deconfinement occurs the additional mass disappears.
The lattice QCD studies have shown, that at baryon chemical density µ B ∼ 0, color deconfinement and chiral symmetry restoration coincide.Hence, in a medium of low baryon density, the mass of the constituent quark vanishes at the deconfinement point T c and also the screening radius of the gluon cloud size vanishes.At low T and high µ, there is no reason to expect a similar behavior and probably there may be an intermediate region of massive dressed quarks between the hadronic phase and the deconfined and chiral symmetry restoration phase.At low T and high µ, other possibilities could exist such as quarkyonic matter and color superconductivity [1,2].A possible phase diagram is shown in Figure 1 [3].
The energy density of quark-gluon matter at a temperature T is given by (g q + g q)]T 4 (1) where g g , g q and g q are the degeneracy numbers of the gluons, quarks and antiquarks.For two flavors the total number of degrees of freedom (DOF) ∼37.
Quantitative predictions have not been evaluated using experimental data from high energy heavy ion collisions.In one of the recent article several big questions regarding the formation of QGP have been posed [4].The two most interesting problems are the determination of temperature and the thermalization in pp and AA collisions.The central challenge of this paper is to obtain the thermodynamics of the QGP from data and compare with the predictions of lattice QCD results.There are other review articles in the literature on the deconfinement both in pp and AA collisions [5][6][7].

High Energy Nuclear Collisions
A diagram of the time sequence of a heavy ion collision is shown in Figure 2. The time unit is f m/c = 3.3 × 10 −24 s.The initial energy density is developed in 1 f m/c.A new form of matter, the QGP, is composed of a dense volume of deconfined Q and Q pairs , which expands in ∼10 f m/c and subsequently forms the freeze out hadrons.To create a system of a large number (∼1000) of color strings, in a nuclear collision requires an accelerator with ∼100×10 9 eV (100 GeV) kinetic energy per nucleon with two colliding beams.The two accelerators used in these nuclear collision studies are the Relativistic Heavy Ion Collider (RHIC) located at Brookhaven National Laboratory on Long Island New York and the Large Hadron Collider (LHC) at the Central European Research Nuclear Laboratory (CERN), in Geneva, Switzerland.
The strong field of a color source can lift a quark-antiquark pair out of the vacuum.This is the Schwinger QED 2 mechanism [9].The interaction of an initial quark Q in one nucleon and Q in the colliding nucleon are connected by the strong color string that subsequently breaks into two strings creating a q and q deconfined pair at a time τ pro given by [9] with m t being the average transverse mass of the q q pair.More details of string fragmentation is given by Lund model which has tried to address the collectivity in pp collisions [10].

The Color String Percolation Model
The thermodynamic quantities are derived from RHIC and LHC data using the Color String Percolation Model (CSPM), as seen in Figure 3.In CSPM strings have a finite area in the transverse space in terms of the color field.The area is restricted to a finite disc of a given radius dictated by confinement.When the string with radii r 0 and string density ξ are uniformly distributed on the much larger transverse nuclear collision area S N they may overlap.ξ is given by N s S 1 /S N where the single disc area is S 1 and N s is the total number of strings.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 the 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.At a certain critical density ξ c = 1.2 a macroscopic percolation cluster appears that marks the percolation phase transition [11][12][13].This is the Color String Percolation Model (CSPM).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 Knowing the color charge Q n one can obtain the multiplicity µ n and the mean transverse momentum squared p 2 t n of the particles produced by a cluster of n strings [13] where µ 1 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µ 1 , p 2 t n = p 2 t 1 .When color strings fully overlap, S n = S 1 and therefore µ n = √ nµ 1 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 is conservation of the total transverse momentum.
In the thermodynamic limit, one obtains an analytic expression [12,13] 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(3) 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 [14].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 [15,16].Saturation results from the overcrowding in impact parameter of low x partons of a boosted hadrons or nucleus, leading to the appearance of a scale, Q 2 s .This is the basic idea of CGC.For example the particle density in CSPM is given by In CGC particle density is related to the coupling constant α s (Q 2 s ), In both case particle density increases with the number of participants N part [16].

Color Suppression Factor
In the present work we have extracted the color suppression factor F(ξ) in high multiplicity events in pp collisions using CMS data from the transverse momentum spectra of pions at √ s = 0.9, 2.76, 7 and 13 TeV [17,18].Figure 4 shows a transverse momentum spectra for two multiplicity cuts at √ s = 7 TeV.For comparison purpose the spectra from pp collisions √ s = 200 GeV is also shown.The spectra becomes harder for higher multiplicity cuts.This is due to the fact that high string density color sources are created in the higher multiplicity events.To evaluate the initial value of F(ξ) from data in high multiplicity events in pp collisions, a parameterization of the experimental data of p t distribution in low energy pp collisions √ s = 200 GeV was used [20].The charged particle spectrum is described by a power law [19] where a is the normalization factor, p 0 and α are fitting parameters with p 0 = 1.98 and α = 12.87 [20].This parameterization is used in high multiplicity events in pp collisions to take into account the interactions of the strings [19].The color suppression factor F(ξ) encodes the effect of the interaction among strings once they overlap.The parameter p 0 in Equation ( 7) is for independent strings and gets modified.
where F(ξ) mult pp is the multiplicity dependent color suppression factor.In pp collisions F(ξ) pp ∼ 1 at low energies due to the low overlap probability.The spectra were fitted using Equation (9) in the softer sector with p t in the range 0.12-1.0GeV/c.Equation ( 9) is similar to Tsallis distribution [21] which can be obtained in the clustering of color sources frame by doing the convolution of the distribution of the different formed clusters size with the thermal distributions of the decay of each cluster.In the thermodynamic limit the color suppression factor F(ξ) is related to the percolation density parameter ξ [19].
Figure 5 shows the extracted value of F(ξ) as a function of N tracks from CMS experiment for √ s = 0.9-13 TeV.It is observed that the F(ξ) has a lower value for high multiplicity events.To compare high energy pp with the heavy ions results, we need to normalized N tracks with the interaction area in pp collisions.The interaction area S ⊥ = πR 2  pp has been computed in the IP-Glasma model, where R pp is the interaction radius [22].This is based on an impact parameter description of pp collisions, combined with an underlying description of particle production based on the theory of Color Glass Condensate [22].The interaction radius R pp is approximately a linear function of the charged particle multiplicity.In the IP-Glasma model R pp is dependent on gluon multiplicity [22] R pp = f pp (dN g /dy) 1/3 , ( 11) where x = (dN g /dy) 1/3 .The gluon multiplicity dN g /dy is related to the number of tracks seen in the CMS experiment by where ∆η ∼ 4.8 units of pseudorapidity.The interaction cross section S ⊥ as a function of N tracks /∆η, using Equation (11), is shown in Figure 6.S ⊥ increases with the multiplicity and for very high multiplicities it is approximately constant.Figure 7 shows the extracted value of F(ξ) as a function of N tracks /∆η scaled by the interaction area S ⊥ from the CMS experiment for √ s = 0.9-13 TeV.N tracks is the total charged particle multiplicity in the region |η| < 2.4 with ∆η ∼ 4.8 units of pseudorapidity.Color Suppression Factor F(ξ) in pp, pp and Au-Au collisions vs. N tracks /∆η scaled by the transverse area S ⊥ .For pp and pp collisions S ⊥ is multiplicity dependent as obtained from IP-Glasma model [22].In case of Au-Au collisions the nuclear overlap area was obtained using the Glauber model [23].
The results from FNAL (Fermi National Accelerator Laboratory) E735 experiment on pp collisions at √ s = 1.8 TeV is also shown in Figure 7 [24].In the E735 experiment the total charged particle multiplicity was 10 < N tracks < 200 in the pseudorapidity range |η| <3.25 with ∆η ∼ 6.5 units of pseudorapidity.It is observed that the E735 results follow the trend seen in CMS data.The decrease in F(ξ) for high multiplicity events is due to the high string density created in theses events.

Temperature
The connection between F(ξ) 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 [9,25] dn/dp 2 t ∼ exp(−π p 2 t /x 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 cannot be constant but changes randomly from point to point [26,27].Such fluctuations lead to a Gaussian distribution of the string tension which gives rise to thermal distribution [26] with x 2 = π p 2 t 1 /F(ξ).The temperature is expressed as [20,28] We will adopt the point of view that the experimentally determined chemical freeze-out temperature is a good measure of the upper end of the cross over phase transition temperature, T c [29].The single string average transverse momentum p 2 t 1 is calculated at ξ c = 1.2 with the universal chemical freeze-out temperature of 167.7 ± 2.6 MeV [30].This gives p 2 t 1 = 207.2± 3.3 MeV which is close to 200 MeV used previously in the calculation of percolation transition temperature [28].
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 deceleration of the colliding nuclei [31].The thermalization in this case is due to the Hawking-Unruh effect [32,33].
It is well known that the black holes evaporates by quantum pair production and behave as if they have an effective temperature of where 1/4GM is the acceleration of gravity at the surface of a black hole of mass M. The rate of pair production in the gravitational background of the black hole can be evaluated by considering the tunneling through the event horizon.Unruh showed that a similar effect arises in a uniformly accelerated frame, where an observer detects the thermal radiation with the temperature T = a/2π, where a is the acceleration.Similarly, in hadronic interactions the probability to produce states of masses M due to the chromoelectric field E and color charge is given by the Schwinger mechanism [9] which is similar to the Boltzmann weight in a heat bath with an effective temperature In CSPM the strong color field inside the large cluster produces deceleration of the primary q q pair which can be seen as a thermal temperature by means of the Hawking-Unruh effect.This implies that the radiation temperature is determined by the transverse extension of the color flux tube/cluster in terms of the string tension [34][35][36].
The string percolation density parameter ξ which characterizes the percolation clusters measures the initial temperature of the system.Since the 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.The temperature obtained using Equation ( 16) was ∼193.6 MeV for Au-Au collisions at √ s NN = 200 GeV in reasonable agreement with T i = 221 ± 19 stat ± 19 sys MeV from the enhanced direct photon experiment measured by the PHENIX Collaboration [37].For Pb-Pb collisions at √ s NN = 2.76 TeV the temperature is ∼262.2MeV for 0-5% centrality, which is expected to be ∼35% higher than the temperature from Au-Au collisions [20].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 [38].The direct photon measurements from ALICE gives the temperature of T i = 304 ± 51 MeV [39].The agreement with the measured temperature shows that the temperature obtained using Equation ( 16) can be termed as the initial temperature of the percolation cluster.
Figure 8 shows a plot of temperature as a function of N tracks /∆η scaled by S ⊥ .Temperature from both hadron-hadron and nucleus-nucleus collisions fall on a universal curve when multiplicity is scaled by the transverse interaction area.The horizontal line at ∼165 MeV is the universal hadronization temperature obtained from the systematic comparison of the statistical thermal model parametrization of hadron abundances measured in high energy e + e − , pp and A-A collisions [30].In Figure 8 for √ s = 7 and 13 TeV higher multiplicity cuts show temperatures above the hadronization temperature and similar to those observed in Au-Au collisions at √ s NN = 200 GeV.The thermalization in pp and A-A collisions occurs through the existence of an event horizon caused by a rapid deceleration of the colliding nuclei [34,35].The thermalization is due to the Hawking-Unruh effect [32][33][34].In CSPM the strong color field inside the large cluster produces deceleration of the primary q q pair which can be seen as a thermal temperature.In case of A-A collisions it was observed that average transverse momentum is twice the temperature <p t > = 2T.This shows that the charged particle transverse momentum spectrum is exponentially distributed and the inverse slope parameter is the thermalized temperature [19].
After the initial temperature T > T c the CSPM perfect fluid may expand according to Bjorken boost invariant 1D hydrodynamics [41] where ε is the energy density, S n nuclear overlap area, and τ pro the production time for a boson (gluon) [25].The τ pro is given by Equation ( 2) ∼1 f m/c.The dynamics of massless particle production has been studied in QE2 quantum electrodynamics.QE2 can be scaled from electrodynamics to quantum chromodynamics using the ratio of the coupling constants [25].Figure 9 shows a plot of ε/T 4 as a function of T/T c that is in excellent agreement with numerical lattice gauge calculations for 0.8 < T/T c < 2. At T/T c ∼ 1.2 (T ∼ 195 MeV) we find that the QGP has ∼37 quark gluon degrees of freedom.
It is important to ask the question if matter created in high multiplicity pp events is nearly a perfect fluid with a low shear viscosity as observed in heavy ion collisions [19].

Shear Viscosity and Trace Anomaly
The shear viscosity to entropy density ratio η/s was obtained in the framework of relativistic kinetic theory and the string percolation [42].In this picture the relevant parameter is the string density ξ.The relativistic kinetic theory relation for the shear viscosity over entropy density ratio, η/s is given by [43,44] η s Tλ m f p where T is the temperature and λ m f p 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 ( 22) can be applied to obtain the shear viscosity.In CSPM the number density is given by the effective number of sources per unit volume where L is the longitudinal extension of the source, L = 1.0-1.1 fm [28].η/s is obtained from Equations ( 22) and (24) as Well below ξ c , as the temperature increases, the string density increases and the area is filled rapidly and λ m f p 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 λ m f p 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 10a shows η/s as a function of the temperature.The lower bound shown in Figure 10a is given by the AdS/CFT conjecture [45].The results from pp collisions from √ s = 13 TeV shows a very small η/s and that is 2.7 times the AdS/CFT conjectured lower bound 1/4π.The theoretical estimates of η/s has been obtained as a function of T/T c for both the weakly interacting (wQGP) and strongly interacting (sQGP) coupled QCD plasma are shown in Figure 5 [43].√ s = 0.9, 2.76, 7 and 13 TeV.The lower bound shown is given by the AdS/CFT [45].For comparison purposes the results from Au-Au and Pb-Pb at √ s NN = 200 GeV and 2.76 TeV respectively are also shown as solid squares for 0-5 % centrality [19].η/s estimates for wQGP and sQGP are shown as dotted blue and green lines respectively [43].(b) The trace anomaly ∆ = (ε − 3p)/T 4 vs. temperature.Blue open squares are from HotQCD Collaboration [46].Black stars are from Wuppertal Collaboration [47].The CSPM results are obtained as ∆ = 1/(η/s) [19].The black dashed line both in (a,b) corresponds to extrapolation from CSPM at higher temperatures.
η/s has also been obtained in several other calculations for pure glue matter [48], in the semi quark gluon plasma [49] and in quasi particle description [50].In pure SU(3) gluodynamics a conservative upper bound for η/s was obtained η/s = 0.134 (33) at T = 1.65T c [51].In the quasi particle approach also low η/s ∼ 0.2 is obtained for T > 1.05T c and rises very slowly with the increase in temperature [52].In CSPM also η/s grows with temperature as 0.16T/T c .
In CSPM η/s is not needed to reproduce the elliptic flow and higher harmonics.The cluster formed by the strings has generally an asymmetric form in the transverse plane and acquires dimensions comparable to the nuclear overlap.This azimuthal asymmetry is at the origin of the elliptic flow in CSPM.The partons emitted at some point inside the cluster have to pass through the strong color field before appearing on the surface.The results of the simulation for different harmonics are in reasonable agreement with the experimental data on the p t and centrality dependencies [53,54].
The trace anomaly (∆) is the expectation value of the trace of the energy-momentum tensor, Θ µ µ = (ε − 3p), which measures the deviation from conformal behavior and thus identifies the interaction still present in the medium [55].We consider the ansatz that inverse of η/s is equal to trace anomaly ∆. η/s is in quantitative agreement with (ε − 3p)/T 4 over a wide range of temperatures [56,57].Figure 10b shows 1/(η/s) and the dimensionless quantity, (ε − 3p)/T 4 , obtained from lattice QCD simulations [46].The minimum in η/s ∼ 0.22 determines the peak of the interaction measure ∼4.5 in agreement with the recent HotQCD values [46].This happens at the temperature of T c ∼ 195 MeV. Figure 10b also shows the results from Wuppertal Collaboration [47].The maximum in ∆ corresponds to the minimum in η/s.Both ∆ and η/s describe the transition from a strongly coupled QGP to a weakly coupled QGP.
We are not aware of any theoretical work which directly relates the trace anomaly with the shear viscosity to entropy density ratio.However, the bulk viscosity ζ is related to both ∆ and η [58].A detailed study based on low energy theorems and the lattice result for ∆ shows that ζ/s rises very fast close to the critical temperature in such a way that its value at temperatures higher than T > 1.1T c is quite negligible [59].It was observed that ζ scales as α 4 s η where α s is the coupling constant.The trace anomaly ∆ is proportional to α 2 s [58].There are many other works in which ∆ and η/s have been obtained separately [48,50,60,61].
In our earlier work F(ξ) was obtained in Au-Au collisions at √ s NN = 200 GeV for various centralities using STAR data [19].The results are shown in Figure 7 along with pp (CMS) and pp (E735) collisions.It is observed that F(ξ) as a function of dN c /∆η scaled by the transverse interaction area falls on a universal scaling curve for both hadron-hadron and nucleus-nucleus collisions.F(ξ) values in high multiplicity events in pp collisions at √ s = 13 TeV are similar to those obtained in most central events in Au-Au collisions at √ s NN = 200 GeV.In Figure 8 the temperatures from Au-Au collisions at √ s NN = 200 GeV are compared with the pp collisions.In high multiplicity events for pp collisions at √ s = 13 TeV the temperature is the same as in most central events in Au-Au.The η/s as a function of temperature from Au-Au at √ s NN = 200 GeV and Pb-Pb at √ s NN = 2.76 TeV collisions is shown in Figure 10a along with pp collisions.The η/s value is similar both in the high multiplicity pp collisions at √ s = 13 TeV and in most central collisions in Au-Au at √ s NN = 200 GeV.This shows that the QGP created in high multiplicity pp collisions is also strongly coupled.
The similarity between high multiplicity pp events and Au-Au suggests that the thermalization in both systems is reached through the stochastic process (Hawking-Unruh) instead of kinetic process [62].

Summary and Conclusions
In the initial collision of pp and A-A nuclei at high energy, the Schwinger QED 2 mechanism scaled from electric to color fields, lifts color neutral q q pairs from the vacuum at the initial time t i = 2.405 h/<m t > ∼1 f m/c.Here <m t > is the average transverse mass of quarks when a color source string breaks.Color sources are formed by the overlapping of n individual color strings with transverse area S 1 .The overall result of color string breaking is the increase of the transverse momentum p t of the q q pairs and the decrease of the multiplicity expected from n individual strings.In the transverse plane the total color string area is reduced by F(ξ) due to the overlapping of individual strings.F(ξ) is determined by comparing the charged particle transverse momentum spectra from a lower energy pp collision (single string) with the various high energy charged particle transverse momentum spectra from pp and A-A collisions(overlapping strings).In the overlap area the strings start interacting to form clusters in the transverse plane very much similar to discs in 2D percolation theory.When the 2D percolation density parameter ξ reaches ξ =1.2 , the spanning cluster is formed and represents a communicating QGP droplet.
The initial temperature is given by Equation ( 16) and is due to the Hawking-Unruh barrier penetration of the event horizon provided by color confinement.The similarity between high multiplicity pp events and Au-Au indicates that the thermalization in both systems reached is due to a stochastic process rather than a kinetic process.We find ∼37 degrees of freedom at highest temperature of 195 MeV reached in Au-Au collisions at √ s NN = 200 GeV.
The percolation clustering begins ∼145 MeV and ends at the hadronization ie.critical temperature (T = 167.8MeV).The center of this crossover phase transition is in agreement with the lattice gauge pseudo-critical temperature T = 155 ± 9 MeV.The temperature dependence of the η/s ratio across the crossover region has a clear minimum value at T = 195 MeV and η/s = 0.22, which is ∼2.7 times the Ads/CFT value of 1/4π.Here the QGP has a low shear viscosity with strong residual interactions.We found that the inverse η/s of closely tracks the lattice QCD values of the trace anomaly ∆ = (ε − 3p)/T 4 for T from 150 MeV to 260 MeV.
The color string percolation model (CSPM) gives a successful description of the experimental data in the soft region.The clustering of color sources is not derived directly from QCD, but has its most fundamental features.The non-Abelian character is reflected in the coherence sum of the color fields.Confinement is responsible for the event horizon and the small transverse size of the strings.The interpretation of data using color string percolation describes the equation of state of the QGP close to and above the critical (hadronization) temperature.

Figure 2 .
Figure 2. Illustration of the dynamical evolution of relativistic heavy-ion collisions and the QGP expansion [8].

Figure 4 .
Figure 4. Transverse momentum spectra of pions from CMS experiment at √ s = 7 TeV for two different multiplicity cuts N track = 131 (red solid circle) and N track = 75 (green solid circle)[17].For comparison purpose the p t spectra from pp at √ s = 200 GeV is also shown (solid blue circle)[19].

Figure 8 .
Figure 8. Temperature vs. N tracks /∆η scaled by S ⊥ from pp and Au-Au collisions.The horizontal line at ∼165 MeV is the universal hadronization temperature [30,40].

Figure 10 .
Figure 10.(a) η/s as a function of temperature T using Equation (25) for √ s = 0.9, 2.76, 7 and 13 TeV.The lower bound shown is given by the AdS/CFT [45].For comparison purposes the results from Au-Au and Pb-Pb at √ s NN = 200 GeV and 2.76 TeV respectively are also shown as solid squares for 0-5 % centrality [19].η/s estimates for wQGP and sQGP are shown as dotted blue and green lines respectively [43].(b) The trace anomaly ∆ = (ε − 3p)/T 4 vs. temperature.Blue open squares are from HotQCD Collaboration[46].Black stars are from Wuppertal Collaboration[47].The CSPM results are obtained as ∆ = 1/(η/s)[19].The black dashed line both in (a,b) corresponds to extrapolation from CSPM at higher temperatures.