# Review of Anisotropic Flow Correlations in Ultrarelativistic Heavy-Ion Collisions

Authors: You Zhou
Journal: Advances in High Energy Physics
 (2016)

---

## Abstract

Anisotropic flow phenomena are a key probe of the existence of Quark-Gluon Plasma. Several new observables associated with correlations between anisotropic flow harmonics are developed, which are expected to be sensitive to the initial fluctuations and transport properties of the created matter in heavy-ion collisions. I review recent developments of correlations of anisotropic flow harmonics. The experimental measurements, together with the comparisons to theoretical model calculations, open up new opportunities of exploring novel QCD dynamics in heavy-ion collisions.


---

## Body

## 1. Introduction

One of the fundamental questions in the phenomenology of Quantum Chromo Dynamics (QCD) is, what are the properties of matter at extreme densities and temperatures where quarks and gluons are in a new state of matter, the so-called Quark-Gluon Plasma (QGP)? [1, 2]. Collisions of high-energy heavy ions, at the Brookhaven Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC), allow us to create and study the properties of the QGP matter in the laboratory. This matter expands under large pressure gradients, which transfer the inhomogeneous initial conditions into azimuthal anisotropy of produced particles in momentum space. This anisotropy of produced particles is one of the probes of the properties of the QGP [3, 4]. It can be characterized by an expansion of the single-particle azimuthal distribution P
(
φ
):(1)Pφ=12π∑n=-∞+∞Vn→e-inφ,where φ is the azimuthal angle of emitted particles, V
n
→ is the nth order flow vector defined as V
n
→
=
v
n
e
i
n
Ψ
n, its magnitude v
n is the nth order anisotropic flow harmonic, and its orientation is symmetry plane (participant plane) angle Ψ
n. Alternatively, this anisotropy can be generally given by the joint probability density function (PDF) in terms of v
n and Ψ
n as(2)Pvm,vn,…,Ψm,Ψn,…=1NeventdNeventvmvn⋯dvmdvn⋯dΨmdΨn⋯. In the last decade, the experimental measurements of anisotropic flow v
n [5–55], combined with theoretical advances from calculations made in a variety of frameworks [56–62], have led to broad and deep knowledge of initial conditions and properties of the created hot/dense QCD matter. In particular, the precision anisotropic flow measurements based on the huge data collected at the LHC experiments and the successful description from hydrodynamic calculations demonstrate that the QGP created in heavy-ion collisions behaves like a strongly coupled liquid with a very small specific shear viscosity η
/
s [63–68], which is close to a quantum limit 1/4π [69].It has been investigated into great details of event-by-event fluctuations of single flow harmonic. Based on the measurements of higher-order cumulants of anisotropic flow [43, 48, 51, 74, 75] and the event-by-event v
n distributions [40], it was realized that the newly proposed Elliptic-Power function [76–78] gives the best description of underlying PDF of single harmonic v
n distributions [72, 79, 80]. On the other hand, it has been known for a while that both the flow harmonic (magnitude) v
n and its symmetry plane (orientation) Ψ
n of the flow vector V
n
→ fluctuate event-by-event [81–83], but only recently p
T and η dependent flow angle (Ψ
n) and magnitude (v
n) were predicted by hydrodynamic calculations [84, 85]. Many indications were quickly obtained in experiments by looking at the deviations from unity of v
n
[
2
]
/
v
n
{
2
} [86] and factorization ratio r
n [52, 55, 86]. These measurements were nicely predicted or reproduced by hydrodynamic calculations and are found to be sensitive to the initial state density fluctuations and/or the shear viscosity of the expanding fireball medium [84, 85, 87]. Most of these above-mentioned studies are focused on the fluctuations of single flow harmonics and their corresponding symmetry planes, as a function of collisions centrality, transverse momentum p
T, and pseudorapidity η. Results of correlations between symmetry planes [28, 41] reveal a new type of correlations between different order flow vectors, which was investigated in the observable of v
2
n
/
Ψ
n before [88–90]. In particular, some of the symmetry planes correlations show quite different centrality dependence from the initial state and final state, and this characteristic sign change during system evolution is correctly reproduced by theoretical calculations [62, 82, 91], thus confirming the validity of hydrodynamic framework in heavy-ion collisions and further yielding valuable additional insights into the fluctuating initial conditions and hydrodynamic response [62, 82, 92].In addition to all these observables, the (anti)correlations between anisotropic flow harmonicsv
m and v
n are found to be extremely interesting [45, 62, 70, 71, 93]. A completely new set of information on the joint probability density function (PDF) can be obtained from the rich correlation pattern observed in experiments. On the other hand, no existent theoretical calculations [62, 70, 71, 93] could provide quantitative descriptions of data [36]. Thus, it is crucial to investigate in depth the relationship between different flow harmonics: whether they are correlated, anticorrelated, or not correlated from both experimental and theoretical points of view.


## 2. Correlations ofv
n and v
m Fluctuations

