# Chemical Potentials of Light Flavor Quarks from Yield Ratios of Negative to Positive Particles in Au+Au Collisions at RHIC

Authors: Ya-Qin Gao; Hai-Ling Lao; Fu-Hu Liu
Journal: Advances in High Energy Physics
 (2018)

---

## Abstract

The transverse momentum spectra ofπ-, π+, K-, K+, p¯, and p produced in Au+Au collisions at center-of-mass energy sNN=7.7, 11.5, 19.6, 27, 39, 62.4, 130, and 200 GeV are analyzed in the framework of a multisource thermal model. The experimental data measured at midrapidity by the STAR Collaboration are fitted by the (two-component) standard distribution. The effective temperature of emission source increases obviously with the increase of the particle mass and the collision energy. At different collision energies, the chemical potentials of up, down, and strange quarks are obtained from the antiparticle to particle yield ratios in given transverse momentum ranges available in experiments. With the increase of logarithmic collision energy, the chemical potentials of light flavor quarks decrease exponentially.


---

## Body

## 1. Introduction

The constructions of the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) have been opening a new epoch for the studies of nuclear and quark matters. One of the major goals of the RHIC and LHC studies is to obtain information on the quantum chromodynamics (QCD) phase diagram [1]. The phase diagram includes at least a fundamental phase transition between the hadron gas and the quark-gluon plasma (QGP) or quark matter and is usually plotted as chemical freeze-out temperature (Tch) versus baryon chemical potential (μbaryon). Nowadays, the detailed characteristics of the phase diagram are not known yet. The experimental and theoretical nuclear physicists have been focusing their attentions on the searching for the critical end point and phase boundary. Lattice QCD calculations show that a system is produced at small μbaryon or high energies through a crossover at the quark-hadron phase transition [2–4]. Based on the lattice QCD [5] and several QCD-based models calculations [6–9], as well as mathematical extensions of lattice techniques [10–13], researchers suggest that the transition at larger μbaryon is the first order and the QCD critical end point is existent.Pinpointing the phase boundary and the critical end point is the central issue to understand the properties of interacting matter under extreme conditions and to map the QCD phase diagram. The matter produced in high-energy heavy-ion collisions provides the opportunity to search for the phase boundary and the critical end point [6, 14]. To this end, the STAR Collaboration at the RHIC has undertaken the first phase of the beam energy scan (BES) program [15–17], starting the second phase from 2018 to 2019 [18]. The program is to vary the collision energy which enables a search for nonmonotonic excitation functions over a broad domain of the phase diagram. Before looking for an evidence for the existence of a critical end point and the phase boundary, it is important to know the (Tch,μbaryon) region of phase diagram one can access. The produced particles spectra and yield ratios allow us only to infer the values of Tch and μbaryon [19]. Furthermore, the bulk properties such as rapidity density dN/dy, mean transverse momentum 〈pT〉, particle ratios, and freeze-out properties may provide an insight into the particle production mechanisms at BES energies. Therefore, it is very important to study these bulk properties systematically, which may reveal the evolution and the changes of the system created in high-energy heavy-ion collisions.As one of the most important measured quantities, the transverse momentum (pT) spectrum includes abundant information which is related to the excitation degree of the collision system. The spectra of identified particles can also provide useful information about temperature, particle ratio, and chemical potential by using thermal and statistical investigations [20]. For any system, one can determine the direction and limitation of mass transfer by comparing the chemical potentials of particles; that is to say, the chemical potential is a sign to mark the direction of spontaneous chemical reaction. The chemical potential can also be a criterion for determining whether thermodynamic equilibrium does exist in the interacting region in high-energy collisions [1]. Generally, a low absolute value of chemical potential corresponds to a high degree of thermodynamic equilibrium. Therefore, the chemical potential is also one of the major solutions for investigating the QGP. One can see that the chemical potentials of quarks are an important subject at high energy. Therefore, we are very interested in measuring the chemical potentials of quarks.In this paper, we extract the chemical potentials of light flavor quarks from the yield ratios of negatively to positively charged particles. By using the (two-component) standard distribution, thepT spectra of π-, π+, K-, K+, p¯, and p produced in Au+Au collisions at center-of-mass energy (per nucleon pair) sNN=7.7, 11.5, 19.6, 27, 39, 62.4, 130, and 200 GeV measured by the STAR Collaboration in midrapidity interval (|y|<0.1) [19, 21] are described. The considered energies stretch across a wide energy range which covers the main range of the RHIC at its BES.


