De-Confinement in small systems: Clustering of color sources in high multiplicity \bar{p}p collisions at \sqrt{s}= 1.8 TeV

De-Confinement in small systems: Clustering of color sources in high multiplicity p collisions at = 1.8 TeV

L. G. Gutay A. S. Hirsch C. Pajares R. P. Scharenberg B. K. Srivastava Department of Physics and Astronomy, Purdue University, West Lafayette, IN-47907, USA Departamento de Fisica de Particulas, Universidale de Santiago de Compos tela and Instituto Galego de Fisica de Atlas Enerxias(IGFAE), 15782 Santiago, de Compostela , Spain
July 14, 2019

It is shown that de-confinement can be achieved in high multiplicity non jet p collisions at = 1.8 TeV Fermi National Accelerator Laboratory(FNAL- E735) experiment. Previously the evidence for de-confinement was the demonstrated by the constant freeze out energy density in high multiplicity events. In this paper we use the same data but analyze the transverse momentum spectrum in the framework of the clustering of color sources. The charged particle pseudorapidities densities in the range 7.0 26.0 are considered. Results are presented for both thermodynamic and transport properties. The initial temperature and energy density are obtained and compared with the Lattice Quantum Chromo Dynamics(LQCD) simulations. The energy density () 11.5 for 25.0 is close to the value for 0-10% central events in Au+Au collisions at = 200 GeV. The shear viscosity to entropy density ratio() is 0.2 at the transition temperature. The result for the trace anomaly is in excellent agreement with LQCD simulations. These results confirm our earlier observation that the de-confined state of matter was created in high multiplicity events in p collisions at =1.8 TeV.

QGP, De-confinement, Shear viscosity, Trace anomaly

1 Introduction

The observation of high total multiplicity, high transverse energy, non-jet, isotropic events led Van Hove hove () and Bjorken boje () to conclude that high energy density events are produced in high energy p collisions larry (). These events have a far greater cross section than the jet production. In these events the transverse energy is proportional to the number of low transverse momentum particles. This basic correspondence can be explored over a wide range of the charged particle pseudorapidity density in p collisions at center of mass energy = 1.8 TeV. The analysis of charged particle transverse momentum data from the E735 experiment exhibits flow velocity of mesons and anti-baryons also indicating the possible evidence of QGP formation muller ().

Collective hydrodynamics flow has been successful in explaining the collisions at RHIC and LHC energies flow (). However collisions have been considered different from the heavy ion collisions. The observation of long range rapidity correlations, the so called “ridge”, similar to that seen in heavy ion collisions in high multiplicity collisions at = 7 TeV by the CMS experiment suggests the evidence of strong radial flow CMS (); shuryak (). There are several theoretical papers which support the view that QGP can be formed in high multiplicity events werner (); andres (); ghosh ().

In our earlier work, published in 2002, the evidence of hadronic de-confinement in p collisions at = 1.8 TeV was presented e735d (). Based on the HBT analysis a constant freeze-out energy density 1 for high multiplicity p events was measured e735d (). The freeze-out energy density was found to be independent of the .

The objective of this work is to further analyze the published E735 data on the transverse momentum spectra of charged particles in the framework of Color String Percolation Model(CSPM). CSPM has been successfully applied to heavy ion data for thermodynamics and transport coefficients eos (); cpod13 (); eos2 (); IS2013 (); eos3 (). After a brief description of the color string percolation model (CSPM)pajares1 (); pajares2 () and the E735 experiment the results are presented for the shear viscosity to entropy density ratio() and the Equation of State (EOS) of the deconfined matter.

2 String Interactions and Percolation

Multiparticle production at high energies is currently described in terms of color strings stretched between the projectile and target. Hadronizing 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. 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 isich (); satz1 (). As increases the fraction of covered by this spanning cluster increases.

We assume that a cluster of strings behaves as a single string with an energy-momentum that corresponds to the sum of energy-momenta of the individual strings and with a higher color field, corresponding to the vectorial sum of the color field of each individual string pajares1 (); pajares2 (). One can obtain the multiplicity and the mean transverse momentum squared of the particles produced by a cluster of strings pajares2 ()


where and are the mean multiplicity and of particles produced from a single string with a transverse area . In the limit of high string density, one obtains pajares1 (); pajares2 ()


where is the color suppression factor. is related to the .


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 were found to be in excellent agreement with experiment IS2013 ().

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 cgc (); perx ().

3 E735 experiment

