Advances in High Energy Physics

Volume 2013 (2013), Article ID 965735, 15 pages

http://dx.doi.org/10.1155/2013/965735

Research Article

Transverse Momentum Distributions of Final-State Particles Produced in Soft Excitation Process in High Energy Collisions

Fu-Hu Liu, Ya-Hui Chen, Hua-Rong Wei, and Bao-Chun Li

Institute of Theoretical Physics, Shanxi University, Taiyuan, Shanxi 030006, China

Received 17 July 2013; Revised 5 September 2013; Accepted 20 September 2013

Academic Editor: Gonang Xie

Copyright © 2013 Fu-Hu Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Transverse momentum distributions of final-state particles produced in soft process in proton-proton (pp) and nucleus-nucleus (AA) collisions at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) energies are studied by using a multisource thermal model. Each source in the model is treated as a relativistic and quantum ideal gas. Because the quantum effect can be neglected in investigation on the transverse momentum distribution in high energy collisions, we consider only the relativistic effect. The concerned distribution is finally described by the Boltzmann or two-component Boltzmann distribution. Our modeling results are in agreement with available experimental data.

1. Introduction

High energy collisions are an important research topic in particle and nuclear physics. The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) did firstly collider experiments on heavy ions [1], and the center-of-mass energy per nucleon pair at the RHIC reached highly 200 GeV [2]. The Large Hadron Collider (LHC) at European Laboratory for Particle Physics (CERN) renovated value of to TeV region [3]. It seems that a new state of matter, namely, Quark-Gluon Plasma (QGP), is possibly formed in heavy ion collisions at RHIC and LHC energies due to high temperature and density [4, 5]. At initial stage of high energy collisions, another possible new state of matter, namely, color glass condensate (CGC), is caused by strong color fields in the low-gluon realm [6, 7], wheredenotes the ratio of quark or gluon momentum to hadron one. A CGC is in fact a region of the nuclear wave function at low-andand exists already before the collisions, wheredenotes the square momentum of virtual photon. On the other hand, the CGC may not be a new state, but more like a model or calculation for initial state hadron behavior. One cannot measure the QGP and CGC directly. However, one can measure final-state particle spectra at freeze-out to extract thermal and other characteristics of interacting system and give a judgment on formation and property of the new matters.

The final-state particle spectra include rapidity (or pseudorapidity ) [8, 9], transverse momentum (or transverse mass) [10, 11], transverse energy [12, 13], and other distributions [14]. It is known that and distributions reflect, respectively, the degrees of longitudinal extension and transverse excitation of interacting system. Especially, for transverse excitation, soft excitation and hard scattering processes can affect, respectively, distributions in low- and high- ranges. The soft and hard processes correspond to different physics mechanisms and distribution laws [15]. In low energy collisions, the soft process is main process, and the hard process can be neglected due to almost zero contribution. In high energy collisions, although the hard process cannot be neglected, the soft process is still main process.

To understand the transverse excitation, we need firstly to study the soft excitation process. A lot of models have been introduced to describe the soft process, although some of them can be used to describe the hard process too [16, 17]. Among the models, the multisource thermal model proposed by us is a very simple one and can be used to describe spectra in both the soft and hard processes if source’s contribution is given by an Erlang distribution [18]. Finally, the considered distribution is described by a multicomponent Erlang distribution [19, 20]. Different from some simulation codes, our model gives directly a few statistical laws by analytical expressions in describing some quantities. In the case of being incapable of analytical expressions, we could use a Monte Carlo method to give a numerical result. Our model is easy to be used by experimental experts.

Due to significances of the considered model and topic, in this paper, based on Boltzmann distribution for a single source, we describe spectra of final-state particles produced in soft process in proton-proton () and nucleus-nucleus () collisions at RHIC and LHC energies. Some interesting results are obtained.

2. The Model

According to the multisource thermal model, many emission sub-sources of final-state particles are assumed to form in high energy collisions [19, 20]. These multiple sub-sources can be different regions in the overlap region or different mechanisms, and these particles can be created/emitted at different times in the collisions. In fact, these sub-sources can be divided into different groups (sources) due to different interacting mechanisms or event samples. Obviously, soft process corresponds to sources with low degree of excitation or to particles with low transverse momentum, and hard process corresponds to sources with high degree of excitation or to particles with high transverse momentum, where the excitation means to create particles through string breaking, direct scattering, recombination, and their hybrid.

