Advances in High Energy Physics

Volume 2013 (2013), Article ID 574871, 22 pages

http://dx.doi.org/10.1155/2013/574871

## On the Higher Moments of Particle Multiplicity, Chemical Freeze-Out, and QCD Critical Endpoint

^{1}Egyptian Center for Theoretical Physics (ECTP), MTI University, Cairo 11571, Egypt^{2}Research Center for Einstein Physics, Freie-University Berlin, Berlin, Germany

Received 6 January 2013; Revised 11 May 2013; Accepted 2 June 2013

Academic Editor: Edward Sarkisyan-Grinbaum

Copyright © 2013 A. Tawfik. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We calculate the first six nonnormalized moments of particle multiplicity within the framework of the hadron resonance gas model. In terms of the lower order moments and corresponding correlation functions, general expressions of higher order moments are derived. Thermal evolution of the first four normalized moments and their products (ratios) are studied at different chemical potentials, so that it is possible to evaluate them at chemical freeze-out curve. It is found that a nonmonotonic behaviour reflecting the dynamical fluctuation and strong correlation of particles starts to appear from the normalized third order moment. We introduce novel conditions for describing the chemical freeze-out curve. Although the hadron resonance gas model does not contain any information on the criticality related to the chiral dynamics and singularity in the physical observables, we are able to find out the location of the QCD critical endpoint at MeV and temperature MeV.

#### 1. Introduction

Recently, the higher order multiplicity moments have gained prominence in high energy physics with a huge hope in pinpointing the QCD critical endpoint (CEP) [1] connecting the first order boundary separating the hadronic from the partonic matter at high density with the cross-over boundary at low density [2, 3]. It is clear that the hadron resonance gas partition function is an approximation to a nonsingular part of the free energy of QCD in the hadronic phase. There are large theoretical uncertainties of its location, and even its existence is not fully confirmed yet [4–6]. The characteristic feature of CEP is critical dynamical fluctuations [7–14]. The higher order moments are conjectured to reflect the large fluctuations associating the hadron-quark phase transition. This was the motivation of a remarkable number of experimental and theoretical studies [1, 15–20]. Recently, various calculations have shown that the higher order moments of the multiplicity distributions of some conserved quantities, such as net-baryon, net-charge, and net-strangeness, are sensitive to the correlation length [21–23], which in turn is related to the higher order moments, themselves. In realistic heavy-ion collisions, the correlation length is found to remain finite.

Apparently, it should not be an exception that the strongly interacting QCD matter undergoes phase transition(s) at extreme conditions, as almost all matter types suffer from such a critical change as the extreme conditions change. The lattice QCD calculations predict that a cross-over takes place between the hadronic phase and the Quark Gluon Plasma (QGP) phase, when the temperature exceeds critical value of MeV. As given in [3, 24, 25], depending on the different parameters (for instance, the quark favours and their masses) and on the lattice configurations, lattice QCD assigns different values to . Furthermore, the value of depends on the baryon chemical potential . With vanishing , the lattice QCD calculations [26] show that the higher order susceptibilities can be related to the higher order moments of the corresponding multiplicity distributions. They also show nonmonotonic behavior near . Apparently, the QCD phase diagram including CEP can be mapped out using this characteristic behaviour. In the literature, one finds that the QCD phase diagram has been determined using various criteria, such as critical energy density [27–29] or line of constant physics [3]. Also the chemical freeze-out curve can be characterized assuming constant physics. For a recent and extensive review, we refer to [30] and the references therein.

At large chemical potential , various calculations based on QCD-based models indicate that the transition from the hadronic phase to the QGP phase is of first order. The endpoint connecting this line with the one of the cross-over (at small ) is likely of second order [4–6, 31, 32]. Therefore, the phase transition should be accompanied with large fluctuations in different physical quantities. Experimentally, we study the QCD phase diagram by varying the colliding energy in heavy-ion collisions and determining the critical temperature [2]. The possibility of finding QCD CEP has motivated the RHIC beam energy scan program [33]. By tuning the collision energy from a center-of-mass energy of GeV down to GeV, we are able to vary the baryon chemical potential from ~2 to ~500 MeV. Nevertheless, there are large theoretical uncertainties of its location and even its existence is not fully confirmed [4–6, 34]. The higher order moments of various physical quantities are conjectured to reflect the nonmonotonic behaviour near the QCD CEP. In the present work, we give a systematic study for the higher order moments of the multiplicity distribution and show their ability to represent the fluctuations along the QCD phase diagram.

Many reasons speak for utilizing the physical resonance gas model (HRG) in predicting the hadron abundances and their thermodynamics. The HRG model seems to provide a good description for the thermal evolution of the thermodynamic quantities in the hadronic matter [3, 7–9, 27–29, 35, 36] and has been successfully utilized to characterize the conditions deriving from the chemical freeze-out at finite densities [37, 38]. In light of this, HRG can be used in calculating the higher order moments of particle multiplicity using a grand canonical partition function of an ideal gas with all experimentally observed states up to a certain large mass as constituents. The HRG grand canonical ensemble includes two important features [3]; the kinetic energies and the summation over all degrees of freedom and energies of resonances. On other hand, it is known that the formation of resonances can only be achieved through strong interactions [39]; *resonances (fireballs) are composed of further resonances (fireballs), which in turn consist of resonances (fireballs) and so on*. In other words, the contributions of the hadron resonances to the partition function are the same as that of free particles with some effective mass. At temperatures comparable to the resonance half-width, the effective mass approaches the physical one [3]. Thus, at high temperatures, the strong interactions are conjectured to be taken into consideration through including heavy resonances. It is found that hadron resonances with masses up to GeV are representing suitable constituents for the partition function [3, 7–9, 27–29, 35, 36]. Such a way, the singularity expected at the Hagedorn temperature [27, 28] can be avoided and the strong interactions are assumed to be considered. Nevertheless, the validity of HRG is limited to temperatures below the critical one, .

