Extracting jet transport coefficient via single hadron and dihadron productions in high-energy heavy-ion collisions

TL;DR

This paper uses a next-to-leading-order perturbative QCD model combined with higher-twist energy loss to analyze jet quenching and extract the jet transport coefficient in high-energy heavy-ion collisions at RHIC and LHC energies.

Contribution

It provides a consistent description of nuclear modification factors for single and dihadron productions and quantitatively extracts the jet quenching parameter using a global analysis.

Findings

Extracted jet quenching parameter $\hat{q}/T^3$ at RHIC and LHC energies.

Provided predictions for dihadron nuclear modification factors at LHC and Xe+Xe collisions.

Results are consistent with JET Collaboration findings.

Abstract

We study the suppressions of high transverse momentum single hadron and dihadron productions in high-energy heavy-ion collisions based on the framework of a next-to-leading-order perturbative QCD parton model combined with the higher-twist energy loss formalism.Our model can provide a consistant description for the nuclear modification factors of single hadron and dihadron productions in central and non-central nucleus-nucleus collisions at RHIC and the LHC energies. We quantitatively extract the value of jet quenching parameter $\hat q$ via a global $\chi^2$ analysis, and obtain ${\hat{q}}/{T^3} = 4.1 \sim 4.4$ at $T = 378$~MeV at RHIC and ${\hat{q}}/{T^3} = 2.6 \sim 3.3$ at $T = 486$~MeV at the LHC, which are consistent with the results from JET Collaboration. We also provide the predictions for the nuclear modification factors of dihadron productions in Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV and in Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}$ = 5.44 TeV.

Equations

$$\frac{d{\sigma }_{pp}^h}{dy{d}^2{p}_T}$$

$$\frac{d{\sigma }_{pp}^{h_1{h}_2}}{dPS}$$

$$\frac{d{N}_{AB}^h}{dy{d}^2{p}_T}$$

$$\frac{d{N}_{AB}^{h_1{h}_2}}{dPS}$$

$$f_{a/A}(x_a,{\mu }^2,\vec{r})$$

$$\displaystyle S_{a/A}(x_{a},\mu^{2},\vec{r})=1+[S_{a/A}(x_{a},\mu^{2})-1]\frac{At_{A}(\vec{r})}{\int{d^{2}}r[t_{A}(\vec{r})]^{2}},\ \ \$$

$${\tilde{D}}_c^h(z_c,{\mu }^2,\Delta E_c)=(1-e^{-⟨N_g⟩})[\frac{z_c^{\prime }}{z_c}{D}_h^c(z_c^{\prime },{\mu }^2)$$

$$+⟨N_g⟩\frac{{z_g}^{\prime }}{z_c}{D}_g^h({z_g}^{\prime },{\mu }^2)]+e^{-⟨N_g⟩}{D}_c^h(z_c,{\mu }^2),$$

$$\frac{\Delta E}{E}$$

$$⟨N_g⟩$$

$$\hat{q}={\hat{q}}_0\frac{T^3}{T_0^3}\frac{p^{\mu }{u}^{\mu }}{p_0},$$

$$R_{AA}=\frac{d{N}_{AB}^h/dy{d}^2{p}_T}{⟨T_{AA}⟩d{\sigma }_{pp}^h/dy{d}^2{p}_T},$$

$$I_{AA}(p_T^{assoc})=\frac{D_{AA}(p_T^{assoc})}{D_{pp}(p_T^{assoc})},$$

$$I_{AA}(z_T)=\frac{D_{AA}(z_T)}{D_{pp}(z_T)},$$

$$\displaystyle D_{AA}(z_{T})=p_{T}^{trig}\frac{d{N}_{AA}^{h_{1}h_{2}}/dy^{trig}dp_{T}^{trig}dy^{assoc}dp_{T}^{assoc}}{\langle T_{AA}\rangle d{\sigma}_{AA}^{h_{1}}/dy^{trig}dp_{T}^{trig}}.\ \ \ \$$

$$\frac{\hat{q}}{T^3}=4.1\sim 4.4,T=378\mathrm{MeV};$$

$$\frac{\hat{q}}{T^3}\approx 3.5,T=469\mathrm{MeV};$$

$$\frac{\hat{q}}{T^3}=2.6\sim 3.3,T=486\mathrm{MeV};$$

$$\frac{\hat{q}}{T^3}\approx 2.5,T=516\mathrm{MeV}.$$

Full text

I Introduction

