Pseudo-rapidity distributions of charged particles in asymmetric collisions using Tsallis thermodynamics

The pseudo-rapidity distributions of the charged particles produced in the asymmetric collision systems p+Al, p+Au and 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3$$\end{document}He+Au at sNN=200\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_\text {NN}} = 200$$\end{document} GeV are evaluated in the framework of a fireball model with Tsallis thermodynamics. The fireball model assumes that the experimentally measured particles are produced by fireballs following the Tsallis distribution and it can effectively describe the experimental data. Our results as well as previous results for d+Au collisions at sNN=200\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_\text {NN}} = 200$$\end{document} GeV and p+Pb collisions at sNN=5.02\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_\text {NN}} = 5.02$$\end{document} TeV validate that the fireball model based on Tsallis thermodynamics can provide a universal framework for pseudo-rapidity distribution of the charged particles produced in asymmetric collision systems. We predict the centrality dependence of the total charged particle multiplicity in the p+Al, p+Au and 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3$$\end{document}He+Au collisions. Additionally, the dependences of the fireball model parameters (y0a\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_{0a}$$\end{document}, y0A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_{0A}$$\end{document}, σa\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{a}$$\end{document} and σA\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{A}$$\end{document}) on the centrality and system size are studied.


Introduction
High-energy heavy-ion collisions provide a unique way to understand the origin of the universe.However, their processes cannot be directly observed in experiments.We can only study the collision process indirectly by analyzing the properties of the final particles produced in the collisions.The pseudo-rapidity distribution of charged particles is one of the important experimental observables.The study of this observable could lead to a better understanding of the properties of the particles produced in the collisions, the particle production mechanism and so on.There have been numerous works in previous studies using different models, such as HIJING [1], AMPT [2][3][4], EPOS-LHC [5], a multi-source thermal model [6,7], a new revised Landau hydrodynamics model [8], a 1 + 1-dimensional hydrodynamics model [9,10], a dynamical initial state model coupled to (3 + 1)D viscous relativistic hydrodynamics [11] and so on, to analyze the existing experimental data of pseudo-rapidity distributions of the charged particles [12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30].Although these models are based on different physical ideas, valuable physical information on the collision process has been extracted and learned.
Recently, a fireball model based on Tsallis thermodynamics was utilized to analyze the pseudo-rapidity distribution of charged particles measured in high-energy heavy-ion collisions [31][32][33].In our previous works [32,33], we used the fireball model to study the pseudo-rapidity distributions of the charged particles produced in p+p(p ) collisions for energies ranging from √ s NN = 23.6 GeV to 13 TeV and A+A collisions at the RHIC and LHC and extended the fireball model to the asymmetric collision systems, i.e., d+Au collisions at √ s NN = 200 GeV and p+Pb collisions at √ s NN = 5.02 TeV, by considering the asymmetric collision geometry configuration.In this paper, we utilize recently published data from the PHINEX Collaboration [30] at RHIC to systematically study the pseudo-rapidity distributions of the charged particles produced in asymmetric collision systems, including p+Al, p+Au and 3 He+Au collisions at √ s NN = 200 GeV.We also predict the total multiplicities of the charged particles from the fireball model and study their centrality dependence.Further, we analyze the centrality and system size dependencies of the fireball model parameters obtained from the pseudo-rapidity distributions of the charged particles.
The paper is organized as follows.In Sect.2, the fireball model with Tsallis thermodynamics is briefly introduced.In Sect.3, the fitting results of the fireball model and the total charged particle multiplicities extracted from the fireball model are shown.The dependences of the model parameters on the centrality and size of the collision systems are also presented.A brief conclusion is drawn in Sect. 4.

Theoretical descriptions
In the self-consistent Tsallis thermodynamics, the Tsallis distribution is proposed as a generalization of the Boltzmann-Gibbs distribution [34].To describe the transverse momentum spectrum of particles, the Tsallis distribution is written as [31][32][33] where g is the particle state degeneracy, V is the volume, m T = √ m 2 0 + p 2 T is the transverse mass and m 0 is the particle rest mass, y is the rapidity, q is the entropic factor, which measures the non-additivity of the entropy [34,35], is the chemical potential and T is the temperature.The Boltzmann distribution is recovered when q = 1 .We take = 0 because the multiplicities of + and − are equal and they are the majority of particles produced in the collision systems considered.For the middle rapidity y ≈ 0 , Eq. ( 1) can be rewrit- ten as (1) The parameters q and T are extracted from the experimental transverse momentum spectrum of the particles.
In the fireball model with Tsallis thermodynamics [31][32][33], the particles measured in the experiment were produced by fireballs following Tsallis distribution Eq. ( 1).The density distribution of these fireballs in the rapidity space is (y f ) , where y f is the rapidity of the fireball.Therefore the transverse momentum spectrum of particles can be written as where N is the total particle multiplicity and A is the normalization constant such that Sometimes, the experimental data are measured in the pseudo-rapidity space.To describe the experimental data dN d , we substitute the relation between rapidity and pseudorapidity [36] into Eq.( 3) and integrate the transverse momentum in the equation to obtain [32,33] where Because of the term of T cosh 2 y , Eq. ( 6) cannot be analytically integrated over p T and it is done numerically.
In this paper the asymmetric collision systems are studied, and the distribution (y f ) is assumed to be the sum of two asymmetric q-Gaussian functions [33], (2) (3) where y 0a(A) and a(A) are the centroid position and width of the fireball distribution in the direction of the light (heavy) nucleus beam, respectively.The normalization of Eq. ( 8) is handled by the normalization constant A in Eq. ( 3).x is the parameter to characterize the extent of asymmetry, which was first proposed in our previous work [33].In this work, we take q � = q as in [31][32][33].A representative figure of Eq. ( 8) is shown in Appendix 1 with the parameters obtained for the p+Al collisions at 0-5% centrality and √ s NN = 200 GeV.