## 2. The Model and Method

To extract the chemical potentials of quarks, we need to know the yield ratios of negatively to positively charged particles. Although we can have the values of yield ratios directly in experiments, they are not complete and comprehensive in some cases. Usually, thepT spectra of charged particles are given in many experiments and we can get the yield ratios by fitting the available data. Then, the values of chemical potentials for the up, down, and strange quarks can be obtained from the yield ratios π-/π+, K-/K+, and p¯/p which are synthetically considered in special ways.In this paper, thepT spectra are analyzed in the framework of a multisource thermal model [22], which assumes that various sources are involved in high-energy collisions. These sources are divided into few groups by different interaction mechanisms, geometrical relations, or event samples. Each group of sources forms a relatively large emission source which stays in a local thermal equilibrium state at the chemical or kinetic freeze-out. Each emission source is considered to emit particles in its rest frame and treated as a thermodynamic system of relativistic and quantum ideal gas. This means that each emission source can be described by the thermal and statistical model or other similar models and distributions. The final-state distribution is attributed to all sources in the whole system, which results in a multicharacteristic emission process [22] if we use the standard distribution [23–26]. This also means that pT spectrum can be described by a multicomponent standard distribution in which each component describes a given emission source.We now structure the multicomponent standard distribution. It is assumed that there arel components to be considered. For the i-th component, the standard Boltzmann, Fermi-Dirac, and Bose-Einstein distributions [23–26] can be uniformly expressed as(1)fipT=1NdNdpT=CipTpT2+m02∫yminymaxcoshy×exppT2+m02coshy-μTi+S-1dy,where Ci is the normalization constant which results in ∫0∞‍fi(pT)dpT=1; N, m0, μ, and Ti denote the particle number, the rest mass of the considered particle, the chemical potential of the considered particle, and the effective temperature for the i-th component, respectively; ymin is the minimum rapidity and ymax is the maximum rapidity; the values of S are 0, +1, and -1, which denote the Boltzmann, Fermi-Dirac, and Bose-Einstein distributions, respectively. We neglect the existence of μ in (1) due to the fact that it has mainly effect on the normalization which can be redone, but not the trend of curve.In the final state,pT spectrum is resulting from l components; that is,(2)fpT=1NdNdpT=∑i=1lwifipT,where wi (i=1,2,…,l) is the relative weight resulting from the i-th component. Because of the probability distribution being acquiescently normalized to 1, the coefficient obeys the normalization condition of ∑‍wi=1. Considering the relative contribution of each component, we have the mean effective temperature to be Teff=∑‍wiTi, which reflects the mean excitation degree of different sources corresponding to different components and can be used to describe the effective temperature of whole interacting system. It should be noted that the effective temperature contains the contributions of transverse flow and thermal motion. It is not the “real" temperature of the interacting system.According to [27, 28], the relation between antiproton to proton yield ratios can be written as(3)p¯p=exp-2μpTch≈exp-2μbaryonTch,where μp denotes the chemical potential of proton. In the framework of the statistical thermal model of noninteracting gas particles with the assumption of standard Maxwell-Boltzmann statistics, there is an empirical expression for Tch [29–32]; one has(4)Tch=Tlim11+exp2.60-lnsNN/0.45,where sNN is in the units of GeV and the “limiting" temperature Tlim=0.164 GeV [29, 30].In a similar way, the yield ratios of antiparticles to particles for other hadrons can be written as(5)kπ≡π-π+=exp-2μπTch,kK≡K-K+=exp-2μKTch,kp≡p¯p=exp-2μpTch,kD≡D-D+=exp-2μDTch,kB≡B-B+=exp-2μBTch,where kj (j=π, K, p, D, and B) denote the yield ratios of negatively to positively charged particles obtained from the normalization constants of pT spectra. The symbols μπ, μK, μD, and μB represent the chemical potentials of π, K, D, and B, respectively. In the above discussion, the symbol of a given particle is used for its yield for the purpose of simplicity. Furthermore, we have(6)μπ=-12Tch·lnkπ,μK=-12Tch·lnkK,μp=-12Tch·lnkp,μD=-12Tch·lnkD,μB=-12Tch·lnkB.Letμq denote the chemical potential for quark flavor, where q=u, d, s, c, and b represent the up, down, strange, charm, and bottom quarks, respectively. In principle, we can use kj to give relations among different μq. The values of μq are then expected from these relations. According to [33, 34], based on the same Tch, kj in terms of μq are(7)kπ=exp-μu-μd/Tchexpμu-μd/Tch=exp-2μu-μdTch,kK=exp-μu-μs/Tchexpμu-μs/Tch=exp-2μu-μsTch,kp=exp-2μu+μd/Tchexp2μu+μd/Tch=exp-22μu+μdTch,kD=exp-μc-μd/Tchexpμc-μd/Tch=exp-2μc-μdTch,kB=exp-μu-μb/Tchexpμu-μb/Tch=exp-2μu-μbTch.Thus, we have(8)μu=-16Tch·lnkπ·kp,μd=-16Tch·lnkπ-2·kp,μs=-16Tch·lnkπ·kK-3·kp,μc=-16Tch·lnkπ-2·kp·kD3,μb=-16Tch·lnkπ·kp·kB-3.As can be seen from (8), μq are obtained from kj. In addition to the yield ratios π-/π+, K-/K+, and p¯/p, other combinations can also give μq if the spectra in the numerator and denominator are under the same experimental conditions.


