Review of anisotropic flow correlations in ultrarelativistic heavy-ion collisions
Anisotropic flow phenomena is a key probe of the existence of Quark-Gluon Plasma. Several new observable 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.
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 :
where is the azimuthal angle of emitted particles, is the -th order flow-vector defined as , its magnitude is the -th order anisotropic flow harmonic and its orientation is symmetry plane (participant plane) angle . Alternatively, this anisotropy can be generally given by the joint probability density function () in terms of and as:
In the last decade, the experimental measurements of anisotropic flow (5); (6); (7); (8); (9); (10); (12); (11); (13); (14); (15); (16); (17); (18); (19); (20); (21); (22); (23); (24); (25); (26); (27); (28); (29); (30); (31); (32); (33); (34); (35); (36); (37); (38); (39); (40); (41); (45); (46); (47); (48); (49); (50); (51); (52); (53); (54); (55); (56); (57); (58), combined with theoretical advances from calculations made in a variety of frameworks (59); (60); (61); (62); (63); (64), have led to a 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 (66); (67); (68); (69); (70); (71), which is close to a quantum limit 1/4 (72).
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 (74); (46); (51); (54); (75) and the event-by-event distributions (40), it was realized that the newly proposed Elliptic-Power function (76); (77); (78) gives the best description of underlying of single harmonic distributions (79); (80); (73). On the other hand, it has been known for a while that both the flow harmonic (magnitude) and its symmetry plane (orientation) of the flow-vector fluctuate event-by-event (85); (86); (81), but only recently the and dependent flow angle () and magnitude () were predicted by hydrodynamic calculations (87); (88). Many indications were quickly obtained in experiments by looking at the deviations from unity of (89) and factorization ratio (89); (55); (58). These measurements were nicely predicted or reproduced by hydrodynamic calculations, and are found to be sensitive to either the initial-state density fluctuations and/or the shear viscosity of the expanding fireball medium (87); (88); (90). Most of these above mentioned studies are focused on the fluctuations of single flow harmonics and its corresponding symmetry planes, as a function of collisions centrality, transverse momentum and pseudorapidity . Results of correlations between symmetry planes (41); (28) reveal a new type of correlations between different order flow-vectors, which was investigated in the observable of before (43); (44); (42). In particular, some of the symmetry planes correlations show quite different centrality dependence from the initial- and final-state, and this characteristic sign change during system evolution is correctly reproduced by theoretical calculations (86); (65); (64), thus confirms the validity of hydrodynamic framework in heavy-ion collisions and further yields valuable additional insights into the fluctuating initial conditions and hydrodynamic response (86); (83); (64).
In addition to all these observables, the (anti-)correlations between anisotropic flow harmonics and are found to be extremely interesting (91); (92); (64); (48); (93). A completely new set of information on the joint is carried by the rich pattern observed in experiments. On the other hand, no existed theoretical calculations (91); (92); (64); (93) could provide quantitative descriptions of data (36). Thus, it’s crucial to investigate in depth of the relationship between different flow harmonics: whether they are correlated, anti-correlated or not correlated, from both experimental and theoretical point of view.
Ii Correlations of and 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 (91); (82); (84); (92); (36). In order to better understand the event-by-event distribution, it’s critical to investigate the relationship between and . Considering the naive ellipsoidal shape of the overlap region in non-central heavy ion collisions generates non-vanishing even flow harmonics , 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 odd flow harmonics for central and non-central collisions (73). A linear correlation function was proposed to study the relationship between and (81). It is defined as:
where is the standard deviation of the quantity , is 1 (or -1) if and are linearly (anti-linearly) correlated, and is 0 if not correlated.
It was shown in Fig 1 that there is an anti-correlationsbetween and , while a correlation was observed between and . In addition, it was demonstrated that depends on both the initial conditions and while is only sensitive to (81). Nevertheless, it cannot be accessible easily in experimental measurements, which rely on two- and multi-particle 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 and , are discussed in the following section.
ii.1 Event Shape Engineering (ESE)
The first experimental attempt was made by ATLAS Collaboration (48), using the Event-Shape Engineering (ESE) (94). This is a technique to select events according to the magnitude of reduced flow vector . Fig. 2 shows the performance of event shape selection on (left) and (right) in ATLAS detector. For each centrality the data sample is divided into several event classes according to the or distributions. Then the and relationship was investigated by measurements of and in each event class from ESE selection. Without using ESE selection, a boomerang-like patten was observed for the centrality dependence of - correlation. This is mainly due to the fact that has a weaker centrality dependence than . By using ESE, it was observed in Fig. 3 (right) that for event class with the same centrality (shown as the same color), decreases as increasing. It suggests that is anti-correlated with . Considering the linear hydrodynamic response of and from eccentricity and triangularity , the anti-correlation between and might reveal the anti-correlation between and of the initial geometry. This indication of initial anti-correlations between and is observed in model calculations (94); (105).
Figure 4 shows the investigation of relationship between and . A boomerang-like patten, although weaker than that for the - relationship shown in Fig. 3 (left), is observed in Fig. 4 (left), prior to the ESE selection. After the ESE selection, it is found in Fig. 4 (right) that increases with increasing . 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 (48). This non-linear contribution of from is further investigated by fitting the correlation pattern using , where and denote the linear and non-linear components. It is found that the linear component has a weak centrality dependence, while the non-linear component, increasing dramatically with collision centrality, becomes the dominant contribution in the most peripheral collisions (48).
These (anti)correlation patten between and 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 non-flow effects, and they require sub-dividing such calculations and modeling resolutions associated with ESE due to finite event-wise multiplicities. Considering the computational constraints, this approach can not be performed easily in hydrodynamic calculations which usually are based on limited statistics compared to experimental data.
ii.2 Symmetric Cumulants ()
A new type of observable for the analyses of flow harmonic correlations, Symmetric Cumulants (originally named Standard Candles (SC) in (91)), was proposed as . If , the isotropic part of the corresponding four-particle cumulant is given by:
For a detector with uniform acceptance in azimuthal direction, the asymmetric terms, e.g. , are averaged to zero. The single event 4-particle correlation could be calculated as:
And the single event 2-particle correlation could be obtained as:
Then, the weights of and are used to get the event-averaged 2- and 4-particle correlations, as introduced in (91). Due to the definition, this new type of 4-particle cumulant is independent of the symmetry planes and , and is expected to be less sensitive to non-flow correlations, which should be strongly suppressed in 4-particle cumulants. This was confirmed by the calculation using HIJING model (106); (107) which does not include anisotropic collectivity but e.g. azimuthal correlations due to jet production. It is observed that both and are non-zero, while are compatible with zero in HIJING simulations (36). This confirms that the measurements are nearly insensitive to non-flow correlations. Therefore, it is believed that is nonzero if there is (anti-)correlations of and . The investigation of will allow us to know whether finding larger than in an event will enhance or reduce the probability of finding larger than in that event, which provides a unique information for the event-by-event simulations of anisotropic flow harmonics.
Figure 5 shows the first calculation of (solid markers) and (open markers) as a function of centrality from AMPT model (91). Non-zero values for both and are observed. The positive suggests a correlation between the event-by-event fluctuations of and , which indicates that finding larger than in an event enhances the probability of finding larger than in that event. On the other hand, the negative results of implies that finding larger than enhances the probability of finding smaller than (91).
Several configurations of the AMPT model have been investigated to better understand the results based on AMPT simulations (91). 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 (109); (96). 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 (109); (96). The and calculations for three different scenarios: (a) 3 mb; (b) 10 mb; (c)10 mb, no rescattering are presented in Fig. 5 (left). 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 decreases. Additionally, the ‘10mb, no rescattering’ setup seems to give slightly smaller magnitudes of and .
Further studies have been performed in AMPT initial conditions, based on the observable of which is defined as (80). The centrality dependence of and are presented as red circles and blue diamonds in Fig. 5 (left bottom). Positive and increasing trend from central to peripheral collisions has been observed for . In contrast, negative and decreasing trend was observed for in the AMPT initial conditions. This shows that finding larger than in an event enhances the probability of finding larger than , while in parallel enhancing the probability of finding smaller than in that event. Same conclusions were obtained using MC-Glauber initial conditions (75).
Based on AMPT calculations, it seems that the signs of (for 2, 3, 4) in the final state are determined by the correlations of in the initial state, while its magnitude also depends on the properties of the created system. This clearly suggests that is a new promising observable to constrain the initial conditions and the transport properties of the system.
The first experimental measurements of centrality dependence of (red squares) and (blue circles) are presented in Fig. 6 (left). Positive values of are observed for all centralities. This confirms a correlation between the event-by-event fluctuations of and . On the other hand, the measured negative results of show the anti-correlation between and magnitudes. The same measurements are performed using the like-sign technique, which is another powerful approach to estimate non-flow effects (27). It was found that the difference between correlations for like-sign and all charged combinations, which might be mainly due to non-flow effects, are much smaller compared to the magnitudes of itself. This further proves that non-zero values of measured in experiments cannot be explained by non-flow effects solely.
In addition, the comparison between experimental data and the event-by-event perturbative-QCDsaturationhydro (“EKRT”) calculations (64), which incorporate both initial conditions and hydrodynamic evolution, is shown in Fig. 6. It was shown that this model can capture quantitatively the centrality dependence of individual , and harmonics in central and mid-central collisions (64). However, it can only qualitatively, but not quantitatively predict the measurements by ALICE. For a given parameterization tuned by individual flow harmonic, the calculation can not describe and simultaneously for any single centrality. Experimental measurements are also compared to the VISH21 model calculations (see Fig. 7), using various combinations of initial conditions (IC) from (a) MC-Glb; (b) MC-KLN and (c) MC-AMPT with 0.08 and 0.20. It is noticed that the one with MC-Glb IC and 0.08 is compatible with measurement and the calculation with MC-AMPT IC and 0.08 can describe the measurement (93). However, just like EKRT calculations, none of these combinations is able to describe and simultaneously. Thus, it is concluded that the new SC observables provide a better handle on the initial conditions and than each of the individual harmonic measurement alone.
After being presented for the first time at Quark Matter 2015 conference, preliminary results of and gained a lot of attention (108). One of the key suggestions was to normalize by dividing with the products , in order to get rid of influences from individual flow harmonics. The results are shown in Fig. 6 (right), with normalized and observables by dividing with the products and , respectively (36). The 2-particle correlations and are obtained with a pseudorapidity gap of to suppress contributions from non-flow effects. It was shown in Fig. 8 (top left) that the normalized observable exhibits a clear sensitivity to different parameterizations and the initial conditions, which provides a unique opportunity to discriminate between various possibilities of the detailed setting of of the produced QGP and the initial conditions used in hydrodynamic calculations. On the other hand, normalized is independent of the setting of . In addition, it was demonstrated in Fig. 5 (right) that the normalized , also named in the following text, is compatible with its corresponding observable in the initial state. Thus, the could be taken as golden observable to directly constrain initial conditions without demands for precise knowledge of transport properties of the system (93). Furthermore, none of existing theoretical calculations can reproduce the data, there is still a long way to go for the development of hydrodynamic calculations.
|Observables||Equations||number of particles||Exp.||Th.|
|4||(36)||(93), (92), (80)|
|4||(36)||(93), (92), (80), (110)|
|4||(93), (92), (110)|
|4||(93), (92), (110)|
Predictions of relationship between other harmonics are provided in (93) and shown in Fig. 8. Besides different sensitivities to IC and 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 is observed in central collision using 0.20, while the maximum value is seen in more peripheral collision if 0.08 is used.
Compared to the previous measurements of relationship between flow harmonics investigated using the ESE technique, the 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, , , which is clearly non trivial. Although one did not use the information of symmetry planes in both and studies, recent study just reveals that flow harmonic correlations might be not completely independent on symmetry plane correlations (110). The proportionality relations between symmetric cumulants involving higher harmonics or and symmetry plane correlations is derived, which seems 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 into the nature of fluctuating initial conditions and of the created QGP in heavy ion collisions.
In the past two decades, the underlying of each single harmonic was investigated in great details. However, it is an open question at the moment how the joint underlying , including different order symmetry planes and harmonics, is described, especially if these correlations between different flow harmonics modify the single harmonics . 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 , 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 is an exciting challenge for the theory community.
The author thanks J.J. Gaardhøje, K. Gajdošová, L. Yan, J.Y. Ollitrault and H. Song for the comments on the manuscripts and fruitful discussions. The author is supported by the Danish Council for Independent Research, Natural Sciences, and the Danish National Research Foundation (Danmarks Grundforskningsfond).
- T. D. Lee, Phys. Rev. D 19, 1802 (1979).
- E. V. Shuryak, Phys. Rept. 61, 71 (1980).
- J. -Y. Ollitrault, Phys. Rev. D 46, 229 (1992).
- S. Voloshin and Y. Zhang, Z. Phys. C 70, 665 (1996)
- C. Alt et al. [NA49 Collaboration], Phys. Rev. C 68, 034903 (2003).
- K. H. Ackermann et al. [STAR Collaboration], Phys. Rev. Lett. 86, 402 (2001).
- C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 87, 182301 (2001).
- C. Adler et al. [STAR Collaboration], Phys. Rev. C 66, 034904 (2002).
- J. Adams et al. [STAR Collaboration], Phys. Rev. C 72, 014904 (2005).
- J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 95, 122301 (2005).
- L. Adamczyk et al. [STAR Collaboration], Phys. Rev. C 88, 014902 (2013).
- B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 77, 054901 (2008).
- K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 89, 212301 (2002).
- S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91, 182301 (2003).
- A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 162301 (2007).
- S. Afanasiev et al. [PHENIX Collaboration], Phys. Rev. C 80, 024909 (2009).
- A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 105, 062301 (2010)
- A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 107, 252301 (2011)
- B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. Lett. 89, 222301 (2002);
- B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. C 72, 051901 (2005)
- B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. Lett. 94, 122303 (2005);
- S. Manly et al. [PHOBOS Collaboration], Nucl. Phys. A 774, 523 (2006)
- B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. Lett. 97, 012301 (2006);
- B. Alver et al. [PHOBOS Collaboration], Phys. Rev. Lett. 98, 242302 (2007);
- B. Alver et al. [PHOBOS Collaboration], Phys. Rev. Lett. 104, 142301 (2010);
- B. Alver et al. [PHOBOS Collaboration], Phys. Rev. C 81, 034915 (2010)
- KAamodt et al. [ALICE Collaboration], Phys. Rev. Lett. 105, 252302 (2010).
- K. Aamodt et al. [ALICE Collaboration], Phys. Rev. Lett. 107, 032301 (2011).
- B. Abelev et al. [ALICE Collaboration], Phys. Lett. B 719, 18 (2013).
- B. Abelev et al. [ALICE Collaboration], Phys. Rev. Lett. 111, no. 23, 232302 (2013);
- B. B. Abelev et al. [ALICE Collaboration], Phys. Rev. C 90, no. 5, 054901 (2014)
- B. B. Abelev et al. [ALICE Collaboration], JHEP 1506, 190 (2015).
- J. Adam et al. [ALICE Collaboration], Phys. Rev. C 93, no. 3, 034916 (2016)
- J. Adam et al. [ALICE Collaboration], Phys. Rev. C 93, no. 4, 044903 (2016).
- J. Adam et al. [ALICE Collaboration], Phys. Rev. Lett. 116, no. 13, 132302 (2016).
- J. Adam et al. [ALICE Collaboration], arXiv:1604.07663 [nucl-ex].
- G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 707, 330 (2012)
- G. Aad et al. [ATLAS Collaboration], Phys. Rev. C 86, 014907 (2012)
- G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 725, 60 (2013)
- G. Aad et al. [ATLAS Collaboration], JHEP 1311, 183 (2013)
- G. Aad et al. [ATLAS Collaboration], Phys. Rev. C 90, no. 2, 024905 (2014)
- J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 92, 062301 (2004)
- A. Andronic et al. [FOPI Collaboration], Nucl. Phys. A 679, 765 (2001)
- P. Chung et al. [E895 Collaboration], Phys. Rev. C 66, 021901 (2002)
- G. Aad et al. [ATLAS Collaboration], Eur. Phys. J. C 74, no. 8, 2982 (2014)
- G. Aad et al. [ATLAS Collaboration], Eur. Phys. J. C 74, no. 11, 3157 (2014)
- G. Aad et al. [ATLAS Collaboration], Phys. Rev. C 90, no. 4, 044906 (2014)
- G. Aad et al. [ATLAS Collaboration], Phys. Rev. C 92, no. 3, 034903 (2015)
- S. Chatrchyan et al. [CMS Collaboration], Eur. Phys. J. C 72, 2012 (2012).
- S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 109, 022301 (2012)
- S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. C 87, no. 1, 014902 (2013)
- S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 110, no. 4, 042301 (2013)
- S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 724, 213 (2013)
- S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. C 89, no. 4, 044906 (2014)
- S. Chatrchyan et al. [CMS Collaboration], JHEP 1402, 088 (2014)
- V. Khachatryan et al. [CMS Collaboration], Phys. Lett. B 742, 200 (2015)
- V. Khachatryan et al. [CMS Collaboration], Phys. Rev. Lett. 115, no. 1, 012301 (2015)
- V. Khachatryan et al. [CMS Collaboration], Phys. Rev. C 92, no. 3, 034911 (2015)
- P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, Phys. Lett. B 503, 58 (2001)
- P. F. Kolb, P. Huovinen, U. W. Heinz and H. Heiselberg, Phys. Lett. B 500, 232 (2001)
- M. Luzum and P. Romatschke, Phys. Rev. C 78, 034915 (2008); Phys. Rev. Lett. 103, 262302 (2009)
- H. Song and U. W. Heinz, Phys. Rev. C 77, 064901 (2008)
- H. Song, S. A. Bass, U. Heinz, T. Hirano and C. Shen, Phys. Rev. Lett. 106, 192301 (2011)
- H. Niemi, K. J. Eskola and R. Paatelainen, Phys. Rev. C 93, no. 2, 024907 (2016)
- D. Teaney and L. Yan, Phys. Rev. C 90, no. 2, 024902 (2014)
- U. Heinz and R. Snellings, Ann. Rev. Nucl. Part. Sci. 63, 123 (2013)
- M. Luzum and H. Petersen, J. Phys. G 41, 063102 (2014)
- P. Huovinen, Int. J. Mod. Phys. E 22, 1330029 (2013)
- E. Shuryak, arXiv:1412.8393 [hep-ph].
- H. Song, Pramana 84, 703 (2015)
- K. Dusling, W. Li and B. Schenke, Int. J. Mod. Phys. E 25, no. 01, 1630002 (2016)
- P. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005)
- J. Jia, J. Phys. G 41, no. 12, 124003 (2014)
- R. Snellings [ALICE Collaboration], J. Phys. G 38, 124013 (2011)
- Y. Zhou, PhD Thesis, Utrecht University (2016), CERN-THESIS-2016-005 , doi:1874/325588.
- L. Yan, J. Y. Ollitrault and A. M. Poskanzer, Phys. Rev. C 90, 024903 (2014)
- L. Yan and J. Y. Ollitrault, Phys. Rev. Lett. 112, 082301 (2014)
- L. Yan, J. Y. Ollitrault and A. M. Poskanzer, Phys. Lett. B 742, 290 (2015)
- L. V. Bravina et al., Eur. Phys. J. C 75, no. 12, 588 (2015)
- Y. Zhou, K. Xiao, Z. Feng, F. Liu and R. Snellings, Phys. Rev. C 93, no. 3, 034909 (2016)
- H. Niemi, G. S. Denicol, H. Holopainen and P. Huovinen, Phys. Rev. C 87, 054901 (2013)
- R. S. Bhalerao, M. Luzum and J. -Y. Ollitrault, Phys. Rev. C 84, 034910 (2011)
- R. S. Bhalerao, J. Y. Ollitrault and S. Pal, Phys. Rev. C 88, 024909 (2013)
- D. Teaney and L. Yan, Phys. Rev. C 86, 044908 (2012)
- H. Petersen, G. Y. Qin, S. A. Bass and B. Muller, Phys. Rev. C 82, 041901 (2010)
- Z. Qiu and U. W. Heinz, Phys. Rev. C 84, 024911 (2011)
- U. Heinz, Z. Qiu and C. Shen, Phys. Rev. C 87, no. 3, 034913 (2013)
- F. G. Gardim, F. Grassi, M. Luzum and J. Y. Ollitrault, Phys. Rev. C 87, no. 3, 031901 (2013)
- Y. Zhou [ALICE Collaboration], Nucl. Phys. A 931, 949 (2014)
- I. Kozlov, M. Luzum, G. S. Denicol, S. Jeon and C. Gale, Nucl. Phys. A 931, 1045 (2014)
- A. Bilandzic, C. H. Christensen, K. Gulbrandsen, A. Hansen and Y. Zhou, Phys. Rev. C 89, 064904 (2014)
- R. S. Bhalerao, J. Y. Ollitrault and S. Pal, Phys. Lett. B 742, 94 (2015)
- X. Zhu, Y. Zhou, H. Xu, H. Song, “Correlations of event-by-event flow harmonics in 2.76A TeV Pb-Pb collisions”, to be published.
- J. Schukraft, A. Timmins and S. A. Voloshin, Phys. Lett. B 719, 394 (2013)
- Z. Qiu and U. Heinz, Phys. Lett. B 717, 261 (2012)
- Z. -W. Lin, C. M. Ko, B. -A. Li, B. Zhang and S. Pal, Phys. Rev. C 72, 064901 (2005)
- X. -N. Wang and M. Gyulassy, Phys. Rev. Lett. 86, 3496 (2001)
- B. Zhang, Comput. Phys. Commun. 109, 193 (1998)
- L. -W. Chen and C. M. Ko, Phys. Lett. B 634, 205 (2006)
- B. -A. Li and C. M. Ko, Phys. Rev. C 52, 2037 (1995)
- J. Xu and C. M. Ko, Phys. Rev. C 84, 014903 (2011)
- N. Borghini, P. M. Dinh and J. -Y. Ollitrault, Phys. Rev. C 64, 054901 (2001)
- A. Bilandžić, R. Snellings and S. Voloshin, Phys. Rev. C 83, 044913 (2011)
- G. Agakishiev et al. [STAR Collaboration], Phys. Rev. C 86, 014904 (2012)
- P. Huo, J. Jia and S. Mohapatra, Phys. Rev. C 90, no. 2, 024910 (2014)
- X. N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991).
- M. Gyulassy and X. N. Wang, Comput. Phys. Commun. 83, 307 (1994)
- Y. Zhou [ALICE Collaboration], arXiv:1512.05397 [nucl-ex].
- Y. Zhou, S. S. Shi, K. Xiao, K. J. Wu and F. Liu, Phys. Rev. C 82, 014905 (2010)
- G. Giacalone, L. Yan, J. Noronha-Hostler and J. Y. Ollitrault, arXiv:1605.08303 [nucl-th].