In this paper, we study the first six nonnormalized higher order moments of the particle multiplicity in the hadron resonance gas (HRG) model, which is introduced in Section 2. In Section 2.1, we suggest an answer to the question: what is needed when going from “trivial” lower order to “sophisticated” higher order moments? Based on results, general expressions for arbitrary higher moments are deduced, so that it is possible to conclude that going from lower to higher order moments is achievable through a series of lower order moments and correlation functions. The thermal evolution of the first four normalized moments and their products (ratios) are studied at different chemical potentials in Section 2.2. Therefore, it is possible to evaluate them at the chemical freeze-out curve, which is characterized by constant at all values of , where and are entropy density and temperature, respectively. Sections 3 and 4 are devoted to introduce novel conditions describing the chemical freeze-out curve and define the location of the QCD critical endpoint.

#### 2. Higher Order Moments of the Particle Multiplicity

In a grand canonical ensemble, it is straightforward to derive an expression for the pressure. In Section 2.1, we give a list of different moments of particle multiplicity. As we move from lower order to higher order multiplicities, certain distribution functions are being added/subtracted. These distribution functions represent higher order correlations. In Section 2.2, we list out various normalized higher order moments, the normalization being based on the variance . The products (ratios) of normalized moments are elaborated in Section 2.3.

The hadron resonances treated as a free gas [3, 27–29, 35] are conjectured to add to the thermodynamic pressure in the hadronic phase (below ). This statement is valid for free as well as for strongly interacting resonances. It has been shown that the thermodynamics of strongly interacting system can also be approximated to an ideal gas composed of hadron resonances with masses ≤2 GeV [3, 40]. Therefore, the confined phase of QCD, the hadronic phase, is modelled as a noninteracting gas of resonances. The grand canonical partition function reads
where is the Hamiltonian of the system and is the temperature. The Hamiltonian is given by the sum of the kinetic energies of relativistic Fermi and Bose particles. The main motivation of using this Hamiltonian is that it contains all relevant degrees of freedom of confined and strongly interacting matter. It includes implicitly the interactions that result in resonance formation. In addition, it has been shown that this model can submit a quite satisfactory description of particle production in heavy-ion collisions. With the above assumptions the dynamics of the partition function can be calculated exactly and be expressed as a sum over *single-particle partition* functions of all hadrons and their resonances. Consider the following:
where is the th particle dispersion relation, is spin-isospin degeneracy factor, and stands for bosons and fermions, respectively.

Before the discovery of QCD, a probable phase transition of a massless pion gas to a new phase of matter was speculated [41]. Based on statistical models like Hagedorn [42, 43] and Bootstrap [44, 45], the thermodynamics of such an ideal pion gas is studied, extensively. After the QCD, the new phase of matter is known as Quark Gluon Plasma (QGP). The physical picture was that at the additional degrees of freedom carried by QGP are to be released resulting in an increase in the thermodynamic quantities like energy and pressure densities. The success of hadron resonance gas model in reproducing lattice QCD results at various quark flavours and masses (below ) changed this physical picture drastically. Instead of releasing additional degrees of freedom at , the interacting system reduces its effective degrees of freedom at . In other words, the hadron gas has much more degrees of freedom than QGP.

At finite temperature and baryon chemical potential , the pressure of the th hadron or resonance species reads As no phase transition is conjectured in HRG, summing over all hadron resonances results in the final thermodynamic pressure in the hadronic phase.

The switching between hadron and quark chemistry is given by the relations between the *hadronic* chemical potentials and the quark constituents; , where being baryon (strange) quantum number. The chemical potential assigned to the light quarks is , and the one assigned to strange quark reads . The strangeness chemical potential is calculated as a function of and under the assumption that the overall strange quantum number has to remain conserved in heavy-ion collisions [3].

##### 2.1. Nonnormalized Higher Order Moments

As given in the introduction, the higher order moments can be studied in different physical quantities. For example, the higher order moments of charged-particle multiplicity distribution have been predicted four decades ago [46]. Recently, the higher order moments of various multiplicity distributions have been reported by STAR measurements [15–19] and lattice QCD calculations [47–49]. The empirical relevance of the higher order moments to the experimental measurements has been suggested in [1]. Accordingly, the measurement of the correlation length seems to be very much crucial. In a future work, we shall discuss the experimental signatures of the higher order moments. Another ingredient would be the experimental sensitivity for the suggested signatures, which are based on singular behavior of thermodynamic functions. For the time being, we just want to mention that the experimental measurements apparently take place at the final state of the collision, which would mean that the signals have to survive the extreme conditions in such high energy collisions. It has been pointed out that the contribution of the critical fluctuations to the higher order moments is proportional to a positive power of . The latter is conjectured to diverge at QCD CEP, and such an assumption is valid in the thermodynamic limit.

In the present work, we first study the contribution of the hadron resonance mass spectrum in clarifying the interrelations between the various nonnormalized and normalized moments. The motivation of using resonance gas is clear [27–29]. As introduced about five decades ago, the resonance gas is only valid below [3, 27–29, 35, 42, 43]. Therefore, the present work is valid in the hadronic phase, only.

For the th particle, the “first” order moment is given by the derivative of with respect to the dimensionless quantity . When taking into account the antiparticles, we add a negative sign to the chemical potential. The first derivative describes the multiplicity distribution or an expectation operator, which is utilized to estimate the number or multiplicity density The “second” order moment is known as the variance. It gives the susceptibility of the measurements. Consider the following: The “third” order moment measures the lopsidedness of the distribution. As given in Section 2.2, the normalization of third order moment is known as skewness or the asymmetry in the distribution. Skewness tells us about the direction of variation of the data set. Consider the following: In general, the normalization of th order moment is obtained by dividing it by , where is the standard deviation. The normalization is assumed to remove the brightness dependence.

The “fourth” order moment compares the tallness and skinniness or shortness and squatness, that is, shape, of a certain measurement to its normal distribution. It is defined as the uncertainty is an uncertainty or *“the location- and scale-free movement of probability mass from the shoulders of a distribution into its center and tails and to recognize that it can be formalized in many ways”* [50]. Consider the following:
The “fifth” order moment measures the asymmetry sensitivity of the “fourth” order moment. Consider the following:
The “sixth” order moment is generally associated with compound options. Consider the following:

