Hadron Spectra Parameters within the Non-Extensive Approach

TL;DR

This paper examines the application of non-extensive statistical mechanics to high-energy hadron spectra, revealing scaling behaviors and correlations of temperature and non-extensive parameter with collision energy.

Contribution

It demonstrates the transferability of non-extensive distributions and identifies explicit scaling laws for temperature and q in different collision systems.

Findings

Non-extensive distributions can be transformed into each other.

Temperature scales linearly with collision energy in proton-proton collisions.

Temperature and q parameter are linearly correlated in nucleus-nucleus collisions.

Abstract

We investigate how the non-extensive approach works in high-energy physics. Transverse momentum ($p_T$) spectra of several hadrons are fitted by various non-extensive momentum distributions and by the Boltzmann--Gibbs statistics.~It is shown that some non-extensive distributions can be transferred one into another.~We find explicit hadron mass and center-of-mass energy {scaling both in the temperature and in the non-extensive parameter, $q$,} in proton--proton and heavy-ion collisions. We find that the temperature depends linearly, but the Tsallis $q$ follows a logarithmic dependence on the collision energy in proton--proton collisions. In the nucleus--nucleus collisions, on~the other hand, $T$ and $q$ correlate linearly, as was predicted in our previous work.

Tables

Table 1. Table 1: The values of χ 2 / d . o . f . formulae-sequence superscript 𝜒 2 𝑑 𝑜 𝑓 \chi^{2}/d.o.f. of spectral fits for pions, kaons, and protons in p ​ p 𝑝 𝑝 pp collisions at 900 GeV and 7 TeV as examples.

Collision Energy ($\sqrt{s}$)	Produced Hadrons	${𝒇}_{𝑩𝑮}$	${𝒇}_{\mathrm{𝟏}}$	${𝒇}_{\mathrm{𝟐}}$	${𝒇}_{\mathrm{𝟑}}$	${𝒇}_{\mathrm{𝟒}}$	${𝒇}_{\mathrm{𝟓}}$
	$\pi$	110.8	0.2814	0.2814	0.4697	0.2814	1.456
900 GeV	$K$	8.047	0.1748	0.1749	0.1698	0.1749	0.6669
	$p$	3.491	0.3724	0.3724	0.3735	0.3724	0.4145
	$\pi$	1316.0	0.9681	0.9681	3.417	0.9681	0.3049
7 TeV	$K$	520.2	0.4202	0.4202	0.4313	0.4202	3.100
	$p$	254.3	0.4481	0.4481	0.4356	0.4481	4.357

Table 2. Table 2: Fitting p T subscript 𝑝 𝑇 p_{T} ranges of spectra for different charged particles in p ​ p 𝑝 𝑝 pp collisions  ALICE ( 2011 , 2014 ); PHENIX ( 2011 ); ALICE ( 2014 ) .

$\sqrt{𝒔}$	$𝝅$ [GeV]	$𝑲$ [GeV]	$𝒑$ [GeV]
62.4 GeV	0.3–2.9	0.4–2	0.6–3.6
200 GeV	0.3–3	0.4–2	0.5–4.6
900 GeV	0.1–2.6	0.2–2.4	0.35–2.4
2.76 TeV	0.1–20	0.2–20	0.3–20
5.02 TeV	0.1–20	0.2–20	0.3–20
7 TeV	0.1–20	0.2–20	0.3–20

Table 3. Table 3: Fitting parameters of Equation ( 13 ) in use within Figure II-2 :

