Jets, Mach cone, hot spots, ridges, harmonic flow, dihadron and
$\gamma$-hadron correlation in high-energy heavy-ion collisions

Show preview

Hide preview

ar X
iv :1
01 1.
52 49
v3 [
nu cl-
th] 2
1 A pr
20 11
Jets, Mach cone, hot spots, ridges, harmonic flow, dihadron and γ-hadron correlation
in high-energy heavy-ion collisions
Guo-Liang Ma1 and Xin-Nian Wang2, 3
1Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China 2Institute of Particle Physics, Huazhong Normal University, Wuhan 430079, China
3Nuclear Science Division MS 70R0319, Lawrence Berkeley National Laboratory, Berkeley, California 94720
Within the AMPT Monte Carlo model, fluctuations in the initial transverse parton density are shown to lead to harmonic flows. The net back-to-back dihadron azimuthal correlation after sub- traction of contributions from harmonic flows still has a double-peak that is independent of the initial geometric triangularity and unique to the jet-induced Mach cone and expanding hot spots distorted by radial flow. The longitudinal structure of hot spots also leads to a near-side ridge in dihadron correlation with a large rapidity gap. By successively randomizing the azimuthal angle of the transverse momenta and positions of initial partons, one can isolate the effects of jet-induced medium excitation and expanding hot spots on the dihadron azimuthal correlation. The double- peaks in the net dihadron and γ-hadron correlation are quantitatively different since the later is caused only by jet-induced Mach cone.
PACS numbers: 25.75.-q, 25.75.Bh,25.75.Cj,25.75.Ld
Dihadron correlation in high-energy heavy-ion exper- iments provides a useful tool to study jet quenching [1] that suppresses not only single inclusive large pT hadrons [2] but also large pT back-to-back correlation [3]. By decreasing pT of the away-side hadrons associated with a high pT triggered hadron, one can further study the medium modification of hadron distribution from an en- ergetic jet as well as the jet-induced medium excitations. Experimental data in Au + Au collisions at the Rela- tivistic Heavy-Ion Collider (RHIC) [4–6] show a double- peaked back-to-back azimuthal dihadron correlation with a maximum opening angle of ∆φ ≈ 1 (rad) relative to the back-side of a high pT trigger. Such a feature in dihadron correlation was suspected to be caused by Mach-cone ex- citation induced by jet-medium interaction [7, 8]. How- ever, hydrodynamic study of jet-induced medium excita- tion with a realistic source term for energy-momentum deposition [9] and parton transport study of jet propaga- tion [10] both show that it is the deflection of jet show- ers and medium excitation by radial flow that produces double-peaked back-to-back dihadron correlation. More- over, recent studies also found that expanding hot spots in the initial parton transverse density distorted by radial flow [11, 12] and the triangular flow [13, 14] all contribute to a double-peaked back-to-back dihadron correlation. If these local fluctuations are extended in the longitudinal direction like a string model [15] or glasma [16, 17], they would also lead to a ridge structure in the near-side di- hadron correlation with large rapidity gap [18, 19]. With these different mechanisms contributing to the dihadron correlation, it is important to explore ways to separate different contributions and study the characteristics of the dihadron correlation from each of them.
In this Letter, we will use a multiphase transport (AMPT) model [20] to study dihadron correlation as a result of harmonic flows, expanding hot spots, jets and jet-induced medium excitation. We first study harmonic flows of hadron spectra due to fluctuations in initial par-
ton density, which all contribute to dihadron azimuthal correlation. The dihadron correlation after subtraction of contributions from harmonic flows should provide in- formation on the unique structures of hadron correlation from medium modified jets, jet-induced medium excita- tion and expanding hot spots distorted by radial flow in high-energy heavy-ion collisions. We will also investigate the longitudinal structures of hot spots and the result- ing dihadron azimuthal correlation with large rapidity gap. By successively randomizing the azimuthal angles of transverse momenta and coordinates of initial jet shower partons, we can isolate the effects of medium modified dijets, jet-induced medium excitation and expanding hot spots. We will also compare dihadron correlation after subtraction of the contributions from harmonic flow with γ-hadron correlation which is only affected by jet-induced medium excitation. The AMPT results shown in this paper are central Au + Au collisions with fixed impact- parameter b = 0 at the RHIC energy
√ s = 200 GeV.
AMPT [20] model combines initial conditions from HI- JING model [22] with parton and hadron cascade for final state interaction. We will use a version of AMPT with string melting and parton recombination for hadroniza- tion which was shown to better describe the collective phenomenon in heavy-ion collisions at RHIC [21]. The parton cascade employed in AMPT includes only elas- tic parton collisions whose cross sections are controlled by the values of strong coupling constant and the Debye screening mass. Within HIJING, Glauber model for mul- tiple nucleon scattering is used to describe the initial par- ton production in heavy-ion collisions. Nucleon-nucleon scatterings contain both independent hard parton scat- tering and coherent soft interaction that is modeled by string formation for each participant nucleon. Strings are then converted into soft partons via string melting scheme in AMPT. Such multiple scatterings lead to fluc- tuation in local parton number density or hot spots from both soft and hard interactions which are proportional