In the rest frame of a source, we consider the source as a thermodynamic system of relativistic and quantum ideal gas. The momentum () distribution of final-state particles in the natural unit system is given by [21]

where is the number of particles, is the normalization constant, is the rest mass of a considered particle, is the chemical potential, is the temperature parameter, denotes fermions, and denotes bosons, respectively. Our calculations show that at RHIC and LHC energies the quantum effect and chemical potential can be neglected compared to the relativistic effect [22]. Then, we have a simple expression for momentum distribution to be [23, 24]

where and is the modified Bessel function of order 2.

The distribution can be written as a Boltzmann distribution [25]:

where is the normalization constant. Because of interactions among different sources, the considered source has a deformation and/or movement in the transverse plane. Letdenote a relative deformation and letdenote an absolute movement of the source; that is, we use instead of in (3). The revised distribution can be given by

In the case of considering multiple sources, we have

or

where , , and denote the contribution ratio, normalization constant, and temperature of theth source, respectively. Because the effects of deformation and movement of the source can be neglected in the calculation of transverse momentum, we take the default values of and in the revised distribution, which results from the Boltzmann distribution or a multicomponent Boltzmann distribution.

We should have a few sources to describe the soft and hard processes. For the soft process, the number of sources is generally 1 or 2. For the hard process, the number of sources is also 1 or 2. The total number of sources will be from 2 to 4 for a wide distribution. In this paper, we pay our attention on the soft process which has a narrow distribution. It is hard to say that what the distribution range is for the soft process. What we can say is that for low energy collisions the distribution range is narrower. In the present work, we regard the distribution range as 0–10 GeV/c. The difference between the single and multisource models is obvious. The former one describes a narrower distribution which corresponds to an equilibrium state with a lower degree of excitation. The latter one describes a wider distribution which corresponds to a few local equilibrium sates with different excitations.

3. Comparisons with Experimental Data

Figure 1 presents the transverse momentum distributions of (a) , and , (b) , and , (c) , (d) , and (e) produced in collision at center-of-mass energy  GeV with different ranges and magnifications shown in the figure. The symbols represent the experimental data of the STAR [26] (Figures 1(a) and 1(b)) and BRAHMS Collaborations [27] (Figures 1(c)–1(e)), and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution. In the calculation, we have used a fitting method to obtain parameter values which are shown in Table 1 with values of per degree of freedom (). To give a short presentation, the values corresponding to “negative/positive” charged particles are given in terms of “the first value/the second value” or “value” in the case of the first value and the second value being the same. We would like to point out that the presenting styles of rapidity ranges for Figures 1(a) and 1(b) as well as for Figures 1(c)–1(e) are different due to different presentations in [26, 27]. One can see that the modeling results with 1 or 2 sources are in agreement with the experimental data. For emissions of , and , the temperature parameter increases with increase of particle mass (Figures 1(a) and 1(b)), which indicates the impact of radial flow and/or the early emission of heavy hadrons. The temperature does depend nonobviously on rapidity range (Figures 1(c)–1(e)).

Table 1: Parameter values corresponding to the curves in Figure 1.

Figure
Collision
Particle
Rapidity
(GeV)

(GeV)

Figure 1(a)/Figure 1(b)
pp 200 GeV

0.165
1.000
—
1.128/1.212

0.187
1.000
—
0.211/0.243

/

0.197/0.212
1.000
—
0.340/0.274

Figure 1(c)
pp 200 GeV

0.175
0.896
0.360
0.849

0.175
0.930
0.360
0.958

0.175
0.912
0.360
0.968

0.175
0.823
0.270
1.415

0.175
0.856
0.270
1.170

0.166
0.888
0.270
1.669

0.205
1.000
—
1.273

Figure 1(d)
pp 200 GeV

0.160
0.700
0.360
0.499

0.160
0.600
0.360
0.786

0.160
0.600
0.360
0.437

0.198
0.800
0.360
1.048

0.176
0.788
0.300
0.758

0.160
0.780
0.273
1.689

0.160
0.667
0.253
0.714

Figure 1(e)
pp 200 GeV

0.160
0.400
0.265
1.028

0.160
0.540
0.290
1.240

0.160
0.540
0.290
0.821

0.156
0.875
0.280
1.457

0.156
0.872
0.280
0.750

0.190
0.859
0.280
1.561

0.176
1.000
—
0.907

(a)

(b)

(c)

(d)

(e)

