Loading [MathJax]/jax/output/HTML-CSS/jax.js

Thermal production of charmonia in Pb-Pb collisions at sNN=5.02 TeV

  • This work uses the Boltzmann transport model to study the thermal production of J/ψ and ψ(2S) in the quark gluon plasma (QGP) produced by sNN=5.02 TeV Pb-Pb collisions. The J/ψ nuclear modification factors are studied in detail alongside the mechanisms of primordial production and the recombination of charm and anti-charm quarks in the thermal medium. The ψ(2S) binding energy is much smaller in the hot medium compared to the ground state; thus, ψ(2S) with middle to low pT can be thermally regenerated in the later stages of QGP expansions, enabling ψ(2S) to inherit larger collective flows from the bulk medium. We quantitatively study the nuclear modification factors of both J/ψ and ψ(2S) in different centralities and transverse momentum bins for sNN=5.02 TeV Pb-Pb collisions.
  • 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 [25] and theoretical models [612]. 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) [1420]. 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 new J/ψs in the QGP [2123], in a process called “regeneration”. As most charm quarks are distributed in low and middle pT regions, the regeneration process dominates the nuclear modification factor and the collective flows of J/ψ in the low and middle pT bins [24]. In the high pT 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 ratio Rψ(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 when 3<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 to J/ψ, ψ(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 than J/ψ, and inherits larger collective flows from the bulk medium [17, 27, 28]. This sequential regeneration can affect the pT dependence of the ratio Rψ(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 temperature Td 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 for J/ψ 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η)τη+vTT]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 αΨ,

    αΨ=1ETd3k(2π)3EgσgΨ(p,k,T)FgΨ(p,k)fg(k,T),

    (2)

    where Eg and fg 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, the J/ψ decay rate αΨ is compared with the quasi-free dissociation rate [27]. Most J/ψs can survive in the relatively low temperature region, and two different transport models of the decay rate present similar final results for the region T<300 MeV, where most QGP and charmonia are located [7, 30].

    Figure 1

    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 at T>Td strongly suppresses the ψ(2S) production, where no excited states can survive inside the QGP at T>Tχc,ψ(2S)d1.1Tc. Similarly to fast cooling systems, charmonium states might survive in the region T>Td if the medium quickly cools down to below Td. 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 of J/ψs due to their large Td. The ψ(2S) decay rate is extracted from the nuclear geometry scale, αψ(2S)=αJ/ψ×r2ψ(2S)/r2J/ψ, and χc is found in a similar way. The mean radius of charmonia in a vacuum is calculated with the potential model rJ/ψ,χ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 [3436]. 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(xTb)cosh(η)τ0,

    (5)

    where TA and TB are the thickness functions of two colliding nuclei, with the definition TA(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 region 2.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)=(n1)π(n2)p2Tpp[1+p2T(n2)p2Tpp]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 rapidity 2.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 p2Tpp|y=0=12.5 (GeV/c)2 and n=3.2. For the charmonium momentum distribution in other rapidities, which have less constraints, p2Tpp(y) is determined by the relation,

    p2TJ/ψpp(y)=p2TJ/ψpp|y=0×[1(yymax)2],

    (8)

    where ymax=cosh1(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 to J/ψ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 mass ms=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 of Tc=170 MeV. The initial maximum temperature of the QGP is calculated as T0(xrmT=0,τ0)=510 MeV in the central rapidity |y|<2.4, and 450 MeV in the forward rapidity 2.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+ˉcJ/ψ+g is plotted with the dashed line and is proportional to the number of charm quark pairs in the QGP, it dominates total J/ψ 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 in J/ψ inclusive production is fitted as fB=NBJ/ψpp/(Npromptpp+NBJ/ψ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 high pT to relatively low pT. This hot medium modification of the non-prompt charmonium (or bottom quark) momentum distribution is characterized with a quench factor RQ. In high pT,RQ is smaller than unity, and is calculated as 0.4 for the non-prompt J/ψ RAA [33]. This value is employed in the entire pT 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 of dσcˉcpp/dy in the transport model, this paper performs two calculations for RAA under a change in dσ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

    Figure 2.  (color online) (left-hand panel) Inclusive nuclear modification factor RAA of J/ψ as a function of the number of participants Np 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 for J/ψ inclusive of RAA. The color band is for the uncertainties of the inputs, with dσcˉcpp/dy changed by ±20%. Experimental data is from the ALICE Collaboration [47]. (right-hand panel) Inclusive RAA 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 larger Np is due to the combined effects of the decreased primordial production and increased regeneration in the final J/ψ yield. The experimental data in the right-hand panel of Fig.2 is plotted for the range 0.3<pT<8 GeV/c. This excludes the contribution from coherent photoproduction, which occurs below 0.3 GeV/c and can make the total RAA larger than unity in ultra-peripheral collisions [31].

    In order to show the contributions of primordial production and regeneration, the pT-differential RAA is also plotted in the left-hand panel of Fig.3. The significant increase of RAA in the low pT region is caused by regeneration, and its large suppression in the high pT 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 with pT due to the leakage effect. The theoretical results for RJ/ψAA in both central and forward rapidities can explain the experimental data well. Note that even the charm pair cross-section dσ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 prompt J/ψ and ψ(2S) to RAA, as a function of rapidity, are also presented.

    Figure 3

    Figure 3.  (color online) (left-hand panel) J/ψ nuclear modification factor RAA as a function of the transverse momentum pT. 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 of dσcˉcpp/dy (see Fig.2). The differences between the solid and dotted lines are mainly due to the contributions of regeneration at low pT and B hadron decay at high pT, respectively. Experimental data is from the ALICE Collaboration [49]. (right-hand panel) The contribution of prompt J/ψ towards RAA as a function of rapidity in the centrality 0%-100%. Prompt ψ(2S)s contribution to RAA(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 bound J/ψ. The ψ(2S) decay rate here is extracted from J/ψ' using the nuclear geometry scale of Sec. II. In Fig.4, charmonia at high pTmainly arise from the primordial production. In moving from peripheral to central collisions, the dissociation rate of ψ(2S) increases, and the ratio of Rψ(2S)AA/RJ/ψAA decreases with Np. 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, the J/ψ and ψ(2S) nuclear modification factors approach unity, which entails that Rψ(2S)AA/RJ/ψAA1 as Np2, (see left-hand panel of Fig.4). The individual RJ/ψAA in the high pT bin is also plotted in the right-hand panel of Fig.4. The decreasing tendencies of J/ψ and ψ(2S) RAA with increasing Np are explained well.

    Figure 4

    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 of 6.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 of Np in the central rapidity region with momentum cut-offs of 6.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, as pT0, there will be significant regeneration of J/ψ. The loosely-bound ψ(2S)can only be thermally produced in the later stages of QGP expansion compared with J/ψ, 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 larger pT (dashed line in right-hand panel of Fig.5) compared with the regenerated J/ψ. Due to the different pT distributions of the thermally produced J/ψ and ψ(2S), the shapes of J/ψ 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 ratio Rψ(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 as pT0 (see the solid line in the left-hand panel of Fig.5.

    Figure 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 momentum pT 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) The J/ψ 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 prompt RJ/ψAA is also plotted for comparison. The experimental data is from [26, 50].

    In the low and middle pT regions, regeneration becomes important for J/ψ and ψ(2S), and enhances their RAA in semi-central and central collisions. In order to show the role of ψ(2S) regeneration on Rψ(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 with Np 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 of Rψ(2S)AA/RJ/ψAA (see the dotted-dashed line). Note that in Ref. [33], predictions about prompt Rψ(2S)AA/RJ/ψAA in 2.76 TeV Pb-Pb collisions have been made. Its value is predicted to be around 0.15 in all pT 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

    Figure 6.  (color online) Ratio of the J/ψ and ψ(2S) prompt nuclear modification factors as a function of Np in rapidity 1.6<|y|<2.4 with momentum cut-offs of3<pT<30 GeV/c. The dotted-dashed line represents the scenario for J/ψ 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-produced J/ψ at pT>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 low pT, 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 of v2 at pT3 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 of J/ψ'. In the high pT region, the momentum anisotropy of ψ(2S) is larger than that of J/ψ', as they are easily dissociated and sensitive to the anisotropy of the bulk medium.

    Figure 7

    Figure 7.  (color online) (left-hand panel) The elliptic flows of the prompt J/ψ and ψ(2S) as a function of the transverse momentum pT in centrality 20%-40% for the forward rapidity 2.5<|y|<4 sNN=5.02 TeV Pb-Pb collisions. Solid and dotted lines are used for the prompt J/ψ and ψ(2S), respectively. (right-hand panel) The elliptic flows of inclusive J/ψ as a function of the transverse momentum pT. The solid and dotted lines represent the inclusive J/ψ 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 inclusive J/ψ, assuming the kinetic equilibrium for bottom quarks as the upper limit [30]. The non-prompt contribution becomes important at high pT and therefore the kinetically thermalized bottom quarks can enhance the inclusive vJ/ψ2 by 2% at pT8 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 low pT and primordial production at high pT, respectively. With different binding energies, J/ψ and ψ(2S) can be sequentially produced in the different stages of QGP anisotropic expansion, resulting in different pT distributions of the regenerated J/ψ and ψ(2S). This explains well both theJ/ψ and ψ(2S) RAAs 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 of J/ψ 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.

    [1] T. Matsui and H. Satz, Phys. Lett. B, 178: 416 (1986) doi: 10.1016/0370-2693(86)91404-8
    [2] A. Adare et al [PHENIX Collaboration], Phys. Rev. Lett., 98: 232301 (2007) doi: 10.1103/PhysRevLett.98.232301
    [3] B. Abelev et al [ALICE Collaboration], Phys. Rev. Lett., 109: 072301 (2012) doi: 10.1103/PhysRevLett.109.072301
    [4] S. Chatrchyan et al [CMS Collaboration], JHEP, 1205: 063 (2012)
    [5] J. Adam et al [ALICE Collaboration], JHEP, 1511: 127 (2015)
    [6] L. Grandchamp and R. Rapp, Phys. Lett. B, 523: 60 (2001) doi: 10.1016/S0370-2693(01)01311-9
    [7] X. Zhao and R. Rapp, Phys. Rev. C, 82: 064905 (2010) doi: 10.1103/PhysRevC.82.064905
    [8] X. Du and R. Rapp, arXiv: 1808.10014[nucl-th]
    [9] X. l. Zhu, P. f. Zhuang, and N. Xu, Phys. Lett. B, 607: 107 (2005) doi: 10.1016/j.physletb.2004.12.023
    [10] B. Chen, Phys. Rev. C, 93(5): 054905(2016); B. Chen, K. Zhou, and P. Zhuang, Phys. Rev. C, 86: 034906(2012)
    [11] J. P. Blaizot, D. De Boni, P. Faccioli et al, Nucl. Phys. A, 946: 49 (2016) doi: 10.1016/j.nuclphysa.2015.10.011
    [12] J. P. Blaizot and M. A. Escobedo, Phys. Rev. D, 98(7): 074007 (2018) doi: 10.1103/PhysRevD.98.074007
    [13] C. Gerschel and J. Hufner, Phys. Lett. B, 207: 253 (1988) doi: 10.1016/0370-2693(88)90570-9
    [14] H. Satz, J. Phys. G, 32: R25 (2006) doi: 10.1088/0954-3899/32/3/R01
    [15] R. Katz and P. B. Gossiaux, Annals Phys., 368: 267 (2016) doi: 10.1016/j.aop.2016.02.005
    [16] J. P. Blaizot and M. A. Escobedo, JHEP, 1806: 034 (2018)
    [17] B. Chen, Phys. Rev. C, 95(3): 034908 (2017) doi: 10.1103/PhysRevC.95.034908
    [18] X. Yao and B. Müller, Phys. Rev. C, 97(1): 014908(2018) Erratum: [Phys. Rev. C, 97: 049903(E) (2018)]
    [19] X. Yao and T. Mehen, arXiv: 1811.07027[hep-ph]
    [20] X. Yao and B. Müller, arXiv: 1811.09644[hep-ph]
    [21] R. L. Thews, M. Schroedter, and J. Rafelski, Phys. Rev. C, 63: 054905 (2001) doi: 10.1103/PhysRevC.63.054905
    [22] A. Andronic, P. Braun-Munzinger, K. Redlich et al, Phys. Lett. B, 571: 36 (2003) doi: 10.1016/j.physletb.2003.07.066
    [23] L. Yan, P. Zhuang, and N. Xu, Phys. Rev. Lett., 97: 232301 (2006) doi: 10.1103/PhysRevLett.97.232301
    [24] X. Zhao and R. Rapp, Nucl. Phys. A, 859: 114 (2011) doi: 10.1016/j.nuclphysa.2011.05.001
    [25] V. Khachatryan et al [CMS Collaboration], Phys. Rev. Lett., 113(26): 262301 (2014) doi: 10.1103/PhysRevLett.113.262301
    [26] A. M. Sirunyan et al [CMS Collaboration], Phys. Rev. Lett., 118(16): 162301 (2017) doi: 10.1103/PhysRevLett.118.162301
    [27] X. Du and R. Rapp, Nucl. Phys. A, 943: 147 (2015) doi: 10.1016/j.nuclphysa.2015.09.006
    [28] J. Zhao and B. Chen, Phys. Lett. B, 776: 17 (2018) doi: 10.1016/j.physletb.2017.11.014
    [29] B. Chen, T. Guo, Y. Liu et al, Phys. Lett. B, 765: 323 (2017) doi: 10.1016/j.physletb.2016.12.021
    [30] K. Zhou, N. Xu, Z. Xu et al, Phys. Rev. C, 89(5): 054911 (2014) doi: 10.1103/PhysRevC.89.054911
    [31] W. Shi, W. Zha, and B. Chen, Phys. Lett. B, 777: 399 (2018) doi: 10.1016/j.physletb.2017.12.055
    [32] Y. Burnier, O. Kaczmarek, and A. Rothkopf, Phys. Rev. Lett., 114(8): 082001 (2015) doi: 10.1103/PhysRevLett.114.082001
    [33] B. Chen, Y. Liu, K. Zhou et al, Phys. Lett. B, 726: 725 (2013) doi: 10.1016/j.physletb.2013.09.036
    [34] B. I. Abelev et al [STAR Collaboration], Phys. Rev. Lett., 98: 192301(2007) Erratum: [Phys. Rev. Lett., 106: 159902(2011)]
    [35] A. Adare et al [PHENIX Collaboration], Phys. Rev. Lett., 98: 172301 (2007) doi: 10.1103/PhysRevLett.98.172301
    [36] B. Abelev et al [ALICE Collaboration], Phys. Rev. Lett., 111: 102301 (2013) doi: 10.1103/PhysRevLett.111.102301
    [37] W. Zhao, H. j. Xu, and H. Song, Eur. Phys. J. C, 77(9): 645 (2017) doi: 10.1140/epjc/s10052-017-5186-x
    [38] J. Zhao, S. Shi, N. Xu et al, arXiv: 1805.10858[hep-ph]
    [39] R. Aaij et al [LHCb Collaboration], JHEP, 1706: 147 (2017)
    [40] B. Abelev et al [ALICE Collaboration], Phys. Lett. B, 718: 295(2012) Erratum: [Phys. Lett. B, 748: 472(2015)]
    [41] B. Abelev et al [ALICE Collaboration], JHEP, 1211: 065 (2012)
    [42] R. Vogt, Phys. Rev. C, 71: 054902 (2005) doi: 10.1103/PhysRevC.71.054902
    [43] K. J. Eskola, V. J. Kolhinen, and C. A. Salgado, Eur. Phys. J. C, 9: 61 (1999) doi: 10.1007/s100529900005
    [44] K. J. Eskola, H. Paukkunen, and C. A. Salgado, JHEP, 0904: 065 (2009)
    [45] J. Sollfrank, P. Huovinen, M. Kataja et al, Phys. Rev. C, 55: 392 (1997) doi: 10.1103/PhysRevC.55.392
    [46] C. Patrignani et al [Particle Data Group], Chin. Phys. C, 40(10): 100001 (2016) doi: 10.1088/1674-1137/40/10/100001
    [47] R. T. Jiménez Bustamante[ALICE Collaboration], Nucl. Phys. A, 967: 576 (2017) doi: 10.1016/j.nuclphysa.2017.06.026
    [48] J. Adam et al [ALICE Collaboration], Phys. Lett. B, 766: 212 (2017) doi: 10.1016/j.physletb.2016.12.064
    [49] B. Paul[ALICE Collaboration], J. Phys. Conf. Ser., 779(1): 012037 (2017)
    [50] A. M. Sirunyan et al [CMS Collaboration], Eur. Phys. J. C, 78(6): 509 (2018) doi: 10.1140/epjc/s10052-018-5950-6
    [51] S. Acharya et al [ALICE Collaboration], Phys. Rev. Lett., 119(24): 242301 (2017) doi: 10.1103/PhysRevLett.119.242301
  • [1] T. Matsui and H. Satz, Phys. Lett. B, 178: 416 (1986) doi: 10.1016/0370-2693(86)91404-8
    [2] A. Adare et al [PHENIX Collaboration], Phys. Rev. Lett., 98: 232301 (2007) doi: 10.1103/PhysRevLett.98.232301
    [3] B. Abelev et al [ALICE Collaboration], Phys. Rev. Lett., 109: 072301 (2012) doi: 10.1103/PhysRevLett.109.072301
    [4] S. Chatrchyan et al [CMS Collaboration], JHEP, 1205: 063 (2012)
    [5] J. Adam et al [ALICE Collaboration], JHEP, 1511: 127 (2015)
    [6] L. Grandchamp and R. Rapp, Phys. Lett. B, 523: 60 (2001) doi: 10.1016/S0370-2693(01)01311-9
    [7] X. Zhao and R. Rapp, Phys. Rev. C, 82: 064905 (2010) doi: 10.1103/PhysRevC.82.064905
    [8] X. Du and R. Rapp, arXiv: 1808.10014[nucl-th]
    [9] X. l. Zhu, P. f. Zhuang, and N. Xu, Phys. Lett. B, 607: 107 (2005) doi: 10.1016/j.physletb.2004.12.023
    [10] B. Chen, Phys. Rev. C, 93(5): 054905(2016); B. Chen, K. Zhou, and P. Zhuang, Phys. Rev. C, 86: 034906(2012)
    [11] J. P. Blaizot, D. De Boni, P. Faccioli et al, Nucl. Phys. A, 946: 49 (2016) doi: 10.1016/j.nuclphysa.2015.10.011
    [12] J. P. Blaizot and M. A. Escobedo, Phys. Rev. D, 98(7): 074007 (2018) doi: 10.1103/PhysRevD.98.074007
    [13] C. Gerschel and J. Hufner, Phys. Lett. B, 207: 253 (1988) doi: 10.1016/0370-2693(88)90570-9
    [14] H. Satz, J. Phys. G, 32: R25 (2006) doi: 10.1088/0954-3899/32/3/R01
    [15] R. Katz and P. B. Gossiaux, Annals Phys., 368: 267 (2016) doi: 10.1016/j.aop.2016.02.005
    [16] J. P. Blaizot and M. A. Escobedo, JHEP, 1806: 034 (2018)
    [17] B. Chen, Phys. Rev. C, 95(3): 034908 (2017) doi: 10.1103/PhysRevC.95.034908
    [18] X. Yao and B. Müller, Phys. Rev. C, 97(1): 014908(2018) Erratum: [Phys. Rev. C, 97: 049903(E) (2018)]
    [19] X. Yao and T. Mehen, arXiv: 1811.07027[hep-ph]
    [20] X. Yao and B. Müller, arXiv: 1811.09644[hep-ph]
    [21] R. L. Thews, M. Schroedter, and J. Rafelski, Phys. Rev. C, 63: 054905 (2001) doi: 10.1103/PhysRevC.63.054905
    [22] A. Andronic, P. Braun-Munzinger, K. Redlich et al, Phys. Lett. B, 571: 36 (2003) doi: 10.1016/j.physletb.2003.07.066
    [23] L. Yan, P. Zhuang, and N. Xu, Phys. Rev. Lett., 97: 232301 (2006) doi: 10.1103/PhysRevLett.97.232301
    [24] X. Zhao and R. Rapp, Nucl. Phys. A, 859: 114 (2011) doi: 10.1016/j.nuclphysa.2011.05.001
    [25] V. Khachatryan et al [CMS Collaboration], Phys. Rev. Lett., 113(26): 262301 (2014) doi: 10.1103/PhysRevLett.113.262301
    [26] A. M. Sirunyan et al [CMS Collaboration], Phys. Rev. Lett., 118(16): 162301 (2017) doi: 10.1103/PhysRevLett.118.162301
    [27] X. Du and R. Rapp, Nucl. Phys. A, 943: 147 (2015) doi: 10.1016/j.nuclphysa.2015.09.006
    [28] J. Zhao and B. Chen, Phys. Lett. B, 776: 17 (2018) doi: 10.1016/j.physletb.2017.11.014
    [29] B. Chen, T. Guo, Y. Liu et al, Phys. Lett. B, 765: 323 (2017) doi: 10.1016/j.physletb.2016.12.021
    [30] K. Zhou, N. Xu, Z. Xu et al, Phys. Rev. C, 89(5): 054911 (2014) doi: 10.1103/PhysRevC.89.054911
    [31] W. Shi, W. Zha, and B. Chen, Phys. Lett. B, 777: 399 (2018) doi: 10.1016/j.physletb.2017.12.055
    [32] Y. Burnier, O. Kaczmarek, and A. Rothkopf, Phys. Rev. Lett., 114(8): 082001 (2015) doi: 10.1103/PhysRevLett.114.082001
    [33] B. Chen, Y. Liu, K. Zhou et al, Phys. Lett. B, 726: 725 (2013) doi: 10.1016/j.physletb.2013.09.036
    [34] B. I. Abelev et al [STAR Collaboration], Phys. Rev. Lett., 98: 192301(2007) Erratum: [Phys. Rev. Lett., 106: 159902(2011)]
    [35] A. Adare et al [PHENIX Collaboration], Phys. Rev. Lett., 98: 172301 (2007) doi: 10.1103/PhysRevLett.98.172301
    [36] B. Abelev et al [ALICE Collaboration], Phys. Rev. Lett., 111: 102301 (2013) doi: 10.1103/PhysRevLett.111.102301
    [37] W. Zhao, H. j. Xu, and H. Song, Eur. Phys. J. C, 77(9): 645 (2017) doi: 10.1140/epjc/s10052-017-5186-x
    [38] J. Zhao, S. Shi, N. Xu et al, arXiv: 1805.10858[hep-ph]
    [39] R. Aaij et al [LHCb Collaboration], JHEP, 1706: 147 (2017)
    [40] B. Abelev et al [ALICE Collaboration], Phys. Lett. B, 718: 295(2012) Erratum: [Phys. Lett. B, 748: 472(2015)]
    [41] B. Abelev et al [ALICE Collaboration], JHEP, 1211: 065 (2012)
    [42] R. Vogt, Phys. Rev. C, 71: 054902 (2005) doi: 10.1103/PhysRevC.71.054902
    [43] K. J. Eskola, V. J. Kolhinen, and C. A. Salgado, Eur. Phys. J. C, 9: 61 (1999) doi: 10.1007/s100529900005
    [44] K. J. Eskola, H. Paukkunen, and C. A. Salgado, JHEP, 0904: 065 (2009)
    [45] J. Sollfrank, P. Huovinen, M. Kataja et al, Phys. Rev. C, 55: 392 (1997) doi: 10.1103/PhysRevC.55.392
    [46] C. Patrignani et al [Particle Data Group], Chin. Phys. C, 40(10): 100001 (2016) doi: 10.1088/1674-1137/40/10/100001
    [47] R. T. Jiménez Bustamante[ALICE Collaboration], Nucl. Phys. A, 967: 576 (2017) doi: 10.1016/j.nuclphysa.2017.06.026
    [48] J. Adam et al [ALICE Collaboration], Phys. Lett. B, 766: 212 (2017) doi: 10.1016/j.physletb.2016.12.064
    [49] B. Paul[ALICE Collaboration], J. Phys. Conf. Ser., 779(1): 012037 (2017)
    [50] A. M. Sirunyan et al [CMS Collaboration], Eur. Phys. J. C, 78(6): 509 (2018) doi: 10.1140/epjc/s10052-018-5950-6
    [51] S. Acharya et al [ALICE Collaboration], Phys. Rev. Lett., 119(24): 242301 (2017) doi: 10.1103/PhysRevLett.119.242301
  • 加载中

Cited by

1. Zhou, A., Zheng, S., Chen, B. Canonical Effect on Quarkonium Enhancement in the Deconfined Medium[J]. Symmetry, 2025, 17(6): 830. doi: 10.3390/sym17060830
2. Zheng, S., Liu, J., Zheng, S. et al. Charmonia production in hot QCD matter and electromagnetic fields[J]. Physical Review C, 2025, 111(5): 054908. doi: 10.1103/PhysRevC.111.054908
3. Zhang, L., Li, X., Xing, W.-J. et al. Nuclear modification of Bc mesons in relativistic heavy-ion collisions based on a linear Boltzmann transport model[J]. Physical Review C, 2025, 111(5): 054915. doi: 10.1103/PhysRevC.111.054915
4. Coquet, M.. Investigating the interplay between initial hard processes and final-state effects measuring prompt and non-prompt J/ψ[J]. EPJ Web of Conferences, 2025. doi: 10.1051/epjconf/202531604009
5. Bai, Y., Chen, B. Probing QGP droplets with charmonium in high-multiplicity proton–proton collisions[J]. European Physical Journal C, 2024, 84(11): 1193. doi: 10.1140/epjc/s10052-024-13599-4
6. Zheng, S., Liu, J. Review of Deep Learning in High-Energy Heavy-Ion Collisions[J]. Symmetry, 2024, 16(11): 1426. doi: 10.3390/sym16111426
7. Chen, G., Chen, B., Zhao, J. Bottomonium evolution with in-medium heavy quark potential from lattice QCD[J]. European Physical Journal C, 2024, 84(8): 869. doi: 10.1140/epjc/s10052-024-13266-8
8. Zhao, J., Chen, B. J/ψ polarization in relativistic heavy ion collisions[J]. European Physical Journal C, 2024, 84(8): 875. doi: 10.1140/epjc/s10052-024-13264-w
9. Bai, X.. Recent charmonium measurements in Pb-Pb collisions with ALICE[J]. EPJ Web of Conferences, 2024. doi: 10.1051/epjconf/202429609008
10. Andronic, A., Gossiaux, P.B., Petreczky, P. et al. Comparative study of quarkonium transport in hot QCD matter[J]. European Physical Journal A, 2024, 60(4): 88. doi: 10.1140/epja/s10050-024-01306-6
11. Acharya, S., Adamová, D., Adler, A. et al. Measurements of inclusive J/ψ production at midrapidity and forward rapidity in Pb–Pb collisions at sNN = 5.02 TeV[J]. Physics Letters Section B Nuclear Elementary Particle and High Energy Physics, 2024. doi: 10.1016/j.physletb.2024.138451
12. Acharya, S., Adamová, D., Aglieri Rinella, G. et al. Prompt and non-prompt J/ψ production at midrapidity in Pb–Pb collisions at √sNN = 5.02 TeV[J]. Journal of High Energy Physics, 2024, 2024(2): 66. doi: 10.1007/JHEP02(2024)066
13. Mattioli, K.R.. Towards Experimental Confirmation of Quarkonia Melting in Quark–Gluon Plasma: A Review of Recent Measurements of Quarkonia Production in Relativistic Heavy-Ion Collisions[J]. Symmetry, 2024, 16(2): 225. doi: 10.3390/sym16020225
14. Shukri, A., Abdulsalam, A., Ilahi, M. et al. Computation of quarkonium dissociation and recombination rates in quark-gluon plasma using Bateman equation[J]. International Journal of Modern Physics E, 2023, 32(12): 2341002. doi: 10.1142/S0218301323410021
15. Pan, C., Zheng, S., Yang, M. et al. Nonthermal distributions of charmed hadrons in relativistic heavy-ion collisions[J]. Physical Review C, 2023, 108(3): 034903. doi: 10.1103/PhysRevC.108.034903
16. Abdulsalam, A., Ilahi, M., Alabbasi, L. et al. Dissociation rates and recombination likelihood of bottomonium states in heavy-ion collisions at s NN = 5. 0 2 TeV[J]. International Journal of Modern Physics A, 2023, 38(24): 2350104. doi: 10.1142/S0217751X2350104X
17. Chen, B., Tong, B. Bottomonium transport in p -Pb and Pb-Pb collisions at energies available at the CERN Large Hadron Collider[J]. Physical Review C, 2023, 108(1): 014901. doi: 10.1103/PhysRevC.108.014901
18. Zhao, J., Shi, S. Detecting nuclear mass distribution in isobar collisions via charmonium[J]. European Physical Journal C, 2023, 83(6): 511. doi: 10.1140/epjc/s10052-023-11657-x
19. Yang, M., Zheng, S., Tong, B. et al. Bottom energy loss and nonprompt J/ψ production in relativistic heavy ion collisions[J]. Physical Review C, 2023, 107(5): 054917. doi: 10.1103/PhysRevC.107.054917
20. Wen, L., Chen, B. Bottomonium sequential suppression and strong heavy-quark potential in heavy-ion collisions[J]. Physics Letters Section B Nuclear Elementary Particle and High Energy Physics, 2023. doi: 10.1016/j.physletb.2023.137774
21. Zhao, J., Chen, B. Charmonium transport in the high- μb medium[J]. Physical Review C, 2023, 107(4): 044909. doi: 10.1103/PhysRevC.107.044909
22. Wen, L., Du, X., Shi, S. et al. Investigating color screening in proton-nucleus collisions with complex potentials*[J]. Chinese Physics C, 2022, 46(11): 114102. doi: 10.1088/1674-1137/ac7fe6
23. Li, C., Chen, B. Anisotropic Flows of Charmonium in the Relativistic Heavy-Ion Collisions[J]. Mathematics, 2022, 10(22): 4284. doi: 10.3390/math10224284
24. Tong, B., Chen, B. Charmonium dissociation at high baryon chemical potential[J]. Physical Review C, 2022, 106(3): 034911. doi: 10.1103/PhysRevC.106.034911
25. Zhao, J., Zhuang, P. Effects of cold and hot nuclear matter on J/ψ production at energies selected for the beam energy scan at the BNL Relativistic Heavy Ion Collider[J]. Physical Review C, 2022, 105(6): 064907. doi: 10.1103/PhysRevC.105.064907
26. Chen, B., Jiang, L., Liu, X.-H. et al. X(3872) production in relativistic heavy-ion collisions[J]. Physical Review C, 2022, 105(5): 054901. doi: 10.1103/PhysRevC.105.054901
27. Zhao, J., Chen, B., Zhuang, P. Charmonium triangular flow in high energy nuclear collisions[J]. Physical Review C, 2022, 105(3): 034902. doi: 10.1103/PhysRevC.105.034902
28. Xie, Z., Han, R. Flavor-dependent photoproduction in heavy-ion collisions[J]. Communications in Theoretical Physics, 2021, 73(12): 125301. doi: 10.1088/1572-9494/ac3231
29. Chen, B., Jiang, L., Liu, Y. Hadronic effects on charmonium elliptic flows in heavy-ion collisions[J]. Physical Review C, 2021, 103(6): 064910. doi: 10.1103/PhysRevC.103.064910
30. Abdulsalam, A.. The rates of charmonium dissociation and recombination in heavy-ion collisions at sNN=5.02 TeV using Bateman equation[J]. Nuclear Physics A, 2021. doi: 10.1016/j.nuclphysa.2020.122130
31. Chen, B., Hu, M., Zhang, H. et al. Probe the tilted Quark-Gluon Plasma with charmonium directed flow[J]. Physics Letters Section B Nuclear Elementary Particle and High Energy Physics, 2020. doi: 10.1016/j.physletb.2020.135271

Figures(7)

Get Citation
Baoyi Chen. Thermal production of charmonia in Pb-Pb collisions at sNN=5.02 TeV[J]. Chinese Physics C, 2019, 43(12): 124101. doi: 10.1088/1674-1137/43/12/124101
Baoyi Chen. Thermal production of charmonia in Pb-Pb collisions at sNN=5.02 TeV[J]. Chinese Physics C, 2019, 43(12): 124101.  doi: 10.1088/1674-1137/43/12/124101 shu
Milestone
Received: 2019-07-09
Article Metric

Article Views(3163)
PDF Downloads(27)
Cited by(31)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Thermal production of charmonia in Pb-Pb collisions at sNN=5.02 TeV

  • 1. Department of Physics, Tianjin University, Tianjin 300350, China
  • 2. Institut für Theoretische Physik, Goethe-Universität Frankfurt, Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany

Abstract: This work uses the Boltzmann transport model to study the thermal production of J/ψ and ψ(2S) in the quark gluon plasma (QGP) produced by sNN=5.02 TeV Pb-Pb collisions. The J/ψ nuclear modification factors are studied in detail alongside the mechanisms of primordial production and the recombination of charm and anti-charm quarks in the thermal medium. The ψ(2S) binding energy is much smaller in the hot medium compared to the ground state; thus, ψ(2S) with middle to low pT can be thermally regenerated in the later stages of QGP expansions, enabling ψ(2S) to inherit larger collective flows from the bulk medium. We quantitatively study the nuclear modification factors of both J/ψ and ψ(2S) in different centralities and transverse momentum bins for sNN=5.02 TeV Pb-Pb collisions.

    HTML

    1.   Introduction
    • 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 [25] and theoretical models [612]. 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) [1420]. 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 new J/ψs in the QGP [2123], in a process called “regeneration”. As most charm quarks are distributed in low and middle pT regions, the regeneration process dominates the nuclear modification factor and the collective flows of J/ψ in the low and middle pT bins [24]. In the high pT 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 ratio Rψ(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 when 3<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 to J/ψ, ψ(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 than J/ψ, and inherits larger collective flows from the bulk medium [17, 27, 28]. This sequential regeneration can affect the pT dependence of the ratio Rψ(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 temperature Td 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 for J/ψ 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.

    2.   Transport model and hydrodynamics
    • 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η)τη+vTT]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 αΨ,

      αΨ=1ETd3k(2π)3EgσgΨ(p,k,T)FgΨ(p,k)fg(k,T),

      (2)

      where Eg and fg 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, the J/ψ decay rate αΨ is compared with the quasi-free dissociation rate [27]. Most J/ψs can survive in the relatively low temperature region, and two different transport models of the decay rate present similar final results for the region T<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 at T>Td strongly suppresses the ψ(2S) production, where no excited states can survive inside the QGP at T>Tχc,ψ(2S)d1.1Tc. Similarly to fast cooling systems, charmonium states might survive in the region T>Td if the medium quickly cools down to below Td. 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 of J/ψs due to their large Td. The ψ(2S) decay rate is extracted from the nuclear geometry scale, αψ(2S)=αJ/ψ×r2ψ(2S)/r2J/ψ, and χc is found in a similar way. The mean radius of charmonia in a vacuum is calculated with the potential model rJ/ψ,χ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 [3436]. 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(xTb)cosh(η)τ0,

      (5)

      where TA and TB are the thickness functions of two colliding nuclei, with the definition TA(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 region 2.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)=(n1)π(n2)p2Tpp[1+p2T(n2)p2Tpp]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 rapidity 2.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 p2Tpp|y=0=12.5 (GeV/c)2 and n=3.2. For the charmonium momentum distribution in other rapidities, which have less constraints, p2Tpp(y) is determined by the relation,

      p2TJ/ψpp(y)=p2TJ/ψpp|y=0×[1(yymax)2],

      (8)

      where ymax=cosh1(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 to J/ψ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 mass ms=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 of Tc=170 MeV. The initial maximum temperature of the QGP is calculated as T0(xrmT=0,τ0)=510 MeV in the central rapidity |y|<2.4, and 450 MeV in the forward rapidity 2.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.

    3.   Numerical results and analysis
    • 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+ˉcJ/ψ+g is plotted with the dashed line and is proportional to the number of charm quark pairs in the QGP, it dominates total J/ψ 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 in J/ψ inclusive production is fitted as fB=NBJ/ψpp/(Npromptpp+NBJ/ψ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 high pT to relatively low pT. This hot medium modification of the non-prompt charmonium (or bottom quark) momentum distribution is characterized with a quench factor RQ. In high pT,RQ is smaller than unity, and is calculated as 0.4 for the non-prompt J/ψ RAA [33]. This value is employed in the entire pT 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 of dσcˉcpp/dy in the transport model, this paper performs two calculations for RAA under a change in dσ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 of J/ψ as a function of the number of participants Np 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 for J/ψ inclusive of RAA. The color band is for the uncertainties of the inputs, with dσcˉcpp/dy changed by ±20%. Experimental data is from the ALICE Collaboration [47]. (right-hand panel) Inclusive RAA 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 larger Np is due to the combined effects of the decreased primordial production and increased regeneration in the final J/ψ yield. The experimental data in the right-hand panel of Fig.2 is plotted for the range 0.3<pT<8 GeV/c. This excludes the contribution from coherent photoproduction, which occurs below 0.3 GeV/c and can make the total RAA larger than unity in ultra-peripheral collisions [31].

      In order to show the contributions of primordial production and regeneration, the pT-differential RAA is also plotted in the left-hand panel of Fig.3. The significant increase of RAA in the low pT region is caused by regeneration, and its large suppression in the high pT 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 with pT due to the leakage effect. The theoretical results for RJ/ψAA in both central and forward rapidities can explain the experimental data well. Note that even the charm pair cross-section dσ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 prompt J/ψ and ψ(2S) to RAA, as a function of rapidity, are also presented.

      Figure 3.  (color online) (left-hand panel) J/ψ nuclear modification factor RAA as a function of the transverse momentum pT. 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 of dσcˉcpp/dy (see Fig.2). The differences between the solid and dotted lines are mainly due to the contributions of regeneration at low pT and B hadron decay at high pT, respectively. Experimental data is from the ALICE Collaboration [49]. (right-hand panel) The contribution of prompt J/ψ towards RAA as a function of rapidity in the centrality 0%-100%. Prompt ψ(2S)s contribution to RAA(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 bound J/ψ. The ψ(2S) decay rate here is extracted from J/ψ' using the nuclear geometry scale of Sec. II. In Fig.4, charmonia at high pTmainly arise from the primordial production. In moving from peripheral to central collisions, the dissociation rate of ψ(2S) increases, and the ratio of Rψ(2S)AA/RJ/ψAA decreases with Np. 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, the J/ψ and ψ(2S) nuclear modification factors approach unity, which entails that Rψ(2S)AA/RJ/ψAA1 as Np2, (see left-hand panel of Fig.4). The individual RJ/ψAA in the high pT bin is also plotted in the right-hand panel of Fig.4. The decreasing tendencies of J/ψ and ψ(2S) RAA with increasing Np 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 of 6.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 of Np in the central rapidity region with momentum cut-offs of 6.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, as pT0, there will be significant regeneration of J/ψ. The loosely-bound ψ(2S)can only be thermally produced in the later stages of QGP expansion compared with J/ψ, 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 larger pT (dashed line in right-hand panel of Fig.5) compared with the regenerated J/ψ. Due to the different pT distributions of the thermally produced J/ψ and ψ(2S), the shapes of J/ψ 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 ratio Rψ(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 as pT0 (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 momentum pT 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) The J/ψ 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 prompt RJ/ψAA is also plotted for comparison. The experimental data is from [26, 50].

      In the low and middle pT regions, regeneration becomes important for J/ψ and ψ(2S), and enhances their RAA in semi-central and central collisions. In order to show the role of ψ(2S) regeneration on Rψ(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 with Np 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 of Rψ(2S)AA/RJ/ψAA (see the dotted-dashed line). Note that in Ref. [33], predictions about prompt Rψ(2S)AA/RJ/ψAA in 2.76 TeV Pb-Pb collisions have been made. Its value is predicted to be around 0.15 in all pT 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 of Np in rapidity 1.6<|y|<2.4 with momentum cut-offs of3<pT<30 GeV/c. The dotted-dashed line represents the scenario for J/ψ 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-produced J/ψ at pT>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 low pT, 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 of v2 at pT3 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 of J/ψ'. In the high pT region, the momentum anisotropy of ψ(2S) is larger than that of J/ψ', 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 momentum pT in centrality 20%-40% for the forward rapidity 2.5<|y|<4 sNN=5.02 TeV Pb-Pb collisions. Solid and dotted lines are used for the prompt J/ψ and ψ(2S), respectively. (right-hand panel) The elliptic flows of inclusive J/ψ as a function of the transverse momentum pT. The solid and dotted lines represent the inclusive J/ψ 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 inclusive J/ψ, assuming the kinetic equilibrium for bottom quarks as the upper limit [30]. The non-prompt contribution becomes important at high pT and therefore the kinetically thermalized bottom quarks can enhance the inclusive vJ/ψ2 by 2% at pT8 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].

    4.   Summary
    • 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 low pT and primordial production at high pT, respectively. With different binding energies, J/ψ and ψ(2S) can be sequentially produced in the different stages of QGP anisotropic expansion, resulting in different pT distributions of the regenerated J/ψ and ψ(2S). This explains well both theJ/ψ and ψ(2S) RAAs 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 of J/ψ 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.

Reference (51)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return