Fitting Functions	${𝑻}_{\mathrm{𝟎}}$	${𝑻}_{\mathrm{𝟏}}$	${𝒒}_{\mathrm{𝟎}}$	${𝒒}_{\mathrm{𝟏}}$
$f_{BG}$	0.2515 $\pm$ 0.0005	0.1335 $\pm$ 0.0002	-	-
$f_1$	0.1343 $\pm$ 0.0003	$-$0.0041 $\pm$ 0.0001	1.135 $\pm$ 0.002	0.009 $\pm$ 0.001
$f_2$	0.1343 $\pm$ 0.0003	$-$0.0041 $\pm$ 0.0001	1.135 $\pm$ 0.002	0.009 $\pm$ 0.001
$f_3$	0.1190 $\pm$ 0.0002	$-$0.0412 $\pm$ 0.0002	1.129 $\pm$ 0.001	0.004 $\pm$ 0.001
$f_4$	0.1083 $\pm$ 0.0003	0.0011 $\pm$ 0.0004	1.135 $\pm$ 0.002	0.009 $\pm$ 0.001
$f_5$	0.1222 $\pm$ 0.0005	0.0007 $\pm$ 0.0001	1.127 $\pm$ 0.002	0.013 $\pm$ 0.002

Table 4. Table 4: Fitting p T subscript 𝑝 𝑇 p_{T} range of different hadron spectra in heavy-ion collisions in this work  ALICE ( 2014 , 2016 , 2013 , 2014 , 2015 ) .

Particles	Mass [GeV]	$𝒑𝑷𝒃$ [GeV]	$𝑷𝒃𝑷𝒃$ [GeV]
$\pi$	0.140	0.11–2.85	0.11–19
$K$	0.494	0.225–2.45	0.225–19
$K_S^0$	0.498	0.05–7	0.45–11
$K^{\ast }$	0.896		0.55–4.5
$p$	0.938	0.325–3.9	0.325–17.5
$\varphi$	1.019		0.65–4.5
$\Lambda$	1.116	0.65–7	0.65–11
$\Xi$	1.321		0.7–7.5
$\Omega$	1.672		1.3–7.5

Table 5. Table 5: Fitting parameters of Equation ( 14 ) between T 𝑇 T and q − 1 = 1 / n 𝑞 1 1 𝑛 q-1=1/n for spectra of π 𝜋 \pi (upper) and K 𝐾 K (lower) in p ​ p 𝑝 𝑝 pp , p ​ P ​ b 𝑝 𝑃 𝑏 pPb , and  P ​ b ​ P ​ b 𝑃 𝑏 𝑃 𝑏 PbPb collisions (note that f B ​ G subscript 𝑓 𝐵 𝐺 f_{BG} is not included because q = 1 𝑞 1 q=1 is a constant).

Particles	Fittings	${𝑻}_{\mathrm{𝟏}}$ in $𝒑𝒑$	${𝑻}_{\mathrm{𝟎}}$ in $𝒑𝒑$	${𝑻}_{\mathrm{𝟏}}$ in $𝒑𝑷𝒃$	${𝑻}_{\mathrm{𝟎}}$ in $𝒑𝑷𝒃$	${𝑻}_{\mathrm{𝟏}}$ in $𝑷𝒃𝑷𝒃$	${𝑻}_{\mathrm{𝟎}}$ in $𝑷𝒃𝑷𝒃$
	$f_1$	$-$0.36 $\pm$ 0.02	0.08 $\pm$ 0.01	-	-	1.40 $\pm$ 0.02	0.335 $\pm$ 0.004
	$f_2$	$-$0.36 $\pm$ 0.02	0.08 $\pm$ 0.01	-	-	1.40 $\pm$ 0.02	0.335 $\pm$ 0.004
$\pi$	$f_3$	$-$0.14 $\pm$ 0.04	0.07 $\pm$ 0.02	-	-	1.22 $\pm$ 0.02	0.249 $\pm$ 0.005
	$f_4$	$-$0.22 $\pm$ 0.01	0.08 $\pm$ 0.01	-	-	1.52 $\pm$ 0.02	0.333 $\pm$ 0.005
	$f_5$	$-$0.31 $\pm$ 0.03	0.08 $\pm$ 0.01	-	-	1.31 $\pm$ 0.01	0.311 $\pm$ 0.007
	$f_1$	$-$1.30 $\pm$ 0.02	0.011 $\pm$ 0.001	1.55 $\pm$ 0.02	0.470 $\pm$ 0.001	1.67 $\pm$ 0.06	0.434 $\pm$ 0.003
	$f_2$	$-$1.30 $\pm$ 0.02	0.011 $\pm$ 0.001	1.55 $\pm$ 0.02	0.470 $\pm$ 0.001	1.67 $\pm$ 0.06	0.434 $\pm$ 0.003
