# Energy-Dependent Chemical Potentials of Light Hadrons and Quarks Based on Transverse Momentum Spectra and Yield Ratios of Negative to Positive Particles

Authors: Xing-Wei He; Feng-Min Wu; Hua-Rong Wei; Bi-Hai Hong
Journal: Advances in High Energy Physics
 (2020)

---

## Abstract

We describe the transverse momentum spectra or transverse mass spectra ofπ±, K±, p, and p¯ produced in central gold-gold (Au-Au), central lead-lead (Pb-Pb), and inelastic proton-proton (pp) collisions at different collision energies range from the AGS to LHC by using a two-component (in most cases) Erlang distribution in the framework of multisource thermal model. The fitting results are consistent with the experimental data, and the final-state yield ratios of negative to positive particles are obtained based on the normalization constants from the above describing the transverse momentum (or mass) spectra. The energy-dependent chemical potentials of light hadrons (π, K, and p) and quarks (u, d, and s) in central Au-Au, central Pb-Pb, and inelastic (pp) collisions are then extracted from the modified yield ratios in which the contributions of strong decay from high-mass resonance and weak decay from heavy flavor hadrons are removed. The study shows that most types of energy-dependent chemical potentials decrease with increase of collision energy over a range from the AGS to LHC. The curves of all types of energy-dependent chemical potentials, obtained from the fits of yield ratios vs. energy, have the maximum at about 3.510 GeV, which possibly is the critical energy of phase transition from a liquid-like hadron state to a gas-like quark state in the collision system. At the top RHIC and LHC, all types of chemical potentials become small and tend to zero at very high energy, which confirms that the high energy collision system possibly changes completely from the liquid-like hadron-dominant state to the gas-like quark-dominant state and the partonic interactions possibly play a dominant role at the LHC.


---

## Body

## 1. Introduction

The critical energy of phase transition [1–4] is important for studying the quantum chromodynamics (QCD) phase diagram [5, 6] and the properties of quark-gluon plasma (QGP) [7–9], so more and more scientists devote to finding the critical energy. The experiments performed on the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), especially the beam energy scan program at the RHIC, deal with a collision energy range from a few to several tens of GeV [1, 7, 10, 11], which may contain the energy of the critical end point of hadron-quark phase transition [1–4, 12]. The STAR Collaboration found that the critical energy may be or below 19.6 GeV (unless otherwise noted, the energy values presented in this paper are in the center-of-mass coordinate system) [1]. One study based on yield ratio (the yield ratio of negative to positive particles) and the correlation between collision energy and transverse momentum indicated that the critical energy maybe range from 11.5 GeV to 19.6 GeV [1, 13–15], while another study based on yield ratio showed that the critical energy may be about 4 GeV [12]. Studies about a striking pattern of viscous damping and an excitation function for (Rout2−Rside2) extracted for central collisions indicated the critical energy may be close to 62.4 GeV [16–18]. It is not hard to see that the value of critical energy has not been determined so far, so finding the critical energy arouses our great interest.Lattice QCD [19–21], a powerful tool to investigate the QGP matter in high-temperature and high-density system, indicates that the critical end point (CEP) of phase transition on QCD phase diagram is a crossover at small chemical potentials or high collision energies [22, 23]. So it is important to study baryon chemical potential for finding the CEP on QCD phase diagram. When collisions occur at high energy, especially at RHIC and LHC, the collision system probably creates the QGP matter [24–26] where the partonic interactions play an important role, and the baryon chemical potential is small, even close to 1 MeV or zero [12, 27–29]. While when energy is not very high, the transition from hadron to quark has not yet taken place in the collision system, where the hadronic interactions play an important role [1, 13–15], and the value of baryon chemical potential is larger. We could predict that the chemical potential corresponding to the CEP should be an inflection point or abrupt change point in chemical potential-energy plane. It is therefore worthwhile to study the trend of chemical potential with energy.The yield ratio of negative to positive particles is an important quantity in high energy study. Generally, one can get yield ratio by many ways. One way is to directly collect the values of yield ratio from the productive international collaborations, which is a rapid and convenient method. Another way needs the aid of the extracted normalization constant in describing the transverse momentum spectra of negative and positive particles with consistent statistical law, but the workload is huge. In this paper, due to the fact that the experiment data of some particles correspond to a narrow range of transverse momentum (pT), we adopt the second method to obtain a relatively accurate result for the normalization constant being extracted from a wider range of transverse momentum distribution. In addition, one can extract the yield ratio as a fitting parameter in many models. One can obtain some information about the very hot and dense nuclear matter by yield ratio. For example, one can analyze some statistical thermal models to extract temperature, baryon chemical potential (μB) at chemical freeze-out, and so on by describing the ratios of hadron yields, and further to establish the “line of chemical freeze-out” [30] by which one can continue to study QGP, QCD phase transition, and QCD phase diagram [7, 31–35]. Except for baryon chemical potential, one can study the chemical potentials of other hadrons and quarks because these chemical potentials are also important in studying collision system evolution and particle production.In the present work, we describe the transverse momentum (pT) or transverse mass (mT) spectra of π±, K±, p, and p¯ produced in central gold-gold (Au-Au), central lead-lead (Pb-Pb), and inelastic proton-proton (pp) collisions in mid-rapidity interval (in most cases) over a center-of-mass energy (sNN) range from the AGS to LHC [36–51] by using a two-component (in most cases) Erlang distribution [52, 53] in the framework of a multisource thermal model [53–55] and obtain the yield ratios, π−/π+, K−/K+, and p¯/p, of negative to positive particles according to the extracted normalization constants. The energy-dependent chemical potentials of light hadrons (π, K, and p) and quarks (u, d, and s) in central Au-Au, central Pb-Pb, and inelastic (pp) collisions are then extracted from the modified yield ratios in which the contributions of strong decay from high-mass resonance and weak decay from heavy flavor hadrons are removed.


## 2. The Model and Formulism

