Multiplicity Dependence in the Non-Extensive Hadronization Model Calculated by the HIJING++ Framework

TL;DR

This paper explores how the parameters of the non-extensive hadronization model depend on multiplicity in high energy collisions, using the HIJING++ framework to compare with experimental data and predict high-multiplicity behavior.

Contribution

It provides a status update on tuning the HIJING++ generator and investigates multiplicity dependence of model parameters with new calculations.

Findings

HIJING++ shows good agreement with existing data

Parameters exhibit significant multiplicity dependence

Predictions made for high-multiplicity regions with limited data

Abstract

The non-extensive statistical description of the identified final state particles measured in high energy collisions is well-known by its wide range of applicability. However, there are many open questions that need to be answered, including but not limited to, the question of the observed mass scaling of massive hadrons or the size and multiplicity dependence of the model parameters. This latter is especially relevant, since currently the amount of available experimental data with high multiplicity at small systems is very limited. This contribution has two main goals: On the one hand we provide a status report of the ongoing tuning of the soon-to-be-released HIJING++ Monte Carlo event generator. On the other hand, the role of multiplicity dependence of the parameters in the non-extensive hadronization model is investigated with HIJING++ calculations. We present cross-check comparisons of HIJING++ with existing experimental data to verify its validity in our range of interest as well as calculations at high-multiplicity regions where we have insufficient experimental data.

Tables

Table 1. Table 1: Main internal parameters in HIJING++ .

Parameter	Description
$p_0$	soft-hard separation scale: minimum $p_T$ transfer of hard or semihard scatterings
${\sigma }_{soft}$	the inclusive cross section for soft interactions
${\sigma }_0$	the cross section that characterizes the geometrical size of a nucleon
${\mu }_0$	the parameter in the scaled eikonal function of nucleon used to calculate total cross-section
$K$	K-factor for the differential jet cross sections in the lowest order pQCD calculation
$\max p_T^{}{}_{cut}{}^{}$	$p_T$ cut for classifying the connected-independent type strings at fragmentation
$m_{inv-cut}$	invariant mass cut-off for the dipole radiation of a string system below which soft gluon radiations are terminated
$m_{min-inv-ex.str.}$	minimum value for the invariant mass of the excited string system in a hadron–hadron interaction
$S_{p_{T_1}}$	the parameter that regularizes the singularity at $p_T=0$ in the distribution of the soft $p_T$ kick
$S_{p_{T_2}}$	the parameter that gives the scale beyond which the $p_T$ kick distribution will be similar to $1/p_T^4$
$F$	the scale in the form factor to suppress the $p_T$ transfer to diquarks in hard scatterings
$v_{q_i}$	phenomenological parameters ($i=1,2,3$) of the soft parton distribution function that yield an $x$ distribution of the valence quarks in a soft interaction
$s_{q_i}$	phenomenological parameters ($i=1,2,3$) of the soft parton distribution function that yield an $x$ distribution of the sea quarks in a soft interaction
StringPT:temperature	the temperature parameter in the Lund fragmentation model as described in PYTHIA
StringPT:tempPreFactor	the temperature prefactor for strange quarks and diquarks in the Lund fragmentation model as described in PYTHIA
StringZ:aExtraSQuark	parameters in the Lund symmetric fragmentation function as described in PYTHIA
StringZ:aExtraDiuark

Table 2. Table 2: Tsallis parameters extracted from “minimum bias” I ​ N ​ E ​ L > 0 𝐼 𝑁 𝐸 𝐿 0 INEL>0 proton–proton collisions at s = 7 𝑠 7 \sqrt{s}=7 TeV, measured by ALICE  ALICE:PID