$K$	$f_3$	$-$0.90 $\pm$ 0.04	0.032 $\pm$ 0.005	1.94 $\pm$ 0.03	0.436 $\pm$ 0.004	1.96 $\pm$ 0.03	0.394 $\pm$ 0.007
	$f_4$	$-$0.81 $\pm$ 0.01	0.010 $\pm$ 0.004	2.03 $\pm$ 0.03	0.470 $\pm$ 0.003	2.11 $\pm$ 0.05	0.427 $\pm$ 0.006
	$f_5$	$-$1.59 $\pm$ 0.02	$-$0.001 $\pm$ 0.0005	2.43 $\pm$ 0.01	0.453 $\pm$ 0.002	0.73 $\pm$ 0.02	0.309 $\pm$ 0.007

Equations

$${\text{e}}_q(x):=[1+(1-q)x{]}^{\frac{1}{1-q}}.$$

$$E\frac{{\text{d}}^{3}N}{{\text{d}}^{3}p}=\frac{{\text{d}}^{3}N}{\text{d}y{p}_T\text{d}{p}_T\text{d}\varphi }=\frac{1}{2\pi p_T}\frac{{\text{d}}^{2}N}{\text{d}y\text{d}{p}_T}.$$

$$f_0$$

$$f_1$$

$$f_2$$

$$f_3$$

$$f_4$$

$$f_5$$

$$A_1=A_2\cdot \frac{(n_2-1)(n_2-2)}{2\pi n_2{T}_2[n_2{T}_2+m(n_2-2)]}=A_2\cdot C_q,\mathrm{and}{n}_1=n_2.$$

$$f_1=A_1\cdot {\left(1-\frac{m}{n_1{T}_1}\right)}^{-n_1}\cdot {\left(1+\frac{m_T}{n_1{T}_1-m}\right)}^{-n_1}.$$

$$A_1\cdot {\left(1-\frac{m}{n_1{T}_1}\right)}^{-n_1}=A_4,n_1=n_4,\mathrm{and}{n}_1{T}_1-m=n_4{T}_4.$$

$$\displaystyle\textrm{\,e}_{q}(-\frac{p_{T}}{T})\longrightarrow\left\{\begin{array}[]{l}~{}\textrm{\,e}^{-p_{T}/T}\qquad~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}p_{T}\to 0\\ \\ \Big{(}(q-1)\frac{p_{T}}{T}\Big{)}^{\frac{1}{1-q}}\qquad p_{T}\to\infty.\end{array}\right.~{}.$$

$$\frac{1}{2\pi p_T}\frac{{\text{d}}^{2}N}{\text{d}y\text{d}{p}_T}=A_2\cdot C_q\cdot {\left(1+\frac{m_T-m}{n_2{T}_2}\right)}^{-n_2}$$

$${\int }_0^{p_{Tmax}}\frac{1}{2\pi p_T}\frac{{\text{d}}^{2}N}{\text{d}y\text{d}{p}_T}2\pi p_T\text{d}{p}_T=\frac{\text{d}N}{\text{d}y}.$$

$$\displaystyle T=T_{0}+T_{1}\cdot\Big{(}\frac{\sqrt{s}}{m}\Big{)}~{},\quad{\rm and}\quad~{}q=q_{0}+q_{1}\cdot\ln\Big{(}\frac{\sqrt{s}}{m}\Big{)}~{}.$$

$$T\approx T_0-(q-1)T_1,$$

Full text

I Introduction

In high-energy nuclear physics, the investigation of transverse momentum ($p_{T}$) spectra is a fundamental measure in statistical approaches. The $p_{T}$ spectrum reveals information on the kinetic properties of the particles produced in high-energy collisions. Strong correlation phenomena were recently observed in proton–proton and heavy-ion collisions ALICE (2017, 2017), their statistical and thermodynamical description points beyond the classical Boltzmann–Gibbs (BG) statistics. It has long been realized that data on single inclusive particle distributions show a power-law behavior in the high-$p_{T}$ region. For these, the Pareto–Hagedorn–Tsallis distribution has been frequently applied Pareto (1896); Hagedorn (1983); Tsallis (1988). Its form coincides with the generalized $q$-exponential function Tsallis (2009):

