arxiv: v2 [nucl-th] 25 Sep 2016

Transcription

1 EPJ manuscrit No. (will be inserted by the editor) Nuclear Modification Factor Using sallis Non-extensive Statistics Sushanta riathy, rambak Bhattacharyya 2, Prakhar Garg,a, Prateek Kumar, Raghunath Sahoo,b, and Jean Cleymans 2 Disciline of Physics, School of Basic Sciences, Indian Institute of echnology Indore, Simrol, M.P , India 2 UC-CERN Research Centre and Deartment of Physics, University of Cae own, Rondebosch 77, South Africa arxiv: v2 [nucl-th] 25 Se 26 October 6, 28 Introduction Abstract. he nuclear modification factor is derived using sallis non-extensive statistics in relaxation time aroximation. he variation of the nuclear modification factor with transverse momentum for different values of the non-extensive arameter, q, is also observed. he exerimental data from RHIC and LHC are analysed in the framework of sallis non-extensive statistics in a relaxation time aroximation. It is shown that the roosed aroach exlains the of all articles over a wide range of transverse momentum but does not seem to describe the rise in at very high transverse momenta. PACS q Relativistic heavy-ion collisions Cj Heavy-quark roduction in heavy-ion collisions One of the major goals of studying heavy-ion collisions at high-energies is to search for a deconfined state of quarks and gluons, also known as quark gluon lasma (QGP), and to study its roerties. he bulk roerties of the QGP are governed by light quarks and gluons. he heavy quarks act as robes for QGP roerties due to the fact that they witness the entire lasma evolution as they are roduced in the initial hard scattering and endure until hadronization. Also, as their time scale of thermalization is longer than that of light quarks, they can retain the entire interaction history more effectively. Similarly, the energy loss of light quarks becomes imortant to study the flavor deendence of the energy loss in heavy-ion collisions. he energy loss is more for high- heavy and light quark flavors [,2,3,4] due to interaction with the medium. Finally they aear as constituents of hadrons. he roagation of energetic quarks through the medium has been treated as Brownian motion which is described by means of Fokker-Planck equation. In this equation, the interaction is encoded in drag and diffusion coefficients. Many theoretical efforts have been made using Fokker-Planck equation to reroduce the exerimentally observed value of of heavy and light quarks [5,6,7,8,9,,,2,3]. Also, there have been studies to elucidate the dominant mode of energy loss [5] or the flavor deendence of energy a Present Address: Deartment of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, NY , USA b Raghunath.Sahoo@cern.ch (corresonding author) loss [4,5]. But, till date the issues are far from being settled [6]. he nuclear modification factor ( ) is a measure of the modification of article roduction. It can be reresented as = f fin f in, () where f in is the distribution of the highly energetic articles immediately after their formation and f fin is the distribution of the articles after the interaction with the medium. is defined as ( ) = (/N evt AA )d2 N AA /dyd ( N coll /σ inel NN ) d2 σ /dyd, (2) where d 2 N AA /dyd is the yield in A+A collisions, N coll is the number of binary nucleon-nucleon collisions averaged over the imact arameter range of the corresonding centrality bin calculated by Glauber Monte-Carlo simulation [7]. σnn inel is the inelastic cross section and d2 σ /dyd is the differential cross section for inelastic + collisions. N evt is the number of events in A+A collisions. If, AA =, this indicates that A+A collisions are mere suerosition of scaled + collisions. A deviation of from unity indicates the medium modification. It has been observed that high- article yields in Au+Au and Pb+Pb collisions at RHIC and LHC are suressed as comared to + collisions [8,9], which suggests the formation of a dense medium. In this work, we reresent the initial distribution of the energetic articles with the hel of sallis ower law