The strongly-interacting quark-gluon plasma (QGP) can be created in high-energy heavy-ion collisions performed at the Large Hadron Collider (LHC) and the Relativistic Heavy-Ion Collider (RHIC). Jet quenching Gyulassy:1990ye ; Wang:1991xy ; Qin:2015srf has been regarded as an extremely useful tool for studying the properties of such hot and dense nuclear matter. When hard quarks or gluons traverse the QGP matter, they interact with the medium via multiple scatterings and medium-induced gluon radiations. The elastic and inelastic interactions between jet and medium may cause the energy loss of hard jet and also change the energy distribution among jet partons. As one of the consequences of jet quenching and energy loss, the yield of high transverse momentum hadrons fragmented from the surviving hard partons is suppressed as compared to that in proton-proton collisions normalized by the number of binary nucleon-nucleon collisions. Phenomenological studies have been performed on various jet quenching observables, such as the nuclear modifications of single hadron productions Bass:2008rv ; Armesto:2009zi ; Burke:2013yra ; Cao:2017hhk , dihadron and photon-hadron correlations Zhang:2007ja ; Majumder:2004pt ; Qin:2009bk ; Renk:2008xq ; Chen:2016vem ; Chen:2017zte , as well as the observables related to fully reconstructed jets in relativistic nuclear collisions Qin:2010mn ; CasalderreySolana:2010eh ; He:2011pd ; Young:2011qx ; Zapp:2012ak ; Wang:2013cia ; Chang:2016gjp ; Tachibana:2017syd .

In recent years, jet quenching studies have entered the quantitative era in that much effort has been devoted to the quantitative extraction of the so-called jet quenching parameter $\hat{q}$. This parameter is defined as the transverse momentum squared per unit length exchanged between the propagating hard parton and the traversed medium, $\hat{q}=d\langle(\Delta p_{T})^{2}\rangle/dt$, and may be directly related to the gluon density of the nuclear medium Baier:1996sk . Jet transport parameter $\hat{q}$ also controls the amount of medium-induced gluon radiation and thus radiative jet energy loss Baier:1996kr ; Baier:1996sk ; Baier:1998kq ; Guo:2000nz ; Wang:2001ifa ; Majumder:2009ge . In addition, the transverse momentum broadening effect as controlled by $\hat{q}$ may lead to significant nuclear modification on back-to-back dijet, dihadron and other jet-related angular correlations Chen:2016vem ; Chen:2017zte . Among many quantitative jet quenching studies, one of the most important steps is performed by JET Collaboration in Ref. Burke:2013yra which has compared five different theoretical jet quenching models with the nuclear modification data on single hadron productions in most central collisions at RHIC and the LHC and quantitatively extracted the temperature dependence of jet quenching parameter $\hat{q}$. The values of $\hat{q}$ temperatures available at RHIC and the LHC have been obtained as: $\hat{q}/T^{3}=4.6\pm 1.2$ at $T\approx 370$ MeV and $\hat{q}/T^{3}=3.7\pm 1.4$ at $T\approx 470$ MeV for a $10$ GeV quark jet Andres:2016iys . Following this direction, Ref Andres:2016uio ; Andres:2017awo has studied the centrality and collision energy dependence of $\hat{q}$ values at both RHIC and the LHC. Also, Refs. Chen:2016vem has utilized the nuclear modification data on back-to-back dihadron and hadron-jet angular correlations to extract the value of $\hat{q}$ at RHIC.

This paper follows closely the above efforts and study the nuclear modifications of both single hadron and dihadron productions at high transverse momenta using a next-to-leading-order (NLO) perturbative QCD model combined with the higher-twist energy loss formalism. In particular, we perform a global $\chi^{2}$ analysis on the nuclear modification data on single hadron and dihadron productions at RHIC Adare:2008qa ; Adare:2012wg ; Adams:2006yt ; STAR:2016jdz and the LHC Abelev:2012hxa ; CMS:2012aa ; Khachatryan:2016odn ; Acharya:2018qsh ; Acharya:2018eaq ; Adam:2016xbp ; Aamodt:2011vg ; Conway:2013xaa and quantitatively extract the values of jet quenching parameter $\hat{q}$. Our analysis yields ${\hat{q}}/{T^{3}}=4.1\sim 4.4$ at $T=378$ MeV at RHIC and ${\hat{q}}/{T^{3}}=2.6\sim 3.3$ at $T=486$ MeV at the LHC. These results are quantitatively consistent with JET Collaboration. We also extract the $\hat{q}$ values for Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV using the single hadron nuclear modification data, and predict the nuclear modification factors for dihadron productions for these collisions.