Experiment E735 was run during the 1988-1989 Tevatron running period, primarily collecting data triggered to enrich high multiplicity events in collisions at = 1.8 TeV e735a (). The E735 was located at the interaction region of the Fermi National Accelerator Laboratory (FNAL) e735a (); e735b (); e735c (). The p interaction region was surrounded by a cylindrical drift chamber which in turn was covered by a single layer hodoscope including endcaps. This system measured the total charged particle multiplicity 10 200 in the pseudorapidity range 3.25. A sidearm magnetic spectrometer with tracking chambers and time of flight counters, provided particle identified momenta spectra in the range 0.1 1.5 GeV/c. The spectrometer covered −0.37 +1.00 with ( is the azimuthal angle around the beam direction).

The multiplicity dependence of the transverse momentum spectra was measured in p collisions at = 1.8 TeVe735c (). The invariant spectra was fitted with a power law ua1 (). Table 1 shows the power law fit parameters for high multiplicity events e735c ().

47 (minbias) 7.3 1.25 8.35
85 13.07 1.052 7.038
105 16.15 1.001 6.743
135 20.76 1.001 6.581
165 25.38 1.061 6.766
Table 1: Number of tracks as measured by the E735 experiment in the pseudorapidity range 3.25, , and the fit parameters and to the invariant distribution e735b (); e735c ().

4 Determination of the color suppression factor F()

The suppression factor is determined by comparing the charged particle spectra from low energy collisions and high multiplicity transverse momentum spectra from p. To evaluate the initial value of from data a parameterization of events at = 200 GeV is also used to compute the distribution


where a is the normalization factor. = 1.98 and = 12.88 are parameters used to fit the data. This parameterization is used for high multiplicity events in p collisions at = 1.8 TeV to take into account the interactions of the strings pajares2 ().


In pp collisions 1 at = 200 GeV due to the low overlap probability eos (). Figure (1) shows a plot of as a function of charged particle multiplicity per unit transverse area for high multiplicity events along with the minbias events. For collisions the overlap area is taken as the inelastic cross section levente (); cdf (); alice (). For p collision at =1.8 60 mb cdf (). The error on is . The transverse overlap area measured in E735 by Hanbury-Brown-Twiss is in excellent agreement with this inelastic cross section e735d (). The results from Au+Au collisions at = 200 GeV eos2 (); eos3 () are also shown in Fig. (1).

Figure 1: The color suppression factor as a function of . The solid green squares are from p at =1.8 TeV. The solid red circles are for Au+Au collisions at = 200 GeV(STAR data) eos2 (). The error is smaller than the size of the symbol. Lines are fit to p and the STAR data.

5 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 color 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 eos (); bialas (); pajares3 ()


with = . The temperature is expressed as


where is the average transverse momentum squared of particles produced from a single string. At the percolation transition = 1.2 and the critical value of the temperature is taken = 167 MeV from the statistical model, which is the universal chemical freeze-out temperature and is a good measure of the phase transition temperature fr1 (); fr2 (). This determines = 207 MeV eos (). The temperature obtained using Eq. (9) for 25 is 188 MeV. This temperature is closed to 193.6 MeV obtained for 0-10 %Au+Au collisions at = 200 GeVeos ().

T (MeV)
7.30(minbias) 0.778 160.31 0.50 0.30
13.07 1.39 170.73 0.86 0.23
16.15 1.42 171.22 1.07 0.23
20.76 1.84 178.06 1.39 0.21
25.38 2.30 185.07 1.75 0.21
Table 2: , the measured percolation density parameter , initial temperature , initial energy density and for p =1.8 TeV

6 Energy Density

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. After the initial temperature the CSPM perfect fluid may expand according to Bjorken boost invariant 1D hydrodynamics bjorken ()


where is the energy density, the transverse overlap area of the colliding nuclei, and the proper time. Above the critical temperature only massless particles are present in CSPM. To evaluate we use the charged pion multiplicity . The factor 3/2 in Eq.(12) accounts for the neutral pions. The average transverse mass is given by , where is the average transverse momentum of pion and the mass of pion. In Schwinger model is given by swinger ().

The pion multiplicity was obtained from the total multiplicity after subtracting the contributions due to kaons and protons e735a (); e735 (). The average for pions are flat in the region 7 26. Table 2 shows the values of and temperature for various values of .

Figure 2: as a function of T/. Solid green squares are from p collisions at = 1.8 TeV. The lattice QCD calculation is shown as dotted blue line hotqcd (). CSPM for Au+Au at = 200 GeV are red solid circles eos2 ().