It is found recently that the relationship between different order flow harmonics can be used to probe the initial state conditions and the hydrodynamic response of the QGP [36, 71, 93–95]. In order to better understand the event-by-event P
(
φ
) distribution, it is critical to investigate the relationship between v
m and v
n. Considering the naive ellipsoidal shape of the overlap region in noncentral heavy-ion collisions generating nonvanishing even flow harmonics v
2
n, the correlations between the even flow harmonics are expected. However, it is not straightforward to use geometrical argument to explain the relationship between even flow harmonics for central collisions, where all the harmonics are driven by fluctuations instead of geometry, and to explain the relationship between even and odd flow harmonics for central and noncentral collisions [80]. A linear correlation function c
(
v
m
,
v
n
) was proposed to study the relationship between v
m and v
n [83]. It is defined as(3)cvm,vn=vm-vmevvn-vnevσvmσvnev,where σ
v
m is the standard deviation of the quantity v
m; c
(
v
m
,
v
n
) is 1 (or −1) if v
m and v
n are linearly (antilinearly) correlated and is 0 if they are not correlated. It was shown in Figure 1 that there is an anticorrelation between v
2 and v
3, while a correlation was observed between v
2 and v
4. In addition, it was demonstrated that c
(
v
2
,
v
4
) depends on both the initial conditions and η
/
s, while c
(
v
2
,
v
3
) is only sensitive to η
/
s [83]. Nevertheless, it cannot be accessible easily in experimental measurements, which rely on two-particle and multiparticle correlations techniques. Thus, it is critical to find an observable which studies the relationship between flow harmonics without contributions from symmetry plane correlations and can be accessed with observable techniques from experiments. Two different approaches, named E
v
e
n
t
S
h
a
p
e
E
n
g
i
n
e
e
r
i
n
g and S
y
m
m
e
t
r
i
c
C
u
m
u
l
a
n
t, are discussed in the following section.Figure 1
p
T dependence of c
(
v
2
,
v
3
) (a) and c
(
v
2
,
v
4
) (b) in centrality 20–30% in Pb–Pb collisions at s
N
N
=
2.76 TeV. Figures are taken from [83].
(a)
(b)


### 2.1. Event Shape Engineering (ESE)

The first experimental attempt was made by ATLAS Collaboration [45], using the Event Shape Engineering (ESE) [96]. This is a technique to select events according to the magnitude of reduced flow vector V
n
→. Figure 2 shows the performance of event shape selection on V
2 (a) and V
3 (b) in ATLAS detector. For each centrality the data sample is divided into several event classes according to V
2 or V
3 distributions. Then v
2 and v
3 relationship was investigated by measurements of v
2 and v
3 in each event class from ESE selection. Without using ESE selection, a boomerang-like pattern was observed for the centrality dependence of v
2-v
3 correlation. This is mainly due to the fact that v
3 has weaker centrality dependence than v
2. By using ESE, it was observed in Figure 3(b) that, for event class with the same centrality (shown as the same color), v
3 decreases as v
2 increases. It suggests that v
2 is anticorrelated with v
3. Considering the linear hydrodynamic response of v
2 and v
3 from eccentricity ε
2 and triangularity ε
3, the anticorrelation between v
2 and v
3 might reveal the anticorrelation between ε
2 and ε
3 of the initial geometry. This indication of initial anticorrelations between ε
2 and ε
3 is observed in model calculations [96, 97].Figure 2
Distributions ofV
2 (a) and V
3 (b) calculated with ATLAS forward calorimeter for centrality interval 0-1%. Figures are taken from [45].
(a)
(b)Figure 3
The correlation ofv
2 (x-axis) with v
3 (y-axis) measured in 0.5 <
p
T
<
2 GeV/c. (a) shows v
2 and v
3 values for fourteen 5% centrality intervals over the centrality range of 0–70% without event shape selection. (b) shows v
2 and v
3 values in 15 q
2 intervals in seven centrality ranges (markers) with larger v
2 value corresponding to larger q
2 value. Figures are taken from [45].
(a)
(b)Figure4 shows the investigation of relationship between v
2 and v
4. A boomerang-like pattern, although weaker than that for v
2-v
3 relationship shown in Figure 3(a), is observed in Figure 4(a), prior to the ESE selection. After the ESE selection, it is found in Figure 4(b) that v
4 increases with increasing v
2. This suggests a correlation between the two harmonics and it can be understood by the interplay between linear and nonlinear collective dynamics in the system evolution [45]. This nonlinear contribution of v
4 from v
2 is further investigated by fitting the correlation pattern using v
4
=
c
0
2
+
(
c
1
v
2
2
)
2, where c
0 and c
1 denote the linear and nonlinear components. It is found that the linear component has weak centrality dependence, while the nonlinear component, increasing dramatically with collision centrality, becomes the dominant contribution in the most peripheral collisions [45].Figure 4
The correlation ofv
2 (x-axis) with v
4 (y-axis) measured in 0.5 <
p
T
<
2 GeV/c. (a) shows v
2 and v
3 values for fourteen 5% centrality intervals over the centrality range of 0–70% without event shape selection. (b) shows v
2 and v
4 values in 15 q
2 intervals in seven centrality ranges (markers) with larger v
2 value corresponding to larger q
2 value. Figures are taken from [45].
(a)
(b)These (anti)correlation patterns betweenv
m and v
n observed in experiments open a new window to the understanding of the collectivity phenomena in heavy-ion collisions. However, it was also noticed that these measurements were based on 2-particle correlations, which might be suffered by nonflow effects, and they require subdividing such calculations and modeling resolutions associated with ESE due to finite event-wise multiplicities. Considering the computational constraints, this approach cannot be performed easily in hydrodynamic calculations which usually are based on limited statistics compared to experimental data.


### 2.2. Symmetric Cumulants (SC)

