Abstract

By using the method of data-driven reanalysis, the midrapidity transverse momentum spectra of charged hadrons (, , and ) produced in central and peripheral gold-gold (Au-Au) collisions from the Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider (RHIC) are fitted by using the blast-wave model with the Boltzmann-Gibbs statistics. The model results are in agreement with the experimental data measured by the STAR Collaboration at the RHIC-BES energies. We observe that the kinetic freeze-out temperature , transverse flow velocity , mean transverse momentum , and initial temperature increase with collision energy as well as with event centrality.

1. Introduction

One of the most fundamental questions in high energy and nuclear physics is to determine the phase structure of the strongly interacting quantum chromodynamics (QCD) matter [13]. The yield ratios, transverse momentum spectra, and other data for various identified particles produced in proton-proton , proton-nucleus (), and nucleus-nucleus () collisions at high energies are important observable quantities for determining the phase structure. The experimental facilities such as the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) provide excellent tools to study the properties of Quark-Gluon Plasma (QGP) [46], which are expected to create collision events with high multiplicities.

The phase diagram of the QCD matter is usually expressed in terms of the chemical freeze-out temperature () and the baryon chemical potential () [7, 8]. Besides, other quantities such as the kinetic freeze-out temperature ( or ) and transverse flow velocity () are useful to understand the phase diagram [9]. To search for the possible critical energy in the phase transition from hadronic matter to QGP in high-energy collisions, the STAR Collaboration has been performing the Beam Energy Scan (BES) program [1013] at the RHIC. Besides, other experiments at similar or lower energies at other accelerators are scheduled [14, 15].

Generally, the processes of high-energy collisions result possibly in three main stages [1618]: (i)The initial stage: at this stage, the collisions are in the beginning. The temperature at this stage is called the initial temperature which is one of the main factors to affect the particle spectra, which is less studied in the community comparatively. After the initial state, the “fireball” leads to a decrease in the temperature and finally to the hadronization(ii)The chemical freeze-out stage: at this stage, the inner collisions among various particles are elastic and the yield ratios of differential types of particles remain invariant. The chemical freeze-out temperature can be obtained from the particle ratios, which is much studied in the community comparatively(iii)The kinetic freeze-out stage: at this stage, the scattering processes stop, the hadrons decouple from the rest of the system, and the hadron’s energy/momentum spectra freeze in time. The temperature at this stage is known as the kinetic freeze-out temperature which can be obtained from the spectra

When one studies from the spectra, the effect of should be eliminated. If the effect of is not eliminated in the temperature, this temperature is called the effective temperature ( or ). At the stage of kinetic freeze-out, and are two important parameters which describe the thermal motion of the produced particles and the collective expansion of the emission source, respectively. The spectra in a low- region () which is mainly contributed by the soft excitation process essentially separate the contribution of the thermal motion and the collective expansion, if one only extracts and . The spectra in a high- region are contributed by the hard scattering process which is not needed in extracting and .

We are very interested in the extraction of and in collisions at the RHIC-BES energies which are very suitable to study the spectra in a low- region, where the spectra in a high- region are not produced due to not too high energies. In this work, the double-differential spectra of charged particle dependences on collision energy and event centrality in gold-gold (Au-Au) collisions are analyzed with the blast-wave model with the Boltzmann-Gibbs statistics by means of data-driven analysis. The model results are compared with the data measured by the STAR Collaboration at the RHIC-BES energies [19, 20].

The remainder of this work consists of The Method and Formalism, Results and Discussion, and Conclusions. We shall describe the remanent parts orderly.

2. The Method and Formalism

Various methods can be used for the extraction of and , e.g., the blast-wave model with the Boltzmann-Gibbs statistics [2123], the blast-wave model with the Tsallis statistics [2426], an alternative method by using the Boltzmann-Gibbs statistics [22, 2733], and the alternative method by using the Tsallis distribution [3339]. In this work, we choose the blast-wave model with the Boltzmann-Gibbs statistics due to its similarity with the ideal gas model in thermodynamics and few parameters. However, these methods only describe the spectra in the low- region. For the spectra in the high- region if available, the Hagedorn function which is known as the inverse power law [40, 41] can be used. We shall discuss these issues in detail as follows.

In general, there are two main processes responsible for the contribution of spectra. They are (i) the soft excitation process which contributes the soft component in the low- region and (ii) the hard scattering process which contributes the hard component in the whole region if one uses the general superposition function or in the high- region if one uses the usual step function.

