On Descriptions of Particle Transverse Momentum Spectra in High Energy Collisions

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.


Introduction
In high energy collisions, one can use different distributions to describe the transverse momentum spectra of finalstate 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][3][4][5][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][14][15][16][17][18][19][20][21][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 2 Advances in High Energy Physics spectra by using two main parameters, a temperature  and a nonequilibrium degree (or an entropic index) , where  > 1 and  = 1 mean a nonequilibrium state and an equilibrium state, respectively [13][14][15][16][17][18][19][20][21][22].The equilibrium indicates a thermal and kinetic equilibrium and does not require a chemical equilibrium.
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][39][40][41][42][43][44], a comparison and example are given in Sections 3 and 4, respectively.Finally, we give our conclusions in Section 5.

The Formalism
The Tsallis distribution has more than one version .In the present work, we adopt the version used in [38].Correspondingly, the Tsallis forms of the three standard distributions have different expressions [39][40][41][42][43][44].Based on [38], we can combine the different expressions for the given version to a uniform one.On the invariant particle momentum () distribution, we have where  is the particle number,  is the degeneracy factor, and  is the volume; ,  0 , and  are the energy, the rest mass, and the chemical potential of the considered particle, respectively;  = 0, +1, and −1 denote the Boltzmann, Fermi-Dirac, and Bose-Einstein distributions, respectively; a ± sign is introduced in the above equation, + is for √ 2 +  2 0 > , and − is for √ 2 +  2 0 ≤  [38].Then, the normalized Tsallis forms of the three standard distributions for momenta can be given by where  1 is the normalization constant.Further, we have the Tsallis forms of the three standard unit-density functions of transverse momentum (  ) and rapidity () to be where a ± sign is introduced in the above equation, + is for √ 2  +  2 0 cosh  > , and − is for √ 2  +  2 0 cosh  ≤ ; the Tsallis forms of the three standard transverse momentum distributions are Advances in High Energy Physics 3 where  max and  min denote the maximum and minimum rapidities, respectively; and the Tsallis forms of the three standard rapidity distributions are where   max 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 where, comparing with (1), the limitations √ 2 +  2 0 >  and  = 0 and the upper index /( − 1) are used in (6).Then, we have the normalized Tsallis distribution of momenta to be where  2 is the normalization constant.The unit-density function of   and  is [38] Further, we have the transverse momentum distribution to be and the rapidity distribution to be Because  is very close to 1, we have the upper index /( − 1) in ( 6)-( 10) to approximate 1/(−1).The upper index 1/(−1) is a result of a mean field approach [45] which gives a smaller  too close to 1.
Particularly, in the limit where  → 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 The normalized form of the three standard distributions for momenta can be given by [46,47] where  3 is the normalization constant.The unit-density function of   and  is [38] The transverse momentum distribution can be given by The rapidity distribution is 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   and Boltzmann temperature (physical temperature)   , Fermi-Dirac temperature  FD and   , Bose-Einstein temperature  BE and   , Tsallis-B temperature  TB and   , Tsallis-FD temperature  TFD and   , and Tsallis-BE temperature  TBE and   are unclear.The simplest relations are approximately linear, especially for light and heavy flavor particles, respectively [37].We have respectively, where  ,FD,BE,TB,TFD,TBE and  ,FD,BE,TB,TFD,TBE are corrected factors in the relations.Naturally, we do not eliminate other possibilities on the relations.We regard   as the physical temperature due to the fact that   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   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 /, Λ, 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].

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  ≪  0 can be neglected as we have known in a previous work [47].Only in the case of  →  0 , one can see some differences between quantum and nonquantum distributions.In the formal calculation, we take  0 = 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  = 160 MeV for heavy ion collisions at the RHIC and LHC.In addition, we take  = 1.05, 1.1, 1.2, and 1.3 for comparisons, which will make a large difference in distributions.The selected values of  cover approximately possible range obtained from experimental data [37,38].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).
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 Advances in High Energy Physics the results for different , can be observed in the case of using the given set of parameters.One can also expect that the twoor 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   = 80 (Figure 3  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][24][25][26][27][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][50][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.

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  beam ,   , 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  beam .We have the transverse momentum of the projectile fragment to be   =    beam sin Θ.Because   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 where and   denotes the rest mass of a nucleon; a ± sign is introduced in the above equation, + is for   √ 2 beam sin 2 Θ +  2  cosh  > , and − is for for the emission angles of projectile fragments can be given by where, comparing with ( 22), the limitations   √ 2 beam sin 2 Θ +  2  cosh  >  and  = 0 and the upper index /( − 1) are used in (24).One can also use the upper index 1/(−1) in (24).The three standard distributions for the emission angles of projectile fragments are where  = 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 due to the fact that the probability distribution of transverse momenta obeys which is a Rayleigh distribution with the dispersion   = √    .The relativistic ideal gas model gives the probability distribution of emission angles for the projectile fragments to be due to the fact that the probability distribution of transverse momenta obeys 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 Θ ≈ Θ and sin 2Θ ≈ 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 which is a Rayleigh distribution with the dispersion  Θ = √  /(   2 beam ).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.

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.

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

Figure 3 :
Figure 3: Comparisons of transverse momentum spectra for different distributions.Figures 3(a)-3(d) correspond to different   indicated in the panels and a given  of 1.012.The meanings of the curves and symbols are the same as those for Figure 1.
(a)), 160 (Figure 3(b)), 240 (Figure 3(c)), and 320 MeV (Figure 3(d)), respectively; and  0 = 138 MeV/c 2 ,  = 0, and  = 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

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

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 .It is difficult to get the same or similar results by different distributions with proper temperatures for a large .