With the determination of and one can compare the energy density expressed as with the available lattice QCD results hotqcd (). Figure 2 shows a plot of as a function of T/. The lattice QCD results are from the HotQCD Collaboration hotqcd (). The result for most central collisions for Au+Au at = 200 GeV is also shown in Fig.2. It is observed that for highest multiplicity 27, is close to the value obtained for Au+Au at 200 GeV. The minimum bias result is below the percolation threshold. This result confirms the formation of the QGP in high multiplicity events in p interactions at = 1.8 TeV.

7 Shear viscosity to entropy density ratio /s

In our earlier work the shear viscosity to entropy density ratio was obtained in the framework of kinetic theory and the string percolation eos2 (). The following expression was obtained for gul2 (); gul1 (); eos2 ().


where T is the temperature and L is the longitudinal extension of the source 1 fm.

Fig. 3 shows as a function of the temperature eos2 (). The lower bound shown in Fig. 3 is given by the AdS/CFT conjecture kss ().

The results from Au+Au at = 200 GeV collisions are also shown in Fig. 3 for comparison purposes. It is seen that for 25 the is equal to the value from the most central collisions in Au+Au at 200 GeV.

Figure 3: as a function of temperature T using Eq. (11) eos2 (). Solid green squares are from p collisions at = 1.8 TeV. CSPM results for Au+Au at = 200 GeV are shown as solid red circles eos2 (). The meson gas value for 0.7 is shown as solid black circle at T 150 MeV prakash (). The lower bound shown is given by the AdS/CFT kss ().

8 and Trace anomaly

The trace anomaly () is the expectation value of the trace of the energy-momentum tensor, , which measures the deviation from conformal behavior and thus identifies the interaction still present in the medium cheng (). We find that the reciprocal of is in quantitative agreement with from LQCD over a wide range of temperatures cpod13 (); IS2013 (). Fig. 4 shows the for p along with the CSPM calculation and LQCD simulations. The minimum in determines the peak of the interaction measure 5 in agreement with the recent HotQCD values lattice12 (). Figure 4 also shows the results from the Wuppertal Collaboration wuppe (). The maximum in corresponds to the minimum in . This is true for the highest multiplicity events in p collisions.

Figure 4: The trace anomaly vs temperature lattice12 (). Solid green squares are from p collisions at = 1.8 TeV. Blue open squares are from HotQCD Collaboration lattice12 (). Black stars are from Wuppertal Collaboration wuppe ().

9 Equation of State: Sound velocity

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 far the sound velocity can be written bjorken ()


In the above equation the entropy density is expressed as . The pressure is related to trace anomaly . We can express in terms of

Figure 5: The speed of sound from CSPM (red circles) and Lattice QCD-p4 (blue dash line) versus hotqcd (). The green solid squares are from p at = 1.8 TeV.

Since there is no direct way to obtain pressure in CSPM, we have made the assumption that the dimensionless ratio . Fig. 5 shows a plot of as a function of . It is observed that the CSPM results are in very good agreement with the lattice calculations hotqcd (). This suggests that the as obtained from the data can be represented by .

10 Summary

The p collisions at = 1.8 TeV from the E735 experiment have been analyzed in using the clustering of color sources phenomenology. The results for temperature, energy density and transport coefficient have been extracted for high multiplicity events. The trace anomaly is also obtained and compared with the LQCD results. The clustering of color sources has shown that the determination of as a function of temperature is an important quantity that relates to the trace anomaly . The main assumption of the present approach is that the inverse of represents the trace anomaly, . The clustering of color sources (percolation) provides us with a microscopic partonic picture that also connects the transport properties of the QGP to its thermodynamics.

These results strongly suggest that even in small systems at high energy and high multiplicity events QGP formation is possible as seen in collisions at = 1.8 TeV.

The analysis of the collisions at = 7 and 14 TeV at the CERN LHC can map event higher temperatures and energy densities. Then it will be possible to obtain the trace anomaly at higher temperature as shown in Fig.(4). The attenuation of high particles in high multiplicity events would be a second signal that the de-confined matter is indeed the QGP.

