# On Productions of Net-Baryons in Central Au-Au Collisions at RHIC Energies

Authors: Ya-Hui Chen; Guo-Xing Zhang; Fu-Hu Liu
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

The transverse momentum and rapidity distributions of net-baryons (baryons minus antibaryons) produced in central gold-gold (Au-Au) collisions at 62.4 and 200 GeV are analyzed in the framework of a multisource thermal model. Each source in the model is described by the Tsallis statistics to extract the effective temperature and entropy index from the transverse momentum distribution. The two parameters are used as input to describe the rapidity distribution and to extract the rapidity shift and contribution ratio. Then, the four types of parameters are used to structure some scatter plots of the considered particles in some three-dimensional (3D) spaces at the stage of kinetic freeze-out, which are expected to show different characteristics for different particles and processes. The related methodology can be used in the analyses of particle production and event holography, which are useful for us to better understand the interacting mechanisms.


---

## Body

## 1. Introduction

In high energy nucleus-nucleus collisions, due to the fast speed and the great mass of the projectile and target nuclei, the collisions can form new substances such as pions, kaons, antiprotons, which do not obviously exist in the nuclei [1, 2]. In the area filled with new and rich phenomena, there are many issues which are difficult to deal with, and the collisions still stay at the phenomenological stage in understanding [3]. To restore the original process of the collisions by analyzing various properties of final-state particles is a major aspect in high energy nuclear physics [4].The relativistic heavy ion collider (RHIC) located in the United States and the large hadron collider (LHC) located in Switzerland are the main sources of experimental data of heavy-ion collisions [5–7]. They play very significant roles in the studies of high energy nucleus-nucleus collisions. The transverse momentum and rapidity distributions of different final-state particles can be measured from the experiments. The two distributions are related to the excitation degree of the interacting system and the penetrating power of the projectile nucleus, respectively. From the two distributions, we can extract some other distributions and correlations. Scientists worked in fields of high energy and nuclear physics are very interested in experiments performed in the colliders.Since the 1980s or earlier days, many phenomenological models have been proposed to explain the large number of experimental data in high energy particle-particle, particle-nucleus, and nucleus-nucleus collisions [4]. In these models, except for the standard (Boltzmann, Fermi-Dirac, and Bose-Einstein) and other distributions, the Tsallis statistics can describe the thermal source (fireball) and transverse momentum distributions [8–10]. An effective temperature parameter which reflects the excitation degree of interacting system and an entropy index (nonextensive parameter) which describes the degree of nonequilibrium of interacting system are used in the Tsallis statistics. Particularly, the Tsallis statistics describes the temperature fluctuations in two- or three-standard distributions.Among abundant experimental data, the transverse momentum and (pseudo)rapidity distributions of identified particles are very important. From the two distributions, we can extract some parameters such as the (effective) temperature, entropy index, rapidity shift, and contribution ratio in the framework of a multisource thermal model [11–14]. These parameters can be then used to extract some scatter plots of the considered particles in some three-dimensional (3D) spaces. These scatter plots are expected to show different characteristics for different particles and processes. We then propose a new method for studies of particle production and event holography. In addition, baryon transport is a fundamental aspect of high energy nucleus-nucleus collisions as it pertains to the initial scatterings where the Quark-Gluon Plasma (QGP) is formed. It is a process that is not very well understood and so further theoretical studies are needed.In this paper, in the framework of a multisource thermal model [11–14] in which each source is described by the Tsallis (transverse) momentum distribution [8–10, 15, 16], we analyze the transverse momentum and rapidity distributions of net-baryons (baryons minus antibaryons) produced in central gold-gold (Au-Au) collisions at the center-of-mass energies of 62.4 GeV and 200 GeV. Based on the descriptions of the two distributions, some scatter plots of the considered particles in some 3D spaces are structured to show partly the event picture in the whole phase space. A methodology related to particle production and event holography is presented to better understand the interacting mechanisms.


## 2. The Model and Method

