![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-9ec99d26-99ec-4d77-b47e-ba4210df3b97.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-f7f1cee1-15af-48df-98cc-a6c046c73374.jpg)

Nuclear Physics A

Nuclear Physics A 00 (2014) 1-4

## Heavy Flavor Dynamics in QGP and Hadron Gas

## Shanshan Cao a , Guang-You Qin b , Ste ff en A Bass a

a

Department of Physics, Duke University, Durham, NC 27708, USA

b Institute of Particle Physics and Key Laboratory of Quark and Lepton Physics (MOE), Central China Normal University, Wuhan, 430079, China

## Abstract

We study heavy flavor evolution in the quark-gluon plasma matter and the subsequent hadron gas created in ultrarelativistic heavyion collisions. The motion of heavy quarks inside the QGP is described using our modified Langevin framework that incorporates both collisional and radiative energy loss mechanisms; and the scatterings between heavy mesons and the hadron gas are simulated with the UrQMD model. We find that the hadronic interaction further suppresses the D meson R AA at high p T and enhances its v 2 . And our calculations provide good descriptions of experimental data from both RHIC and LHC. In addition, we explore the heavy-flavor-tagged angular correlation functions and find them to be a potential candidate for distinguishing di ff erent energy loss mechanisms of heavy quarks inside the QGP.

Keywords: heavy flavor, energy loss, hadronic interaction, correlation function

## 1. Introduction

Heavy flavor hadrons serve as excellent probes of the quark-gluon plasma (QGP) matter created in ultrarelativistic heavy-ion collisions. Experimental observations at both RHIC and LHC have revealed many interesting data on heavy flavor mesons and their decay electrons especially their small values of R AA and large values of v 2 which are comparable to those of light hadrons [1, 2, 3, 4]. This seems inconsistent with earlier expectation of the mass hierarchy of parton energy loss and thus requires a thorough understanding of heavy flavor dynamics.

Various theoretical frameworks have been developed to study the heavy flavor such as the Boltzmann-based parton cascade model [5], the linearized Boltzmann transport [6], and the Langevin approach [7]. In this work, we construct a framework that describes the full-time evolution of heavy flavor in heavy-ion collisions. The motion of heavy quarks inside the QGP is described with our improved Langevin equation that incorporates both collisional and radiative energy loss [8, 7], the hadronization process is calculated with our hybrid model of fragmentation plus coalescence [7], and the hadronic interaction is simulated with the UrQMD model [9]. With this framework, we not only provide good descriptions of single heavy flavor spectra, but also explore the angular correlation functions of heavy flavor pairs and find them sensitive to di ff erent energy loss mechanisms of heavy quarks inside the QGP.

## 2. Evolution of Heavy Flavor in Heavy-ion Collisions

Since most heavy quarks are produced in the primordial stage of collisions, we use a Monte-Carlo Glauber model to initialize their positions and pQCD calculations for their momenta. The CTEQ parametrization is adopted for the parton distribution functions which are further modified by the EPS09 parametrization to include the nuclear shadowing e ff ect. The bulk matter (QGP) - is initialized with either the Glauber or the KLN model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-d5117141-e152-4753-89e4-f029a965d78f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-f7975ab0-47e5-443e-9f16-a22d1fd0fe69.jpg)

T

Figure 1. (Color online) D meson suppression (left) and flow (right) in the LHC experiment.

In the QGP stage, the space-time evolution of the bulk matter is simulated by a (2 + 1)-D viscous hydrodynamic model [10] and the in-medium evolution of heavy quarks is described using our improved Langevin equation [7]:

$$\frac { d \vec { p } } { d t } = - \eta _ { D } ( p ) \vec { p } + \vec { \xi } + \vec { f } _ { g }.$$