2 x (fm)-8 -6 -4 -2 0 2 4 6 8
y (fm
)
-8
-6
-4
-2
0
2
4
6
8 =0.022∈ =0.023∈
x (fm)-8 -6 -4 -2 0 2 4 6 8 -6
-4
-2
0
2
4
6 -1 < y < 1 fm η
-8
FIG. 1: (Color online) Contour plot of initial parton density (in arbitrary unit) dN/dxdy in transverse plane (left panel) and dN/dxdη (right panel) in x-η (pseudorapidity) plane in a typical AMPT central Au+ Au event (b = 0) at
√ s = 200
GeV, with ellipticity ǫ2 = 0.02 and triangularity ǫ3 = 0.02 of the transverse parton distribution.
to local transverse density of participant nucleons and binary collisions, respectively.
Shown in Fig. 1 are contour plots of initial parton den- sity in the transverse plane dN/dxdy (left panel) and in x − η(pseudorapidity) plane dN/dxdη within a slice |y| < 1 fm (right panel) of a typical AMPT event. We assume the transverse position of both hard and soft partons are confined to the size of their parent nucle- ons rN ∼ 1 fm. Therefore, fluctuations in local density of participant nucleons and binary nucleon-nucleon col- lisions lead to fluctuations in parton transverse density or hot spots with the smallest transverse size of about 1 fm. These hot spots and valleys give rise to finite values of initial geometric irregularities ǫn defined as
ǫn =
√ 〈r2 cos(nϕ)〉2 + 〈r2 sin(nϕ)〉2
〈r2〉 , (1)
where (r, ϕ) are polar coordinates of each parton and the average 〈· · · 〉 is density weighted. Normally, ǫ2 is referred to as eccentricity and ǫ3 as triangularity.
Hot spots in the fluctuating initial parton density distribution are also extended in the longitudinal di- rection as shown in the right panel of Fig. 1. Such extended longitudinal distribution in pseudorapidity η is partially from soft partons via the materialization of strings. Partons from initial state radiation associ- ated with hard scatterings have a distribution dN/dy = dN/d log(1/x) ∼ 1 which is also extended in rapidity. Collective expansion due to parton rescattering will
translate the initial geometric irregularities into har- monic flows in momentum space [23]. Shown in Fig. 2 are the harmonic flows of final hadron spectra vn = 〈cosn(φ− ψn)〉 from AMPT calculations, where the
(GeV/c) T
p 0 0.5 1 1.5 2 2.5 3 3.5 4
n v
0
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
1 v
2 v
3 v
4 v
5 v
6 v
Au+Au 200 GeV (b = 0 fm)
FIG. 2: (Color online) Azimuthal anisotropies of hadron spec- tra vn(pT ) (n = 1 − 6) in central (b = 0) Au + Au collisions at √ s = 200 GeV from AMPT model calculation.
event plane angle for each harmonics is given by
ψn = 1
n
[ arctan
⟨ r2 sin(nϕ)
⟩ 〈r2 cos(nϕ)〉 + π
] . (2)
Even in central collisions at fixed impact-parameter b = 0, the geometrical fluctuation of initial parton density leads to anisotropic collective expansion which translates the initial geometric irregularities into significant values of harmonic flows in momentum space in the central ra- pidity region. It is interesting to observe that all vn(pT ) decrease at high pT but the turning points shift to higher pT for higher harmonics. The higher harmonic flows vn for n > 6 are insignificant due to viscous diffusion. In the study of dihadron correlation to search for effects
of medium modification of jet structure and jet-induced medium excitations, it is important to isolate and sub- tract contributions from harmonic flows, especially the triangular flow which contributes the most to the double- peak structure of back-to-back dihadron correlation. We will use the ZYAM (zero yield at minimum) scheme [4] to subtract contributions to dihadron correlation,
f(∆φ) = B
( 1 +
∞∑ n=1
2〈vtrign vasson 〉 cosn∆φ ) , (3)
from harmonic flows, where B is a normalization fac- tor determined by the ZYAM scheme, vtrign and v
asso n
are harmonic flow coefficients for trigger and associated hadrons. Shown in Fig. 3 are dihadron correlations be- fore (dot-dashed) and after (solid) the removal of contri-
butions from harmonic flows for ptrigT > 2.5 GeV/c and 1 < passoT < 2 GeV/c. Also shown are contributions from each harmonic flow n = 2-6 (dashed). These contribu- tions are significant for up to n = 5 harmonics. After subtraction of contributions from harmonic
flows, the dihadron correlation still has a double-peak feature on the away-side of the trigger which should re- flect the azimuthal structure from deflection of medium