For the soft component, according to Refs. [2123], the probability density function of the spectra in the blast-wave model with the Boltzmann-Gibbs statistics results in where is the number of particles, is the normalization constant, is the transverse mass, is the rest mass of the considered particle, and are the radial position and the maximum radial position, respectively, and are the modified Bessel functions of the first and second kinds, respectively, is the boost angle, is a self-similar flow profile, is the flow velocity on the surface, and is used in the original form [21]. Particularly, . The parameter is used in different works, e.g., or noninteger in Refs. [22, 24, 42], which corresponds to the centrality from the center to the periphery.

Equation (1) and similar or related functions are not enough to describe the whole spectra. In particular, the maximum reaches up to in collisions at the LHC [43]. Then, one needs other functions such as the Tsallis–Lévy- [44, 45] or Tsallis–Pareto-type function [44, 46] and the Hagedorn function [40, 41] or inverse power law [4749] to the spectra in high- and very high- regions. In this work, the hard component is simply represented by the inverse power law. That is, where and are free parameters and is the normalization constant which is related to the free parameters.

However, the structure of spectra is very complex. In fact, several regions have been observed and analyzed in Ref. [50]. These regions include the first one with , the second one with , and the third one with . Different regions may correspond to different mechanisms. The first region in our discussion is regarded as the region of the soft excitation process, while the second and third regions are regarded as the regions of the hard and the very hard excitation process, respectively. In particular, a special region with is considered due to the resonant production in some cases, and it is regarded as the region of the very soft excitation process.

Generally, the whole region discussed above can be uniformly superposed by two methods: (i) the general superposition in which the contribution regions of different components overlap each other and (ii) the Hagedorn model (the usual step function) [40] in which there is no overlapping of different regions of different components.

Considering , , , and which denote the probability density functions by the soft, hard, very soft, and very hard components, respectively, where and are assumed to be in the form of and , respectively, the unified superposition according to the first method is where is the contribution fraction of the very soft component, while and denote the contributions of the soft and very hard components, respectively.

The step function can be used to structure the superposition according to the Hagedorn model [40]; i.e., where , , , and are the constants which make the interfacing components link to each other perfectly.

Particularly, if the contributions of the very soft and very hard components can be neglected, Equations (3) and (4) are, respectively, simplified to be Further, if the contribution of the hard component at the RHIC-BES energies can be neglected, Equations (5) and (6) are simplified to be the same form:

This work deals with Au-Au collisions at the RHIC-BES energies, for which Equation (7), i.e., Equation (1), is suitable. In the following section, we shall use Equation (1) to fit the experimental data measured by the STAR Collaboration at the RHIC-BES energies [19, 20].

In particular, the mean and the root-mean-square can be expressed, respectively, as due to where denotes the maximum considered by us. In this work, we take .

It should be noted that although only Equation (1) is used in the analysis, we would like to continue to have the statement and formalism for other functions or distributions such as the inverse power law and its superposition with thermal distribution and the discussions on the very soft, hard, and very hard components. In fact, due to the existence of other functions or distributions, the mentioned method of data-driven reanalysis can be used in the spectra in wide coverage, which is not the case in this work. In addition, it is possible to use simultaneously the (very) soft and (very) hard components in other cases which are more universal.

3. Results and Discussion

Figure 1 presents the event centrality-dependent double-differential spectra, , of , , and produced in the midrapidity interval in Au-Au collisions at the center-of-mass energy per nucleon pair at the RHIC-BES, where denotes the rapidity. The symbols represent the experimental data measured by the STAR Collaboration [19], and the curves are our fitting results by using the blast-wave model with the Boltzmann-Gibbs statistics, Equation (1) [2123]. The spectra in centrality classes 0–5%, 5–10%, 10–20%, 20–30%, 30–40%, 40–50%, 50–60%, 60–70%, and 70–80% are scaled by 1, 1/2, 1/4, 1/6, 1/8, 1/10, 1/12, 1/14, and 1/16, respectively. In the fit, the least-square method is used to determine the best values of parameters. The related parameters along with and degree of freedom (dof) are listed in Table 1, where the centrality classes are listed together. One can see that Equation (1) fits well the data in Au-Au collisions at 7.7 GeV at the RHIC.

Figure 2 is the same as Figure 1, but it shows the spectra at . One can see that Equation (1) fits well the data in Au-Au collisions at 11.5 GeV at the RHIC-BES.

Figure 3 is also the same as Figure 1, but it shows the spectra at , where the data are cited from Ref. [20]. Once again, Equation (1) fits well the data in Au-Au collisions at 14.5 GeV at the RHIC-BES.

Figures 46 are also the same as Figure 1, but they show the spectra at , 27, and 39 GeV, respectively. Once more, Equation (1) fits well the data in Au-Au collisions at other RHIC-BES energies.