11 Acknowledgment

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. CP thanks Ministry of Economy of Spain and FEDER funds for financial support.


  • (1) L. V. Hove, Phys. Lett. B 118 (1982) 138.
  • (2) J. D. Bjorken, Fermilab Pub. 82/44-THY.
  • (3) L. McLerran, Rev. Mod. Phys. 58 (1986) 1021.
  • (4) P. Levai, B. Muller, Phys. Rev. Lett. 67 (1991) 1519.
  • (5) U. Heinz, R. Snellings, Ann. Rev. Nucl. Part. Sci 63 (2013) 123.
  • (6) V. Khachatryan et al. (CMS Collaboration), JHEP 1009 (2010) 091.
  • (7) T. Kalaydzhyan, E. Shuryak, arXiv:1503.0521.
  • (8) K. Werner, B. Guiot, I. Karpenko, T. Piergo, Nucl. Phys. A 931 (2014) 83.
  • (9) C. Andres, A. Moscoso, C. Pajares, Phys. Rev. C 90 (2014) 054902.
  • (10) P. Ghosh, S. Muhuri, J. K. Nayak, R. Varma, J. Phys. G 41 (2014) 03516.
  • (11) T. Alexopoulos et al. (E735 Collaboration), Phys. Lett. B 528 (2002) 43.
  • (12) R. P. Scharenberg, B. K. Srivastava, A. S. Hirsch, Eur. Phys. J. C 71 (2011) 1510.
  • (13) R. P. Scharenberg, PoS (CPOD 2013) 017.
  • (14) J. D. de Deus, A. S. Hirsch, C. Pajares, R. P. Scharenberg, B. K. Srivastava, Eur. Phys. J. C 72 (2012) 2123.
  • (15) B. K. Srivastava, Nucl. Phys. A 926 (2014) 142.
  • (16) M. A. Braun, J. D. de Deus, A. S. Hirsch, C. Pajares, R. P. Scharenberg, B. K. Srivastava, arXiv:1501.01524.
  • (17) M. A. Braun, C. Pajares, Eur. Phys. J. C 16 (2000) 349.
  • (18) M. A. Braun, F. del Moral, C. Pajares, Phys. Rev. C 65 (2002) 024907.
  • (19) J. Schwinger, Phys. Rev. 128 (1962) 2425.
  • (20) A. Bialas, Phys. Lett. B 466 (1999) 301.
  • (21) M. B. Isichenko, Rev. Mod. Phys. 64 (1992) 961.
  • (22) H. Satz, Rep. Prog. Phys. 63 (2000) 1511.
  • (23) L. McLerran, R. Venugopalan, Phys. Rev. D 49 (1994) 2233.
  • (24) J. D. de Deus, C. Pajares, Phys. Lett. B 695 (2011) 211.
  • (25) T. Alexopoulos et al. (E735 Collaboration), Phys. Rev. D 48 (1993) 984.
  • (26) T. Alexopoulos et al. (E735 Collaboration), Phys. Lett. B 336 (1994) 599.
  • (27) T. Alexopoulos et al. (E735 Collaboration), Phys. Lett. B 435 (1998) 453.
  • (28) C. Albajar et al. (UA1 Collaboration), Nucl. Phys. B 335 (1990) 261.
  • (29) B. I. Abelev et al. (STAR Collaboration), Phys. Rev. C 79 (2009) 34909.
  • (30) S. Klimenko, J. Konigsberg, T. M. Liss, Fermilab-FN-0741(2003).
  • (31) B. Abelev et al. (ALICE Collaboration), Eur. Phys. J. C 73 (2013) 2456.
  • (32) C. Y. Wong, Introduction to high energy heavy ion collisions, 1994.
  • (33) J. D. de Deus, C. Pajares, Phys. Lett. 642 (2006) 455.
  • (34) P. Braun-Munzinger, J. Stachel, C. Wetterich, Phys. Lett. B 596 (2004) 61.
  • (35) F. Becattini, P. Castorina, A. Milov, H. Satz, Eur. Phys. J. C 66 (2010) 377.
  • (36) J. D. Bjorken, Phys. Rev. D 27 (1983) 140.
  • (37) T. Alexopoulos et al. (E735 Collaboration), Phys. Rev. Lett. 64 (1990) 991.
  • (38) A. Bazavov et al., Phys. Rev. D 80 (2009) 014504.
  • (39) P. Danielewicz, M. Gyulassy, Phys. Rev. D 31 (1985) 53.
  • (40) T. Hirano, M. Gyulassy, Nucl. Phys. A 769 (2006) 71.
  • (41) P. K. Kovtun, D. T. Son, A. O. Starinets, Phys. Rev. Lett. 94 (2005) 111601.
  • (42) M. Prakash, M. Prakash, R. Venugopalan, G. Wleke, Phys. Rep. 227 (1993) 321.
  • (43) M. Cheng et al., Phys. Rev. D 80 (2010) 054504.
  • (44) P. Petreczky, Lattice 2012 Cairns, Australia (24 - 30 June).
  • (45) S. Borsanyi et al., Phys. Rev. D 80 (2009) 014504.
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