2 2 Sushanta riathy et al.: Nuclear Modification Factor... distribution arameterized by the sallis q arameter and the sallis temerature, remembering the fact that their genesis is due to very hard scatterings. We lug the initial distribution (f in ) in Boltzmann ransort Equation (BE) and solve it with the hel of Relaxation ime Aroximation (RA) of the collision term to find out the final distribution (f fin ). Hence, the ratio in Eq. exressible in terms of q, and relaxation time τ can be comuted and comared with the exerimentally observed values. he aer is organized as follows. In section 2, the nuclear modification factor is derived using RA of the BE and the exression for nuclear modification factor is derived. In section 3, fits to the exerimental data using the roosed model along with results and discussions are resented; and lastly, we summarize our findings in section 4. 2 Nuclear Modification Factor in Relaxation time aroximation (RA) he evolution of the article distribution owing to its interaction with the medium articles can be studied through Boltzmann transort equation, df(x,, t) dt = f t + v. xf + F. f = C[f], (3) where f(x,, t) is the distribution of articles which deends on osition, momentum and time. v is the velocity and F is the external force. x and are the artial derivatives with resect to osition and momentum, resectively. C[f] is the collision term which encodes the interaction of the robe articles with the medium. he Boltzmann ransort Equation has earlier also been used in relaxation time aroximation to study the time evolution of temerature fluctuation in a non-equilibrated system [2]. Assuming homogeneity of the system ( x f = ) and absence of external force (F = ), the second and third terms of the above equation become zero and Eq. 3 becomes, df(x,, t) dt = f t = C[f] (4) In relaxation time aroximation [2, 22], the collision term can be exressed as, C[f] = f f eq τ (5) where f eq is Boltzmann local equilibrium distribution characterized by a temerature eq. τ is the relaxation time, the time taken by a non-equilibrium system to reach equilibrium. With the ansatz in Eq. 5, Eq. 4 becomes, f t = f f eq τ (6) Solving the above equation in view of initial conditions i.e. at t =, f = f in and at t = t f, f = f fin ; leads to, f fin = f eq + (f in f eq )e t f τ, (7) where t f is the freeze-out time. Using Eq.7, the nuclear modification factor can be exressed as, = f fin = f ( eq + f ) eq e t f τ (8) f in f in f in Eq. 8 is the derived nuclear modification factor after incororating relaxation time aroximation, which is the basis of our analysis in the resent aer. It involves the (ower law-like) initial distribution and the equilibrium distribution. In this analysis, the initial distribution is arameterized using the thermodynamically consistent sallis distribution [23]. sallis statistics is widely used to analyse the detected article sectra in high-energy collisions starting from e + e, + to heavy-ions [24,25,26,27,28,29,3,3,32,33,34,35, 36,37,38,39,4,4]. Also, as the system stays away from thermal equilibrium during the formation of highly energetic articles immediately after the collision, the article distribution can be arameterized with the hel of sallis distribution. A thermodynamically consistent nonextensive sallis distribution function, to be used as the initial distribution, is given by, [23] f in = gv [ (2π) 2 m + (q ) m ] q q, (9) which is used for studying the article distribution stemming from the roton-roton collisions as discussed in Ref. [23]. he Boltzmann equilibrium distribution is given by f eq = gv (2π) 2 m e m eq () Here, V is the system volume, m = 2 + m2 is the transverse mass and q is the non-extensive arameter, which measures the degree of deviation from equilibrium. Using Eqs. 9 and (both for mid-raidity and for zero chemical otential) nuclear modification factor can be exressed as, e m eq = ( + (q ) m [ e m eq ) q q ( + (q ) m ) q q ] + e t f τ () his will be comared to exerimental results in the next section. 3 Results and Discussion o illustrate the formula given in Eq. we take as an examle the case of the J/ψ article. he variation of nuclear modification factor with transverse momentum for