Thus, from (4)–(9), a general expression for the th order moment can be deduced where the coefficients read where and vanishes for all . It is obvious that the coefficients of a certain moment are to be determined from a long chain of all previous ones. Such a conclusion dates back to about four decades [46], where it has been shown that the coefficients are related to high order correlation functions. Should this assumption is proven to be valid, then expression (11) gets a novel interpretation. It seems to sum up the correlation functions up to the th order. According to [46] and when neglecting the two-particle correlations , then the higher order moments read where , and is the th particle. To the second order moment , we have to add the effects of the two particle correlation function, . The third order moment gets approximately three times this amount. The three and four particle correlations, (14) and (15) appear first in third and fourth moments (13) and (14), respectively.

When ignoring these high order correlations, we are left with the widely used two particle correlations. In the present work, we restrict the discussion to this type, only. It is of great interest, as it is accessible experimentally and achievable numerically. The two particle correlations are suggested as a probe for the bulk QCD medium, energy loss, medium response, jet properties, and intensity interferometry [51–57]. In addition to this list of literature, the comprehensive review [58, 59] can be recommended.

Taking into account the particle multiplicities, then expressions (12) and (13) can be rewritten as follows: As mentioned in Section 2.1, the experimental signals for the higher order moments shall be discussed in a future work.

As given in (10), the dependence of first six nonnormalized moments on lower ones can explicitly be deduced as follows: Naively spoken, we conclude that raising lower order moments to higher ones is achievable through a series of all lower order moments and an additional term reflecting the correlations, themselves [46]. From (17), the additional terms are proportional to . Apparently, they are essential in order to judge whether going from lower to higher order moments would make the signatures of dynamical fluctuations clearer than when excluding them.

The first and second terms can be generalized. Then, a general expression would read where is the coefficient of the th moment. The last term in (18) can be elaborated when higher order moments are calculated. The latter are essential in order to find out a clear pattern.

Figure 1 presents the results of the higher order moments as calculated in the HRG model. As introduced previously, raising the orders of multiplicity moments results in new coefficients and new integrals. The earlier are partly characterized in (18). The dimensionless nonnormalized moments are given depending on the temperature in the left panel of Figure 1. The new coefficients that resulted in when going from lower to higher order moments are depicted in the right panel. The integrals can be compared with the correlation functions [46]. We find that these integrals are related to , where is the order of the moment. It is obvious that successive moments have a difference of about one order of magnitude. Therefore, higher order can be disregarded with reference to their vanishing contributions.

In order to draw a picture about the contributions of the integrals appearing in (17), the phase space integral can be replaced by a series representation. For the th particle species, the pressure can be expressed as where is modified Bessel function of the second kind and the factor represents fermions and bosons, respectively. In deriving previous expression, we expand the functions appearing in (4), for instance. The resulting expression is applicable as long as is valid. Then, the derivatives with respect to leads to which is valid for all order moments of order .

For completeness, we mention here that taking into consideration one particle and its antiparticle is not rare in the literature. In this case, the expression (19) can be rewritten as and accordingly, Comparing (20) with (19), we conclude that calculating higher order moments is associated with a gradual change in . Expression (10) reflects almost the same behaviour, where higher order moments are associated with change in the coefficients.

##### 2.2. Normalized Higher Order Moments

The normalization of higher order moments which can be deduced through derivatives with respect to the chemical potential of given charges, apparently gives additional insights about the properties of higher order moments. From statistical point of view, the normalization is done with respect to the standard deviation , which is be related to . Therefore, it provides a tool to relate moments with various orders to the experimental measurement. The susceptibility of the distribution gives a measure for . Should we are interested on the multiplicity, then the susceptibility is simply given by the derivative of first order moments with respect to . It has been shown that the susceptibility is related to [1]. The results of in hadronic, bosonic, and fermionic resonances are calculated at different and given in Figure 2. As discussed earlier, the strange quantum number has to remain conserved in high energy collisions. As per the standard model, this is one of the global symmetries in strong interactions. The procedure of keeping strange degrees of freedom conserved in HRG is introduced in [3]. This is the origin of the dependence, especially in the bosonic gas. Although baryon chemical potential vanishes per definition, the chemical potential associated with strange quark remains finite. The dependence of on is depicted in Figure 16. Another feature in these calculations is the assumption that the freeze-out boundary is determined by constant , where is the entropy density [37]. At chemical freeze-out boundary, the dependence of on is given in Figure 5.

The “third” order moment normalized to is known as “skewness.” For standard Gaussian distribution, the skewness is obviously vanishing. Therefore, the skewness is an ideal quantity probing the non-Gaussian fluctuation feature as expected near . The QCD CEP is conjectured to be sensitive to skewness. Experimentally, it has been shown that the skewness is related to [21]. The skewness for bosonic and fermionic resonance gas, respectively, reads At different chemical potentials, the skewness is calculated depending on and given in Figure 3. It is worthwhile to notice that the values of in the fermionic resonance gas are much larger than the values in the bosonic resonance gas, especially at small chemical potentials. But at large chemical potential, this observation is exactly the opposite. The value of , at which the large (or small) skewness is flipped, seems to reflect the nature of phase transition, where the fermionic and/or bosonic degrees of freedom turn out to be undistinguishable. Such an observation shall be utilized in positioning the QCD CEP, Section 3.

The normalization of fourth order moment is known as, heteroskedacity or kurtosis. It means varying volatility or more accurately, varying variance. Actually, the kurtosis is given by normalized fourth order moment minus . The subtraction of , which arises from the Gaussian distribution, is usually omitted [50, 60, 61]. Therefore, the kurtosis is an ideal quantity for probing the non-Gaussian fluctuation feature as expected near and CEP. A sign change of skewness or kurtosis is conjectured to indicate that the system crosses the phase boundary [26, 62, 63]. As HRG is valid below , the sign change is not accessible. It has been shown that kurtosis is related to [21]. Consider the following:

##### 2.3. Products of Higher Order Moments

