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

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Charm hadron azimuthal angular correlations in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from parton scatterings

Hai Wang，

Jin-Hui Chen，

Yu-Gang Ma，

Song Zhang

Nuclear Science and Techniques

Vol.30, No.12

Article number 185

Published in print 01 Dec 2019

Available online 23 Nov 2019

DOI：10.1007/s41365-019-0706-z

41200

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₀ p_T spectrum in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV by triggering the $c \bar{c}$ 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 plasmaHeavy flavorTwo particle correlations

1 Introduction

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-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, 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-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-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-11]. Experimentally, this scenario was first found in measurements of the semi-leptonic electron decay spectra at RHIC [12, 13], followed by its confirmation by charm quark collectivity measurements at both RHIC and LHC [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]. 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-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, 21]. In this study, we used a multi-phase transport (AMPT) model [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 $c \bar{c}$ 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 $\sqrt{s_{NN}}$ = 200 GeV. Then, the charm hadron azimuthal angular correlations study was used to investigate the collision energy loss dynamics.

2 Method and results

2.1 The AMPT model

The AMPT model [22, 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], which is an extension of the PYTHIA model [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] is introduced. Scatterings among partons are modeled by Zhang’s parton cascade (ZPC) [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] including baryon–baryon, baryon–meson, and meson–meson elastic and inelastic scatterings. Details of the AMPT model can be found in Ref. [22].

2.2 Quark phase space and charm hadron p_T distributions in the AMPT model

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]. In this version of ZPC, two partons undergo scattering every time when they approach each other with a distance smaller than $\sqrt{σ / π}$ and with a total parton elastic scattering cross section of $σ \approx 9 π α_{s}^{2} / (2 μ^{2})$ . This transport process in AMPT simulates the parton energy loss in hot QGP medium [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]. 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 $\sqrt{s_{NN}}$ = 200 GeV from AMPT model. Fig. 1(a) shows the quark rapidity density distributions while Fig. 1(b) shows their transverse momentum spectrum distributions. The charm quarks density dN/dy at mid-rapidity is approximately 0.4 and the rapidity density of light quarks is close to 500.

Fig. 1.

(Color online) (a) Rapidity density of charm quarks

(c + \bar{c})

and strange quarks

(s + \bar{s})

as well as light quarks

(u + \bar{u} + d + \bar{d})

at parton freeze-out in Au + Au collisions at

\sqrt{s_{NN}}

= 200 GeV. (b) and (c) shows p_T distributions, for quarks (b) and for the corresponding charm hadrons (c) compared with experimental data [18, 28].

We found that this charm quark density is significantly lower than the result at RHIC energies [18, 28]. To improve the prediction power, we enabled the $c \bar{c}$ trigger in HIJING [24], to enhance the total rate of $c \bar{c}$ 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 $\sqrt{s_{NN}}$ = 200 GeV. The AMPT model with similar direction is currently under development [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). With the enhanced production of $c \bar{c}$ 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 Ds and J/ψs [cf. 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] was strictly followed, but with a focus on the charm quarks. 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 $\sqrt{s_{NN}}$ = 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), 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) 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].

Fig. 2.

(a) Normalized probability distributions of charm quark freezing-out after N_scattering collisions in AMPT with cross sections of 3 mb and 10 mb, respectively. (b) Corresponding average energy loss of charm quarks with initial transverse momentum of 1.6 lt; p_T lt; 2.0 GeV/c as a function of their number of collisions in ZPC. (c) Average energy loss of charm quarks in different momentum windows.

2.3 Two-particle angular correlation

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-34], and we briefly introduce the method in this section. Selected particles from each event are paired for correlations as

S (Δ η, Δ ϕ) = \frac{d^{2} N_{pairs}^{signal}}{d Δ η d Δ ϕ},

(1)

while they are combined with particles from different events to build the background distribution

B (Δ η, Δ ϕ) = \frac{d^{2} N_{pairs}^{mixed}}{d Δ η d Δ ϕ},

(2)

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

C (Δ η, Δ ϕ) = \frac{S (Δ η, Δ ϕ)}{B (Δ η, Δ ϕ)} \times \frac{N_{pairs}^{mixed}}{N_{pairs}^{signal}} .

(3)

A one-dimensional Δϕ correlation function can be constructed from the C(Δη,Δϕ) by integrating over Δ≤ta as

