![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0495v2-figure-252d9781-451e-4625-84e6-f773946ab68a.jpg)

Nuclear Physics A

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0495v2-figure-0762b46b-996b-401b-88c6-94ba0dd641eb.jpg)

Nuclear Physics A 00 (2014) 1-5

## Baseline for the cumulants of net-proton distributions at STAR

## Xiaofeng Luo a,1 , Bedangadas Mohanty b , Nu Xu a,c

a Institute of Particle Physics and Key Laboratory of Quark &amp; Lepton Physics (MOE), Central China Normal University, Wuhan, 430079, China. b School of Physical Sciences, National Institute of Science Education and Research, Bhubaneswar 751005, India c Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA.

## Abstract

We present a systematic comparison between the recently measured cumulants of the net-proton distributions by STAR for 0-5% central Au + Au collisions at √ s NN = 7.7-200 GeV and two kinds of possible baseline measure, the Poisson and Binomial baselines. These baseline measures are assuming that the proton and anti-proton distributions independently follow Poisson statistics or Binomial statistics. The higher order cumulant net-proton data are observed to deviate from all the baseline measures studied at 19.6 and 27 GeV. We also compare the net-proton with net-baryon fluctuations in UrQMD and AMPT model, and convert the net-proton fluctuations to net-baryon fluctuations in AMPT model by using a set of formula.

Keywords: QCD Critical Point, Higher Moments, Net-proton, Heavy-ion Collision, Quantum Chromodynamics

## 1. Introduction

In the phase diagram of Quantum ChromoDynamics (QCD-theory of strong interactions), it is conjectured on the basis of theoretical calculations that there will be a QCD Critical Point (CP) at high temperature ( T ) and non-zero baryonic chemical potential ( µ B ) [1]. Since the ab initio Lattice QCD calculation meet the notorious sign problem, there are still large uncertainties in theoretically determining the location of the CP in the QCD phase diagram [2, 3]. Di ff erent QCD based models also give very di ff erent results [4]. Finding the existence of a CP experimentally will be an excellent test of QCD theory in the non-perturbative region and a milestone of exploring the QCD phase diagram [5]. This is one of the main goals of the Beam Energy Scan (BES) Program at the Relativistic Heavy Ion Collider (RHIC). By tuning the colliding energy of gold nucleus from √ sNN = 200 GeV down to 7.7 GeV, one can access a broad region of the QCD Phase Diagram (20 &lt; µ B &lt; 420 MeV) [6]. Due to the high sensitivity to the correlation length ( ξ ) of the dynamical system [7, 8, 9] and direct connection to the susceptibilities in theoretical calculations, for example, the Lattice QCD calculations [5, 10], higher moments of multiplicity distributions of conserved quantities, such as net-baryon, net-charge and net-strangeness, have been applied to search for the QCD critical point in heavyion collision experiments [11]. As the correlation length will diverge near the CP, the non-monotonic variation of the moments of multiplicity distribution with respected to the colliding energy is the golden signature of the CP. To extract the CP induced experimental signal, it is crucial to understand the non-CP physics e ff ects in heavy-ion collisions on the experimental observable, such as the e ff ects of conservations for charges (electric, baryon number and strangeness number), finite size, resonance decay and hadronic scattering. Moreover, proper baseline needs to be constructed for experimental observables to search for the CP.

1 xfluo@mail.ccnu.edu.cn

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0495v2-figure-ec8c5a44-ceaf-4510-a897-d2ecef4c7e17.jpg)

part

Figure 1. (Color Online) E ffi ciency corrected mean values of the net-proton, proton and anti-proton distributions as a function of the average number of participant nucleons ( &lt; Npart &gt; ) in Au + Au collisions at √ sNN = 7.7-200 GeV [18]. The dashed lines are used to guide eyes.

In this paper, we will make the comparison between the baselines (Poisson, Binomial) and recently measured cumulants of net-proton distribution published by the STAR Collaboration and discuss the deviations of the data from the baselines. In addition, the results from the AMPT [12] and UrQMD [13] models will be discussed.

## 2. Results and Discussion