There are several techniques to scale the correlation functions. The survey system’s optional statistics module represents the most common technique, that is, Pearson or product moment correlation. This module includes the so-called partial correlation which seems to be useful when the relationship between two variables is to be highlighted, while effect of one or two other variables can be removed. In the present work, we study the products of higher order moments of the distributions of conserved quantities. The justification of this step is that certain products can be directly connected to the corresponding susceptibilities in lattice QCD simulation and related to long range correlations [26, 64, 65]. Seeking for simplicity, we start with the Boltzmann approximation.

###### 2.3.1. Modified Bessel Function: Boltzmann Statistics

When relativistic momentum integrals are replaced by summation over modified Bessel functions, the products of higher order moments in Boltzmann approximation read The origin of the second term in all expressions containing is obvious. Expressions (25) justify the conclusions in [66] that are in Boltzmann approximation while These expressions are valid in the final state, which can be characterized by chemical and thermal freeze-out. In other words, these expressions are functions of the chemical potential. Using relativistic momentum integrals (see Section 2.3.2) shows that the products of moments result in a constant dependence on [66].

###### 2.3.2. Relativistic Momentum Integrals: Quantum Statistics

The fluctuations of conserved quantities are assumed to be sensitive to the structure of the hadronic system in its final state. As mentioned above, crossing the phase boundary or passing through critical endpoint is associated with large fluctuations. The most proposed fluctuation of observables is variations of second order moments of the distribution, such as particle ratio [11–14] and charged dynamical measurement [67]. Then, the fluctuations are approximately related to [23].

The ratio of standard deviation and the mean multiplicity for fermions and bosons read where stands for fermions and bosons, respectively. The results are given in Figure 4. We notice that the bosonic resonance gas results in smaller values than the fermionic one, especially in final state. Furthermore, we notice that decreases with increasing of the bosons at middle temperature. This is exactly the opposite in the fermionic resonance gas. This would give an explanation for the observation that the results in hadron resonances are not spread as in the other two sectors. Figure 5 shows the dependence of on at the chemical freeze-out boundary, which is characterized by . The ratio is equivalent to , (32). The dependence of bosonic and fermionic on is given in Figure 5, as well. We notice that is smaller than , especially at small . At large both quantities are almost equal.

The multiplication of skewness by the standard deviation is directly related to the thermodynamics of the number susceptibility of the lattice QCD. In HRG, the bosonic and fermionic products read The product is equivalent to (33). The results are given in Figure 6. It is obvious that for either bosons or fermions. Then, for hadrons, . Nevertheless, the fine structure seems to reveal interesting features.

The dependence of and on is illustrated in Figures 2 and 4, respectively. It is obvious that both quantities have a monotonic behavior. Their dependences on are given in Figure 5. Also, this type of dependences seems to be monotonic. As given in Figure 6, the product has a characteristic dependence on . Regardless of the tiny change, we notice that increasing increases the bosonic product but decreases the fermionic product. In both cases, it forms a folding fan. The hadronic product makes an amazing bundle in the middle (at a characteristic ). This will be discussed in details in Section 3. When is studied as a function of and given in Figure 15, we find that remains almost constant, while raises with increasing . We also notice that both curves cross at a certain point. This is not the case of the dependence of and , Figure 5. So far, we conclude that starting from the third normalized moment, a nonmonotonic behavior appears. Thus, the kurtosis and its products should be expected to highlight such a nonmonotonic behavior [26, 64, 68].

The multiplication of kurtosis by called [69] is apparently equivalent to the ratio of third order moment to second order moment (34). In lattice QCD and QCD-like models, is found to diverge near the critical endpoint [26, 64]. In HRG, the bosonic and fermionic products read When ignoring the constant term in (24), then the second terms in the previous expressions entirely disappear. The results are given in Figure 7. In the hadronic sector, the dependence of on the temperature at different values is given in Figure 7(a). We notice that increasing is accompanied with a drastic declination in . Also, we find that flips its sign at large . When comparing the thermal evolution of with the one of , Figure 6, we simply find that the latter is much drastically changing than the former. Also, when comparing their dependence on the chemical potentials at the freeze-out boundary, Figures 9 and 14, it is apparent that the dependence increases when the normalized fourth order moment is included. The product calculated at the freeze-out boundary leads to some interesting findings. First, almost vanishes or even flips its sign. Second, the and corresponding to vanishing are coincident with the phenomenologically measured freeze-out parameters. Third, the freeze-out boundaries of bosons and fermions are crossing at a point located very near to the one assumed by the lattice QCD calculations to be the QCD CEP, Section 4.

The skewness and kurtosis can be combined through the standard deviation as follows: Also when ignoring the constant term in (24), the second terms in previous expressions disappear. The results of are given in Figure 8. In the fermionic sector, decreases with increasing . Increasing makes the decrease much faster. Despite the drastic change in , this behavior can be compared with , Figure 6, qualitatively. In the bosonic sector, there is a very slow decrease with increasing . The same sector in , Figure 6, shows a very slow but an increasing dependence on . Therefore, in the hadronic sector, is overall diminishing with increasing . The declination becomes faster for larger . Figure 9 illustrates the dependence, that is, is estimated at the freeze-out curve. We notice that the bosons result in an almost unvarying with growing . Therefore, the overall decrease in the hadronic sector is originated in the fermionic degrees of freedom.

From (28)–(31), we summarize that Obviously, other products would complete missing ratios, for example, relating third and fourth order moments to the first one, respectively. As discussed above, the susceptibility is equivalent to . These products seem to be sensitive to the volume independent multiplicities.

The results of are given in Figure 10. At the freeze-out boundary, the dependence of seems to be much stronger than the dependence on , Figure 11.

The thermal evolution of is illustrated in Figure 12. Increasing is accompanied by an almost overall increase of values. At low , fluctuations appear. In bosonic sector, increasing leads to an obvious decrease in so that the curves are distributed over a wide range. This is not the case for the fermionic sector. Figure 12(a) shows the results in the hadronic sector. At the freeze-out boundary, the values of at various -values are given in Figure 13. At large , both bosonic and fermionic ratios vanish. With decreasing , there is a bounce. It reaches its maximum over ~0.25 GeV. Then both curves come closer to each other. They cross again over ~0.2 GeV. Continue decreasing results in another bounce off. At smaller , the fermionic sector switches again to positive values, while the bosonic sector results in large negative values of . For hadrons, remains almost negative. At very low , it turns out to small positive values.