A new type of observables for the analyses of flow harmonic correlations,symmetric cumulants (originally namedStandard Candles (SC) in [93]), was proposed as S
C
(
m
,
n
)
=
⟪
cos
⁡
(
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
)
⟫
c. If m
≠
n, the isotropic part of the corresponding four-particle cumulant is given by (4)⟪cos⁡mφ1+nφ2-mφ3-nφ4⟫c=⟪cos⁡mφ1+nφ2-mφ3-nφ4⟫-⟪cos⁡mφ1-φ2⟫⟪cos⁡nφ1-φ2⟫=vm2vn2-vm2vn2.For a detector with uniform acceptance in azimuthal direction, the asymmetric terms, f
o
r
e
x
a
m
p
l
e, ⟪
cos
⁡
(
m
φ
1
-
n
φ
2
)
⟫, are averaged to zero. The single event 4-particle correlation ⟪
c
o
s
(
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
)
⟫ could be calculated as(5)
c
o
s
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
=
1
M
M
-
1
M
-
2
M
-
3
V
m
2
V
n
2
-
2
R
e
V
m
+
n
V
m
∗
V
n
∗
-
2
R
e
V
m
V
m
-
n
∗
V
n
∗
+
V
m
+
n
2
+
V
m
-
n
2
-
M
-
4
V
m
2
+
V
n
2
+
M
M
-
6
.And the single event 2-particle correlation ⟪
c
o
s
[
m
(
φ
1
-
φ
2
)
]
⟫ could be obtained as(6)cosmφ1-φ2=1MM-1Vm2-M.Then, the weights of M
(
M
-
1
) and M
(
M
-
1
)
(
M
-
2
)
(
M
-
3
) are used to get the event-averaged 2-particle and 4-particle correlations, as introduced in [93]. Due to the definition, this new type of 4-particle cumulants S
C
(
m
,
n
) is independent of the symmetry planes Ψ
m and Ψ
n and is expected to be less sensitive to nonflow correlations, which should be strongly suppressed in 4-particle cumulants. This was confirmed by S
C
(
m
,
n
) calculation using HIJING model [98, 99] which does not include anisotropic collectivity but, for example, azimuthal correlations due to jet production. It is observed that both ⟪
c
o
s
(
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
)
⟫ and ⟪
c
o
s
[
m
(
φ
1
-
φ
2
)
]
⟫
⟪
c
o
s
[
n
(
φ
1
-
φ
2
)
]
⟫ are nonzero, while S
C
(
m
,
n
) are compatible with zero in HIJING simulations [36]. This confirms that S
C
(
m
,
n
) measurements are nearly insensitive to nonflow correlations. Therefore, it is believed that S
C
(
m
,
n
) is nonzero if there is (anti)correlations of v
n and v
m. The investigation of S
C
(
m
,
n
) will allow us to know whether finding v
m larger than v
m in an event will enhance or reduce the probability of finding v
n larger than v
n in that event, which provides unique information for the event-by-event simulations of anisotropic flow harmonics.Figure5 shows the first calculation of S
C
(
4,2
) (solid markers) and S
C
(
3,2
) (open markers) as a function of centrality from AMPT model [93]. Nonzero values for both S
C
(
4,2
) and S
C
(
3,2
) are observed. Positive S
C
(
4,2
) suggests a correlation between the event-by-event fluctuations of v
2 and v
4, which indicates that finding v
2 larger than 〈
v
2
〉 in an event enhances the probability of finding v
4 larger than 〈
v
4
〉 in that event. On the other hand, the negative results of S
C
(
3,2
) imply that finding v
2 larger than 〈
v
2
〉 enhances the probability of finding v
3 smaller than 〈
v
3
〉 [93].Figure 5
The centrality dependence of symmetric cumulantsS
C
(
4,2
) and S
C
(
3,2
) at s
N
N. Figures are taken from [72, 93].
(a)
(b)Several configurations of the AMPT model have been investigated to better understand the results based on AMPT simulations [93]. Partonic interactions can be tweaked by changing the partonic cross section: the default value is 10 mb, while using 3 mb generates weaker partonic interactions in ZPC [100, 101]. One can also change the hadronic interactions by controlling the termination time in ART. Setting NTMAX = 3, where NTMAX is a parameter which controls the number of time steps in ART (rescattering time), will effectively turn off the hadronic interactions [100, 101]. For S
C
(
4,2
) and S
C
(
3,2
) calculations for three different scenarios, (a) 3 mb, (b) 10 mb, and (c) 10 mb, no rescattering is presented in Figure 5. It is found that when the partonic cross section is decreasing from 10 mb (lower shear viscosity) to 3 mb (higher shear viscosity), the strength of S
C
(
4,2
) decreases. Additionally, the “10 mb, no rescattering” setup seems to give slightly smaller magnitudes of S
C
(
4,2
) and S
C
(
3,2
).Further studies have been performed in AMPT initial conditions, based on the observable ofS
C
(
m
,
n
)
ε which is defined as ε
m
2
ε
n
2
-
ε
m
2
ε
n
2 [72]. The centrality dependence of S
C
(
4,2
)
ε and S
C
(
3,2
)
ε is presented as red circles and blue diamonds in Figure 5(b). Positive and increasing trend from central to peripheral collisions has been observed for S
C
(
4,2
)
ε. In contrast, negative and decreasing trend was observed for S
C
(
3,2
)
ε in the AMPT initial conditions. This shows that finding ε
2 larger than 〈
ε
2
〉 in an event enhances the probability of finding ε
4 larger than 〈
ε
4
〉, while in parallel enhancing the probability of finding ε
3 smaller than 〈
ε
3
〉 in that event. Same conclusions were obtained using MC-Glauber initial conditions [75].Based on AMPT calculations, it seems that the signs ofS
C
(
m
,
n
)
v (for m
,
n
=
2,3
,
4) in the final state are determined by the correlations of S
C
(
m
,
n
)
ε in the initial state, while its magnitude also depends on the properties of the created system. This clearly suggests that S
C
(
m
,
n
)
v is a new promising observable to constrain the initial conditions and the transport properties of the system.The first experimental measurements of centrality dependence ofS
C
(
4,2
) (red squares) and S
C
(
3,2
) (blue circles) are presented in Figure 6(a). Positive values of S
C
(
4,2
) are observed for all cases of centrality. This confirms a correlation between the event-by-event fluctuations of v
2 and v
4. On the other hand, the measured negative results of S
C
(
3,2
) show the anticorrelation between v
2 and v
3 magnitudes. The same measurements are performed using the like-sign technique, which is another powerful approach to estimate nonflow effects [27]. It was found that the difference between correlations for like-sign and all charged combinations, which might be mainly due to nonflow effects, is much smaller compared to the magnitudes of S
C
(
m
,
n
) itself. This further proves that nonzero values of S
C
(
m
,
n
) measured in experiments cannot be explained by nonflow effects solely.Figure 6
The centrality dependence of symmetric cumulantsS
C
(
4,2
) (red markers) and S
C
(
3,2
) (blue markers) at s
N
N
=
2.76 TeV Pb–Pb collisions. Figures are taken from [36].
(a)
(b)In addition, the comparison between experimental data and the event-by-event perturbative-QCD+saturation+hydro (“EKRT”) calculations [62], which incorporate both initial conditions and hydrodynamic evolution, is shown in Figure 6. It was shown that this model can capture quantitatively the centrality dependence of individual v
2, v
3, and v
4 harmonics in central and mid-central collisions [62]. However, it can only qualitatively but not quantitatively predict S
C
(
m
,
n
) measurements by ALICE. For given η
/
s
(
T
) parameterization tuned by individual flow harmonic, the calculation cannot describe S
C
(
4,2
) and S
C
(
3,2
) simultaneously for any single centrality. Experimental measurements are also compared to the VISH2+1 model calculations (see Figure 7), using various combinations of initial conditions (IC) from (a) MC-Glb, (b) MC-KLN, and (c) MC-AMPT with η
/
s
= 0.08 and 0.20. It is noticed that the one with MC-Glb IC and η
/
s
= 0.08 is compatible with S
C
(
4,2
) measurement and the calculation with MC-AMPT IC and η
/
s
= 0.08 can describe S
C
(
3,2
) measurement [70]. However, just like EKRT calculations, none of these combinations is able to describe S
C
(
4,2
) and S
C
(
3,2
) simultaneously. Thus, it is concluded that the new S
C
(
m
,
n
) observables provide better handle on the initial conditions and η
/
s
(
T
) than each of the individual harmonic measurements alone.Figure 7
The centrality dependence of symmetric cumulantsS
C
(
4,2
) (red markers) and S
C
(
3,2
) (blue markers) at s
N
N
=
2.76 TeV Pb–Pb collisions by VISH2+1 simulations. Figures are taken from [70].After being presented for the first time at Quark Matter 2015 conference, preliminary results ofS
C
(
4,2
) and S
C
(
3,2
) gained a lot of attention [102]. One of the key suggestions was to normalize S
C
(
m
,
n
) by dividing with the products v
m
2
v
n
2 in order to get rid of influences from individual flow harmonics. The results are shown in Figure 6(b), with normalized S
C
(
3,2
) and S
C
(
4,2
) observables by dividing with the products v
3
2
v
2
2 and v
4
2
v
2
2, respectively [36]. The 2-particle correlations v
m
2 and v
n
2 are obtained with a pseudorapidity gap of Δ
η
>
1.0 to suppress contributions from nonflow effects. It was shown in Figure 8(a) that the normalized S
C
(
4,2
) observable exhibits clear sensitivity to different η
/
s parameterizations and the initial conditions, which provides a unique opportunity to discriminate between various possibilities of the detailed setting of η
/
s
(
T
) of the produced QGP and the initial conditions used in hydrodynamic calculations. On the other hand, normalized S
C
(
3,2
) is independent of the setting of η
/
s
(
T
). In addition, it was demonstrated in Figure 9 that the normalized S
C
(
3,2
), also named N
S
C
v
(
3,2
) in the following text, is compatible with its corresponding observable S
C
ε
(
3,2
) in the initial state. Thus, N
S
C
v
(
3,2
) could be taken as golden observable to directly constrain initial conditions without demands for precise knowledge of transport properties of the system [70]. Furthermore, none of existing theoretical calculations can reproduce the data; there is still a long way to go for the development of hydrodynamic calculations.Figure 8
The centrality dependence of normalized symmetric cumulantsN
S
C
(
m
,
n
) at s
N
N
=
2.76 TeV Pb–Pb collisions by VISH2+1 simulations. Figures are taken from [70].Figure 9
The centrality dependence ofN
S
C
(
3,2
) (a, b, and c) and C
(
v
3
2
,
v
2
2
) (d, e, and f) and the corresponding observables in the initial conditions at s
N
N
=
2.76 TeV Pb–Pb collisions from VISH2+1. Figures are taken from [70].
(a)
(b)
(c)
(d)
(e)
(f)Predictions of relationship between other harmonics are provided in [70] and shown in Figure 8. Besides different sensitivities to IC and η
/
s as seen above, the centrality dependence of the relationship between flow harmonics seems quite different. For instance, despite the differences in the initial conditions, a maximum value of S
C
(
5,3
) is observed in central collision using η
/
s
=
0.20, while the maximum value is seen in more peripheral collision if η
/
s
=
0.08 is used.Compared to the previous measurements of relationship between flow harmonics investigated using the ESE technique,S
C
(
m
,
n
) observable provides a quantitative measure of these correlation strengths. Further investigations on relationship between flow harmonics using list of observables in Table 1 could be performed as a function of centrality, transverse momentum, and pseudorapidity et al., which is clearly nontrivial. Although one did not use the information of symmetry planes in both E
S
E and S
C studies, recent study just reveals that flow harmonic correlations might not be completely independent on symmetry plane correlations [73]. The proportionality relations between symmetric cumulants involving higher harmonics v
4 or v
5 and symmetry plane correlations are derived, which seem to build the bridge between flow harmonic correlations and flow angle correlations (symmetry plane correlations). This might point out to a new direction of investigations of correlations between flow vectors and will shed a new light on the nature of fluctuating initial conditions and η
/
s of the created QGP in heavy-ion collisions.Table 1
List of observables for correlations of flow harmonics, including all combinations of symmetric 2-harmonic 4-particle cumulants (up tov
6).
Observables
Equations
Number of particles
Exp.
Th.
cos
⁡
2
φ
1
+
3
φ
2
-
2
φ
3
-
3
φ
4
c
〈
v
2
2
v
3
2
〉
-
v
2
2
v
3
2
4
[36]
[70–72]
cos
⁡
2
φ
1
+
4
φ
2
-
2
φ
3
-
4
φ
4
c
〈
v
2
2
v
4
2
〉
-
v
2
2
v
4
2
4
[36]
[70–73]
cos
⁡
2
φ
1
+
5
φ
2
-
2
φ
3
-
5
φ
4
c
〈
v
2
2
v
5
2
〉
-
v
2
2
v
5
2
4
[70, 71, 73]
cos
⁡
2
φ
1
+
6
φ
2
-
2
φ
3
-
6
φ
4
c
〈
v
2
2
v
6
2
〉
-
v
2
2
v
6
2
4
cos
⁡
3
φ
1
+
4
φ
2
-
3
φ
3
-
4
φ
4
c
〈
v
3
2
v
4
2
〉
-
v
3
2
v
4
2
4
[70]
cos
⁡
3
φ
1
+
5
φ
2
-
3
φ
3
-
5
φ
4
c
〈
v
3
2
v
5
2
〉
-
v
3
2
v
5
2
4
[70, 71, 73]
cos
⁡
3
φ
1
+
6
φ
2
-
3
φ
3
-
6
φ
4
c
〈
v
3
2
v
6
2
〉
-
v
3
2
v
6
2
4
cos
⁡
4
φ
1
+
5
φ
2
-
4
φ
3
-
5
φ
4
c
〈
v
4
2
v
5
2
〉
-
v
4
2
v
5
2
4
cos
⁡
4
φ
1
+
6
φ
2
-
4
φ
3
-
6
φ
4
c
〈
v
4
2
v
6
2
〉
-
v
4
2
v
6
2
4
cos
⁡
5
φ
1
+
6
φ
2
-
5
φ
3
-
6
φ
4
c
〈
v
5
2
v
6
2
〉
-
v
5
2
v
6
2
4
⋯
⋯
6


