Clustering of Color Sources and the Equation of State of the QGP

Clustering of Color Sources and the Equation of State of the QGP

Brijesh K Srivastava Department of Physics, Purdue University, West Lafayette, IN-47906, USA

Possible phase transition of strongly interacting matter from hadron to a quark-gluon plasma state have in the past received considerable interest. It has been suggested that this problem might be treated by percolation theory. The clustering of color sources with percolation (CSPM) is used to determine the equation of state (EOS) and the transport coefficient of the Quark-Gluon Plasma (QGP) produced in central A-A collisions at RHIC and LHC energies.


-titleInternational Conference on New Frontiers in Physics(ICFP) 2012

1 Introduction

One of the main goal of the study of relativistic heavy ion collisions is to study the deconfined matter, known as Quark-Gluon Plasma (QGP), which is expected to form at large densities. It has been suggested that the transition from hadronic to QGP state can be treated by the percolation theory celik (). The formulation of percolation problem is concerned with elementary geometrical objects placed on a random d-dimensional lattice. The objects have a well defined connectivity radius , and two objects can communicate if the distance between them is less than . Several objects can form a cluster of communication. At certain density of the objects a infinite cluster appears which spans the entire system. This is defined by the dimensionless percolation density parameter isich (). Percolation theory has been applied to several areas ranging from clustering in spin system to the formation of galaxies. Figure 1 shows the transition from disconnected to connected system at high densities.

Figure 1: Left: Disconnected system. Right: Connected system

In nuclear collisions there is indeed, as a function of parton density, a sudden onset of large scale color connection. There is a critical density at which the elemental objects form one large cluster, loosing their independent existence. Percolation would correspond to the onset of color deconfinement and it may be a prerequisite for any subsequent QGP formation.

The determination of the EOS of hot, strongly interacting matter is one of the main challenges of strong interaction physics. Recent lattice QCD (LQCD) calculations for the bulk thermodynamic observables, e.g. pressure, energy density, entropy density and for the sound velocity have been reported bazavov (). Recently, attention has been focused on the shear viscosity to entropy density ratio as a measure of the fluidity teaney (); teaney1 (); lacey (); rom (). The observed temperature averaged , based on viscous hydrodynamics analyses of RHIC data, are suggestive of a strongly coupled plasma gul1 (); larry (). In this talk a percolation model coupled with thermodynamical relations has been utilized to obtain EOS and the from the experimental data.

2 String Interactions and Percolation

Multiparticle production at high energies is currently described in terms of color strings stretched between the projectile and target. Hardonizing these strings produce the observed hadrons. The strings act as color sources of particles through the creation of pairs from the sea. At low energies only valence quarks of nucleons form strings that then hadronize. The number of strings grows with the energy and with the number of nucleons of participating nuclei. Color strings may be viewed as small discs in the transverse space filled with the color field created by colliding partons. Particles are produced by the Schwinger mechanisms swinger (). With growing energy and size of the colliding nuclei the number of strings grow and start to overlap to form clusters pajares1 (); pajares2 (). At a critical density a macroscopic cluster appears that marks the percolation phase transition. 2D percolation is a non-thermal second order phase transition, but in CSPM the Schwinger barrier penetration mechanism for particle production and the fluctuations in the associated string tension due to the strong string interactions make it possible to define a temperature. Consequently the particle spectrum is ”born” with a thermal distribution bialas (). With an increasing number of strings there is a progression from isolated individual strings to clusters and then to a large cluster which suddenly spans the area. In two dimensional percolation theory the relevant quantity is the dimensionless percolation density parameter given by pajares1 (); pajares2 ()


where N is the number of strings formed in the collisions and is the transverse area of a single string and is the transverse nuclear overlap area. The critical cluster which spans , appears for 1.2 satz1 (). As increases the fraction of covered by this spanning cluster increases.

The percolation theory governs the geometrical pattern of string clustering. It requires some dynamics to describe the interaction of several overlapping strings. It is assumed that a cluster behaves as a single string with a higher color field corresponding to the vectorial sum of the color charge of each individual string. Knowing the color charge, one can calculate the multiplicity and the mean transverse momentum squared of the particles produced by a cluster of strings. One finds pajares1 (); pajares2 ()


where and are the mean multiplicity and average transverse momentum squared of particles produced by a single string . In the saturation limit, all the strings overlap into a single cluster that approximately occupies the whole interaction area, one gets the following universal scaling law


In the limit of high density one obtains