C (Δ ϕ) = A \times \frac{\int S (Δ η, Δ ϕ) d Δ η}{\int B (Δ η, Δ ϕ) d Δ η},

(4)

where the normalization constant A is given by $N_{pairs}^{mixed} / N_{pairs}^{signal}$ . The distribution of pairs in Δϕ can be expanded to a Fourier series,

\frac{d N_{p a i r s}}{d Δ ϕ} \propto 1 + 2 \sum_{n = 1}^{\infty} v_{n, n} (p_{T}^{1}, p_{T}^{2}) \cos (n Δ ϕ) .

(5)

The coefficients v_n,n can be directly calculated by

v_{n, n} = 〈 \cos (n Δ ϕ) 〉 = \frac{\sum_{m = 1}^{N} \cos (n Δ ϕ_{m}) C (Δ ϕ_{m})}{\sum_{m = 1}^{N} C (Δ ϕ)},

(6)

where n = 2,3,4 and N = 200 is the number of Δϕ bins. The harmonic flow coefficients vn (n = 2,3,4) can be calculated as $v_{n} = v_{n, n} / \sqrt{| v_{n, n} |}$ .

2.4 Azimuthal correlation of charm quarks and charm hadrons

The azimuthal anisotropy of charm and light hadrons can be obtained from the two-particle angular correlations as described above. Fig. 3 shows the elliptic flow v₂ of mid-rapidity hadrons as a function of p_T in 0–80% Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, for πs (Fig. 3(a)), $K_{S}^{0}$ s (Fig. 3(b)), protons (Fig. 3(c)), and charm hadrons (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. 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, 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₂ values. This is also the case for charm hadrons [cf. 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₀ data. Calculations on Λc and Ds 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. Further measurements of the elliptic flow of Λc and Ds at RHIC energies can help to distinguish the parton cross section dependence of v₂(p_T) behavior.

Fig. 3.

(Color online) Elliptic flow v₂ of mid-rapidity hadrons as a function of p_T in 0–80% Au + Au collisions at

\sqrt{s_{NN}}

= 200 GeV. The calculations were performed with total parton cross sections of 3 and 10 mb. The experimental data were obtained from Refs.[4, 35].

Then the charm–charm azimuthal correlations were studied in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. 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)), 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 and discussed in Ref. [36]. The increase in the total parton cross section enhances the collision probability among quarks [cf. 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), with a larger total parton cross section calculation, where the correlations between the charm quarks are nearly the same as light quarks.

Fig. 4.

(Color online) Two-particle azimuthal correlations between charm quarks (open symbols), charm quarks and light quarks (blue stars), and light quarks (red squares) in 0–80% Au + Au collisions at

\sqrt{s_{NN}}

= 200 GeV from AMPT: the total parton cross section of 3 mb (a) and 10 mb (b). A scenario of charm–charm azimuthal correlations without parton interaction is also shown in (a).

These charm quarks coalesce into a charm hadron with the nearby quarks and undergo hadronic interaction in the AMPT. 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)), the correlations between charm hadrons and between light hadrons are the same. 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.

Fig. 5.

(Color online) Two-particle azimuthal correlations of D–D, D–π, and π-π in 0–80% Au + Au collisions at

\sqrt{s_{NN}}

= 200 GeV from AMPT: total parton cross section of 3 mb (a) and 10 mb (b). A scenario of D–D correlations without parton interaction is plotted for reference in panel (a).

We also studied the high p_T charm hadron azimuthal correlations in QGP medium. 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), 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. 6.

(Color online) The two-particle azimuthal correlations of D-D, D-h^±, and h^±-h^± in 0–80% Au + Au collisions at

\sqrt{s_{NN}}

= 200 GeV, with a trigger particle of p_T gt; 5 GeV/c and associated particle of p_T 2 GeV/c. Experimental data on charge hadron azimuthal correlation is obtained from Ref. [37].

3 SUMMARY

We studied the charm quark evolution in QGP medium created in ultrarelativistic heavy ion collisions with the AMPT model. By including an additional $c \bar{c}$ production in the AMPT to reproduce the open charm hadron p_T spectrum in 0–80% Au + Au collisions at $\sqrt{s_{NN}}$ = 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 Ds, Λ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

[1]

E. Shuryak,

Strongly coupled quark-gluon plasma in heavy ion collisions

. Rev. Mod. Phys. 89, 035001 (2017). doi: 10.1103/RevModPhys.89.035001