On Descriptions of Particle Transverse Momentum Spectra in High Energy Collisions

Abstract

The transverse momentum spectra obtained in the frame of an isotropic emission source are compared in terms of Tsallis, Boltzmann, Fermi-Dirac, and Bose-Einstein distributions and the Tsallis forms of the latter three standard distributions. It is obtained that, at a given set of parameters, the standard distributions show a narrower shape than their Tsallis forms which result in wide and/or multicomponent spectra with the Tsallis distribution in between. A comparison among the temperatures obtained from the distributions is made with a possible relation to the Boltzmann temperature. An example of the angular distributions of projectile fragments in nuclear collisions is given.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Academic Editor:Edward Sarkisyan-Grinbaum

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

Received 6 May 2014; Revised 20 July 2014; Accepted 14 August 2014; 28 August 2014

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. The publication of this article was funded by SCOAP³ .

1. Introduction

In high energy collisions, one can use different distributions to describe the transverse momentum spectra of final-state particles. For example, in the framework of a multisource thermal model which was proposed by us some years ago [1], one can use the Rayleigh distribution [2-6], Boltzmann distribution [7, 8], Fermi-Dirac distribution [9], and Erlang distribution [10, 11] to describe the transverse momentum spectrum contributed by a given isotropic emission source. Particularly, in the classical version of the multisource thermal model, one can use the Rayleigh distribution to describe the transverse momentum spectrum, which is the same as that in the classical Maxwell ideal gas model. In the relativistic (and quantum) multisource thermal model, we can use the Boltzmann (or Fermi-Dirac or Bose-Einstein) distribution to describe the related spectrum.

In collisions at very high energy, it is important to find a proper distribution for the description of momenta or transverse momenta. Generally, the transverse momentum spectrum of final-state particles is a two-, three-, or multicomponent Boltzmann, Fermi-Dirac, Erlang, or other distributions. It shows that the types of sources produced in high energy collisions are two, three, or multiple. The rapidity spectrum of particles produced in the rest sources can be obtained by the spectrum of transverse momenta and assumption of isotropic emission, and the rapidity spectrum of final-state particles is the sum of different sources with different rapidity shifts [12].

In transverse momentum distribution, except for the contribution ratio, the temperature parameters are different for different components in the mentioned two-, three-, or multicomponent distributions. These two-, three-, or multicomponent distributions render temperature fluctuations in high energy collisions, which can be described by the nonextensive Tsallis distribution [13-22] or the Tsallis forms of the standard (Boltzmann, Fermi-Dirac, and Bose-Einstein) distributions, and the temperature fluctuations can be described by a parameter (an entropic index) which characterizes the degree of nonequilibrium. In fact, the Tsallis distribution describes the final-state transverse momentum spectra by using two main parameters, a temperature T and a nonequilibrium degree (or an entropic index) q , where q>1 and q=1 mean a nonequilibrium state and an equilibrium state, respectively [13-22]. The equilibrium indicates a thermal and kinetic equilibrium and does not require a chemical equilibrium.

In experiments at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) [23-28], the measured transverse momentum spectra are found to obey the Tsallis distribution [13-22] or similar distribution such as the Lévy distribution [11, 29-31], where the latter one is a special expression of the former one. The abundant data obtained at the RHIC and LHC have attracted many theoretical or modelling analyses [29-38], and some analyses show two or multiple processes' emission [34, 37]. In some cases [30, 31], the Lévy distribution is also called the Tsallis distribution. The Tsallis distribution has more than one version used in literature [13-38], in which the power laws are somewhat different.

Different distributions can be adopted in the multisource thermal model [1]. Generally, for a given isotropic emission source, the temperature parameter should be a given value. In the case of using the same temperature parameter for different distributions, different results can be obtained. It is an important issue to choose a proper distribution for the rest sources in the multisource thermal model. In most cases, there is more than one distribution to be chosen. We are interested in studying the differences and similarities among these distributions of transverse momentum spectra. These differences and similarities are useful for the descriptions of the rest sources in the multisource thermal model [1].