Our paper is organized as follows. In Sec. II, we briefly introduce our framework to study the productions of single hadrons and dihadrons at high transverse momenta in proton-proton and nucleus-nucleus collisions. In Sec. III, we perform a global $\chi^{2}$ analysis and extract jet quenching parameter $\hat{q}$ from the nuclear modification data on single hadron and dihadron productions at RHIC and the LHC. We also provide our predictions for the nuclear modification factors of dihadron productions in central and non-central Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV at the LHC. Sec. IV contains our summary.

II Framework

In high-energy proton-proton collisions, the production cross section of high transverse momentum hadrons can be factorized into a convolution of parton distribution functions (PDFs), the cross section of hard partonic scatterings, and fragmentation functions (FFs),

[TABLE]

Here, $f_{a}(x_{a},\mu^{2})$ and $f_{b}(x_{b},\mu^{2})$ are parton distribution functions which we take from CT14 Hou:2016nqm ; $D_{c}^{h}(z_{c},\mu^{2})$ is fragmentation function which we take from Refs. Kretzer:2000yf ; Wang:2004yv ; $d\sigma_{ab\rightarrow cd}/d\hat{t}$ is the tree-level $2\to 2$ partonic scattering cross section. The NLO correction at $\mathcal{O}(\alpha_{s}^{3})$ contains $2\rightarrow 2$ virtual diagrams and $2\rightarrow 3$ tree diagrams, and has been included in our calculation. It has been shown in Ref. Zhang:2007ja that NLO perturbative QCD calculation for single $\pi^{0}$ production in proton-proton collisions agrees well with the experimental data at RHIC.

Similarly, the production cross section for high transverse momentum dihadrons in high-energy proton-proton collisions can be written as,

[TABLE]

where the phase space is $dPS=dy^{h_{1}}d^{2}p_{T}^{h_{1}}dy^{h_{2}}d^{2}p_{T}^{h_{2}}$.

In relativistic nucleus-nucleus collisions, one has to consider both cold nuclear matter effect in the initial state and hot nuclear matter effect in the final state. The yield of single hadron production at high transverse momentum may be obtained as Zhang:2007ja ; Zhang:2009rn ,

[TABLE]

Similarly, the yield of dihadron production at high transverse momentum in nucleus-nucleus collisions may be calculated as Zhang:2007ja ; Zhang:2009rn ; Wang:2003mm

[TABLE]

In the above two equations, $t_{A}(\vec{r})$ is the nuclear thickness function, normalized as $\int d^{2}rt_{A}(\vec{r})=A$, with $A$ the mass number of the nucleus. Here we use the Woods-Saxon form for the nuclear density distribution. $f_{a/A}(x_{a},\mu^{2},\vec{r})$ is the nuclear modified PDF, which we calculate as follows Wang:1996yf ; Li:2001xa :

[TABLE]

where $Z$ is the proton number of the nucleus. Here, $S_{a/A}(x_{a},\mu^{2},\vec{r})$ is called the nuclear shadowing factor and denotes the nuclear modification to the PDF in a free proton $f_{a/p}(x_{a},\mu^{2})$. The shadowing factor $S_{a/A}(x_{a},\mu^{2},\vec{r})$ is calculated using the following form (Emelyanov:1999pkc, ; Hirano:2003pw, ),

[TABLE]

where $S_{a/A}(x_{a},\mu^{2})$ is taken from the EPPS16 Eskola:2016oht . $\tilde{D}_{c}^{h}(z_{c},\Delta E_{c})$ is the medium-modified fragmentation function and is calcualted as follows Wang:2004yv ; Zhang:2007ja ; Zhang:2009rn :

[TABLE]

where $\Delta E_{c}$ is the energy loss of parton $c$, $z_{c}=p_{T}/p_{Tc}$, $z^{\prime}_{c}=p_{T}/(p_{Tc}-\Delta E_{c})$, $z^{\prime}_{g}=\langle{N_{g}}\rangle p_{T}/\Delta{E_{c}}$ and $\langle N_{g}\rangle$ is the average number of gluons radiated by parton $c$. In this work, we use the higher twist formalism Wang:2009qb ; Wang:2001cs ; Wang:2002ri to calculate medium-induced gluon radiation and parton energy loss. For a quark with initial energy $E$, the total energy loss $\Delta E$ can be calculated as,

[TABLE]