## 3. Results and Discussion

The energy dependent double-differentialpT spectra of π-, π+, K-, K+, p¯, and p produced in central Au+Au collisions at sNN=7.7, 11.5, 19.6, 27, 39, 62.4, 130, and 200 GeV at the midrapidity |y|<0.1 are presented in Figure 1, where the centrality interval at 130 GeV is 0–6% and at other energies is 0–5%. The different symbols represent the data measured by the STAR Collaboration [19, 21], and the curves are the results fitted here by the (two-component) standard distribution. Generally, the standard distribution is firstly used in the fit process. If it does not fit the data, the two-component standard distribution is used. It is because of the quality of the measurements that (two-component) standard distribution is used. In the case of using the two-component standard distribution, the first component results in narrow pT region and the second component results in wide pT regions. That is, in low pT region both components contribute to the spectra, and in high pT region only the second component contributes to the spectra. In the calculation, the values of the free parameters (T1, w1, and T2), the normalization constant (N0), and χ2 obtained by fitting the data are listed in Table 1 including the degrees of freedom (dof). One can see that the data are well fitted by the (two-component) standard distribution. From the parameter values, one can see that the effective temperature increases with the increase of the particle mass and the collision energy for emissions of the six types of particles.Table 1
Values ofT1, w1, T2, N0, χ2, and dof corresponding to the curves in Figure 1. The values for positively and negatively charged particles are given using the slash (/), where the values for positive particles are shown before the slash and the values for negative particles after the slash. The values of N0 are obtained due to the comparisons between the experimental (1/2πpT)d2N/dpTdy and the calculated (1/2πpT)N0f(pT)/dy, where dy=0.2 and f(pT) is presented in (2).
s N N (GeV) Particles T 1 (MeV) w 1 T 2 (MeV) N 0 χ 2 dof 7.7 π ± 87 ± 7 / 89 ± 10 0.65 ± 0.04 / 0.63 ± 0.07 181 ± 2 / 177 ± 4 81.0 ± 4.0 / 81.0 ± 9.0 0.38/0.45 22/22 K ± 185 ± 6 / 170 ± 8 1.00/1.00 - 9.5 ± 0.8 / 3.8 ± 0.5 0.75/0.98 21/21 p/p- 224 ± 7 / 257 ± 20 1.00/1.00 - 19.0 ± 0.2 / 0.13 ± 0.01 1.15/0.52 27/13  11.5 π ± 102 ± 10 / 105 ± 10 0.70 ± 0.04 / 0.69 ± 0.03 193 ± 4 / 190 ± 4 92.0 ± 8.0 / 93.0 ± 8.0 0.25/0.48 22/22 K ± 190 ± 5 / 178 ± 4 1.00/1.00 - 11.0 ± 0.8 / 5.7 ± 0.5 0.49/0.97 23/21 p/p- 217 ± 10 / 216 ± 20 1.00/1.00 - 16.0 ± 3.0 / 0.54 ± 0.10 1.70/3.07 26/21  19.6 π ± 118 ± 10 / 116 ± 10 0.79 ± 0.02 / 0.83 ± 0.02 215 ± 8 / 219 ± 7 109.0 ± 10.0 / 119.0 ± 10.0 0.19/0.40 22/22 K ± 181 ± 7 / 181 ± 5 0.88 ± 0.04 / 0.89 ± 0.05 260 ± 20 / 239 ± 20 13.0 ± 1.4 / 8.7 ± 0.8 0.55/1.61 22/22 p/p- 234 ± 20 / 237 ± 10 1.00/1.00 - 11.5 ± 2.0 / 1.4 ± 0.1 0.69/19.81 27/20  27 π ± 117 ± 8 / 118 ± 7 0.79 ± 0.03 / 0.81 ± 0.03 219 ± 7 / 221 ± 7 120.0 ± 10.0 / 125.0 ± 10.0 0.19/0.27 22/22 K ± 178 ± 7 / 180 ± 6 0.85 ± 0.03 / 0.84 ± 0.06 262 ± 20 / 236 ± 10 14.0 ± 1.0 / 10.0 ± 1.0 0.27/0.92 22/21 p/p- 239 ± 7 / 247 ± 10 1.00/1.00 - 10.0 ± 2.0 / 2.0 ± 0.2 0.84/1.28 21/20  39 π ± 114 ± 10 / 116 ± 10 0.78 ± 0.02 / 0.78 ± 0.02 222 ± 7 / 222 ± 6 129.0 ± 12.0 / 128.0 ± 10.0 0.39/0.20 22/22 K ± 189 ± 7 / 189 ± 5 0.95 ± 0.01 / 0.96 ± 0.02 332 ± 18 / 347 ± 23 14.0 ± 1.5 / 11.0 ± 1.0 0.20/0.36 22/22 p/p- 250 ± 20 / 252 ± 20 1.00/1.00 - 8.3 ± 1.0 / 2.8 ± 0.3 0.81/2.42 20/21  62.4 π ± 139 ± 6 / 137 ± 8 1.00/1.00 - 146.0 ± 10.0 / 150.0 ± 10.0 1.14/1.07 8/8 K ± 210 ± 25 / 214 ± 25 1.00/1.00 - 15.8 ± 1.0 / 13.6 ± 0.8 0.02/0.37 8/8 p/p- 335 ± 20 / 346 ± 40 1.00/1.00 - 8.3 ± 0.4 / 3.8 ± 0.3 1.97/1.60 13/14  130 π ± 136 ± 8 / 137 ± 7 1.00/1.00 - 181.0 ± 13.0 / 185.0 ± 11.0 1.98/3.68 8/8 K ± 204 ± 12 / 210 ± 13 1.00/1.00 - 18.5 ± 1.4 / 17.3 ± 1.3 0.39/0.20 11/11 p/p- 373 ± 15 / 385 ± 22 1.00/1.00 - 7.5 ± 0.7 / 5.5 ± 0.3 1.73/0.47 11/11  200 π ± 116 ± 8 / 115 ± 5 0.76 ± 0.03 / 0.76 ± 0.03 263 ± 23 / 262 ± 25 209.0 ± 14.0 / 215.0 ± 12.0 0.21/0.17 7/7 K ± 235 ± 30 / 239 ± 30 1.00/1.00 - 19.3 ± 1.1 / 18.5 ± 1.0 0.06/0.06 8/8 p/p- 382 ± 40 / 393 ± 35 1.00/1.00 - 8.3 ± 0.9 / 6.4 ± 0.6 0.10/0.25 14/15Figure 1
Midrapidity (|y|<0.1) double-differential pT spectra for π-, π+, K-, K+, p-, and p in central Au+Au collisions at sNN=7.7, 11.5, 19.6, 27, 39, 62.4, 130, and 200 GeV, where the centrality interval at 130 GeV is 0–6% and at other energies is 0–5%. The different symbols represent the measurements done by the STAR experiment [19, 21] and the curves represent the results fitted by the (two-component) standard distribution. The values of parameters can be found in Table 1.Based on the above successful fits of thepT spectra of antiparticles and particles, we can use (8) and the pT spectra in Figure 1 to study the dependence of μq on sNN. This is done by integrating the yield over the given pT ranges available in experiments at different energies. Figure 2 shows the correlations between μq and sNN, where μ is used on the vertical axis to replace μq which are marked in the panel for different styles of symbols. With the increase of logarithmic sNN, an exponential decrease of μq is observed. Corresponding to the solid, dashed, and dotted curves which fit to the dependences of μu, μd, and μs on sNN, respectively, we have(9)μu=820.1±0.1sNNGeV-0.914±0.025MeV,μd=681.1±0.1sNNGeV-0.834±0.031MeV,μs=420.7±0.2sNNGeV-1.004±0.063MeVwith the χ2/dof=1.67/6, 2.73/6, and 1.04/6, respectively.Figure 2
Correlations betweenμq and sNN for central Au+Au collisions at RHIC. The symbols represent μq obtained from the ratios by integrating the yield over the given pT ranges available in experiments in Figure 1. The solid, dashed, and dotted curves are fitted results corresponding to μu, μd, and μs, respectively.The similarity in up and down quark masses renders the similarity in their chemical potentials. The difference between the chemical potentials of up (or down) and strange quarks is caused by the difference between their masses. At the lowest BES energy the difference between the chemical potentials is dozens of MeV, while at the highest RHIC energy these quantities are around a few MeV. The decrease inμq is obvious, which indicates the change of mean free path of produced quarks in the middle state. If the produced quarks at the lowest BES energy have a small mean free path which looks as if a liquid-like middle state is formed, the produced quarks at the highest RHIC energy should have a large mean free path which looks as if a gas-like middle state is formed. The main difference at different energies is different mean free paths of the produced quarks. To search for the critical energy at which the change from a liquid-like middle state to a gas-like middle state had happenned is beyond the focus of the present work.From (9) we can obtain a linear relation between lnμq and lnsNN,(10)lnμq=a-b·lnsNN,where the intercept a and slope -|b| can be obtained from the parameters in (9). In particular, -|b| is close to -1. The large negative slope shows an obvious anticorrelation between lnμq and lnsNN. It is expected that lnμq will be smaller at higher energy or larger at lower energy. In particular, at the LHC energies, lnμq will be negative since μq will be less than 1 MeV. The limiting value of μq is close to 0 at the LHC [35], which results in an obvious negative lnμq.The main conclusion observed from Figure2 is that μq is high (from dozens of MeV to ~100 MeV) at the BES and close to 0 at the LHC [35]. This is consistent with the trend of μbaryon (~100–300 MeV at the BES and ~1 MeV at the LHC) obtained from other works [1, 29–32, 36]. This is natural due to the fact that baryon consists of valence quarks. If we regard μbaryon=Σμq, where Σ denotes the sum over all valence quarks in baryon, the present work is consistent with the models which study μbaryon [1, 29–32, 36].We would like to point out that although we have used the (two-component) standard distribution in the fits ofpT spectra and T1 (T2) has been used, the values of μq obtained by us are independent of models and parameters. In fact, μq is only related to kj if Tch is known. We can use directly the yield ratios of data to obtain μq. The reason why we use the function form instead of data is to extend pT spectrum in intermediate region to low and high regions where the data are not available. In our opinion, the function that fitted the data in intermediate pT region can predict approximately the trends in low and high pT regions.


## 4. Conclusions

In summary, we found a good fit of the transverse momentum spectra of charged particles produced in central Au+Au collisions at the RHIC at its BES energies. It is shown that the (two-component) standard distribution successfully fitted the data measured at midrapidity by the STAR Collaboration, though other distributions are also acceptable. The effective temperature parameter increases with the increase of the particle mass and the collision energy.At BES energies, the chemical potentials of light flavor quarks were obtained from the yield ratios of negatively to positively charged particles in given transverse momentum ranges available in experiments. At low energy, the chemical potentials of up and down quarks are consistent but differ from that of strange quark. At high energy, the three chemical potentials seem to deviate from each other, and they finally approach zero at very high energy.From the lowest BES energy to the highest RHIC energy, with the increase of logarithmic collision energy, an exponential decrease of the chemical potentials of light flavor quarks is observed. The similarity in up and down quark masses renders the similarity in their chemical potentials. The difference between the chemical potentials of up (or down) and strange quarks is caused by their different masses. The difference between the chemical potentials changes from dozens of MeV to a few MeV. The decrease in chemical potential indicates that the mean free path of produced quarks changes from a small value to a large one.



---

*Source: 6047960-2018-10-11.xml*