In this paper, we compare the Tsallis, Boltzmann, Fermi-Dirac, Bose-Einstein, Tsallis form of Boltzmann (Tsallis-B), Tsallis form of Fermi-Dirac (Tsallis-FD), and Tsallis form of Bose-Einstein (Tsallis-BE) distributions expected to be observed in high energy collisions. Following the formalism of these distributions in Section 2 which is based on [38-44], a comparison and example are given in Sections 3 and 4, respectively. Finally, we give our conclusions in Section 5.

2. The Formalism

The Tsallis distribution has more than one version [13-38]. In the present work, we adopt the version used in [38]. Correspondingly, the Tsallis forms of the three standard distributions have different expressions [39-44]. Based on [38], we can combine the different expressions for the given version to a uniform one. On the invariant particle momentum (p ) distribution, we have [figure omitted; refer to PDF] where N is the particle number, g is the degeneracy factor, and V is the volume; E , _m0 , and μ are the energy, the rest mass, and the chemical potential of the considered particle, respectively; S=0 , +1 , and -1 denote the Boltzmann, Fermi-Dirac, and Bose-Einstein distributions, respectively; a ± sign is introduced in the above equation, + is for ^p2 +^m02 >μ , and - is for ^p2 +^m02 ...4;μ [38]. Then, the normalized Tsallis forms of the three standard distributions for momenta can be given by [figure omitted; refer to PDF] where _C1 is the normalization constant. Further, we have the Tsallis forms of the three standard unit-density functions of transverse momentum (_pT ) and rapidity (y ) to be [figure omitted; refer to PDF] where a ± sign is introduced in the above equation, + is for ^pT2 +^m02 cosh...y>μ , and - is for ^pT2 +^m02 cosh...y...4;μ ; the Tsallis forms of the three standard transverse momentum distributions are [figure omitted; refer to PDF] where _ymax... and _ymin... denote the maximum and minimum rapidities, respectively; and the Tsallis forms of the three standard rapidity distributions are [figure omitted; refer to PDF] where _pTmax... denotes the maximum transverse momentum.

In particular, some special distributions can be obtained due to some limitations. According to [38], the Tsallis distribution which describes well experimental data of invariant particle momenta in high energy collisions can be given by [figure omitted; refer to PDF] where, comparing with (1), the limitations ^p2 +^m02 >μ and S=0 and the upper index q/(q-1) are used in (6). Then, we have the normalized Tsallis distribution of momenta to be [figure omitted; refer to PDF] where _C2 is the normalization constant. The unit-density function of _pT and y is [38] [figure omitted; refer to PDF] Further, we have the transverse momentum distribution to be [figure omitted; refer to PDF] and the rapidity distribution to be [figure omitted; refer to PDF] Because q is very close to 1, we have the upper index q/(q-1) in (6)-(10) to approximate 1/(q-1) . The upper index 1/(q-1) is a result of a mean field approach [45] which gives a smaller q too close to 1.

Particularly, in the limit where q[arrow right]1 , the Tsallis distribution reduces to the standard Boltzmann distribution [38]. Considering the standard Fermi-Dirac distribution and the standard Bose-Einstein distribution, we uniformly have the three standard distributions to be [figure omitted; refer to PDF] The normalized form of the three standard distributions for momenta can be given by [46, 47] [figure omitted; refer to PDF] where _C3 is the normalization constant. The unit-density function of _pT and y is [38] [figure omitted; refer to PDF] The transverse momentum distribution can be given by [figure omitted; refer to PDF] The rapidity distribution is [figure omitted; refer to PDF]

In the above discussions, we have used the same temperature value in different distributions and obtained different results. For a given transverse momentum spectrum, different distributions can be used at different temperatures to give the same or similar results. The relations between the Tsallis temperature _TT and Boltzmann temperature (physical temperature) _TB , Fermi-Dirac temperature _TFD and _TB , Bose-Einstein temperature _TBE and _TB , Tsallis-B temperature _TTB and _TB , Tsallis-FD temperature _TTFD and _TB , and Tsallis-BE temperature _TTBE and _TB are unclear. The simplest relations are approximately linear, especially for light and heavy flavor particles, respectively [37]. We have [figure omitted; refer to PDF] respectively, where _{aT,FD,BE,TB,TFD,TBE} and _{CT,FD,BE,TB,TFD,TBE} are corrected factors in the relations. Naturally, we do not eliminate other possibilities on the relations.