where $C_{A}=3$, and $l_{T}$ is the transverse momentum of radiated gluon. We assume the energy loss of a gluon is simply $9/4$ times that of a quark Wang:2009qb . The average number of radiated gluons from the propagating hard parton is calculated as Chang:2014fba ,

[TABLE]

The parton energy loss is controlled by jet transport parameter $\hat{q}$ (Baier:1996sk, ), for which we take the following form:

[TABLE]

where $T$ is the local temperature of the medium, $T_{0}$ is a reference temperature which is usually taken as the temperature at the center of the medium at the hydrodynamics initial time $\tau_{0}=0.6$ fm in central nucleus-nucleus collisions, and $u^{\mu}$ is the four flow velocity of the fluid. In our calculation, the dynamical evolution of the QGP medium is obtained using the OSU (2+1)-dimensional viscous hydrodynamics model (VISH2+1) Song:2007fn ; Song:2007ux ; Qiu:2011hf ; Qiu:2012uy .

III Numerical results

In this section, we present our numerical results for single hadron and dihadron nuclear modification factors in Au+Au collisions at $\sqrt{s_{\rm{NN}}}=0.2$ TeV, Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV and $5.02$ TeV, and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV. A global $\chi^{2}$ analysis is performed to extract the jet quenching parameter $\hat{q}$ in different collision systems and different collision energies at RHIC and the LHC. Based on our analysis, we also provide the predictions for the nuclear modification factors of dihadron productions in Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV.

The nuclear modification factor $R_{AA}$ for single hadron production in heavy-ion collisions is defined as Wang:2004yv ,

[TABLE]

where $T_{AA}(\vec{b})=\int d^{2}rt_{A}(\vec{r})t_{B}(\vec{r}+\vec{b})$ is the overlap function of two colliding nuclei and the average in the equation is taken for a given centrality class.

As for dihadron production at high transverse momentum in heavy-ion collisions, the nuclear modification factor $I_{AA}$ can be defined either as a function of $p_{T}^{assoc}$ or as a function of $z_{T}=p_{T}^{assoc}/p_{T}^{trig}$ Zhang:2007ja

[TABLE]

where $D_{AA}(z_{T})=p_{T}^{trig}D_{AA}(p_{T}^{assoc})$ is called hadron-triggered fragmentation function Wang:2003aw ,

[TABLE]

III.1 Au+Au collisions at $\sqrt{s_{\rm{NN}}}=0.2$ TeV at RHIC

Fig. 1 shows our calculations for single hadron and dihadron nuclear modification factors in central ($0-10\%$) Au+Au collisions at $\sqrt{s_{\rm{NN}}}=0.2$ TeV at RHIC compared with the experimental data taken from PHENIX Adare:2008qa ; Adare:2012wg and STAR STAR:2016jdz Collaborations. In each plot, different lines represent our model calculations for $R_{AA}$ or $I_{AA}$ using different values of jet quenching parameter $\hat{q}_{0}$. The solid line in the middle denotes the result using the best value of $\hat{q}_{0}$ obtained from our global $\chi^{2}$ analysis, which is shown in Fig. 2. In the figure, we also show $\chi^{2}/{\rm d.o.f}$ as a function of $\hat{q}_{0}$ using only $R_{AA}$ data or only $I_{AA}$ data. We can see that two fitting results are consistent with each other. This means that with the similar value of $\hat{q}_{0}$, both single hadron and dihadron nuclear suppression factors can be described consistently within our jet energy loss model. Our global $\chi^{2}$ analysis renders: $\hat{q}_{0}=1.1\sim 1.2$ GeV${}^{2}/$fm at $T_{0}=378$ MeV. In terms of the scaled dimensionless jet quenching parameter, it reads, $\hat{q}/T^{3}=4.1\sim 4.4$ at $T=378$ MeV. These values are consistent with the results obtained by JET Collaboration Burke:2013yra .

To test the goodness of our approach, we use the same $\hat{q}_{0}$ value obtained above to calculate the nuclear modification factors $R_{AA}$ and $I_{AA}$ in mid-central Au+Au collision at $\sqrt{s_{\rm{NN}}}=0.2$ TeV at RHIC. The result is shown in Fig. 3, where the solid lines in the middle denote the results using the best $\hat{q}_{0}$ value (i.e, $\hat{q}_{0}=1.2$ GeV2/fm at $T_{0}=378$ MeV), while the other two lines represent the uncertainty for the extracted $\hat{q}_{0}$ value ($\hat{q}_{0}=1.1$ or $1.3$  GeV2/fm for the two lines). We can see that with the similar $\hat{q}_{0}$ value, our model can provide a good description of experimental data on single and dihadron nuclear modification in both central and non-central Au+Au collisions at $\sqrt{s_{\rm{NN}}}=0.2$ TeV at RHIC.

