# Clustering of Color sources and the Equation of State in Heavy Ion Collisions at RHIC and LHC Energies

###### Abstract

The initial temperature , energy density , and formation time of the initial state of the QGP formed in the heavy ion collisions at RHIC and LHC energies are determined using the data driven Color String Percolation Model (CSPM). Multiparticle production by interacting strings stretched between projectile and target form a spanning cluster at the percolation threshold. The relativistic kinetic theory relation for is evaluated as a function of and the mean free path () using data and CSPM. (, ) describes the transition from a strongly interacting QGP at to a weakly coupled QGP at . We find that the reciprocal of is equal to the trace anomaly which also describes the transition. We couple this initial state of the QGP to a 1D Bjorken expansion to determine the sound velocity of the QGP for 0.85 . The bulk thermodynamic quantities and the equation of state are in excellent agreement with LQCD results.

## 1 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) [1, 2]. 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, . We assume that a cluster of strings that occupies an area of behaves as a single color source with a higher color field corresponding to the vectorial sum of the color charges of each individual string . The resulting color field covers the area of the cluster. As , and the individual string colors may be oriented in an arbitrary manner respective to each other , the average is zero, and .

Knowing the color charge one can obtain the multiplicity and the mean transverse momentum squared of the particles produced by a cluster of strings [2]

(1) |

where and are the mean multiplicity and of particles produced from a single string with a transverse area . This implies a simple relation between the multiplicity and transverse momentum , which means conservation of the total transverse momentum produced.

In the thermodynamic limit, one obtains an analytic expression [1, 2]

(2) |

where is the color suppression factor and is the percolation density parameter. Eq.(1) can be written as and . The critical cluster which spans , appears for 1.2 [3]. It is worth noting that CSPM is a saturation model similar to the Color Glass Condensate (CGC), where plays the same role as the saturation momentum scale in the CGC model [4, 5].

## 2 Experimental Determination of the Color Suppression Factor

The suppression factor is determined by comparing the and A+A transverse momentum spectra. To evaluate the value of from data for Au+Au collisions, a parameterization of events at 200 GeV is used to compute the distribution

(3) |

where a is the normalization factor. 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 [2]

(4) |

The color suppression factor is related to the percolation density parameter .

(5) |

## 3 Energy Density

The initial energy density above is given by [6]

(6) |

To evaluate we use the charged pion multiplicity at midrapidity and values from STAR for 0-10% central Au-Au collisions with 200 GeV [7]. We can calculate using the CSPM thermal distribution Eqs.(9) and (10). For , GeV, adding the extra energy required for the rest mass of pions at hadronization GeV. The error on is due to the error on [8].

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 [9]. The production time for a boson (gluon) is [10]

(7) |

For Au-Au collisions at 200 GeV Eqs. (6) and (7) gives = 2.270.16 GeV/ at =2.88. In CSPM the energy density is proportional to . From the measured value of and it is found that is proportional to for the range , GeV/ [7, 11] as shown in Fig.1.

## 4 Determination of the Temperature

The connection between the measured and the temperature involves the Schwinger mechanism (SM) for particle production. The Schwinger distribution for massless particles is expressed in terms of [10, 9]

(8) |

where the average value of the string tension is . 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 [12]. Such fluctuations lead to a Gaussian distribution of the string tension for the cluster, which transforms SM into the thermal distribution [12]

(9) |

with = . The temperature is expressed as [8, 13]

(10) |

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 [14]. The thermalization in this case is due to the Hawking-Unruh effect [15, 16]. In CSPM the strong color field inside the large cluster produces de-acceleration of the primary 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 = 1.2 the connectivity percolation transition at models the thermal deconfinement transition.

We adopt the point of view that the experimentally determined universal chemical freeze-out temperature () is a good measure of the phase transition temperature, . is evaluated using Eq.(10) at = 1.2 with = 167.7 2.6 MeV [17]. This gives = 207.2 3.3 MeV. This calibrates the CSPM temperature scale.

Figure 2 shows a plot of as a function of T/. The lattice QCD results are from HotQCD Collaboration [18]. 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/. Table I gives the CSPM values , , and at = 0.88, 1, 1.16 and 1.57.

## 5 The Shear Viscosity to Entropy Density ratio

The relativistic kinetic theory relation for the shear viscosity over entropy density ratio, is given by [19, 20]

(11) |

where T is the temperature and is the mean free path given by

(12) |

is the number density of an ideal gas of quarks and gluons and the transport cross section for these constituents.

System | T (MeV) | |||
---|---|---|---|---|

Meson Gas | 0.22 | 150.0 | - | 0.76 |

Hadron to QGP | 1.2 | 167.7 | 0.94 | 0.240 |

Au+Au | 2.88 | 193.6 | 2.27 | 0.204 |

Pb+Pb | 10.56 | 262.2 | 8.32 | 0.260 |