Figure 1: The transverse momentum distributions of (a) , , and , (b) , , and , (c) , (d) , and (e) produced in collision at  GeV with differentranges and magnifications are shown in the figure. The symbols represent the experimental data of the STAR [26] and BRAHMS Collaborations [27] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution.

Figure 2 shows the distributions of produced in (a) and (b) collision at  GeV, (c) d-Au collisions at  GeV, and (d) p-Pb collisions at beam energy being 400 GeV with different y ranges, global scale uncertainty (GSU) (or invariant mass ()) ranges, and magnifications shown in the figure, where anddenote the cross section and dilepton branching ratio, respectively. The symbols represent the experimental data of the PHENIX [28] (Figures 2(a) and 2(c)) [29, 30] (Figure 2(b)) and NA50 Collaborations [31] (Figure 2(d)), and the curves are our results calculated by the Boltzmann or two-component Boltzmann distributions. The values of parameters and are given in Table 2. We see again that the model with 1 or 2 sources describes the experimental data. For emission of the temperature parameter does depend nonobviously on rapidity range (Figures 2(a) and 2(c)).

Table 2: Parameter values corresponding to the curves in Figure 2.

Figure
Collision
Particle
Rapidity
(GeV)

(GeV)

Figure 2(a)
pp 200 GeV
in pp

0.462
1.000
—
0.663

0.462
1.000
—
1.157

0.462
1.000
—
1.383

Figure 2(b)
pp 200 GeV
in pp

0.160
0.175
0.541
1.371

0.160
0.175
0.654
1.638

Figure 2(c)
d-Au 200 GeV
in d-Au

0.525
1.000
—
1.332

0.525
1.000
—
1.125

0.525
1.000
—
1.845

Figure 2(d)
400 GeV p-Pb
in p-Pb

0.205
0.650
0.368
1.041

(a)

(b)

(c)

(d)

Figure 2: The transverse momentum distributions of produced in (a) and (b) collision at  GeV, (c) d-Au collisions at  GeV, and (d) p-Pb collisions at beam energy being 400 GeV with differentranges, GSU (or ) ranges, and magnifications are shown in the figure. The symbols represent the experimental data of the PHENIX [28–30] and NA50 Collaborations [31] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution.

The distributions of identified charged particles produced in-Au collisions at  GeV, Cu-Cu collisions at  GeV, Au-Au collisions at , 130, and 200 GeV with different centrality classes are displayed in Figures 3–7, respectively. The symbols represent the experimental data of the STAR [26] (Figures 3, 5, 6, and 7) and PHENIX Collaboration [32] (Figure 4), and the curves are our results calculated by the Boltzmann or two-component Boltzmann distributions. Correspondingly, the values of parameters and are given in Tables 3–7, respectively. Once more, the model with 1 or 2 sources describes the experimental data. From Tables 5, 6, and 7 we see clearly that for emissions of , , and (Figures 5(a), 6(a), and 7(a)), as well as , , and (Figures 5(b), 6(b), and 7(b)), the temperature parameter increases with increases of particle mass, impact centrality, and , where we would like to point out that a large centrality (a small percentage) corresponds to a small impact parameter. The similar conclusions can be obtained from Tables 3 and 4.

Table 3: Parameter values corresponding to the curves in Figure 3(a)/Figure 3(b). The collisions are d + Au at 200 GeV.

Particle
Centrality
(GeV)

0–20%
0.178
0.564/0.704

20–40%
0.170/0.175
1.053/0.779

40–100%
0.170
0.879/0.905

0–100%
0.173
0.796/0.771

0–20%
0.206/0.228
0.119/0.052

20–40%
0.206/0.218
0.112/0.114

40–100%
0.206/0.213
0.216/0.074

0–100%
0.206/0.217
0.137/0.131

0–20%
0.243/0.239
0.036/0.032

20–40%
0.232/0.229
0.047/0.036

40–100%
0.215/0.212
0.091/0.082

0–100%
0.229
0.032/0.013

Table 4: Parameter values corresponding to the curves in Figure 4.

Figure
Collision
Particle
Centrality
(GeV)

(GeV)

Figure 4(a)/Figure 4(b)
Cu-Cu

22.5 GeV

0–10%
0.183
0.934/0.912
0.310
1.079/1.008

10–30%
0.183
0.934/0.917
0.310
1.195/0.848

30–60%
0.173/0.183
0.934/0.931
0.300/0.310
0.992/1.810