#### 3. Chemical Freeze-Out

In the rest frame of produced particle, the hadronic matter can be determined by constant degrees of freedom, for instance, [70]. The chemical freeze-out is related to the particle creation. Therefore, the abundance of different particle species is controlled by chemical potential, which obviously depends on . With the beam energy, is increasing, while the baryon densities at midrapidity are decreasing. The estimation of the macroscopic parameters of the chemical freeze-out can be extracted from particle ratio. These parameters collected over the last two decades seem to follow regular patterns as the beam energy increases [30, 37]. The higher order moments have been suggested to control the chemical freeze-out, so that several conditions have been proposed [71–73].

As introduced in Section 2.3.2, the thermal evolution of declines. This result seems to update previous ones [48, 74], where is assumed to remain finite and positive with increasing . In the present work, we find that the sign of is flipped at high [71–73]. It has been found that the and parameters, at which the sign is flipped, amazingly coincide with the ones of the chemical freeze-out, Figure 14. Also, we find that the freeze-out boundaries of bosons and fermions are crossing at one point located at the hadronic curve. This point is close to the one that lattice QCD calculations suggest for QCD CEP, Section 4. The results are given in Figure 14. Vanishing leads to The rhs and lhs in both expressions can be rewritten as which is valid for bosons and fermions. Then, the chemical freeze-out is defined, if the condition is fulfilled. At chemical freeze-out curve, a naive estimation leads to fm. In doing this, it is assumed that the proportionality coefficients of and , are equal. In the heavy-ion collisions, has been measured [75]. Near a critical point, the experimental value ~2-3 fm (only factor 3 larger) agrees well with our estimation.

At the chemical freeze-out curve, the intensive parameters and which are related to the extensive properties entropy and particle number, respectively, have to be determined over a wide range of beam energies. Left panel of Figure 14 collects a large experimental data set. For a recent review, we refer to [30] and the references therein. The double-dotted curve represents a set of and , at which vanishes. It is obvious that this curve reproduces very well the experimental data. As given above, at this curve the normalized fourth order moment is equal to three times the squared susceptibility . This new condition seems to guarantee the condition introduced in [37]; , right panel of Figure 14, at least over the range GeV. When excluding all degrees of freedom but the fermionic ones, decreases with increasing (dashed curve). An opposite dependence is accompanied with the bosonic degrees of freedom (dotted curve).

#### 4. QCD Critical Endpoint

In the right panel of Figure 14, it is interesting to notice that both fermionic and bosonic curves intersect at the hadronic curve at one point. Also, in the freeze-out diagram, all curves crossing at one point, left panel of Figure 14. Furthermore, based on nonperturbative convergence radius, the critical endpoint as calculated in lattice QCD [5, 6] given by solid square in Figure 14 is located very near to the crossing point. It is clear that the HRG model does not contain any information on criticality related with the chiral dynamics and singularity in physical observables required to locate the CP. The HRG partition function is an approximation to a nonsingular part of the free energy of QCD in the hadronic phase. It can be used only as the reference for LGT calculations or HIC to verify the critical behavior, but it can not be used as an origin to search for the chiral critical structure in QCD medium. This is the motivation of Figure 15.

In Figure 15, the skewness and kurtosis of bosons and fermions calculated at the freeze-out curve are given depending on (left panel) and center-of-mass energy (right panel). It seems that the evolution of fermionic and bosonic skewness and kurtosis coincide at one point, marked with the vertical band. This fact would reflect the nature of the phase transition. It would be the critical endpoint connecting cross-over with the first order deconfinement phase transitions. At QCD CEP, the phase transition is conjectured to be of second order. It is worthwhile to mention that the crossing point is amazingly coincident with the QCD CEP measured in lattice QCD [5, 6], Figure 14, regardless of its uncertainties.

From (23) and (24), the following expressions have to be solved in , individually and/or dependently, in order to determine of the crossing point: Due to the mathematical difficulties in dealing with these expressions, the integral over phase space has to be simplified. A suitable simplification is given in Section 2.3.1. where and being baryon (strange) quantum number and is the baryon (strange) chemical potential of light and strange quarks, respectively. The quarks chemistry is introduced in Section 2. Accordingly, the difference between baryons and fermions is originated in the exponential function. For simplicity, we consider one fermion and one boson particle. Then, the baryon chemical potential at the chemical freeze-out curve at which the fermionic and bosonic skewness (or kurtosis) curves of these two particles cross with each other can be given as In the relativistic limit, while in the nonrelativistic limit . It is obvious that the bosonic and fermionic degrees of freedom play an essential role in determining (42). Furthermore, it seems that the chemical potential of strange quark has no effect at the crossing point.

The dependence of on as calculated in HRG is given in Figure 16. As mentioned above, is calculated to guarantee strange number conservation in heavy-ion collisions. At small , has a linear dependence, (Hooke’s limit). At large , the dependence is no longer linear.

As discussed earlier, the critical behavior and the existence of QCD CEP can be identified by means of signatures sensitive to singular parts of the free energy, especially the ones reflecting dynamical fluctuations of conserved charges, such as baryon number and charge density [76]. The reason that the fermionic and bosonic higher order moments are crossing on a point very near the QCD CEP shall be elaborated in a forthcoming work [77].

#### 5. Conclusion