Hadron	$𝒏$	$𝒒$	$𝑻$ (GeV)	$𝑨$	$𝒗$	${𝝌}^{\mathrm{𝟐}}/𝒏𝒅𝒇$
${\pi }^{\pm }$	7.415 $\pm$ 0.033	1.135 $\pm$ 0.005	0.089 $\pm$ 0.010	73.188 $\pm$ 9.700	0.000 $\pm$ 0.119	174.225 / 54
$K^{\pm }$	7.539 $\pm$ 0.086	1.133 $\pm$ 0.013	0.155 $\pm$ 0.010	0.915 $\pm$ 0.095	0.000 $\pm$ 0.066	20.274 / 47
$p\left(\overline{p}\right)$	8.805 $\pm$ 0.184	1.114 $\pm$ 0.023	0.191 $\pm$ 0.012	0.124 $\pm$ 0.013	0.000 $\pm$ 0.054	18.462 / 45

Table 3. Table 3: Multiplicity classes used in HIJING++ runs.

Class	I	II	III	IV	V	VI	VII	VIII	IX	X	XI	XII	XIII	XIV	XV
	0	10	15	20	25	30	35	40	45	50	55	60	70	80	90
	10	15	20	25	30	35	40	45	50	55	60	70	80	90	100

Table 4. Table 4: The multiplicity averaged temperature parameters for charge averaged pions, kaons, and protons.

Hadron	(GeV)
${\pi }^{\pm }$	$0.063\pm 0.003$
$K^{\pm }$	$0.092\pm 0.001$
$p\left(\overline{p}\right)$	$0.106\pm 0.002$

Equations

$$\left.\frac{d^{2}N}{p_{T}dp_{T}dy}\right|_{INEL>0}=A\cdot m_{T}\cdot\left(1+\frac{E}{nT}\right)^{-n}\ \ \$$

$$\text{Class}={⟨d{N}_{ch}/d\eta ⟩}_{min}<⟨d{N}_{ch}/d\eta ⟩\le {⟨d{N}_{ch}/d\eta ⟩}_{max}.$$

Full text

\pubvolume

xx \issuenum1 \articlenumber5

\historyReceived: 30 April 2019; Accepted: 23 May 2019; Published: date

\TitleMultiplicity Dependence in the Non-Extensive Hadronization Model Calculated by the HIJING++ Framework \AuthorGábor Bíró 1,2,∗,†\orcidA, Gergely Gábor Barnaföldi 1,∗,†\orcidC, Gábor Papp 2,∗,†\orcidB and Tamás Sándor Biró 1,∗,†\orcidD

\AuthorNamesGábor Bíró, Gergely Gábor Barnaföldi, Gábor Papp, Tamás Sándor Biró

\corresCorrespondence: biro.gabor@wigner.mta.hu (G.B.); barnafoldi.gergely@wigner.mta.hu (G.G.B); pg@ludens.elte.hu (G.P.); biro.tamas@wigner.mta.hu (T.S.B.)

\firstnoteThese authors contributed equally to this work.

1 Introduction

The transverse momentum ($p_{T}$) distribution of identified hadrons stemming from high-energy proton–proton, proton–nucleus, and nucleus–nucleus collisions is one of the most fundamental observables in high-energy physics. In recent years, the Tsallis–Pareto-like distributions, motivated from non-extensive statistical physics, have received close attention because their applicability in this field BIROG:TSALLIS2 ; GRIGORYAN:TSALLIS ; ZHU:TSALLIS ; CLEYMANS:TSALLIS ; WONG:TSALLIS ; TRIPATHY:TSALLIS ; SHEN:TSALLIS2 . With the appearance of high precision experimental data spanning from low- to high-$p_{T}$, neither the thermal models with a bare Boltzmann–Gibbs exponential distribution nor perturbative Quantum Chromodynamics (pQCD)-motivated power-law distributions are able to describe the whole spectrum. On the other hand, the Tsallis–Pareto distributions combine these two regions perfectly (see, e.g., BIROG:TSALLIS1 ; BIROG:TSALLIS2 ; BIROG:TSALLIS3 ; BIROG:TSALLIS4 ; TAKACS:TSALLIS ; GRIGORYAN:TSALLIS ; URMOSSY:TSALLIS ; ZHU:TSALLIS ; VAN:TSALLIS ; KHUNTIA:TSALLIS ; WILK:TSALLIS ; CLEYMANS:TSALLIS ; SHEN:TSALLIS ; WONG:TSALLIS ; TRIPATHY:TSALLIS ; SHEN:TSALLIS2 and references therein). During the investigation of the parameters, we showed that they possess non-trivial relations such as mass- and energy scaling BIROG:TSALLIS1 ; BIROG:TSALLIS2 . There are also implications that for larger systems a soft-hard extension is needed BIROG:TSALLIS4 ; SHEN:TSALLIS . These studies indicate that increasing the size of the colliding system (roughly speaking, the volume of the quark-gluon plasma) may also reflect in the parameters. Our goal is therefore to systematically explore the parameter space as the function of the event multiplicity.