60–100%
0.162/0.160
0.934/0.900
0.293/0.270
1.394/0.959

0–100%
0.183
0.934/0.920
0.310
1.229/1.084

Figure 4(c)/Figure 4(d)
Cu-Cu

22.5 GeV

0–10%
0.183
0.860/0.747
0.310
0.636/0.593

10–30%
0.183
0.840/0.732
0.310
0.632/0.818

30–60%
0.183
0.900/0.755
0.310
1.087/1.440

60–100%
0.183
0.975/0.747
0.310
1.677/1.311

0–100%
0.183
0.860/0.747
0.310
0.632/0.515

Figure 4(e)/Figure 4(f)
Cu-Cu

22.5 GeV

0–10%
0.183/0.200
0.760/0.510
0.310
1.898/0.904

10–30%
0.183/0.200
0.860/0.567
0.310
0.825/1.159

30–60%
0.183/0.200
0.883/0.580
0.310
1.849/1.747

60–100%
0.160/0.200
1.000/0.888
—/0.310
1.953/1.855

0–100%
0.183/0.200
0.862/0.566
0.310
1.090/1.466

Table 5: Parameter values corresponding to the curves in Figure 5(a)/Figure 5(b). The collisions are Au + Au at 62.4 GeV.

Particle
Centrality
(GeV)

0–5%
0.185
0.184/0.184

5–10%
0.185
0.192/0.156

10–20%
0.182
0.224/0.246

20–30%
0.180
0.230/0.228

30–40%
0.180
0.276/0.302

40–50%
0.175
0.112/0.128

50–60%
0.170
0.118/0.146

60–70%
0.170
0.138/0.234

70–80%
0.165
0.254/0.292

0–5%
0.275
0.056/0.010

5–10%
0.275
0.024/0.018

10–20%
0.275
0.040/0.018

20–30%
0.253
0.064/0.046

30–40%
0.243
0.102/0.154

40–50%
0.235
0.082/0.132

50–60%
0.215
0.166/0.084

60–70%
0.205
0.100/0.080

70–80%
0.205
0.362/0.130

0–5%
0.512
0.114/0.166

5–10%
0.475
0.068/0.120

10–20%
0.442
0.048/0.096

20–30%
0.393
0.040/0.136

30–40%
0.355
0.076/0.082

40–50%
0.305
0.104/0.080

50–60%
0.270
0.070/0.086

60–70%
0.250
0.038/0.058

70–80%
0.220
0.020/0.084

Table 6: Parameter values corresponding to the curves in Figure 6(a)/Figure 6(b). The collisions are Au + Au at 130 GeV.

Particle
Centrality
(GeV)

0–6%
0.185
0.210/0.268

6–11%
0.182
0.290/0.232

11–18%
0.180
0.430/0.328

18–26%
0.180
0.394/0.310

26–34%
0.175
0.222/0.172

34–45%
0.170
0.126/0.144

45–58%
0.170
0.192/0.234

58–85%
0.165
0.334/0.354

/
0–6%
0.275
0.088/0.146

6–11%
0.275
0.268/0.087

11–18%
0.273/0.263
0.151/0.584

18–26%
0.283
0.189/0.584

26–34%
0.255/0.235
0.101/0.598

34–45%
0.260/0.245
0.846/0.293

45–58%
0.245/0.205
0.639/0.448

58–85%
0.205
0.311/0.543

/
0–6%
0.580
0.177/0.177

6–11%
0.530
0.157/0.128

11–18%
0.450
0.123/0.104

18–26%
0.410
0.105/0.118

26–34%
0.410
0.103/0.068

34–45%
0.363
0.071/0.121

45–58%
0.303
0.076/0.065

58–85%
0.265
0.029/0.055

Table 7: Parameter values corresponding to the curves in Figure 7(a)/Figure 7(b). The collisions are Au + Au at 200GeV.

Particle
Centrality
(GeV)

0–5%
0.193
0.940/0.746

5–10%
0.193
0.865/0.907

10–20%
0.193/0.190
0.780/0.444

20–30%
0.190/0.188
0.628/0.661

30–40%
0.188
1.432/1.449

40–50%
0.188
0.570/0.577

50–60%
0.180
0.974/0.914

60–70%
0.175
1.269/1.295

70–80%
0.170
1.269/1.311

0–5%
0.327
0.019/0.020

5–10%
0.327
0.019/0.032

