-
Due to their large masses, heavy flavors have unique advantages in both experimental and theoretical studies of quantum chromodynamics (QCD). Since
J/ψ was initially proposed as a probe of the deconfined matter known as “quark-gluon plasma” (QGP) [1], its anomalous yield suppression through parton scattering and enhancement from the recombination of charm and anti-charm quarks in QGP have been widely studied in experiments [2–5] and theoretical models [6–12]. Charmonium is produced by an initial hard-scattering process referred to as “primordial production”, which occurs in hadronic collisions with spectator nucleons [13]. Charmonium states are assumed to form before QGP reaches local equilibrium, they then undergo inelastic scatterings and color screening effects when moving through QGP, which results in dissociations and also transitions between different eigenstates (J/ψ ,ψ′ ,χc ) [14–20]. These eigenstates are finally detected in experiments through their decay into dileptons. At LHC energies, abundant charm pairs are produced in nuclear collisions, this significantly enhances the probability of c andˉc combining to generate newJ/ψ s in the QGP [21–23], in a process called “regeneration”. As most charm quarks are distributed in low and middlepT regions, the regeneration process dominates the nuclear modification factor and the collective flows ofJ/ψ in the low and middlepT bins [24]. In the highpT bin, charmonia are mainly produced from the initial hadronic collisions [10].More experimental data about the charmonium excited state
ψ(2S) has been measured in√sNN=2.76 TeV [25] and 5.02 TeV [26] Pb-Pb collisions, across different centralities and transverse momentum bins. The ratioRψ(2S)AA/RJ/ψAA has been measured for these two collision energies, and there is considerable discrepancy between the values reported. At 2.76 TeV,Rψ(2S)AA/RJ/ψAA becomes larger than unity in the most central collisions when3<pT<30 GeV/c. At 5.02 TeV, the ratio is∼0.5 in a similar centrality and momentum bin. Both sets of experimental data feature large error bars, which prevents any solid conclusions being drawn. In contrast toJ/ψ ,ψ(2S) is a loosely bound state. Its wavefunction is significantly modified by the high-temperature medium, which can obscure its dissociation and regeneration rates. Having a smaller binding energy,ψ(2S) is thermally produced more frequently in the lower temperature region thanJ/ψ , and inherits larger collective flows from the bulk medium [17, 27, 28]. This sequential regeneration can affect thepT dependence of the ratioRψ(2S)AA/RJ/ψAA .In this paper, a two-component transport model is employed to study both
J/ψ andψ(2S) production in different centralities and momentum bins for√sNN=5.02 TeV Pb-Pb collisions. The decay rates of excited states are updated with a more realistic formula instead of employing a survival temperatureTd above which no charmonia can survive. This improvement can explain effectively the ratio ofψ(2S)/J/ψ in 5.02 TeV Pb-Pb collisions. The rest of this paper proceeds as follows. In Sec. 2, the details of an improved Boltzmann transport model for charmonium evolution, and the hydrodynamic equations for QGP expansion, are introduced. In Sec. 3, realistic calculations forJ/ψ andψ(2S) in√sNN=5.02 TeV Pb-Pb collisions are performed and the results are compared with experimental data. A final summary is given in Sec. 4. -
Heavy quarkonium evolution in phase-space has, though the use of Boltzmann transport models, been well-studied for hot deconfined mediums in the SPS [9] and LHC [29, 30], for both p-Pb and Pb-Pb collisions. Focusing on hot-medium effects, one can start quarkonium evolutions after their production in hard-scattering processes. The three-dimensional transport equation for charmonium evolution is simplified as,
[cosh(y−η)∂∂τ+sinh(y−η)τ∂∂η+vT⋅∇T]fΨ=−αΨfΨ+βΨ,
(1) fΨ is the Ψ phase-space density. y and η are the rapidities in momentum and coordinate space, respectively.vT=pT/ET=pT/√m2Ψ+p2T is the transverse velocity of charmonium, which represents the leakage effects for a cooling system with finite size, i.e., charmonia with large velocities tend to escape from the thermal medium instead of dissociating. Primordially-produced charmonia in the initial hadronic collisions suffer color screening effects and parton inelastic scatterings, which are both included in the decay rateαΨ ,αΨ=1ET∫d3k(2π)3EgσgΨ(p,k,T)FgΨ(p,k)fg(k,T),
(2) where
Eg andfg are gluon energy and density in the thermal medium, respectively.FgΨ is the flux factor. In the expanding QGP,uμ represents the four-velocity of the fluid. The gluonic Ψ cross-section in a vacuum is extracted from the perturbative calculation with a Coulomb potential approximation. For the thermal medium, this paper takes a similar approach to Ref. [9], using a reduced charmonium binding energy of the form,σgΨ(w)=A0(w/ϵΨ−1)3/2(w/ϵΨ)5
(3) with
A0=(211π/27)(m3cϵΨ)−1/2 andϵΨ representing the binding energy of Ψ. The charm quark mass is taken as the mass of a D meson, to fit the binding energy of charmonium in a vacuum.w=pμΨpgμ/mΨ is the gluon energy in the Ψ rest-frame. In Fig.1, theJ/ψ decay rateαΨ is compared with the quasi-free dissociation rate [27]. MostJ/ψ s can survive in the relatively low temperature region, and two different transport models of the decay rate present similar final results for the regionT<300 MeV, where most QGP and charmonia are located [7, 30].Figure 1.
J/ψ decay rate in the thermal medium as a function of temperature T. The decay rate from quasi-free dissociation is shown for comparison. The solid line is from the improved version of the transport model, developed by TSINGHUA Group [9, 10, 31], and the dotted line is from the calculations of TAMU Group [6, 7, 27].The heavy quark potential
V(r,T) can be partially screened by the thermal medium, especially at the large distances and high temperatures suggested by lattice QCD calculations [32]. Charmonium bound states may disappear sequentially in the static medium. The maximum survival temperature of a certain bound state is called the "dissociation temperature"Td , above this the bound state disappears. In nuclear collisions, the assumption that no bound states survive atT>Td strongly suppresses theψ(2S) production, where no excited states can survive inside the QGP atT>Tχc,ψ(2S)d≈1.1Tc . Similarly to fast cooling systems, charmonium states might survive in the regionT>Td if the medium quickly cools down to belowTd . In this work, the approximation ofαΨ(T>Td)=+∞ employed in [33] is replaced by a large but finite value, shown in Fig.1. The new decay rate increases the survival probability of excited states, and weakly affects the production ofJ/ψ s due to their largeTd . Theψ(2S) decay rate is extracted from the nuclear geometry scale,αψ(2S)=αJ/ψ×⟨r⟩2ψ(2S)/⟨r⟩2J/ψ , andχc is found in a similar way. The mean radius of charmonia in a vacuum is calculated with the potential model⟨r⟩J/ψ,χc,ψ(2S)=(0.5,0.72,0.9) fm [14].At LHC energies, many charm pairs are also produced in Pb-Pb collisions, which can significantly increase the recombination of uncorrelated charm and anti-charm quarks in the QGP. This process is included in Eq.(1) by the term
βΨ . The regeneration rate depends on both charm and anti-charm quark densities in the QGP, and also their recombination probability. At high temperatures, charmonium binding energies are reduced significantly, which suppresses the regeneration probability of charmonium. Asψ(2S) s are loosely bound, they are thermally produced in the hadronization of QGP. Charm quarks with color charges strongly couple with the QGP and lose energy. In relativistic heavy ion collisions, large quench factors and collective flows for charmed mesons have been observed [34–36]. Therefore, one can approximately take the kinetically thermalized phase-space distribution for charm quarks as being atτ⩾τ0 , whereτ0 is the time-scale of the QGP local equilibrium [37]. As heavy quarks are rarely produced in the thermal medium due to their large masses, the total number of charm pairs is conserved with spatial diffusions inside the QGP [38]. The spatial density is controlled by the conservation equation.∂μ(ρcuμ)=0.
(4) The initial charm quark density at
τ0 is obtained from nuclear geometry,ρc(xT,η,τ0)=dσcˉcppdηTA(xT)TB(xT−b)cosh(η)τ0,
(5) where
TA andTB are the thickness functions of two colliding nuclei, with the definitionTA(B)(xT)= ∫∞−∞dzρA(B)(xT,z) .ρA(B)(xT,z) is understood to be the Woods-Saxon nuclear density. The rapidity distribution of charm pairs in√sNN=5.02 TeV pp collisions is obtained by interpolation of the experimental data for 2.76 TeV and 7 TeV collisions,dσcˉcpp/dy=0.86 mb in the central rapidity region|y|<0.9 , and 0.56 mb in the forward rapidity region2.5<|y|<4 [39].The momentum distribution of charmonium primordially produced in A-A collisions is scaled from its distribution in pp collisions. The parametrization of the charmonium initial distribution at 5.02 TeV is similar to the form it takes at 2.76 TeV,
d2σJ/ψppdy2πpTdpT=fNormJ/ψ(pT|y)⋅dσJ/ψppdy,
(6) fNormJ/ψ(pT|y)=(n−1)π(n−2)⟨p2T⟩pp[1+p2T(n−2)⟨p2T⟩pp]−n.
(7) The charmonium rapidity differential cross-section at 5.02 TeV is
dσJ/ψpp/dy= 5.0 µb in the central rapidity|y|<1 , and 3.25 µb in the forward rapidity2.5<|y|<4 , these were found through an interpolation between the experimental results for 2.76 TeV [40] and 7 TeV [41] pp collisions.fNormJ/ψ(pT|y) is the normalized transverse momentum distribution of charmonium with rapidity y. The mean squared transverse momentum⟨p2T⟩ and the parameter n are calculated as⟨p2T⟩pp|y=0=12.5 (GeV/c)2 andn=3.2 . For the charmonium momentum distribution in other rapidities, which have less constraints,⟨p2T⟩pp(y) is determined by the relation,⟨p2T⟩J/ψpp(y)=⟨p2T⟩J/ψpp|y=0×[1−(yymax)2],
(8) where
ymax=cosh−1(√sNN/(2ET)) , and is the maximum rapidity of charmonium in pp collisions with zero transverse momentum. As the masses of charmonium excited states (χc,ψ(2S) ) are similar toJ/ψ s, their initial momentum distributions are approximated to be the same as those in Eqs.(6,7).In nuclear collisions, the charmonium initial distribution is also modified by shadowing effects in the nucleus [42, 43]. This paper employs the EPS09 NLO model [44] to generate the modification factors for primordially-produced charmonium in
√sNN=5.02 TeV Pb-Pb collisions. This suppression factor is∼0.8 depending on the impact parameter. For the regeneration, shadowing effects reduce the number of charm pairs by around 20%, and suppress the regeneration by a factor of∼0.82 .The expanding QGP background of charmonium evolution is simulated with the (2+1)-dimensional ideal hydrodynamic equations in the transverse plane, under the assumption of Bjorken expansion in the longitudinal direction.
∂μTμν=0,
(9) Tμν=(e+p)uμuν−gμνp , and is the energy-momentum tensor. e and p are the energy density and pressure, respectively.uμ is the four-velocity of the QGP fluid, which can affect charm quark spatial diffusions through Eq.(4) as well as charmonium regeneration. It also determines the collective flows of light hadrons, charmed mesons and regenerated charmonia. The deconfined matter is treated as an ideal gas of massless gluons, u and d quarks and strange quarks with massms=150 MeV [45]. Hadron gas is an ideal gas made up of all known hadrons and resonances with masses of up to 2 GeV [46]. The two phases are connected by a first-order phase transition with a critical temperature ofTc=170 MeV. The initial maximum temperature of the QGP is calculated asT0(xrmT=0,τ0)=510 MeV in the central rapidity|y|<2.4 , and 450 MeV in the forward rapidity2.5<|y|<4 . Hereτ0=0.6 fm/c, and is the time-scale for the QGP to reach local equilibrium [37]. The lifetime of QGP is∼10 fm/c for the most central Pb-Pb collisions at√sNN=5.02 TeV. -
With the transport model for charmonium evolution and the hydrodynamic equations for QGP collective expansion, one can obtain realistic nuclear modification factors for charmonia in heavy ion collisions. In the left-hand panel of Fig.2, primordially-produced charmonia undergo dissociations from peripheral to central Pb-Pb collisions, shown with the dotted line. The regeneration process
c+ˉc→J/ψ+g is plotted with the dashed line and is proportional to the number of charm quark pairs in the QGP, it dominates totalJ/ψ production in central collisions. The experimental data in the left-hand panel of Fig.2 is inclusive production; it includes the non-prompt contribution from B hadron decays, which contributes around 10% of the final inclusive yield. The detailed momentum dependence of the non-prompt fraction inJ/ψ inclusive production is fitted asfB=NB→J/ψpp/(Npromptpp+ NB→J/ψpp)=0.04+0.023pT/(GeV/c) [33], with a weak dependence on the rapidity and colliding energy√sNN . In nuclear collisions, bottom quarks suffer strong energy losses in the thermal medium. B hadrons and non-prompt charmonia are shifted from highpT to relatively lowpT . This hot medium modification of the non-prompt charmonium (or bottom quark) momentum distribution is characterized with a quench factorRQ . In highpT ,RQ is smaller than unity, and is calculated as 0.4 for the non-promptJ/ψ RAA [33]. This value is employed in the entirepT region. With both prompt and non-prompt charmonia, one can obtain charmonium's inclusive nuclear modification factors, which are shown in Fig.2. Taking into account the large uncertainties ofdσcˉcpp/dy in the transport model, this paper performs two calculations forRAA under a change indσcˉcpp/dy of ±20% (see the color band in Fig.2) . In most central collisions, primordial production is strongly suppressed and regeneration dominates the total yield.Figure 2. (color online) (left-hand panel) Inclusive nuclear modification factor
RAA ofJ/ψ as a function of the number of participantsNp in central rapidity for Pb-Pb collisions at√sNN=5.02 TeV. The dotted and dashed lines represent primordial production and regeneration, respectively. The solid black line in the middle of the band is forJ/ψ inclusive ofRAA . The color band is for the uncertainties of the inputs, withdσcˉcpp/dy changed by ±20%. Experimental data is from the ALICE Collaboration [47]. (right-hand panel) InclusiveRAA in forward rapidity. The lines and band are similar to the left panel. Experimental data is from ALICE Collaboration [48].In the forward rapidity, both initial conditions of the hydrodynamic equations and transport model are updated by comparison to the central rapidity collisions. In the right panel of Fig.2, the
J/ψ nuclear modification factor from primordial production (dotted line), regeneration (dashed line), and inclusive production (color band) are plotted separately. The flat tendency of the experimental data at largerNp is due to the combined effects of the decreased primordial production and increased regeneration in the finalJ/ψ yield. The experimental data in the right-hand panel of Fig.2 is plotted for the range0.3<pT<8 GeV/c. This excludes the contribution from coherent photoproduction, which occurs below 0.3 GeV/c and can make the totalRAA larger than unity in ultra-peripheral collisions [31].In order to show the contributions of primordial production and regeneration, the
pT -differentialRAA is also plotted in the left-hand panel of Fig.3. The significant increase ofRAA in the lowpT region is caused by regeneration, and its large suppression in the highpT region is due to the effects of color screening and the inelastic collisions of partons. The dotted line represents the initial production, and it increases slightly withpT due to the leakage effect. The theoretical results forRJ/ψAA in both central and forward rapidities can explain the experimental data well. Note that even the charm pair cross-sectiondσcˉcpp/dy value is larger in a central rapidity than in a forward one,RAA is similar across the two rapidities, because in a central rapidity with a hotter medium, the QGP strong expansion “blows” charm quarks to a larger volume, which suppresses the charm quark spatial density and the regeneration rate of charmonium. Meanwhile, the elliptic flows of regenerated charmonia become larger in the central rapidity. These will be discussed in detail below. In the right-hand panel of Fig.3, the contributions of promptJ/ψ andψ(2S) toRAA , as a function of rapidity, are also presented.Figure 3. (color online) (left-hand panel)
J/ψ nuclear modification factorRAA as a function of the transverse momentumpT . The dotted line represents initial production, and the solid line is the inclusive production comprising initial production, regeneration and B hadron decay. The color band is due to the uncertainties ofdσcˉcpp/dy (see Fig.2). The differences between the solid and dotted lines are mainly due to the contributions of regeneration at lowpT and B hadron decay at highpT , respectively. Experimental data is from the ALICE Collaboration [49]. (right-hand panel) The contribution of promptJ/ψ towardsRAA as a function of rapidity in the centrality 0%-100%. Promptψ(2S) s contribution toRAA(y) is also predicted. The experimental data is from [50].The situation is more complicated for
ψ(2S) in the hot medium, owing to its dissociation rate compared with the tightly boundJ/ψ . Theψ(2S) decay rate here is extracted fromJ/ψ ' using the nuclear geometry scale of Sec. II. In Fig.4, charmonia at highpT mainly arise from the primordial production. In moving from peripheral to central collisions, the dissociation rate ofψ(2S) increases, and the ratio ofRψ(2S)AA/RJ/ψAA decreases withNp . In the most central collisions at LHC, the ratio of charmonium nuclear modification factors is proportional to their decay rates. In peripheral collisions, the charmonium path length in the hot medium becomes smaller. Under weak suppression, theJ/ψ andψ(2S) nuclear modification factors approach unity, which entails thatRψ(2S)AA/RJ/ψAA→1 asNp→2 , (see left-hand panel of Fig.4). The individualRJ/ψAA in the highpT bin is also plotted in the right-hand panel of Fig.4. The decreasing tendencies ofJ/ψ andψ(2S) RAA with increasingNp are explained well.Figure 4. (color online) (left-hand panel) The ratio of
J/ψ andψ(2S) prompt nuclear modification factors in the central rapidity region with momentum cut-offs of6.5<pT<30 GeV/c for√sNN=5.02 TeV Pb-Pb collisions. Experimental data is from the CMS Collaboration [26]. (right-hand panel)J/ψ andψ(2S) prompt nuclear modification factors as a function ofNp in the central rapidity region with momentum cut-offs of6.5<pT<30 GeV/c for√sNN=5.02 TeV Pb-Pb collisions. The experimental data is from [50].The transverse momentum dependence of
Rψ(2S)AA/RJ/ψAA is also studied via the sequential regeneration mechanism. In Fig.5, aspT→0 , there will be significant regeneration ofJ/ψ . The loosely-boundψ(2S) can only be thermally produced in the later stages of QGP expansion compared withJ/ψ , which leads to the production of regeneratedψ(2S) with larger velocities and collective flows inherited from the bulk medium, because of the strong couplings between the charm quarks and the deconfined medium. Therefore, the regeneratedψ(2S) are distributed at a largerpT (dashed line in right-hand panel of Fig.5) compared with the regeneratedJ/ψ . Due to the differentpT distributions of the thermally producedJ/ψ andψ(2S) , the shapes ofJ/ψ andψ(2S) RAA(pT) s (solid black and blue lines in right-hand panel of Fig.5) are different. There is a "peak" in the ratioRψ(2S)AA/RJ/ψAA in the left-hand panel of Fig.5 due to the sequential regeneration ofψ(2S) . In the absence ofψ(2S) regeneration, the ratio will decrease to zero aspT→0 (see the solid line in the left-hand panel of Fig.5.Figure 5. (color online) (left-hand panel) The ratio of the
J/ψ andψ(2S) prompt nuclear modification factors as a function of the transverse momentumpT in the minimum bias (corresponding to the impact parameter b = 8.4 fm) for√sNN=5.02 TeV Pb-Pb collisions. Dotted-dashed and solid lines are with and withoutψ(2S) regeneration, respectively. Experimental data is from the CMS Collaboration [26, 50]. (Right-hand panel) TheJ/ψ andψ(2S) nuclear modification factors as a function of transverse momentum. Dotted, dashed and solid black lines are for initial, regenerated and promptψ(2S) , respectively. The promptRJ/ψAA is also plotted for comparison. The experimental data is from [26, 50].In the low and middle
pT regions, regeneration becomes important forJ/ψ andψ(2S) , and enhances theirRAA in semi-central and central collisions. In order to show the role ofψ(2S) regeneration onRψ(2S)AA/RJ/ψAA , this study presents two calculations with and withoutψ(2S) regeneration, respectively, in Fig.6. In Fig.6, neglecting the regeneration ofψ(2S) (solid line),Rψ(2S)AA/RJ/ψAA continues to decrease withNp due to the stronger QGP suppression of charmonium excited states. The contribution of regeneration increasesψ(2S) production, most notably in the central collisions, and enhances the value ofRψ(2S)AA/RJ/ψAA (see the dotted-dashed line). Note that in Ref. [33], predictions about promptRψ(2S)AA/RJ/ψAA in 2.76 TeV Pb-Pb collisions have been made. Its value is predicted to be around∼0.15 in allpT bins. The binding energy and regeneration rate forψ(2S) in Ref. [33] are smaller, which suppresses theψ(2S) production. In this study these calculations are extended from 2.76 TeV [33] to 5.02 TeV, and both calculations are consistent with the experimental data at 5.02 TeV. The difference between the solid and dotted-dashed lines in Fig.6 is due to theψ(2S) regeneration component.Figure 6. (color online) Ratio of the
J/ψ andψ(2S) prompt nuclear modification factors as a function ofNp in rapidity1.6<|y|<2.4 with momentum cut-offs of3<pT<30 GeV/c. The dotted-dashed line represents the scenario forJ/ψ andψ(2S) with both primordial production and regeneration, and the solid line is the scenario without regeneration forψ(2S) (J/ψ regeneration is included in both lines). Experimental data at√sNN=2.76 TeV and 5.02 TeV is from the CMS Collaboration [26].Furthermore, the anisotropies of the
J/ψ andψ(2S) momentum distributions are shown in Fig.7. Charmonia moving inside the QGP are likely to be in a color-neutral bound state, thus they are weakly coupled with the bulk medium, and less affected by the collective expansion of the QGP compared with unbound charm quarks. The non-zero elliptic flow of primordially-producedJ/ψ atpT>6 GeV/c is mainly due to the effects of path length differences in the transverse plane (see the solid line in the left-hand panel of Fig.7). At lowpT ,J/ψ production is dominated by the regeneration component. These heavy quarks are strongly coupled with the thermal medium and inherit collective flows, which results in a peak ofv2 atpT∼3 GeV/c. The elliptic flows of the promptψ(2S) are also displayed as a dotted line. Asψ(2S) are regenerated in the later stages of QGP anisotropic expansion, their elliptic flow is larger than that ofJ/ψ '. In the highpT region, the momentum anisotropy ofψ(2S) is larger than that ofJ/ψ ', as they are easily dissociated and sensitive to the anisotropy of the bulk medium.Figure 7. (color online) (left-hand panel) The elliptic flows of the prompt
J/ψ andψ(2S) as a function of the transverse momentumpT in centrality 20%-40% for the forward rapidity2.5<|y|<4 √sNN=5.02 TeV Pb-Pb collisions. Solid and dotted lines are used for the promptJ/ψ andψ(2S) , respectively. (right-hand panel) The elliptic flows of inclusiveJ/ψ as a function of the transverse momentumpT . The solid and dotted lines represent the inclusiveJ/ψ in forward and central rapidities, respectively. Experimental data is from the ALICE Collaboration [51].In the right-hand panel of Fig.7 the experimental data is for the inclusive
J/ψ , and includes a non-prompt contribution from the B hadron decay. The solid line represents the inclusiveJ/ψ , assuming the kinetic equilibrium for bottom quarks as the upper limit [30]. The non-prompt contribution becomes important at highpT and therefore the kinetically thermalized bottom quarks can enhance the inclusivevJ/ψ2 by∼2 % atpT∼8 GeV/c.vJ/ψ2 in central rapidity is also calculated and shown as a dotted line for comparison. The situation for inclusiveψ(2S) is connected with the energy loss of bottom quarks in QGP, and has been comprehensively studied in Ref. [33]. -
In summary, this work employs an improved transport model to study the thermal production of
J/ψ andψ(2S) in Pb-Pb collisions at√sNN=5.02 TeV. Charmonium nuclear modification factors are dominated by regeneration at lowpT and primordial production at highpT , respectively. With different binding energies,J/ψ andψ(2S) can be sequentially produced in the different stages of QGP anisotropic expansion, resulting in differentpT distributions of the regeneratedJ/ψ andψ(2S) . This explains well both theJ/ψ andψ(2S) RAA s and their ratio in 5.02 TeV Pb-Pb collisions, and clearly shows how the open charm quark evolution can affect the final charmonium production. The sequential regeneration ofJ/ψ andψ(2S) describes the charm quark diffusion histories and QGP expansion in heavy-ion collisions.I acknowledge instructive discussions with Prof. Pengfei Zhuang and Jiaxing Zhao. I am also grateful to Prof. Carsten Greiner for the kind hospitality during this study.