[TABLE]

Hadron spectra can be described by the Lorentz-invariant particle spectra. These were successfully fitted by the non-extensive distributions in a wide center-of-mass energy and $p_{T}$ range Osada (2008); STAR (2006); BRAHMS (2005); CMS (2010, 2010); PHENIX (2004); STAR (2006, 2005); Wilk (2000); ALICE (2011, 2014); PHENIX (2011); ALICE (2014); Tang (2009). In the following, we focus on the most often used formulas from Osada (2008); STAR (2006); BRAHMS (2005); CMS (2010, 2010); PHENIX (2004); STAR (2006, 2005); Wilk (2000) for representing identified particle spectra in various collisions. This work explores differences between ($m_{T}-m$) and $m_{T}$-dependent, as well as simple $p_{T}$ functions:

[TABLE]

Different research groups used various kinds of expressions of it in order to describe $p_{T}$ spectra. We consider functions of $m_{T}-m$ and $p_{T}$ in the non-extensive approach, after applying the normalized functions and the thermodynamically motivated ones Shen (2019). Our aim is to find the best-fitting functions among these, while assigning a physical interpretation to their parameters. We investigate the following distribution forms:

[TABLE]

There are relations among the distributions defined above. It is easy to realize that $f_{1}$ and $f_{2}$ coincide whenever their amplitudes satisfy the relation

[TABLE]

Accounting for the differences between ($m_{T}-m$) and $m_{T}$ dependencies, we re-cast $f_{1}$ and $f_{4}$ described in Equation (3) as follows:

[TABLE]

Comparing this with $f_{4}$, we arrive at the relations

[TABLE]

These comments are important for the comparison of different approaches. They also demonstrate that no inconsistency occurs by applying different fit formulas. However, differences arise from the statistical physical motivations behind these formulas Osada (2008); STAR (2006); BRAHMS (2005); CMS (2010, 2010); PHENIX (2004); Shen (2019, 2018). The corresponding results and discussions are investigated next. Note that for all the physical quantities, we use the natural units, $c=1$, for convenience in this paper.

II Results and Discussions

In this section, we analyze the transverse momentum distributions of identified pions and kaons stemming from the elementary ($pp$) and heavy-ion ($pPb$ and $PbPb$) collisions fitted by the functions listed in Equation (3). All the relevant parameters are then analyzed in order to investigate further the non-extensive physics behind these collisions.

II.1 Analysis of the $pp$ Spectra

In high-energy physics, even the smallest hadron–hadron ($pp$) collisions are rather complicated processes. One usually separates two main regimes of hadron production: one is a soft multiparticle production, dominant at low transverse momenta, where the spectra can also be fitted by an exponential behavior Hagedorn (1995), cf. the curve $f_{BG}$ in Figure II-1. We realize that $f_{BG}$ describes well this part of the spectra even in $pp$ collisions. As $p_{T}$ gets higher ($p_{T}>$3 GeV), the spectrum displays a power-law tail. They are predicted by perturbative QCD, owing to the hard scattering of current quarks and gluons. In a number of publications ALICE (2011, 2014); PHENIX (2011); ALICE (2014); Tang (2009), the Tsallis statistical distribution was successfully applied to describe data for $pp$ collisions over a wide range of the transverse momenta because of its two limits: the exponential shape at small $p_{T}$ and the power-like distribution at large $p_{T}$,

[TABLE]

