Charm hadron azimuthal angular correlations in Au + Au collisions at sNN = 200 GeV from parton scatterings

Charm hadron azimuthal angular correlations in Au + Au collisions at <math id="M1"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV from parton scatterings 1 2 first-author 3 corresp Corresponding author, chenjinhui@fudan.edu.cn Corresponding author, chenjinhui@fudan.edu.cn chenjinhui@fudan.edu.cn 3 1 3 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China University of Chinese Academy of Sciences, Beijing 100049, China Institute of Modern Physics and Key Laboratory of Nuclear Physics and Ion-beam Application (MOE), Fudan University, Shanghai 200433, China Institute of Modern Physics and Key Laboratory of Nuclear Physics and Ion-beam Application (MOE), Fudan University, Shanghai 200433, China The dynamic evolution of the charm hadron in hot quark matter was studied in the framework of a multiphase transport (AMPT) model. We first reproduced the open charm hadron D 0 p T spectrum in Au + Au collisions at <math id="M2"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV by triggering the <math id="M3"><mrow><mtext>c</mtext><mover accent="true"><mtext>c</mtext><mo>¯</mo></mover></mrow></math> production in AMPT, and then the elliptic flow of charm hadrons was described with different parton cascade cross sections. Charm hadron azimuthal angular correlations were proposed, and they are affected by the different parton cross section parameter applied in the model, which can facilitate our understanding of the loss of collision energy of charm quarks in hot quark medium and can stimulate further experimental studies. Quark gluon plasma Heavy flavor Two particle correlations 1 Introduction 1 Heavy flavor hadrons are important in the studies of the properties of quark–gluon plasma (QGP), which is believed to be created in ultrarelativistic heavy-ion collisions [ 1 1 - 3 3 1-3 ]. The energy–momentum scale of hard probes is sufficiently large to enable the calculation of their initial production rate, the medium modification of the final spectra, and the correlations at high transverse momentum p T based on perturbative quantum chromodynamics (pQCD). These hard probes can provide valuable information on the hot QCD medium. It is suggested that in heavy-ion collisions at the relativistic heavy ion collider (RHIC) and at the large hadron collider (LHC), the massive heavy quarks undergo Brownian-like motion at low momentum, which provide information on the hadronization mechanisms at intermediate momentum, and merge into a radiative energy loss regime at high momentum [ 4 4 , 5 5 ]. Therefore, to study heavy quarks the spectral modification can be used, which results from the interactions between the heavy and light quarks and gluons propagating in their dynamically evolving QCD medium [ 5 5 - 7 7 5-7 ]. The heavy flavor quarks have mass effects, which enable systematic investigations of the variation of the prevalent processes in different p T regions. Compared with light quarks and gluons, the large mass of heavy flavor quarks suppresses the small angle gluon radiation, which results in smaller loss of radiated energy [ 5 5 - 7 7 5-7 ]. The large mass slows down the equilibration speed of the heavy flavor quarks through the medium relative to their light counterparts. Therefore, the non-equilibrated heavy flavor quarks in the final state can provide valuable information on their interaction with the medium throughout their propagation in the QGP medium. For example, charm quarks with significantly larger mass compared to that of light quarks, can achieve a strong collective flow when they diffuse through the QGP [ 8 8 - 11 11 8-11 ]. Experimentally, this scenario was first found in measurements of the semi-leptonic electron decay spectra at RHIC [ 12 12 , 13 13 ], followed by its confirmation by charm quark collectivity measurements at both RHIC and LHC [ 12 12 - 14 14 12-14 ]. The interactions of heavy flavor with the medium can be illustrated by scattering on the medium partons. When hard-scattered partons diffuse through the QGP, they lose energy in the QGP medium. At high p T , the mass effect diminishes and heavy flavor observables degenerate to light flavor. At low p T , a dead-cone effect is expected to arise, which suggests an inverse mass dependence of the energy loss from heavy quarks to light quarks and to gluons [ 15 15 ]. However, the experimental data do not fully support this mechanism. The observations of the heavy flavor nuclear modification factor and elliptic flow are nearly comparable to those of light hadrons [ 16 16 - 19 19 16-19 ], which suggest the importance of the loss of elastic energy in the QGP. Heavy flavor hadrons are ideal for the systematic investigation of the relationship of the radiative and collisional energy loss mechanisms over a broad momentum region and for the identification of the transition between the two [ 20 20 , 21 21 ]. In this study, we used a multi-phase transport (AMPT) model [ 22 22 ] to investigate the dynamics of elastic scattering among partons to understand their propagation into the heavy flavor production and their evolution in the QGP medium. We introduced the <math id="M4"><mrow><mtext>c</mtext><mover accent="true"><mtext>c</mtext><mo>¯</mo></mover></mrow></math> trigger to enhance the charm quark rate in AMPT, in order to reproduce the open charm hadron p T distributions in Au + Au collisions at <math id="M5"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV. Then, the charm hadron azimuthal angular correlations study was used to investigate the collision energy loss dynamics. 2 Method and results 1 2.1 The AMPT model 2 The AMPT model [ 22 22 , 23 23 ] is a transport model consisting of four main components: initial conditions, partonic interactions, conversion from partonic to hadronic matter, and hadronic interactions. The initial conditions, with spatial and momentum distributions of minijet partons and soft string excitations included, are obtained from the heavy ion jet interaction generator (HIJING) model [ 24 24 ], which is an extension of the PYTHIA model [ 25 25 ]. A Woods–Saxon radial shape is used for the colliding gold nuclei and a parametrized nuclear shadowing function that depends on the impact parameter of the collision [ 24 24 ] is introduced. Scatterings among partons are modeled by Zhang’s parton cascade (ZPC) [ 26 26 ], which at present, includes only two-body elastic scatterings with cross sections obtained from pQCD with screening mass. In the default version of the AMPT model, after the partons stop interacting, they recombine with their parent strings, which are produced from initial soft nucleon–nucleon interactions. The resulting strings are converted to hadrons using the Lund string fragmentation model. However, in the case of the string melting version of AMPT model, the hadrons produced from string fragmentation are converted to their valence quarks and antiquarks. The subsequent partonic interactions are modeled by ZPC. Following the freeze-out of the partons, they are recombined into hadrons by a quark coalescence process. The dynamic evolution of the hadronic phase is subsequently described by an extended relativistic transport (ART) model [ 27 27 ] including baryon–baryon, baryon–meson, and meson–meson elastic and inelastic scatterings. Details of the AMPT model can be found in Ref. [ 22 22 ]. 2.2 Quark phase space and charm hadron p T distributions in the AMPT model 2 As we focus on the parton scattering effect on charm hadron azimuthal angular correlation study, in the following, the string melting version of the AMPT model (v2.26t5) is employed [ 23 23 ]. In this version of ZPC, two partons undergo scattering every time when they approach each other with a distance smaller than <math id="M6"><mrow><msqrt><mrow><mi>σ</mi><mo>/</mo><mi>π</mi></mrow></msqrt></mrow></math> and with a total parton elastic scattering cross section of <math id="M7"><mrow><mi>σ</mi><mo>≈</mo><mn>9</mn><mi>π</mi><msubsup><mi>α</mi><mi>s</mi><mn>2</mn></msubsup><mo>/</mo><mo stretchy="false">(</mo><mn>2</mn><msup><mi>μ</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow></math> . This transport process in AMPT simulates the parton energy loss in hot QGP medium [ 22 22 ]. We set the strong coupling constant as α s = 0.33 and screening mass as μ = 2.265/fm or 1.241/fm for σ = 3 mb or 10 mb, respectively [ 22 22 ]. Fig. 1 Fig. 1 shows the momentum space density of charm quarks together with light quarks and the corresponding charm hadrons in 0–80% Au + Au collisions at <math id="M8"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV from AMPT model. Fig. 1(a) Fig. 1(a) shows the quark rapidity density distributions while Fig. 1(b) Fig. 1(b) shows their transverse momentum spectrum distributions. The charm quarks density d N /d y at mid-rapidity is approximately 0.4 and the rapidity density of light quarks is close to 500. We found that this charm quark density is significantly lower than the result at RHIC energies [ 18 18 , 28 28 ]. To improve the prediction power, we enabled the <math id="M13"><mrow><mtext>c</mtext><mover accent="true"><mtext>c</mtext><mo>¯</mo></mover></mrow></math> trigger in HIJING [ 24 24 ], to enhance the total rate of <math id="M14"><mrow><mtext>c</mtext><mover accent="true"><mtext>c</mtext><mo>¯</mo></mover></mrow></math> in our study. The physical process of charm quark evolution is identical to that of other charm quarks in the model. The trigger rate is determined by the experimentally measured cross section. For example, to match the rate at RHIC, we randomly triggered half of the full event sample, and the mid-rapidity density of charm quark reached ∼ 3.6 in the 0–80% central Au + Au collisions at <math id="M15"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV. The AMPT model with similar direction is currently under development [ 29 29 ]. As the inelastic scatterings among partons are not included in the current version of the AMPT model, the quarks produced from the melting of strings scatter uninterruptedly and then freeze out. Those freeze-out charm quarks coalescence with the nearby partons to form hadrons, which continue their evolution in the hadronic medium modeled by ART, and then freeze out, as shown in Fig. 1(c) Fig. 1(c) . With the enhanced production of <math id="M16"><mrow><mtext>c</mtext><mover accent="true"><mi>c</mi><mo>¯</mo></mover></mrow></math> in the early stage of the AMPT process, the result of the D meson p T spectrum from the AMPT model can reasonably reproduce the experimental data. However, the results from the J/ψ p T spectrum are lower than the data. This can be attributed to the fact, that there is no initial production and dissociation of J/ψ particles in the medium. A prediction for Λ c is also presented, which can be found between D s and J/ψ s [cf. Fig. 1(c) Fig. 1(c) ]. Then, the history of charm quark interactions was analyzed in AMPT, by tracing the dynamics of parton cascading by the number of collisions ( N scattering ) of one charm quark with other partons. The procedure described in Ref. [ 30 30 ] was strictly followed, but with a focus on the charm quarks. Fig. 2(a) Fig. 2(a) shows the probability distributions of the charm quarks freeze-out after N scattering collisions. In average, partons are subjected to 〈 N scattering 〉 = 4-5 Au + Au collisions at <math id="M17"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV with σ = 3 mb and impact parameter b = 8 fm, while the number becomes 8–9 for σ = 10 mb. As can be seen in Fig. 2(b) Fig. 2(b) , the average elastic energy loss of a charm quark with p T is close to 2 GeV/ c . The charm quarks continuously lose energy due to the two-body elastic scattering in ZPC. The 〈Δ E 〉 is the largest in the first collision, it drops gradually as more collisions occur, and apparently, it stabilizes with larger 〈 N scattering 〉. It is understood that the 〈Δ E 〉 is smaller in the 10-mb scenario because it typically experiences a larger number of scatterings. Fig. 2(c) Fig. 2(c) shows the average energy loss of charm quarks as a function of their momentum. These values are expected to be similar in the range of σ = 3–10 mb, and their difference is due to the probability distribution of N scattering . The 〈ΔE〉 increases as the momentum of the charm quark increases, and its value is close to that of the calculation based on a linearized Boltzmann transport model [ 20 20 ]. 2.3 Two-particle angular correlation 2 To study the azimuthal angle dependence of the charm evolution in hot dense medium simulated in AMPT, two-particle angular correlations are used, which is a powerful tool to determine the interaction of the lower energy jet (or parton) with the surrounding medium. The analysis process is described in detail in Refs. [ 31 31 - 34 34 31-34 ], and we briefly introduce the method in this section. Selected particles from each event are paired for correlations as <math display="block" id="M18"><mrow><mi>S</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>η</mi><mn>,</mn><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msup><mtext>d</mtext><mn>2</mn></msup><msubsup><mi>N</mi><mrow><mtext>pairs</mtext></mrow><mrow><mtext>signal</mtext></mrow></msubsup></mrow><mrow><mtext>d</mtext><mi>Δ</mi><mi>η</mi><mtext>d</mtext><mi>Δ</mi><mi>ϕ</mi></mrow></mfrac><mn>,</mn></mrow></math> while they are combined with particles from different events to build the background distribution <math display="block" id="M19"><mrow><mi>B</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>η</mi><mn>,</mn><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msup><mtext>d</mtext><mn>2</mn></msup><msubsup><mi>N</mi><mrow><mtext>pairs</mtext></mrow><mrow><mtext>mixed</mtext></mrow></msubsup></mrow><mrow><mtext>d</mtext><mi>Δ</mi><mi>η</mi><mtext>d</mtext><mi>Δ</mi><mi>ϕ</mi></mrow></mfrac><mn>,</mn></mrow></math> where Δϕ is the relative azimuthal angle and Δ≤ta is the relative pseudorapidity between the particle pair. < > In our analysis, each event is mixed with 10 other events to improve the statistical power of the background estimation, while the direction of the impact parameter of the collisions in the AMPT events is rotated randomly in the transverse plane to calculate B (Δ≤ta,Δϕ). %We also carry out another check by mixing event with similar event plane direction, and the difference between different background reconstructions is found to be negligible. Then, the two-particle correlation function can be obtained as <math display="block" id="M20"><mrow><mi>C</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>η</mi><mn>,</mn><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mi>S</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>η</mi><mn>,</mn><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo></mrow><mrow><mi>B</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>η</mi><mn>,</mn><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo></mrow></mfrac><mo>×</mo><mfrac><mrow><msubsup><mi>N</mi><mrow><mtext>pairs</mtext></mrow><mrow><mtext>mixed</mtext></mrow></msubsup></mrow><mrow><msubsup><mi>N</mi><mrow><mtext>pairs</mtext></mrow><mrow><mtext>signal</mtext></mrow></msubsup></mrow></mfrac><mn>.</mn></mrow></math> A one-dimensional Δϕ correlation function can be constructed from the C (Δ η ,Δ ϕ ) by integrating over Δ≤ta as <math display="block" id="M21"><mrow><mi>C</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>A</mi><mo>×</mo><mfrac><mrow><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mi>S</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>η</mi><mn>,</mn><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>Δ</mi><mi>η</mi></mrow><mrow><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mi>B</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>η</mi><mn>,</mn><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>Δ</mi><mi>η</mi></mrow></mfrac><mn>,</mn></mrow></math> where the normalization constant A is given by <math id="M22"><mrow><msubsup><mi>N</mi><mrow><mtext>pairs</mtext></mrow><mrow><mtext>mixed</mtext></mrow></msubsup><mo>/</mo><msubsup><mi>N</mi><mrow><mtext>pairs</mtext></mrow><mrow><mtext>signal</mtext></mrow></msubsup></mrow></math> . The distribution of pairs in Δϕ can be expanded to a Fourier series, <math display="block" id="M23"><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>N</mi><mrow><mtext>p</mtext><mi>a</mi><mi>i</mi><mi>r</mi><mi>s</mi></mrow></msub></mrow><mrow><mtext>d</mtext><mi>Δ</mi><mi>ϕ</mi></mrow></mfrac><mo>∝</mo><mn>1</mn><mo>+</mo><mn>2</mn><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>∞</mi></munderover></mstyle><msub><mi>v</mi><mrow><mi>n</mi><mn>,</mn><mi>n</mi></mrow></msub><mo stretchy="false">(</mo><msubsup><mi>p</mi><mtext>T</mtext><mn>1</mn></msubsup><mn>,</mn><msubsup><mi>p</mi><mtext>T</mtext><mn>2</mn></msubsup><mo stretchy="false">)</mo><mi>cos</mi><mo stretchy="false">(</mo><mi>n</mi><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mn>.</mn></mrow></math> The coefficients v n , n can be directly calculated by <math display="block" id="M24"><mrow><msub><mi>v</mi><mrow><mi>n</mi><mn>,</mn><mi>n</mi></mrow></msub><mo>=</mo><mo>〈</mo><mi>cos</mi><mo stretchy="false">(</mo><mi>n</mi><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo><mo>〉</mo><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><mi>cos</mi><mo stretchy="false">(</mo><mi>n</mi><mi>Δ</mi><msub><mi>ϕ</mi><mi>m</mi></msub><mo stretchy="false">)</mo><mi>C</mi><mo stretchy="false">(</mo><mi>Δ</mi><msub><mi>ϕ</mi><mi>m</mi></msub><mo stretchy="false">)</mo></mrow><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><mi>C</mi><mo stretchy="false">(</mo><mi>Δ</mi><mi>ϕ</mi><mo stretchy="false">)</mo></mrow></mfrac><mn>,</mn></mrow></math> where n = 2,3,4 and N = 200 is the number of Δϕ bins. The harmonic flow coefficients v n ( n = 2,3,4) can be calculated as <math id="M25"><mrow><msub><mi>v</mi><mi>n</mi></msub><mo>=</mo><msub><mi>v</mi><mrow><mi>n</mi><mn>,</mn><mi>n</mi></mrow></msub><mo>/</mo><msqrt><mrow><mo>|</mo><msub><mi>v</mi><mrow><mi>n</mi><mn>,</mn><mi>n</mi></mrow></msub><mo>|</mo></mrow></msqrt></mrow></math> . 2.4 Azimuthal correlation of charm quarks and charm hadrons 2 The azimuthal anisotropy of charm and light hadrons can be obtained from the two-particle angular correlations as described above. Fig. 3 Fig. 3 shows the elliptic flow v 2 of mid-rapidity hadrons as a function of p T in 0–80% Au + Au collisions at <math id="M26"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV, for π s ( Fig. 3(a) Fig. 3(a) ), <math id="M27"><mrow><msubsup><mi>K</mi><mi>S</mi><mn>0</mn></msubsup></mrow></math> s ( Fig. 3(b) Fig. 3(b) ), protons ( Fig. 3(c) Fig. 3(c) ), and charm hadrons ( Fig. 3(d) Fig. 3(d) ). Systematically, the AMPT model calculations with a parton cross section of 3 mb generate a smaller elliptic flow than those from 10 mb, and they describe the low p T data better than results from 10 mb, as shown in Fig. 3 Fig. 3 . When p T increases to 1.5 GeV/ c and above, results from the AMPT model with a parton cross section of 3 mb underestimate the data. From this p T range, results from the AMPT model with a parton cross section of 10 mb are closer to the data. As shown in Fig. 2 Fig. 2 , results in the range of 3––10 mb represent the different number of parton collisions and they are responsible for the difference in the v 2 values. This is also the case for charm hadrons [cf. Fig. 2(d) Fig. 2(d) ]. For the D hadron result, as data is only available from p T gt; 1.0 GeV/ c , the AMPT result with a cross section of 10 mb gives better description of the D 0 data. Calculations on Λ c and D s are also available. They follow a similar p T and parton cross section dependence behavior as those other particles presented in panels (a-d) of Fig. 3 Fig. 3 . Further measurements of the elliptic flow of Λ c and D s at RHIC energies can help to distinguish the parton cross section dependence of v 2 ( p T ) behavior. Then the charm–charm azimuthal correlations were studied in Au + Au collisions at <math id="M29"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV. Fig. 4 Fig. 4 shows the quarks azimuthal correlations, which was performed at the stage of parton freeze-out in ZPC, prior to the hadronization in AMPT. Partons with p T lt; 2.5 GeV/ c were chosen to cover the majority of quarks in the collisions, which required a η gap of (0.8,2.4) between quarks. It should be noted that the kinematic window needs to be adjusted with different acceptance in different experiments. In the case of the total parton cross section of 3 mb ( Fig. 4(a) Fig. 4(a) ), the charm–charm correlation is suppressed at the near side and enhanced at the far side, which is different from the distributions of light quark azimuthal correlations. This is possibly due to a stronger effect of the initial production on the massive charm quarks, as indicated by the open triangles in Fig. 4 Fig. 4 and discussed in Ref. [ 36 36 ]. The increase in the total parton cross section enhances the collision probability among quarks [cf. Fig. 2 Fig. 2 ]; thus, it reduces the difference from the effect of initial production between the charm and light quarks. This can also be seen in Fig. 4(b) Fig. 4(b) , with a larger total parton cross section calculation, where the correlations between the charm quarks are nearly the same as light quarks. These charm quarks coalesce into a charm hadron with the nearby quarks and undergo hadronic interaction in the AMPT. Fig. 5 Fig. 5 shows the corresponding charm hadron azimuthal correlations together with light hadrons. In the case of the total parton cross section of 3 mb, the correlations between D – π are nearly identical to those of π – π , suggesting hadronization and hadronic scattering effect on the evolution of charm hadrons, while the correlations between D – D are slightly different from those of light hadrons. They have a slightly higher distribution on the far side and a lower distribution on the near side. According to the study with the larger total parton cascade cross section of 10 mb (panel Fig. 4(b) Fig. 4(b) ), the correlations between charm hadrons and between light hadrons are the same. Fig. 5(a) Fig. 5(a) shows a calculation of D – D correlations without parton interaction (the 0-mb cross section scenario), which gives a flat distribution along Δϕ. This is a scenario with hadronization and hadronic interaction only, which is similar to the default version of the AMPT model. The comparison of D – D azimuthal correlations among different parton cross section parameters applied in AMPT suggests that the number of parton collisions affects the evolution of charm quarks in QGP medium. We also studied the high p T charm hadron azimuthal correlations in QGP medium. Figure 6 Figure 6 shows the two-particle azimuthal correlations with a trigger particle p T gt; 5 GeV/ c and its associated particle at p T 2 GeV/ c . The correlations are normalized with respect to the number of triggers. The results of charge hadron azimuthal correlations from AMPT partially describe the experimental data, which are slightly different in parton cross sections of 3 and 10 mb, because only elastic scattering is included in the current calculations. For the D – D correlations with small parton cross section, similar to the low p T azimuthal correlation results, shown in of Fig. 5(a) Fig. 5(a) , the high p T D – D correlations exhibit different behavior than the light flavor charge hadrons, in this case with lower yield on both the near and far sides. For the case of 10-mb parton cross section, correlations between D – D are similar to light charge hadrons on the far side, and lower on the near side, as shown in Fig. 5(b) Fig. 5(b) . 3 SUMMARY 1 We studied the charm quark evolution in QGP medium created in ultrarelativistic heavy ion collisions with the AMPT model. By including an additional <math id="M33"><mrow><mtext>c</mtext><mover accent="true"><mi>c</mi><mo>¯</mo></mover></mrow></math> production in the AMPT to reproduce the open charm hadron p T spectrum in 0–80% Au + Au collisions at <math id="M34"><mrow><msqrt><mrow><msub><mi>s</mi><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV, we obtained a reasonable description of the elliptic flow of the D meson, and then predictions were provided for other charm hadrons including D s , Λ c , and J/ψ . We also studied the azimuthal angular correlations between D mesons. We found that the total parton cross section presents a clear effect on the D meson azimuthal correlations. By combining the different charm quark average energy loss with different parton cascade cross section parameters, our studies provide an effective method to understand the collisional energy loss of charm quarks in hot QGP medium. Therefore, further experimental measurements from LHC and sPHENIX are required. References E. Shuryak , Strongly coupled quark-gluon plasma in heavy ion collisions . Rev. Mod. Phys . 89 , 035001 ( 2017 ). doi: 10.1103/RevModPhys.89.035001 http://doi.org/10.1103/RevModPhys.89.035001 Strongly coupled quark-gluon plasma in heavy ion collisions J.H. Chen , D. Keane and Y.G. Ma et al ., Antinuclei in heavy-ion collisions . Phys. Rep . 760 , 1 ( 2018 ). doi: 10.1016/j.physerp.2018.07.002 http://doi.org/10.1016/j.physerp.2018.07.002 Antinuclei in heavy-ion collisions H.C. Song , Y. Zhou , K. Gajdosova , Collective flow and hydrodynamics in large and small systems at the LHC . NUCL SCI TECH 28 , 99 ( 2017 ). doi: 10.1007/s41365-017-0245-4 http://doi.org/10.1007/s41365-017-0245-4 Collective flow and hydrodynamics in large and small systems at the LHC L. Adamczyk , J.K. Adkins and G. Agakishiev et al ., Measurement of D0 azimuthal anisotropy at midrapidity in Au+Au collisions at <math id="M35"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 118 , 212301 ( 2017 ). doi: 10.1103/PhysRevLett.118.212301 http://doi.org/10.1103/PhysRevLett.118.212301 Measurement of D0 azimuthal anisotropy at midrapidity in Au+Au collisions at sNN = 200 GeV X. Dong , Y.J. Lee and R. Rapp , Open Heavy-Flavor Production in Heavy-Ion Collisions . arXiv:1903.07709 Open Heavy-Flavor Production in Heavy-Ion Collisions S.S. Cao , G. Coci and S.K. Das et al ., Toward the determination of heavy-quark transport coefficients in quark-gluon plasm . Phys. Rev. C 99 , 054907 ( 2019 ). doi: 10.1103/PhysRevC.99.054907 http://doi.org/10.1103/PhysRevC.99.054907 Toward the determination of heavy-quark transport coefficients in quark-gluon plasm A. Andronic , F. Arleo and R. Arnaldi et al ., Heavy-flavour and quarkonium production in the LHC era: from proton–proton to heavy-ion collisions . Eur. Phys. J. C 76 , 107 ( 2016 ). doi: 10.1140/epjc/s10052-015-3819-5 http://doi.org/10.1140/epjc/s10052-015-3819-5 Heavy-flavour and quarkonium production in the LHC era: from proton–proton to heavy-ion collisions H. Hees , V. Greco and R. Rapp , Heavy-quark probes of the quark-gluon plasma and interpretation of recent data taken at the BNL Relativistic Heavy Ion Collider . Phys. Rev. C 73 , 034913 ( 2006 ). DOI: 10.1103/PhysRevC.73.034913 http://doi.org/10.1103/PhysRevC.73.034913 Heavy-quark probes of the quark-gluon plasma and interpretation of recent data taken at the BNL Relativistic Heavy Ion Collider P.B. Gossiaux , R. Bierkandt and J. Aichelin , Tomography of quark gluon plasma at energies available at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC) . Phys. Rev. C 79 , 044906 ( 2016 ). doi: 10.1103/PhysRevC.79.044906 http://doi.org/10.1103/PhysRevC.79.044906 Tomography of quark gluon plasma at energies available at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC) M. He , R.J. Fries and R. Rapp , Heavy-quark diffusion and hadronization in quark-gluon plasma . Phys. Rev. C 86 , 014903 ( 2012 ). doi: 10.1103/PhysRevC.86.014903 http://doi.org/10.1103/PhysRevC.86.014903 Heavy-quark diffusion and hadronization in quark-gluon plasma S. Plumari , V. Minissale and S.K. Das et al ., Charmed hadrons from coalescence plus fragmentation in relativistic nucleus-nucleus collisions at RHIC and LHC . Eur. Phys. J. C 78 , 348 ( 2018 ). doi: 10.1140/epjc/s10052-018-5828-7 http://doi.org/10.1140/epjc/s10052-018-5828-7 Charmed hadrons from coalescence plus fragmentation in relativistic nucleus-nucleus collisions at RHIC and LHC S.S. Adler , S. Afanasiev and C. Aidala et al ., Nuclear modification of electron spectra and implications for heavy quark energy loss in Au+Au collisions at <math id="M36"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 96 , 032301 ( 2006 ). doi: 10.1103/PhysRevLett.96.032301 http://doi.org/10.1103/PhysRevLett.96.032301 Nuclear modification of electron spectra and implications for heavy quark energy loss in Au+Au collisions atsNN = 200 GeV A. Adare , S. Afanasiev and C. Aidala et al ., Energy loss and flow of heavy quarks in Au+Au collisions at <math id="M37"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 98 , 172301 ( 2007 ). doi: 10.1103/PhysRevLett.98.172301 http://doi.org/10.1103/PhysRevLett.98.172301 Energy loss and flow of heavy quarks in Au+Au collisions at sNN = 200 GeV B.I. Abelev , M.M. Aggarwal and Z. Ahammed et al ., Transverse momentum and centrality dependence of high-pT nonphotonic electron suppression in Au+Au collisions at <math id="M38"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 98 , 192301 ( 2007 ). doi: 10.1103/PhysRevLett.98.192301 http://doi.org/10.1103/PhysRevLett.98.192301 Transverse momentum and centrality dependence of high-pT nonphotonic electron suppression in Au+Au collisions at sNN = 200 GeV A. Buzzatti and M. Gyulassy , Jet flavor tomography of quark gluon plasmas at RHIC and LHC . Phys. Rev. Lett . 108 , 022301 ( 2012 ). doi: 10.1103/PhysRevLett.108.022301 http://doi.org/10.1103/PhysRevLett.108.022301 Jet flavor tomography of quark gluon plasmas at RHIC and LHC B. Abelev , J. Adam and D. Adamova et al ., Suppression of high transverse momentum D mesons in central Pb-Pb collisions at <math id="M39"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 2.76 TeV . J. High Energy Phys . 2012 , 112 (2012) . doi: 10.1007/JHEP09(2012)112 http://doi.org/10.1007/JHEP09(2012)112 Suppression of high transverse momentum D mesons in central Pb-Pb collisions at sNN = 2.76 TeV B. Abelev , J. Adam and D. Adamova et al ., D Meson Elliptic Flow in Noncentral Pb-Pb Collisions at <math id="M40"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 2.76 TeV . Phys. Rev. Lett . 111 , 102301 ( 2013 ). doi: 10.1103/PhysRevLett.111.102301 http://doi.org/10.1103/PhysRevLett.111.102301 D Meson Elliptic Flow in Noncentral Pb-Pb Collisions at sNN = 2.76 TeV L. Adamczyk , J.K. Adkins and G. Agakishiev et al ., Observation of D0 meson nuclear modifications in Au+Au collisions at <math id="M41"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 113 , 142301 ( 2014 ). doi: 10.1103/PhysRevLett.113.142301 http://doi.org/10.1103/PhysRevLett.113.142301 Observation of D0 meson nuclear modifications in Au+Au collisions at sNN = 200 GeV A. Adare , C. Aidala and N.N. Ajitanand et al ., Heavy-quark production and elliptic flow in Au+Au collisions at <math id="M42"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 62.4 GeV . Phys. Rev. C 91 , 044907 ( 2015 ). doi: 10.1103/PhysRevC.91.044907 http://doi.org/10.1103/PhysRevC.91.044907 Heavy-quark production and elliptic flow in Au+Au collisions at sNN = 62.4 GeV S.S. Cao , T. Luo and G.Y. Qing et al ., Linearized Boltzmann transport model for jet propagation in the quark-gluon plasma: Heavy quark evolution . Phys. Rev. C 94 , 014909 ( 2016 ). doi: 10.1103/PhysRevC.94.014909 http://doi.org/10.1103/PhysRevC.94.014909 Linearized Boltzmann transport model for jet propagation in the quark-gluon plasma: Heavy quark evolution M. He , R. J. Fries and R. Rapp , Heavy flavor at the large hadron collider in a strong coupling approach . Phys. Lett. B 735 , 445 ( 2014 ). doi: 10.1016/j.physletb.2014.05.050 http://doi.org/10.1016/j.physletb.2014.05.050 Heavy flavor at the large hadron collider in a strong coupling approach Z. W. Lin , C. M. Ko and B. A. Li et al ., Multiphase transport model for relativistic heavy ion collisions . Phys. Rev. C 72 , 064901 ( 2005 ). doi: 10.1103/PhysRevC.72.064901 http://doi.org/10.1103/PhysRevC.72.064901 Multiphase transport model for relativistic heavy ion collisions AMPT source codes are available at http://myweb.ecu.edu/linz/ampt/ http://myweb.ecu.edu/linz/ampt/ AMPT source codes are X. N. Wang and M. Gyulassy , HIJING: A Monte Carlo model for multiple jet production in pp, p A, and A A collisions . Phys. Rev. D 44 , 3501 ( 1991 ). doi: 10.1103/PhysRevD.44.3501 http://doi.org/10.1103/PhysRevD.44.3501 HIJING: A Monte Carlo model for multiple jet production in pp, p A, and A A collisions http://home.thep.lu.se/torbjorn/Pythia.html http://home.thep.lu.se/torbjorn/Pythia.html B. Zhang , Erratum to "ZPC 1.0.1" a parton cascade for ultrarelativistic heavy ion collisions" . Comput. Phys. Commun . 109 , 193 - 206 ( 1998 ). doi: 10.1016/S0010-4655(98)00010-1 http://doi.org/10.1016/S0010-4655(98)00010-1 Erratum to "ZPC 1.0.1" a parton cascade for ultrarelativistic heavy ion collisions" B. A. Li and C. M. Ko , Formation of superdense hadronic matter in high energy heavy-ion collisions . Phys. Rev. C 52 , 2037 ( 1995 ). doi: 10.1103/PhysRevC.52.2037 http://doi.org/10.1103/PhysRevC.52.2037 Formation of superdense hadronic matter in high energy heavy-ion collisions L. Adamczyk , J.K. Adkins and G. Agakishiev et al ., J/Ψ production at low pT in Au+Au and Cu+Cu collisions at <math id="M43"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV with the STAR detector . Phys. Rev. C 90 , 024906 ( 2014 ). doi: 10.1103/PhysRevC.90.024906 http://doi.org/10.1103/PhysRevC.90.024906 J/Ψ production at low pT in Au+Au and Cu+Cu collisions at sNN = 200 GeV with the STAR detector C. Zhang , L. Zheng and F. Liu et al ., Update of a multiphase transport model with modern parton distribution functions and nuclear shadowing . Phys. Rev. C 99 , 064906 ( 2019 ). doi: 10.1103/PhysRevC.99.064906 http://doi.org/10.1103/PhysRevC.99.064906 Update of a multiphase transport model with modern parton distribution functions and nuclear shadowing L. He , T. Edmonds and Z.W. Lin et al ., Anisotropic parton escape is the dominant source of azimuthal anisotropy in transport models . Phys. Lett. B 753 506 ( 2016 ). doi: 10.1016/j.physletb.2015.12.051 http://doi.org/10.1016/j.physletb.2015.12.051 Anisotropic parton escape is the dominant source of azimuthal anisotropy in transport models G. Aad , B. Abbott and J. Abdallah et al ., Measurement of the azimuthal anisotropy for charged particle production in <math id="M44"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 2.76 TeV lead-lead collisions with the ATLAS detector . Phys. Rev. C 86 , 014907 ( 2012 ). doi: 10.1103/PhysRevC.86.014907 http://doi.org/10.1103/PhysRevC.86.014907 Measurement of the azimuthal anisotropy for charged particle production in sNN = 2.76 TeV lead-lead collisions with the ATLAS detector S. Chatrchyan , V. Khachatryan and A.M. Sirunyan et al ., Multiplicity and transverse momentum dependence of two- and four-particle correlations in pPb and PbPb collisions . Phys. Lett. B 724 213 ( 2013 ). doi: 10.1016/j.physletb.2013.06.028 http://doi.org/10.1016/j.physletb.2013.06.028 Multiplicity and transverse momentum dependence of two- and four-particle correlations in pPb and PbPb collisions L. Y. Zhang , J. H. Chen and Z. W. Lin et al ., Two-particle angular correlations in pp and p-Pb collisions at energies available at the CERN Large Hadron Collider from a multiphase transport model . Phys. Rev. C 98 , 034912 ( 2018 ). doi: 10.1103/PhysRevC.98.034912 http://doi.org/10.1103/PhysRevC.98.034912 Two-particle angular correlations in pp and p-Pb collisions at energies available at the CERN Large Hadron Collider from a multiphase transport model L. Y. Zhang , J. H. Chen and Z. W. Lin et al ., Two-particle angular correlations in heavy ion collisions from a multiphase transport model . Phys. Rev. C 99 , 054904 ( 2019 ). doi: 10.1103/PhysRevC.99.054904 http://doi.org/10.1103/PhysRevC.99.054904 Two-particle angular correlations in heavy ion collisions from a multiphase transport model L. Adamczyk , J.K. Adkins and G. Agakishiev et al ., Centrality and transverse momentum dependence of elliptic flow of multistrange hadrons and ϕ meson in Au+Au collisions at <math id="M45"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. Lett . 116 , 062301 ( 2016 ). doi: 10.1103/PhysRevLett.116.062301 http://doi.org/10.1103/PhysRevLett.116.062301 Centrality and transverse momentum dependence of elliptic flow of multistrange hadrons and ϕ meson in Au+Au collisions at sNN = 200 GeV X. Zhu , N. Xu and P. Zhuang , Effect of partonic wind on charm quark correlations in high-energy nuclear collisions . Phys. Rev. Lett 100 , 152301 ( 2008 ). doi: 10.1103/PhysRevLett.100.152301 http://doi.org/10.1103/PhysRevLett.100.152301 Effect of partonic wind on charm quark correlations in high-energy nuclear collisions A. Adare et al ., Measurement of two-particle correlations with respect to second- and third-order event planes in Au+Au collisions at <math id="M46"><mrow><msqrt><mrow><msub><mtext>s</mtext><mrow><mtext>NN</mtext></mrow></msub></mrow></msqrt></mrow></math> = 200 GeV . Phys. Rev. C 99 , 054903 ( 2019 ). doi: 10.1103/PhysRevC.99.054903 http://doi.org/10.1103/PhysRevC.99.054903 Measurement of two-particle correlations with respect to second- and third-order event planes in Au+Au collisions at sNN = 200 GeV We are grateful to Dr. Chen Zhong, who maintains the high-performance computer center, where the calculations of this study were performed. 10.1007/s41365-019-0706-z This work was partly supported by the Key Research Program of the Chinese Academy of Science (No. XDPB09) and the National Natural Science Foundation of China (Nos. 11890714, 11775288, 11421505, and 11520101004). 2019-09-19 2019-10-18 2019-10-28 2019-11-23 2019-12-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2019 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-019-0706-z http://dx.doi.org/10.1007/s41365-019-0706-z . 2019 Nuclear Science and Techniques 12 30