3 Sushanta riathy et al.: Nuclear Modification Factor... 3 different values of non-extensive arameter is lotted in Fig.. For this figure, m = 3.96 GeV, eq =.6 GeV, =.7 GeV and t f /τ =.6 are taken. It is observed that for higher values of q, decreases for all the values of. his suggests that when the initial distribution remains closer to equilibrium, the suression becomes less. he model resented here fails to describe the increase in at larger transverse momenta, which is seen in exerimental data for light flavor hadrons q=. q=. q=. Fig. : Nuclear modification factor versus for different values of the non-extensive arameter using sallis Boltzmann distribution as shown in Eq.. Here m = 3.96 GeV, eq =.6 GeV, =.7 GeV and t f /τ =.6. We now roceed to the more detailed analysis of the exerimental data with the model roosed above. Keeing all the arameters free, we fit the sectra for different articles in different centralities for Pb+Pb and Au+Au collisions using Minuit class available in ROO library [42] to get a convergent solution. he convergent solution is obtained by χ 2 minimization technique. Here, q and t f /τ are the fitting arameters for the exerimental data. he equilibrium temerature eq is fixed to 6 MeV throughout the analysis. In our analysis, it is observed that the fitting of low for light flavor articles fails due to the reason that it involves different hysical rocesses such as regeneration, coalescence, shadowing etc., which are out for the scoe of resent formalism. But in heavy flavor articles, no such rocesses are involved at the discussed energies. hus our roosed model exlains successfully the heavy flavor data. Also, for very high-, our model fails to exlain the increase in (most rominent for K ± in Fig. 3). Fig. 2 shows the fitting of exerimental data using Eq. for π meson in most central Au+Au collisions at snn = 2 GeV. he derived exression for fits the data in intermediate to high- range. Also, in Fig. 2 we show the fitting of exerimental data for π + + π in most central Pb+Pb collisions at s NN = 2.76 ev. he model considered here, fits the data from intermediate to high- range. he fitting arameters are shown in table along with the χ 2 /ndf values Pb+Pb@2.76 ev (-5%) Au+Au@2 GeV (-5%) Fig. 2: (Color Online) Fitting of exerimental data for nuclear modification factor with our roosed model (Eq.) for π [43] (blue triangles) in Au+Au collisions at s NN = 2 GeV and π + + π [44] in Pb+Pb collisions at s NN = 2.76 ev (blue dots). he solid red line shows the fitting for blue triangles and the dotted red line shows the fitting for blue dots. Similarly, in Fig. 3 we fit the exerimental data for K + +K in most central Pb+Pb collisions at s NN = 2.76 ev. he roosed model fits the data for intermediate range. he derived exression of could not fit the data in high- range as an enhancement is observed in high-. he fitting arameters are shown table along with the χ 2 /ndf values. Fig. 4 shows the fitting of exerimental data for KS in most central Pb+Pb collisions at s NN = 2.76 ev. he roosed model fits the data for intermediate to high- range. he fitting arameters are shown table along with the χ 2 /ndf values. Fig.5 shows the fitting of exerimental data for + in most central Pb+Pb collisions at s NN = 2.76 ev. he roosed model fits the data for intermediate to high- range. χ 2 /ndf value and the fitting arameters are shown table. Fig. 6 shows the fitting of exerimental data using Eq. for Λ+ Λ in most central Pb+Pb collisions at s NN = 2.76 ev. he roosed model fits the data from intermediate to high- range starting from 2 GeV to 4 GeV. χ 2 /ndf value and the fitting arameters are shown table. In Fig. 7 we fit the exerimental data for D meson in most central, (-)% Au+Au collisions at s NN = 2 GeV. As less data oints are available, the fitting arameters cannot be established very well. hus the χ 2 /ndf value is very high comared to other articles, which can π

4 4 Sushanta riathy et al.: Nuclear Modification Factor ev + - K +K.8 Pb+Pb@2.76 ev (-5%) Fig. 3: (Color Online) Fitting of exerimental data for nuclear modification factor with our roosed model(eq.) for K + + K [44] in most central Pb+Pb collisions at snn = 2.76 ev (blue dots). he solid red line shows the fitting for blue dots Fig. 5: (Color Online) Fitting of exerimental data for nuclear modification factor with our roosed model (Eq.) for + [44] in most central Pb+Pb collisions at s NN = 2.76 ev (blue dots). he solid red line shows the fitting for blue dots Pb+Pb@2.76 ev (-5%) K s Fig. 4: (Color Online) Fitting of exerimental data for nuclear modification factor with our roosed model (Eq.) for K S in most central Pb+Pb collisions at s NN = 2.76 ev(blue dots). he solid red line shows the fitting for blue dots. be seen table. Also in Fig. 7 we show the fitting of exerimental data for D meson for (3-5)% central Pb+Pb collisions at s NN = 2.76 ev. he roosed model fits the data accurately for all the ranges as the enhancement is not involved which originate from regeneration through coalescence mechanism. χ 2 /ndf value and the fitting arameters are shown table. In Fig. 8 we show the fitting of exerimental data for J/ψ in minimum bias Pb+Pb collisions (-9% centrality) at s NN = 2.76 ev. he roosed model fits the data accurately for all ranges as the enhancement in is not observed in exerimental J/ψ data. he χ 2 /ndf value Pb+Pb@2.76 ev (-5%) Λ+Λ Fig. 6: (Color Online) Fitting of exerimental data for nuclear modification factor with our roosed model (Eq.) for Λ+ Λ in most central Pb+Pb collisions at s NN = 2.76 ev (blue dots). he solid red line shows the fitting for blue dots. and the fitting arameters are shown table. We have exlicitly checked that the roosed model also exlains the J/ψ for central Pb+Pb collisions at s NN = 5.2 ev.