The HIJING++ framework is a soon-to-be-published general purpose Monte Carlo event generator, currently in the final phase of development HIJING ; HIJING2 ; HPP1 ; HPP2 ; HPP3 ; HPP4 . It will serve as the successor of the FORTRAN HIJING, completely rewritten in modern C++. With the flexibility gained by using modular C++ structures, HIJING++ also utilizes several external packages PYTHIA ; LHAPDF ; GSL ; VEGAS ; ROOT ; PDF1 ; PDF2 . Currently the internal parameters of HIJING++ are being tuned to main experimental observables using Professor PROFESSOR ; PARALLEL . This provides an excellent opportunity to test the capabilities of HIJING++ and calculate high-multiplicity events.

In the next section, we briefly summarize the progress of tuning in HIJING++ and present the current status. In Section 3, a theoretical description of the transverse momentum spectra is given, and the HIJING++ calculations are given in Section 4.

2 Tuning of HIJING++ Parameters

A typical general purpose Monte Carlo event generator, e.g., the HIJING++ framework, developed to be able to simulate high-energy heavy-ion collisions, has parameters that are not determined by theory and need to be tuned to reproduce measured experimental data with the highest possible precision. One of the main features of the HIJING++ framework is that very few input parameters are needed to fully define a run, such as the species of the projectile and target beam and center-of-mass energy. Given this information, all of the other intrinsic parameters are calculated automatically.

Since HIJING++ is based on the convolution of sequential collisions of nucleon–nucleon pairs in each nucleus–nucleus interaction, it is highly important to have a solid proton–proton collisions baseline. In this section, we present the up-to-date result of the tuning process using the following $\sqrt{s}=7$ TeV proton–proton experimental data:

• •

$p_{T}$ spectra of identified $\pi^{\pm}$, $K^{\pm}$, and $p(\bar{p})$ hadrons with $INEL>0$ normalization (at least one charged particle in the $|\eta|<1.0$ region is required) up to $p_{T}=20$ GeV/c ALICE:PID ;

• •

charged hadron multiplicity distribution in the range of $\left<dN_{ch}/d\eta\right>=0-70$, where $N_{ch}$ is the number of charged particles ALICE:MULTPICITY1 ; ALICE:MULTPICITY2 ;

• •

charged hadron $\eta=\frac{1}{2}\ln{\frac{p+p_{Z}}{p-p_{Z}}}$ pseudorapidity distribution at mid-pseudorapidity $|\eta|<1.0$ ALICE:MULTPICITY1 .

The tuning process is performed iteratively utilizing the Professor tool PROFESSOR ; PARALLEL . In Table 1, we list the main tunable parameters.

The detailed process of tuning and the parameters values will be described in the technical release paper of HIJING++. Here we present only the tuning status regarding the experimental data listed above.