## 2.1. Event Shape Engineering (ESE)

The first experimental attempt was made by ATLAS Collaboration [45], using the Event Shape Engineering (ESE) [96]. This is a technique to select events according to the magnitude of reduced flow vector V
n
→. Figure 2 shows the performance of event shape selection on V
2 (a) and V
3 (b) in ATLAS detector. For each centrality the data sample is divided into several event classes according to V
2 or V
3 distributions. Then v
2 and v
3 relationship was investigated by measurements of v
2 and v
3 in each event class from ESE selection. Without using ESE selection, a boomerang-like pattern was observed for the centrality dependence of v
2-v
3 correlation. This is mainly due to the fact that v
3 has weaker centrality dependence than v
2. By using ESE, it was observed in Figure 3(b) that, for event class with the same centrality (shown as the same color), v
3 decreases as v
2 increases. It suggests that v
2 is anticorrelated with v
3. Considering the linear hydrodynamic response of v
2 and v
3 from eccentricity ε
2 and triangularity ε
3, the anticorrelation between v
2 and v
3 might reveal the anticorrelation between ε
2 and ε
3 of the initial geometry. This indication of initial anticorrelations between ε
2 and ε
3 is observed in model calculations [96, 97].Figure 2
Distributions ofV
2 (a) and V
3 (b) calculated with ATLAS forward calorimeter for centrality interval 0-1%. Figures are taken from [45].
(a)
(b)Figure 3
The correlation ofv
2 (x-axis) with v
3 (y-axis) measured in 0.5 <
p
T
<
2 GeV/c. (a) shows v
2 and v
3 values for fourteen 5% centrality intervals over the centrality range of 0–70% without event shape selection. (b) shows v
2 and v
3 values in 15 q
2 intervals in seven centrality ranges (markers) with larger v
2 value corresponding to larger q
2 value. Figures are taken from [45].
(a)
(b)Figure4 shows the investigation of relationship between v
2 and v
4. A boomerang-like pattern, although weaker than that for v
2-v
3 relationship shown in Figure 3(a), is observed in Figure 4(a), prior to the ESE selection. After the ESE selection, it is found in Figure 4(b) that v
4 increases with increasing v
2. This suggests a correlation between the two harmonics and it can be understood by the interplay between linear and nonlinear collective dynamics in the system evolution [45]. This nonlinear contribution of v
4 from v
2 is further investigated by fitting the correlation pattern using v
4
=
c
0
2
+
(
c
1
v
2
2
)
2, where c
0 and c
1 denote the linear and nonlinear components. It is found that the linear component has weak centrality dependence, while the nonlinear component, increasing dramatically with collision centrality, becomes the dominant contribution in the most peripheral collisions [45].Figure 4
The correlation ofv
2 (x-axis) with v
4 (y-axis) measured in 0.5 <
p
T
<
2 GeV/c. (a) shows v
2 and v
3 values for fourteen 5% centrality intervals over the centrality range of 0–70% without event shape selection. (b) shows v
2 and v
4 values in 15 q
2 intervals in seven centrality ranges (markers) with larger v
2 value corresponding to larger q
2 value. Figures are taken from [45].
(a)
(b)These (anti)correlation patterns betweenv
m and v
n observed in experiments open a new window to the understanding of the collectivity phenomena in heavy-ion collisions. However, it was also noticed that these measurements were based on 2-particle correlations, which might be suffered by nonflow effects, and they require subdividing such calculations and modeling resolutions associated with ESE due to finite event-wise multiplicities. Considering the computational constraints, this approach cannot be performed easily in hydrodynamic calculations which usually are based on limited statistics compared to experimental data.