Recently, the STAR Collaboration has published the beam energy dependence of higher moments of the netproton (as proxy of net-baryon [14]) multiplicity distributions from RHIC BES Au + Au collision data [15]. The protons and anti-protons are identified with ionization energy loss in the Time Projection Chamber (TPC) of the STAR detector within transverse momentum range 0 4 . &lt; pT &lt; 0 8 GeV c and at mid-rapidity . / | y | &lt; 0 5. . Several analysis techniques [16, 17] have been applied, such as defining a new collision centrality with charged particle multiplicity within large pseudo-rapidity ( | η | &lt; 1) and excluding the particle of interest used in the analysis, centrality bin width correction (CBWC) and e ffi ciency correction. These address the e ff ects of auto-correlation, volume fluctuation and finite tracking e ffi ciency in the moments calculation. Figure 1 shows the mean of the net-proton, proton and antiproton distributions are almost linearly increasing with the average number of participant nucleons for each energy. The di ff erent mean values of the proton and anti-proton distributions at each energy are determined by the interplay between baryon number pair production and baryon stopping e ff ects. At high energies, the pair production process dominates the production of protons and anti-protons at mid-rapidity, while at low energies the e ff ect of baryon stopping is more important than at high energies. In the following, we discuss some expectations for cumulants of net-proton multiplicity distributions from some basic distributions.

- 1.
- Poisson : If the protons and anti-protons are independently distributed as Poissonian distributions. Then the net-proton multiplicity will follow the Skellam distribution, which is expressed as: P N ( ) = ( Mp Mp ) N / 2 I N (2 √ MpMp ) exp[ -( Mp + Mp )] , where the N is the net-proton number, I N ( x ) is a modified Bessel function, Mp and Mp are the mean number of protons and anti-protons, as shown in Fig. 1. The various order cumulants ( Cn ) are closely connected with the moments, e.g., C 1 = &lt; N &gt; = M C , 2 = &lt; ( ∆ N ) 2 &gt; = σ , 2 C 3 = &lt; ( ∆ N ) 3 &gt; = S σ , 3 C 4 = &lt; ( ∆ N ) 4 &gt; -3 &lt; ( ∆ N ) 2 &gt; 2 = κσ 4 , where the ∆ N = N -&lt; N &gt; , the σ 2 , S and κ are variance, skewness and kurtosis, respectively. Then we construct, S σ = C 3 / C 2 = ( Mp -Mp ) / ( Mp + Mp ) and κσ 2 = C 4 / C 2 = 1, which provides the Poisson expectations for the various order cumulants moments of net-/ proton distributions. The only input parameters of the Poisson baseline for cumulants of net-proton distributions are the mean values of the proton and anti-proton distributions.
- 2. Binomial Negative Binomial: / If the protons and anti-protons are independently distributed as Binomial or Negative Binomial distributions (BD NBD). Then various order cumulants of the net-proton distributions can be /

expressed in term of cumulants of the proton and anti-proton distributions: C net -p n = C p n + ( -1) n C p ¯ n . The first four order cumulants can be written as: C x 2 = σ 2 x = ε µ , x x C x 3 = S x σ 3 x = ε µ x x (2 ε x -1) , C x 4 = κ σ x 4 x = ε µ x x (6 ε 2 x -6 ε x + 1) , where ε x = σ /µ 2 x x , µ x = Mx , Mx is the mean values of protons or anti-protons distributions, x = p or ¯. p ε x &lt; 1 means the underlying distributions of protons or anti-protons are Binomial distributions, while ε x &gt; 1 gives Negative Binomial distributions. The input parameters for BD NBD expectations are the measured mean and / variance of the proton and anti-proton distributions.

Figure 2. Energy dependence of cumulants ( C 1 -C 4) of net-proton distributions for 0 -5% Au + Au collisions [18]. The error bars are statistical only. The dashed lines in the left panels are the Poisson expectations, while those are the Binomial expectations in the right panels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0495v2-figure-76a003bb-63d9-40d1-a775-cb19de506f6c.jpg)