Results and discussion
Because the data of the transverse momentum spectra of the charged particles produced in p+Al, p+Au and 3 He+Au collisions at √ s NN = 200 GeV have not been released yet, the data of the transverse momentum spectra of 0 produced from these collisions are obtained from [37] in this study.Using Eq. ( 2) by taking g = 1 and V as a free parameter, the parameters T and q for the Tsallis distribution are extracted and listed in Table 2 in property of the collision system.We take the parameters T and q from the closest centrality when the centrality of the particle transverse momentum spectrum and centrality of the charged particle pseudo-rapidity distribution are not the same.These two parameters and the fireball model with Tsallis thermodynamics, Eqs. ( 6) and ( 8), are then utilized to study the pseudo-rapidity distribution of the charged particles produced in the collisions.The corresponding values of x in Eq. ( 8) are also listed in Table 3   The curves are the results from Eqs. ( 6) and ( 8) √ s NN = 200 GeV shown in [33], i.e., in terms of larger errors, a fewer number of data points as well as a lower pseudo-rapidity coverage, which leads to larger uncertainties to the fireball model parameters and affects our analyses of the fireball model parameters versus collision centrality and the collision system size to some extent later in the following.The pseudo-rapidity distribution of the charged particles for centrality 5-10% is lower than the case for centrality 10-20% in some pseudo-rapidity regions for the 3 He+Au collisions, which is observed in Fig. 3.A larger x at centrality 5-10% compared with the others for the 3 He+Au collisions is also observed in Table 3 in Appendix 2. We emphasize that the same fitting protocol is applied for all the pseudo-rapidity distribution data of the charged particles.Because the collision system is asymmetric, the pseudo-rapidity distribution of the charged particles has significant forward/backward asymmetry.Fewer particles are produced in the direction of the light nucleus (p, 3 He) beam compared to the heavy nucleus (Al, Au) beam.As the d+Au collision system at √ s NN = 200 GeV and the p+Pb collision system at √ s NN = 5.02 TeV we studied in [33], the pseudo-rapidity distributions of the charged particles produced by these collision systems also become more symmetric from the central to peripheral collisions.This is because the peripheral collisions for asymmetric collision systems are more similar to the symmetric p+p collisions according to collision geometry.
We then evaluate the centrality dependence of the total multiplicities of the charged particles produced in these collision systems.Integrating Eq. ( 6) over the ∈ [−10, 10] we obtain the total multiplicity of the charged particles for each centrality from the fireball model.Because the corresponding experimental data are not yet available, we only analyze the results extracted from the fireball model and treat them as predictions.Figure 4 shows the total multiplicities of the charged particles calculated from the fireball model versus the collision centrality c. c = 0 represents the most central collisions, and c = 1 represents the most peripheral colli- sions.It can be observed that the fitting function taken from [14] can effectively describe the centrality dependence of the total multiplicities of the charged particles.As the centrality changes from the central to peripheral collisions, fewer charged particles are produced.
We also analyze both the centrality and system size dependence of the parameters ( y 0a , y 0A , a and A ) of the fireball model.In Fig. 5 the dependence of the fireball model parameters on the collision centrality in the p+Al, p+Au and of our previous work [33], the linear fittings are performed to guide the eyes in Fig. 5 and the fitting functions are listed in Table 1.The negative and positive slopes of the linear fittings for the fireball model parameters (y 0a , y 0A and A ) versus centrality are similar to those of the d+Au and p+Pb collisions in [33].In the following discussion, the results for d+Au and p+Pb are obtained from [33].It can be observed that the slopes of the linear fittings for parameter y 0a versus centrality are positive and the corresponding slopes of the linear fittings for parameter y 0A versus centrality are nega- tive for the p+Al, p+Au and d+Au collisions at √ s NN = 200 GeV.However, the slopes of the linear fittings for parameters (y 0a , y 0A ) reverse the signs, respectively, for the 3 He+Au collision at √ s NN = 200 GeV compared to the above-men- tioned cases.The slopes of the linear fittings for parameters (y 0a , y 0A ) are positive for p+Pb at √ s NN = 5.02 TeV.It can also be observed that there is a universal trend with increasing centrality for parameter a(A) except the lightest collision system p+Al, i.e., a increases with increasing centrality in the direction of the light nucleus beam and A decreases with The fitting function is from [14] and specified in the legend increasing centrality in the direction of the heavy nucleus beam.In the p+Al collision system, a and A have opposite trends with increasing centrality to their counterparts in the other collision systems.These different patterns indicate the complex dynamics in the asymmetric collisions relevant to GeV.The parameters for the p+p, d+Au and Au+Au collisions are taken from [33] Table 1 Results of the linear fits are shown in Fig. the combinations of the projectile and target as well as the collision energy, which needs more investigations.Figure 6 shows the collision system size dependence of the fireball model parameters at √ s NN = 200 GeV.For collision systems other than p+p, the parameters of the most central collisions are considered.In the p+p and Au+Au collisions, the parameters y 0a = y 0A = y 0 , a = A = , where y 0 and are the rapidity centroid and width of fireball distribution in the symmetric collision system, as detailed in [33].It can be deduced that when the light nucleus is p, y 0a decreases as the size of the heavy nucleus increases, whereas y 0A shows the oppo- site trend.This indicates that a larger heavy nucleus has stronger stopping power for p.When the heavy nucleus is Au, y 0a increases as the size of the light nucleus increases, whereas y 0A shows the opposite trend.This means that a larger light nucleus is more difficult to stop by Au but has a stronger stopping power for Au.For parameters a and A , no conclusive patterns are observed.We expect that more discussions can be added when the data quality of the pseudo-rapidity distributions of the charged particles is improved by experimentalists.These phenomena manifest the complex dynamics in the asymmetric collisions.