10–20%
0.285
0.044/0.034

20–30%
0.285
0.029/0.022

30–40%
0.243
0.112/0.107

40–50%
0.243
0.082/0.079

50–60%
0.243
0.007/0.021

60–70%
0.215
0.081/0.071

70–80%
0.200
0.235/0.143

0–5%
0.570
0.046/0.031

5–10%
0.516
0.071/0.013

10–20%
0.516
0.037/0.045

20–30%
0.435
0.047/0.035

30–40%
0.396
0.039/0.036

40–50%
0.338/0.345
0.035/0.032

50–60%
0.312/0.320
0.025/0.031

60–70%
0.270
0.086/0.057

70-80%
0.246
0.045/0.022

(a)

(b)

Figure 3: The transverse momentum distributions of (a) , , and as well as (b) , , and produced in d-Au collisions at  GeV with and different centrality classes and magnifications are shown in the figure. The symbols represent the experimental data of the STAR Collaboration [26] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4: The transverse momentum distributions of (a) , (b) , (c) , (d) , (e) , and (f) produced in Cu-Cu collisions at  GeV with different centrality classes and magnifications are shown in the figure. The symbols represent the experimental data of the PHENIX Collaboration [32] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution.

(a)

(b)

Figure 5: The transverse momentum distributions of (a) , , and , and (b) , , and produced in Au-Au collisions at  GeV with different centrality classes and magnifications are shown in the figure. The symbols represent the experimental data of the STAR Collaboration [26] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution.

(a)

(b)

Figure 6: The same as Figure 5, but showing the results at  GeV.

(a)

(b)

Figure 7: The same as for Figure 5, but showing the results at  GeV.

Figure 8 gives the distributions of (a)–(c) identified particles and (d) charged particles in range of in nonsingle diffraction (NSD) produced in collision at  GeV, where and denote numbers of events and charged particles, respectively. The symbols represent the experimental data of the ALICE Collaboration [33, 34] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution. The values of parameters and are given in Table 8. We see that the model with 1 or 2 sources describes the experimental data. For emissions of charged hadrons, the temperature parameter increases with increase of particle mass.

Table 8: Parameter values corresponding to the curves in Figure 8.

Figure
Collision
Particle
(GeV)

(GeV)

Figure 8(a)/Figure 8(b)
pp 900 GeV

0.148/0.142
0.803/0.765
0.349/0.341
0.153/0.369

/
0.185
0.670
0.427
0.158/0.375

/
0.185
0.700/0.670
0.390/0.427
0.103/0.156

Figure 8(c)
pp 900 GeV

0.132
0.505
0.369
0.375

0.323
1.000
—
1.026

0.311
1.000
—
1.375

0.362
1.000
—
0.890

0.325
1.000
—
0.664

Figure 8(d)
pp 900 GeV
charged
0.259
0.985
0.754
1.903

(a)

(b)

(c)

(d)

Figure 8: The transverse momentum distributions of (a) , , and , (b) , , and , (c) ,, , , and , and (d) charged particles in range of in NSD produced in pp collisions at  GeV with different magnifications are shown in the figure. The symbols represent the experimental data of the ALICE Collaboration [33, 34] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution.

The distributions of identified particles produced in (a) central Pb-Pb collisions at  TeV, (b) Pb-Pb collisions with different centralities at  TeV and inelastic collision at  TeV, (c) central rapidity region in Pb-Pb collisions with different centralities at  TeV, and (d) central rapidity region in collision at  TeV are presented in Figure 9. The symbols represent the experimental data of the ALICE Collaboration [35–38] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution. The values of parameters and are given in Table 9. We see that in most cases the model with 1 or 2 sources describes the experimental data. Especially, for emissions of , , and (Figure 9(a)), the temperature parameter increases with increase of particle mass; for emission of (Figure 9(b)), the temperature parameter does not depend on impact centrality; and for emission of (Figure 9(c)), the temperature parameter increases with increase of impact centrality (or with decrease of impact parameter).

Table 9: Parameter values corresponding to the curves in Figure 9.

Figure
Collision
Particle
Type
(GeV)

(GeV)

Figure 9(a)
Pb-Pb 2.76 TeV

0–5%
0.169
0.595
0.364
0.858

0–5%
0.355
0.900
0.460
0.141

0–5%
0.697
1.000
—
0.832

Figure 9(b)
Pb-Pb 2.76 TeV

