Abstract

The Tsallis distribution has been tested to fit all the particle spectra at mid-rapidity from central events produced in d + Au, Cu + Cu, and Au + Au collisions at RHIC and p + Pb, Pb + Pb collisions at LHC. Even though there are strong medium effects in Cu + Cu and Au + Au collisions, the results show that the Tsallis distribution can be used to fit most of particle spectra in the collisions studied except in Au + Au collisions where some deviations are seen for proton and at low . In addition, as the Tsallis distribution can only fit part of the particle spectra produced in Pb + Pb collisions where is up to 20 GeV/c, a new formula with one more fitting degree of freedom is proposed in order to reproduce the entire region.

1. Introduction

The heavy-ion collision experiments at RHIC and LHC give us the opportunity to study the phase transition from nuclear matter to quark gluon plasma (QGP), the collective motion, the nuclear medium effects, and so on. The particle spectrum is one of the basic quantities measured in experiments to address the questions raised in such studies. Recently, the Tsallis distribution has attracted the attention of many theorists and experimentalists in high energy heavy-ion collisions [125]. It has been applied to particle spectra produced in different reaction systems, from pp, pA to AA, to understand the particle production mechanism and extract physical quantities, for example, temperature [24, 7, 1214, 1721, 24] and chemical potential [26]. In pp collisions, the excellent ability to fit the spectra of identified hadrons and charged particles in a large range of up to 200 GeV/c is quite impressive [2, 2325]. A systematic investigation of particle spectra in p + p collisions at RHIC and LHC has been conducted in [2]. The results show that the Tsallis distribution can fit all the particle spectra at different energies in p + p collisions. A possible cascade particle production mechanism is proposed. Recently, a Tsallis distribution scaling function was found for charged hadron spectra in p + p and collisions [1]. Comparing to nucleus-nucleus collisions, the pp collision is very simple. It has been used as a baseline for nucleus-nucleus collisions. A nuclear modification factor or was proposed to show the nuclear medium effects in pA or AA collisions referring to pp collisions [5, 1012, 2734]. The nuclear modification factor different from unity is a manifestation of medium effects. Many authors have successfully applied Tsallis distribution to fit particle spectra in pA and AA collisions even though the spectra were affected by nuclear medium modification [38, 10, 11, 18, 21]. We also notice that many works only show the small part of the particle spectra, while the exponential distribution also can fit the low region [11, 35]. It should cover all regions of particle spectra, available in experiment, in order to show the advantage and/or the fitting power of the Tsallis distribution. In recent years, the experimental groups at RHIC and LHC have published the wide range of particle spectra for different particles in different reaction systems. Such data allow us to conduct the systematic study of particle spectra in heavy-ion collisions at RHIC and LHC using the Tsallis distribution, as we have done for p + p collisions [2].

In this work, we would like to test whether the Tsallis distribution can fit all the particle spectra produced at RHIC and LHC, which can help us to understand the particle production mechanism. Before we start to conduct our investigation, we can get some clues to estimate whether it can fit the particle spectrum or not from the nuclear modification factor. If or is flat for the whole region, according to its definition, this means that the particle spectrum is similar in shape and only differs in magnitude to the one in p + p collisions. Based on the previous studies [1, 2], we are sure that the Tsallis distribution can fit the particle spectrum since it can fit all the particle spectra produced in p + p collisions, especially up to extremely high [2, 2325]. In the pA reactions, are flat and very close to 1 for most of the produced particles at different centralities [5, 10, 11, 27, 29, 34], while in the AA collisions the nuclear medium effects play an important role. increases from the most central collisions to peripheral collisions [11, 12, 28, 3033]. Since the nuclear modification factors of different particles are almost 1 in peripheral heavy-ion collisions [5, 1012, 2734], as we discussed, Tsallis distribution should be able to fit the particle spectra. Therefore we will only focus on the particle spectra at the most central collisions where the Tsallis distribution may not fit all of them. We have collected data of particle spectra from d + Au, Cu + Cu, and Au + Au collisions at RHIC and p + Pb, Pb + Pb collisions at LHC and select the data for most of the particles with the highest  GeV to conduct this study.

The paper is organized as follows. In Section 2, we show different versions of the Tsallis distribution used in the literature. More details can be found in [2]. We also give the form of Tsallis distribution used in our analysis. In Section 3, we show our results of particle spectra from d + Au, p + Pb, Cu + Cu, Au + Au, and Pb + Pb. Another distribution is proposed to fit the particle spectra in Pb + Pb collisions since Tsallis distribution can only fit part of the particle spectra in the case. A brief conclusion is given in Section 4.

2. Tsallis Distributions

In the literature, several versions of Tsallis distribution with different arguments can be found [225]. The asymptotic behaviors of these distributions at low and high limits can be found in [2]. We only briefly show them here.

The STAR [9] and PHENIX [5, 12] Collaborations at RHIC along with ALICE [13, 14] and CMS [15] Collaborations at LHC adopted this form of Tsallis distribution:where is the transverse mass and is the mass of the particle. , , and are fitting parameters.

In [8, 1922, 36], the following Tsallis form is used:based on thermodynamic consistency arguments. Where is the degeneracy of the particle, is the volume, is the rapidity, is the chemical potential, is the temperature, and is a parameter. In (2), there are four parameters , , , and . was assumed to be  0 in [1921, 36] which is a reasonable assumption because the energy is high enough and the chemical potential is small compared to temperature. In the mid-rapidity region, (2) is reduced toIn [4, 21], (2) has been rewritten asto take into account the width of the corresponding rapidity distribution of the particles.

