Investigation of Particle Distributions in Xe-Xe Collision at ffiffiffiffiffiffiffi s NN p = 5 : 44 TeV with the Tsallis Statistics

The distribution characteristic of final-state particles is one of the significant parts in high-energy nuclear collisions. The transverse momentum distribution of charged particles carries essential evolution information about the collision system. The Tsallis statistics is used to investigate the transverse momentum distribution of charged particles produced in Xe-Xe collisions at ffiffiffiffiffiffiffi sNN p = 5:44TeV. On this basis, we reproduce the nuclear modification factor of the charged particles. The calculated results agree approximately with the experimental data measured by the ALICE Collaboration.


Introduction
One of the major goals of high-energy nucleus-nucleus (AA) collisions is to study quark-gluon plasma (QGP) at high energy density and high temperature. The Large Hadron Collider (LHC) has performed different species of collisions at one or more energies, such as lead-lead, protonlead, and proton-proton collisions. The Xe-Xe ion collision [1, 2] at ffiffiffiffiffiffiffi s NN p = 5:44 TeV is a new collision experiment and is an intermediate-size collision system at the LHC. Since the mass number value of xenon is between proton and lead, it helps us to understand the system-scale effect of the final-state particle properties in ion collisions at high energy [3][4][5][6]. Compared with the sphere of the Pb nucleus, the deformation of the Xe nucleus is long and flattened in collisions. The deformed shape of Xe will provide us with different kinds of collision configurations. The deformed Xe nucleus will affect the initial condition of the reaction. How much impact does the deformation have on particle production and distribution? Many charged particles are produced and measured in the AA collisions. The investigation of the particle spectra is of great interest and is very helpful for comprehending the collision reaction mechanism and the particle production process in the different species of collision systems at different center-of-mass energies [7][8][9][10][11][12][13].
With respect to the final-state observations, the experimental transverse momentum p T spectrum is of great significance in understanding the production process of the moving particles. In past years, theoretical efforts have been carried out in statistical models to analyze the particle spectra over a broad range of collision energies [14-18]. At RHIC and LHC energies, the p T spectra have been investigated intensively in various collision systems like Au+Au, Pb+Pb, and pp at different energies. A statistical model can achieve some features in treating the multiparticle system in RHIC and LHC. Recently, the ALICE Collaboration reported the p T spectra and nuclear modification factors of charged particles produced in Xe-Xe collisions at ffiffiffiffiffiffiffi s NN p = 5:44 TeV [1]. The nuclear modification factor R AA is also an important observation and can provide information about the dynamics of QGP matter at extreme densities and temperatures [19][20][21][22][23][24][25][26].
In this paper, we discuss the p T spectra and the nuclear modification factor R AA in the Tsallis statistics. By the investigation of the p T spectra, we extract the parameters, which provide the calculation foundation for the nuclear modification factor R AA .

Description of the Particle Distribution in the Tsallis Statistics
The Tsallis statistics has been widely used to study the properties of final-state particles produced in nucleus-nucleus and proton-proton collisions at high energy [27][28][29][30]. In The Tsallis statistics, more than one version of the Tsallis distribution is used to investigate particle distributions. According to the Tsallis statistics, the number of the particles is where g and μ are the degeneracy factor and the chemical potential of the multiparticle system, respectively. T and q are the Tsallis temperature and the degree parameter of deviation from equilibrium, respectively. The first equation and second equation are two versions. The second equation (equation (1b)) can naturally meet the thermodynamic consistency [31][32][33]. At μ = 0, the transverse momentum distribution is The nuclear modification factor R AA acts as a probe to understand the nuclear medium effect in the AA collision and is a measure of the particle production modification. It is typically expressed as a ratio of the particle p T spectra in AA collisions to that in pp collisions: where N AA is the production yield in AA collisions and σ pp is the production cross-section in pp collisions. The average nuclear overlap function hT AA i is estimated via a Glauber model of nuclear collisions. The R AA is also expressed as where f in is the distribution of the initial particles produced at an early time of the hadronization. Then, these particles interact with the medium system. The function f fin is the distribution of the final-state particles, which no longer interact with each other.
According to the Boltzmann transport equation, the distribution of the particles f ðx, p, tÞ is The evolution of the particle distribution is attributed to its interaction with the medium particles. The terms v and F are the velocity and the external force, respectively. In relaxation time approximation, the collision term C½ f is given by where τ is the relaxation time. The Boltzmann local equilibrium distribution f eq is where T eq is the equilibrium temperature of the QCD phase transition. Considering A solution of the equation is where t f is the freeze-out time. The initial distribution is taken as the Tsallis distribution, i.e., equation (2). Therefore, the final-state distribution is Then, the nuclear modification factor R AA is obtained as Advances in High Energy Physics The equation is the calculation basis of the nuclear modification factor. In the relaxation time approximation, the R AA is derived in the Tsallis statistics.