In Figure 1, the multiplicity and the pseudorapidity distribution of charged hadrons are presented, while in Figure 2 the $p_{T}$ spectrum of identified $\pi^{\pm}$, $K^{\pm}$, and $p(\bar{p})$ hadrons calculated and measured at $\sqrt{s}=7$ TeV proton–proton collisions can be seen.

The results above show that HIJING++ reproduces the event multiplicity excellently. In Figure 1, the agreement between the HIJING++ results and the experimental data is $\sim$15% for the multiplicity and $\sim$1% for the pseudorapidity distribution. The charged pion and kaon spectra also show a good agreement above $p_{T}=2$ GeV/c, but the production is slightly overestimated at lower $p_{T}$ values. The best agreement for the $\pi^{\pm}$ results is $\sim$1% between 2 and 15 GeV/c. For kaons, the yield is slightly underestimated above 2 GeV/c, where the agreement is $\sim$15–20%. On the other hand, the proton yield is overestimated in the large $p_{T}$ region, the agreement is $\sim$20–30%.

3 The Non-Extensive Hadronization Model

It is a well-known and an intensively studied phenomenon that the transverse momentum spectra of hadrons stemming from high-energy particle collisions can be described by Tsallis–Pareto type distributions BIROG:TSALLIS1 ; BIROG:TSALLIS2 ; BIROG:TSALLIS3 ; BIROG:TSALLIS4 ; TAKACS:TSALLIS ; GRIGORYAN:TSALLIS ; ZHU:TSALLIS ; KHUNTIA:TSALLIS ; WILK:TSALLIS ; CLEYMANS:TSALLIS . Although this observation itself has further consequences, the theory has even more subtle details because of the observed non-trivial dependence on the center-of-mass energy and hadron mass. In the following sections, we show that the parameters also depend on the event multiplicity, i.e., on the size of the system.

We adopted the usual blast-wave assumptions regarding the system, namely that the fireball is azimuthally symmetric and is expanding with a $v$ radial flow velocity (in units of $c=1$). Moreover, the freeze-out occurs instantly on a hypersurface according to the Cooper – Frye formulation at a given freeze-out temperature GRIGORYAN:TSALLIS ; URMOSSY:TSALLIS . With these assumptions, we used the following simple form of the invariant yield:

[TABLE]

where $A$ is the amplitude incorporating the irrelevant spin degeneracy and constant factors as well as the invariant volume, $m_{T}=\sqrt{p_{T}^{2}+m^{2}}$ is the transverse mass, $E=\gamma(m_{T}-vp_{T})-m$ is the one-particle energy in the co-moving coordinate system, $\gamma=1/\sqrt{1-v^{2}}$ is the Lorentz factor, $T$ is a parameter with a temperature unit, and finally $n=\frac{1}{q-1}$ is the non-extensivity parameter, characterizing the temperature fluctuations. We note that $T$ is not necessarily the freeze-out temperature and therefore is not necessarily the same for all hadron species ZHU:TSALLIS ; VAN:TSALLIS . The notation $INEL>0$ means that only those events where there is at least one charged particle in the $|\eta|<1.0$ region are considered. This choice is in agreement with the experimental definitions described in the previous section.

As a reference, in Table 2 and in Figure 3, we show the parameters and curves fitted on the experimental “minimum bias” (in the sense that there is no event multiplicity classification) data. These results are consistent with our previous observations BIROG:TSALLIS1 : the heaviest proton has the largest temperature and the smallest $q$. We note that, in the lower part of Figure 3, a periodic oscillation is visible. This is an effect in addition to the scaling. This has been investigated, for example, in WILK:OSCILLATIONS ; RYBCZYNSKI:OSCILLATIONS .

4 The Multiplicity Dependence of the Non-Extensive Model

In Section 2, we showed that the HIJING++ framework is able to reproduce the main experimental observables such as event multiplicity distribution and the $p_{T}$ spectra of various identified hadrons. In Section 3, we briefly summarized the main features of the blast-wave motivated non-extensive hadronization model. In this section, we take advantage of the power of HIJING++ and extract the Tsallis parameters from a wide range of event multiplicity classes. The event classes of the HIJING++ run are classified as