In the present work, the first six nonnormalized order moments of the particle multiplicity are calculated in the hadron resonance gas model. General expressions for arbitrary higher order moments are deduced in terms of lower ones. We concluded that going from lower to higher order moments is possible through adding up a series that consists of all lower order moments plus correlation functions. We studied the thermal evolution of the first four normalized order moments and their products (ratios) at different chemical potentials . By doing that, it was possible to evaluate first four normalized moments at the chemical freeze-out curve. The freeze-out curve is characterized by constant at all values of , where and are entropy density and temperature, respectively. It has been found that nonmonotonic behavior reflecting dynamical fluctuation and strong correlations appears starting from the normalized third order moment (skewness ). Furthermore, nonmonotonicity is observed in the normalized fourth order moment, the kurtosis , and its products. These are novel observations. Although HRG is exclusively applicable below , that is, it does not include deconfinement phase transition, it is apparent that the higher order moments are able to give signatures for the critical phase transition. The accuracy of our calculations made it possible to have clear evidence on the phase transition. Based on these findings, we introduced novel conditions characterizing the chemical freeze-out curve and the location of the QCD critical endpoint as follows. The chemical freeze-out curve is described by , where is the susceptibility in particle number, that is, the second order moment. The location of QCD critical endpoint (at and axes) is positioned when the condition or is fulfilled. The subscripts and refer to bosons and fermions, respectively. Accordingly, the QCD endpoint is positioned at MeV and MeV. We are able to estimate these two quantities, although the hadron resonance gas model basically does not contain information on criticality related with the chiral dynamics and singularity in physical observables required to locate the critical endpoint. After submitting this work, a new paper was posted in arXiv, in which the authors used the second order moment (susceptibility) and the fourth order one (kurtosis) that are apparently sensitive to the phase transition [78].

#### Acknowledgment

This work is partly supported by the German-Egyptian Scientific Projects (GESP ID: 1378).

#### References