We focus on the fittings of the produced charged particle spectra in elementary collisions with the non-extensive functions in Equation (3). Data were taken for pions, kaons, and protons in $pp$ collisions at $\sqrt{s}=62.4$ GeV, $200$ GeV from the PHENIX Collaboration PHENIX (2011) and at 900 GeV ALICE (2011), 2.76 TeV ALICE (2014), 5.02 TeV, and 7 TeV ALICE (2014) from the ALICE Collaboration. We restrict our analysis to the midrapidity region $|y|<0.5$ within the $p_{T}$ ranges, as shown in Table 2. Note that in the following, $\pi$, $K$, and $p$ mark the spectra of $\frac{\pi^{+}+\pi^{-}}{2}$, $\frac{K^{+}+K^{-}}{2}$, and $\frac{p+\bar{p}}{2}$, respectively.

Figure II-1 shows that all of the different non-extensive functions we used fit the pion and kaon spectra very well for various kinds of beam energies at midrapidity. The ratios of $\chi^{2}/d.o.f.$ of the relevant fits are given in Table 1. Specifically, the first two distributions ($f_{1}$ and $f_{2}$) of $m_{T}-m$ and $f_{4}$ of $m_{T}$ show close-fitting results. The distribution, $f_{3}$, derived thermodynamically, does not display large differences in the goodness of fit either. Checking the fitting parameters $A$, $T$, and $q=1+1/n$, we observe that, as we expected and introduced in the previous section, all these functions share the same Tsallis parameter $n$. The two $m_{T}-m$ functions ($f_{1}$ and $f_{2}$) lead to fitting values of the temperature $T$, which are different from the pure $m_{T}$ fit ($f_{4}$). This indicates that the normalization constant does not affect the fitted $T$ and $q$ parameters but the integrated yield $\textrm{\,d}N/\textrm{\,d}y$. Namely, by normalizing the momentum spectrum

[TABLE]

with the $C_{q}$ normalization constant and the condition of $A_{2}=\textrm{\,d}N/\textrm{\,d}y$, we obtain the integral over $p_{T}$ from 0 to its maximal values $p_{Tmax}$:

[TABLE]

Moving towards physical interpretation issues, we investigate the temperature, $T$, and the non-extensive parameter, $q$. Investigations in PHENIX (2011); Cleymans (2013) showed that both of them express $\sqrt{s}$ dependence. In this paper, we found that they are also dependent on the hadron mass, $m$. The $\sqrt{s}/m$ dependence, as a result, is studied in order to analyze hadron spectra parameters within the non-extensive approach. Following the phenomenological observations in Barnafoldi (2011); Takacs (2019), a QCD-like evolution can be introduced for both the parameters $T$ and $q$. While analyzing data, we found that the temperature $T$ had a weak logarithmic $\sqrt{s}/m$ dependence. Thus, here we assume a linear $\sqrt{s}/m$ dependence to analyze the temperature $T$, but the non-extensive parameter $q$ is kept with the stronger logarithmic distribution:

[TABLE]

In summary, our work indicates that the BG distribution is not suitable for describing the hadron spectra over a wide range of $p_{T}$. Comparisons of their corresponding fitting errors $\chi^{2}/d.o.f.$ show that both $m_{T}-m$ and $m_{T}$ functions share the same goodness between $f_{1}$ and $f_{2}$, cf. Equation (3). Together with the thermodynamically derived $f_{3}$, all the non-extensive approaches ($f_{1}\sim f_{4}$) follow the experimental data accurately. The fitting temperature, $T$, is nearly constant when changing the ratio of the collision energy to hadron mass, $\sqrt{s}/m$. Specifically, distributions of $f_{1}$, $f_{2}$, $f_{4}$, and $f_{5}$ are described best with such a connection, as shown in the left panel of Figure II-2. From Table 3, we also see that the slope parameters in these four cases are almost zero, which means that they are constant around some values. The non-extensive parameter $q$, on the other hand, follows a logarithmic dependence, agreeing with a pQCD-based motivation, cf. Shen (2019). Note that our results on $T$ and $q$ are different from the work by Cleymans et al. Cleymans (2013). Those authors parameterized this relation as a power-law.

II.2 Analysis of the $pPb$ and $PbPb$ Results