III.2 Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV at the LHC

Now we present our numerical results for Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV at the LHC. Fig. 4 shows our calculations for single hadron and dihadron nuclear modification factors in central ($0-10\%$) Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV compared with the experimental data from ALICE Abelev:2012hxa ; Aamodt:2011vg ; Adam:2016xbp and CMS CMS:2012aa ; Conway:2013xaa Collaborations. In each plot, different lines represent our model calculations for $R_{AA}$ and (or) $I_{AA}$ using different $\hat{q}_{0}$ values. The solid line in the middle denotes the results using the best $\hat{q}_{0}$ value obtained from our global $\chi^{2}$ analysis. Also show in Fig. 5 is the $\chi^{2}$ analysis of $\hat{q}_{0}$ value using only $R_{AA}$ data or $I_{AA}$ data. Although there is some small difference between two fitting results, they are quantitatively consistent with each other within the uncertainties. Such consistency implies that with the similar values of $\hat{q}_{0}$, our jet energy loss model can provide a consistent description of both single hadron and dihadron nuclear suppression factors in Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV. From Fig. 5, we obtain: $\hat{q}_{0}=1.5\sim 1.9$ GeV${}^{2}/$fm at $T_{0}=486$ MeV, which translates into the scaled jet quenching parameter, $\hat{q}/T^{3}=2.6\sim 3.3$ at $T=486$ MeV. This values are also consistent with JET Collaboration Burke:2013yra .

We also test our approach by using the same $\hat{q}_{0}$ value obtained above to calculate the nuclear modification factors $R_{AA}$ and $I_{AA}$ in the non-central ($50-60\%$) Pb+Pb collision at $\sqrt{s_{\rm{NN}}}=2.76$ TeV. The result is shown in Fig. 6: the solid lines in the middle denote the results using the best $\hat{q}_{0}$ value (i.e., $\hat{q}_{0}=1.6$ GeV2/fm at $T_{0}=486$ MeV), while the other two lines (using $\hat{q}_{0}=1.5$ and $1.9$ GeV2/fm) represent the uncertainty for our extracted $\hat{q}_{0}$ value. We can see that with the same $\hat{q}_{0}$ value, our jet energy loss model can also describe the experimental data on single and dihadron nuclear modification in non-central Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV. Another interesting result is that for both Au+Au collisions at RHIC and Pb+Pb collision at the LHC, the nuclear modification factors $I_{AA}$ for dihadron productions are typically larger than single hadron suppression factors $R_{AA}$ given the same nucleus-nucleus collision conditions. One of the main reasons for such difference is the dominance of tangential emissions in dijet (dihadron) events, as has been been pointed out in Ref. Zhang:2007ja .

III.3 Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV at the LHC

Recently, ALICE Khachatryan:2016odn and CMS Acharya:2018qsh Collaborations have published their measurements on the nuclear modification factor $R_{AA}$ for single hadron productions in Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV. These new results provide a good opportunity for studying the collision energy and system size dependences of jet quenching in relativistic heavy-ion collisions. Since no experimental data on dihadron nuclear modification factor $I_{AA}$ are available for these collisions, we will extract the $\hat{q}_{0}$ values only using the available $R_{AA}$ data. Given that our model can provide a consistent description of both single hadron and dihadron nuclear modifications in Au+Au collisions at $\sqrt{s_{\rm{NN}}}=0.2$ TeV and Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV, we then use the extracted $\hat{q}_{0}$ values to predict dihadron nuclear modification factor $I_{AA}$ in Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV.