Summary
In this paper, we studied the pseudo-rapidity distributions of the charged particles produced in p + Al, p + Au and 3 He + Au collisions at √ s NN = 200 GeV using the fireball model with Tsallis thermodynamics.The model can well fit the experimental data from the asymmetric collisions.We also extracted the total multiplicities of the charged particles from the fireball model as predictions and analyzed their dependence on collision centrality.Notably, the data quality of the pseudo-rapidity distributions of the charged particles produced in the p + Al, p + Au and 3 He + Au collisions at √ s NN = 200 GeV affected our results to some extent.Combining our previous results of d+Au collisions at √ s NN = 200 GeV and p+Pb collision at √ s NN = 5.02 TeV, we analyzed the centrality and system size dependence of the fireball model parameters ( y 0a , y 0A , a and A ). Interesting patterns were revealed, which indicated the complex dynamics in the asymmetric collisions.Our results confirmed the conclusion made previously in [33] that the fireball model with Tsallis thermodynamics as a universal framework could also describe the pseudo-rapidity distribution of charged particles produced in asymmetric collision systems.
This work was supported in part by the National Natural Science Foundation of China (Nos.11905120, 11947416), the Natural Science Foundation of Sichuan Province (No. 2023NSFSC1322), the Natural Science Basic Research Plan in Shaanxi Province of China (No. 2023-JC-YB-012), the United States Department of Energy (# DE-FG02-93ER40773) and the NNSA (No. DENA0003841 (CENTAUR)).172 Page 2 of 8 Appendix 2. A representative figure of the transverse momentum spectra of 0 for p+Al collisions at √ s NN = 200GeV is shown in Appendix 3.It is worth noting that the transverse momentum spectrum of 0 is very similar to that of ± at √ s NN = 200 GeV for a given collision centrality, and the temperature parameter T extracted from 0 should reasonably characterize the(8) in Appendix 2. In Figs. 1, 2 and 3, the results of the pseudo-rapidity distributions of the charged particles from the fireball model with Tsallis thermodynamics for different centrality bins in p+Al, p+Au and 3 He+Au collisions at √ s NN = 200 GeV are shown.The fireball model effectively describes the experimental data within the errors.Notably, the data quality of the pseudo-rapidity distributions of the charged particles is not as good as that for the d+Au collisions at

Fig. 1 (
Fig. 1 (Color online) The pseudo-rapidity distributions of the charged particles produced in p+Al collisions at √ s NN = 200 GeV for different centralities.The symbols are experimental data taken from [30].

3
He+Au collisions at √ s NN = 200 GeV is shown.Inspired by the linear relation of the centrality dependence of the fireball model parameters for d+Au collisions at √ s NN = 200 GeV and p+Pb collisions at √ s NN = 5.02 TeV shown in Fig. 12

Fig. 4
Fig. 4 Total charged particle multiplicities produced in the p+Al, p+Au and 3 He+Au collisions at √ s NN = 200 GeV versus the collision centrality c.The squares are the fireball model results.The lines are the fitting results.The fitting function is from [14] and specified in the legend

Fig. 5 (Fig. 6 (
Fig. 5 (Color online) Centrality dependence of model parameters y 0a , y 0A , a and A in p+Al, p+Au and 3 He+Au collision at √ s NN = 200 GeV.The lines are the linear fit results to guide the eyes

Table 2
Parameters q and T for the p+Al, p+Au and 3 He+Au collisions at √ s NN = 200 GeV for different centralities

Table 3
Parameter x for the p+Al, p+Au and 3 He+Au collisions at √ s NN = 200 GeV for different centralities