[TABLE]

The multiplicity ranges of each class used in this study are listed in Table 3.

Using this event classification, we calculated the mid-rapidity transverse momentum spectra of charge averaged pions, kaons, and protons in $INEL>0$ events. We generated $200$M events. To avoid superfluous overcrowding of the available space, we show only the low, moderate, and high multiplicity spectra along with the fitted Tsallis–Pareto curves defined by Equation (1) in Figure 4.

In Figure 4, we can see that the best fits occurred at the high multiplicity events. For the pions, the fitted curves follow the points well at the low-$p_{T}$ region, while for the kaons and protons with higher mass at the low-$p_{T}$ region the model overpredicts the yield. We note, however, that for the HIJING++ run we investigated the same $0.1$ GeV/c$<p_{T}<20$ GeV/c region for all hadrons, while the low-$p_{T}$ part for the kaons and protons in the case of the experimental results is missing. At the high-$p_{T}$ region, the fit breaks down because of the low statistics.

In Figure 5, the fitted parameters in the function of the event multiplicity are shown. Using the distribution form Equation (1), we observe that with increasing multiplicity the $q$ parameter increases for each hadron but with different slopes: the increase of $q$ (or the decrease of $n$) is the largest for the heaviest hadron. On the other hand, the temperature decreases slowly with the increasing pseudorapidity density. Here, the previously observed $T_{\pi^{\pm}}<T_{K^{\pm}}<T_{p(\bar{p})}$ mass hierarchy stays valid with the multiplicity averaged values, as can be seen in Table 4.

The radial flow velocity also increases with the multiplicity, but also with different rates. While at low multiplicity the lightest pions have the smallest $v$, it increases rapidly with the increasing multiplicity. On the other hand, the rate of increase in the case of protons and kaons are approximately the same. These observations require further investigation. Finally, the amplitudes are increasing for each hadron species with the multiplicity. The value of pions is much higher than those of the heavier hadrons, which indicates that with increasing multiplicity the number of the produced pions grows faster than the number of kaons and protons.

5 Conclusions

In this contribution, we investigated the multiplicity dependence of the parameters of the non-extensive hadronization model in proton–proton collisions using HIJING++ calculations. We presented the current status of the tuning process of HIJING++ and showed that it is able to reproduce the main high-energy physics observables such as multiplicity and $p_{T}$ distributions. We presented the non-extensive hadronization model that we used to describe the transverse momentum distribution of identified hadrons. In accordance with our previous results, we showed that a mass hierarchy emerges in the Tsallis parameters. Utilizing the tuned HIJING++ calculations, we also extracted the parameters from $\sqrt{s}=7$ TeV proton–proton collisions Monte Carlo calculations with various event multiplicity classifications. Our study showed that the $q$ non-extensivity parameter increases with increasing multiplicity, while the $T$ temperature has only a slight decrease. On the other hand, all hadrons result in a non-zero, increasing radial flow velocity. All parameters show the earlier observed mass hierarchy. Our findings suggest that these parameters are sensitive to the event size and may serve as a thermometer of the collision.

\authorcontributions

Software and formal analysis: G.B. and G.P.; investigation: G.B. and G.G.B.; writing—original draft preparation: G.B.; writing—review and editing: G.G.B.; supervision: G.G.B., T.S.B., and G.P.; funding acquisition: G.G.B. and T.S.B.

\funding

This research was funded by Hungarian-Chinese cooperation grant No. MOST 2014DFG02050, the Wigner HAS-OBOR-CCNU grant, and OTKA grants K120660 and K123815. G.B. acknowledges the support of the Wigner Data Center and Wigner GPU Laboratory.

\conflictsofinterest

The authors declare no conflict of interest.

\reftitle

References