0–20%
0.332
0.985
0.900
1.321

20–40%
0.332
0.985
0.900
1.760

40–60%
0.332
0.985
0.900
1.324

60–80%
0.332
0.975
0.900
0.583

pp 2.76 TeV

Inelastic
0.332
0.965
0.900
1.816

Figure 9(c)
Pb-Pb 2.76 TeV

0–5%,
0.632
1.000
—
0.517

5–10%,
0.600
1.000
—
0.324

10–20%,
0.580
1.000
—
0.275

20–30%,
0.540
1.000
—
0.200

30–40%,
0.532
1.000
—
0.275

40–50%,
0.461
1.000
—
0.432

50–60%,
0.450
1.000
—
0.296

60–70%,
0.385
1.000
—
0.371

Figure 9(d)
pp 7 TeV

0.350
0.700
0.730
0.510

0.565
1.000
—
0.666

(a)

(b)

(c)

(d)

Figure 9: The transverse momentum distributions of (a) , , and , (b) , (c) , and (d) and produced in Pb-Pb collisions at  TeV ((a), (b), and (c)) and pp collision at and 7 TeV ((b) and (d)) with differentranges, centrality classes, and magnifications are shown in the figure. The symbols represent the experimental data of the ALICE Collaboration [35–38] and the curves are our results calculated by the Boltzmann or two-component Boltzmann distribution.

From Tables 1–9 we see that some values of are too low pointing to overestimated errors of the experimental points. In fact, in the case of errors being not available in related references, we have used a half size of the experimental points to give the errors. This treatment may cause larger errors in some cases.

4. Discussions and Conclusions

From the above discussions we see that the model used in the present work is just a simple phenomenology which does not contain other processes such as parton-hadron string dynamics, hydrodynamic flows, and resonances. The successful description renders that the mentioned processes should contribute a higher transverse momentum at multi-GeV energy or a refined structure in distribution curve. In the concerned transverse momentum region and for the concerned distribution curves, we just need to consider the Maxwell-Boltzmann thermal law.

The present work is justified to compare fits in a low transverse momentum region (<10 GeV/c) for different particles by the same thermal law. Although the difference for charged and neutral particles is unlikely due to Coulomb effects which are important for very soft charged particles only, both the charged and neutral particles obey the same thermal law. The transverse momentum can extend to more than 100 GeV/c at multi-GeV energy. The distribution in the low transverse momentum region is mainly contributed by the soft processes. The hard processes which contribute high transverse momentums can be partly described by the thermal law.

To conclude, we have used the multisource thermal model to describe the transverse momentum distributions of particles produced in the soft process in and collisions at RHIC and LHC energies. For single source, the relativistic ideal gas model is applied in description of particle behavior. The concerned distribution is finally described by single source or two sources which result from a Boltzmann or two-component Boltzmann distributions. The modeling results are in agreement with available experimental data, which renders that an equilibrium or two local equilibriums are reached in high energy collisions. Because of the evolvement time of interesting system in collisions being very short, the particles should reach rapidly to the state of equilibrium.

The present work can be used to extract nuclear temperature for soft process. For emissions of charged hadrons, the temperature parameter increases with increases of particle mass, impact centrality, and center-of-mass energy and does depend nonobviously on rapidity range. That the temperature increases with particle mass indicates the impact of radial flow and/or the early emission of heavy hadrons. The temperature parameter for emission of does depend nonobviously on rapidity range too, which is consistent with charged hadrons. However, for emission of the temperature parameter does not depend on impact centrality, which is inconsistent with charged hadrons. Different behaviors for and render different production mechanisms. Especially, there are Coulomb corrections for emissions of charged particles, which affects the extraction of temperature [39].

The values of temperature parameter for emissions of are about 160–190 MeV which reaches the temperature (166–172 MeV) of creating QGP at zero baryon-chemical potential, where 172 MeV is the equilibrium phase transition temperature and 166 MeV is due to finite hadron size [40]. In most cases the temperature for emission of heavy hadrons is greater than that for pions, which renders the impact of radial flow and/or the early emission of heavy hadrons in collisions.

Conflict of Interests

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant no. 10975095 and no. 11247250, the China National Fundamental Fund of Personnel Training under Grant no. J1103210, the Open Research Subject of the Chinese Academy of Sciences Large-Scale Scientific Facility under Grant no. 2060205, and the Shanxi Scholarship Council of China.

References