being the color suppression factor. It follows that


A similar scaling is found in the Color Glass Condensate approach (CGC)cgc (); perx (). The saturation scale in CGC corresponds to in CSPM. The net effect due to is the reduction in hadron multiplicity and increase in the average transverse momentum of particles. The CSPM model calculation for hadron multiplicities and momentum spectra was found to be in excellent agreement with experiment diasde (); diasde2 ().

3 Color Suppression Factor

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


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 pajares2 ()


In pp collisions 1 at these energies due to the low overlap probability. Figure 2 shows a plot of as a function of charged particle multiplicity per unit transverse area for Au+Au collisions at 200 GeV for various centralities for the STAR data nucleo (); levente (); eos (). decreases in going from peripheral to central collisions. The value is obtained using Eq. (6), which increases with the increase in centrality. These experimental values of are used to obtain temperature, EOS and .

Figure 2: Color suppression factor as a function of . The solid red circles are for Au+Au collisions at 200 GeV(STAR data) nucleo (). 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.

4 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 swinger (); wong ()


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 bialas (). Such fluctuations lead to a Gaussian distribution of the string tension for the cluster, which gives rise to the thermal distribution bialas (); eos ()


with = . The temperature is expressed as


5 Equation of State

Among the most important and fundamental problems in finite-temperature QCD are the calculation of the bulk properties of hot QCD matter and characterization of the nature of the QCD phase transition. The QGP according to CSPM is born in local thermal equilibrium because the temperature is determined at the string level. We use CSPM and thermodynamical relations to calculate energy density, entropy density and sound velocity (EOS). After the initial temperature the CSPM perfect fluid may expand according to Bjorken boost invariant 1D hydrodynamics bjorken ()


where is the energy density, s the entropy density, the proper time, and the sound velocity. Above the critical temperature only massless particles are present in CSPM. The initial energy density above is given by bjorken ()


To evaluate we use the charged pion multiplicity at midrapidity and values from STAR for 0-10% central Au+Au collisions at 200 GeV levente (). The factor 3/2 in Eq.(17) accounts for the neutral pions. The average transverse mass is given by , where is the transverse momentum of pion and being the mass of pion. The dynamics of massless particle production has been studied in two-dimensional quantum electrodynamics (QED2). QED2 can be scaled from electrodynamics to quantum chromodynamics using the ratio of the coupling constants wong (). The production time for a boson (gluon) is swinger ().


From the measured value of and it is found that is proportional to for the range . Figure 3 shows a plot of energy density as a function of . GeV/ nucleo (); levente ().

Figure 3: Energy density as a function of the percolation density parameter . The value for LHC energy is shown as blue square.

Figure 4 shows a plot of as a function of T/. The lattice QCD results are from HotQCD Collaboration bazavov (). CSPM is in excellent agreement with the LQCD calculations in the phase transition region for 1.5. The sound velocity and entropy density results are also in good agreement with LQCD calculations eos ().

Figure 4: as a function of T/.The lattice QCD calculation is shown as dotted blue line bazavov ().

6 Shear Viscosity

The relativistic kinetic theory relation for the shear viscosity over entropy density ratio, is given by gul1 (); gul2 ()


where T is the temperature and is the mean free path . is the number density of an ideal gas of quarks and gluons and the transport cross section for these constituents. 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 . The effective no. of sources is given by the total area occupied by the strings divided by the effective area of the string eos2 ().


is obtained from and the temperature


Figure 5 shows a plot of as a function of T/. The lower bound shown in Fig. 5 is given by AdS/CFT kss (). The theoretical estimates of obtained for both the weakly (wQGP) and strongly (sQGP) coupled QCD plasma are shown in Fig. 5 gul1 (). 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 5.8.

Figure 5: as a function of T/. Au+Au at 200 GeV for 0-10 centrality is shown as solid black square. wQGP and sQGP values are shown as dotted blue and green lines respectively gul1 (). 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 0.7 is shown as solid black circle at T/0.88 meson ().


The STAR results for Au+Au collisions at = 200 GeV have been extrapolated to estimate F values for Pb+Pb collisions at =2.76 TeV. The STAR points in Fig.2 has the functional form


Recently, the ALICE experiment at LHC published the charged-particle multiplicity density data as a function of centrality in Pb+Pb collisions at = 2.76 TeV alice1 (). The ALICE data points are shown in Fig.2. For central 0-5 in Pb+Pb collisions = 10.56 as compared to = 2.88 for central Au+Au collisions at 200 GeV. The extrapolated value of for central Pb+Pb collision at 2.76 TeV is 8.32 as shown in Fig.3. For Pb+Pb collisions the temperature is 262.2 MeV for 0-5 centrality, which is 35 higher than the temperature from Au+Au collisions eos ().