## 2.2. Symmetric Cumulants (SC)

A new type of observables for the analyses of flow harmonic correlations,symmetric cumulants (originally namedStandard Candles (SC) in [93]), was proposed as S
C
(
m
,
n
)
=
⟪
cos
⁡
(
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
)
⟫
c. If m
≠
n, the isotropic part of the corresponding four-particle cumulant is given by (4)⟪cos⁡mφ1+nφ2-mφ3-nφ4⟫c=⟪cos⁡mφ1+nφ2-mφ3-nφ4⟫-⟪cos⁡mφ1-φ2⟫⟪cos⁡nφ1-φ2⟫=vm2vn2-vm2vn2.For a detector with uniform acceptance in azimuthal direction, the asymmetric terms, f
o
r
e
x
a
m
p
l
e, ⟪
cos
⁡
(
m
φ
1
-
n
φ
2
)
⟫, are averaged to zero. The single event 4-particle correlation ⟪
c
o
s
(
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
)
⟫ could be calculated as(5)
c
o
s
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
=
1
M
M
-
1
M
-
2
M
-
3
V
m
2
V
n
2
-
2
R
e
V
m
+
n
V
m
∗
V
n
∗
-
2
R
e
V
m
V
m
-
n
∗
V
n
∗
+
V
m
+
n
2
+
V
m
-
n
2
-
M
-
4
V
m
2
+
V
n
2
+
M
M
-
6
.And the single event 2-particle correlation ⟪
c
o
s
[
m
(
φ
1
-
φ
2
)
]
⟫ could be obtained as(6)cosmφ1-φ2=1MM-1Vm2-M.Then, the weights of M
(
M
-
1
) and M
(
M
-
1
)
(
M
-
2
)
(
M
-
3
) are used to get the event-averaged 2-particle and 4-particle correlations, as introduced in [93]. Due to the definition, this new type of 4-particle cumulants S
C
(
m
,
n
) is independent of the symmetry planes Ψ
m and Ψ
n and is expected to be less sensitive to nonflow correlations, which should be strongly suppressed in 4-particle cumulants. This was confirmed by S
C
(
m
,
n
) calculation using HIJING model [98, 99] which does not include anisotropic collectivity but, for example, azimuthal correlations due to jet production. It is observed that both ⟪
c
o
s
(
m
φ
1
+
n
φ
2
-
m
φ
3
-
n
φ
4
)
⟫ and ⟪
c
o
s
[
m
(
φ
1
-
φ
2
)
]
⟫
⟪
c
o
s
[
n
(
φ
1
-
φ
2
)
]
⟫ are nonzero, while S
C
(
m
,
n
) are compatible with zero in HIJING simulations [36]. This confirms that S
C
(
m
,
n
) measurements are nearly insensitive to nonflow correlations. Therefore, it is believed that S
C
(
m
,
n
) is nonzero if there is (anti)correlations of v
n and v
m. The investigation of S
C
(
m
,
n
) will allow us to know whether finding v
m larger than v
m in an event will enhance or reduce the probability of finding v
n larger than v
n in that event, which provides unique information for the event-by-event simulations of anisotropic flow harmonics.Figure5 shows the first calculation of S
C
(
4,2
) (solid markers) and S
C
(
3,2
) (open markers) as a function of centrality from AMPT model [93]. Nonzero values for both S
C
(
4,2
) and S
C
(
3,2
) are observed. Positive S
C
(
4,2
) suggests a correlation between the event-by-event fluctuations of v
2 and v
4, which indicates that finding v
2 larger than 〈
v
2
〉 in an event enhances the probability of finding v
4 larger than 〈
v
4
〉 in that event. On the other hand, the negative results of S
C
(
3,2
) imply that finding v
2 larger than 〈
v
2
〉 enhances the probability of finding v
3 smaller than 〈
v
3
〉 [93].Figure 5
The centrality dependence of symmetric cumulantsS
C
(
4,2
) and S
C
(
3,2
) at s
N
N. Figures are taken from [72, 93].
(a)
(b)Several configurations of the AMPT model have been investigated to better understand the results based on AMPT simulations [93]. Partonic interactions can be tweaked by changing the partonic cross section: the default value is 10 mb, while using 3 mb generates weaker partonic interactions in ZPC [100, 101]. One can also change the hadronic interactions by controlling the termination time in ART. Setting NTMAX = 3, where NTMAX is a parameter which controls the number of time steps in ART (rescattering time), will effectively turn off the hadronic interactions [100, 101]. For S
C
(
4,2
) and S
C
(
3,2
) calculations for three different scenarios, (a) 3 mb, (b) 10 mb, and (c) 10 mb, no rescattering is presented in Figure 5. It is found that when the partonic cross section is decreasing from 10 mb (lower shear viscosity) to 3 mb (higher shear viscosity), the strength of S
C
(
4,2
) decreases. Additionally, the “10 mb, no rescattering” setup seems to give slightly smaller magnitudes of S
C
(
4,2
) and S
C
(
3,2
).Further studies have been performed in AMPT initial conditions, based on the observable ofS
C
(
m
,
n
)
ε which is defined as ε
m
2
ε
n
2
-
ε
m
2
ε
n
2 [72]. The centrality dependence of S
C
(
4,2
)
ε and S
C
(
3,2
)
ε is presented as red circles and blue diamonds in Figure 5(b). Positive and increasing trend from central to peripheral collisions has been observed for S
C
(
4,2
)
ε. In contrast, negative and decreasing trend was observed for S
C
(
3,2
)
ε in the AMPT initial conditions. This shows that finding ε
2 larger than 〈
ε
2
〉 in an event enhances the probability of finding ε
4 larger than 〈
ε
4
〉, while in parallel enhancing the probability of finding ε
3 smaller than 〈
ε
3
〉 in that event. Same conclusions were obtained using MC-Glauber initial conditions [75].Based on AMPT calculations, it seems that the signs ofS
C
(
m
,
n
)
v (for m
,
n
=
2,3
,
4) in the final state are determined by the correlations of S
C
(
m
,
n
)
ε in the initial state, while its magnitude also depends on the properties of the created system. This clearly suggests that S
C
(
m
,
n
)
v is a new promising observable to constrain the initial conditions and the transport properties of the system.The first experimental measurements of centrality dependence ofS
C
(
4,2
) (red squares) and S
C
(
3,2
) (blue circles) are presented in Figure 6(a). Positive values of S
C
(
4,2
) are observed for all cases of centrality. This confirms a correlation between the event-by-event fluctuations of v
2 and v
4. On the other hand, the measured negative results of S
C
(
3,2
) show the anticorrelation between v
2 and v
3 magnitudes. The same measurements are performed using the like-sign technique, which is another powerful approach to estimate nonflow effects [27]. It was found that the difference between correlations for like-sign and all charged combinations, which might be mainly due to nonflow effects, is much smaller compared to the magnitudes of S
C
(
m
,
n
) itself. This further proves that nonzero values of S
C
(
m
,
n
) measured in experiments cannot be explained by nonflow effects solely.Figure 6
The centrality dependence of symmetric cumulantsS
C
(
4,2
) (red markers) and S
C
(
3,2
) (blue markers) at s
N
N
=
2.76 TeV Pb–Pb collisions. Figures are taken from [36].
(a)
(b)In addition, the comparison between experimental data and the event-by-event perturbative-QCD+saturation+hydro (“EKRT”) calculations [62], which incorporate both initial conditions and hydrodynamic evolution, is shown in Figure 6. It was shown that this model can capture quantitatively the centrality dependence of individual v
2, v
3, and v
4 harmonics in central and mid-central collisions [62]. However, it can only qualitatively but not quantitatively predict S
C
(
m
,
n
) measurements by ALICE. For given η
/
s
(
T
) parameterization tuned by individual flow harmonic, the calculation cannot describe S
C
(
4,2
) and S
C
(
3,2
) simultaneously for any single centrality. Experimental measurements are also compared to the VISH2+1 model calculations (see Figure 7), using various combinations of initial conditions (IC) from (a) MC-Glb, (b) MC-KLN, and (c) MC-AMPT with η
/
s
= 0.08 and 0.20. It is noticed that the one with MC-Glb IC and η
/
s
= 0.08 is compatible with S
C
(
4,2
) measurement and the calculation with MC-AMPT IC and η
/
s
= 0.08 can describe S
C
(
3,2
) measurement [70]. However, just like EKRT calculations, none of these combinations is able to describe S
C
(
4,2
) and S
C
(
3,2
) simultaneously. Thus, it is concluded that the new S
C
(
m
,
n
) observables provide better handle on the initial conditions and η
/
s
(
T
) than each of the individual harmonic measurements alone.Figure 7
The centrality dependence of symmetric cumulantsS
C
(
4,2
) (red markers) and S
C
(
3,2
) (blue markers) at s
N
N
=
2.76 TeV Pb–Pb collisions by VISH2+1 simulations. Figures are taken from [70].After being presented for the first time at Quark Matter 2015 conference, preliminary results ofS
C
(
4,2
) and S
C
(
3,2
) gained a lot of attention [102]. One of the key suggestions was to normalize S
C
(
m
,
n
) by dividing with the products v
m
2
v
n
2 in order to get rid of influences from individual flow harmonics. The results are shown in Figure 6(b), with normalized S
C
(
3,2
) and S
C
(
4,2
) observables by dividing with the products v
3
2
v
2
2 and v
4
2
v
2
2, respectively [36]. The 2-particle correlations v
m
2 and v
n
2 are obtained with a pseudorapidity gap of Δ
η
>
1.0 to suppress contributions from nonflow effects. It was shown in Figure 8(a) that the normalized S
C
(
4,2
) observable exhibits clear sensitivity to different η
/
s parameterizations and the initial conditions, which provides a unique opportunity to discriminate between various possibilities of the detailed setting of η
/
s
(
T
) of the produced QGP and the initial conditions used in hydrodynamic calculations. On the other hand, normalized S
C
(
3,2
) is independent of the setting of η
/
s
(
T
). In addition, it was demonstrated in Figure 9 that the normalized S
C
(
3,2
), also named N
S
C
v
(
3,2
) in the following text, is compatible with its corresponding observable S
C
ε
(
3,2
) in the initial state. Thus, N
S
C
v
(
3,2
) could be taken as golden observable to directly constrain initial conditions without demands for precise knowledge of transport properties of the system [70]. Furthermore, none of existing theoretical calculations can reproduce the data; there is still a long way to go for the development of hydrodynamic calculations.Figure 8
The centrality dependence of normalized symmetric cumulantsN
S
C
(
m
,
n
) at s
N
N
=
2.76 TeV Pb–Pb collisions by VISH2+1 simulations. Figures are taken from [70].Figure 9
The centrality dependence ofN
S
C
(
3,2
) (a, b, and c) and C
(
v
3
2
,
v
2
2
) (d, e, and f) and the corresponding observables in the initial conditions at s
N
N
=
2.76 TeV Pb–Pb collisions from VISH2+1. Figures are taken from [70].
(a)
(b)
(c)
(d)
(e)
(f)Predictions of relationship between other harmonics are provided in [70] and shown in Figure 8. Besides different sensitivities to IC and η
/
s as seen above, the centrality dependence of the relationship between flow harmonics seems quite different. For instance, despite the differences in the initial conditions, a maximum value of S
C
(
5,3
) is observed in central collision using η
/
s
=
0.20, while the maximum value is seen in more peripheral collision if η
/
s
=
0.08 is used.Compared to the previous measurements of relationship between flow harmonics investigated using the ESE technique,S
C
(
m
,
n
) observable provides a quantitative measure of these correlation strengths. Further investigations on relationship between flow harmonics using list of observables in Table 1 could be performed as a function of centrality, transverse momentum, and pseudorapidity et al., which is clearly nontrivial. Although one did not use the information of symmetry planes in both E
S
E and S
C studies, recent study just reveals that flow harmonic correlations might not be completely independent on symmetry plane correlations [73]. The proportionality relations between symmetric cumulants involving higher harmonics v
4 or v
5 and symmetry plane correlations are derived, which seem to build the bridge between flow harmonic correlations and flow angle correlations (symmetry plane correlations). This might point out to a new direction of investigations of correlations between flow vectors and will shed a new light on the nature of fluctuating initial conditions and η
/
s of the created QGP in heavy-ion collisions.Table 1
List of observables for correlations of flow harmonics, including all combinations of symmetric 2-harmonic 4-particle cumulants (up tov
6).
Observables
Equations
Number of particles
Exp.
Th.
cos
⁡
2
φ
1
+
3
φ
2
-
2
φ
3
-
3
φ
4
c
〈
v
2
2
v
3
2
〉
-
v
2
2
v
3
2
4
[36]
[70–72]
cos
⁡
2
φ
1
+
4
φ
2
-
2
φ
3
-
4
φ
4
c
〈
v
2
2
v
4
2
〉
-
v
2
2
v
4
2
4
[36]
[70–73]
cos
⁡
2
φ
1
+
5
φ
2
-
2
φ
3
-
5
φ
4
c
〈
v
2
2
v
5
2
〉
-
v
2
2
v
5
2
4
[70, 71, 73]
cos
⁡
2
φ
1
+
6
φ
2
-
2
φ
3
-
6
φ
4
c
〈
v
2
2
v
6
2
〉
-
v
2
2
v
6
2
4
cos
⁡
3
φ
1
+
4
φ
2
-
3
φ
3
-
4
φ
4
c
〈
v
3
2
v
4
2
〉
-
v
3
2
v
4
2
4
[70]
cos
⁡
3
φ
1
+
5
φ
2
-
3
φ
3
-
5
φ
4
c
〈
v
3
2
v
5
2
〉
-
v
3
2
v
5
2
4
[70, 71, 73]
cos
⁡
3
φ
1
+
6
φ
2
-
3
φ
3
-
6
φ
4
c
〈
v
3
2
v
6
2
〉
-
v
3
2
v
6
2
4
cos
⁡
4
φ
1
+
5
φ
2
-
4
φ
3
-
5
φ
4
c
〈
v
4
2
v
5
2
〉
-
v
4
2
v
5
2
4
cos
⁡
4
φ
1
+
6
φ
2
-
4
φ
3
-
6
φ
4
c
〈
v
4
2
v
6
2
〉
-
v
4
2
v
6
2
4
cos
⁡
5
φ
1
+
6
φ
2
-
5
φ
3
-
6
φ
4
c
〈
v
5
2
v
6
2
〉
-
v
5
2
v
6
2
4
⋯
⋯
6


## 3. Summary

In the past two decades, the underlying PDF of each single harmonicP
(
v
n
) was investigated in great detail. However, at the moment, how the joint underlying PDF, including different order symmetry planes and harmonics, is described is an open question, especially if these correlations between different flow harmonics modify the single harmonics P
(
v
n
). New observables discussed here begin to answer these open questions. Nevertheless, many more investigations between different flow harmonics, including higher-order cumulants and higher harmonics, are necessary to reasonably constrain the joint PDF and ultimately lead to new insights into the nature of fluctuation of the created matter in heavy-ion collisions. How to turn the multitude of measured and possibly measurable in future relationships between anisotropic flow harmonics into a focused search for correct initial conditions and detailed setting of η
/
s is an exciting challenge for the theory community.



---

*Source: 9365637-2016-09-20.xml*