Net-proton measurements at RHIC and the QCD phase diagram

Net-proton measurements at RHIC and the QCD phase diagram


Two measurements related to the proton and anti-proton production near midrapidity in = 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV Au+Au collisions using the STAR detector at the Relativistic Heavy Ion Collider (RHIC) are discussed. At intermediate impact parameters the slope parameter of the directed flow versus rapidity () for the net-protons shows a non-monotonic variation as a function of the beam energy. This non-monotonic variation is characterized by the presence of a minimum in between = 11.5 and 19.6 GeV and a change in the sign of twice between = 7.7 and 39 GeV. At small impact parameters the product of the moments of net-proton distribution, kurtosis variance () and skewness standard deviation () are observed to be significantly below the corresponding measurements at large impact parameter collisions for = 19.6 and 27 GeV. The and values at these beam energies deviate from the expectations from Poisson statistics and that from a Hadron Resonance Gas model. Both these measurements have implications towards the understanding of the QCD phase structures, the first order phase transition and the critical point in the high baryonic chemical potential region of the phase diagram.

QCD phase diagram, Critical Point, First Order Phase Transition, Heavy-Ion Collisions, Quark Gluon Plasma, Net-proton

ain]Bedangadas Mohanty

1 Introduction

The formation of a hot and dense medium of deconfined quarks and gluons (QGP) has been established in high energy heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) facility at Brookhaven National Laboratory and the Large Hadron Collider (LHC) facility at CERN [1]. The transition from QGP to a hadron gas has been shown to be a cross over [2]. The focus of research in this field has now shifted towards two aspects, (a) characterizing the transport properties of QGP and (b) establishing the QCD phase structures at high baryonic chemical potential () region of the QCD phase diagram.

A rigorous phenomenological analysis of the precision data from relativistic heavy-ion collisions and theoretical advances over 14 years has led to quantitative estimates for some of the transport properties of a strongly interacting de-confined state of quarks and gluons. The estimated shear viscosity to entropy density ratio () is found to reflect the inviscid liquid property of the QGP and has a value between (1-2)/4 [3, 4]. That reflecting the stopping power or the opacity of QGP has been estimated by obtaining the square of the momentum transferred by the parton to the QGP per unity length () and is found to lie between 2-10 GeV/fm [5, 6].

Figure 1: Conjectured QCD phase diagram [7] .

On the other hand a dedicated program called the Beam Energy Scan (BES) to establish the phase diagram of QCD was launched at RHIC in the year 2010 to unravel the QCD phase structure at large . A range of from 20 MeV to 400 MeV of the phase diagram was covered by varying the from 200 to 7.7 GeV. The rich phase structure in the high region can be seen from the conjectured QCD phase diagram shown in Fig. 1 [7]. The two distinct features of the phase diagram at large are the first order phase boundary and the critical point (CP). In this paper we concentrate on the status of the experimental search for these two phase structures through the measurement of proton and antiproton production in heavy-ion collisions. Specifically we discuss the observable related to the azimuthal and multiplicity distributions for net-proton (the difference in number of protons and anti-protons) in Au+Au collisions at midrapidity for = 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV. We find the results to be very intriguing, which must to be further quantified by having a high event statistics second phase of BES program in near future at RHIC.

In the next section we discuss the two observables related to the search for the first order phase transition and the critical point. In section 3 we present the experimental results on the directed flow measurements of net-protons, an observable for first order phase transition. In section 4 we discuss the measurements related to the product of higher moments of net-proton multiplicity distribution, observable for critical point search. Finally in section 5 we summarize the findings.

2 Observables

The experimental results presented here are from the data recorded in the STAR detector at RHIC in the year 2010 and 2011.

2.1 Directed flow