According to the multisource thermal model [11–14], in high energy nucleus-nucleus collisions, the target nucleus and the projectile nucleus penetrate each other. We can assume that many emission sources are formed in the interacting system. In a given frame of reference such as the center-of-mass reference frame, these sources with different rapidity shifts yx in the rapidity space can be divided into four parts: a leading target nucleon cylinder (LT) in rapidity interval [yLTmin,yLTmax], a leading projectile nucleon cylinder (LP) in rapidity interval [yLPmin,yLPmax], a target cylinder (TC) in rapidity interval [yTmin,0], and a projectile cylinder (PC) in rapidity interval [0,yPmax]. It is expected that the LT and TC are mainly contributed by the target nucleus, and the LP and PC are mainly contributed by the projectile nucleus. Because the TC and PC are linked at midrapidity, they can be regarded as a whole central cylinder. For symmetric collisions such as Au-Au collisions which are studied in the present work, the equations yTmin=-yPmax, yLTmax=-yLPmin, and yLTmin=-yLPmax can be used. The free rapidity shifts are then yPmax, yLPmin, and yLPmax.In the rest frame of each source, the source is assumed to emit final-state particles isotropically, and the rapidity of a considered particle isy′. Then, the rapidity of the considered particle in the center-of-mass reference frame is(1)y=yx+y′.According to different origins of different final-state particles, yx can be in [yLTmin,yLTmax], [yLPmin,yLPmax], [yTmin,0], or [0,yPmax] corresponding to the LT, LP, TC, and PC, respectively. Generally, y′ can be obtained by the definition of rapidity, and yx can be obtained due to different contribution ratios of the four sources.In the Tsallis statistics [8–10, 15, 16], let T denote the temperature parameter of the considered source, q denote the nonextensive parameter (entropy index), N denote the number of particles, and p′, m0, and μ denote the momentum, rest mass, and chemical potential of the considered particle, respectively. We have the final-state particle momentum distribution in the rest frame of the considered source to be(2)fp′=1NdNdp′=Cp′21+q-1Tp′2+m02-μ-q/q-1,where C is the normalization constant which gives (2) to be a probability distribution. Generally, the chemical potential μ can be neglected due to its small value at high energy. Then, (2) can be simplified as(3)fp′=1NdNdp′=Cp′21+q-1Tp′2+m02-q/q-1.In the Monte Carlo method, we can obtain the momentump′, the energy E′=p′2+m02, the polar angle θ′, the azimuth ϕ′=ϕ, the transverse momentum pT′=p′sin⁡θ′=pT, the x-component of momentum px′=pTcos⁡ϕ=px, the y-component of momentum py′=pTsin⁡ϕ=py, the longitudinal momentum pz′=p′cos⁡θ′, the rapidity y′=0.5ln⁡[(E′+pz′)/(E′-pz′)], the rapidity shift yx, the energy E=pT2+m02cosh⁡y, the longitudinal momentum pz=pT2+m02sinh⁡y, the transverse velocity βT=pT/E, the x-component of the velocity βx=βTcos⁡ϕ, the y-component of the velocity βy=βTsin⁡ϕ, the longitudinal velocity βz=pz/E, and so forth, where the quantities with no mark denote those in the center-of-mass reference frame.By using the above method, we can get many distributions and scatter plots for a given collision process. Generally,T and q are sensitive to the transverse momentum distribution and insensitive to the rapidity distribution. We can use the transverse momentum distribution to determine values of T and q. The rapidity shifts and contribution ratios of different sources are usually determined by the rapidity distribution. The contributions of the TC and PC are mainly in the central rapidity region, and the contributions of the LT and LP are mainly in the backward and forward rapidity regions, respectively, where the latter two regions are also called the fragmentation regions. In the present work, to describe the rapidity distribution as accurate as possible, we may use T and q obtained from pT distribution at y=0 for the TC and PC and T and q obtained from pT distribution at y=3.0 (or 2.9) for the LT and LP.Particularly, for the considered Au-Au collisions at 62.4 and 200 GeV in the present work, we can describepT distribution obtained by the above Monte Carlo method by the analytic expression [10](4)fpTpT=1NdNdpT=CTpTpT2+m02∫ymin⁡ymax⁡cosh⁡y·1+q-1TpT2+m02cosh⁡y-q/q-1dy,where ymin and ymax denote the minimum and maximum rapidities, respectively, and CT is the normalization constant which gives (4) to be a (normalized) probability distribution. At the same time, px (py) distribution obtained by the above Monte Carlo method can be approximately parameterized to(5)fpx,ypx,y=1NdNdpx,y=Cx,ypx,y2+m02∫ymin⁡ymax⁡cosh⁡y·1+q-1Tpx,y2+m02cosh⁡y-q/q-1dy,where Cx,y is the normalization constant which gives (5) to be a probability distribution. In the present calculation, we can fit a given pT distribution by using the Monte Carlo method or analytic expression to obtain values of free parameters T and q. It should be noted that the given central y in (4) and (5) should be shifted to 0 when we extract T and q because the effect of vectored longitudinal motion has to be subtracted.