Discussions and Conclusions
In this section, we discuss the transverse momentum spectra and the nuclear modification factor of the charged particles produced in Xe-Xe collisions at ffiffiffiffiffiffiffi s NN p = 5:44 TeV. The transverse momentum contributes significantly to the characterization of the matter formed in high energy collisions because p T is sensitive to the matter properties at an early time. The transverse momentum spectra in the kinematic range 0:15 < p T < 50 GeV/c and jηj < 0:8 are presented for nine centrality classes in Figure 1 Table 1. The nonequilibrium degree q is a constant value. The freeze-out time t increases with increasing collision centrality. The final-state transverse momentum spectra for different centralities are determined by the temperature T, at which there are no interactions between the final-state particles. By the analysis of the p T spectra, the thermodynamics parameters are extracted.
The dotted lines are the results of the Boltzmann statistics, which can agree with the experimental data in the low p T range.
The nuclear modification factor is also an important observation and is a measure of the particle-production modification. In Figure 1, we compare the p T spectra of the model results and the experiment data, and can extract the parameters, which are required in the calculation of the nuclear modification factor R AA . Figure 2 (11). The parameters used in the calculation are determined by the model results in Figure 1. The nuclear modification factor R AA depends strongly on the collision centrality. The R AA rises linearly at low p T (about below 2.2 GeV). At high p T , the R AA first declines linearly and then rises slowly. The model can approximately describe the nuclear modification factor at the high p T region, as shown in Figure 3. The dotted lines are the results of the Boltzmann statistics. Same as the above description of the transverse momentum spectra, they agree with the experimental data at low p T .
Both experimentally and theoretically, the study of the particle spectra can contribute to our understanding of the particle production and the evolution dynamics in the collision system. The Tsallis statistics has attracted extensive attention due to the investigation of final-state particles produced in nuclear collisions at high energies. Compared with Levy-Tsallis, Boltzmann, and Blast wave, the Tsallis distribution can describe the transverse momentum spectra at a large range. It can extract the temperature and the nonequilibrium degree, which provide the requirements of the R AA calculation. It is suc-cessful in explaining the experimental data of the transverse momentum spectra and can obtain some thermodynamics information, such as the temperature and the chemical potential. In our previous work [34][35][36][37], the statistics model is only used to study the transverse momentum spectra of particles produced in one or more collision systems at different energies. The present work is a new attempt. The model is improved by the Tsallis statistics in relaxation time approximation. Considering relaxation time approximation of the collision term, we achieve the final-state distribution by solving the Boltzmann transport equation, where the initial distribution is inserted consistently. And, the expression of the R AA calculation in the Tsallis statistics is derived. In our previous work [31][32][33][34], the Tsallis distribution can describe the p T distributions of particles produced in one or more collision systems, such as p, Cu, Au, and Pb collisions at various energies. Compared with these collision systems, the Xe nucleus has a moderate prolate deformation. But, p T distributions in Xe-Xe collisions can also be described well by the Tsallis distribution. The improved model can not only describe transverse momentum spectra but also

Data Availability
The used data in the model calculation are available and have been listed in Table 1.

Conflicts of Interest
The authors declare that they have no conflicts of interest. [7] F. Bellini and Alice Collaboration, "Testing the system size dependence of hydrodynamical expansion and thermal particle production with π, K, p, and ϕ in Xe -Xe and Pb-Pb collisions with ALICE," Nuclear Physics A, vol. 982, pp. 427-430, 2019.