5 Sushanta riathy et al.: Nuclear Modification Factor ev (3-5%) Au+Au@2 GeV (-%) Fig. 7: (Color Online) Fitting of exerimental data for nuclear modification factor with our roosed model (Eq.) for D meson [45] (blue triangles) in most central Au+Au collisions at s NN = 2 GeV and D meson [46] in (3-5)% central Pb+Pb collisions at s NN = 2.76 ev (blue dots). he solid red line shows the fitting for blue triangles and the dotted red line shows the fitting for blue dots. D And hence, t f /τ becomes less and the degrees of freedom are shared with other arameters used in the fit. t f /τ Pb+Pb@2.76 ev (Most central) mass (GeV/c ) Fig. 9: (Color Online) t f /τ as a function of article mass for most central Pb+Pb collisions at s NN = 2.76 ev Pb+Pb@2.76 ev (-2)% Fig. 8: (Color Online) Fitting of exerimental data for nuclear modification factor with our roosed model (Eq.) for J/ψ [47] in Pb+Pb collisions at s NN = 2.76 ev (blue dots) with centrality (-9)%. he solid red line shows the fitting to exerimental data. Finally, we have shown the variation of t f /τ with mass in Fig. 9. It is observed that t f /τ decreases with increasing the article mass, which suggests that the heavy articles have more relaxation time comared to lighter articles. Note here, that t f is the freeze-out time, and τ is the relaxation time, which differs from article to article, as the equilibration deends on the article secies and their interaction with the rest of the medium. Intuitively, the heavier the article is, the more is the relaxation time. J/ψ 4 Summary and Conclusion In this work, we reresent the initial distribution of the energetic articles with the hel of sallis ower law distribution arameterized by the sallis q arameter and the sallis temerature, remembering the fact that their genesis is due to very hard scatterings. We lug the initial distribution (f in ) in Boltzmann ransort Equation (BE) and solve it with the hel of Relaxation ime Aroximation (RA) of the collision term to find out the final distribution (f fin ). Hence, the ratio in Eq. exressible in terms of q, and relaxation time, τ can be comuted and comared with the exerimentally observed values. he variation of nuclear modification factor with transverse momentum for different values of non-extensive arameter, is also observed. he suression is found to be higher for a system with higher degree of deviation from equilibrium. Also, we analyse the exerimental data from RHIC and LHC with calculated nuclear modification factor. It is observed that the calculated exlains accurately for heavy flavor articles in all range, but it can only exlain for light flavor articles in intermediate to high- range with an excetion for kaons, where the enhancement at high- can not be exlained with the resent model. he relaxation time, as found from data, is higher for heavy articles. Acknowledgements S acknowledges the financial suort by DS INSPIRE rogram of Govt. of India. B acknowledges the discussion with Dr. Santosh K. Das.

