J/psi Production in Quark-Gluon Plasma
Li Yan,Pengfei Zhuang,N. Xu
2006-01-01
Abstract:J/ψ Production in Quark-Gluon Plasma Li Yan, 1 Pengfei Zhuang, 1 and Nu Xu 2 Physics Department, Tsinghua University, Beijing 100084, China Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA (Dated: October 30, 2006) We study J/ψ production at RHIC and LHC energies with both initial production and regener- ation. We solve the coupled set of transport equation for the J/ψ distribution in phase space and the hydrodynamic equation for evolution of quark-gluon plasma. At RHIC, continuous regeneration is crucial for the J/ψ momentum distribution while the elliptic flow is still dominated by initial production. At LHC energy, almost all the initially created J/ψs are dissociated in the medium and regeneration dominates the J/ψ properties. PACS numbers: 25.75.-q, 12.38.Mh, 24.85.+p The goal of high-energy nuclear collisions is to iden- tify and study the equation of state of the Quark-Gluon Plasma (QGP) which is believed to exist at the early stage of our universe [1]. In s N N = 200 GeV Au + Au collisions, the observations of jet-quenching [2] and collective flow have demonstrated the formation of hot and dense matter with partonic collectivity [3, 4]. Local thermalization of the system created in heavy ion colli- sions is yet to be tested. Heavy flavors including charm and bottom quarks are powerful tools [5] for studying the early collision dynamics because their masses are much larger than possible excitation, temperature for example, of the system created in collisions at RHIC. Recent results have shown that heavy flavors are all pro- duced in the initial collisions [6, 7]. Therefore, study- ing open charm and charmonium production will yield important information on the properties of QGP. The J/ψ is a particularly sensitive probe of the early stages because its survial probability depends on the environ- ment. Lattice gauge theory calculations have indicated that J/ψs do exist above the critical temperature [8, 9]. The J/ψ suppression was first proposed as a direct signature to identify the QGP formation twenty years ago [10]. Besides the normal suppression induced by nuclear absorption, the J/ψs initially produced by hard processes are anomalously suppressed [11–17] by inter- actions in the hot medium. While charm quark produc- tion at the SPS is expected to be small, there are more than 10 c¯ pairs produced in a central Au + Au colli- c sion at RHIC, and is probably more than 200 pairs at the LHC [18]. These uncorrelated charm quark pairs in the QGP can be recombined to form J/ψs. Obvi- ously, regeneration will enhance the J/ψ yield and al- ter its momentum spectra. Recently, the regeneration approach for J/ψ production at RHIC has been widely discussed with different models, such as thermal pro- duction on the hadronizaton hypersurface according to statistic law [19–21], the coalescence mechanism [22] and the kinetic model [23, 24] which considers contin- uous J/ψ regeneration in a QGP. The medium created in high-energy nuclear collisions evolves dynamically. In order to extract information about the medium by analyzing the J/ψ distributions, both the hot and dense medium and the J/ψ produc- tion processes must be treated dynamically. In this Letter, we treat continuous regeneration of the J/ψ in a QGP self-consistently, including hydrodynamic evo- lution of the QGP itself. Both regenerated and initially produced J/ψs are simultaneously taken into account. Including both processes is important for RHIC since, unlike the collisions at the LHC, not all the initially produced J/ψs are destroyed. The massive J/ψs are unlikely fully thermalized with the medium. Thus their phase space distribution should be governed by transport rather than the kinetic equation of Refs. [23, 24]. To comprehensively treat the J/ψ distribution, we determine the J/ψ transport equation including both initial production and anoma- lous suppression as well as regeneration. The transport equation is then solved together with the hydrodynamic equation which characterizes the space-time evolution of the QGP. In order to compare with the experimen- tal measurements, we calculate the average transverse momentum squared, p 2 , elliptic flow, v 2 , and nuclear t modification factor, R AA , for J/ψ at both RHIC and LHC. In p + p collisions, about 30-40% of final state J/ψs are from the feed-down [25] of ψ and χ c . In order to simplify the numerical calculation, we neglect the ψ contribution and assume 40% of the final state J/ψ’s come from χ c [16]. Since the Ψ (=J/ψ, χ c ) is heavy, we use a classical Boltzmann-type transport equation to describe its evolution. The distribution function f Ψ (p t , x t , τ |b) in the central rapidity region and in the transverse phase space (p t , x t ) at fixed impact param- eter b is controlled by the equation ∂f Ψ /∂τ + v Ψ · f Ψ = −α Ψ f Ψ + β Ψ . The second term on the left-hand side arises from free-streaming of Ψ with transverse velocity v Ψ = p t / p 2 + m 2 , which leads to the “leakage” effect and t is needed to explain the averaged transverse momen- tum squared at the SPS [15]. The anomalous suppres- sion and regeneration mechanisms are reflected in the loss term α Ψ and gain term β Ψ , respectively. Suppose the medium locally equilibrates at time τ 0 , when nuclear absorption of the initially produced J/ψs has ceased. Absorption effects can be included in the initial distribution [15, 16], f Ψ (p t , x t , τ 0 |b), of the transport equation. To determine the initial condition at RHIC, we take τ 0 = 0.6 fm, the initial J/ψ produc- tion cross section σ pp = 2.61 ± 0.20 µb [26] and nuclear absorption cross section σ abs = 3 mb [27]. We also as- sume that a gN = 0.076 GeV 2 /fm and p 2 pp = 4.31 t GeV 2 [26] to describe the initial p t broadening due to multiple gluon scattering.