## 3. Comparisons and Extractions

Figure1 shows the transverse momentum distributions of net-baryons produced in central Au-Au collisions at (a)–(d) 62.4 and (e)–(h) 200 GeV. From Figures 1(a)–1(d), the corresponding rapidities are y = 0, 0.65, 2.3, and 3.0, respectively. From Figures 1(e)–1(h), the corresponding rapidities are y = 0, 0.9, 1.9, and 2.9, respectively. The symbols represent the experimental data of the BRAHMS Collaboration [17–19] measured at the RHIC and collected together in [20] in which the longitudinal axis should not be d2N/dyd2pT, but 2πpT-1d2N/dydpT. The solid curves are the Tsallis fitting results by us. The values of free parameters T and q obtained by using the analytic expression, normalization constant N0 obtained by fitting the experimental data, and χ2/ndf (number of degree of freedom) obtained in the fit are given in Table 1, where the last three points in Figure 1(h) have not been included in the calculation of χ2/ndf due to the low statistics.Table 1
Values of parameters andχ2/ndf corresponding to the fits in Figure 1, where Figures 1(a)–1(d) are for central Au-Au collisions at 62.4 GeV and Figures 1(e)–1(h) are for central Au-Au collisions at 200 GeV. The unit of T is GeV. The last three points in Figure 1(h) have not been included in the calculation of χ2/ndf due to the low statistics.
Figure
y
T
q
N
0
χ
2
/
ndf
Figure1(a)
0
0.270 ± 0.013
1.023 ± 0.002
3.45 ± 0.17
0.757
Figure1(b)
0.65
0.230 ± 0.012
1.023 ± 0.002
5.13 ± 0.26
0.169
Figure1(c)
2.3
0.181 ± 0.009
1.023 ± 0.002
10.80 ± 0.54
0.546
Figure1(d)
3.0
0.158 ± 0.008
1.023 ± 0.002
14.18 ± 0.71
0.209
Figure1(e)
0
0.282 ± 0.014
1.023 ± 0.002
1.61 ± 0.08
0.764
Figure1(f)
0.9
0.239 ± 0.012
1.023 ± 0.002
1.46 ± 0.07
0.464
Figure1(g)
1.9
0.223 ± 0.011
1.023 ± 0.002
3.39 ± 0.17
0.581
Figure1(h)
2.9
0.210 ± 0.011
1.023 ± 0.002
7.79 ± 0.39
0.571Figure 1
Transverse momentum distributions of net-baryons produced in central Au-Au collisions at different rapidities shown in the panels. (a)–(d) and (e)-(f) are for 62.4 and 200 GeV, respectively, where the symbols represent the experimental data of the BRAHMS Collaboration [17–19] measured at the RHIC and collected together in [20], and the curves are our modeling results.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)From Figure1 and Table 1, one can see that T decreases obviously with the increase of y, and q has a constant value with y. This reflects nonequivalent excitations of the sources of net-baryons with different rapidities. The net-baryons with low rapidity are produced in the sources located in the central rapidity region, which are in fact in the TC/PC, for which the degree of excitation is high and the longitudinal speed is low. The net-baryons with high rapidity are produced in the sources located in the fragmentation regions, which are in fact in the LT/LP, for which the degree of excitation is low and the longitudinal speed is high. The fact that the effective temperature decreases with increasing of rapidity renders also that the velocity of radial flow in the fragmentation region is less than that in the central region. In the same region, compared with 62.4 GeV, the interacting system at 200 GeV has higher excitation degree and larger radial flow. In all rapidity region, both the interacting systems at 62.4 and 200 GeV are close to equilibrium state due to small values of q.Figure2 presents the rapidity distributions of net-baryons produced in central Au-Au collisions at (a) 62.4 and (b) 200 GeV. The solid symbols represent the experimental data of the BRAHMS Collaboration [17–19] and collected together in [20], and the open symbols are reflected from the data around midrapidity, where [19] is a conference proceeding which means that most points in Figure 2(b) right at y>1 are preliminary only. The solid curves are our fitting results. From left to right, the dotted or dashed curves presented in each panel are orderly the contributions of the LT, TC, PC, and LP. The maximum values of the net-baryons calculated by us are around the beam rapidities (4.2 at 62.4 GeV and 5.4 at 200 GeV) shown in the figure by the arrows. A few net-baryons extend beyond the beam rapidities due to the forward excitation of nucleons in the beam nuclei and the statistical fluctuations in the Monte Carlo calculation. Other particles have probabilities to extend beyond the beam rapidity [21]. In the calculation, we have used T and q obtained from pT distribution at y = 0 for the TC/PC and T and q obtained from pT distribution at y = 3.0 (or 2.9) for the LT/LP. Other free parameters such as yPmax, yLPmin, and yLPmax and the contribution ratio k of LP are shown in Table 2 with the normalization constant N0 and χ2/ndf. By comparing the modeling results with the experimental data, we find that the experimental rapidity distributions of net-baryons can be well described by the model. Both of the rapidity spans yPmax-0 and yLPmax-yLPmin in central Au-Au collisions at 200 GeV are greater than those at 62.4 GeV. These trends are natural results.Table 2
Values of parameters andχ2/ndf corresponding to the curves in Figures 2(a) and 2(b), where T and q are taken from Table 1 for the central cylinders at y=0 and for the leading nucleons at y=3.0 (or 2.9). See [14] on the discussions for the number of free parameters and the ndf.
Figure
y
Pmax
⁡
y
LPmin
⁡
y
LPmax
⁡
k
N
0
χ
2
/
ndf
Figure3(a)
3.60 ± 0.18
1.90 ± 0.10
3.55 ± 0.18
0.21 ± 0.01
297 ± 15
0.339
Figure3(b)
4.35 ± 0.22
2.32 ± 0.12
4.55 ± 0.23
0.32 ± 0.02
256 ± 13
0.957Figure 2
Rapidity distributions of net-baryons in central Au-Au collisions at (a) 62.4 and (b) 200 GeV. The solid symbols represent the experimental data of the BRAHMS Collaboration [17–19] measured at the RHIC and collected together in [20], and the squares and diamonds are for 0–5% and 0–10% centralities, respectively. The open symbols are reflected from the data around midrapidity. The curves are our calculated results, where the dotted and dashed curves from left to right are the contributions of LT, TC, PC, and LP, respectively, and the solid curves are the sum of the four sources. Arrows indicate the beam rapidities. See text on the discussion of a few net-baryons to extend beyond the beam rapidities.
(a)
(b)We would like to point out that the normalization constant for Figure2 should be known (it should be the number of participant baryons Npart), but not a free parameter due to baryon number conservation. In fact, the situation is complex, and we have to regard it as a free parameter. As we know, except for the limitation of experimental acceptance range, some baryons in the participant region are not real participants due to the limited scattering cross-section. These no-participants are in fact the spectator baryons. In addition, some baryons in 0–5% and 0–10% centralities discussed in the present work are real spectator baryons. These spectator baryons move straight forward along beam direction and should be excluded in the analysis. Because we do not know exactly the number of spectator baryons, we have to regard the normalization constant as a free parameter. Generally, N0<Npart=357±8 [18]. Table 2 confirms this situation.The scatter plots of 1000 net-baryons in 3D spacepx-py-y in central Au-Au collisions at (a) 62.4 and (b) 200 GeV are given in Figure 3. The black, magenta, olive, and blue balls represent orderly the contributions of LT, TC, PC, and LP from low to high in the rapidity space due to different contribution ratios. Similar to Figure 3, the scatter plots of 1000 net-baryons in 3D (momentum) space px-py-pz and in 3D (velocity) space βx-βy-βz in central Au-Au collisions at (a) 62.4 and (b) 200 GeV are shown in Figures 4 and 5, respectively. The meanings of balls with different colors are the same as those in Figure 3.Figure 3
Scatter plots of 1000 net-baryons in 3Dpx-py-y space in central Au-Au collisions at (a) 62.4 and (b) 200 GeV. The black, magenta, olive, and blue balls represent orderly the contributions of LT, TC, PC, and LP from low to high in the rapidity space due to different contribution ratios.
(a)
(b)Figure 4
Same as Figure3, but showing the scatter plots in 3D momentum space.
(a)
(b)Figure 5
Same as Figure3, but showing the scatter plots in 3D velocity space.
(a)
(b)One can see different scatter regions of the four sources at the stage of kinetic freeze-out. The neighbor scatter regions overlap partly iny direction in 3D space px-py-y, in pz direction in 3D space px-py-pz, and in βz direction in 3D space βx-βy-βz. The distribution range of pz is much larger than that of px and py. The distribution pattern in βz direction is very different from that in βx and βy directions. The particle scatter plots are bumpy cylinders in 3D spaces px-py-y and px-py-pz. The contributions of TC/PC are mainly in the regions with small y and pz, while the contributions of LT/LP are mainly in the regions with large y and pz. The density in small px and py region is greater than that in other regions. Meanwhile, the particle scatter plot is a bumpy sphere in the velocity space. The densities close to the maximum βz are greater than those in other regions. Particularly, in 3D spaces px-py-y and px-py-pz, the higher the center-of-mass energy is, the broader the scatter plots are; and in the velocity space, the higher the center-of-mass energy is, the higher the densities close to the maximum βz are.