6 6 Sushanta riathy et al.: Nuclear Modification Factor... able : Centrality, χ 2 /ndf and different extracted arameters after fitting Eq. to the data of different articles for Pb+Pb collisions and Au+Au collisions at s NN = 2.76 ev and s NN = 2 GeV, resectively. Pb+Pb 2.76 ev Particle Centrality(%) χ 2 /ndf t f /τ q (GeV) π + + π ± ± ±.43 K + + K ± ± ±.33 Ks ±.73.4 ± ± ± ± ±. Λ + Λ ± ± ±.52 D ± ± ±.568 J/ψ ± ± ±.48 Au+Au 2 GeV π ±.56.5 ± ±.55 D ± ± ±.688 References. I. Arsene et al. [BRAHMS Collaboration], Phys. Lett. B 65, 29 (27). 2. J. D. Bjorken, Phys. Rev. D 27, 4 (983). 3. M. Pluemer, M. Gyulassy and X. N. Wang, Nucl. Phys. A 59, 5C (995). 4. R. Baier, Y. L. Dokshitzer, S. Peigne and D. Schiff, Phys. Lett. B 345, 277 (995). 5. J. e. Alam, P. Roy and A. K. Dutt-Mazumder, arxiv: heh/ B. Svetitsky, Phys. Rev. D 37, 2484 (988). 7. M. Golam Mustafa, D. Pal and D. K. Srivastava, Phys. Rev. C 57, 889 (998) 3499(E) (998). 8. H. van Hees, V. Greco and R. Ra, Phys. Rev. C 73, 3493 (26). 9. H. van Hees, M. Mannarelli, V. Greco and R. Ra, Phys. Rev. Lett., 923 (28).. S. K. Das, J. e. Alam and P. Mohanty, Phys. Rev. C 8, 5496 (29).. S. K. Das, F. Scardina, S. Plumari and V. Greco, Phys. Lett. B 747, 26 (25). 2. W. M. Alberico, et al., Eur. Phys. J. C 73, 248 (23). 3. S. Mazumder,. Bhattacharyya, J. e. Alam and S. K. Das, Phys. Rev. C 84, 449 (2). 4. B. G. Zakharov, JEP Lett. 96, 66 (23). 5. B. Z. Koeliovich, I. K. Potashnikova and I. Schmidt, Phys. Rev. C 82, 379 (2). 6. I. Kolbe and W. A. Horowitz, arxiv:5.933 [he-h]. 7. R. J. Glauber and G. Matthiae, Nucl. Phys. B 2, 35 (97). 8. S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. Lett. 9, 723 (23). 9. K. Aamodt et al. (ALICE Collaboration), Phys. Lett. B 696, 3 (2). 2.. Bhattacharyya, P.Garg, R. Sahoo, P. Samantray, Eur. Phys. J. A 52, 283 (26). 2. R. Balescu, Equilibrium and Non-Equilibrium Statistical Mechanics, John Wiley and Sons, USA (975). 22. W. Florkowski and R. Ryblewski, Phys. Rev. C 93, 6493 (26) 23. J. Cleymans, D. Worku, J. Phys. G 39, 256 (22). 24. D. hakur, S. riathy, P. Garg, R. Sahoo and J. Cleymans, arxiv: [he-h]. 25. I. Bediaga, E.M.F. Curado, J.M. de Miranda, Physica A 286, 56 (2). 26. G. Wilk and Z. W lodarczyk, Acta Phys. Polon. B 46, 3 (25). 27. K. Ürmössy, G.G. Barnaföldi,.S. Biró, Phys. Lett. B 7, (2). 28. K. Ürmössy, G.G. Barnaföldi,.S. Biró, Phys. Lett. B 78, 25 (22). 29. P. K. Khandai, P. Sett, P. Shukla, V. Singh, Int. Jour. Mod. Phys. A 28, 3566 (23). 3. B.-C. Li, Y.-Z. Wang and F.-H. Liu, Phys. Lett. B 725, 352 (23). 3. L. Marques, J. Cleymans and A. Deman Phys. Rev. D 9, 5425 (25). 32. B. I. Abelev et al. (SAR collaboration), Phys. Rev. C 75, 649 (27). 33. A. Adare et al. (PHENIX collaboration), Phys. Rev. D 83, 524 (2). 34. A. Adare et al. (PHENIX collaboration), Phys. Rev. C 83, 6493 (2). 35. K. Aamodt et al. (ALICE collaboration), Phys. Lett. B 693, 53 (2). 36. K. Aamodt et al. (ALICE collaboration), Eur. Phys. J C 7, 655 (2). 37. V. Khachatryan et al. (CMS collaboration), J. of High Eng. Phys. 2, 4 (2). 38. V. Khachatryan et al. (CMS collaboration), Phys. Rev. Lett. 5, 222 (2). 39. G. Aad, et al. (ALAS collaboration), New J. Phys. 3, 5333 (2). 4. B. Abelev et al. (ALICE collaboration), Phys. Rev. Letts. 9, 2523 (22). 4.. Bhattacharyya, J. Cleymans, A. Khuntia, P. Pareek and R. Sahoo, Eur. Phys. J. A 52, 3 (26). 42. CERN ROO V.5.34/32 (June 23, 25) Package: htt://root.cern.ch. 43. A. Adare et al. (PHENIX Collaboration), Phys. Rev. Lett., 2323 (28). 44. B. B. Abelev et al. (ALICE Collaboration), Phys. Lett. B 736, 96 (24). 45. L. Adamczyk et al. (SAR Collaboration), Phys. Rev. Lett. 3, 423 (24). 46. B. B. Abelev et al. (ALICE Collaboration), Phys. Rev. C 9, 3494 (24). 47. B. B. Abelev et al. (ALICE Collaboration), Phys. Lett. B 734, 34 (24).