According to our method, to obtain the normalization constants, we need firstly to describe thepT spectra of π±, K±, p, and p¯ with a multicomponent Erlang distribution [52, 53] which is in the framework of a multisource thermal model [53–55]. The model assumes that many emission sources are formed in high energy collisions and are classified into a few groups due to the existent of different interacting mechanisms in the collisions and different event samples in experiment measurements. The sources in the same group have the same excitation degree and stay at a common local equilibrium state, which can be described by an Erlang pT distribution. All emission sources in different groups result in the final-state distribution, which can be described by a multicomponent Erlang pT distribution.The multicomponent Erlang distribution based on the above multisource thermal model has the following form. According to thermodynamic system, particles generated from one emission source obey to an exponential distribution of transverse momentum,(1)fijptij=1ptijexp−ptijptij,where ptij is the transverse momentum of the ith source in the jth group, and ptij is the mean value of ptij. We assume that the source number in the jth group and the transverse momentum of the mj sources are denoted by mj and pT, respectively. All the sources in the jth group then result in the folding result of exponential distribution
(2)fjpT=pTmj−1mj−1!ptijmjexp−pTptij,which is the normalized Erlang distribution. The contribution of the l group of sources can be expressed as
(3)fpT=∑j=1lkjfjpT,where kj denotes the relative weight contributed by the jth group and meets the normalization ∑j=1lkj=1. This is the multicomponent Erlang distribution.In fact, in the present work, we describe the transverse momentum spectra of final-state light flavor particles by using a two-component Erlang distribution, where one component reflects the soft excitation process, while the other one reflects the hard scattering process. The soft process corresponding to low-pT region is regarded as the contribution of the interactions among a few sea quarks and gluons, and the hard process corresponding to high-pT region is regarded as originating from a harder head-on scattering between a few valent quarks. Due to the fact that the experimental data of some particles correspond to a narrow range of pT, we adopt one-component Erlang distribution to fit these data.Some experimental data we collect are about transverse mass distribution, notpT distribution, so we give the transformational relation between pT distribution and mT distribution based on the relation between pT and mT (mT=pT2+m02, where m0 is the rest mass of particle), i.e.,
(4)dNNdmT=mTpTdNNdpT.The same as in our previous work [28], in the present work, we only calculate the chemical potentials of some light hadrons (π, K, and p) and some light quarks (u, d, and s). For the hadrons containing c or b quark, considering that there is a lack of the experimental data of pT spectra continuously varying with energy, we do not calculate the chemical potentials of the hadrons containing c or b quark and c and b quarks. In addition, due to the lifetimes of the hadrons containing t quark being too short to measure, we also can not obtain the chemical potentials of the hadrons containing t quark. According to the statistical arguments based on the chemical and thermal equilibrium within the thermal and statistical model [56], we can get the relations between antiparticle to particle (negative to positive particle) yield ratios and chemical potentials of hadrons to be [42, 56, 57]
(5)kπ=exp−2μπTch,kK=exp−2μKTch,kp=exp−2μpTch,where kπ, kK, and kp denote the yield ratios of antiparticles, π−, K−, and p¯ to particles, π+, K+, and p, respectively, and μπ, μK, and μp represent the chemical potentials of π, K, and p, respectively. In addition, Tch represents the chemical freeze-out temperature of interacting system and can be empirically obtained by the following formula
(6)Tch=Tlim11+exp2.60−lnsNN/0.45within the framework of a statistical thermal model of noninteracting gas particles with the assumption of standard Maxwell-Boltzmann statistics [7, 8, 58], where the “limiting” temperature Tlim is 0.164 GeV, and sNN is in the unit of GeV [58, 59].Assuming thatμu, μd, and μs represent the chemical potentials of u, d, and s quarks, respectively, and according to Equation (5) and references [12, 57, 60] under the same value of chemical freeze-out temperature, the yield ratios in terms of quark chemical potentials can be written as
(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−2μu−μdTch.Based on Equations (5) and (7), one can obtain the chemical potentials of hadrons and quarks in terms of yield ratios, respectively,(8)μπ=−12Tch·lnkπ,μK=−12Tch·lnkK,μp=−12Tch·lnkp,(9)μu=−16Tch·lnkπ·kp,μd=−16Tch·lnkπ−2·kp,μs=−16Tch·lnkπ·kK−3·kp.In the present work, by describing thepT (or mT) spectra of some light particles, π±, K±, p, and p¯ in central Au-Au, central Pb-Pb, and inelastic pp collisions in mid-rapidity interval at collision energy from the AGS to LHC with a two-component (in most cases) Erlang distribution, we obtain the yield ratios, π−/π+, K−/K+, and p¯/p based on the extracted normalization constants, and the chemical potentials of light hadrons (π, K, and p) and light quarks (u, d, and s) based on the yield ratios modified by removing the contributions of strong decay from high-mass resonance and weak decay from heavy flavor hadrons. Then the dependencies of chemical potentials on sNN are analyzed.


## 3. Results and Discussion

Figure1 shows the transverse mass distributions of π± (Figures 1(a) and 1(b)) and K± (Figures 1(c) and 1(d)) produced in central (0–5%) Au-Au collisions at mid-rapidity in the center-of-mass energy (sNN) range from 2.67 to 4.84 GeV, where dN/dy on axis denotes the rapidity density. The experimental data represented by different kinds of symbols were measured by the E895 Collaboration [36] for π± at 2.67, 3.31, 3.81, and 4.28 GeV and the E866 and E917 Collaborations [37, 38] for K± at 3.31, 3.81, 4.28, and 4.84 GeV. The data at each energy are scaled by suitable factors for clarity. The plotted error bars include both statistical and systematic uncertainties for π± and only statistical uncertainty for K±. The solid curves are our results calculated by using the two-component Erlang distribution. The values of free parameters (m1, pti1, k1, m2, and pti2), normalization constant (N0), and χ2 per degree of freedom (χ2/dof), as well as p values corresponding to the two-component Erlang distribution are listed in Table 1, where the normalization constant is for comparison between curve and data. Here, χ2 is calculated according to the following formula of
(10)χ2=∑Niexp−Nitheo2σ2where Niexp, Nitheo, and σ denote experimental value, theoretical value, and error value, respectively. One can see that the two-component Erlang distribution can well describe the experimental data of the considered particles in Au-Au collisions at the AGS. The values of m2 corresponding to high-pT region for different particles are 2, which reflects that the hard process origins from a hard head-on scattering between two valent quarks, while the values of m1 corresponding to low-pT region for different particles are 3, which reflects that the soft process origins from the interaction among a few sea quarks and gluons. The values of weight factor k1 of soft excitation process are more than 50%, which shows that soft excitation is the main excitation process, and the normalization constant N0 increases with increase of energy. It should be noted that the particle yield ratio is represented by N0 from the spectrum of negative or positive particles. The relative value of N0 is enough to obtain the particle yield ratio.Figure 1
Transverse mass spectra for positive (π+ and K+) and negative (π− and K−) particles produced in central Au-Au collisions at mid-rapidity over an energy range from 2.67 to 4.84 GeV. The experimental data represented by the symbols are measured by the E895 Collaboration [36] for π± at 2.67, 3.31, 3.81, and 4.28 GeV and the E866 and E917 Collaborations [37, 38] for K± at 3.31, 3.81, 4.28, and 4.84 GeV. The data at each energy are scaled by successive powers of 2 for clarity. The plotted error bars include both statistical and systematic uncertainties for π± and only statistical uncertainty for K±. The solid curves are our results calculated by using the two-component Erlang distribution.Table 1
Values of free parameters, normalization constant, andχ2/dof, as well as p values corresponding to two-component Erlang pT (or mT) distribution for Au-Au collisions in Figures 1 and 2.
FiguresNN (GeV)Particlem1pti1 (GeV/c)k1m2pti2 (GeV/c)N0χ2/dofp valueFigure1(a)2.67π+30.078±0.0010.74±0.0120.182±0.01212.051±0.3657.043/170.983Figure1(b)π−30.061±0.0010.74±0.0220.163±0.01020.811±0.36511.690/240.983Figure1(a)3.31π+30.083±0.0020.68±0.0220.251±0.01027.819±0.7473.768/221Figure1(c)K+30.149±0.0030.91±0.0620.162±0.0202.344±0.1722.144/60.906Figure1(b)π−30.067±0.0010.58±0.0120.176±0.00438.257±0.77719.052/340.982Figure1(d)K−30.149±0.0080.91±0.0620.162±0.0600.186±0.0178.857/70.263Figure1(a)3.81π+30.084±0.0020.60±0.0220.280±0.01039.084±0.7702.654/171Figure1(c)K+30.190±0.0110.80±0.0620.202±0.0234.829±0.2850.431/50.994Figure1(b)π−30.068±0.0010.51±0.0120.187±0.00450.197±0.79736.777/340.341Figure1(d)K−30.174±0.0070.88±0.0620.185±0.0260.572±0.0073.598/50.608Figure1(a)4.28π+30.081±0.0030.51±0.0220.267±0.01649.619±1.1035.465/140.978Figure1(c)K+30.197±0.0110.94±0.0420.235±0.0307.717±0.0311.304/30.728Figure1(b)π−30.072±0.0030.51±0.0220.192±0.00262.627±0.74440.248/340.213Figure1(d)K−30.150±0.0090.80±0.0720.235±0.0321.355±0.1071.738/40.784Figure1(c)4.84K+30.187±0.0060.87±0.0720.267±0.04010.798±0.3835.328/50.377Figure1(d)K−30.180±0.0030.91±0.0620.247±0.0302.025±0.0695.463/50.362Figure2(a)7.7π+20.172±0.0040.63±0.0620.233±0.00496.122±3.32613.970/200.832K+30.197±0.0030.93±0.0720.300±0.03320.070±0.6626.218/170.991p40.215±0.0030.89±0.0820.270±0.05452.211±2.2035.070/231Figure2(d)π−20.149±0.0060.52±0.0320.219±0.003107.122±3.5789.297/200.979K−30.186±0.0030.92±0.0820.285±0.0577.208±0.26814.213/170.652p¯40.232±0.0130.71±0.1420.334±0.0660.412±0.0202.045/90.991Figure2(b)11.5π+20.153±0.0070.52±0.0420.236±0.003125.208±4.5832.186/201K+30.201±0.0030.87±0.1320.262±0.04824.436±0.7181.486/191p40.211±0.0030.90±0.0820.244±0.04842.924±1.8717.386/220.998Figure2(e)π−20.146±0.0080.51±0.0420.230±0.003135.170±5.7851.593/201K−30.191±0.0030.92±0.0820.221±0.04412.017±0.3821.338/171p¯40.209±0.0040.92±0.0820.234±0.0461.374±0.05917.977/170.390Figure2(c)19.6π+20.156±0.0080.57±0.0520.249±0.005165.077±7.9570.737/201K+30.198±0.0050.67±0.1320.294±0.00929.706±0.8791.316/201p40.224±0.0040.79±0.0720.278±0.04034.690±1.3744.885/231Figure2(f)π−20.145±0.0080.56±0.0420.246±0.004176.077±8.3810.682/201K−30.192±0.0040.64±0.0720.292±0.00718.620±0.5852.950/201p¯40.222±0.0040.93±0.0720.247±0.0493.937±0.1697.144/160.970Figure2 presents the transverse momentum spectra of π±, K±, p, and p¯ in central (0–5%) Au-Au collisions at center-of-mass energy 7.7 (Figures 2(a) and 2(d)), 11.5 (Figures 2(b) and 2(e)), and 19.6 (Figures 2(c) and 2(f)) GeV. The symbols represent the experimental data recorded by the STAR Collaboration in the mid-rapidity range ∣y∣<0.1 [39]. The uncertainties are statistical and systematic added in quadrature. The curves are our results fitted by using the two-component Erlang distribution. The values of m1, pti1, k1, m2, pti2, N0, and χ2/dof, as well as p values corresponding to the two-component Erlang distribution are given in Table 1. It is not hard to see that the experimental data can be well fitted by the two-component Erlang distribution. Similarly, the values of m2 are 2, and the values of m1 are 2, 3, and 4. The values of weight factor k1 are more than 50%, and N0 in most cases increases with increase of collision energy.Figure 2
Mid-rapidity transverse momentum spectra for positive (π+, K+, and p) and negative (π−, K−, and p¯) particles produced in central Au-Au collisions at (a, d) 7.7, (b, e) 11.5, and (c, f) 19.6 GeV. The symbols represent the experimental data recorded by the STAR Collaboration [39]. The errors are the combined statistical and systematic ones, and the curves are our results by the two-component Erlang distribution.Figure3 gives the same as Figure 2 but for Au-Au collisions at 27 (Figures 3(a) and 3(d)), 39 (Figures 3(b) and 3(e)), and 62.4 (Figures 3(c) and 3(f)) GeV. All the experimental data were recorded by the STAR Collaboration [39, 40]. The results calculated by using the two-component Erlang distribution are shown in the solid curves, where the values of corresponding free parameters, normalization constant, and χ2/dof, as well as p values are shown in Table 2. Obviously, the calculation results by the two-component Erlang distribution are in good agreement with the experimental data of the considered particles in Au-Au collisions. Once more, the values of m2 are 2, and the values of m1 are 2, 3, and 4. The values of weight factor k1 are more than 50%, and N0 in most cases increases with increase of collision energy.Figure 3
The same as Figure2 but for Au-Au collisions at (a, d) 27, (b, e) 39, and (c, f) 62.4 GeV.Table 2
Values of free parameters, normalization constant, andχ2/dof, as well as p values corresponding to two-component Erlang pT (or mT) distribution for Au-Au collisions in Figures 3 and 4.
FiguresNN (GeV)Particlem1pti1 (GeV/c)k1m2pti2 (GeV/c)N0χ2/dofp valueFigure3(a)27π+20.152±0.0070.51±0.0520.249±0.003182.402±6.2022.449/201K+30.205±0.0040.97±0.0220.575±0.11529.993±1.0462.393/201p40.226±0.0050.85±0.1220.284±0.05630.191±1.1654.566/170.999Figure3(d)π−20.164±0.0060.66±0.0420.264±0.004186.402±6.7102.325/201K−30.198±0.0020.97±0.0120.531±0.05221.872±0.8198.210/190.984p¯40.228±0.0040.92±0.0820.343±0.0685.877±0.2395.235/160.994Figure3(b)39π+20.155±0.0080.54±0.0520.265±0.004185.159±7.2583.265/201K+30.211±0.0030.94±0.0620.359±0.04531.219±1.0245.650/200.999p40.239±0.0050.79±0.0920.293±0.05526.115±1.0813.078/161Figure3(e)π−20.153±0.0100.53±0.0420.258±0.004191.409±7.1591.088/201K−30.206±0.0030.86±0.0620.352±0.02024.658±0.7155.446/200.999p¯40.233±0.0040.92±0.0820.280±0.0568.086±0.3195.566/170.996Figure3(c)62.4π+20.172±0.0030.65±0.0520.274±0.012232.461±1.7200.261/40.992K+30.246±0.0040.75±0.0420.256±0.01539.598±2.0320.463/40.977p40.252±0.0010.92±0.0320.436±0.02728.457±0.1542.891/90.968Figure3(f)π−20.175±0.0020.67±0.0520.269±0.008234.954±1.2690.825/40.935K−30.223±0.0150.90±0.1020.239±0.04732.071±1.8446.873/40.143p¯40.281±0.0020.77±0.0320.435±0.01015.011±0.11711.697/100.306Figure4(a)130π+20.128±0.0030.58±0.0220.262±0.020288.154±9.35019.320/80.013K+30.166±0.0070.65±0.0320.382±0.03046.172±1.88311.031/70.137p40.224±0.0080.61±0.0420.420±0.01529.451±0.55212.703/110.313Figure4(c)π−20.115±0.0100.60±0.0220.269±0.012291.470±9.10025.039/80.002K−30.143±0.0070.59±0.0420.326±0.01644.407±1.93311.880/70.104p¯40.240±0.0070.79±0.0420.466±0.03018.813±0.64417.215/110.102Figure4(b)200π+20.162±0.0080.62±0.0220.291±0.004314.469±3.50037.817/190.006K+30.208±0.0030.60±0.0320.409±0.00647.468±0.79328.553/100.002p40.266±0.0030.92±0.0220.568±0.00715.345±0.26619.924/160.224Figure4(d)π−20.179±0.0040.65±0.0120.297±0.001293.451±2.50038.653/220.016K−30.223±0.0020.79±0.0120.480±0.00742.471±0.4569.724/100.465p¯40.270±0.0030.98±0.0220.650±0.00611.195±0.22526.068/160.053ThepT spectra of π±, K±, p, and p¯ in central (0–5%) Au-Au collisions at 130 (Figures 4(a) and 4(c)) and 200 (Figures 4(b) and 4(d)) GeV are displayed in Figure 4. The symbols also denote the experimental data recorded by the PHENIX Collaboration [41, 42]. The data for each type of particle are divided by suitable factors for clarity. The error bars indicate the combined uncorrelated statistical and systematic uncertainties for 130 GeV and are statistical only for 200 GeV. The curves are the two-component Erlang model fits to the spectra. The values of all free parameters, normalization constant, and χ2/dof, as well as p values corresponding to the two-component Erlang distribution are listed in Table 2. Similarly, our calculation results with the two-component Erlang model are consistent with the experimental data. The values of m2 are 2, and the values of m1 are 2, 3, and 4. The values of weight factor k1 are more than 50%, and N0 in most cases increases with increase of collision energy.Figure 4
The same as Figure2 but for Au-Au collisions at (a, c) 130 and (b, d) 200 GeV.Figure5 exhibits the mT spectra of π± at 0<y<0.2, K± at ∣y∣<0.1, p, and p¯ produced in central Pb-Pb collisions at 6.3 (Figures 5(a) and 5(d)), 7.7 (Figures 5(b) and 5(e)), and 8.8 (Figures 5(c) and 5(f)) GeV. The experimental data, represented by symbols, were taken by the NA49 Collaboration [43–45], where p and p¯ were done near mid-rapidity and covered the rapidity intervals of 1.5<y<2.2 (yc.m.=1.88) for 6.3 GeV, 1.6<y<2.3 (yc.m.=2.08) for 7.7 GeV, and 1.9<y<2.3 (yc.m.=2.22) for 8.8 GeV. The error bars on the spectra points are statistical only. The curves are fits of two-component Erlang function to the spectra. The values of free parameters, normalization constant, and χ2/dof, as well as p values are summarized in Table 3. We can see that the experimental data for all hadrons and energies are well described by the fit function. The values of m2 are 2, and the values of m1 are 2, 3, and 4. The values of weight factor k1 are more than 50%, and N0 increases with increase of collision energy.Figure 5
Transverse mass spectra forπ±, K±, p, and p¯ in central Pb-Pb collisions at (a, d) 6.3, (b, e) 7.7, and (c, f) 8.8 GeV. The symbols represent the experimental data taken by the NA49 Collaboration [43–45]. The errors are statistical only. The curves are fits of two-component Erlang function to the spectra.Table 3
Values of free parameters, normalization constant, andχ2/dof, as well as p values corresponding to two-component Erlang pT (or mT) distribution for Pb-Pb collisions in Figures 5 and 6.
FiguresNN (GeV)Particlem1<pti1> (GeV/c)k1m2<pti2> (GeV/c)N0χ2/dofp valueFigure5(a)6.3π+30.095±0.0030.51±0.0320.228±0.00572.088±2.3799.764/100.461K+30.194±0.0040.90±0.1020.295±0.05916.508±0.6446.875/40.143p40.209±0.0021--2.889±0.07827.988/110.003Figure5(d)π−30.082±0.0030.51±0.0320.225±0.00683.773±2.7654.884/100.899K−30.179±0.0050.86±0.0920.328±0.0815.644±0.18613.617/40.009p¯40.234±0.0181--0.004±0.0011.420/40.841Figure5(b)7.7π+20.172±0.0030.90±0.0220.436±0.07987.611±2.54110.081/100.433K+30.202±0.0041--20.584±0.6386.746/40.456p40.214±0.0040.87±0.1220.420±0.0555.278±0.3062.421/80.965Figure5(e)π−20.151±0.0040.79±0.0320.317±0.02298.711±3.2574.493/100.922K−30.205±0.0060.82±0.0820.220±0.0388.170±0.30218.012/140.206p¯40.233±0.0280.75±0.1520.580±0.1160.021±0.0040.4253/10.514Figure5(c)8.8K+30.203±0.0040.91±0.0520.205±0.04121.556±0.6255.397/40.249p40.221±0.0030.92±0.0820.270±0.0548.605±0.4227.546/80.479Figure5(f)K−40.149±0.0040.75±0.1520.280±0.0318.062±0.28213.337/40.010p¯40.219±0.0100.92±0.0720.370±0.0740.084±0.0103.408/40.492Figure6(a)12.3K+30.204±0.0040.90±0.0720.220±0.02724.872±0.6966.261/80.618p40.224±0.0030.85±0.0720.290±0.04613.871±0.5698.960/80.346Figure6(d)K−30.192±0.0040.83±0.0720.327±0.06511.307±0.3395.198/120.951p¯30.312±0.0100.94±0.0620.370±0.0740.348±0.0223.947/40.413Figure6(b)17.3π+20.199±0.0040.91±0.0420.330±0.03817.753±0.78144.048/200.001K+40.162±0.0050.51±0.0620.289±0.0144.246±0.23819.806/240.708p30.292±0.0080.96±0.0220.322±0.0200.403±0.0146.545/120.886Figure6(e)π−20.188±0.0030.88±0.0320.307±0.01418.926±0.66224.600/200.217K−40.150±0.0050.51±0.1020.286±0.0112.378±0.1319.820/240.995p¯30.304±0.0110.93±0.0320.170±0.0340.028±0.0027.744/120.805Figure6(c)2760π+20.177±0.0070.58±0.0320.364±0.004118.422±3.79014.385/350.999K+30.283±0.0090.64±0.1020.473±0.02117.258±0.5522.657/301p40.360±0.0050.82±0.0520.470±0.0645.311±0.1759.093/361Figure6(f)π−20.180±0.0070.58±0.0320.364±0.004118.191±3.78212.640/350.999K−30.277±0.0100.59±0.0920.467±0.01617.198±0.6022.851/301p¯40.360±0.0060.81±0.0720.470±0.0765.253±0.17310.028/361Figure6 presents the mT and pT spectra of π±, K±, p, and p¯ in central Pb-Pb collisions at 12.3 (Figures 6(a) and 6(d)), 17.3 (Figures 6(b) and 6(e)), and 2760 (Figures 6(c) and 6(f)) GeV, where σtrig on the vertical axis denotes the interaction cross section satisfying a T0 centrality trigger. The symbols represent the experimental data reported by the NA49 Collaboration for 12.3 GeV at mid-rapidity [∣y∣<0.1 for K± [45], 2.2<y<2.6 (yc.m.=2.57) for p and p¯ [44]], the NA44 Collaboration for 17.3 GeV near mid-rapidity (2.4<y<3.1 for π±, 2.4<y<3.5 for K±, and 2.3<y<2.9 for p and p¯) [46], and the ALICE Collaboration for 2760 GeV at mid-rapidity ∣y∣<0.5 [47]. Some data for different particles are divided by suitable factors for clarity. The errors are statistical for 12.3 GeV, are systematic for 17.3 GeV, and are quadratic sum of statistical errors and systematic errors for 2760 GeV. The curves represent the two-component Erlang fits. The values of free parameters, normalization constant, and χ2/dof, as well as p values are summarized in Table 3. Obviously, the experimental data for all particles at all energies are in good agreement with the fits. The values of m2 are 2, and the values of m1 are 2, 3, and 4. The values of weight factor k1 are more than 50%, and N0 increases with increase of collision energy.Figure 6
Transverse mass and momentum spectra forπ±, K±, p, and p¯ at mid-rapidity in central Pb-Pb collisions at (a, d) 12.3, (b, e) 17.3, and (c, f) 2760 GeV. The symbols represent the experimental data reported by the NA49 Collaboration for 12.3 GeV at mid-rapidity [∣y∣<0.1 for K± [45], 2.2<y<2.6 (yc.m.=2.57) for p and p¯ [44]], the NA44 Collaboration for 17.3 GeV near mid-rapidity (2.4<y<3.1 for π±, 2.4<y<3.5 for K±, and 2.3<y<2.9 for p and p¯) [46], and the ALICE Collaboration for 2760 GeV at mid-rapidity ∣y∣<0.5 [47]. Some data for different particles are divided by suitable factors for clarity. The errors are statistical for 12.3 GeV, are systematic for 17.3 GeV, and are quadratic sum of statistical errors and systematic errors for 2760 GeV. The curves are fits of two-component Erlang function to the spectra.Figure7 shows the pT spectra of π± (Figures 7(a) and 7(b)) and K± (Figures 7(c) and 7(d)) produced in mid-rapidity y≈0 inelastic pp collisions at 6.3, 7.7, 8.8, 12.3, and 17.3 GeV. The measurements were performed at the CERN-Super Proton Synchrotron (SPS) by the large acceptance NA61/SHINE hadron spectrometer [48]. Spectra at different energies are scaled by appropriate factors for better visibility. The error bars on data points correspond to combined statistical and systematic uncertainties. The curves are our fitting results by using the one- or two-component Erlang function. For some curves, we use one-component Erlang function because the number of corresponding experimental data points is small. Due to the proportion of the second component is small, it has little effect on the calculated particle ratio, despite the absence of the second component. The values of free parameters, normalization constant, and χ2/dof, as well as p values are given in Table 4. As can be seen, the fits for all hadrons at all energies are in good agreement with the experimental data. The values of m2 and m1 are 2 and 3, respectively. The values of weight factor k1 are more than 50%. It should be noted that the dof for π− at 12.3 GeV in Table 4 is zero, which means the dash curve in Figure 7(b) is drawn to guide the eye.Figure 7
Transverse momentum spectra for (a, b)π± and (c, d) K± in y≈0 inelastic pp collisions at SPS energies (6.3, 7.7, 8.8, 12.3, and 17.3 GeV). The symbols represent the experimental data reported by the NA61/SHINE Collaboration [48]. Spectra at different energies are scaled by appropriate factors for better visibility. The errors are quadratic sum of statistical errors and systematic errors. The curves are fits of one- or two-component Erlang function to the spectra, where the dash curve in the figure is drawn to guide the eye.Table 4
Values of free parameters, normalization constant, andχ2/dof, as well as p values corresponding to one- or two-component Erlang pT distribution for inelastic pp collisions in Figures 7 and 8.
FiguresNN (GeV)Particlem1<pti1> (GeV/c)k1m2<pti2> (GeV/c)N0χ2/dofp valueFigure7(a)6.3π+30.103±0.0060.77±0.0920.163±0.0110.538±0.0512.193/60.901Figure7(c)K+30.143±0.0121--0.056±0.00112.005/20.002Figure7(b)π−30.106±0.0020.90±0.0620.116±0.0100.358±0.0381.874/40.759Figure7(d)K−30.135±0.0080.80±0.0720.184±0.0200.016±0.0032.566/40.632Figure7(a)7.7π+30.102±0.0080.69±0.0820.202±0.0110.698±0.0860.246/50.998Figure7(c)K+30.145±0.0131--0.066±0.0060.862/20.650Figure7(b)π−30.110±0.0420.73±0.1220.170±0.0120.482±0.0780.084/30.993Figure7(d)K−30.158±0.0171--0.028±0.0010.487/20.784Figure7(a)8.8π+30.107±0.0040.80±0.0420.232±0.0150.711±0.0550.370/40.985Figure7(c)K+30.144±0.0040.75±0.0720.190±0.0600.067±0.0015.585/20.061Figure7(b)π−30.113±0.0050.75±0.0620.203±0.0180.461±0.0380.653/20.722Figure7(d)K−30.162±0.0070.75±0.0820.193±0.0900.030±0.0011.179/20.411Figure7(a)12.3π+30.098±0.0090.51±0.0420.198±0.0080.796±0.0680.212/40.995Figure7(c)K+30.155±0.0050.76±0.0720.180±0.0300.076±0.0034.164/20.125Figure7(b)π−30.117±0.0090.80±0.0920.202±0.0160.631±0.0890.165/00Figure7(d)K−30.159±0.0040.64±0.0620.246±0.0300.044±0.0014.512/20.105Figure7(a)17.3π+30.108±0.0080.61±0.0820.198±0.0170.813±0.0490.054/10.817Figure7(c)K+30.154±0.0070.80±0.0620.178±0.0300.077±0.0062.300/30.512Figure7(b)π−30.112±0.0060.77±0.0620.227±0.0200.688±0.0040.032/10.859Figure7(d)K−30.164±0.0080.80±0.0720.178±0.0300.051±0.0033.057/30.382Figure8(a)62.4π+20.160±0.0030.88±0.0120.310±0.0050.965±0.0395.886/261K+30.175±0.0030.84±0.0220.455±0.0190.085±0.0033.026/161p40.186±0.0040.83±0.0420.404±0.0160.040±0.0026.489/271Figure8(c)π−20.160±0.0030.89±0.0120.310±0.0040.966±0.0358.208/261K−30.168±0.0040.76±0.0320.362±0.0140.076±0.0031.856/161p¯40.186±0.0030.95±0.0220.552±0.0460.022±0.0017.113/271Figure8(b)200π+20.170±0.0030.91±0.0120.390±0.0041.071±0.04111.062/270.997K+30.162±0.0040.73±0.0220.434±0.0120.120±0.0051.123/161p40.194±0.0040.88±0.0120.532±0.0080.043±0.00217.121/340.993Figure8(d)π−20.170±0.0030.89±0.0120.370±0.0041.028±0.0408.247/271K−30.156±0.0050.64±0.0320.380±0.0090.119±0.0053.264/161p¯40.192±0.0040.86±0.0220.497±0.0060.036±0.00210.592/341Figure8 presents the pT spectra of π±, K±, p, and p¯ produced in inelastic pp collisions at 62.4 (Figures 8(a) and 8(c)) and 200 (Figures 8(b) and 8(d)) GeV, where E and σ on the vertical axis denote the particle energy and cross section, respectively. The data measured by the PHENIX Collaboration in the mid-pseudorapidity range ∣η∣<0.35 [49] are represented in different panels by different symbols. Spectra for different particles are scaled by appropriate factors for better visibility. The error bars are statistical only. The curves are our results fitted by using the two-component Erlang distribution. The values of free parameters, normalization constant, and χ2/dof, as well as p values are given in Table 4. One can see that all the fitting results by using the two-component Erlang function are consistent with the experimental data. The values of m2 are 2, and the values of m1 are 2, 3, and 4. The values of weight factor k1 are more than 50%, and N0 increases with increase of collision energy.Figure 8
Transverse momentum spectra for positive (π+, K+, and p) and negative (π−, K−, and p¯) particles produced in inelastic pp collisions at (a, c) 62.4 and (b, d) 200 GeV. The experimental data represented by the symbols are measured by the PHENIX Collaboration in ∣η∣<0.35 [49]. The errors are statistical only. The curves are fits by the two-component Erlang distribution.Figure9 exhibits the pT spectra of π±, K±, p, and p¯ produced in inelastic pp collisions at 900 GeV (Figures 9(a) and 9(c)) and 2.76 TeV (Figures 9(b) and 9(d)). The symbols also denote the experimental data recorded by the CMS Collaboration in the range ∣y∣<1 [50]. The data for different particles are scaled by suitable factors for clarity. The error bars indicate the combined uncorrelated statistical and systematic uncertainties, and the fully correlated normalization uncertainty is 3.0%. The curves are fits of the two-component Erlang function to the spectra. The values of free parameters, normalization constant, and χ2/dof, as well as p values are summarized in Table 5. We can see that the experimental data are well described by the fit function. The values of m2 are 2, and the values of m1 are 2, 3, and 4. The values of weight factor k1 are more than 50%, and N0 increases with increase of collision energy.Figure 9
Transverse momentum spectra for positive (π+, K+, and p) and negative (π−, K−, and p¯) particles produced in inelastic pp collisions at (a, c) 900 GeV and (b, d) 2.76 TeV. The symbols represent the experimental data recorded by the CMS Collaboration in ∣y∣<1 [50]. The errors are the combined uncorrelated statistical and systematic ones, and the fully correlated normalization uncertainty is 3.0%. The curves are our results fitted by the two-component Erlang distribution.Table 5
Values of free parameters, normalization constant, andχ2/dof, as well as p values corresponding to one- or two-component Erlang pT distribution for inelastic pp collisions in Figures 9 and 10.
FiguresNN (GeV)Particlem1<pti1> (GeV/c)k1m2<pti2> (GeV/c)N0χ2/dofp valueFigure9(a)900π+20.160±0.0020.80±0.0120.390±0.0173.924±0.04314.393/220.887K+30.165±0.0040.52±0.0320.465±0.0180.480±0.0061.769/171p40.171±0.0030.51±0.0420.530±0.0200.212±0.00411.337/270.996Figure9(c)π−20.160±0.0020.80±0.0120.420±0.0183.884±0.04717.291/220.747K−30.165±0.0040.51±0.0320.460±0.0150.472±0.0062.163/171p¯40.168±0.0030.51±0.0220.536±0.0180.202±0.00319.958/270.833Figure9(b)2760π+20.156±0.0020.73±0.0120.384±0.0114.944±0.05911.582/220.965K+30.190±0.0040.51±0.0420.370±0.0120.578±0.0086.135/170.992p40.179±0.0040.51±0.0420.670±0.0270.276±0.00430.949/270.273Figure9(d)π−20.159±0.0020.74±0.0120.387±0.0114.824±0.05816.848/220.772K−30.201±0.0060.51±0.0720.330±0.0150.560±0.0089.175/170.935p¯40.181±0.0030.51±0.0320.610±0.0240.262±0.00432.664/270.208Figure10(a)7000π+20.144±0.0020.60±0.0220.321±0.0076.084±0.06718.074/220.702K+30.236±0.0050.62±0.0520.260±0.0120.700±0.0093.234/171p30.261±0.0060.51±0.0520.600±0.0250.346±0.00512.963/270.989Figure10(c)π−20.150±0.0020.63±0.0220.336±0.0075.984±0.07219.714/220.600K−30.240±0.0050.63±0.0620.251±0.0140.696±0.0086.961/170.984p¯30.260±0.0050.51±0.0520.600±0.0300.342±0.00519.058/270.868Figure10(b)13000π+20.159±0.0040.70±0.0220.441±0.0275.724±0.0928.045/220.997K+30.219±0.0080.71±0.1220.266±0.0310.600±0.0133.199/171p30.264±0.0090.51±0.0620.569±0.0370.314±0.00611.583/260.993Figure10(d)π−20.168±0.0020.73±0.0320.456±0.0395.604±0.07828.200/220.169K−30.230±0.0100.69±0.1020.250±0.0260.596±0.0152.820/171p¯30.270±0.0080.51±0.0620.569±0.0280.306±0.00610.087/260.998Figure10 shows the same as Figure 9, but for s=7 TeV (Figures 10(a) and 10(c)) and s=13 TeV (Figures 10(b) and 10(d)). The symbols also denote the experimental data recorded by the CMS collaboration in the range ∣y∣<1 [50, 51]. The error bars indicate the combined uncorrelated statistical and systematic uncertainties, and the fully correlated normalization uncertainty is 3.0%. The curves are our results fitted by using the two-component Erlang distribution. The values of free parameters, normalization constant, and χ2/dof, as well as p values are summarized in Table 5. It is not hard to see that the experimental data can be well fitted by the two-component Erlang distribution. Similarly, the values of m2 are 2, and the values of m1 are 2 and 3. The values of weight factor k1 are more than 50%, and N0 increases with increase of collision energy.Figure 10
The samepT spectra as Figure 9, but for (a, c) s=7 TeV and (b, d) s=13 TeV. The symbols also represent the experimental data recorded by the CMS Collaboration in ∣y∣<1 [50, 51]. Similarly, the errors are the combined uncorrelated statistical and systematic ones, and the fully correlated normalization uncertainty is 3.0%. The curves are our results fitted by the two-component Erlang distribution.According to the extracted normalization constants from the above comparisons, the three types of energy-dependent yield ratios,π−/π+, K−/K+, and p¯/p, of negative to positive particles from different collision systems are obtained and are shown in Figure 11. The black, red, and blue circles denote, respectively, the calculated results from inelastic pp, central Au-Au, and central Pb-Pb collisions. For comparison, the black, red, and blue triangles correspondingly denote the experimental results [36–48, 50, 51, 61–68] from inelastic or NSD pp, central Au-Au, and central Pb-Pb collisions, respectively. One can see that our calculation results based on transverse momentum (or mass) spectra are consistent with the experimental data. That means the differences between our calculation results and experimental data are very small, which indicates that our calculation method is correct.Figure 11
sNN-dependent yield ratios: (a) π−/π+, (b) K−/K+, and (c) p¯/p. The black, red, and blue circles denote, respectively, the calculated results from inelastic pp, central Au-Au, and central Pb-Pb collisions according to normalization constants, and the triangles correspondingly denote the experimental results for comparison.In fact, the yield ratios we calculate on the basis of transverse momentum (or mass) spectra, are at the stage of kinetic freeze-out and are affected by strong decay from high-mass resonance and weak decay from heavy flavor hadrons. In order to obtain the yield ratios at the stage of chemical freeze-out, one needs to remove the contributions of strong decay and weak decay from the above yield ratios. According to reference [69], we remove the contributions of strong and weak decays and obtain the modified (primary) yield ratios, kπ, kK, and kp. The calculation result shows that strong decay affects mainly kπ and kK, and weak decay affects mainly kp. Strong decay can pull down kπ and lift kK, and weak decay can lift kp, although these two decays do not significantly affect the primary yield ratios as a whole.The three types of yield ratios show regular trends with increase of collision energy.kπ from pp collisions, kK, and kp increase with increase of collision energy, and kπ from central Au-Au and central Pb-Pb collisions decreases with increase of collision energy. To see more clearly the dependences of the three types of yield ratios on collision energy, we show the logarithms of the three yield ratios, lnkπ, lnkK, and lnkp, with 1/sNN in Figure 12. The black squares, red circles, and blue triangles denote the calculated results from inelastic pp, central Au-Au, and central Pb-Pb collisions at mid-rapidity, respectively. One can see that lnkK and lnkp show obviously linear dependence on 1/sNN, which should be fitted by linear functions for clarity. lnkK and lnkp from all collision systems mentioned above, decrease monotonously with increase of 1/sNN and can be described by the below linear functions of
(11)lnkK=−8.690±0.187sNN+−0.021±0.030,lnkp=−45.034±0.637sNN+−0.035±0.030,with χ2/dof to be 7.589/26 and 3.675/17, respectively. lnkπ displays different behavior from the above two yield ratios. With the increase of sNN, lnkπ from inelastic pp collisions increases obviously and that from nucleus-nucleus (Au-Au and Pb-Pb) collisions slightly decrease. The dependence of lnkπ on 1/sNN can also be empirically described by the following functions of
(12)lnkπpp=−2.859±0.814sNN+−0.016±0.045,lnkπNN=2.890±1.084sNN2+0.205±0.142sNN+−0.010±0.016,with χ2/dof to be 0.543/8 and 56.886/13, respectively, where kπpp and kπNN represent the kπ from pp and nucleus-nucleus (Au-Au and Pb-Pb) collisions, respectively. In fact, we also can use a linear function to fit the energy dependent kπNN roughly, but in order to more accurately describe the data points in low energy region, we adopt the above polynomial function. The fitting results are represented by the solid and dotted curves, where the solid curves correspond to the data points in the energy range mentioned above, and the dotted curves show the changed trends of data points. It is noticed that the values of intercepts of the above four curves are asymptotically 0, which means the limiting values of the three yield ratios are one at very high energy.Figure 12
Energy-dependent and modified yield ratios. The symbols denote the yield ratios modified by removing strong decay from high-mass resonance and weak decay from heavy flavor hadrons. The curves are the fits according to Equations (11)–(14), respectively.Based on the corrected yield ratios of negative to positive particles and Equations (8) and (9), the energy-dependent chemical potentials, μπ, μK, and μp, of light hadrons, π, K, and p, and the chemical potentials, μu, μd, and μs, of light quarks, u, d, and s, are obtained and are shown in Figure 13 with different symbols. The black squares, red circles, and blue triangles denote the calculated results from inelastic pp, central Au-Au, and central Pb-Pb collisions at mid-rapidity, respectively. The curves are the derivative results according to Equations (11)–(14) corresponding to the fitted curves in Figure 12. The red curves in Figures 13(a), 13(d), 13(e), and 13(f) are the derivative results related to kπNN from central nucleus-nucleus collisions, and the black curves are the derivative results related to kπpp from inelastic pp collisions or other yield ratios. The solid and dotted curves in Figure 13 correspond to the solid and dotted curves in Figure 13, respectively. In fact, the μp in this work is close to the μB extracted from thermal fits [34, 70, 71], which indicates to some extent that our calculation of the chemical potential is correct. One can see that with the increase of sNN from AGS to LHC, μπ from central nucleus-nucleus collisions increases obviously and that from inelastic pp collisions decreases obviously, while μK, μp, μu, μd, and μs from both central nucleus-nucleus and inelastic pp collisions decrease obviously. At the same energy, μK is larger than μπ but less than μp, and μu is almost as large as μd but larger than μs due to the difference in mass between different particles. The limiting values of the six types of chemical potentials from central nucleus-nucleus and inelastic pp collisions are zero at very high energy.Figure 13
Energy-dependent chemical potentials, (a)μπ, (b) μK, and (c) μp, of light hadrons, π, K, and p, and (d) μu, (e) μd, and (f) μs, of light quark, u, d, and s. The symbols denote the calculated results according to the modified yield ratios and Equations (8) and (9), and the curves are the derivative results based on the fits of Equations (11)–(14) in Figure 12.From Figure12, one can see that in central Au-Au and Pb-Pb collisions, lnkπNN>0 increases and lnkK<0 and lnkp<0 decrease obviously with the increase of energy. The difference between lnkπNN and lnkK (lnkp) is caused by different production mechanisms. In the processes of producing pions, kaons, and protons, the difference of cross section of absorption, content of primary proton in nuclei, and so on can result in the difference in the yield of these particles. lnkπNN>0 and lnkπpp<0 result in μπNN>0 and μπpp<0. As the energy increases to the LHC, in both central Au-Au (Pb-Pb) collisions and pp collisions, kπ, kK, and kp approach to one, and μπ, μK, and μp approach to zero. These same limiting values indicate that hard scattering process possibly plays an important role, mean-free-path of particles becomes large, and the collision system possibly changes completely from the hadron-dominant state to the quark-dominant state.In Figure13, it should be noted that the derived curves of hadron and quark chemical potentials from the fits of the energy-dependent yield ratios in Figure 12 simultaneously show a maximum at around 4 GeV, which is not observed from the fits of yield ratios. In order to figure out the accurate energies at these maximums, we make the following calculation according to Equations (6), (8), (9), (11)–(14). For these black curves, whose derivation does not involve Equation (14), we figure out the analytical solutions for the energies at these maximums. The derivation is as follows. From the fits of the logarithms of the three types of yield ratios with 1/sNN in Equations (11)–(13), one can see that all the intercepts on the vertical axes approximate to zero. For simplicity of calculation, we assume that all the intercepts are zero; then, Equations (11)–(13) have the following form of
(13)lnki=AisNNi=π,K,andp,where Ai is the slope. Then, according to Equations (8) and (9), the chemical potential μj can be given by
(14)μj=TchBjsNNj=π,K,p,u,d,ands,where Bj is a constant. In consideration of Tch given by Equation (6), μj can be written as
(15)μj=TlimBj1+exp2.60−lnsNN/0.45SNN−1.Then,(16)dμjdSNN=TlimBj1+exp2.60−lnsNN/0.45SNN−21−119exp2.60−lnsNN/0.45.Letdμj/dsNN=0, then,
(17)1−119exp2.60−lnsNN/0.45=0.Finally, we obtain the energy value at the maximum,sNN=3.526 GeV. One can see that all the energy values at these maximums are the same, which means that in the case of linear fitting (Equations (11)–(13)), the energy is independent of the slope parameter of linear equation. It needs to be emphasized that in Figure 12, due to the lack of data in low energy region, we can only make prediction about them by the linear fits that can well describe the data points in relatively high energy region. That means the obtained energy value at the maximum is based on the above linear fits. For these red curves, whose derivation involves Equation (14), due to the calculation being complicated, we only give the numerical solutions for the energies at the maximums of these curves. From the data of these derived curves obtained by Equations (6), (8), (9), (11)–(13) in Figure 13, we give the numerical solutions for the energies at the maximums of all curves. From Figures 13(a)–13(f), the values of the energies at these maximums one by one are 3.584 GeV (for black curve), 3.398 GeV (for red curve), 3.310 GeV, 3.534 GeV, 3.537 GeV (for black curve), 3.590 GeV (for red curve), 3.527 GeV (for black curve), 3.428 GeV (for red curve), 3.521 GeV (for black curve), and 3.666 GeV (for red curve). One can see that the energies at these maximum range from 3.310 GeV to 3.666 GeV, and the average value of these energies is 3.510 GeV.The special energy (around 3.510 GeV) at the maximum possibly is the critical energy of phase transition from a liquid-like hadron state to a gas-like quark state in the collision system, where the liquid-like state and the gas-like state are the states in which the mean-free-path of interacting particles are relatively short and relatively long, respectively. In other words, at this special energy, the collision system starts to change initially its state from the liquid-like nucleons and hadrons to the gas-like quarks, and many properties of the system also change. The curve of proton chemical potential having the maximum at this special energy, indicates that the density of baryon number in nucleus-nucleus collisions has the largest value and the mean-free-path of particles has the smallest value at this special energy, which means that the hadronic interactions play an important role at this stage [72]. When sNN>3.510GeV, the chemical potential gradually decreases with increasing sNN, implying that the density of baryon number gradually decreases [72], shear viscosity over entropy density gradually weakens [73], and mean-free-path gradually increases. Meanwhile, the hadronic interactions gradually fade and the partonic interactions gradually become greater. When sNN increases to the top RHIC, especially the LHC, all types of chemical potentials approach to zero, when the high energy collision system possibly changes completely from the hadron-dominant state to the quark-dominant state and signifies that the partonic interactions possibly play a dominant role at the top RHIC and LHC [1, 15], and the strongly coupled QGP (sQGP) has been observed [24–26]. It should be pointed out that, due to the lack of experimental data in low energy range, the existence of this maximum is not actually certain. The maximum is only a calculated result according to some empirical formulas, so the energy of critical point is large fluctuations. Although the trend of the chemical potential in low energy region (sNN<3.510GeV) is unlikely, the maximum point is at least a turning point, which implies the possibility of phase transition. These results are consistent with our previous work [28]. Our result (around 3.510 GeV) of the critical energy of phase transition is consistent with that (below 19.6 GeV) by the STAR Collaboration [1] and less than the result (between 11.5 GeV and 19.6 GeV) of a study based on the correlation between collision energy and transverse momentum [13–15] and the result (around 62.4 GeV) of the study based on a striking pattern of viscous damping and an excitation function [16]. One can see that although there are many study results, none of them have been confirmed to be reliable. Therefore, we need to continue to study and confirm the exact critical energy of the phase transition.


## 4. Summary and Conclusion

The transverse momentum (or mass) spectra of final-state light flavor hadrons,π±, K±, p, and p¯, produced in central Au-Au, central Pb-Pb, and inelastic pp collisions at mid-rapidity over an energy range from AGS to LHC, are described by a two- or one-component Erlang distribution in the frame of multisource thermal model. The fitting results are in agreement with the experimental data recorded by the E866, E917, E895, NA49, NA44, NA61/SHINE, PHENIX, STAR, and ALICE Collaborations.From the fitting parameters, in most cases, the experimental data ofpT (or mT) spectra are suitable for the two-component Erlang distribution, where the first component corresponding to a narrow low-pT (or mT) region is contributed by the soft excitation process in which a few sea quarks and gluons take part in, and the second component corresponding to a wide high-pT (or mT) region is contributed by the hard scattering process which is a more violent collision between two valent quarks in incident nucleons. The study shows that the values of the contribution ratio of soft excitation process are more than 50%, which means that the excitation degrees of these collision systems are mainly contributed by the soft excitation processes.Based on the normalization constants in fitting the transverse momentum or mass spectra of final-state light flavor particles, the final-state yield ratios of negative to positive particles are obtained. The energy-dependent chemical potentials of light hadrons,μπ, μK, and μp, and quarks, μu, μd, and μs, are extracted from the modified yield ratios in which the contributions of strong decay from high-mass resonance and weak decay from heavy flavor hadrons are removed. With the increase of sNN over a range from a few GeV to more than 10 TeV, the μK, μp, μu, μd, and μs decrease obviously in central Au-Au, central Pb-Pb, and inelastic pp collisions, while μπ increases in central Au-Au and Pb-Pb collisions and it decreases in inelastic pp collisions. When collision energy increases to the top RHIC and LHC, all types of chemical potentials are small and the limiting values of them are zero in central Au-Au, central Pb-Pb, and inelastic pp collisions at very high energy.The logarithms of the yield ratios,lnkπpp, lnkK, and lnkp, show obviously linear dependences on 1/sNN in mentioned energy range, and lnkπNN increases as a polynomial function of 1/sNN. Base on the above relationships, we find that at about 3.510 GeV, the derived curves of hadron and quark chemical potentials simultaneously show the maximum. This special energy possibly is the critical energy of phase transition from a liquid-like hadron state to a gas-like quark state in high energy collision system, where the density of baryon number in nucleus-nucleus collisions has a large value and the hadronic interactions play an important role at this stage. When collision energy increases to the top RHIC, especially the LHC, all types of chemical potentials approach to zero, which indicates that high energy collision system possibly changes completely from the hadron-dominant liquid-like state to the quark-dominant gas-like state and the partonic interactions possibly play a dominant role at the LHC.



---

*Source: 1265090-2020-11-05.xml*