Figure 2 shows the comparison between the cumulants of net-proton distributions and the Poisson Binomial base-/ lines. For the Poisson case, when the order of cumulant increases, the deviations of the data from the Poisson expectations for net-proton and proton increase. Largest deviations are found for C 4 at 19.6 and 27 GeV. For the Binomial case, the agreements persist up to the third order, while it fails to describe the net-proton and proton fourth order cumulants at 19.6 and 27 GeV. It is understandable that the Binomial baseline describes the data better, as it uses the variance as an input parameter, in addition to the mean values, and as the variance is correlated with higher order cumulants. The cumulants of anti-proton distributions can be described by the Poisson and Binomial baselines very well. More baseline discussions from Hadronic Resonance Gas model and transport model UrQMD can be found in [19]. The findings are consistent with ours. As we discussed that the binomial distributions can be used as baseline comparison for net-proton multiplicity distributions, an alternative interpretation, based on the negative binomial distribution, was given in [20]. The negative binomial distributions used there is not appropriate for the net-proton multiplicity distributions.

The STAR experiment measures net-proton fluctuations instead of net-baryon fluctuations and one may want to know to what extend they can reflect the net-baryon fluctuations in heavy-ion collisions. Therefore, fig. 3 demonstrates the comparison between moments of net-proton and net-baryon distributions from AMPT and UrQMD model calculations. We can find that the κσ 2 of net-baryon distributions are systematically lower than the net-proton results. The di ff erences are even bigger for low energies than high energies. There are two possible e ff ects for the di ff erence between net-proton and net-baryon fluctuations, one is the non inclusion of neutrons in the net-proton fluctuations, and the other one is the nucleon isospin exchanging process due to ∆ resonance formation via p π and n π interaction, the so called isospin randomization, which will modify the net-proton fluctuations after the chemical freezeout. Asakawa and Kitazawa have derived a set of formulas [21] to convert the measured net-proton cumulants to the net-baryon cumulants by taking into account the above two e ff ects. The converting formulas for various order net-baryon cumulants can be written as:

$$C _ { 1 } ^ { n e t - B } & = 2 C _ { 1 } ^ { n e t - p }, \quad C _ { 2 } ^ { n e t - B } = 4 C _ { 2 } ^ { n e t - p } - 2 C _ { 1 } ^ { l o t - p }, \\ C _ { 3 } ^ { n e t - B } & = 8 C _ { 3 } ^ { n e t - p } - 1 2 ( C _ { 2 } ^ { p } - C _ { 2 } ^ { \flat } ) + 6 C _ { 1 } ^ { n e t - p }, \\ C _ { 4 } ^ { n e t - B } & = 1 6 C _ { 4 } ^ { n e t - p } + 1 6 C _ { 3 } ^ { l o t - p } - 6 4 ( C _ { 3 } ^ { p } + C _ { 3 } ^ { \flat } ) + 4 8 C _ { 2 } ^ { n e t - p } + 1 2 C _ { 2 } ^ { l o t - p } - 2 6 C _ { 1 } ^ { l o t - p },$$

where tot -p means proton number plus anti-proton number. The right side of fig. 3 shows the net-baryon κσ 2 ( C 4 / C 2 ) results, converted from the net-protons fluctuations. Within large uncertainties, they are consistent with the net-baryon results directly calculated with the AMPT model.

Figure 3. Energy dependence of κσ 2 of net-proton and net-baryon distributions for 0-5% Au + Au collisions from the UrQMD (left) and the AMPT string melting model (right). The results marked as solid black stars are based on theoretical calculations using Asakawa and Kitazawa's formula. The error calculation is based on the Bootstrap method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0495v2-figure-427b500a-8472-472f-8adb-76067f17b97e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0495v2-figure-09c1c403-d0be-42c4-ace8-25d7050deb61.jpg)

## 3. Summary

We have compared the energy dependence of cumulants of net-proton distributions in 0-5% central Au + Au collisions with Poisson and Binomial baselines. The Binomial baseline describes the data better than the Poisson baseline. The largest deviation is observed for the fourth order net-proton cumulant ( C 4 ) from Poisson and Binomial baselines in the most central collisions at 19.6 and 27 GeV . A proper comparison of the experimental measurements to QCD calculations is needed to extract the exact physics process that leads to deviation of the data from the baselines presented. Di ff erences between net-proton and net-baryon fluctuations are observed in the AMPT and UrQMD models. Finally, we use the A&amp;K's formulas to convert the net-proton fluctuations in the AMPT model to net-baryon fluctuations, which is consistent with the net-baryon results directly calculated from the AMPT model.