One way to verify the validity of extrapolation from RHIC to LHC energy is to compare the energy density expressed as with the available lattice QCD results. The LHC point is shown in Fig. 4. 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/.

The estimated value of for Pb+Pb is shown in Fig. 5 at T/ = 1.57. These results from STAR and ALICE data show that the value is 2.5 and 3.3 times the KSS bound kss ().

8 Summary

Two central objectives in the experimental characterization of the QGP are the Equation Of State (EOS) and the shear viscosity to entropy ratio . The percolation analysis of the color sources applied to STAR data at RHIC provides a compelling argument that the QGP is formed in central Au+Au collisions at 200 GeV. It also suggests that the QGP is produced in all soft high energy high multiplicity collisions when the string density exceeds the percolation transition. We found = 0.204 at = 1.15 ( RHIC ) and =0.260 at = 1.57 (LHC). In the phase transition region 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. Thus Clustering and percolation can provide a conceptual basis for the QCD phase diagram which is more general than the symmetry breaking satzx ().

9 Acknowledgement

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.


  • (1) T. Celik, F. Karsch, H. Satz, Phys. Lett. B97, (1980) 128.
  • (2) M. B. Isichenko, Rev. Mod. Phys.64, (1992) 961.
  • (3) A. Bazavov et al., Phys. Rev. D80, (2009) 014504.
  • (4) D. Teaney, Phys. Rev. C68, (2003) 034913.
  • (5) T. Schafer and D. Teaney, Rep. Prog. Phys. C72, (2009) 126001.
  • (6) R. A. Lacey et al., Phys. Rev. Lett. 98, (2007) 092301.
  • (7) P. Romatschke and U. Romatschke, Phys. Rev. Lett. 99, (2007) 272301.
  • (8) T. Hirano and M. Gyulassy, Nucl. Phys. A769, (2006) 71.
  • (9) L. P. Csernai, J. I. Kapusta and L. D. McLerran Phys. Rev. Lett. 97, (2006) 152303.
  • (10) J. Schwinger, Phys. Rev. 128, (1962) 2425.
  • (11) M. A. Braun and C. Pajares, Eur. Phys. J. C16, (2000) 349.
  • (12) M. A. Braun, F. del Moral and C. Pajares, Phys. Rev. C65, (2002) 024907.
  • (13) A. Bialas, Phys. Lett. B466, (1999) 301.
  • (14) H. Satz, Rep. Prog. Phys. 63, (2000) 1511.
  • (15) L. McLerran and R. Venugopalan, Phys. Rev. D49, (1994) 2233;49, (1994) 3352.
  • (16) J. Dias de Deus and C. Pajares, Phys. Lett. B695, (2011) 455.
  • (17) J. Dias de Deus, E. G. Ferreiro, C. Pajares, R. Ugoccioni, Eur. Phys. J. C40, (2005) 229.
  • (18) P. Brogueira, J. Dias de Deus, J. G. Milhano, Nucl. Phys. A. 832, (2010) 76.
  • (19) B. K. Srivastava, R. P. Scharenberg and T. Tarnowsky, (STAR Collaboration), Nukleonika 51, (2006) s109.
  • (20) B. I. Abelev et al., (STAR Collaboration), Phys. Rev. C79, (2009) 34909.
  • (21) R. P. Scharenberg, B. K. Srivastava and A. S. Hirsch, Eur. Phys. J. C71, (2011) 1510.
  • (22) J. D. Bjorken, Phys. Rev. D27, (1983) 140.
  • (23) C Y. Wong, Introduction to high energy heavy ion collisions, 289 (1994).
  • (24) P. Danielewicz and M. Gyulassy, Phys. Rev. D31, (1985) 53.
  • (25) J. Dias de Deus et al., Eur. Phys. J. C72, (2012) 2123.
  • (26) D. Teaney in QGP4, page 207, Eds. R. C. Hwa and X. N. Wang, World scientific (2010).
  • (27) P. K. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. C94, (2005) 111601.
  • (28) K. Aamodt et al., (ALICE Collaboration), Phys. Rev. Lett. 106, (2011) 032301.
  • (29) H. Satz, Quantum fiels theory in extreme environments, Paris, April 2009.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description