Hydrodynamic description of D meson production in high-energy heavy-ion collisions

  • The large values and constituent-quark-number scaling of the elliptic flow of low-pT D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low-pT charm quarks, we compare the D0 meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. D0 mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D0 mesons at pT<3 GeV/c as measured in Au+Au collisions at sNN=200 GeV. Though the experimental data on D0 spectra are consistent with the hydrodynamic result at small pT1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum,pT>2 GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Our comparative study indicates that charm quarks only approach local thermal equilibrium at small pT, even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate pT.
  • The large collective flow and momentum anisotropy as manifested in the final hadron spectra in high-energy heavy-ion collisions indicate a strong collectivity of the produced dense and hot nuclear matter during its dynamic evolution [1, 2]. The observed approximate constituent quark number (NCQ) scaling of the elliptic flow for light quark hadrons [3-6] and strong jet quenching [7-10] in high-energy heavy-ion collisions at both the Relativistic Heavy-ion Collider (RHIC) and the Large Hadron Collider (LHC) suggest the formation of a hot, deconfined, and opaque quark-gluon plasma (QGP). Similar phenomena are also observed for heavy quark mesons [11-17]. Experimental data from the RHIC and LHC show a large elliptic flow of D0 mesons at low-pT, which also approximately obeys the NCQ scaling of light quark hadrons [13, 18]. Heavy quarks are predominantly produced through initial hard processes and have mass scales much larger than the typical temperature of the QGP medium. Therefore, it still remains an interesting question as to whether and to what degree charm quarks become thermalized [2] and flow with the QGP due to their strong interaction with the hot medium.

    Past studies usually describe the thermalization processes of heavy quarks with transport equations. Their low-momentum (pM) dynamics in the QGP medium are usually treated as Brownian motion [19-21], considering the masses of heavy quarks to be much larger than the typical temperatures in the medium MT. For the transport of heavy quarks with intermediate momentum pM, both the collisional and radiative energy loss have to be considered. For high-momentum heavy quarks, most transport models [22-24] consider the gluon radiations and interactions between heavy quarks and the medium as perturbative processes or use transport coefficients with weakly-coupled assumptions as inputs.

    It was found that transport models need a large heavy-quark momentum-diffusion parameter (ˆq=0.7...2.7 GeV2/fm at T=300 MeV and p=30 GeV/c [25]) to describe the observed medium modification of the heavy-quark meson spectra and elliptic flow coefficients [25, 26]. Several studies have attempted to evaluate the heavy quark diffusion coefficient at zero momentum non-perturbatively from lattice QCD [27-29]. The estimated interaction strength is comparable to the values from phenomenological determination [24, 30] and much larger than the natural expectation of the weakly-coupled theory at the leading order. This poses questions regarding the weak-coupling assumption on the nature of the interaction between heavy quarks and the QGP medium. It is possible that, with a strong coupling, the dynamics of the low-pT charm quarks in the QGP are better described by a hydrodynamic approach.

    In this study, we compute the production of low-pT D mesons in an extreme limit of hydrodynamic evolution with the following assumptions. First, the number of (anti)charm quarks is conserved during the lifetime of the fireball, i.e., thermal production of heavy quark pairs is negligible. Second, the coupling is so strong that the initially produced charm quarks quickly diffuse into the QGP and reach local kinetic equilibrium. Finally, the differential yield of D mesons is computed using the Cooper-Frye formula [31] on the hydrodynamic freeze-out hypersurface, with an effective charm chemical potential to guarantee the charm yield is the same as measured in the experiment. We compare the resultant pT-dependent spectra and elliptic flow of low-pT D mesons, which are considered as the manifestation of the extreme limit of complete thermalization of heavy quarks, to the experimental data as well as transport calculations to quantify the degree of heavy-quark thermalization in heavy-ion collisions. We use the (3+1)D hydrodynamic model (CLVisc) [32] and a linearized Boltzmann-Langevin model [24] for heavy quark transport.

    The paper is organized as follows. In Section II, we introduce the relativistic hydrodynamic model and the Boltzmann-Langevin transport model for the calculations of pT spectra and elliptic flow of charmed mesons. In Section III, we first discuss the calibration of the hydro model to the experimental data on light quark hadron spectra and elliptic flow in Au+Au collisions at the RHIC energy. We then compare the hydro results on the pT spectra and elliptic flow of charmed mesons at low pT to the experimental data and transport calculations and discuss the implications on the degree of thermalization for heavy quarks in the QGP formed in high-energy heavy-ion collisions.

    In the usual paradigm for charmed meson production in high-energy heavy-ion collisions, charm quarks are produced through the initial hard scatterings [33]. These heavy quarks experience energy loss and momentum diffusion in the QGP through both elastic and inelastic collisions with the medium, which can be modeled as the drag and diffusion coefficients in the Boltzmann-Langevin equations. The final D mesons are formed through charm quark fragmentation at high pT [34-36] or charm-light quark recombination at low and intermediate pT [37]. Such transport models for initial heavy quark and final meson production can describe the charm meson spectra and elliptic flow well in high-energy heavy-ion collisions [25, 26]. However, it is still interesting to investigate whether and to what degree the heavy quarks achieve kinematic equilibrium. If these heavy quarks indeed reach local thermal equilibrium, one should expect that the hydrodynamic model can also describe heavy meson spectra and elliptic flow. Furthermore, a hydrodynamic picture overcomes the difficulty of the transport equation when the coupling becomes large between charm quarks or mesons and the medium.

    In this study, we consider the extreme scenario and investigate the charmed meson spectra and elliptic flow in the limit of the fully thermalized low-pT charm quarks. These low-pT charm quarks are still produced initially through hard processes; however, the interactions are assumed to be strong enough that they quickly diffuse into the medium and lose the memory of their initial distribution in phase space, both in terms of the spatial location and momentum of the initial production through hard processes. In the momentum space, they are assumed to reach full kinetic equilibrium and comove with the medium, flowing with the strongly coupled QGP as light quarks and gluons. Finally, the charm quarks transit to charm hadrons with a phase transition of the bulk medium, where we assume the interaction in the hadronic phase is still strong enough to maintain the kinetic equilibrium and equilibrium ratios of different species of charm mesons. In this situation, the low-pT D mesons are produced in the same way as other light hadrons on the freeze-out hypersurface in a relativistic hydrodynamic model. The temperature and fluid velocity profiles on the freeze-out hypersurface are crucial for a reasonable estimate of the D-meson spectra and elliptic flow in the limit of a complete heavy-quark thermalization. We will use the rapidity distribution, pT spectra, and elliptic flow of charged pions to calibrate the hydrodynamic model.

    For completeness, we have to consider D0 mesons from the feed-down of D. Using the feed-down tables from [37, 38],

    D(2007)064.7%D0(1865)+π0,

    (1)

    D(2007)035.3%D0(1865)+γ,

    (2)

    D(2010)+68%D0(1865)+π+

    (3)

    Note that the D's are vector mesons with spin 1. Therefore, the spin degeneracy ratio for D/D is 3:1.

    The hydrodynamic model we use, CLVisc, is a (3+1)D viscous hydrodynamic model parallelized on GPU using OpenCL [32]. The program is well tested against several analytical solutions and can describe the bulk hadron spectra and anisotropic flow in high-energy heavy-ion collisions at top RHIC and LHC energy. It simulates the fluid dynamic evolution of the strongly interacting QCD matter created in high-energy heavy-ion collisions by solving the fluid dynamic equations

    μTμν=0,  with  Tμν=εuμuνPΔμν+πμν,

    (4)

    where ε is the energy density, P is the pressure as a function of energy density given by the equation of state (EoS), Δμν=gμνuμuν is a projection operator, uμ is the fluid four-velocity obeying uμuμ=1 and πμν is the shear stress tensor.

    The initial condition for entropy density distribution in the transverse plane is provided by the Trento Monte Carlo model [39, 40]. An envelope function is used to approximate the longitudinal distribution along the space-time rapidity,

    H(ηs)=exp[(|ηs|ηw)22σ2ηθ(|ηs|ηw)],

    (5)

    where ση=1.5 and σw=1.3 are used for Au+Au collisions at sNN=200 GeV.

    We have assumed an initial time for the hydrodynamics τ0=0.6 fm. In the present study, we use the partial chemical equilibrium EoS with chemical freeze-out temperature 165 MeV and a smooth crossover between a QGP at high temperature and hadron resonance gas (HRG) EoS at low temperature [41] as inspired by the lattice QCD study.

    Baryons and mesons passing through the freeze-out hyper-surface are assumed to obey Fermi-Dirac and Bose-Einstein distributions, respectively. Their momentum distributions are given by the Cooper-Frye formula [31],

    dNidypTdpTdϕ=gi(2π)3pμdΣμf(pu)(1+δf),

    (6)

    where gi=2spin+1 is the spin degeneracy, pμ is the four-momenta of produced particles in the laboratory frame, Σμ is the freeze-out hyper-surface, f(pu) is the Fermi-Dirac/Bose-Einstein distribution function,

    f(pu)=1exp[(puμi)/Tfrz]±1,

    (7)

    and δf is the non-equilibrium correction,

    δf=(1feq)pμpνπμν2T2frz(ε+P).

    (8)

    We have chosen the freeze-out temperature Tfz=137 MeV for light flavor hadrons. The freeze-out temperatures for D mesons are different, and we consider several values to provide an estimate of the uncertainties.

    The elliptic flow of D mesons is defined as the second coefficient of the Fourier decomposition of their azimuthal angle distributions with respect to the event plane of light hadrons,

    d3NpTdpTdydϕ=d2N2πpTdpTdy[1+n=12vncos(n(ϕΨEP))].

    (9)

    We compare the charm meson spectra and flow calculated from the hydrodynamic freeze-out model to that from a transport model [24]. The transport approach assumes that heavy quarks, including those with low momentum in the comoving frame of the medium, remain good quasi-particles in the QGP. Therefore, the dynamics of low-pT heavy flavors can be described by a Boltzmann-type transport equation,

    (t+v)fQ(t,x,p)=[dRdq3(p+q,q)fQ(p+q)dRdq3(p,q)fQ(p)]dq3.

    (10)

    Here, fQ(t,x,p) is the phase-space density of heavy flavors: heavy quarks or heavy mesons (heavy baryons are omitted in this study). dR(p,q)/dq3 is the differential rate for a heavy flavor particle with momentum p to transfer three-momentum q to the local medium with flow velocity u and temperature T.

    The hydrodynamic and the transport model have both overlapping and distinct regimes of application. Both models contain the equilibrium situation. Compared to the hydrodynamic approach, a transport model also governs the far-from-equilibrium dynamics of heavy flavor particles. The traces of off-equilibrium effects can be important in a finite and expanding plasma with moderate coupling between heavy quarks and the medium. In the meantime, hydrodynamics can describe the evolution with large couplings and in the non-quasiparticle regime, which is beyond the applicability of the transport approach. Therefore, both models provide complementary pictures to understand the experimental data.

    The initial charm quark spectrum that initializes the transport equation is obtained from the perturbative QCD based FONLL [42, 43] calculation with EPPS16 nuclear parton distribution function [44]. In the medium above the pseudo-critical temperature (T>Tc), we assume the heavy flavors exist as deconfined heavy quarks. The interaction rate consists of two parts with a comparable contribution to heavy quark energy loss at intermediate pT region: collisional processes and medium-induced gluon radiations.

    Collisional process in the QGP Collisional processes are mainly modeled by two-body scatterings between the heavy quark and medium partons,

    dRi=q,gdq3=¯|M2QiQi|dif(pi)16π2EQEQd3pi2Eid3pi2Ei×δ(4)(pQ+pipQpi),

    (11)

    where Q (i) and Q (i) labels the initial and final state heavy quark (medium quark or gluon). pQ=pQ+q is the three-momentum transfer to the heavy quark. ¯|M2QiQi| is the squared amplitude of the two-body collision at leading order, averaging over the initial-stage quantum number and summed over the final-state quantum number. The ˆt-channel divergence is screened using the QCD non-perturbative scale Λ=0.2 GeV and the Debye screening mass mD=6παsT of a three-flavor plasma,

    ¯|M2QiQi|1ˆt21(ˆtΛ2)(ˆtm2D).

    (12)

    Finally, f(pi)=epiu/T is the classical thermal distribution function of the medium parton in plasma with local four-velocity u and temperature T, obtained from hydrodynamic simulations. The degeneracy factors for medium quarks and gluons are dg=2dA=16 and dq=4NfNc=36, respectively. Besides perturbative scatterings, an additional effective heavy-quark diffusion constant of the form

    ΔκDT3=A+BET

    (13)

    is introduced to mimic a possible non-perturbative contribution that peaks at low energy and low temperature. The effective heavy-quark-medium coupling parameter and the parameters in the non-perturbative diffusion constant are tuned in [24] to the suppression and momentum anisotropy of the production of open heavy flavor particles.

    Medium-induced radiative process in the QGP phase Heavy quarks can radiate additional gluons in the collision with medium partons. This is treated similarly to the rate in Eq. (11) for collisional processes, using two-to-three-body matrix-elements. One complication is that energetic gluon radiation in the rest frame of the medium is suppressed due to the QCD Landau-Pomeranchuk-Migdal (LPM) effect. In [24], this was included by restricting the phase-space of the radiated gluon using an interference factor motivated by the Higher-Twist approach 2[1cos(t/τf)], where t is the time since the last gluon radiation, while τf=2x(1x)EQ/(k2+x2M2) is the gluon formation time. Despite the LPM suppression, the radiative energy loss was found to be equally important in the intermediate pT region.

    Heavy quark hadronization and hadronic rescattering At the pseudo-critical temperature, charm quarks hadronize to D-mesons through both fragmentation and recombination mechanisms [45]. In the hadronic phase, D-mesons continued to interact with the light mesons via D-π and D-ρ scatterings [46] as implemented in Ultra-relativistic Quantum Molecular Dynamics (UrQMD) [47, 48].

    To calibrate the relativistic hydrodynamic model, we have computed the pseudorapidity density of charged particles dNch/dη as shown in Fig. 1, and the pT spectra and the elliptic flow v2 of π+ as shown in Fig. 2 for central Au+Au collisions with centrality range 0-5% at sNN=200 GeV. The calibrated parameters are the scale factor of 53 multiplied to the initial entropy density from the Trento model [39, 40], the starting time for hydro τ0=0.6 fm, the shear viscosity over entropy density ratio η/s=0.15, and the freeze-out temperature for light hadrons Tfz=137 MeV. These values have been tuned such that the predicted dNch/dη, pT spectra and v2(pT) of charged pions from relativistic hydrodynamics agree with the experimental data. Relativistic hydrodynamics with the same set of parameters also describes data well at other centralities and we refer readers to Ref. [32] for a more detailed description of CLVisc and the parameters. These parameters are fixed in the following calculations of the heavy quark meson spectra, except the freeze-out temperatures for D mesons. We will vary the value of the D meson freeze-out temperature and study its effect on the final D meson spectra and the elliptic flow.

    Figure 1

    Figure 1.  (color online) Pseudorapidity density of charged hadrons from integrating the Cooper-Frye formula over the freeze-out hyperspace obtained in the CLVisc calculation. It agrees with the data of the PHOBOS experiment for Au+Au collisions at sNN=200GeV [49].

    Figure 2

    Figure 2.  (color online) Left: transverse momentum spectra of π+ in Au+Au collisions. Right: anisotropic flow coefficient v2 in Au+Au collisions using the event-plane method. In both panels, CLVisc calculations (lines) are compared to data from the PHENIX experiment [50] and STAR experiment [51].

    The effects of freeze-out temperature and resonance decays on D0 meson production are studied in this section. In the left panel of Fig. 3, we show the ratio between D0 pT spectra with and without contributions from the resonance decays. In the right panel of Fig. 3, we compare the D0 meson elliptic flow as a function of pT with (dot-dashed) and without (solid) resonance decays for three different values of the D meson freeze-out temperatures, Tfz=100, 137, and 155 MeV. We see that the resonance decay has a much larger effect on the pT spectra of D0 (left panel) than that on the elliptic flow v2 of D0 (right panel). For a given D meson freeze-out temperature, resonance decays from D contribute more D0 mesons at low pT than that at high pT. The ratio between D0 with and without decay decreases as pT increases for all three freeze-out temperatures considered here. As one increases the freeze-out temperature from 100 MeV to 137 MeV and 155 MeV, not only the D0 yield but also the ratio between D0 with and without decay increases, from approximately 2.35 to 3.0 and 3.3 at pT<1 GeV. The elliptic flow without decay (black solid), in contrast, almost overlaps that with decays (red dash-dotted) for all three different freeze-out temperatures; therefore, it is insensitive to the contributions from the resonance decays. The freeze-out temperature, however, has a large effect on the shape and magnitude of the elliptic flow as a function of transverse momentum. At small transverse momentum (pT < 2 GeV), the elliptic flow increases with the freeze-out temperature. At large transverse momentum (pT > 2 GeV), the trend goes the opposite way – the elliptic flow decreases as the freeze-out temperature increases.

    Figure 3

    Figure 3.  (color online) Left: freeze-out temperature dependence of the effect of resonance decay on the D0 pT spectra. Right: freeze-out temperature dependence of the effect of resonance decay on the D0 flow coefficient v2. Calculations in both panels use one-shot hydrodynamic simulation of 20%-30% central Au+Au collisions at sNN=200 GeV.

    In Fig. 4, we show the hydrodynamic calculations of the transverse momentum (pT) spectra (left) and the pT-dependent elliptic flow (right) of D0 mesons in Au+Au collisions at sNN=200 GeV for different centralities as compared with the results of the transport approach and experimental data from the STAR experiment [12]. As the thermal production of charm quark pairs is negligible in the QGP at the temperatures reached in heavy-ion collisions at the RHIC and LHC colliding energies, we assume the charm quark number from the initial hard processes is conserved via an introduction of charm chemical potential at the freeze-out. The value of the chemical potential μc at the freeze-out temperature Tfrz is adjusted to fit the magnitude of the experimental data on the pT spectra at low pT in the most central collisions. For Tfrz=150 MeV, the effective charm chemical potential μc=700 MeV is found to fit the experimental data on D0 spectra (red) at low pT (pT<1 GeV) for 0-10% (solid), 10%-40% (dot-dashed), and 0-80% (dashed) centrality. With the same freeze-out temperature and the effective charm chemical potential in the most central collisions, the hydrodynamic model can describe the low pT spectra well at other centralities. However, the hydro calculations overpredict the pT spectra at pT>1 GeV/c.

    Figure 4

    Figure 4.  (color online) Left: hydrodynamic calculations (lines) of D0 pT spectra in Au+Au collisions at sNN=200 GeV for three different centrality bins, using two different sets of freeze-out temperatures and effective charm chemical potentials. They are compared to transport calculations in the same collision system (blue bands) and data from the STAR experiment [12]. We also compare the D-meson production in p+p collisions from FONLL calculation (green band) [42, 43] to the STAR measurement [52]. Right: hydrodynamic calculations of D0 meson v2 are compared to transport calculations and STAR data [13].

    Different hadrons can freeze-out at different temperatures in the hydrodynamic model of hadron production. For example, it was shown in the UrQMD model studies that protons freeze-out earlier than pions and kaons [53]. In the hydro calculation of D0 spectra (red in the left panel), we have assumed that D0 mesons freeze out (Tfrz=150 MeV) earlier than the light quark hadrons (Tfrz=137 MeV). If we set the freeze-out temperature of D0 the same as that of the light quark hadrons at 137 MeV with a re-adjusted charm chemical potential μc=800 MeV, the slope of the pT spectra becomes slightly larger (blue) due to increased radial expansion, further away from experimental data. Increasing the freeze-out temperature slightly decreases the slope of the spectra at high pT. However, the change is too small to describe the data at high pT.

    In the left panel of Fig. 4, we further compare the transport calculations (blue bands) of the D0 spectra in Au+Au collisions. Before discussing the results in Au+Au collisions, it is necessary to present the FONLL calculation of the baseline D0 spectra (green band) in p+p collisions. The D0 meson spectra in p+p collisions are computed using the FONLL program [42, 43], for which major uncertainty stems from the variation of the renormalization scale mT/2<μR<2mT, where mT=M2+p2T is the transverse mass of the heavy quark. From Fig. 4, the central prediction from FONLL (green line) systemically underestimates the STAR data on D0 production in p+p collisions at low pT, while the upper band of the prediction with μR=mT/2 is found to be closer to the measurements. Moving to results for nuclear collisions, we can see that the transport approach, which includes effects of non-equilibrium evolution, describes the experimental data on pT spectra at large pT>2 GeV/c. At low pT, the calculations are still consistent with experiments within error bands; however, the central points of the band under-predict the experimental data. We comment that the production of low-pT D mesons in nuclear collisions is sensitive to the small-x region of the nuclear parton distribution function (nPDF), which is less well-constrained than the parton distribution function of the proton. Despite the EPPS16 nPDF publishing its uncertainties, we have only shown calculations using the central fit of the EPPS16 nPDF.

    Comparison of the high-pT part of the spectra between the hydrodynamic results and experimental data shows that charm quarks have not reached complete equilibrium for pT>1 GeV/c, as the spectra from experiments are still below the hydrodynamic results, assuming a fully equilibrated system of charm quarks that flow with the QGP fluid. The enhanced ratios of Λc/D0 and Ds/D0 observed at the RHIC [54, 55] for Au+Au collisions and LHC [56, 57] for Pb+Pb collisions as compared with p+p collisions may partially contribute to this D0 suppression [58-61].

    Both the relativistic hydro and transport models can describe the experimental data on the elliptic flow of D0 mesons within the error bars as shown in the right panel of Fig. 4. One can see the mass ordering of the elliptic flow by comparing it to the elliptic flow of charged pions (blue) from the same hydrodynamic calculations. The hydrodynamic results on v2(pT) at low pT are insensitive to the value of the freeze-out temperature.

    We stress that even though the hydrodynamic model fails to describe the pT spectra of D0 beyond pT>1 GeV/c, its results on the elliptic flow agree with the experimental data and transport model calculation (solid blue) for pT up to 4 GeV/c. This implies that the experimental data on the elliptic flow cannot provide a stringent constraint on the kinetic equilibration of heavy quarks. The strong interaction between heavy quarks and the medium, however, can "drag" the heavy quarks along collective flow developed for the bulk medium even though the interaction might not drive the heavy quarks to full kinetic equilibration.

    We have calculated the D0 meson spectra and elliptic flow in Au+Au collisions at the RHIC energy within a relativistic viscous hydrodynamic model assuming that charm quarks initially produced through hard processes become fully kinetically equilibrated in the QGP. We neglect the thermal charm quark pair production in the QGP and use an effective charm chemical potential at the freeze-out, which is tuned to describe the number of charm quarks produced in the initial hard processes. With the charm freeze-out temperature Tfrz=137150 MeV, the hydrodynamic model can describe the experimental data on D0 elliptic flow for pT<4 GeV/c as well as the elliptic flow of light quark hadrons. The mass ordering of the elliptic flow for pions and D mesons in the experimental data is also observed for a freeze-out temperature between 137 and 150 MeV. We have also compared the hydrodynamic results with those from a transport model with both energy loss and momentum diffusion, which can describe the observed elliptic flow of D0 mesons at both low and high pT.

    The hydrodynamic model, however, fails to describe the pT spectra of D0 for pT>1 GeV/c, significantly over-predicting the D0 pT spectra at large pT. The transport model, in contrast, can describe the pT spectra well at large pT> 2 GeV/c due to parton energy loss. However, its central values under-predict the spectra at low pT due to the baseline spectra of the initial charm production through hard processes in p+p collisions.

    Our comparative study indicates that only the low-pT D mesons in the experiments might have reached kinetic equilibrium while charm quarks in the intermediate region of pT are partially thermalized. Parton energy loss and momentum diffusion in the transport model can describe well the non-equilibrium behavior of D mesons spectra and elliptic flow at large pT. The interaction of charm quarks in partial equilibrium with the medium in the intermediate pT, however, can develop an elliptic flow as large as that of fully equilibrated charm quarks.

    [1] B. A. Gelman, E. V. Shuryak, and I. Zahed, Phys. Rev. C 74, 044908 (2006), arXiv:nuclth/0601029 doi: 10.1103/PhysRevC.74.044908
    [2] U. W. Heinz, J. Phys. A 42, 214003 (2009), arXiv:0810.5529 [nucl-th doi: 10.1088/1751-8113/42/21/214003
    [3] S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 182301 (2003), arXiv:nucl-ex/0305013 doi: 10.1103/PhysRevLett.91.182301
    [4] J. Adams et al. (STAR), Phys. Rev. Lett. 92, 052302 (2004), arXiv:nucl-ex/0306007 doi: 10.1103/PhysRevLett.92.052302
    [5] B. B. Abelev et al. (ALICE), JHEP 06, 190 (2015), arXiv:1405.4632 [nucl-ex
    [6] C. Adler et al. (STAR), Phys. Rev. Lett. 90, 082302 (2003), arXiv:nucl-ex/0210033 doi: 10.1103/PhysRevLett.90.082302
    [7] K. Adcox et al. (PHENIX), Phys. Rev. Lett. 88, 022301 (2002), arXiv:nucl-ex/0109003
    [8] S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 072301 (2003), arXiv:nucl-ex/0304022 doi: 10.1103/PhysRevLett.91.072301
    [9] C. Adler et al. (STAR), Phys. Rev. Lett. 89, 202301 (2002), arXiv:nucl-ex/0206011 doi: 10.1103/PhysRevLett.89.202301
    [10] S. Chatrchyan et al. (CMS), Eur. Phys. J. C 72, 1945 (2012), arXiv:1202.2554 [nucl-ex doi: 10.1140/epjc/s10052-012-1945-x
    [11] L. Adamczyk et al. (STAR), Phys. Rev. Lett. 113, 142301 (2014), [Erratum: Phys.Rev.Lett. 121, 229901 (2018)], arXiv: 1404.6185 [nucl-ex]
    [12] J. Adam et al. (STAR), Phys. Rev. C 99, 034908 (2019), arXiv:1812.10224 [nucl-ex doi: 10.1103/PhysRevC.99.034908
    [13] L. Adamczyk et al. (STAR), Phys. Rev. Lett. 118, 212301 (2017), arXiv:1701.06060 [nucl-ex doi: 10.1103/PhysRevLett.118.212301
    [14] B. Abelev et al. (ALICE), JHEP 09, 112 (2012), arXiv:1203.2160 [nucl-ex
    [15] B. Abelev et al. (ALICE), Phys. Rev. Lett. 111, 102301 (2013), arXiv:1305.2707 [nucl-ex doi: 10.1103/PhysRevLett.111.102301
    [16] A. M. Sirunyan et al. (CMS), Phys. Lett. B 782, 474 (2018), arXiv:1708.04962 [nucl-ex doi: 10.1016/j.physletb.2018.05.074
    [17] A. M. Sirunyan et al. (CMS), Phys. Rev. Lett. 120, 202301 (2018), arXiv:1708.03497 [nucl-ex doi: 10.1103/PhysRevLett.120.202301
    [18] A. M. Sirunyan et al. (CMS), Phys. Rev. Lett. 121, 082301 (2018), arXiv:1804.09767 [hep-ex doi: 10.1103/PhysRevLett.121.082301
    [19] M. C. Wang and G. E. Uhlenbeck, Rev. Mod. Phys. 17, 323 (1945) doi: 10.1103/RevModPhys.17.323
    [20] R. Rapp and H. van Hees (2010) pp. 111–206, arXiv: 0903.1096 [hep-ph]
    [21] G. D. Moore and D. Teaney, Phys. Rev. C 71, 064904 (2005), arXiv:hep-ph/0412346 doi: 10.1103/PhysRevC.71.064904
    [22] S. Cao, T. Luo, G.-Y. Qin et al., Phys. Rev. C 94, 014909 (2016), arXiv:1605.06447 [nucl-th doi: 10.1103/PhysRevC.94.014909
    [23] M. Nahrgang, J. Aichelin, P. B. Gossiaux et al., Phys. Rev. C 93, 044909 (2016), arXiv:1602.03544 [nucl-th doi: 10.1103/PhysRevC.93.044909
    [24] W. Ke, Y. Xu, and S. A. Bass, Phys. Rev. C 98, 064901 (2018) doi: 10.1103/PhysRevC.98.064901
    [25] S. Cao et al., Phys. Rev. C 99, 054907 (2019), arXiv:1809.07894 [nucl-th doi: 10.1103/PhysRevC.99.054907
    [26] A. Beraudo et al., Nucl. Phys. A 979, 21 (2018), arXiv:1803.03824 [nucl-th doi: 10.1016/j.nuclphysa.2018.09.002
    [27] D. Banerjee, S. Datta, R. Gavai et al., Phys. Rev. D 85, 014510 (2012), arXiv:1109.5738 [hep-lat
    [28] H. Ding, A. Francis, O. Kaczmarek et al., Phys. Rev. D 86, 014509 (2012), arXiv:1204.4945 [hep-lat
    [29] A. Francis, O. Kaczmarek, M. Laine et al., Phys. Rev. D 92, 116003 (2015), arXiv:1508.04543 [hep-lat doi: 10.1103/PhysRevD.92.116003
    [30] Y. Xu, J. E. Bernhard, S. A. Bass et al., Phys. Rev. C 97, 014907 (2018), arXiv:1710.00807 [nucl-th doi: 10.1103/PhysRevC.97.014907
    [31] F. Cooper and G. Frye, Phys. Rev. D 10, 186 (1974)
    [32] L.-G. Pang, H. Petersen, and X.-N. Wang, Phys. Rev. C 97, 064918 (2018), arXiv:1802.04449 [nucl-th doi: 10.1103/PhysRevC.97.064918
    [33] J. C. Collins, D. E. Soper, and G. F. Sterman, Nucl. Phys. B 263, 37 (1986) doi: 10.1016/0550-3213(86)90026-X
    [34] E. Braaten, K.-m. Cheung, S. Fleming et al., Phys. Rev. D 51, 4819 (1995), arXiv:hep-ph/9409316 doi: 10.1103/PhysRevD.51.4819
    [35] M. Cacciari and P. Nason, JHEP 09, 006 (2003), arXiv:hep-ph/0306212
    [36] M. Cacciari, S. Frixione, N. Houdeau et al., JHEP 10, 137 (2012), arXiv:1205.6344 [hep-ph
    [37] R. Rapp and E. Shuryak, Phys. Rev. D 67, 074036 (2003), arXiv:hep-ph/0301245 doi: 10.1103/PhysRevD.67.074036
    [38] C. Patrignani et al. (Particle Data Group), Chin. Phys. C 40, 100001 (2016) doi: 10.1088/1674-1137/40/10/100001
    [39] J. S. Moreland, J. E. Bernhard, and S. A. Bass, Phys. Rev. C 92, 011901 (2015), arXiv:1412.4708 [nucl-th doi: 10.1103/PhysRevC.92.011901
    [40] J. E. Bernhard, J. S. Moreland, S. A. Bass et al., Phys. Rev. C 94, 024907 (2016), arXiv:1605.03954 [nucl-th doi: 10.1103/PhysRevC.94.024907
    [41] P. Huovinen and P. Petreczky, Nucl. Phys. A 837, 26 (2010), arXiv:0912.2541 [hep-ph doi: 10.1016/j.nuclphysa.2010.02.015
    [42] M. Cacciari, M. Greco, and P. Nason, JHEP 05, 007 (1998), arXiv:hep-ph/9803400
    [43] M. Cacciari, S. Frixione, and P. Nason, JHEP 03, 006 (2001), arXiv:hep-ph/0102134
    [44] K. J. Eskola, P. Paakkinen, H. Paukkunen et al., Eur. Phys. J. C 77, 163 (2017), arXiv:1612.05741 [hep-ph doi: 10.1140/epjc/s10052-017-4725-9
    [45] S. Cao, G.-Y. Qin, and S. A. Bass, Phys. Rev. C 88, 044907 (2013), arXiv:1308.0617 [nucl-th doi: 10.1103/PhysRevC.88.044907
    [46] Z.-w. Lin, T. G. Di, and C. M. Ko, Nucl. Phys. A 689, 965 (2001), arXiv:nucl-th/0006086 [nucl-th doi: 10.1016/S0375-9474(00)00611-4
    [47] S. A. Bass et al., Prog. Part. Nucl. Phys. 41, 255 (1998), arXiv:nucl-th/9803035 [nucl-th doi: 10.1016/S0146-6410(98)00058-1
    [48] M. Bleicher et al., J. Phys. G 25, 1859 (1999), arXiv:hep-ph/9909407 [hep-ph doi: 10.1088/0954-3899/25/9/308
    [49] B. Back et al., Phys. Rev. Lett. 91, 052303 (2003), arXiv:nucl-ex/0210015 doi: 10.1103/PhysRevLett.91.052303
    [50] S. Adler et al. (PHENIX), Phys. Rev. C 69, 034909 (2004), arXiv:nucl-ex/0307022 doi: 10.1103/PhysRevC.69.034909
    [51] B. I. Abelev et al. (STAR), Phys. Rev. C 77, 054901 (2008), arXiv:0801.3466 [nucl-ex doi: 10.1103/PhysRevC.77.054901
    [52] L. Adamczyk et al. (STAR), Phys. Rev. D 86, 072013 (2012), arXiv:1204.4244 [nucl-ex doi: 10.1103/PhysRevD.86.072013
    [53] H. van Hecke, H. Sorge, and N. Xu, Phys. Rev. Lett. 81, 5764 (1998), arXiv:nucl-th/9804035 doi: 10.1103/PhysRevLett.81.5764
    [54] L. Zhou (STAR), Nucl. Phys. A 967, 620 (2017), arXiv:1704.04364 [nucl-ex doi: 10.1016/j.nuclphysa.2017.05.114
    [55] J. Adam et al. (STAR), Phys. Rev. Lett. 124, 172301 (2020), arXiv:1910.14628 [nucl-ex doi: 10.1103/PhysRevLett.124.172301
    [56] S. Acharya et al. (ALICE), JHEP 10, 174 (2018), arXiv:1804.09083 [nucl-ex
    [57] S. Acharya et al. (ALICE), Phys. Lett. B 793, 212 (2019), arXiv:1809.10922 [nucl-ex doi: 10.1016/j.physletb.2019.04.046
    [58] V. Greco, C. M. Ko, and P. Levai, Phys. Rev. C 68, 034904 (2003), arXiv:nucl-th/0305024 doi: 10.1103/PhysRevC.68.034904
    [59] S. Cao, K.-J. Sun, S.-Q. Li et al., Phys. Lett. B 807, 135561 (2020), arXiv:1911.00456 [nucl-th doi: 10.1016/j.physletb.2020.135561
    [60] M. He and R. Rapp, Phys. Rev. Lett. 124, 042301 (2020), arXiv:1905.09216 [nucl-th doi: 10.1103/PhysRevLett.124.042301
    [61] R.-Q. Wang, J. Song, F.-L. Shao et al., Phys. Rev. C 101, 054903 (2020), arXiv:1911.00823 [hep-ph doi: 10.1103/PhysRevC.101.054903
  • [1] B. A. Gelman, E. V. Shuryak, and I. Zahed, Phys. Rev. C 74, 044908 (2006), arXiv:nuclth/0601029 doi: 10.1103/PhysRevC.74.044908
    [2] U. W. Heinz, J. Phys. A 42, 214003 (2009), arXiv:0810.5529 [nucl-th doi: 10.1088/1751-8113/42/21/214003
    [3] S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 182301 (2003), arXiv:nucl-ex/0305013 doi: 10.1103/PhysRevLett.91.182301
    [4] J. Adams et al. (STAR), Phys. Rev. Lett. 92, 052302 (2004), arXiv:nucl-ex/0306007 doi: 10.1103/PhysRevLett.92.052302
    [5] B. B. Abelev et al. (ALICE), JHEP 06, 190 (2015), arXiv:1405.4632 [nucl-ex
    [6] C. Adler et al. (STAR), Phys. Rev. Lett. 90, 082302 (2003), arXiv:nucl-ex/0210033 doi: 10.1103/PhysRevLett.90.082302
    [7] K. Adcox et al. (PHENIX), Phys. Rev. Lett. 88, 022301 (2002), arXiv:nucl-ex/0109003
    [8] S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 072301 (2003), arXiv:nucl-ex/0304022 doi: 10.1103/PhysRevLett.91.072301
    [9] C. Adler et al. (STAR), Phys. Rev. Lett. 89, 202301 (2002), arXiv:nucl-ex/0206011 doi: 10.1103/PhysRevLett.89.202301
    [10] S. Chatrchyan et al. (CMS), Eur. Phys. J. C 72, 1945 (2012), arXiv:1202.2554 [nucl-ex doi: 10.1140/epjc/s10052-012-1945-x
    [11] L. Adamczyk et al. (STAR), Phys. Rev. Lett. 113, 142301 (2014), [Erratum: Phys.Rev.Lett. 121, 229901 (2018)], arXiv: 1404.6185 [nucl-ex]
    [12] J. Adam et al. (STAR), Phys. Rev. C 99, 034908 (2019), arXiv:1812.10224 [nucl-ex doi: 10.1103/PhysRevC.99.034908
    [13] L. Adamczyk et al. (STAR), Phys. Rev. Lett. 118, 212301 (2017), arXiv:1701.06060 [nucl-ex doi: 10.1103/PhysRevLett.118.212301
    [14] B. Abelev et al. (ALICE), JHEP 09, 112 (2012), arXiv:1203.2160 [nucl-ex
    [15] B. Abelev et al. (ALICE), Phys. Rev. Lett. 111, 102301 (2013), arXiv:1305.2707 [nucl-ex doi: 10.1103/PhysRevLett.111.102301
    [16] A. M. Sirunyan et al. (CMS), Phys. Lett. B 782, 474 (2018), arXiv:1708.04962 [nucl-ex doi: 10.1016/j.physletb.2018.05.074
    [17] A. M. Sirunyan et al. (CMS), Phys. Rev. Lett. 120, 202301 (2018), arXiv:1708.03497 [nucl-ex doi: 10.1103/PhysRevLett.120.202301
    [18] A. M. Sirunyan et al. (CMS), Phys. Rev. Lett. 121, 082301 (2018), arXiv:1804.09767 [hep-ex doi: 10.1103/PhysRevLett.121.082301
    [19] M. C. Wang and G. E. Uhlenbeck, Rev. Mod. Phys. 17, 323 (1945) doi: 10.1103/RevModPhys.17.323
    [20] R. Rapp and H. van Hees (2010) pp. 111–206, arXiv: 0903.1096 [hep-ph]
    [21] G. D. Moore and D. Teaney, Phys. Rev. C 71, 064904 (2005), arXiv:hep-ph/0412346 doi: 10.1103/PhysRevC.71.064904
    [22] S. Cao, T. Luo, G.-Y. Qin et al., Phys. Rev. C 94, 014909 (2016), arXiv:1605.06447 [nucl-th doi: 10.1103/PhysRevC.94.014909
    [23] M. Nahrgang, J. Aichelin, P. B. Gossiaux et al., Phys. Rev. C 93, 044909 (2016), arXiv:1602.03544 [nucl-th doi: 10.1103/PhysRevC.93.044909
    [24] W. Ke, Y. Xu, and S. A. Bass, Phys. Rev. C 98, 064901 (2018) doi: 10.1103/PhysRevC.98.064901
    [25] S. Cao et al., Phys. Rev. C 99, 054907 (2019), arXiv:1809.07894 [nucl-th doi: 10.1103/PhysRevC.99.054907
    [26] A. Beraudo et al., Nucl. Phys. A 979, 21 (2018), arXiv:1803.03824 [nucl-th doi: 10.1016/j.nuclphysa.2018.09.002
    [27] D. Banerjee, S. Datta, R. Gavai et al., Phys. Rev. D 85, 014510 (2012), arXiv:1109.5738 [hep-lat
    [28] H. Ding, A. Francis, O. Kaczmarek et al., Phys. Rev. D 86, 014509 (2012), arXiv:1204.4945 [hep-lat
    [29] A. Francis, O. Kaczmarek, M. Laine et al., Phys. Rev. D 92, 116003 (2015), arXiv:1508.04543 [hep-lat doi: 10.1103/PhysRevD.92.116003
    [30] Y. Xu, J. E. Bernhard, S. A. Bass et al., Phys. Rev. C 97, 014907 (2018), arXiv:1710.00807 [nucl-th doi: 10.1103/PhysRevC.97.014907
    [31] F. Cooper and G. Frye, Phys. Rev. D 10, 186 (1974)
    [32] L.-G. Pang, H. Petersen, and X.-N. Wang, Phys. Rev. C 97, 064918 (2018), arXiv:1802.04449 [nucl-th doi: 10.1103/PhysRevC.97.064918
    [33] J. C. Collins, D. E. Soper, and G. F. Sterman, Nucl. Phys. B 263, 37 (1986) doi: 10.1016/0550-3213(86)90026-X
    [34] E. Braaten, K.-m. Cheung, S. Fleming et al., Phys. Rev. D 51, 4819 (1995), arXiv:hep-ph/9409316 doi: 10.1103/PhysRevD.51.4819
    [35] M. Cacciari and P. Nason, JHEP 09, 006 (2003), arXiv:hep-ph/0306212
    [36] M. Cacciari, S. Frixione, N. Houdeau et al., JHEP 10, 137 (2012), arXiv:1205.6344 [hep-ph
    [37] R. Rapp and E. Shuryak, Phys. Rev. D 67, 074036 (2003), arXiv:hep-ph/0301245 doi: 10.1103/PhysRevD.67.074036
    [38] C. Patrignani et al. (Particle Data Group), Chin. Phys. C 40, 100001 (2016) doi: 10.1088/1674-1137/40/10/100001
    [39] J. S. Moreland, J. E. Bernhard, and S. A. Bass, Phys. Rev. C 92, 011901 (2015), arXiv:1412.4708 [nucl-th doi: 10.1103/PhysRevC.92.011901
    [40] J. E. Bernhard, J. S. Moreland, S. A. Bass et al., Phys. Rev. C 94, 024907 (2016), arXiv:1605.03954 [nucl-th doi: 10.1103/PhysRevC.94.024907
    [41] P. Huovinen and P. Petreczky, Nucl. Phys. A 837, 26 (2010), arXiv:0912.2541 [hep-ph doi: 10.1016/j.nuclphysa.2010.02.015
    [42] M. Cacciari, M. Greco, and P. Nason, JHEP 05, 007 (1998), arXiv:hep-ph/9803400
    [43] M. Cacciari, S. Frixione, and P. Nason, JHEP 03, 006 (2001), arXiv:hep-ph/0102134
    [44] K. J. Eskola, P. Paakkinen, H. Paukkunen et al., Eur. Phys. J. C 77, 163 (2017), arXiv:1612.05741 [hep-ph doi: 10.1140/epjc/s10052-017-4725-9
    [45] S. Cao, G.-Y. Qin, and S. A. Bass, Phys. Rev. C 88, 044907 (2013), arXiv:1308.0617 [nucl-th doi: 10.1103/PhysRevC.88.044907
    [46] Z.-w. Lin, T. G. Di, and C. M. Ko, Nucl. Phys. A 689, 965 (2001), arXiv:nucl-th/0006086 [nucl-th doi: 10.1016/S0375-9474(00)00611-4
    [47] S. A. Bass et al., Prog. Part. Nucl. Phys. 41, 255 (1998), arXiv:nucl-th/9803035 [nucl-th doi: 10.1016/S0146-6410(98)00058-1
    [48] M. Bleicher et al., J. Phys. G 25, 1859 (1999), arXiv:hep-ph/9909407 [hep-ph doi: 10.1088/0954-3899/25/9/308
    [49] B. Back et al., Phys. Rev. Lett. 91, 052303 (2003), arXiv:nucl-ex/0210015 doi: 10.1103/PhysRevLett.91.052303
    [50] S. Adler et al. (PHENIX), Phys. Rev. C 69, 034909 (2004), arXiv:nucl-ex/0307022 doi: 10.1103/PhysRevC.69.034909
    [51] B. I. Abelev et al. (STAR), Phys. Rev. C 77, 054901 (2008), arXiv:0801.3466 [nucl-ex doi: 10.1103/PhysRevC.77.054901
    [52] L. Adamczyk et al. (STAR), Phys. Rev. D 86, 072013 (2012), arXiv:1204.4244 [nucl-ex doi: 10.1103/PhysRevD.86.072013
    [53] H. van Hecke, H. Sorge, and N. Xu, Phys. Rev. Lett. 81, 5764 (1998), arXiv:nucl-th/9804035 doi: 10.1103/PhysRevLett.81.5764
    [54] L. Zhou (STAR), Nucl. Phys. A 967, 620 (2017), arXiv:1704.04364 [nucl-ex doi: 10.1016/j.nuclphysa.2017.05.114
    [55] J. Adam et al. (STAR), Phys. Rev. Lett. 124, 172301 (2020), arXiv:1910.14628 [nucl-ex doi: 10.1103/PhysRevLett.124.172301
    [56] S. Acharya et al. (ALICE), JHEP 10, 174 (2018), arXiv:1804.09083 [nucl-ex
    [57] S. Acharya et al. (ALICE), Phys. Lett. B 793, 212 (2019), arXiv:1809.10922 [nucl-ex doi: 10.1016/j.physletb.2019.04.046
    [58] V. Greco, C. M. Ko, and P. Levai, Phys. Rev. C 68, 034904 (2003), arXiv:nucl-th/0305024 doi: 10.1103/PhysRevC.68.034904
    [59] S. Cao, K.-J. Sun, S.-Q. Li et al., Phys. Lett. B 807, 135561 (2020), arXiv:1911.00456 [nucl-th doi: 10.1016/j.physletb.2020.135561
    [60] M. He and R. Rapp, Phys. Rev. Lett. 124, 042301 (2020), arXiv:1905.09216 [nucl-th doi: 10.1103/PhysRevLett.124.042301
    [61] R.-Q. Wang, J. Song, F.-L. Shao et al., Phys. Rev. C 101, 054903 (2020), arXiv:1911.00823 [hep-ph doi: 10.1103/PhysRevC.101.054903
  • 加载中

Cited by

1. Chen, X., Wu, X.-Y., Cao, S. et al. System-size and shape dependencies of collective-flow fluctuations in relativistic nuclear collisions[J]. Physical Review C, 2024, 109(6): 064915. doi: 10.1103/PhysRevC.109.064915
2. Wang, S., Dai, W., Wang, E. et al. Heavy-Flavour Jets in High-Energy Nuclear Collisions[J]. Symmetry, 2023, 15(3): 727. doi: 10.3390/sym15030727
3. Wu, X.-Y., Qin, G.-Y., Pang, L.-G. et al. (3+1)-D viscous hydrodynamics at finite net baryon density: Identified particle spectra, anisotropic flows, and flow fluctuations across energies relevant to the beam-energy scan at RHIC[J]. Physical Review C, 2022, 105(3): 034909. doi: 10.1103/PhysRevC.105.034909
4. Habashy, D.M., El-Bakry, M.Y., Tawfik, A.N. et al. Particles multiplicity based on rapidity in Landau and artificial neural network (ANN) models[J]. International Journal of Modern Physics A, 2022, 37(2): 2250002. doi: 10.1142/S0217751X22500026

Figures(4)

Get Citation
Chi Ding, Wei-Yao Ke, Long-Gang Pang and Xin-Nian Wang. Hydrodynamic description of D meson production in high-energy heavy-ion collisions[J]. Chinese Physics C. doi: 10.1088/1674-1137/abf645
Chi Ding, Wei-Yao Ke, Long-Gang Pang and Xin-Nian Wang. Hydrodynamic description of D meson production in high-energy heavy-ion collisions[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abf645 shu
Milestone
Received: 2021-01-05
Article Metric

Article Views(1395)
PDF Downloads(40)
Cited by(4)
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:

Hydrodynamic description of D meson production in high-energy heavy-ion collisions

  • 1. Key Laboratory of Quark & Lepton Physics (MOE) and Institute of Particle Physics, Central China Normal University, Wuhan 430079, China
  • 2. Physics Department, University of California, Berkeley, CA 94720, USA
  • 3. Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA

Abstract: The large values and constituent-quark-number scaling of the elliptic flow of low-pT D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low-pT charm quarks, we compare the D0 meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. D0 mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D0 mesons at pT<3 GeV/c as measured in Au+Au collisions at sNN=200 GeV. Though the experimental data on D0 spectra are consistent with the hydrodynamic result at small pT1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum,pT>2 GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Our comparative study indicates that charm quarks only approach local thermal equilibrium at small pT, even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate pT.

    HTML

    I.   INTRODUCTION
    • The large collective flow and momentum anisotropy as manifested in the final hadron spectra in high-energy heavy-ion collisions indicate a strong collectivity of the produced dense and hot nuclear matter during its dynamic evolution [1, 2]. The observed approximate constituent quark number (NCQ) scaling of the elliptic flow for light quark hadrons [3-6] and strong jet quenching [7-10] in high-energy heavy-ion collisions at both the Relativistic Heavy-ion Collider (RHIC) and the Large Hadron Collider (LHC) suggest the formation of a hot, deconfined, and opaque quark-gluon plasma (QGP). Similar phenomena are also observed for heavy quark mesons [11-17]. Experimental data from the RHIC and LHC show a large elliptic flow of D0 mesons at low-pT, which also approximately obeys the NCQ scaling of light quark hadrons [13, 18]. Heavy quarks are predominantly produced through initial hard processes and have mass scales much larger than the typical temperature of the QGP medium. Therefore, it still remains an interesting question as to whether and to what degree charm quarks become thermalized [2] and flow with the QGP due to their strong interaction with the hot medium.

      Past studies usually describe the thermalization processes of heavy quarks with transport equations. Their low-momentum (pM) dynamics in the QGP medium are usually treated as Brownian motion [19-21], considering the masses of heavy quarks to be much larger than the typical temperatures in the medium MT. For the transport of heavy quarks with intermediate momentum pM, both the collisional and radiative energy loss have to be considered. For high-momentum heavy quarks, most transport models [22-24] consider the gluon radiations and interactions between heavy quarks and the medium as perturbative processes or use transport coefficients with weakly-coupled assumptions as inputs.

      It was found that transport models need a large heavy-quark momentum-diffusion parameter (ˆq=0.7...2.7 GeV2/fm at T=300 MeV and p=30 GeV/c [25]) to describe the observed medium modification of the heavy-quark meson spectra and elliptic flow coefficients [25, 26]. Several studies have attempted to evaluate the heavy quark diffusion coefficient at zero momentum non-perturbatively from lattice QCD [27-29]. The estimated interaction strength is comparable to the values from phenomenological determination [24, 30] and much larger than the natural expectation of the weakly-coupled theory at the leading order. This poses questions regarding the weak-coupling assumption on the nature of the interaction between heavy quarks and the QGP medium. It is possible that, with a strong coupling, the dynamics of the low-pT charm quarks in the QGP are better described by a hydrodynamic approach.

      In this study, we compute the production of low-pT D mesons in an extreme limit of hydrodynamic evolution with the following assumptions. First, the number of (anti)charm quarks is conserved during the lifetime of the fireball, i.e., thermal production of heavy quark pairs is negligible. Second, the coupling is so strong that the initially produced charm quarks quickly diffuse into the QGP and reach local kinetic equilibrium. Finally, the differential yield of D mesons is computed using the Cooper-Frye formula [31] on the hydrodynamic freeze-out hypersurface, with an effective charm chemical potential to guarantee the charm yield is the same as measured in the experiment. We compare the resultant pT-dependent spectra and elliptic flow of low-pT D mesons, which are considered as the manifestation of the extreme limit of complete thermalization of heavy quarks, to the experimental data as well as transport calculations to quantify the degree of heavy-quark thermalization in heavy-ion collisions. We use the (3+1)D hydrodynamic model (CLVisc) [32] and a linearized Boltzmann-Langevin model [24] for heavy quark transport.

      The paper is organized as follows. In Section II, we introduce the relativistic hydrodynamic model and the Boltzmann-Langevin transport model for the calculations of pT spectra and elliptic flow of charmed mesons. In Section III, we first discuss the calibration of the hydro model to the experimental data on light quark hadron spectra and elliptic flow in Au+Au collisions at the RHIC energy. We then compare the hydro results on the pT spectra and elliptic flow of charmed mesons at low pT to the experimental data and transport calculations and discuss the implications on the degree of thermalization for heavy quarks in the QGP formed in high-energy heavy-ion collisions.

    II.   MODEL DESCRIPTIONS

      A.   Charm meson production in the limit of local kinetic equilibrium

    • In the usual paradigm for charmed meson production in high-energy heavy-ion collisions, charm quarks are produced through the initial hard scatterings [33]. These heavy quarks experience energy loss and momentum diffusion in the QGP through both elastic and inelastic collisions with the medium, which can be modeled as the drag and diffusion coefficients in the Boltzmann-Langevin equations. The final D mesons are formed through charm quark fragmentation at high pT [34-36] or charm-light quark recombination at low and intermediate pT [37]. Such transport models for initial heavy quark and final meson production can describe the charm meson spectra and elliptic flow well in high-energy heavy-ion collisions [25, 26]. However, it is still interesting to investigate whether and to what degree the heavy quarks achieve kinematic equilibrium. If these heavy quarks indeed reach local thermal equilibrium, one should expect that the hydrodynamic model can also describe heavy meson spectra and elliptic flow. Furthermore, a hydrodynamic picture overcomes the difficulty of the transport equation when the coupling becomes large between charm quarks or mesons and the medium.

      In this study, we consider the extreme scenario and investigate the charmed meson spectra and elliptic flow in the limit of the fully thermalized low-pT charm quarks. These low-pT charm quarks are still produced initially through hard processes; however, the interactions are assumed to be strong enough that they quickly diffuse into the medium and lose the memory of their initial distribution in phase space, both in terms of the spatial location and momentum of the initial production through hard processes. In the momentum space, they are assumed to reach full kinetic equilibrium and comove with the medium, flowing with the strongly coupled QGP as light quarks and gluons. Finally, the charm quarks transit to charm hadrons with a phase transition of the bulk medium, where we assume the interaction in the hadronic phase is still strong enough to maintain the kinetic equilibrium and equilibrium ratios of different species of charm mesons. In this situation, the low-pT D mesons are produced in the same way as other light hadrons on the freeze-out hypersurface in a relativistic hydrodynamic model. The temperature and fluid velocity profiles on the freeze-out hypersurface are crucial for a reasonable estimate of the D-meson spectra and elliptic flow in the limit of a complete heavy-quark thermalization. We will use the rapidity distribution, pT spectra, and elliptic flow of charged pions to calibrate the hydrodynamic model.

      For completeness, we have to consider D0 mesons from the feed-down of D. Using the feed-down tables from [37, 38],

      D(2007)064.7%D0(1865)+π0,

      (1)

      D(2007)035.3%D0(1865)+γ,

      (2)

      D(2010)+68%D0(1865)+π+

      (3)

      Note that the D's are vector mesons with spin 1. Therefore, the spin degeneracy ratio for D/D is 3:1.

    • B.   Relativistic hydrodynamics

    • The hydrodynamic model we use, CLVisc, is a (3+1)D viscous hydrodynamic model parallelized on GPU using OpenCL [32]. The program is well tested against several analytical solutions and can describe the bulk hadron spectra and anisotropic flow in high-energy heavy-ion collisions at top RHIC and LHC energy. It simulates the fluid dynamic evolution of the strongly interacting QCD matter created in high-energy heavy-ion collisions by solving the fluid dynamic equations

      μTμν=0,  with  Tμν=εuμuνPΔμν+πμν,

      (4)

      where ε is the energy density, P is the pressure as a function of energy density given by the equation of state (EoS), Δμν=gμνuμuν is a projection operator, uμ is the fluid four-velocity obeying uμuμ=1 and πμν is the shear stress tensor.

      The initial condition for entropy density distribution in the transverse plane is provided by the Trento Monte Carlo model [39, 40]. An envelope function is used to approximate the longitudinal distribution along the space-time rapidity,

      H(ηs)=exp[(|ηs|ηw)22σ2ηθ(|ηs|ηw)],

      (5)

      where ση=1.5 and σw=1.3 are used for Au+Au collisions at sNN=200 GeV.

      We have assumed an initial time for the hydrodynamics τ0=0.6 fm. In the present study, we use the partial chemical equilibrium EoS with chemical freeze-out temperature 165 MeV and a smooth crossover between a QGP at high temperature and hadron resonance gas (HRG) EoS at low temperature [41] as inspired by the lattice QCD study.

      Baryons and mesons passing through the freeze-out hyper-surface are assumed to obey Fermi-Dirac and Bose-Einstein distributions, respectively. Their momentum distributions are given by the Cooper-Frye formula [31],

      dNidypTdpTdϕ=gi(2π)3pμdΣμf(pu)(1+δf),

      (6)

      where gi=2spin+1 is the spin degeneracy, pμ is the four-momenta of produced particles in the laboratory frame, Σμ is the freeze-out hyper-surface, f(pu) is the Fermi-Dirac/Bose-Einstein distribution function,

      f(pu)=1exp[(puμi)/Tfrz]±1,

      (7)

      and δf is the non-equilibrium correction,

      δf=(1feq)pμpνπμν2T2frz(ε+P).

      (8)

      We have chosen the freeze-out temperature Tfz=137 MeV for light flavor hadrons. The freeze-out temperatures for D mesons are different, and we consider several values to provide an estimate of the uncertainties.

      The elliptic flow of D mesons is defined as the second coefficient of the Fourier decomposition of their azimuthal angle distributions with respect to the event plane of light hadrons,

      d3NpTdpTdydϕ=d2N2πpTdpTdy[1+n=12vncos(n(ϕΨEP))].

      (9)
    • C.   Transport approach of heavy flavor evolution

    • We compare the charm meson spectra and flow calculated from the hydrodynamic freeze-out model to that from a transport model [24]. The transport approach assumes that heavy quarks, including those with low momentum in the comoving frame of the medium, remain good quasi-particles in the QGP. Therefore, the dynamics of low-pT heavy flavors can be described by a Boltzmann-type transport equation,

      (t+v)fQ(t,x,p)=[dRdq3(p+q,q)fQ(p+q)dRdq3(p,q)fQ(p)]dq3.

      (10)

      Here, fQ(t,x,p) is the phase-space density of heavy flavors: heavy quarks or heavy mesons (heavy baryons are omitted in this study). dR(p,q)/dq3 is the differential rate for a heavy flavor particle with momentum p to transfer three-momentum q to the local medium with flow velocity u and temperature T.

      The hydrodynamic and the transport model have both overlapping and distinct regimes of application. Both models contain the equilibrium situation. Compared to the hydrodynamic approach, a transport model also governs the far-from-equilibrium dynamics of heavy flavor particles. The traces of off-equilibrium effects can be important in a finite and expanding plasma with moderate coupling between heavy quarks and the medium. In the meantime, hydrodynamics can describe the evolution with large couplings and in the non-quasiparticle regime, which is beyond the applicability of the transport approach. Therefore, both models provide complementary pictures to understand the experimental data.

      The initial charm quark spectrum that initializes the transport equation is obtained from the perturbative QCD based FONLL [42, 43] calculation with EPPS16 nuclear parton distribution function [44]. In the medium above the pseudo-critical temperature (T>Tc), we assume the heavy flavors exist as deconfined heavy quarks. The interaction rate consists of two parts with a comparable contribution to heavy quark energy loss at intermediate pT region: collisional processes and medium-induced gluon radiations.

      Collisional process in the QGP Collisional processes are mainly modeled by two-body scatterings between the heavy quark and medium partons,

      dRi=q,gdq3=¯|M2QiQi|dif(pi)16π2EQEQd3pi2Eid3pi2Ei×δ(4)(pQ+pipQpi),

      (11)

      where Q (i) and Q (i) labels the initial and final state heavy quark (medium quark or gluon). pQ=pQ+q is the three-momentum transfer to the heavy quark. ¯|M2QiQi| is the squared amplitude of the two-body collision at leading order, averaging over the initial-stage quantum number and summed over the final-state quantum number. The ˆt-channel divergence is screened using the QCD non-perturbative scale Λ=0.2 GeV and the Debye screening mass mD=6παsT of a three-flavor plasma,

      ¯|M2QiQi|1ˆt21(ˆtΛ2)(ˆtm2D).

      (12)

      Finally, f(pi)=epiu/T is the classical thermal distribution function of the medium parton in plasma with local four-velocity u and temperature T, obtained from hydrodynamic simulations. The degeneracy factors for medium quarks and gluons are dg=2dA=16 and dq=4NfNc=36, respectively. Besides perturbative scatterings, an additional effective heavy-quark diffusion constant of the form

      ΔκDT3=A+BET

      (13)

      is introduced to mimic a possible non-perturbative contribution that peaks at low energy and low temperature. The effective heavy-quark-medium coupling parameter and the parameters in the non-perturbative diffusion constant are tuned in [24] to the suppression and momentum anisotropy of the production of open heavy flavor particles.

      Medium-induced radiative process in the QGP phase Heavy quarks can radiate additional gluons in the collision with medium partons. This is treated similarly to the rate in Eq. (11) for collisional processes, using two-to-three-body matrix-elements. One complication is that energetic gluon radiation in the rest frame of the medium is suppressed due to the QCD Landau-Pomeranchuk-Migdal (LPM) effect. In [24], this was included by restricting the phase-space of the radiated gluon using an interference factor motivated by the Higher-Twist approach 2[1cos(t/τf)], where t is the time since the last gluon radiation, while τf=2x(1x)EQ/(k2+x2M2) is the gluon formation time. Despite the LPM suppression, the radiative energy loss was found to be equally important in the intermediate pT region.

      Heavy quark hadronization and hadronic rescattering At the pseudo-critical temperature, charm quarks hadronize to D-mesons through both fragmentation and recombination mechanisms [45]. In the hadronic phase, D-mesons continued to interact with the light mesons via D-π and D-ρ scatterings [46] as implemented in Ultra-relativistic Quantum Molecular Dynamics (UrQMD) [47, 48].

    III.   RESULTS

      A.   Model calibration with experimental data

    • To calibrate the relativistic hydrodynamic model, we have computed the pseudorapidity density of charged particles dNch/dη as shown in Fig. 1, and the pT spectra and the elliptic flow v2 of π+ as shown in Fig. 2 for central Au+Au collisions with centrality range 0-5% at sNN=200 GeV. The calibrated parameters are the scale factor of 53 multiplied to the initial entropy density from the Trento model [39, 40], the starting time for hydro τ0=0.6 fm, the shear viscosity over entropy density ratio η/s=0.15, and the freeze-out temperature for light hadrons Tfz=137 MeV. These values have been tuned such that the predicted dNch/dη, pT spectra and v2(pT) of charged pions from relativistic hydrodynamics agree with the experimental data. Relativistic hydrodynamics with the same set of parameters also describes data well at other centralities and we refer readers to Ref. [32] for a more detailed description of CLVisc and the parameters. These parameters are fixed in the following calculations of the heavy quark meson spectra, except the freeze-out temperatures for D mesons. We will vary the value of the D meson freeze-out temperature and study its effect on the final D meson spectra and the elliptic flow.

      Figure 1.  (color online) Pseudorapidity density of charged hadrons from integrating the Cooper-Frye formula over the freeze-out hyperspace obtained in the CLVisc calculation. It agrees with the data of the PHOBOS experiment for Au+Au collisions at sNN=200GeV [49].

      Figure 2.  (color online) Left: transverse momentum spectra of π+ in Au+Au collisions. Right: anisotropic flow coefficient v2 in Au+Au collisions using the event-plane method. In both panels, CLVisc calculations (lines) are compared to data from the PHENIX experiment [50] and STAR experiment [51].

    • B.   Effects of the freeze-out temperature and resonance decay on D spectra

    • The effects of freeze-out temperature and resonance decays on D0 meson production are studied in this section. In the left panel of Fig. 3, we show the ratio between D0 pT spectra with and without contributions from the resonance decays. In the right panel of Fig. 3, we compare the D0 meson elliptic flow as a function of pT with (dot-dashed) and without (solid) resonance decays for three different values of the D meson freeze-out temperatures, Tfz=100, 137, and 155 MeV. We see that the resonance decay has a much larger effect on the pT spectra of D0 (left panel) than that on the elliptic flow v2 of D0 (right panel). For a given D meson freeze-out temperature, resonance decays from D contribute more D0 mesons at low pT than that at high pT. The ratio between D0 with and without decay decreases as pT increases for all three freeze-out temperatures considered here. As one increases the freeze-out temperature from 100 MeV to 137 MeV and 155 MeV, not only the D0 yield but also the ratio between D0 with and without decay increases, from approximately 2.35 to 3.0 and 3.3 at pT<1 GeV. The elliptic flow without decay (black solid), in contrast, almost overlaps that with decays (red dash-dotted) for all three different freeze-out temperatures; therefore, it is insensitive to the contributions from the resonance decays. The freeze-out temperature, however, has a large effect on the shape and magnitude of the elliptic flow as a function of transverse momentum. At small transverse momentum (pT < 2 GeV), the elliptic flow increases with the freeze-out temperature. At large transverse momentum (pT > 2 GeV), the trend goes the opposite way – the elliptic flow decreases as the freeze-out temperature increases.

      Figure 3.  (color online) Left: freeze-out temperature dependence of the effect of resonance decay on the D0 pT spectra. Right: freeze-out temperature dependence of the effect of resonance decay on the D0 flow coefficient v2. Calculations in both panels use one-shot hydrodynamic simulation of 20%-30% central Au+Au collisions at sNN=200 GeV.

    • C.   Model and data comparisons of the pT spectra and v2 of D0 mesons

    • In Fig. 4, we show the hydrodynamic calculations of the transverse momentum (pT) spectra (left) and the pT-dependent elliptic flow (right) of D0 mesons in Au+Au collisions at sNN=200 GeV for different centralities as compared with the results of the transport approach and experimental data from the STAR experiment [12]. As the thermal production of charm quark pairs is negligible in the QGP at the temperatures reached in heavy-ion collisions at the RHIC and LHC colliding energies, we assume the charm quark number from the initial hard processes is conserved via an introduction of charm chemical potential at the freeze-out. The value of the chemical potential μc at the freeze-out temperature Tfrz is adjusted to fit the magnitude of the experimental data on the pT spectra at low pT in the most central collisions. For Tfrz=150 MeV, the effective charm chemical potential μc=700 MeV is found to fit the experimental data on D0 spectra (red) at low pT (pT<1 GeV) for 0-10% (solid), 10%-40% (dot-dashed), and 0-80% (dashed) centrality. With the same freeze-out temperature and the effective charm chemical potential in the most central collisions, the hydrodynamic model can describe the low pT spectra well at other centralities. However, the hydro calculations overpredict the pT spectra at pT>1 GeV/c.

      Figure 4.  (color online) Left: hydrodynamic calculations (lines) of D0 pT spectra in Au+Au collisions at sNN=200 GeV for three different centrality bins, using two different sets of freeze-out temperatures and effective charm chemical potentials. They are compared to transport calculations in the same collision system (blue bands) and data from the STAR experiment [12]. We also compare the D-meson production in p+p collisions from FONLL calculation (green band) [42, 43] to the STAR measurement [52]. Right: hydrodynamic calculations of D0 meson v2 are compared to transport calculations and STAR data [13].

      Different hadrons can freeze-out at different temperatures in the hydrodynamic model of hadron production. For example, it was shown in the UrQMD model studies that protons freeze-out earlier than pions and kaons [53]. In the hydro calculation of D0 spectra (red in the left panel), we have assumed that D0 mesons freeze out (Tfrz=150 MeV) earlier than the light quark hadrons (Tfrz=137 MeV). If we set the freeze-out temperature of D0 the same as that of the light quark hadrons at 137 MeV with a re-adjusted charm chemical potential μc=800 MeV, the slope of the pT spectra becomes slightly larger (blue) due to increased radial expansion, further away from experimental data. Increasing the freeze-out temperature slightly decreases the slope of the spectra at high pT. However, the change is too small to describe the data at high pT.

      In the left panel of Fig. 4, we further compare the transport calculations (blue bands) of the D0 spectra in Au+Au collisions. Before discussing the results in Au+Au collisions, it is necessary to present the FONLL calculation of the baseline D0 spectra (green band) in p+p collisions. The D0 meson spectra in p+p collisions are computed using the FONLL program [42, 43], for which major uncertainty stems from the variation of the renormalization scale mT/2<μR<2mT, where mT=M2+p2T is the transverse mass of the heavy quark. From Fig. 4, the central prediction from FONLL (green line) systemically underestimates the STAR data on D0 production in p+p collisions at low pT, while the upper band of the prediction with μR=mT/2 is found to be closer to the measurements. Moving to results for nuclear collisions, we can see that the transport approach, which includes effects of non-equilibrium evolution, describes the experimental data on pT spectra at large pT>2 GeV/c. At low pT, the calculations are still consistent with experiments within error bands; however, the central points of the band under-predict the experimental data. We comment that the production of low-pT D mesons in nuclear collisions is sensitive to the small-x region of the nuclear parton distribution function (nPDF), which is less well-constrained than the parton distribution function of the proton. Despite the EPPS16 nPDF publishing its uncertainties, we have only shown calculations using the central fit of the EPPS16 nPDF.

      Comparison of the high-pT part of the spectra between the hydrodynamic results and experimental data shows that charm quarks have not reached complete equilibrium for pT>1 GeV/c, as the spectra from experiments are still below the hydrodynamic results, assuming a fully equilibrated system of charm quarks that flow with the QGP fluid. The enhanced ratios of Λc/D0 and Ds/D0 observed at the RHIC [54, 55] for Au+Au collisions and LHC [56, 57] for Pb+Pb collisions as compared with p+p collisions may partially contribute to this D0 suppression [58-61].

      Both the relativistic hydro and transport models can describe the experimental data on the elliptic flow of D0 mesons within the error bars as shown in the right panel of Fig. 4. One can see the mass ordering of the elliptic flow by comparing it to the elliptic flow of charged pions (blue) from the same hydrodynamic calculations. The hydrodynamic results on v2(pT) at low pT are insensitive to the value of the freeze-out temperature.

      We stress that even though the hydrodynamic model fails to describe the pT spectra of D0 beyond pT>1 GeV/c, its results on the elliptic flow agree with the experimental data and transport model calculation (solid blue) for pT up to 4 GeV/c. This implies that the experimental data on the elliptic flow cannot provide a stringent constraint on the kinetic equilibration of heavy quarks. The strong interaction between heavy quarks and the medium, however, can "drag" the heavy quarks along collective flow developed for the bulk medium even though the interaction might not drive the heavy quarks to full kinetic equilibration.

    IV.   SUMMARY
    • We have calculated the D0 meson spectra and elliptic flow in Au+Au collisions at the RHIC energy within a relativistic viscous hydrodynamic model assuming that charm quarks initially produced through hard processes become fully kinetically equilibrated in the QGP. We neglect the thermal charm quark pair production in the QGP and use an effective charm chemical potential at the freeze-out, which is tuned to describe the number of charm quarks produced in the initial hard processes. With the charm freeze-out temperature Tfrz=137150 MeV, the hydrodynamic model can describe the experimental data on D0 elliptic flow for pT<4 GeV/c as well as the elliptic flow of light quark hadrons. The mass ordering of the elliptic flow for pions and D mesons in the experimental data is also observed for a freeze-out temperature between 137 and 150 MeV. We have also compared the hydrodynamic results with those from a transport model with both energy loss and momentum diffusion, which can describe the observed elliptic flow of D0 mesons at both low and high pT.

      The hydrodynamic model, however, fails to describe the pT spectra of D0 for pT>1 GeV/c, significantly over-predicting the D0 pT spectra at large pT. The transport model, in contrast, can describe the pT spectra well at large pT> 2 GeV/c due to parton energy loss. However, its central values under-predict the spectra at low pT due to the baseline spectra of the initial charm production through hard processes in p+p collisions.

      Our comparative study indicates that only the low-pT D mesons in the experiments might have reached kinetic equilibrium while charm quarks in the intermediate region of pT are partially thermalized. Parton energy loss and momentum diffusion in the transport model can describe well the non-equilibrium behavior of D mesons spectra and elliptic flow at large pT. The interaction of charm quarks in partial equilibrium with the medium in the intermediate pT, however, can develop an elliptic flow as large as that of fully equilibrated charm quarks.

Reference (61)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return