- C. Athanasiou, K. Rajagopal, and M. Stephanov, “Using higher moments of fluctuations and their ratios in the search for the QCD critical point,”
*Physical Review D*, vol. 82, no. 7, Article ID 074008, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. Adams, M. M. Aggarwal, Z. Ahammed et al., “Experimental and theoretical challenges in the search for the quark gluon
plasma: the STAR Collaboration’s critical assessment of the evidence from RHIC collisions,”
*Nuclear Physics A*, vol. 757, p. 102, 2005. View at Google Scholar - A. Tawfik, “QCD phase diagram: a comparison of lattice and hadron resonance gas model calculations,”
*Physical Review D*, vol. 71, Article ID 054502, 2005. View at Google Scholar - M. A. Stephanov, “QCD phase diagram and the critical point,”
*International Journal of Modern Physics A*, vol. 20, p. 4387, 2005. View at Google Scholar - Z. Fodor and S. D. Katz, “Critical point of QCD at finite T and
*μ*, lattice results for physical quark masses,”*Journal of High Energy Physics*, vol. 0404, p. 50, 2004. View at Google Scholar - R. V. Gavai and S. Gupta, “QCD at finite chemical potential with six time slices,”
*Physical Review D*, vol. 78, no. 11, Article ID 114503, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. Tawfik, “Phase space and dynamical fluctuations of kaon-to-pion ratios,”
*Progress of Theoretical Physics*, vol. 126, no. 2, pp. 279–292, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. Tawfik, “Matter—antimatter asymmetry in heavy-ion collisions,”
*International Journal of Theoretical Physics*, vol. 51, pp. 1396–1407, 2012. View at Publisher · View at Google Scholar - A. Tawfik, “Antiproton-to-proton ratios for ALICE heavy-ion collisions,”
*Nuclear Physics A*, vol. 859, no. 1, pp. 63–72, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. Tawfik, “Event-by-event fluctuations of particle ratios in heavy-ion collisions,”
*FIZIKA B*, vol. 18, pp. 141–150, 2009. View at Google Scholar - A. Tawfik, “Dynamical fluctuations in baryon-meson ratios,”
*Journal of Physics G*, vol. 40, Article ID 055109, 2013. View at Google Scholar - A. Tawfik, “Fluctuations of particle yield ratios in heavy-ion collisions,”
*Indian Journal of Physics*, vol. 86, pp. 1139–1146, 2012. View at Google Scholar - A. Tawfik, “Black box QGP,” http://arxiv.org/abs/hep-ph/0604037.
- M. I. Adamovich, M. M. Aggarwal, N. P. Andreeva et al., “Factorial moments of Si-28 induced interactions with Ag(Br) nuclei,”
*Acta Physica Hungarica A*, vol. 13, no. 4, pp. 213–221, 2001. View at Publisher · View at Google Scholar - X. Luo, “Search for the QCD critical point by higher moments of net-proton multiplicity distributions at RHIC,”
*Acta Physica Polonica B*, vol. 5, no. 2, p. 497, 2012. View at Publisher · View at Google Scholar - X. Luo, “Probing the QCD critical point with higher moments of net-proton multiplicity distributions,”
*Journal of Physics: Conference Series*, vol. 316, Article ID 012003, 2011. View at Publisher · View at Google Scholar - M. M. Aggarwal, Z. Ahammed, A. V. Alakhverdyants et al., “Higher moments of net proton multiplicity distributions at RHIC,”
*Physical Review Letters*, vol. 105, Article ID 022302, 2010. View at Publisher · View at Google Scholar - T. J. Tarnowsky, “Searching for the QCD critical point using particle ratio fluctuations and higher moments of multiplicity distributions,”
*Journal of Physics G*, vol. 38, Article ID 124054, 2011. View at Publisher · View at Google Scholar - T. K. Nayak, “Study of the fluctuations of net-charge and net-protons using higher order moments,”
*Nuclear Physics A*, vol. 830, pp. 555C–558C, 2009. View at Publisher · View at Google Scholar - M. A. York and G. D. Moore, “Exploring the phase diagram with Taylor series: epic voyage or just another bad trip,” http://arxiv.org/abs/1106.2535.
- M. A. Stephanov, “Non-Gaussian fluctuations near the QCD critical point,”
*Physical Review Letters*, vol. 102, no. 3, Article ID 032301, 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. Zhou, S. S. Shi, K. Xiao, K. J. Wu, and F. Liu, “Higher moments of net-baryon distribution as probes of QCD critical point,”
*Physical Review C*, vol. 82, Article ID 014905, 2010. View at Publisher · View at Google Scholar - M. Stephanov, K. Rajagopal, and E. Shuryak, “Event-by-event fluctuations in heavy ion collisions and the QCD critical point,”
*Physical Review D*, vol. 60, no. 11, Article ID 114028, 32 pages, 1999. View at Publisher · View at Google Scholar · View at Scopus - Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz, and K. K. Szabó, “The order of the quantum chromodynamics transition predicted by the standard model of particle physics,”
*Nature*, vol. 443, no. 7112, pp. 675–678, 2006. View at Publisher · View at Google Scholar · View at Scopus - M. Cheng, N. H. Christ, S. Datta et al., “The transition temperature in QCD,”
*Physical Review D*, vol. 74, Article ID 054507, 2006. View at Publisher · View at Google Scholar - M. Cheng, P. Hegde, C. Jung et al., “Baryon number, strangeness and electric charge fluctuations in QCD at high temperature,”
*Physical Review D*, vol. 79, Article ID 074505, 2009. View at Publisher · View at Google Scholar - F. Karsch, K. Redlich, and A. Tawfik, “Hadron resonance mass spectrum and lattice QCD thermodynamics,”
*European Physical Journal C*, vol. 29, no. 4, pp. 549–556, 2003. View at Publisher · View at Google Scholar · View at Scopus - F. Karsch, K. Redlich, and A. Tawfik, “Thermodynamics at non-zero baryon number density: a comparison of lattice and hadron resonance gas model calculations,”
*Physics Letters B*, vol. 571, no. 1-2, pp. 67–74, 2003. View at Publisher · View at Google Scholar · View at Scopus - K. Redlich, F. Karsch, and A. Tawfik, “Heavy-ion collisions and lattice QCD at finite baryon density,”
*Journal of Physics G*, vol. 30, no. 8, pp. S1271–S1274, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, “Comparison of chemical freeze-out criteria in heavy-ion collisions,”
*Physical Review C*, vol. 73, no. 3, Article ID 034905, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Ejiri, “Canonical partition function and finite density phase transition in lattice QCD,”
*Physical Review D*, vol. 78, Article ID 074507, 2008. View at Publisher · View at Google Scholar - E. S. Bowman and J. I. Kapusta, “Critical points in the linear
*σ*model with quarks,”*Physical Review D*, vol. 79, Article ID 015202, 2009. View at Publisher · View at Google Scholar - B. I. Abelev, M. M. Aggarwal, Z. Ahammed et al., “Identified particle production, azimuthal anisotropy, and interferometry
measurements in Au + Au collisions at $\sqrt{{S}_{NN}}$ = 9.2-GeV,”
*Physical Review C*, vol. 81, Article ID 024911, 2010. View at Publisher · View at Google Scholar - S.-X. Qin, L. Chang, H. Chen, Y.-X. Liu, and C. D. Roberts, “Phase diagram and critical end point for strongly interacting quarks,”
*Physical Review Letters*, vol. 106, Article ID 172301, 2011. View at Publisher · View at Google Scholar - A. Tawfik, “Influence of strange quarks on the QCD phase diagram and chemical freeze-out,”
*Journal of Physics G*, vol. 31, no. 6, pp. S1105–S1110, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. Tawfik, “In-medium modifications of hadron properties,”
*Indian Journal of Physics*, vol. 85, no. 5, pp. 755–766, 2011. View at Google Scholar · View at Scopus - A. Tawfik, “A universal description for the freezeout parameters in heavy-ion collisions,”
*Nuclear Physics A*, vol. 764, no. 1–4, pp. 387–392, 2006. View at Publisher · View at Google Scholar · View at Scopus - A. Tawfik, “Condition driving chemical freeze-out,”
*Europhysics Letters*, vol. 75, no. 3, pp. 420–426, 2006. View at Publisher · View at Google Scholar · View at Scopus - R. Hagedorn, “Statistical thermodynamics of strong interactions at high-energies,”
*Nuovo Cimento, Supplemento*, vol. 3, pp. 147–186, 1965. View at Google Scholar - R. Venugopalan and M. Prakash, “Thermal properties of interacting hadrons,”
*Nuclear Physics A*, vol. 546, no. 4, pp. 718–760, 1992. View at Google Scholar · View at Scopus - M. Gell-Mann and M. Lévy, “The axial vector current in beta decay,”
*Il Nuovo Cimento*, vol. 16, no. 4, pp. 705–726, 1960. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - R. Hagedorn, “Statistical thermodynamics of strong interactions at high energies. 3. Heavy-pair (quark) production rates,”
*Nuovo Cimento, Supplemento*, vol. 6, pp. 311–354, 1968. View at Google Scholar - R. Hagedorn, “Hadronic matter near the boiling point,”
*Il Nuovo Cimento A*, vol. 56, no. 4, pp. 1027–1057, 1968. View at Publisher · View at Google Scholar - R. J. Eden, P. V. Landshoff, D. I. Olive, and J. C. Polkinghorne,
*The Analytic S-Matrix*, Cambridge University Press, London, UK, 1966. View at MathSciNet - J. Letessier and J. Rafelski,
*Hadrons and Quark-Gluon Plasma*, Cambridge University Press, 2002. - V. Gupta, “Predictions for the higher moments of the charged-particle multiplicity distribution,”
*Lettere al Nuovo Cimento*, vol. 5, no. 1, pp. 108–111, 1972. View at Publisher · View at Google Scholar · View at Scopus - B.-J. Schaefer and M. Wagner, “QCD critical region and higher moments for three flavor models,”
*Physical Review D*, Article ID 034027, 2012. View at Google Scholar - C. Miao, “Higher moments of charge fluctuations in QCD at high temperature,”
*Nuclear Physics A*, vol. 830, pp. 705C–708C, 2009. View at Google Scholar - P. Hegde, “Fluctuations and higher moments of conserved charges from the lattice,” in
*Proceedings of the 19th Particles and Nuclei International Conference (PANIC '11)*, vol. 1441 of*AIP Conference Proceedings*, pp. 898–900, July 2011. View at Publisher · View at Google Scholar · View at Scopus - K. P. Ballanda and H. L. MacGillivray, “Kurtosis: a critical review,”
*American Statistician*, vol. 42, pp. 111–119, 1988. View at Google Scholar - R. Singer, Y. Cho, T. Fields et al., “Two-particle rapidity correlations in 205 GeV/
*c*pp interactions,”*Physics Letters B*, vol. 49, p. 481, 1974. View at Publisher · View at Google Scholar - U. Heinz and B. V. Jacak, “Two-particle correlations in relativistic heavy-ion collisions,”
*Annual Review of Nuclear and Particle Science*, vol. 49, pp. 529–579, 1999. View at Publisher · View at Google Scholar - S. Chapman, J. R. Nix, and U. Heinz, “Extracting source parameters from Gaussian fits to two-particle correlations,”
*Physical Review C*, vol. 52, no. 5, pp. 2694–2703, 1995. View at Publisher · View at Google Scholar · View at Scopus - E. Braidot, “Two-particle correlations at forward rapidity in STAR,”
*Nuclear Physics A*, vol. 854, pp. 168–174, 2011. View at Publisher · View at Google Scholar - C. Adler, Z. Ahammed, C. Allgower et al., “Elliptic flow from two and four particle correlations in Au + Au collisions at $\sqrt{{S}_{\text{NN}}}$ = 130-GeV,”
*Physical Review C*, vol. 66, Article ID 034904, 2002. View at Publisher · View at Google Scholar - G. A. Kozlov, “New features of two-particle correlations,”
*Physics of Particles and Nuclei*, vol. 36, no. 2, pp. S190–S197, 2005. View at Google Scholar · View at Scopus - R. Perez-Ramos, “Single inclusive distribution and two-particle correlations inside one jet at “modified leading logarithmic approximation” of quantum chromodynamics II. Steepest descent evaluation at small x,”
*Journal of High Energy Physics*, vol. 2006, no. 9, article 014, 2006. View at Publisher · View at Google Scholar · View at Scopus - E. A. De Wolf, I. M. Dremin, and W. Kittel, “Scaling laws for density correlations and fluctuations in multiparticle dynamics,”
*Physics Report*, vol. 270, no. 3, pp. 1–141, 1996. View at Google Scholar · View at Scopus - W. Kittel and E. A. De Wolf,
*Soft Multihadron Dynamics*, World Scientific, Singapore, 2005. - N. L. Johnson, S. Kotz, and N. Balakkrishan,
*Continuous Univariate Distribution*, Houghton Mifflin, Boston, MA, 1970. - C. T. Hsu and D. N. Lawley, “The derivation of the fifth and sixth moments of the distribution of ${b}_{2}$ in samples from a normal population,”
*Biometrika*, vol. 31, pp. 238–248, 1940. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Asakawa, S. Ejiri, and M. Kitazawa, “Third moments of conserved charges as probes of QCD phase structure,”
*Physical Review Letters*, vol. 103, no. 26, Article ID 262301, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. Kitazawa, M. Asakawa, and S. Ejiri, “Third moments of conserved charges in QCD phase diagram,”
*PoS LAT2009*, vol. 174, 2009. View at Publisher · View at Google Scholar - B. Stokic, B. Friman, and K. Redlich, “Kurtosis and compressibility near the chiral phase transition,”
*Physics Letters B*, vol. 673, pp. 192–196, 2009. View at Publisher · View at Google Scholar - R. V. Gavai and S. Gupta, “Lattice QCD predictions for shapes of event distributions along the freezeout curve in heavy-ion collisions,”
*Physics Letters B*, vol. 696, no. 5, pp. 459–463, 2011. View at Publisher · View at Google Scholar · View at Scopus - F. Katsch and K. Redlich, “Probing freeze-out conditions in heavy ion collisions with moments of charge fluctuations,”
*Physics Letters B*, vol. 695, pp. 136–142, 2010. View at Google Scholar - B. I. Abelev, M. M. Aggarwal, Z. Ahammed et al., “Beam-energy and system-size dependence of dynamical net charge fluctuations,”
*Physical Review C*, vol. 79, Article ID 024906, 2009. View at Publisher · View at Google Scholar - X. Luo, B. Mohanty, H. G. Ritter, and N. Xu, “Higher moments of net-proton multiplicity distributions,”
*Physics of Atomic Nuclei*, vol. 75, no. 6, pp. 676–678, 2012. View at Publisher · View at Google Scholar - T. Schuster, “Net-baryon-, net-proton-, and net-charge kurtosis in heavy-ion collisions within a relativistic transport approach,”
*The European Physical Journal C*, vol. 72, p. 2143, 2012. View at Google Scholar - J. Cleymans, M. Stankiewicz, P. Steinberg, and S. Wheaton, “The origin of the difference between multiplicities in ${e}^{+}{e}^{-}$ annihilation and heavy ion collisions,” submitted to
*Physics Letters B*. - B. Friman, F. Karsch, K. Redlich, and V. Skokov, “Fluctuations as probe of the QCD phase transition and freeze-out in heavy ion collisions at LHC and RHIC,”
*The European Physical Journal C*, vol. 71, p. 1694, 2011. View at Google Scholar - R. V. Gavai and S. Gupta, “Lattice QCD predictions for shapes of event distributions along the freezeout curve in heavy-ion collisions,”
*Physics Letters B*, vol. 696, pp. 459–463, 2011. View at Publisher · View at Google Scholar - F. Karsch, “Moments of charge fluctuations, pseudo-critical temperatures and freeze-out in heavy-ion collisions,”
*Journal of Physics G*, vol. 38, no. 12, Article ID 124098, 2011. View at Publisher · View at Google Scholar · View at Scopus - F. Karsch, B.-J. Schaefer, M. Wagner, and J. Wambach, “Towards finite density QCD with Taylor expansions,”
*Physics Letters B*, vol. 698, no. 3, pp. 256–264, 2011. View at Publisher · View at Google Scholar · View at Scopus - B. Berdnikov and K. Rajagopal, “Slowing out of equilibrium near the QCD critical point,”
*Physical Review D*, vol. 61, Article ID 105017, 9 pages, 2000. View at Publisher · View at Google Scholar - E. Nakano, B.-J. Schaefer, B. Stokic, B. Friman, and K. Redlich, “Fluctuations and isentropes near the chiral critical endpoint,”
*Physics Letters B*, vol. 682, no. 4-5, pp. 401–407, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Tawfik, “Chemical freeze-out and higher order multiplicity moments,” http://arxiv.org/abs/1306.1025.
- A. Nakamura and K. Nagata, “Probing QCD phase structure by baryon multiplicity distribution,” http://arxiv.org/abs/1305.0760.