Non-extensive Motivated Parton Fragmentation Functions

A new form of fragmentation function is presented here, motivated by earlier non-extensive studies of jet fragmentation. We parametrized our Tsallis-like function on pion spectra and compared it to the most common fragmentation function parametrizations. It is shown that the new form is in agreement with earlier parametrizations, furthermore, its scale evolution overlap better with the experimental data.


Introduction
In high energy collisions the energy distribution of the outgoing particles with low momentum ( 3 GeV/c) follows a thermal-like, exponential distribution. Particles with higher momentum ( 5 GeV/c), where pQCD works better, follow power law distributions. The whole spectrum cannot be described by one exponential or power law function, but the Tsallis distribution fits it perfectly [1].
Statistical systems, like high-energy proton-proton (pp) or electron-positron (e − e + ) collisions are small. The created ∼ 100 particles are far from the thermodynamic limit, highly fluctuating and correlated. Therefore, they could be beyond the validity of statistical mechanics, which was built for weakly interacting systems [2][3][4]. The q entropy is a generalization of Boltzmann's one, which could lead to a better description of this fluctuating collision system [5]. The theory has a temperature-like parameter, T and an additional q parameter describing the non-additivity of the entropy for independent systems. Tsallis is the characteristic distribution, derived from this so-called non-extensive entropy, analogue to the Boltzmann distribution. In the q → 1 limit, Tsallis' entropy and distribution restores Boltzmann's.
In Ref. [6], we built a model for jet fragmentation resulting the so-called microcanonical Tsallis distribution, derived from the observed Gamma shape of the multiplicity distribution. This result is extended here, by applying the form as a fragmentation function (FF). Our paper has two aims: (i) to present the validity of our FF assumption by comparing it to widely applied FFs and experimental data. (ii) To investigate the non-extensive parameters of FFs by comparing them to earlier studies [1].

The Formalism
In the parton model assuming the factorization, perturbative (e.g. partonic cross sections) and non-perturbative (e.g. parton ditribution and hadronization) processes are treated separately. The normalized cross section of the outgoing h hadrons in e − e + annihilation is defined by in [7,8],

arXiv:1811.01974v1 [hep-ph] 5 Nov 2018
where Q is the momentum scale, z = 2E h / √ s is the momentum fraction, where hadron energy is E h and √ s is the center of mass energy. The σ tot (Q) is the total cross section. The inclusive differential cross section of hadron production is where C i functions are calculated from perturbative QCD. The non-perturbative D h i (z, Q) is the fragmentation function, characterizing the probability of an i parton fragmenting into a hadron, h at a given z and Q. The Q evolution of the FF is described by DGLAP evolution equations [7], where i, j runs for all (anti)quarks and the gluon. The P ij are the splitting functions from pQCD.
FFs are arbitrary functions satisfying Eq. (3) and The later is required for the momentum conservation and so the probability interpretation. There is no unified theory to determine FFs, thus phenomenological ansatz is used for the form and then it is parametrized with experimental data [8,9]. From Eq. (2) one can see that the shape of the hadron spectrum is the result of the convolution integral with the arbitrary shaped FF. Therefore, the spectra is not a pure power law despite the frequently referred statement about pQCD.
The most common FF ansatz is polynomial, with Q dependent parameters, Previous studies have determined the parameters of this FF and it uncertainties with experimental data [10][11][12]. This form above works perfectly above 5 GeV/c, however its parameters differ from theoretical expectations, based on the parton splittings [8].
In our previous studies, we determined the parameters of the Tsallis distribution in cases of various spectra for e − e + , pp, pA, AA collisions, different hadrons, and kinematic regimes [1]. We determined q ≈ 1.05 − 1.25 for pp at LHC energies, increasing with the collision energy. The temperature-like parameter T/Q ≈ 0.05 − 0.2, weakly depends on the collision energy and in contrast to q, it is independent from the number of valence quarks. We developed a non-extensive jet fragmentation function, the so-called microcanonical Tsallis distribution in Ref. [6] and we postulate the result as an FF ansatz, where N h i , q h i , T h i are Q dependent parameters, thus it has the same number of parameters as Eq.(4).

Results
To obtain the parameters of our FF, we defined Eq. (5) at a given Q 0 , then we evolved to Q using Eq. (3) to minimize the χ 2 red function, by changing the parameters of the FF, where y and δy is the measured spectrum value and its uncertainty, and N do f = 15 is the number of parameters. We restricted our parameters to be positive, but no further constrains were added, counter to the studies using Eq. (4). In this paper, we use e − e + → π ± data at √ s = 91.2 GeV, from Ref. [10]. We follow all notations and definitions from Ref. [10] in leading order (LO), e.g. independent channels or initial scales Q 0 = 1 GeV for u, d, s and g, Q 0 = m c , m b for c and b respectively. We summarize our first results from minimizing Eq (6). The left panel of Fig. 1 shows the measured spectrum and the fit we got. As a comparison, we plotted calculations with the HKNS [10] and KKP [11] parametrizations in LO using Eq. (4). The smallness of the relative differences highlight that our result is in a good agreement with the data and that there is almost no substantial difference between the various models. The right panel of Fig. 1 shows the probability densities, zD h i (z, Q) of the fitted microcanonical Tsallis FFs, compared to the HKNS with uncertainties and KKP results at a given Q. The different partonic channels agree mostly, however our formula does not diverge at z = 0 and its slopes are different. For simplicity, we did not use flavor tagged data in the fit, causing the differences seen in the c and b channels.  The microcanonical Tsallis FF parametrization in pion spectrum from Ref. [10]. HKNS and KKP labels two LO calculations using polynomial FF from Ref. [10,11]. Agreement between different models is clear. Right: The microcanonical Tsallis FF parametrization compared to HKNS with uncertainty and KKP FFs for different partonic channels at a given Q.  Fig. 2 shows FF parameters α π + i , β π + i from Eq. (4) for HKNS (left) and KKP (middle) parametrization, and q π + i , T π + i parameters from Eq. (5) microcanonical Tsallis FF (right) for different Q between √ s = 2 − 10 5 GeV. If α h i < −2, then dz zD h i → ∞ diverges, therefore FF looses the probability interpretation (gray area on Fig. 2). With Eq. (3), the energy evolution of the pion spectrum from Eq. (1) was determined for all three cases. We fitted these spectra with a single function using Eq. (4) for HKNS and KKP and Eq. (5) for microcanonical Tsallis. On the HKNS and KKP panels one can see that there is a weak correlation between the partonic parameters and parameters obtained from the spectrum fits. However, in the case of the microcanonical Tsallis, the correlation is stronger. Moreover, in all cases, the FF of u-quark agrees best with the spectra fit, which was expected because of their dominating contribution in pions. For more comparison, we fitted some further experimental data with the same way between √ s = 4.8 − 91.2 GeV from Ref. [10], which overlapped with the microcanonical Tsallis description.

Conclusion
In this paper we performed a non-extensive motivated fragmentation function parametrization and we calculated the LO spectra of pions in high-energy e − e + collision. We compared and verified our FF with widely used other FF parametrizations and experimental data as well. We obtained that our microcanonical Tsallis FF differs from the HKNS [10] and KKP [11] FFs at lower energy scales and momenta, which could affect the low momentum calculations of proton-proton collisions. The attractor behaviour of the DGLAP evolution could ensure the same behaviour at high Q [13], supporting the power law tailed pQCD observations. Our fragmentation suggests that there could be a more suitable description for the fragmentation, which overlaps better with the experimental data.