We regard _TB as the physical temperature due to the fact that _TB can be obtained from the relativistic ideal gas model and it degenerates to the temperature in the classical ideal gas model. To give a detailed comparison for _TB and other temperatures in experiments goes beyond the focus of the present work. The relations will be investigated in more detail using the data from RHIC and LHC. The above linear relations are expected to be approximately correct for light particles such as pions, kaons, and nucleons or they are expected to be approximately correct for heavy particles such as J/ψ , Λ , and Δ . The slopes corresponding to light and heavy particles are different. For totally light and heavy particles, the relations seem to be two-linear [37].

3. Comparison

Before we give formally the figures, we have done some calculations for purposes of testing and checking. It is found that, in high energy collisions, the effect with μ...a;_m0 can be neglected as we have known in a previous work [47]. Only in the case of μ[arrow right]_m0 , one can see some differences between quantum and nonquantum distributions. In the formal calculation, we take _m0 =138 MeV/c and μ=130 MeV at a given temperature in the present work. According to the results of the thermal dynamic model [48], the critical temperature for phase transition from hadron matter to quark matter is 156 MeV. We take T=160 MeV for heavy ion collisions at the RHIC and LHC. In addition, we take q=1.05 , 1.1, 1.2, and 1.3 for comparisons, which will make a large difference in distributions. The selected values of q cover approximately possible range obtained from experimental data [37, 38].

Figure 1 presents the transverse momentum spectra calculated by the above several distributions. The solid curves, open circles, dashed and dotted curves which neared the open circles, open squares, and dashed and dotted curves which neared the open squares correspond to the Tsallis, Boltzmann, Fermi-Dirac, Bose-Einstein, Tsallis-B, Tsallis-FD, and Tsallis-BE distributions, respectively. In the calculation, we take _m0 =138 MeV/c, μ=130 MeV, T=160 MeV, and q=1.05 (Figure 1(a)), 1.1 (Figure 1(b)), 1.2 (Figure 1(c)), and 1.3 (Figure 1(d)), respectively, as discussed above. Figure 2 is the same as that for Figure 1, but it shows the calculated results for rapidity spectra. (To check the validity of the numerical results in Figures 1 and 2, we have performed some additional Monte Carlo tests for the Tsallis distribution and found that the two results are in agreement with each other).

Transverse momentum spectra described by several distributions for _m0 =138 MeV/c , μ=130 MeV, T=160 MeV, and q=1.05 (a), 1.1 (b), 1.2 (c), and 1.3 (d), respectively. The solid curves, open circles, dashed and dotted curves which neared the open circles, open squares, and dashed and dotted curves which neared the open squares correspond to the Tsallis, Boltzmann, Fermi-Dirac, Bose-Einstein, Tsallis-B, Tsallis-FD, and Tsallis-BE distributions, respectively.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

It is the same as for Figure 1, but it shows the rapidity spectra.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

One can see from Figures 1 and 2 that, at a given temperature, the three standard distributions result in narrow spectra; the Tsallis-B, Tsallis-FD, and Tsallis-BE distributions result in wide spectra; and the Tsallis distribution results belong to the middle spectra. The quantum effect, that is, the difference in the results for different S , can be observed in the case of using the given set of parameters. One can also expect that the two- or three-temperature standard distributions fit approximately the Tsallis distribution, and the multitemperature standard distribution and the two- or three-temperature Tsallis distribution fit approximately the Tsallis forms of the standard distributions. Particularly, if the Tsallis distribution describes the temperature fluctuations in wide transverse momentum spectrum which can be fitted by a two- or three-component Boltzmann distribution, the Tsallis forms of the standard distributions describe large temperature fluctuations in wider transverse momentum spectrum which can be fitted by a multicomponent Boltzmann distribution.