After the cluster is formed it behaves like a free gas of constituents. Eq. (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 [21]

(13) |

L is the longitudinal extension of the source, L = 1 [13]. The area occupied by the strings is related to through the relation . Thus the effective no. of sources is given by the total area occupied by the strings divided by the effective area of the string .

(14) |

In general is smaller than the number of single strings. The number density of sources from Eqs. (12) and (13) becomes

(15) |

In CSPM the transport cross section is the transverse area of the effective string . Thus is directly proportional to and hence to . The mean free path is given by

(16) |

For a large value of the reaches a constant value. is obtained from and the temperature

(17) |

The behavior of is dominated by the fractional area covered by strings for . This is not surprising because is the ability to transport momenta at large distances and that has to do with the density of voids in the matter.

Figure 3 shows a plot of as a function of T/. The estimated value of for Pb+Pb is shown in Fig. 3 at T/ = 1.57. The lower bound shown in Fig. 3 is given by AdS/CFT [23]. These results from STAR and ALICE data show that the value is 2.5 and 3.3 times the KSS bound [23].

The theoretical estimates of has been obtained as a function of T/ for both the weakly (wQGP) and strongly (sQGP) coupled QCD plasma are shown in Fig. 3 [19]. It is seen that at the RHIC top energy is close to the sQGP. Even at the LHC energy it follows the trend of the sQGP. By extrapolating the CSPM values to higher temperatures it is clear that could approach the weak coupling limit near 6. The CSPM value for the hadron gas is in agreement with the calculated value using measured elastic cross sections for a gas of pions and kaons [22].

## 6 The trace Anomaly and the reciprocal of

The and both describe the transition from a strongly coupled QGP to a weakly coupled QGP. We find that the reciprocal of is in quantitative agreement with the trace anomaly over this wide temperature range. This result is shown in Fig.4. The minimum in =0.20 at 1.15 determines the peak of the interaction measure 5 in agreement with the recent HotQCD values[24].

## 7 The Equation of State

We use CSPM coupled to a 1D Bjorken expansion. The input parameters are the initial temperature T, the initial energy density , and the trace anomaly are determined by data. The Bjorken 1D expansion gives the sound velocity

(18) | |||

(19) | |||

(20) |

where is the energy density, s the entropy density, the proper time, and the sound velocity. Since and one gets

(21) |

From above equations can be expressed in terms of

(22) |

## 8 Conclusions and open questions

Two central objectives in the experimental characterization of the QGP are the Equation Of State (EOS) and the shear viscosity to entropy ratio . We have found that is equal to the trace anomaly . We determine the bulk thermodynamics, , and the EOS of the QGP as a function of (CSPM). The results are in excellent agreement with LQCD numerical simulations as a function of (LQCD). CSPM predicts that the QGP will be formed in high multiplicity p-p collisions at LHC. The CSPM ratio can be used to test the quasi-particle hypothesis where the jet transport parameter is , which can be compared with measured in jet quenching using the the ALICE EMCAL detector at LHC.

## 9 Acknowledgments

This research was supported by the Office of Nuclear Physics within the U.S. Department of Energy Office of Science under Grant No. DE-FG02-88ER40412.

## References

## References

- [1] M. A. Braun and C. Pajares 2000 Eur. Phys. J. C 16 349.
- [2] M. A. Braun, F. del Moral and C. Pajares 2002 Phys. Rev. C 65 024907.
- [3] H. Satz 2000 Rep. Prog. Phys. 63 1511.
- [4] L. McLerran and R. Venugopalan 1994 Phys. Rev. D 49 2233.
- [5] J. Dias de Deus and C. Pajares 2011 Phys. Lett. B 695 455.
- [6] J. D. Bjorken 1983 Phys. Rev. D 27 140.
- [7] B. I. Abelev et al. (STAR Collaboration) 2009 Phys. Rev. C 79 34909.
- [8] R. P. Scharenberg, B. K. Srivastava and A. S. Hirsch 2011 Eur. Phys. J. C 71 1510.
- [9] C. Y. Wong, Introduction to high energy heavy ion collisions (World Scientific,1994).
- [10] J. Schwinger 1962 Phys. Rev. 128 2425.
- [11] B. K. Srivastava, R. P. Scharenberg, T. Tarnowsky, (STAR Collaboration) 2006 Nukleonika 51, s109.
- [12] A. Bialas 1999 Phys. Lett. B 466 301.
- [13] J. Dias de Deus and C. Pajares 2006 Phys. Lett. B 642 455.
- [14] D. Kharzeev, E. Levin and K. Tuchin 2007 Phys. Rev. C 75 044903.
- [15] S. W. Hawking, Commun. 1975 Math. Phys. 43 199.
- [16] W. G. Unruh 1976 Phys. Rev. D 14 870.
- [17] F. Becattini, P. Castorina, A. Milov and H. Satz 2010 Eur. Phys. J. C 66 377.
- [18] A. Bazavov et al. 2009 Phys. Rev. D 80, 014504.
- [19] T. Hirano and M. Gyulassy 2006 Nucl. Phys. A 769 71.
- [20] P. Danielewicz and M. Gyulassy 1985 Phys. Rev. D 31 53.
- [21] I. Bautista, J. Dias de Deus and C. Pajares arXiv:hep-ph/1102:3837.
- [22] M. Prakash et al.1993 Phys. Rept. 227 321.
- [23] P. K. Kovtun, D. T. Son and A. O. Starinets 2005 Phys. Rev. Lett. C 94 111601.
- [24] P. Petreczky, Lattice 2012, Cairns, Australia, 24 - 30 June 2012.