B. B. Back, M. D. Baker, D. S. Barton et al., “Charged-particle pseudorapidity density distributions from Au + Au collisions at  GeV,” Physical Review Letters, vol. 87, no. 10, Article ID 102303, 2001.

I. G. Bearden, D. Beavis, C. Besliu, et al., “Pseudorapidity distributions of charged particles from Au+Au collisions at the maximum RHIC energy, GeV,” Physical Review Letters, vol. 88, no. 20, Article ID 202301, 2001.

A. Toia and ALICE Collaboration, “Bulk properties of Pb-Pb collisions at TeV measured by ALICE,” Journal of Physics G, vol. 38, no. 12, Article ID 124007, 2011.

H. Song, S. A. Bass, U. Heinz, T. Hirano, and C. Shen, “200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid,” Physical Review Letters, vol. 106, no. 19, Article ID 192301, 2011.

G. Roland, “QCD studies with CMS at LHC,” in INT Program 10-2a Quantifying the Properties of Hot QCD Matter, Seattle, Wash, USA, May 2010.

F. Gelis, “Initial stages of heavy-ion collisions from the CGC,” in INT Program 10-2a Quantifying the Properties of Hot QCD Matter, Seattle, Wash, USA, May 2010.

A. Dumitru, “CGC in heavy-ion coll,” in INT Program 10-2a Quantifying the Properties of Hot QCD Matter, Seattle, Wash, USA, May 2010.

B. Alver, B. B. Back, M. D. Baker, et al., “Charged-particle multiplicity and pseudorapidity distributions measured with the PHOBOS detector in Au+Au, Cu+Cu, d+Au, and p+p collisions at ultrarelativistic energies,” Physical Review C, vol. 83, no. 2, Article ID 024913, 2011.

F. Abe, D. Amidei, G. Apollinari, et al., “Pseudorapidity distributions of charged particles produced in interactions as and 1800 GeV,” Physical Review D, vol. 41, no. 7, pp. 2330–2333, 1990.

I. G. Bearden, D. Beavis, C. Besliu et al., “Nuclear stopping in Au + Au collisions at  GeV,” Physical Review Letters, vol. 93, no. 10, Article ID 102301, 2004.

C. Alt, T. Anticic, B. Baatar, et al., “Energy dependence of meson production in central Pb+Pb collisions at to 17 GeV,” Physical Review C, vol. 78, no. 4, Article ID 044907, 2008.

S. Chatrchyan, V. Khachatryan, A. M. Sirunyan, et al., “Missing transverse energy performance of the CMS detector,” Journal of Instrumentation, vol. 6, no. 9, Article ID P09001, 2011.

The CMS Collaboration, “Missing transverse energy performance in minimum-bias and jet events from proton-proton collisions at ,” Preprint CMS-PAS-JME-10-004, 2010

G. Aad, B. Abbott, J. Abdallah, et al., “Measurement of the jet fragmentation function and transverse profile in proton-proton collisions at a center-of-mass energy of 7 TeV with the ATLAS detector,” European Physical Journal C, vol. 71, no. 11, article 1795, 2011.

J. Adams, M. M. Aggarwal, Z. Ahammed, et al., “K(892)* resonance production in Au+Au and p+p collisions at GeV,” Physical Review C, vol. 71, no. 6, Article ID 064902, 2005.

S. Abreu, N. Borghini, S. Jeon, et al., “Heavy-ion collisions at the LHC—last call for predictions,” Journal of Physics G, vol. 35, no. 5, Article ID 054001, 2008.

K. Urmossy, “Multiplicity dependence of hadron spectra in proton-proton collisions at LHC energies and super-statistics,” http://arxiv.org/abs/1212.0260

H. R. Wei, Y. H. Chen, L. N. Gao, and F. H. Liu, “On multi-component Erlang distribution and Levy distribution of transverse momentum spectra in high-energy collisions,” submitted to. Advances in High Energy Physics

F. H. Liu, “Unified description of multiplicity distributions of final-state particles produced in collisions at high energies,” Nuclear Physics A, vol. 810, no. 1–4, pp. 159–172, 2008.

F. H. Liu and J. S. Li, “Isotopic production cross section of fragments in ⁵⁶Fe+p and ¹³⁶Xe(¹²⁴Xe)+Pb reactions over an energy range from 300A to 1500A MeV,” Physical Review C, vol. 78, no. 4, Article ID 044602, 2008.