Our numerical results are shown in Fig. 7 and Fig. 8, in which the left panels show the result for Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and the right for Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV. Fig. 7 shows the nuclear modification factor $R_{AA}$ (in the upper pannels) together with the $\chi^{2}$ analysis (in the lower pannels). Again, the solid lines are the results using the best fit $\hat{q}_{0}$ values. From our $\chi^{2}$ analysis, we obtain: $\hat{q}_{0}\approx 1.7$ GeV${}^{2}/$fm at $T_{0}=516$ MeV ($\hat{q}/T^{3}\approx 2.5$) for central Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and $\hat{q}_{0}\approx 1.8$ GeV2/fm at $T_{0}=469$ MeV ($\hat{q}/T^{3}\approx 3.5$) for central Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV. Using the extracted $\hat{q}_{0}$ values from fitting $R_{AA}$ data, we present in Fig. 8 our predictions for dihadron nuclear modification factor $I_{AA}$ in central Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV (left) and central Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV (right). Different panels are the results with different transverse momenta for trigger hadrons. In each plot, the solid lines in the middle are the results using the best $\hat{q}_{0}$ values, while the other two lines represent the theoretical which we take $\pm 0.1$ GeV2/fm around the best $\hat{q}_{0}$ values fitted from $R_{AA}$ data. One interesting observation is that as the values of $I_{AA}$ also increase as one increases the trigger hadron transverse momentum. This can be understood since dihadrons with larger transverse momenta are more likely produced by tangential emissions, and thus have smaller nuclear modification effect.

III.4 $\hat{q}$ from single hadron and dihadron nuclear suppressions at RHIC and the LHC

In previous subsections, we have quantitatively extracted the jet quenching parameter $\hat{q}$ by performing a detailed $\chi^{2}$ analysis on the comparison of our jet energy loss model calculations to single hadron and dihadron nuclear modification data at RHIC and the LHC. Here we summarize the main results for the extracted $\hat{q}$ values, in terms of the scaled jet quenching parameter $\hat{q}/T^{3}$:

[TABLE]

One can see that the scaled jet quenching parameter $\hat{q}/T^{3}$ has some temperature dependence: it decreases as one increase the temperature, which may be understood as decreasing jet-medium interaction strength at higher temperature regimes. For better visualization, we also plot the above values in Fig. 9, where the results from JET Collaboration on $\hat{q}/T^{3}$ for Au+Au collisions at $\sqrt{s_{\rm{NN}}}=0.2$ TeV and Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV are also shown. We can see that our extracted values for the scaled jet quenching parameter $\hat{q}/T^{3}$ are consistent with the JET Collaboration results.

IV SUMMARY

In this work, we have studied the nuclear suppressions of single hadron and dihadron productions at high transverse momentum regimes in high-energy heavy-ion collisions at RHIC and the LHC. We compute the cross section of single hadron and dihadron productions in relativistic nuclear collisions based on the NLO perturbative QCD framework. For hadron production in heavy-ion collisions, we include both initial-state cold nuclear matter effect and final-state hot nuclear matter effect. The effect of jet energy loss in hot QGP medium is taken into account using medium-modified fragmentation functions, which are calculated based on the higher-twist formalism. The numerical results from our jet energy loss model calculations show consistent descriptions of the nuclear modifications of single hadron and dihadron productions in central and non-central nucleus-nucleus collisions at RHIC and the LHC.

We have further performed a detailed $\chi^{2}$ analysis by comparing our jet energy loss model calculations with the experimental data on single hadron and dihadron nuclear modifications at RHIC and the LHC. From the global $\chi^{2}$ analysis, we have quantitatively extracted the values of $\hat{q}_{0}$ for different collision systems and collision energies. For Au+Au collisions at $\sqrt{s_{\rm{NN}}}=0.2$ TeV at RHIC, we obtain $\hat{q}_{0}=1.1\sim 1.2$ GeV${}^{2}/$fm at $T_{0}=$ 378 MeV (i.e., ${\hat{q}}/{T^{3}}=4.1\sim 4.4$). For Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV at the LHC, we obtain $\hat{q}_{0}=1.5\sim 1.9$ GeV2/fm at $T_{0}=$ 486 MeV (i.e., ${\hat{q}}/{T^{3}}=2.6\sim 3.3$). These results are consistent with the previous JET Collaboration results. As for Pb+Pb collisions at $\sqrt{s_{\rm{NN}}}=5.02$ TeV and Xe+Xe collisions at $\sqrt{s_{\rm{NN}}}=5.44$ TeV, we have used single hadron $R_{AA}$ data to extract the $\hat{q}$ values. These extracted values are then used to predict the nuclear modification effects in dihadron productions in these collisions. Our work provides an important contribution to our quantitative extraction of the temperature dependence of jet quenching parameter by using multiple jet quenching observables from different collision systems and energies, and is helpful to achieve a consistent understanding of jet quenching in relativistic heavy-ion collisions.

ACKNOWLEDGMENTS

This work is supported in part by Natural Science Foundation of China (NSFC) under grant Nos. 11435004, 11775095, 11890711 and 11375072.