In fact, our previous work [49] shows that one can use the multisource thermal model [1] to describe the transverse momentum spectra of particles produced in soft process in proton-proton and nucleus-nucleus collisions at RHIC and LHC energies. The standard Boltzmann distribution can be used to describe the particle behavior for single source. Then, the experimental transverse momentum spectrum is described by single, two-, or three-component Boltzmann distribution. If we consider the hard process which results in a wider transverse momentum spectrum, the multicomponent Boltzmann distribution is needed. In [37] and by us in [22], the Tsallis distribution is shown to describe successfully the transverse momentum (mass) spectra in high energy collisions. One can use the Tsallis distribution to describe uniformly the single, two-, and three-component Boltzmann distributions. A multicomponent Boltzmann distribution describes a wider transverse momentum spectrum which can be fitted by the Tsallis forms of the standard distributions or by a two- or three-component Tsallis distribution.

Figure 3 shows the comparisons of different distributions with Boltzmann distribution on the same or similar transverse momentum spectrum. In the comparisons, we take _TB =80 (Figure 3(a)), 160 (Figure 3(b)), 240 (Figure 3(c)), and 320 MeV (Figure 3(d)), respectively; and _m0 =138 MeV/c² , μ=0 , and q=1.012 for the four cases. The meanings of different curves and symbols are the same as those for Figure 1. To give the same or similar result as the Boltzmann distribution, the selected temperatures for other distributions are listed in Table 1. One can see that, for the same or similar transverse momentum distribution, comparing with the Boltzmann temperature (physical temperature), the standard Bose-Einstein distribution measures a high temperature, and other distributions measure low temperatures, for a small value of q . It is difficult to get the same or similar results by different distributions with proper temperatures for a large q .

Table 1: Values of parameters corresponding to different distributions in Figure 3. The temperatures are in the units of MeV.

Figure	T B	T T	T FD	T BE	T TB	T T FD	T T BE	q
Figure 3(a)	80	75	79	81	75	74	75	1.012
Figure 3(b)	160	150	156	163	150	147	150	1.012
Figure 3(c)	240	229	236	247	224	223	230	1.012
Figure 3(d)	320	305	310	330	302	298	310	1.012

Comparisons of transverse momentum spectra for different distributions. Figures 3(a)-3(d) correspond to different _TB indicated in the panels and a given q of 1.012. The meanings of the curves and symbols are the same as those for Figure 1.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

Different temperatures extracted from the same or similar transverse momentum spectrum result in different thermodynamics and statistics. The most final-state particles are fermions, and the Fermi-Dirac statistics is expected to be the most appropriate description for single source. In high energy collisions, the quantum effect is small and the Boltzmann distribution is expected to be the most appropriate refill. Different thermodynamics and statistics give different expressions of transverse momentum spectra, which extract different temperatures. For a given system, the present work shows that the temperature found from Boltzmann distribution is smaller when compared with that from Bose-Einstein distribution and is larger when compared with those from other distributions.

Although comparing the distributions to some data in proton-proton, proton-nucleus, or nucleus-nucleus collisions would lead to some interesting conclusions and we have not compared the present results with experimental data [23-28], the experimental groundwork of the present work is reliable. In fact, we have used the single or two-component Boltzmann distribution and the Tsallis distributions to describe experimental data in our previous work [22, 49-51]. The extracted temperatures are in the range used in the present work. For a narrow spectrum, all the distributions are very similar for a specific choice of temperatures and a small value of entropic index. We can conclude that all the distributions are equally good to describe the narrow spectrum. Because high energy collisions contain abundant contents such as entropy density [52], phase transition [53], and flow effects [54], more careful and broad investigations are needed in future.

4. An Example: The Angular Distribution of Projectile Fragments in Nuclear Collisions

In the above discussions, we have used the rest frame of the considered source and studied the behaviors of particles produced in violent interacting region such as the participant part in nucleus-nucleus collisions. In fact, these formulas can be also used in nonviolent interacting region such as the spectator parts in nucleus-nucleus collisions. As an example, we study here the space angular distributions of projectile fragments in fixed target experiments in the laboratory reference frame.

Let _pbeam , _AF , and Θ denote the beam momentum per nucleon, the nucleon number of a considered projectile fragment, and the emission angle of the projectile fragment, respectively. The momentum per nucleon of the projectile fragment is assumed to be _pbeam . We have the transverse momentum of the projectile fragment to be _pT =_AFpbeam sinΘ . Because _pT can obey the distributions discussed in Section 2, we have the Θ distributions as in the following.