The patterns of azimuthal anisotropy in particle production, often termed as flow, in heavy-ion collisions can be obtained by studying the Fourier expansion of the azimuthal angle () distribution of produced particles with respect to the reaction plane angle ([8]. Where the reaction plane is defined as the plane subtended by the impact parameter direction and the beam direction. The various (order ) coefficients in this expansion are defined as:


The angular brackets in the definition denote an average over many particles and events. Directed flow is quantified by the first coefficient (). On the other hand the elliptic flow is given by the second coefficient ().

The , which is sensitive to early collision dynamics, is proposed as a signature of first order phase transition based on a hydrodynamic calculation [9, 10, 11]. These calculations whose equation of state incorporates a first-order phase transition from hadronic matter to QGP predict a non-monotonic variation of the slope of directed flow of baryons (and net-baryons) around midrapidity as a function of beam energy, also has a prominent minimum, and a double sign change in the slope, which is not seen in the same hydrodynamic model without a first-order phase transition.

The results discussed in this paper are for the most abundantly measured baryons, antiproton and proton, detected using their energy loss in STAR Time Projection Chamber and by the time-of-flight information from the Time Of Flight detector [12]. Protons and antiprotons have a transverse momentum between 0.4 and 2.0 GeV/ and pseudorapidity () between 1 unit. The first order event plane for 62.4 GeV is constructed using the information from the two Beam Beam Counters at 3.3 5.0. That for = 62.4 and 200 GeV uses the information from STAR ZDC-SMD detectors. All results discussed here are corrected for event plane resolution and proton track reconstruction efficiency.

2.2 Higher moments

Non-monotonic variations of observables related to the moments of the distributions of conserved quantities such as net-baryon, net-charge, and net-strangeness [13] number with are believed to be good signatures of a CP. The moments are related to the correlation length () of the system [14]. Finite size and time effects in heavy-ion collisions put constraints on the significance of the desired signals. A theoretical calculation suggests a non-equilibrium 2-3 fm for heavy-ion collisions [15]. Hence, it is proposed to study the higher moments (like skewness, = and kurtosis, = [] – 3 with = ) of distributions of conserved quantities due to a stronger dependence on  [14]. Further, products of the moments can be related to susceptibilities associated with the conserved numbers. The product of the net-baryon number distribution is related to the ratio of fourth order () to second order () baryon number susceptibilities [16, 17]. The ratio / is expected to deviate from unity near the CP. It has different values for the hadronic and partonic phases [17].

The higher moments of the net-proton multiplicity ( = ) distributions from Au+Au collisions discussed in this paper are for protons and anti-protons detected at midrapidity ( 0.5) in the range 0.4 0.8 GeV/. A good purity of the proton sample (better than 98%) for all beam energies is obtained in the momentum range studied. All results presented are corrected for proton reconstruction efficiency.

3 Search for the first order phase transition

Figure 2: (Color online) Directed flow slope () near mid-rapidity as a function of beam energy () for intermediate-centrality (10-40%) Au+Au collisions [12]. Panels (a), (b) and (c) shows the STAR experiment measurement for antiprotons, protons, and net-protons, respectively, along with corresponding calculations from the UrQMD model obtained with the same cuts and fit conditions. The systematic uncertainties on the measurements are shown as shaded bars. The dashed curves are a smooth fit to guide the eye.

Figure 2 shows the slope of directed flow versus rapidity () near midrapidity as a function of for antiprotons (panel (a)), protons (panel (b)) and net-proton (panel (c)) [12]. The results are for 10-40% Au+Au collision centrality. The antiproton has negative values and shows a monotonic increase with . The proton has positive values for = 7.7 GeV and negative values for rest of the energies studied. The proton dependence on is non-monotonic with a minimum around = 11.5 and 19.6 GeV. The energy dependence of proton involves an interplay between the directed flow of protons associated with baryon number transported from the initial beam rapidity to the vicinity of mid-rapidity, and the directed flow of protons from particle-antiparticle pairs produced near mid-rapidity. The importance of the pair production mechanism increases strongly with beam energy. It is important to distinguish between the two mechanisms before further conclusions can be drawn. The net-proton is expected to provide the contribution from protons associated with baryon number transport. Assuming antiproton directed flow as a proxy for the directed flow of pair produced protons, the proposed net-proton slope can be constructed from

where is the observed rapidity dependence of the ratio of antiprotons to protons at each beam energy [12]. The net-proton slope is shown as a function of in Fig. 2(c). The data shows a non-monotonic dependence on with a minimum around = 11.5 and 19.6 GeV. The values of slope changes sign twice, goes from positive at 7.7 GeV to negative at = 11.5 - 27 GeV, then again becomes positive for 39 GeV.

The corresponding UrQMD results [18, 19], which does not include any first order phase transition effects, in all the three cases of slope of antiproton, proton and net-proton shows a monotonic variation with . The slope values also does not agree with the measurements.

A possible interpretation of the changing sign of the slope is that it reflects a change in equation-of-state. At a given energy where the system undertakes a first order quark–hadron phase transition, one expect the formation of a mixed-phase for which the pressure gradient is small. The softest pressure will naturally produces the observed minimum in slope parameter  [9, 10, 11]. The alternate proposal is that at higher energies pair production is dominant at mid-rapidity and transported baryons have relatively small influence. As there is no preferred direction for pair produced hadrons, so the slope parameter becomes close to zero. While at lower beam energies the observed baryons are strongly influenced by the transported baryons so they are aligned with them hence the slope parameter is positive. The most intriguing is the dependence in the intermediate energies, a mean field model study shows that the energy dependent baryon potential plays an important role in this region [20].

In the search for the signature of a first-order phase transition in the high region of the QCD phase diagram, the findings from from the RHIC beam energy scan program are very compelling and strongly motivates further measurements. To better understand the possible role and relevance of stopping in interpretation of existing data on net-proton directed flow, new higher event statistics measurements of the centrality dependence of at = 7.7 to 19.6 GeV is needed.

4 Search for the critical point

Figure 3: (Color online) Collision energy and centrality dependence of the net-proton and from Au+Au and + collisions at RHIC [22]. Crosses, open squares and filled circles are for the efficiency corrected results of +, 70-80%, and 0-5% Au+Au collisions, respectively. Skellam distributions for corresponding collision centralities are shown in the top panel. Shaded hatched bands are the results from UrQMD. The hadron resonance gas model (HRG) values for and /Skellam are unity. The error bars are statistical and caps are systematic errors. For clarity, + and 70-80% Au+Au results are slightly displaced horizontally.

Figure 3 shows the energy dependence of and for for Au+Au collisions for two collision centralities (0-5% and 70-80%), corrected for () reconstruction efficiency [21, 22]. The Skellam expectations for are calculated using the data as . The values normalized to the corresponding Skellam expectations are shown in the bottom panel of Fig. 3. The Skellam expectations reflect a system of totally uncorrelated, statistically random particle production. The central collision data shows deviation from Skellam expectation with maximum deviation occurring for = 19.6 and 27 GeV. The corresponding results from + collisions at = 62.4 and 200 GeV are also shown and found to be similar to peripheral Au+Au collisions within the statistical errors. For = 19.6 and 27 GeV, differences are observed between the 0-5% central Au+Au collisions and the peripheral collisions. The results are closer to unity for = 7.7 GeV. Higher statistics data for 19.6 GeV will help in quantitatively understanding the suggestive non-monotonic energy dependence of and .

The data also show deviations from the hadron resonance gas model [23, 24] which predict and /Skellam to be unity. The effect of decay is less than 2% as per the HRG calculations in Ref. [24]. To understand the effects of baryon number conservation [25] and experimental acceptance, UrQMD model calculations (a transport model which does not include a CP) [18, 19] for 0-5% Au+Au collisions are shown in the middle and bottom panels of Fig. 3. The UrQMD model shows a monotonic decrease with decreasing beam energy [26]. The centrality dependence of the and from UrQMD [26] (not shown in the figures) closely follow the data at the lower beam energies of 7.7 and 11.5 GeV. Their values are in general larger compared to data for the higher beam energies. The observed dependence for central Au+Au collisions is not explained by UrQMD model.

The current data provide the most relevant measurements over the widest range in (20 to 450 MeV) to date for the CP search, and for comparison with the baryon number susceptibilities computed from QCD to understand the various features of the QCD phase structure [16, 17]. The deviations of and below Skellam expectation are qualitatively consistent with a QCD based model which includes a CP [27]. However, conclusions on the existence of CP can be made only after a high event statistics measurement for 27 GeV in the second phase of the beam energy scan program and comparison to QCD calculations with CP behavior.

5 Summary

We have discussed two striking observations from the RHIC beam energy scan program related to the first order quark–hadron phase transition and the critical point. The measurements uses the produced protons and antiprotons in Au+Au collisions at midrapidity for = 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV.

The slope of the directed flow of protons and net-protons in mid-central collisions (10-40% centrality) at midrapidity () shows a clear non-monotonic variation with respect to (). The minimum value of lies somewhere between () = 27 GeV(160 MeV) to 11.5 GeV(315 MeV). The net-proton changes sign twice in the beam energy range studied. This observable which is driven by the pressure gradients developed in the system is sensitive to first order phase transition effects. The energy dependence of the measured is consistent with a theoretical hydrodynamic model calculation with first order phase transition [9, 10, 11].

Deviations of and for net-proton distribution in 0-5% centrality is observed at = 19.6 and 27 GeV from: (a) 70-80% peripheral collisions, (b) Poisson and hadron resonance gas expectation value of close to unity and (c) transport model based UrQMD calculation within the experimental acceptance. The deviations of and below Poisson expectation are qualitatively consistent with a QCD based model which includes a CP [27]. Higher statistics data set at 20 GeV in the second phase of the beam energy scan program is needed to clarify whether the energy dependence of the observable will follow a non-monotonic variation with a minimum around = 27 GeV to 11.5 GeV, as observed for the net-proton , or a monotonic variation with .

Acknowledgement: This work is supported by the DST Swarnajayanti Fellowship of Government of India.


  1. I. Arsene et al., Nucl. Phys. A 757, 1 (2005); B. B. Back et al., Nucl. Phys. A 757, 28 (2005); J. Adams et al., Nucl. Phys. A 757, 102 (2005); K. Adcox et al., Nucl. Phys. A 757, 184 (2005); S. Chatrchyan et al., JHEP 05, 063 (2012); B. Abelev et al., Phys. Rev. Lett. 109, 072301 (2012).
  2. Y. Aoki et al., Nature 443, 675 (2006).
  3. M. Luzum and P. Romatschke, Phys. Rev. C 78 (2008) 034915 [Erratum-ibid. C 79 (2009) 039903] [arXiv:0804.4015 [nucl-th]].
  4. H. Song, S. A. Bass, U. Heinz, T. Hirano and C. Shen, Phys. Rev. Lett. 106 (2011) 192301 [Erratum-ibid. 109 (2012) 139904] [arXiv:1011.2783 [nucl-th]].
  5. I. Vitev and M. Gyulassy, Phys. Rev. Lett. 89 (2002) 252301 [hep-ph/0209161].
  6. S. A. Bass, C. Gale, A. Majumder, C. Nonaka, G. -Y. Qin, T. Renk and J. Ruppert, Phys. Rev. C 79 (2009) 024901 [arXiv:0808.0908 [nucl-th]].
  7. K. Fukushima and T. Hatsuda, Rept. Prog. Phys. 74 (2011) 014001 [arXiv:1005.4814 [hep-ph]].
  8. A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58 (1998) 1671 [nucl-ex/9805001].
  9. D. H. Rischke, Y. Pursun, J. A. Maruhn, H. Stoecker and W. Greiner, Heavy Ion Phys. 1 (1995) 309 [nucl-th/9505014].
  10. J. Brachmann, S. Soff, A. Dumitru, H. Stoecker, J. A. Maruhn, W. Greiner, L. V. Bravina and D. H. Rischke, Phys. Rev. C 61 (2000) 024909 [nucl-th/9908010].
  11. H. Stoecker, Nucl. Phys. A 750 (2005) 121 [nucl-th/0406018].
  12. L. Adamczyk et al. [STAR Collaboration], arXiv:1401.3043 [nucl-ex].
  13. V. Koch, A. Majumder and J. Randrup, Phys. Rev. Lett. 95, 182301 (2005); M. Asakawa, U. W. Heinz and B. Muller, Phys. Rev. Lett. 85, 2072 (2000).
  14. M. A. Stephanov, Phys. Rev. Lett. 102, 032301 (2009); C. Athanasiou et al., Phys. Rev. D 82, 074008 (2010).
  15. B. Berdnikov and K. Rajagopal, Phys. Rev. D 61, 105017 (2000).
  16. S. Ejiri, F. Karsch and K. Redlich, Phys. Lett. B 633, 275 (2006); M. Cheng et al., Phys. Rev. D 79, 074505 (2009); B. Stokic et al., Phys. Lett. B 673, 192 (2009); S. Gupta et al., Science 332, 1525 (2011).
  17. R. V. Gavai and S. Gupta, Phys. Lett. B 696, 459 (2011).
  18. S. A. Bass, M. Belkacem, M. Bleicher, M. Brandstetter, L. Bravina, C. Ernst, L. Gerland and M. Hofmann et al., Prog. Part. Nucl. Phys. 41 (1998) 255 [Prog. Part. Nucl. Phys. 41 (1998) 225] [nucl-th/9803035].
  19. M. Bleicher, E. Zabrodin, C. Spieles, S. A. Bass, C. Ernst, S. Soff, L. Bravina and M. Belkacem et al., J. Phys. G 25 (1999) 1859 [hep-ph/9909407].
  20. M. Isse, A. Ohnishi, N. Otuka, P. K. Sahu and Y. Nara, Phys. Rev. C 72 (2005) 064908 [nucl-th/0502058].
  21. M. M. Aggarwal et al. [STAR Collaboration], Phys. Rev. Lett. 105 (2010) 022302 [arXiv:1004.4959 [nucl-ex]].
  22. L. Adamczyk et al. [STAR Collaboration], Phys. Rev. Lett. 112 (2014) 032302 arXiv:1309.5681 [nucl-ex].
  23. F. Karsch and K. Redlich, Phys. Lett. B 695, 136 (2011).
  24. P. Garg, D. K. Mishra, P. K. Netrakanti, B. Mohanty, A. K. Mohanty, B. K. Singh and N. Xu, Phys. Lett. B 726, 691 (2013).
  25. A. Bzdak, V. Koch and V. Skokov, Phys. Rev. C 87, 014901 (2013).
  26. X. Luo, J. Xu, B. Mohanty and N. Xu, J. Phys. G 40, 105104 (2013).
  27. M. A. Stephanov, J. Phys. G 38, 124147 (2011).
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Comments 0
Request comment
The feedback must be of minumum 40 characters
Add comment
Loading ...

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