It is noteworthy to point out that Equation (1) for the blast-wave model in the system is assumed to be in local thermodynamic equilibrium, and therefore, a single and should be obtained by the weight average of different particle species. To see clearly the trends of weight average parameters, Figures 7(a) and 7(b) show the dependences of weight averages and on for different event centralities. The symbols represent the parameter values averaged by weighting the yields of different particles which are listed in Table 1. One can see that and increase with the increase in from 7.7 to 39 GeV. Meanwhile, and increase with the increase in the event centrality from the periphery to the center.

In addition, the variation of weight averages on for different collision energies and event centralities is displayed in Figure 7(c), where the symbols represent the parameter values averaged by weighting the yields of different particles. One can see that increases with the increase in . At higher energy and in central collisions, one sees larger and . There is a positive correlation between and .

The dependences of mean transverse momentum and initial temperature () on for different event centralities obtained by weighting the yields of different particles are shown in Figures 8(a) and 8(b), respectively, where according to the color string percolation model [5153]. One can see that and increase with the increase in from 7.7 to 39 GeV. Meanwhile, and increase with the increase in event centrality from the periphery to the center.

We notice that which is quite high for the considered collision energies. Because we obtain from the spectra of particles with nonzero masses, it is possible to have a high value. If we obtain from the spectra of photons, the value will be small. High renders that the excitation degree of the emission source at the stage of the initial state is high. Meanwhile, one of the key issues is whether the transverse flow should be also considered at the initial stage. Naturally, after considering the transverse flow at the initial stage, will be small. It is regretful that we have no clear idea on the extraction of transverse flow at the initial stage. As an alternative method, if we redefine as the initial temperature and as the initial transverse flow velocity, where is to be determined, we may obtain a small and a nonzero .

The reason for the increase in and with the increase in collision energy is due to the fact that more energies are deposited in collisions at higher energy in the considered RHIC-BES energy range. Meanwhile, the system size at higher energy decreases due to a relativistic constriction effect, which results in a smaller volume and then a larger energy density and larger . Meanwhile, at higher energy, the squeeze is more violent, which results in a rapider expansion and larger .

The reason for the increase in and with the increase in event centrality is due to the fact that the central collisions contain more nucleons than the peripheral collisions; then, more energies are deposited in central collisions. Meanwhile, a rapider expansion appears due to more violent squeeze in central collisions, compared to peripheral collisions. As a result, both and in central collisions are larger than those in peripheral collisions.

Because of and being positive correlation with and , with increasing collision energy and event centrality, the increasing trend of and results naturally in the increasing trend of and . This work shows that the two free parameters and and the two derived parameters and appear to be similar law on the dependences of collision energy and event centrality. In particular, and are model-independent, though we obtain them from model-dependent free parameters and in this work. In fact, and can be obtained by the data themselves if the data are across the possible range.

It seems that there are nonmonotonous changes at 11.5 GeV in the excitation functions of , , and in the most central Au-Au collisions. These nonmonotonous changes reflect the minimum or maximum point of equation of state (EoS) of the matter formed in collisions. At a few GeV to about 10 GeV, the matter formed in collisions is baryon-dominant. At above 10 GeV, the matter formed in collisions is meson-dominant. At around 10 GeV, the baryon number density is the largest [54] due to the competition between projectile/target penetrating/stopping and longitudinal contraction.

It is hard to say whether the minimum or maximum point of EoS of the matter formed in the most central Au-Au collisions at 11.5 GeV is related to the search for the QCD critical end point (CEP) which is the main objective of the BES program performed by the STAR Collaboration. Generally, large nonmonotonous changes or saturations or a slight increase should appear in the excitation functions of some quantities at the critical energy which is the energy corresponding to the CEP. The excitation functions considered in this paper change slightly. Although there is no value in the energy range of less than 7.7 GeV, it is expected that the excitation function increases quickly in the energy range of a few GeV while the onset stage of a slight increase appears at around 10 GeV in the excitation functions of and .

It should be noted that there is entanglement in the extraction of and . In fact, if one uses smaller and larger for central collisions, a decreasing trend for from peripheral to central collisions can be obtained. Meanwhile, a negative correlation between and can also be obtained. Thus, this situation is in agreement with some current references [19, 55, 56]. If one even uses almost invariant or slightly larger and properly larger for central collisions, an almost invariant or slightly increased trend for from peripheral to central collisions can be obtained [57]. To show the flexibility in the extraction of and , this work has reported an increasing trend for from peripheral to central collisions and a positive correlation between and .