The Tsallis forms of the three standard distributions for the emission angles of projectile fragments are [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and _mN denotes the rest mass of a nucleon; a ± sign is introduced in the above equation, + is for _AF^pbeam2sin2 Θ+^mN2 cosh...y>μ , and - is for _AF^pbeam2sin2 Θ+^mN2 cosh...y...4;μ . The Tsallis distribution for the emission angles of projectile fragments can be given by [figure omitted; refer to PDF] where, comparing with (22), the limitations _AF^pbeam2sin2 Θ+^mN2 cosh...y>μ and S=0 and the upper index q/(q-1) are used in (24). One can also use the upper index 1/(q-1) in (24). The three standard distributions for the emission angles of projectile fragments are [figure omitted; refer to PDF] where S=0 or +1 for projectile fragments.

Particularly, in the whole phase space, the classical ideal gas model gives the probability distribution of emission angles for the projectile fragments to be [figure omitted; refer to PDF] due to the fact that the probability distribution of transverse momenta obeys [figure omitted; refer to PDF] which is a Rayleigh distribution with the dispersion _σp =_AFmN T . The relativistic ideal gas model gives the probability distribution of emission angles for the projectile fragments to be [figure omitted; refer to PDF] due to the fact that the probability distribution of transverse momenta obeys [figure omitted; refer to PDF] which is the simplest Boltzmann distribution, where the chemical potential and quantum effects are neglected due to small values.

In high energy nucleus-nucleus collisions, we have sinΘ[approximate]Θ and sin2Θ[approximate]2Θ for the projectile fragments. Because the spectators are usually in a low excitation state, the classical ideal gas model can be used for the projectile fragment. Then, the Θ distribution of the projectile fragments can be approximately expressed by [figure omitted; refer to PDF] which is a Rayleigh distribution with the dispersion _σΘ =_mN T/(_AF^pbeam2 ) .

In some cases, we have to use a two- or three-temperature picture to give descriptions for the projectile fragments [5]. The angular distribution is then the sum of two or three Rayleigh distributions with different weights and dispersions (temperatures), and the transverse momentum distribution is also the sum of two or three Rayleigh distributions with different weights and dispersions (temperatures). For the temperature fluctuations which appeared in the projectile fragment distributions, we can also use the Tsallis statistics to give descriptions.

5. Conclusions

From the above discussions, we obtain the following conclusions.

(a) In the rest frame of an isotropic emission source, the spectra described by Tsallis, Boltzmann, Fermi-Dirac, Bose-Einstein, Tsallis form of Boltzmann, Tsallis form of Fermi-Dirac, and Tsallis form of Bose-Einstein distributions are compared in the present work.

(b) The quantum effect can be observed in the case of using the given set of parameters for Figures 1 and 2. At a given temperature, the three standard (Boltzmann, Fermi-Dirac, and Bose-Einstein) distributions result in narrow spectra, the Tsallis forms of the three standard distributions (Tsallis-B, Tsallis-FD, and Tsallis-BE distributions) result in wide spectra, and the Tsallis distribution results in between.

(c) For a given transverse momentum spectrum, the Bose-Einstein distribution shows higher temperature than that by the Boltzmann distribution (physical temperature), while other distributions show lower temperatures. The possible relations between the Boltzmann temperature and others are given by using the corrected factors.

(d) The related distributions can be also used in the descriptions of projectile fragments. Generally, the transverse momentum spectra and angular distributions of projectile fragments obey the Rayleigh distribution. One needs two or three Rayleigh distributions to describe the projectile fragment distributions. The temperature fluctuations which appeared in the projectile fragment distributions can be also described by the Tsallis statistics.

Acknowledgments

The authors acknowledge the referees and editor for their constructive comments. This work was partly finished at the State University of New York at Stony Brook, USA. One of the authors (Fu-Hu Liu) thanks Professor Dr. Roy A. Lacey and the members of the Nuclear Chemistry Group of Stony Brook University for their hospitality. The authors also acknowledge the supports of the National Natural Science Foundation of China (under Grant no. 10975095), the Open Research Subject of the Chinese Academy of Sciences Large-Scale Scientific Facility (under Grant no. 2060205), the Shanxi Scholarship Council of China, and the Overseas Training Project for Teachers at Shanxi University.

Conflict of Interests

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

References

[1] F.-H. Liu, "Particle production in Au-Au collisions at RHIC energies," Physics Letters B , vol. 583, no. 1-2, pp. 68-72, 2004.

[2] N. P. Andreeva, Z. V. Anzon, V. I. Bubnov, "Fragmentation of ²² Ne relativistic nuclei on photoemulsion nuclei," Soviet Journal of Nuclear Physics , vol. 47, no. 1, pp. 102-108, 1988.

[3] S. Kamel, "Transverse-momentum distribution of ultrarelativistic protons in electromagnetic dissociation of ³² S into p + ³¹ P," Nuovo Cimento della Societa Italiana di Fisica A , vol. 112, no. 4, pp. 327-334, 1999.

[4] F.-H. Liu, "Pseudorapidity (rapidity) and transverse momentum distributions of negatively charged hadrons in high energy nucleus-nucleus collisions," Chinese Journal of Physics , vol. 38, no. 1, pp. 42-48, 2000.

[5] F.-H. Liu, "Emission of relativistic light fragments in nucleus-emulsion collisions at high energy," Chinese Journal of Physics , vol. 38, no. 6, pp. 1063-1073, 2000.

[6] F.-H. Liu, "Transverse momentum distribution of strange particles in high-energy Pb-Pb collisions," Chinese Journal of Physics , vol. 41, no. 4, pp. 357-371, 2003.

[7] K. B. Bhalla, M. Chaudhry, S. Lokanathan, R. K. Grover, I. K. Daftari, G. L. Kaul, L. K. Magotra, Y. Prakash, N. K. Rao, S. Garpman, I. Otterlund, "Relativistic α-particles emitted in Fe-emulsion interactions at 1.7A GeV," Nuclear Physics A , vol. 367, no. 3, pp. 446-458, 1981.

[8] M. M. Aggarwal, K. B. Bhalla, G. Das, P. L. Jain, "Angular distributions of relativistic alpha particles in heavy-ion collisions," Physical Review C , vol. 27, no. 2, pp. 640-649, 1983.

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

[10] E. Q. Wang, H. R. Wei, M. A. Rahim, S. Fakhraddin, F. H. Liu, "Multi-component Erlang distribution of final-state particles produced in high energy collisions," Indian Journal of Physics , vol. 87, no. 2, pp. 185-193, 2013.

[11] H.-R. Wei, Y.-H. Chen, L.-N. Gao, F.-H. Liu, "Comparing multicomponent Erlang distribution and Lévy distribution of particle transverse momentum," Advances in High Energy Physics , vol. 2014, 2014.

[12] F.-H. Liu, Y.-Q. Gao, T. Tian, B.-C. Li, "Transverse momentum and pseudorapidity distributions of charged particles and spatial shapes of interacting events in Pb-Pb collisions at 2.76 TeV," Advances in High Energy Physics , 2014.

[13] C. Tsallis, "Possible generalization of Boltzmann-Gibbs statistics," Journal of Statistical Physics , vol. 52, no. 1-2, pp. 479-487, 1988.

[14] C. Tsallis, "Nonadditive entropy and nonextensive statistical mechanics--an overview after 20 years," Brazilian Journal of Physics , vol. 39, no. 2A, pp. 337-356, 2009.

[15] C. Tsallis, R. S. Mendes, A. R. Plastino, "The role of constraints within generalized nonextensive statistics," Physica A , vol. 261, no. 3-4, pp. 534-554, 1998.

[16] W. M. Alberico, A. Lavagno, "Non-extensive statistical effects in high-energy collisions," The European Physical Journal A , vol. 40, no. 3, pp. 313-323, 2009.

[17] W. M. Alberico, P. Czerski, A. Lavagno, M. Nardi, V. Somá, "Signals of non-extensive statistical mechanics in high-energy nuclear collisions," Physica A , vol. 387, no. 2-3, pp. 467-475, 2008.

[18] T. S. Biró Is There a Temperature? Conceptual Challenges at High Energy, Acceleration and Complexity , vol. 171, of Fundamental Theories of Physics, Springer, Berlin, Germany, 2011.

[19] T. S. Biró, G. Purcsel, K. Ürmössy, "Non-extensive approach to quark matter," The European Physical Journal A , vol. 40, no. 3, pp. 325-340, 2009.

[20] G. Wilk, Z. Wlodarczyk, "Power laws in elementary and heavy-ion collisions : a story of fluctuations and nonextensivity?," The European Physical Journal A , vol. 40, no. 3, pp. 299-312, 2009.

[21] G. Wilk, Z. Wlodarczyk, "Multiplicity fluctuations due to the temperature fluctuations in high-energy nuclear collisions," Physical Review C , vol. 79, no. 5, 2009.

[22] B.-C. Li, Y.-Z. Wang, F.-H. Liu, "Formulation of transverse mass distributions in Au-Au collisions at _sNN =200 GeV/nucleon," Physics Letters B , vol. 725, no. 4-5, pp. 352-356, 2013.

[23] A. Adare, S. Afanasiev, C. Aidala, "Identified charged hadron production in p+p collisions at s=200 and 62.4 GeV," Physical Review C , vol. 83, no. 6, 2011.

[24] B. I. Abelev, J. Adams, M. M. Aggarwal, "Strange particle production in p+p collisions at s=200 GeV," Physical Review C , vol. 75, no. 6, 2007.

[25] K. Aamodt, N. Abel, U. Abeysekara, "Transverse momentum spectra of charged particles in proton-proton collisions at s=900 GeV with ALICE at the LHC," Physics Letters B , vol. 693, no. 2, pp. 53-68, 2010.

[26] G. Aad, B. Abbott, J. Abdallah, "Charged-particle multiplicities in pp interactions measured with the ATLAS detector at the LHC," New Journal of Physics , vol. 13, no. 5, 2011.

[27] V. Khachatryan, A. M. Sirunyan, A. Tumasyan, "Transverse-momentum and pseudorapidity distributions of charged hadrons in pp collisions at s=0.9 and 2.36 TeV," Journal of High Energy Physics , vol. 2010, no. 2, 2010.

[28] V. Khachatryann, A. M. Sirunyan, A. Tumasyan, "Transverse-momentum and pseudorapidity distributions of charged hadrons in pp collisions at s=7 TeV," Physical Review Letters , vol. 105, no. 2, 2010.

[29] G. Wilk, Z. Wlodarczyk, "Interpretation of the nonextensivity parameter q in some applications of Tsallis statistics and Levy distributions," Physical Review Letters , vol. 84, no. 13, 2000.

[30] J. Cleymans, G. I. Lykasov, A. S. Parvan, A. S. Sorin, O. V. Teryaev, D. Worku, "Systematic properties of the Tsallis distribution: energy dependence of parameters in high energy p-p collisions," Physics Letters B , vol. 723, no. 4-5, pp. 351-354, 2013.

[31] J. Cleymans, "The Tsallis distribution for p -p collisions at the LHC," Journal of Physics: Conference Series , vol. 455, 2013.

[32] M. D. Azmi, J. Cleymans, "Transverse momentum distributions in proton-proton collisions at LHC energies and Tsallis thermodynamics," Journal of Physics G: Nuclear and Particle Physics , vol. 41, no. 6, 2014.

[33] M. D. Azmi, J. Cleymans, "Transverse momentum distributions in p-Pb collisions and Tsallis thermodynamics," http://arxiv.org/abs/1311.2909

[34] C. Wong, G. Wilk, "Tsallis fits to p_T spectra and multiple hard scattering in pp collisions at the LHC," Physical Review D , vol. 87, no. 11, 2013.

[35] C. Y. Wong, G. Wilk, "Tsallis fits to p_T spectra for pp collisions at the LHC," Acta Physica Polonica B , vol. 43, no. 11, pp. 2047-2054, 2012.

[36] B.-C. Li, Y.-Z. Wang, F.-H. Liu, X.-J. Wen, Y.-E. Dong, "Particle production in relativistic pp (p- ) and AA collisions at RHIC and LHC energies with Tsallis statistics using the two-cylindrical multisource thermal model," Physical Review D , vol. 89, no. 5, 2014.

[37] F.-H. Liu, Y.-Q. Gao, B.-C. Li, "Comparing two-Boltzmann distribution and Tsallis statistics of particle transverse momentums in collisions at LHC energies," The European Physical Journal A , vol. 50, no. 8, 2014.

[38] J. Cleymans, D. Worku, "Relativistic thermodynamics: transverse momentum distributions in high-energy physics," European Physical Journal A , vol. 48, article 160, 2012.

[39] F. Büyükkiliç, D. Demirhan, "A fractal approach to entropy and distribution functions," Physics Letters A , vol. 181, no. 1, pp. 24-28, 1993.

[40] F. Pennini, A. Plastino, A. R. Plastino, "Tsallis entropy and quantal distribution functions," Physics Letters. A , vol. 208, no. 4-6, pp. 309-314, 1995.

[41] J.-C. Chen, Z.-P. Zhang, G.-Z. Su, L.-X. Chen, Y.-G. Shu, " q -generalized Bose-Einstein condensation based on Tsallis entropy," Physics Letters A , vol. 300, no. 1, pp. 65-70, 2002.

[42] A. M. Teweldeberhan, A. R. Plastino, H. G. Miller, "On the cut-off prescriptions associated with power-law generalized thermostatistics," Physics Letters A , vol. 343, no. 1-3, pp. 71-78, 2005.

[43] J. M. Conroy, H. G. Miller, "Color superconductivity and Tsallis statistics," Physical Review D , vol. 78, no. 5, 2008.

[44] J. M. Conroy, H. G. Miller, A. R. Plastino, "Thermodynamic consistency of the q-deformed Fermi-Dirac distribution in nonextensive thermostatics," Physics Letters A , vol. 374, no. 45, pp. 4581-4584, 2010.

[45] A. Lavagno, "Hot and dense hadronic matter in an effective mean-field approach," Physical Review C , vol. 81, no. 4, 2010.

[46] C. D. Dermer, "The production spectrum of a relativistic Maxwell-Boltzmann gas," The Astrophysical Journal , vol. 280, no. 1, part 1, pp. 328-333, 1984.

[47] F.-H. Liu, C.-X. Tian, M.-Y. Duan, B.-C. Li, "Relativistic and quantum revisions of the multisource thermal model in high-energy collisions," Advances in High Energy Physics , vol. 2012, 2012.

[48] P. Braun-Munzinger, "Hadron production in nuclear collisions and the QCD phase boundary," Talk given at 45 years of nuclear theory at Stony Brook: a tribute to Gerald E. Brown, Stony Brook University, Stony Brook, NY, USA, (will be published as a special issue in Nuclear Physics A), 2013, http://tonic.physics.sunysb.edu/gerrybrown/program2.html

[49] F.-H. Liu, Y.-H. Chen, H.-R. Wei, B.-C. Li, "Transverse momentum distributions of final-state particles produced in soft excitation process in high energy collisions," Advances in High Energy Physics , vol. 2013, 2013.

[50] F.-H. Liu, T. Tian, X.-J. Wen, "Transverse mass and rapidity distributions and space dispersion plots of (net-)protons in Pb-Pb collisions at SPS energies," The European Physical Journal A , vol. 50, 2014.

[51] F.-H. Liu, T. Tian, H. Zhao, B.-C. Li, "Extracting chemical potentials of quarks from ratios of negatively/positively charged particles in high-energy collisions," The European Physical Journal A , vol. 50, no. 3, 2014.

[52] B.-C. Li, L.-L. Wang, "Entropy in scalar O (N ) model with a spontaneously broken symmetry," International Journal of Modern Physics A , vol. 24, no. 30, pp. 5725-5736, 2009.

[53] B.-C. Li, M. Huang, "Strongly coupled matter near phase transition," Journal of Physics G , vol. 36, no. 6, 2009.

[54] B.-C. Li, Y.-Y. Fu, L.-L. Wang, F.-H. Liu, "Dependence of elliptic flows on transverse momentum and number of participants in Au+Au collisions at sNN=200 GeV," Journal of Physics G , vol. 40, no. 2, 2013.