In $pPb$ ALICE (2014) collisions at 5.02 TeV and in $PbPb$ ALICE (2016, 2013, 2014, 2015) collisions at 2.76 TeV, more kinds of hadron spectra are analyzed within the formulas of Equation (3). Data are taken from the ALICE Collaboration within wide $p_{T}$ ranges, as seen in Table 4. We observe that all of them present good fittings over the whole range of $p_{T}$ for each hadron at various kinds of centrality bins. On the other hand, similar to the $pp$ cases, the BG formula can still perform well just in the low $p_{T}$ region ($p_{T}<3$ GeV).

In this work, as an example, we analyzed the fitting results of $p_{T}$ spectra of pions and kaons produced in all kinds of collisions mentioned above. It is instructive to plot the relationship between the fitting temperature $T$ and the Tsallis parameter $q$ for the same hadron spectra for different centralities in the same heavy-ion collisions. The results of pions and kaons in $pp$ collisions are also analyzed as comparisons. In Figure II-3, we show the linear correlating appearances for both $\pi$ and $K$ in $pPb$ at 2.76 TeV ALICE (2014) and in $PbPb$ at 5.02 TeV ALICE (2016, 2013) as well as the $pp$ results in all kinds of collision energies ALICE (2011, 2014); PHENIX (2011); ALICE (2014) in this paper. In fact, whatever kinds of particle we study, all these non-extensive fittings give a similar dependence of $T$ on the parameter $q$:

[TABLE]

which agrees with our previous work Shen (2018, 2019) and that of others Wilk (2015).

Note that the slope parameter $T_{1}$ in Table 5 turns negative and $T_{0}$ is nearly zero for the $pp$ case, as discussed in Shen (2018). Results of fittings on pion spectra, typically in $pPb$ collisions at 5.02 TeV, fail in the obvious linear combinations probably due to the small mass of pions and high multiplicities. It is found that all forms of non-extensive distributions feature a similar relation between the temperature $T$ and non-extensive parameter $q$. This, in turn, hopefully promotes a better understanding of the meaning of the non-extensive parameter $q$.

III Summary

In this work, we analyzed various fitting formulas of the hadron spectra in order to explore their sensitivity to different fitting parameters in use within the non-extensive approaches, cf. Equation (3). The hadronization, as well as the $p_{T}$ distributions in high-energy physics (in proton–proton, proton–nucleus, and nucleus–nucleus collisions) are being studied here. For more details, see Shen (2019).

Our results reveal that normalization parameters have no major effect on the shape of these functions. In other words, the fitting formulas of either $m_{T}-m$ or $m_{T}$ lead to the same fit quality. As shown in Table 1, they obtained similar fitting values of $\chi^{2}/d.o.f.$ Finally, we investigated the relationship between the fitting parameters, $T$ and $q$. In $pp$ collisions, the temperature values were fitted by the linear relation of $\sqrt{s}/m$, while the non-extensive parameter $q$ had a logarithmic $\sqrt{s}/m$ dependence, motivated by the QCD-like evolution Barnafoldi (2011); Takacs (2019). All kinds of approaches led to linear relations between the temperature, $T$, and the non-extensive parameter, $q-1$, in heavy-ion collisions at different centralities. This agrees well with our previous results Shen (2018, 2019) and others in Wilk (2015).

Summarizing, based on the Tsallis $q$-exponential, five types of non-extensive formulas in Equation (3) were investigated in parallel to the usual BG distribution. Results showed that the BG statistics failed in describing the hadronization in the whole $p_{T}$ range. Within the non-extensive approaches, $m_{T}-m$ functions obtained similar fitting results to the $m_{T}$ ones. This provides a free choice between the functions $m_{T}-m$ and $m_{T}$ when analyzing the hadron spectra. On the other hand, it does not make any differences with regards to the normalization. Nevertheless, the normalized function, $f_{2}$, is the best choice since it is also connected to the particle yield per unit rapidity, $\textrm{\,d}N/\textrm{\,d}y$, by its normalization, $A_{2}$.

Acknowledgments

This work has been supported by the Hungarian National Research, Development and Innovation Office (NKFIH) under the contract numbers K120660 and K123815 and THOR COST CA 15213.