In addition, we have taken in this work, which closely resembles the hydrodynamic profile as mentioned in Ref. [21]. Although Ref. [58] shows that is the closest approximation to hydrodynamics at freeze-out, Ref. [36] shows that or 1 does not affect obviously the fit curve and free parameters and . If we consider that decays quickly from the surface to the center of the emission source, we are inclined to use . Anyhow, we are not inclined to regard as a free parameter which is too mutable and debatable in our opinion. In current analysis with the blast-wave model [22], not only is mutable (from to ) but also the coverage is narrow and particle-dependent ( for , for , and for ), which uses a single kinetic freeze-out scenario and results in different trends of versus from this work. If we also regard as a free parameter and use narrow and particle-dependent coverage, consistent result with current analysis [22] can be naturally obtained by us.

Indeed, there are too much uncertainties arising from the choice of fit function and flow profile and from the well-known ambiguity in the fit results—in a single spectrum, it is always possible to trade against . That is, and is negatively correlative for a given spectrum. It is possible that we may use suitable and for a set of spectra and obtain a positive or negative correlation. In a positive correlation, decreasing and increasing will result in a negative correlation. Contrarily, in a negative correlation, increasing and decreasing will result in a positive correlation. Indeed, there is an influence if we use a changeable coverage and/or choice on the extraction of the two free parameters. In our opinion, to reduce the uncertainties, one should use a fixed flow profile () and wide and fixed coverage for different particles. In fact, we have used and for different particles in this work and used a multiple freeze-out scenario such as that used in Ref. [59].

Our result (Table 1) shows that the heavier the particle is, the higher and the smaller correspond. This result is in agreement with the hydrodynamic-type behavior [4]. The final and are averaged by weighting different particle yields, which shows a positive correlation between and (Figure 7). Our result is in agreement with the alternative method [36, 37] in which is regarded as the intercept in the linear relation of versus and is regarded as the slope in the linear relation of versus , where is the mean Lorentz factor in the source rest frame. Our result is also in agreement with a very recent work [60] which uses the same method as ours. If the negative correlation can be explained as the result of a longer lifetime (lower excitation degree) which corresponds to lower and a quicker expansion (stronger squeeze) which corresponds to larger , the positive correlation can be explained as the result of a high excitation degree which corresponds to high and quick expansion (strong squeeze) which corresponds to large .

This whole phenomenal analysis results in degrees of thermal motion and collective expansion that are reflected by and , respectively. With the increasing collision energy, the system may undergo different evolution processes. In the considered RHIC-BES energy range, the violent degree of collisions increases with increasing collision energy. The trends of and show approximately monotonous increase in which large fluctuation does not appear, though there are nonmonotonous changes at 11.5 GeV in some cases. The evolution processes at the considered six energies show similar behaviors to each other.

4. Conclusions

The main observations and conclusions are summarized as follows. (a)By using the method of data-driven reanalysis, the blast-wave model with the Boltzmann-Gibbs statistics is used to analyze the collision energy-dependent and event centrality-dependent double-differential transverse momentum spectra of charged particles (, , and ) produced in the midrapidity interval in Au-Au collisions at the RHIC-BES energies. The contribution of soft excitation is considered in this work, but the contribution of the hard process is not excluded if available(b)As the free parameters, the kinetic freeze-out temperature and transverse flow velocity are extracted with the blast-wave model. Both and increase with the increase in collision energy due to more violent collisions at higher energy. The two parameters also increase with the increase in event centrality, as the central collisions contain more nucleons which means more energy deposited and more violent collisions and squeeze, compared with peripheral collisions(c)As the derived parameters, the mean transverse momentum and initial temperature appear to be similar law to the free parameters and when we study the dependences of parameters on collision energy and event centrality. Although and are model-dependent, and are generally model-independent. There is no large fluctuation in the excitation function of the considered parameters at the RHIC-BES, which means a similar collision mechanism

Data Availability

The data used to support the findings of this study are included within the article and are cited at relevant places within the text as references.

Ethical Approval

The authors declare that they are in compliance with ethical standards regarding the content of this paper.

Disclosure

The funding agencies have no role in the design of the study; in the collection, analysis, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

We thank Dr. Muhammad Usman Ashraf for his kind help. This work was supported by the National Natural Science Foundation of China under Grant No. 11575103, the Chinese Government Scholarship (China Scholarship Council), the Scientific and Technological Innovation Programs (STIP) of Higher Education Institutions in Shanxi under Grant No. 201802017, the Shanxi Provincial Natural Science Foundation under Grant No. 201701D121005, and the Fund for Shanxi “1331 Project” Key Subjects Construction.