## 4. Discussions

In the central rapidity region, baryons are produced together with antibaryons and the net-baryon number is small, where these baryons are produced mainly from gluons and sea quarks. We can use the TC/PC to describe baryon productions in the central rapidity region. In the backward/forward rapidity regions, there are almost only baryons and no antibaryons, where these baryons are produced from the valence quarks of the target/projectile nuclei. We can use the LT/LP to describe baryon productions in the backward/forward regions. The effective temperature parameter extracted from the Tsallis statistics decreases with increasing the rapidity, which makes the interactions of gluons and sea quarks in the central rapidity region lose more energy compared with the interactions of valence quarks in the forward rapidity region. The lost energy is expected to transform mainly to the energy of thermal motions. The more the energy lost, the higher the temperature was.We would like to point out that the discussion of transverse momentum distributions in Figure1 does not address radial flow. It means that we have obtained the effective temperature, but not the source temperature [22]. It is known that radial flow plays a huge role for transverse momentum distributions, especially for heavy hadrons such as protons. Normally, to differentiate between the source temperature and radial flow, one does a combined blast wave fit of pions, kaons, and protons. The source temperature describes the thermal motion, and the radial flow describes the blast wave effect. Although the present work does not address particularly the radial flow, the extracted parameter T which we call effective temperature contains together the thermal motion and radial flow (blast wave) effect. This means that the succeeding analyses (Figures 2–5) have not been affected, whether we address particularly the radial flow or not. The only regretful thing is that we do not know the source temperature and radial flow velocity, respectively.In our model, to extract the source temperatureT0 [22] and average transverse velocity uT of radial flow, we can extract the effective temperatures T and mean transverse momentums pT from transverse momentum distributions of pions, kaons, protons, and other light flavor hadrons, respectively. On the planes T-m0γ-  (or T-m0 [22]) and pT-m0γ-, reasonable linear fittings should give T=T0+am0γ- (or T=T0+a′m0) and pT=pT0+uTm0γ-, where γ- is the mean Lorentz factor for a given particle, a (or a′) is a coefficient related to uT, pT0 is the mean transverse momentum for massless particle, and T0 and uT can be obtained by us from the linear relations. It is obvious that we need the transverse momentum distributions for more than three types of particles to do the linear fittings. It is beyond the focus of the present work.The entropy index shows the same value from the central to forward rapidity regions. This renders the nonequilibrium degrees for different sources of net-baryons in the interacting system nearly the same. Although the central rapidity region corresponds to higher source temperature and the forward rapidity region corresponds to higher baryon density, the central rapidity region does not show more advantages or disadvantages compared with the forward rapidity region in the aspect of nonequilibrium degree for the productions of net-baryons. The higher temperature in the central rapidity region is in fact affected by produced particles which carry more thermal motion energies, and the higher baryon density in the forward rapidity region is caused by the leading nucleons which contribute directly to the baryons. Because the values of entropy index are small, the interacting system is close to the equilibrium state.In Figures3 and 4, we have shown px and py separately. Although the data we used only have pT, it does not naturally mean that px and py have the same distribution. In fact, the anisotropic flow is generally covered in pT distribution. If we consider the anisotropic flow such as the elliptic flow, px and py will have different distributions. To extract the anisotropic flow, we need azimuthal information. It is regretful that we neglected this issue in the present work for the purpose of convenience. Although we neglected the anisotropic flow, we still separate px and py and make 3D plots instead of 2D plots for presenting the methodology. Generally, px distribution can be measured in the direction in the reaction plane, and py distribution can be measured in the direction out of the reaction plane. In the case of considering the anisotropic flow, px and py distributions will have different values of parameters. This also explains the reason that the 3D plots in momentum space are different from the actual measurements, d2N/dydpT, that are reported in experiments in 3D momentum space, and the anisotropic flow is covered in them.In Figure5, we have shown a bumpy sphere in the velocity space for presenting partly the interacting event in a whole space image. It seems that there is symmetry between transverse and longitudinal directions. In fact, there is no symmetry due to the density distributions of scatter points in the two directions being totally different. Two small high density areas appear in the regions of large |βz| and small |βx,y|, and a large low density area appears in the region of small |βz| and various |βx,y|. It is normal that a few scatter points have large |βx,y| at RHIC energies due to some particles in Figure 1 having large pT (~2 GeV/c). These large pT result in βT being around 0.9c at y=0. Early Lund string model and experimental data gave a pT of around a few hundred MeV/c in hadronic collisions at low and intermediate energies [23], which results in βT being around 0.5c if pT=600 MeV/c and y=0.Although Lambda decays affect the rapidity distribution somewhere, we have analyzed the net-baryons but not the net-protons measured by the BRAHMS Collaboration [17–19] due to more abundant data of net-baryons being collected in [20]. It is more convenient for us to use [20] instead of [17–19]. In addition, the typical motivation to study net-baryons is that one expects it to be transported in the initial collisions and the baryon number is conservative in the final state. This is why the study of net-baryons is more interesting. Generally, by using the method used in the present work, we can analyze scatter plots for different particles. As a premise, the transverse momentum and rapidity distributions have to be known. If the identified particles such as π+, π-, K+, K-, p, p-, and others produced in collisions at different center-of-mass energies are satisfied the premise, we can structure the corresponding scatter plots for them, respectively. Then, we can compare these scatter plots and related root-mean-square quantities carefully. It is expected that some dependences of scatter plots and root-mean-square quantities on particle mass and center-of-mass energy can be obtained, and one or two points of inflections in the dependences can be found, which may imply the phase transition from the hadronic matter to QGP. The present work is only a presentation for the methodology, but not the study of dependent relation itself.As a thermal and statistical model, the multisource thermal model is successful in the descriptions of transverse momentum and rapidity distributions of net-baryons produced in nucleus-nucleus collisions at high energies. From these descriptions, we can extract some parameters and get some scatter plots of net-baryons in some spaces. Except for the transverse momentum and rapidity distributions, the multisource thermal model is successful in the descriptions of multiplicity distributions [13]. In addition, we can use other distributions to replace the Tsallis distribution. These distributions include, but are not limited to, the Boltzmann distribution, the Fermi-Dirac (Bose-Einstein) distribution, the Erlang distribution, and their combination of multiple components. Particularly, these multicomponent distributions reflect temperature fluctuations which can be described by the Tsallis distribution which uses the entropy index to describe the degree of nonequilibrium.More deeply, a phenomenological study on the net-baryon production in the Color Glass Condensate (CGC) formalism is performed in references [20, 24]. The transverse momentum and rapidity distributions of net-baryons are studied by the CGC formalism and related parameters are extracted [20, 24]. According to the related parameters, it is expected that the scatter plots of net-baryons in some spaces can be obtained. We notice that the maximum rapidities of net-baryons produced in central Au-Au collisions at RHIC energies obtained in the present work are slightly greater than those by the CGC formalism [20]. The shapes of the two distributions have also some differences.In other works [25, 26] which are also in the CGC framework and based on small-coupling quantum chromodynamics (QCD), the rapidity distributions of net-baryons are studied by a gluon saturation model. It is shown that geometric scaling is exhibited in the distributions [25, 26]. We notice that the maximum rapidities of net-baryons obtained in the present work are slightly greater than those by the saturation model [25, 26], too. In addition, the shapes of the two distributions have also some differences.