The first two terms on the right are the drag and the thermal forces describing the multiple scatterings of heavy quarks o ff light partons in the medium background, and are related by the fluctuation-dissipation theorem η D ( p ) = κ/ (2 TE ). Here the momentum space di ff usion coe ffi cient κ is defined in 〈 ξ i ( ) t ξ j ( t ′ ) 〉 = κδ i j δ ( t -t ′ ). The third term /vector f g = -dpg /vector / dt is introduced as the recoil force exerted on the heavy quark while it radiates gluons. The probability of gluon radiation and the energy-momentum distribution of the emitted gluons are calculated according to the following gluon distribution function taken from the Higher-Twist energy loss formalism [11]:

$$\frac { d N _ { g } } { d x d k _ { \perp } ^ { 2 } d t } & = \frac { 2 \alpha _ { s } ( k _ { \perp } ) } { \pi } P ( x ) \frac { \hat { q } } { k _ { \perp } ^ { 4 } } \sin ^ { 2 } \left ( \frac { t - t _ { i } } { 2 \tau _ { f } } \right ) \left ( \frac { k _ { \perp } ^ { 2 } } { k _ { \perp } ^ { 2 } + x ^ { 2 } M ^ { 2 } } \right ) ^ { 4 }, \\ \text{$\text{$\text{$non$confract $\nu$ is the fractional energy carried by radiated $\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\math\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\math\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\mathrm{\math\mathrm{\mathrm{\mathrm{\mathrm{\math\mathrm{\math\mathrm{\math\mathrm{\mathrm{\mathrm{\math\mathrm{\math\mathrm{\math\mathrm{\mathrm{\math\mathrm{\mathrm{\math\mathrm{\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\math\$ is its transverse}$}$}$}$}$}$}$}$}$}$}$}$$

where ˆ is the gluon transport coe q ffi cient, x is the fractional energy carried by radiated gluon, and k ⊥ is its transverse momentum. We relate the spatial di ff usion coe ffi cient of heavy quark D with κ and ˆ via q D = 2 T 2 /κ and ˆ q = 2 κ CA / CF , so that in end only one free parameter exists. It has been shown in [7] that while the collisional energy loss dominates the low p T region of heavy quark, gluon radiation dominates at high p T.

After the heavy quarks traverse the QGP medium, they hadronize based on our hybrid model of fragmentation plus coalescence. And it has been shown in [7] that while the fragmentation mechanism dominates heavy meson formation at high p T, the heavy-light quark coalescence significantly enhances its production at medium p T. The decay of the bulk matter is calculated according to the Cooper-Frye formula.

Finally, both heavy mesons formed from heavy quarks and soft hadrons decayed from the QGP are fed into the UrQMD model for the hadronic interactions. One of the most important ingredients of the UrQMD model are the hadronic scattering cross sections. To simulate the interactions of D mesons with the hadron gas, we introduce the scattering cross sections of charm mesons with pions and rho mesons calculated in Ref. [12] into UrQMD.

We show our calculations of the D meson R AA in the left panel of Fig.1 for the LHC experiments. We observe that due to additional energy loss that the D meson su ff ers inside the hadron gas, its R AA is further decreased at large p T. The error band corresponds to the uncertainty in choosing the cuto ff parameter of the hadron form factor for calculating scattering cross sections [12]. Our calculations provide a good description of the observed D meson R AA with an extraction of the gluon transport coe ffi cient ˆ q of around 2.6 GeV 2 / fm. In the right panel, we find that additional scatterings of D mesons in an anisotropic hadron gas further enhance the v 2 value by over 30%. Two di ff erent initializations of hydrodynamical evolution are also compared which may cause another 30% di ff erence in the D meson v 2 . In Fig.2, we present our calculations for the RHIC experiments, the above conclusions remain and our results are consistent with the experimental data. Centrality and participant number dependences of D meson R AA can be found in our presentation as well.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-6edead6a-3fff-429d-8401-bc0634bb4c9e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-6a228959-db16-4b91-9b94-911ecd29cbd1.jpg)

Figure 2. (Color online) D meson suppression (left) and flow (right) in the RHIC experiment.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-da47f68d-a813-416a-9f15-e2d843a929fc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-a73994a9-715a-4d20-86e8-90fe17b8d074.jpg)

∆φ

Figure 3. (Color online) Fitting D meson R AA with di ff erent energy loss mechanisms by tuning the di ff usion coe ffi cient.

Figure 4. (Color online) Comparison of the angular correlation functions of cc ¯pairs between di ff erent energy loss mechanisms.

## 3. Angular Correlation Functions of Heavy Flavor Pairs

We have shown good descriptions of D meson suppression by including both collisional and radiative energy losses. However, in Fig.3, one sees that by tuning the transport coe ffi cient, collisional or radiative energy loss alone also provide reasonable description of the data within the current experimental uncertainties. This implies the insu ffi -ciency of single particle spectra and motivate us to explore exclusive spectra of heavy flavor [13, 14, 15].

In Fig.4, we investigate the angular correlation functions of cc ¯pairs after they traverse a QGP medium. We start with the leading order approximation of initial c and ¯pairs and find that after they travel through the medium, the c correlation function still peaks around π if only gluon radiation is considered. However, a peak near 0 is observed if one only includes the collisional energy loss. This indicates that unlike the back-to-back initial state, low energy cc ¯ pairs tend to move collinearly in the end because of the boost by the radial flow of the medium.

In Fig.5, we use an improved initialization for cc ¯pairs - the Monte-Carto next-to-leading-order production plus Herwig vacuum radiation [13]. Similar to Fig.4, pure gluon radiation does not a ff ect the angular correlation function, but pure scattering leads to a peak around 0. A more realistic analysis is implemented in Fig.6 where each D meson is looped over all ¯ within an event. We observe that the shapes of D D -¯ correlation functions resemble those of D cc ¯ pairs except for the existence of a large background of uncorrelated D and ¯ . D One may also find D -hadron correlation functions in our presentation though they may not as clean as those directly between heavy flavor pairs. Although these functions are dependent on di ff erent models, they may provide a deeper insight into heavy quark energy loss mechanisms if comparisons can be made between theory and future experiments.

Figure 5. (Color online) Angular correlation functions of cc ¯ pairs with MCNLO + Herwig initial conditions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-8cc1fbc7-ee36-42bb-9afc-1f351aea1f3b.jpg)

Figure 6. (Color online) Angular correlation functions of D D - ¯ with MCNLO Herwig initial conditions. +

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0503v1-figure-c6ab8a12-13ad-4ada-9fbf-60cd7b79747c.jpg)

## 4. Conclusions

We have developed a framework to describe the full-time evolution of heavy flavor in heavy-ion collisions. Within our framework, we have provided the D meson suppression and flow consistent with most existing data, and explored the heavy-flavor-tagged angular correlation functions and found them useful in distinguishing between di ff erent energy loss mechanisms of heavy quarks inside the QGP.

We thank the cooperations within the JET Collaboration. This work was supported by the U.S. Department of Energy Grant No. DE-FG02-05ER41367 and Natural Science Foundation of China (NSFC) under grant No. 11375072.

## References

- [1] L. Adamczyk, et al., Observation of D 0 meson nuclear modifications in Au + Au collisions at √ s NN = 200 GeV. arXiv:1404.6185 .
- [3] B. Abelev, et al., Suppression of high transverse momentum D mesons in central Pb-Pb collisions at √ s NN = 2 76 TeV, JHEP 1209 (2012) . 112. arXiv:1203.2160 , doi:10.1007/JHEP09(2012)112 .
- [2] A. Adare, et al., Heavy Quark Production in p + p and Energy Loss and Flow of Heavy Quarks in Au + Au Collisions at √ s NN = 200 GeV, Phys. Rev. C84 (2011) 044905. arXiv:1005.1627 , doi:10.1103/PhysRevC.84.044905 .

[4] B. Abelev, et al., D meson elliptic flow in non-central Pb-Pb collisions at √ s NN = 2 76TeV, . Phys. Rev. Lett. 111 (2013) 102301. arXiv:1305.2707 , doi:10.1103/PhysRevLett.111.102301 .

- [5] J. Upho ff , O. Fochler, Z. Xu, C. Greiner, Open Heavy Flavor in Pb + Pb Collisions at √ s = 2 76 TeV within a Transport Model, Phys. Lett. . B717 (2012) 430-435. arXiv:1205.4945 , doi:10.1016/j.physletb.2012.09.069 .

[6] P. Gossiaux, J. Aichelin, T. Gousset, V . Guiho, Competition of Heavy Quark Radiative and Collisional Energy Loss in Deconfined Matter, J. Phys. G37 (2010) 094019. arXiv:1001.4166 , doi:10.1088/0954-3899/37/9/094019 .

- [7] S. Cao, G.-Y. Qin, S. A. Bass, Heavy quark dynamics and hadronization in ultra-relativistic heavy-ion collisions: collisional versus radiative energy loss, Phys. Rev. C88 (2013) 044907. arXiv:1308.0617 , doi:10.1103/PhysRevC.88.044907 .
- [8] S. Cao, G.-Y. Qin, S. A. Bass, B. M¨ller, Collisional vs.Radiative Energy Loss of Heavy Quark in a Hot and Dense Nuclear Matter, Nucl. u Phys. A904-905 (2013) 653c-656c. arXiv:1209.5410 , doi:10.1016/j.nuclphysa.2013.02.100 .

[9] S. A. Bass, et al., Microscopic models for ultrarelativistic heavy ion collisions, Prog. Part. Nucl. Phys. 41 (1998) 225-370. arXiv:nucl-th/9803035 .

[10] Z. Qiu, C. Shen, U. Heinz, Hydrodynamic elliptic and triangular flow in Pb-Pb collisions at sqrt(s) = 2.76ATeV, Phys. Lett. B707 (2012) 151-155. arXiv:1110.3033 , doi:10.1016/j.physletb.2011.12.041 .

[11] B.-W. Zhang, E. Wang, X.-N. Wang, Heavy quark energy loss in nuclear medium, Phys. Rev. Lett. 93 (2004) 072301. arXiv:nucl-th/0309040 , doi:10.1103/PhysRevLett.93.072301 .

[12] Z.-W. Lin, T. Di, C. Ko, Charm meson scattering cross-sections by pion and rho meson, Nucl. Phys. A689 (2001) 965-979. arXiv:nucl-th/0006086 , doi:10.1016/S0375-9474(00)00611-4 .

[13] M. Nahrgang, J. Aichelin, P. B. Gossiaux, K. Werner, Azimuthal correlations of heavy quarks in Pb + Pb collisions at LHC ( √ s = 2 76 TeV). . arXiv:1305.3823 .

[14] S. Cao, G.-Y. Qin, S. A. Bass, B. M¨ller, Heavy quark ene rgy loss and angular de-correlation in a quark-gluon plasma matter, J. Phys. Conf. u Ser. 446 (2013) 012035. doi:10.1088/1742-6596/446/1/012035 .

[15] S. Cao, G.-Y. Qin, S. A. Bass, Dynamical Evolution, Hadronization and Angular De-correlation of Heavy Flavor in a Hot and Dense QCD Medium. arXiv:1404.1081 .
Shanshan Cao Guang-You Qin Steffen Bass
Heavy Flavor Dynamics in QGP and Hadron Gas
We study heavy flavor evolution in the quark-gluon plasma matter and the
subsequent hadron gas created in ultrarelativistic heavy-ion collisions. The
motion of heavy quarks inside the QGP is described using our modified Langevin
framework that incorporates both collisional and radiative energy loss
mechanisms; and the scatterings between heavy mesons and the hadron gas are
simulated with the UrQMD model. We find that the hadronic interaction further
suppresses the $D$ meson $R_\mathrm{AA}$ at high $p_\mathrm{T}$ and enhances
its $v_2$. And our calculations provide good descriptions of experimental data
from both RHIC and LHC. In addition, we explore the heavy-flavor-tagged angular
correlation functions and find them to be a potential candidate for
distinguishing different energy loss mechanisms of heavy quarks inside the QGP.