P. Z. Ning, L. Li, and D. F. Min, Foundation of Nuclear Physics: Nucleons and Nuclei, Higher Education Press, Beijing, China, 2003.

F. H. Liu, C. X. Tian, M. Y. Duan, and B. C. Li, “Relativistic and quantum revisions of the multisource thermal model in high-energy collisions,” Advances in High Energy Physics, vol. 2012, Article ID 287521, 9 pages, 2012.

C. D. Dermer, “The production spectrum of a relativistic Maxwell-Boltzmann gas,” The Astrophysical Journal, vol. 280, no. 1, pp. 328–333, 1984.

J. L. Synge, The Relativistic Gas, North-Holland, Amsterdam, The Netherlands, 1957.

C. R. Meng, “Transverse momentum and rapidity distributions of ϕ Mesons produced in Pb-Pb collisions at SPS energies,” Chinese Physics Letters, vol. 26, no. 10, Article ID 102501, 2009.

B. I. Abelev, M. M. Aggarwal, Z. Ahammed, et al., “Systematic measurements of identified particle spectra in pp, d+Au, and Au+Au collisions at the STAR detector,” Physical Review C, vol. 79, no. 3, Article ID 034909, 2009.

H. Yang and BRAHMS Collaboration, “Identified particle production in p+p and d+Au collisions at RHIC,” Journal of Physics G, vol. 34, no. 8, pp. S619–S622, 2007.

A. Adare, S. S. Adler, S. Afanasiev, et al., “Cold nuclear matter effects on J/ψ production as constrained by deuteron-gold measurements at  GeV,” Physical Review C, vol. 77, no. 2, Article ID 024912, 2008.

A. Adare, S. Afanasiev, C. Aidala, et al., “J/ψ Production versus transverse momentum and rapidity in p+p collisions at  GeV,” Physical Review Letters, vol. 98, no. 23, Article ID 232002, 2007.

A. Bickley and PHENIX Collaboration, “Heavy quarkonia production in p+p collisions from the PHENIX experiment,” Journal of Physics G, vol. 34, no. 8, pp. S779–S782, 2007.

B. Alessandro, C. Alexa, R. Arnaldi, et al., “J/ψ and ψ′ production and their normal nuclear absorption in proton-nucleus collisions at 400 GeV,” European Physical Journal C, vol. 48, no. 2, pp. 329–341, 2006.

J. T. Mitchell and PHENIX Collaboration, “The PHENIX potential in the search for the QCD critical point,” in Proceedings of the 3rd International Workshop on the Critical Point and Onset of Deconfinement, Florence, Italy, July 2006, http://arxiv.org/abs/nucl-ex/0701079.

M. Kowalski and ALICE Collaboration, “First results on charged particle production in ALICE experiment at LHC,” Acta Physica Polonica B, vol. 42, no. 3-4, pp. 859–866, 2011.

K. Aamodt, A. Abrahantes Quintana, D. Adamova', et al., “Strange particle production in proton-proton collisions at  TeV with ALICE at the LHC,” European Physical Journal C, vol. 71, no. 3, Article ID 1594, 2011.

M. van Leeuwen and ALICE Collaboration, “High- results from ALICE,” in Proceedings of the Hadron Collider Physics Symposium, Paris, France, November 2012, http://arxiv.org/abs/1201.5205.

M. Floris and ALICE Collaboration, “Identified particles in pp and Pb-Pb collisions at LHC energies with the ALICE detector,” Journal of Physics G, vol. 38, no. 12, Article ID 124025, 2011.

H. Appelshäuser and ALICE Collaboration, “Particle production at large transverse momentum with ALICE,” Journal of Physics G, vol. 38, no. 12, Article ID 124014, 2011.

R. Preghenella and ALICE Collaboration, “Transverse momentum spectra of identified charged hadrons with the ALICE detector in Pb—Pb collisions at  TeV,” in Proceedings of the Europhysics Conference on High Energy Physics, Rhône-Alpes, France, July 2011, http://arxiv.org/abs/1111.0763.

H. Zheng, G. Giuliani, and A. Bonasera, “Coulomb corrections to density and temperature of bosons in heavy ion collisions,” http://arxiv-web3.library.cornell.edu/abs/1306.5741

J. Letessier, G. Torrieri, S. Hamieh, and J. Rafelski, “Quark-gluon plasma fireball explosion,” Journal of Physics G, vol. 27, no. 3, pp. 427–437, 2001.