In [17], Sena and Deppman applied the nonextensive formalism to obtain the probability of particle with momentum aswhere is the normalization constant, is a parameter, , and is the mass of particle. With the approximation very large compared to and [37], (5) can be rewritten aswhere and is the beta-function.

In [24], Wong and Wilk proposed a new form of the Tsallis distribution function to take into account the rapidity cut: wherewith where is assumed to be a constant.

In [2], we have obtainedwhere , , and are the fitting parameters. This is equivalent to (1) but in a simpler form. We adopt (10) to do the analysis here. We notice that (10) has been used by the CMS Collaboration [16, 38] and by Wong et al. in their recent paper [25]. The STAR Collaboration has also applied a formula which is very close to (10) [29].

3. Results

We have selected the data of particle spectra from the most central collisions with the highest  GeV for most of the particles in d + Au, p + Pb, Cu + Cu, A u + Au, and Pb + Pb at RHIC and LHC. We fit the center values of the experimental points. The fit metric used is defined by

As we discussed before, the Tsallis distribution should be able to fit the particle spectra from d + Au and p + Pb. One good example has been shown in [4, 5] for and in d + Au at GeV where their spectra can be obtained by multiplying the particle spectra in p + p collisions with . In Figures 1 and 2, our results for d + Au at GeV and p + Pb at  TeV using (10) have been shown. In order to see the agreement between the data and the Tsallis distribution in linear scale, a ratio data/fit is defined. As shown in Figures 1 and 2, the fits for all particles are good. For the left collision systems, we also do the same comparisons. We would like to emphasize that of charged particle spectrum is up to 45 GeV/c.

Now let us turn to the AA collisions. First we use Tsallis distribution Equation (10) to fit the particle spectra in Au + Au collisions at GeV in Figure 3. All the particle spectra are well fitted except the proton spectrum at GeV/c. This makes a little difference of AA collisions from p + p collisions. We want to check whether this deviation will become larger at higher colliding energy in AA collisions. We considered the particle spectra from Cu + Cu collisions at GeV. The results are shown in Figure 4. The fitting with (10) for different particle spectra is very well. But we do not know whether there is deviation or not for proton at low since the data for GeV/c are not available. Fortunately, the data for different particle spectra at low in Au + Au collisions at GeV are given. In Figure 5, one can see the deviations of particle spectra of proton and at low from the Tsallis distribution Equation (10). While a deviation is observed for proton at  GeV/c which becomes a little larger than the one in Au + Au at GeV, all other particle spectra are well fitted. This makes us curious to fit the particle spectra in Pb + Pb collisions at  TeV. With the successful running at LHC, the identified hadron particle spectra data in Pb + Pb collisions at  TeV are available up to 20 GeV/c. The data satisfy two criteria. One is that there are strong nuclear medium effects in Pb + Pb collisions which can be seen from and the other is that the transverse momenta of the particles reach high values. This gives us an opportunity to test the fitting power of the Tsallis distribution. When we use Tsallis distribution Equation (10) to fit pion spectrum, we find that (10) can only fit part of it. If we choose to fit the low region, (10) can fit the particle spectrum up to 10 GeV/c, as shown in Figure 6 with red dashed line. The blue dotted-dashed line shows the fit for high region which starts from 4 GeV/c.

The exponential form equation was used to fit the particle spectra at RHIC when only the low data are available [11, 35]. With the upgrade of detectors, we have a better ability to measure the particle spectra. When the intermediate data are available, the Tsallis distribution is used to understand the particle spectra and extract physical information. The two-Boltzmann distribution was also used in [7]. But in both cases the number of free fitting parameters increases from 2 to 3. One fitting degree of freedom is increased in this transition. In [6], a double Tsallis formula was proposed to fit particle spectra obtained from central events in Pb + Pb collisions. In this case, three fitting degrees of freedom are increased. Here we will follow the same logic of the transition from exponential distribution to Tsallis distribution to propose a new form equation to fit the particle spectra in Pb + Pb collisions at  TeV by only increasing one fitting degree of freedom. We increase the number of free fitting parameters from 3 to 4 and the proposed formula isThere are four parameters , , , and . We are inspired by the solution of Fokker-Planck equation [56]. We change the power from  2 in [56] to 4 in (12) in order to fit well all the particle spectra with one equation. Figure 6 shows that the fits with (12) are excellent. We would like to mention that when , (12) becomes and when ,Equation (12) has the same asymptotic behaviors as (10).

4. Conclusions

In this paper, we have tested the fitting ability of Tsallis function by fitting different particle spectra produced at the most central collisions in d + Au, p + Pb, Cu + Cu, Au + Au, and Pb + Pb at RHIC and LHC. The Tsallis distribution is able to fit all the particle spectra in d + Au and p + Pb collisions where the medium effects are very weak. This information can be obtained by the nuclear modification factor. In the AA collisions, the Tsallis distribution can fit all the particle spectra very well at RHIC energies except the little deviation observed for proton and at low . However the Tsallis distribution can only fit part of the particle spectra in Pb + Pb at  TeV, either in the low or in the high region. We have proposed a new formula in order to fit all the particle spectra in Pb + Pb by increasing one fitting degree of freedom from Tsallis distribution. This follows the same idea of the transition from the exponential distribution to Tsallis distribution when intermediate data are available in experiments.

According to the results in this paper and [2], we conclude that we can do the systematic analysis of particle spectra with Tsallis distribution in p + p, pA at RHIC and LHC. In the AA collisions at GeV/c, we can do the same analysis as in p + p and pA at RHIC and LHC. But when we consider Pb + Pb collisions, the Tsallis distribution fails.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Dr. J. Mabiala for reading carefully their paper. This work was supported, in part, by the NSFC of China under Grant no. 11205106.