Arxiv:2006.15333V4 [Hep-Ph]
Total Page:16
File Type:pdf, Size:1020Kb
Excitation Functions of Tsallis-like Parameters in High-Energy Nucleus-Nucleus Collisions 1 1, 2, Li-Li Li , Fu-Hu Liu ∗ and Khusniddin K. Olimov ∗ 1 Institute of Theoretical Physics & Collaborative Innovation Center of Extreme Optics & State Key Laboratory of Quantum Optics and Quantum Optics Devices, Shanxi University, Taiyuan 030006, China; shanxi [email protected] (L.L.L.) 2 Laboratory of High Energy Physics, Physical-Technical Institute of SPA “Physics-Sun” of Uzbek Academy of Sciences, b Chingiz Aytmatov str. 2 , 100084 Tashkent, Uzbekistan ∗ Correspondence: [email protected] or [email protected] (F.H.L.); [email protected] (K.K.O.) Abstract: The transverse momentum spectra of charged pions, kaons, and protons produced at mid-rapidity in central nucleus-nucleus (AA) collisions at high energies are analyzed by considering particles to be created from two participant partons which are assumed to be contributors from the collision system. Each participant (contributor) parton is assumed to contribute to the transverse momentum by a Tsallis-like function. The contributions of the two participant partons are regarded as the two components of transverse momentum of the identified particle. The experimental data measured in high-energy AA collisions by international collaborations are studied. The excitation functions of kinetic freeze-out temperature and transverse flow velocity are extracted. The two param- eters increase quickly from 3 to 10 GeV (exactly from 2.7 to 7.7 GeV) and then slowly at above 10 GeV with ≈ ≈ the increase of collision energy. In particular, there is a plateau from near 10 GeV to 200 GeV in the excitation function of kinetic freeze-out temperature. Keywords: Excitation functions of related parameters, participant parton, kinetic freeze-out temperature, transverse flow velocity PACS: 12.40.Ee, 13.85.Hd, 24.10.Pa 1 Introduction freeze-out, a phase transition from QGP to hadrons oc- curred in the system, so the composition and ratio of High-energy collider experiments are designed to study various particles remain unchanged. In the stage of ki- the strongly interacting matter at high temperatures netic freeze-out, elastic collisions among particles stop, and densities [1]. The deconfinement of colliding so their pT and then mT spectra are unchanged [8, 10]. hadrons into quark-gluon plasma (QGP), which then Therefore, by studying the pT (mT ) spectra, we can ob- rapidly expands and cools down [2], is conjectured to be tain some useful information, such as the effective tem- created at such extreme collision energies [3, 4, 5, 6]. In perature (T ), the chemical freeze-out temperature (Tch), high energy and nuclear physics, the study of transverse and the kinetic freeze-out temperature (T0 or Tkin) of arXiv:2006.15333v4 [hep-ph] 18 Apr 2021 [momentum (pT ) or mass (mT )] spectra of charged par- the system, as well as the transverse flow velocity (βT ) of ticles produced in nucleus-nucleus (AA) collisions is very the final state particles. The temperature in which we do important. In particular, the AA collision process at the not exclude the contribution of transverse flow is called Relativistic Heavy Ion Collider (RHIC) and the Large the effective temperature which is related to the kinetic Hadron Collider (LHC) provides a good opportunity to freeze-out temperature. The temperatures in the stages study the signals and characteristics of QGP generation, of chemical and kinetic freeze-outs are called the chem- so as to indirectly study the system evolution and the ical and kinetic freeze-out temperatures respectively. reaction mechanism of particle generation. It is very important to study the behavior of T0 and During the time evolution of collision system [7, 8, 9], βT due to their relation to map the phase diagram of the stages of kinetic freeze-out and chemical freeze-out Quantum Chromodynamics (QCD), though Tch is usu- are two important processes. In the stage of chemical ally used [11, 12, 13, 14, 15, 16] in the phase diagram. 1 In order to extract T0 and βT , and study their depen- The formalism and method are shortly described in Sec- dence on energy, we can analyze the pT (mT ) spectra tion 2. Results and discussion are given in Section 3. of particles using different models. These models in- In Section 4, we summarize our main observations and clude, but are not limited to, the blast-wave model with conclusions. Boltzmann-Gibbs statistics [17, 18] or Tsallis statis- tics [19, 20, 21], as well as other alternative meth- ods [22, 23, 24, 25, 26] based on the standard distribu- 2 Formalism and method tion or Tsallis distribution. Here, the standard distribu- tion denotes together the Boltzmann, Fermi-Dirac, and The Tsallis distribution has different forms or revi- Bose-Einstein distributions. The alternative method re- sions [44, 45, 46, 47], we have the Tsallis-like distribution gards the intercept of T versus m0 as T0, and the slope of pT at mid-y to be of pT versus m as βT , where m0, pT , and m denote h i h i 2 1/(q 1) the rest mass, mean p , and mean energy of the given d N dN (q 1)(mT µ m0) − − T mT 1+ − − − , particles, respectively. dydpT ∝ dy T (1) In our recent work [27, 28], the blast-wave model with Boltzmann-Gibbs statistics or Tsallis statistics and where N denotes the number of particles, the standard distribution have been used to analyze 2 2 mT = pT + m0 (2) the spectra of particles produced in high-energy proton- q proton (pp) and AA collisions. The related parameters can be obtained using pT , were extracted and their excitation functions were ob- 1 q =1+ (3) tained. Not only the blast-wave model [17, 18, 19, 20, 21] n but also the alternative method [22, 23, 24, 25, 26] can is an entropy index that characterizes the degree of equi- be used to extract T0 and βT , though an effective tem- librium or non-equilibrium, n is a parameter related to perature T is used in the latter. The alternative method q, and µ is the chemical potential. In particular, in is partly a new one, in which the extractions of both T0 the expression of m µ m , m is simplified from T − − 0 T and βT are based on T [22, 23, 29] and the related de- m cosh y because cosh y 1 at mid-y. T ≈ rived quantities such as p and m. We have the probability density function of p at h T i T mid-y to be Due to the importance of T0 and βT and their excita- 1/(q 1) tion functions, we use a new method in the framework of 1 dN (q 1)(mT µ m0) − − mT 1+ − − − . multisource thermal model [30] to describe the pT (mT ) N dpT ∝ T spectra of identified particles in this work. Considering (4) the contributions of two participant (contributor) par- Empirically, to fit the spectra of p at mid-y in this tons to p of a given particle, we regard the two con- T T work, Eq. (4) can be revised as tributions as the two components of pT . The pT (mT ) 1/(q 1) spectra of identified particles (concretely charged pions, (q 1)(m µ m ) − − f(p ,T )= Cma0 1+ − T − − 0 , kaons, and protons) produced at mid-rapidity (mid-y) T T T in central AA collisions which include gold-gold (Au- (5) Au) collisions at the Alternating Gradient Synchrotron where C is the normalization constant, a0 is a new non- (AGS), lead-lead (Pb-Pb) collisions at the Super Pro- dimensional parameter that describes the bending de- ton Synchrotron (SPS), Au-Au collisions at the RHIC, gree of the distribution in low-p region (p = 0 1 T T ∼ and Pb-Pb and xenon-xenon (Xe-Xe) collisions at the GeV/c), which is introduced artificially and tested in our LHC are studied. The center-of-mass energy per nu- a0 recent work [48, 49], and mT is revised from mT due to cleon pair, √sNN , considered by us is from 2.7 GeV the introduction of the revised index a0. Because of the to 5.44 TeV. After fitting the experimental data mea- limitation of the normalization, changing the bending sured by the E866 [31], E895 [32, 33], E802 [34, 35], degree in low-pT region will change the slope in high-pT NA49 [36, 37], STAR [38, 39, 40], and ALICE Collabo- a0 region. Although writing CmT in Eq. (5) is not ideal, rations [41, 42, 43], we analyze the tendency of param- as it yields a fractional power unit in C, we have no eters. suitable method to scale out the unit by e.g. m0 due to The remainder of this paper is structured as follows. the nonlinear relationship between mT and m0 shown in 2 Eq. (2). In Eq. (5), the other parameters such as q and where f1,2(pt1,pt2,T ) denotes the united probability a0 do not appear in the function name for the purpose density function of pt1 and pt2. Further, we have the of convenience. In this work, we call Eq. (5) the revised probability density function of pT to be Tsallis-like function. 2π In the framework of the multisource thermal fpT (pT ,T )= fpT ,φ(pT ,φ,T )dφ Z0 model [30], we assume that two participant partons take 2π part in the collisions. Let pt1 and pt2 denote the com- = pT f1(pT cos φ, T )f2(pT sin φ, T )dφ. ponents contributed by the first and second participant Z0 (9) parton to pT respectively, where pt1 (pt2) is less than the transverse momentum of the participant parton.