## 5. Conclusions

From the above calculations and discussions, we obtain following conclusions.(a)
The transverse momentum distributions of net-baryons produced in central Au-Au collisions at 62.4 and 200 GeV are described by the Tsallis distribution which is naturally considered in the multisource thermal model. The effective temperature parameter decreases with the increase of the rapidity. The entropy index does not show an obvious change trend with the increase of the rapidity.(b)
The rapidity distributions of net-baryons produced in central Au-Au collisions at 62.4 and 200 GeV are studied in the framework of the multisource thermal model in which the contributions of leading target/projectile nucleon sources and target/projectile cylinders are included. Each source is described by the Tsallis statistics with the parameter values obtained from the transverse momentum distributions. The rapidity spans for the cylinders and the rapidity shifts for the leading nucleons in central Au-Au collisions at 200 GeV are greater than those at 62.4 GeV.(c)
Based on the parameter values obtained from the transverse momentum and rapidity distributions, a series of quantities is extracted by using the Monte Carlo method. As examples, the particle scatter plots at the stage of kinetic freeze-out in three types of 3D spaces are obtained. In 3D spacespx-py-y, px-py-pz, and βx-βy-βz, one can see different scatter regions of the four sources. The neighbor scatter regions overlap partly in beam direction. The particle scatter plot is a bumpy cylinder in the momentum space and a bumpy sphere in the velocity space.(d)
Particularly, the distribution range ofpz is much larger than that of px and py, and the density in small px and py region is greater than that in the other regions. The distribution pattern in βz direction is very different from that in βx and βy directions, and the densities close to the maximum βz are greater than those in the other regions. In 3D spaces px-py-y and px-py-pz, the higher the center-of-mass energy is, the broader the scatter plots are; and in the velocity space, the higher the center-of-mass energy is, the higher the densities close to the maximum βz are.



---

*Source: 614090-2015-06-29.xml*