## Acknowledgement

The work was supported in part by the MoST of China 973-Project No.2015CB856901, NSFC under grant No. 11205067, 11221504 and 11228513. China Postdoctoral Science Foundation (2012M511237, 2013T60732).

## References

- [1] K. Rajagopal and F. Wilczek. 'At the Frontier of Particle Physics / Handbook of QCD' , volume 3. World Scientific, 2001.
- [2] Z. Fodor, S. Katz, JHEP 0203, 014 (2002); 0404, 050 (2004); R. V. Gavai, S. Gupta, Phys. Rev. D 71 , 114014 (2005); 78 , 114503 (2008).
- [3] R. V. Gavai, arXiv:1404.6615.
- [4] M. A. Stephanov, Int. J. Mod. Phys. A 20 , 4387 (2005).
- [5] S. Gupta, X. Luo, B. Mohanty, H. G. Ritter, N. Xu, Science 332 , 1525 (2011).
- [6] M. M. Aggarwal et al. (STAR Collaboration), arXiv: 1007.2613.
- [7] M. A. Stephanov, Phys. Rev. Lett. 102 , 032301 (2009).
- [8] C. Athanasiou et al. , Phys. Rev. D 82 , 074008 (2010).
- [9] M. A. Stephanov, Phys. Rev. Lett. 107 , 052301 (2011).
- [10] R. V. Gavai and S. Gupta, Phys. Lett. B 696 , 459 (2011).
- [11] L. Adamczyk et al. , (STAR Collaboration), Phys. Rev. Lett. 105 , 022302 (2010); Phys. Rev. Lett. 112 , 032302 (2014); arXiv:1402.1558.
- [12] Z.W. Lin et al. , Phys. Rev. C 72 , 064901 (2005).
- [13] M. Bleicher et al. , J. Phys. G: Nucl. Part. Phys. 25 , 1859 (1999).
- [14] Y. Hatta and M. A. Stephanov, Phys. Rev. Lett. 91 , 102003 (2003).

- [15] L. Adamczyk, et al. , (STAR Collaboration), Phys. Rev. Lett. 112 , 032302 (2014).
- [16] X. Luo, J. Xu, B. Mohanty, N. Xu, J. Phys. G 40 , 105104 (2013); X. Luo, J. Phys.: Conf. Ser. 316 , 012003 (2011).
- [17] X. Luo, J. Phys. G: Nucl. Part. Phys. 39 , 025008 (2012).
- [18] STAR Public Webpage for Data Table: https: // drupal.star.bnl.gov STAR files starpublications 205 data.html. / / / / /
- [19] P. K. Netrakanti, X. Luo, D. K. Mishra, B. Mohanty, A. Mohanty, N. Xu , arXiv: 1405.4617.
- [20] T. J. Tarnowsky and G. D. Westfall, Phys. Lett. B 724 , 51 (2013).
- [21] M.Kitazawa and M. Asakawa, Phys. Rev. C 85 , 021901 (2012); Phys. Rev. C 86 ,024904 (2012) [Erratum-ibid. C 86 (2012) 069902].
Xiaofeng Luo Bedangadas Mohanty Nu Xu
Baseline for the cumulants of net-proton distributions at STAR
We present a systematic comparison between the recently measured cumulants of
the net-proton distributions by STAR for 0-5% central Au+Au collisions at
$\sqrt{s_{NN}}$=7.7-200 GeV and two kinds of possible baseline measures, the
Poisson and Binomial baselines. These baseline measures are assuming that the
proton and anti-proton distributions independently follow Poisson statistics or
Binomial statistics. The higher order cumulant net-proton data are observed to
deviate from all the baseline measures studied at 19.6 and 27 GeV. We also
compare the net-proton with net-baryon fluctuations in UrQMD and AMPT model,
and convert the net-proton fluctuations to net-baryon fluctuations in AMPT
model by using a set of formula.