3 (rad)φ∆ -1 0 1 2 3 4 5
φ∆ )d
N/ d
tr ig
(1/ N
-1
0
1
2
3
4
> 2.5 GeV/c,trig T
p Au+Au 200 GeV (b=0 fm)
constant ZYAM BG removed flow ZYAM BG removed
contributions (n=2-6)nv
< 2 GeV/casso T
1 < p
FIG. 3: (Color online) AMPT results on dihadron azimuthal correlation before (dot-dashed) and after (solid) subtraction
of contribution from harmonic flow vn(n = 2− 6) for ptrigT > 2.5 GeV/c and 1 < passoT < 2 GeV/c.
(rad)φ∆-1 0 1 2 3 4 5
φ∆ )d
N/ d
tr ig
(1/ N
0
0.2
0.4
0.6
0.8
1 < 0.053∈0.00 < < 0.153∈0.05 < < 0.303∈0.15 <
Au+Au 200 GeV (b=0 fm) < 2 GeV/casso
T > 2.5 GeV/c, 1 < ptrig
T p
FIG. 4: (Color online) Dihadron correlations after subtraction of harmonic flow with different values of geometric triangu- larity ǫ3 for p
trig
T > 2.5 GeV/c and 1 < passoT < 2 GeV/c
.
modified dijets and jet-induced medium excitations and hadrons from expanding hot spots under strong radial flow. The structure therefore should be unique and insen- sitive to the fluctuation of the initial geometry of dense matter at a fixed impact-parameter. As shown in Fig. 4, the dihadron correlations after subtraction of contribu- tions from harmonic flows become independent on the initial geometric triangularity ǫ3. In order to study the structure of dihadron azimuthal
correlation from jets (including jet-induced medium ex- citation) and hot spots separately, we successively switch off each mechanism and calculate the dihadron corre- lation within AMPT model. We first randomize the azimuthal angle of each jet shower parton in the ini- tial condition from HIJING simulations. This effectively switches off the initial back-to-back correlation of di-
(rad)φ∆-1 0 1 2 3 4 5
φ∆ )/d
N/ d
tr ig
(1/ N
0
0.2
0.4
0.6
0.8
1 full hot spots soft hot spots smoothed
> 2.5 GeV/ctrig T
p
< 2 GeV/casso T
1 < p
Au+Au 200 GeV (b=0 fm)
FIG. 5: (Color online) Dihadron correlation (with harmonic flow subtracted) from AMPT with different initial conditions
for ptrig T
> 2.5 GeV/c and 1 < passoT < 2 GeV/c. See text for details on the different initial conditions.
(rad)φ∆-1 0 1 2 3 4 5
η∆
-4
-3
-2
-1
0
1
2
3
4
0
0.1
0.2
0.3
0.4
0.5
FIG. 6: (Color online) The contour plot of dihadron cor- relation (with harmonic flow subtracted) from AMPT in
azimuthal and pseudorapidity for ptrig T
> 2.5 GeV/c and 1 < passoT < 2 GeV/c.
jets. The dihadron correlation (dashed) denoted as “hot spots” in Fig. 5 still exhibits a double-peak on the away- side that comes only from hot spots. It has roughly the same opening angle ∆φ ∼ 1 (rad) as in the “full” simula- tion (solid). However, the magnitude of the double-peak in the away-side correlation is reduced by about 40%, which can be attributed to dihadrons from medium mod- ified dijets and jet-induced medium excitation. Without knowing the relative yields of hadrons from jets and hot spots, it is difficult to extract dihadron correlation from dijets (and jet-induced medium excitation) alone. When

4 (rad)φ∆ -1 0 1 2 3 4 5
φ∆ )d
N/ d
tr ig
(1/ N
0
0.5
1
1.5
2 di-hadron (flow BG removed)
-hadronγ
< 2 GeV/casso T
1 < p
Au+Au 200 GeV (b=0 fm)
(h) > 4 GeV/ctrig T
p
) ~ 15 GeV/cγ(trig T
p
FIG. 7: Dihadron correlation (solid) compared with γ-hadron
correlation (dashed) from AMPT for ptrig T
(h) > 4 GeV/c,
ptrig T
(γ) ≥ 15 GeV/c and 1 < passoT < 2 GeV/c .
jet production is turned off in the HIJING initial con- dition, fluctuation in soft partons from strings can still form what we denote as “soft hot spots” that lead to a back-to-back dihadron correlation (dot-dashed) with a weak double-peak. It is clear that jet shower partons in- crease the local parton density in “hot spots” and lead to a stronger double-peaked dihadron correlation than that of “soft hot spots”. Such “soft hot spots” are the likely candidate mechanism for the observed back-to-back di- hadron correlation in heavy-ion collisions at the SPS en- ergy [24] where jet production is insignificant. Without jets in AMPT, one can further randomize the polar an- gle of transverse coordinates of soft partons and therefore eliminate the “soft hot spots”. The dihadron correlation from such “smoothed” initial condition becomes almost flat (dotted). Because of the elongated shape of hot spots in the lon-
gitudinal direction as shown in Fig. 1 (left panel), the structure of dihadron azimuthal correlation should re-
main similar with a large range of rapidity gaps as shown in Fig. 6. Since dihadron correlation from a single jet is only restricted to small rapidity gap δη ∼ 1, the near-side azimuthal correlation with a large rapidity gap, referred to as ridge, can only come from a hot spot. This is con- firmed by our observation that dihadron correlation in AMPT with the “smoothed” initial conditions does not have any ridge structure. In addition, jet production is approximately proportional to binary nucleon-nucleon collisions and should also be correlated with hot spots. Therefore, the ridge in dihadron correlation can occur for both hard and soft trigger hadrons. In non-central collisions, the effects of jet-medium interaction and hot spots are shown to give only a broad single peak in di- hadron correlation on the away-side [25]. With a large rapidity gap, it amounts to a dihadron correlation due to momentum conservation [26].
Since hard direct photons are produced uniformly in azimuthal angle, any structure in γ-hadron azimuthal correlation with large pγtrig can only come from γ-
triggered jets and jet-induced medium excitation [10]. Shown in Fig. 7 are dihadron correlations (solid) after subtraction of harmonic flows as compared with γ-hadron correlation (dashed). The two correlations are compara- ble in magnitude but dihadron has a more pronounced double-peak which can be attributed to additional di- hadrons from hot spots and the geometric bias toward surface and tangential emission that enhances deflection of jet showers and jet-induced medium excitation [10] by radial flow. Such difference is important to measure in experiments that will provide critical information on jet- induced medium excitation and hydrodynamic evolution of hot spots in high-energy heavy-ion collisions.
This work is supported by the NSFC of China under Projects Nos. 10705044, 10975059, 11035009, the Knowl- edge Innovation Project of Chinese Academy of Sciences under Grant No. KJCX2-EW-N01 and by the U.S. DOE under Contract No. DE-AC02-05CH11231 and within the framework of the JET Collaboration.
[1] X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 (1992).
[2] S. S. Adler et al. , Phys. Rev. Lett. 91, 072301 (2003). [3] C. Adler et al., Phys. Rev. Lett. 90, 082302 (2003). [4] J. Adams et al. , Phys. Rev. Lett. 95, 152301 (2005). [5] S. S. Adler et al., Phys. Rev. Lett. 97, 052301 (2006). [6] B. I. Abelev et al., Phys. Rev. Lett. 102, 052302 (2009). [7] J. Casalderrey-Solana, E. V. Shuryak and D. Teaney,
Nucl. Phys. A 774, 577 (2006). [8] H. Stoecker, Nucl. Phys. A 750, 121 (2005). [9] B. Betz, et al., Phys. Rev. C 79, 034902 (2009). [10] H. L. Li et al., Phys. Rev. Lett. 106, 012301 (2011). [11] J. Takahashi et al., Phys. Rev. Lett. 103, 242301 (2009). [12] E. Shuryak, Phys. Rev. C80, 054908 (2009). [13] B. Alver and G. Roland, Phys. Rev. C 81, 054905 (2010).
[14] B. Schenke, S. Jeon, C. Gale, arXiv:1009.3244 [hep-ph]. [15] K. Werner et al., Phys. Rev. C82, 044904 (2010). [16] F. Gelis, T. Lappi, R. Venugopalan, Phys. Rev. D78,
054020 (2008). [17] S. Gavin, L. McLerran, G. Moschelli, Phys. Rev. C79,
051902 (2009). [18] J. Adams et al., Phys. Rev. C 73, 064907 (2006). [19] B. I. Abelev et al., Phys. Rev. C 80, 064912 (2009). [20] B. Zhang et al., Phys. Rev. C 61, 067901 (2000). [21] Z. W. Lin et al., Phys. Rev. C 72, 064901 (2005). [22] X.-N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501
(1991). [23] G. Y. Qin et al., arXiv:1009.1847 [nucl-th]. [24] M. Ploskon, Acta Phys. Hung. A 27, 255 (2006). [25] J. Xu, C. M. Ko, arXiv:1011.3750 [nucl-th].

5 [26] M. Luzum, arXiv:1011.5773 [nucl-th].