# Third Harmonic Flow of Charged Particles in Au+Au Collisions at GeV

###### Abstract

In this proceedings, we report measurements of the third harmonic coefficient of the azimuthal anisotropy, , known as triangular flow. The analysis is for charged particles near midrapidity in Au+Au collisions at = 200 GeV, based on data from the STAR experiment at the Relativistic Heavy Ion Collider. Triangular flow as a function of centrality, pseudorapidity and transverse momentum are reported using various methods, including a study of the signal for particle pairs as a function of their pseudorapidity separation. Results are compared with other experiments and model predictions.

## 1 Introduction

The study of azimuthal anisotropy, based on Fourier coefficients, is widely recognized as an important tool to probe the hot, dense matter created in heavy ion collisions [1].

(1) |

where denotes the azimuthal angle of the outgoing particles, is the orientation of the reaction plane and denotes the harmonic.

The first harmonic coefficient , called directed flow, and the second harmonic coefficient , called elliptic flow, have been extensively studied both experimentally [2] and theoretically, while higher even-order harmonics have also garnered some attention [3]. In contrast, odd harmonics of order three and above were overlooked until recently, because in a picture with smooth initial overlap geometry, it had been assumed that higher-order odd harmonics are required to be zero by symmetry. Now it is realized that event-by-event fluctuations break this symmetry [4, 5, 6]. As a consequence, higher-order odd harmonics carry valuable information about the initial state of the colliding system [7, 8, 9, 10, 11, 12, 13]. The third harmonic coefficient — related to a collective motion known as triangular flow — is thus a new tool to study initial-state fluctuations, and the subsequent evolution of the collision system. Theoretical studies suggest that is more sensitive to viscous effects than because finer details of the higher harmonics are smoothed more by viscosity. We present measurements of vs. the pseudorapidity separation = between the two particles, fit with narrow and wide Gaussians, and present results for the wide Gaussian as a function of and centrality.

## 2 Methods and Analysis

Azimuthal correlations not related to the initial geometry and reaction plane orientation which can arise from resonances, jets, strings, quantum statistics effects, final state interactions (particularly Coulomb effects), and momentum conservation are refereed as non-flow correlations. Two-particle correlations which are sensitive to non-flow correlations, and multiparticle correlations which are less sensitive to non-flow correlations, along with event plane methods have been extensively used in elliptic flow measurements for the past couple of decades. These methods can be extended to measure triangular flow. The event plane method and the two-particle cumulant method are used in the present study.

### 2.1 Event Plane Methods

In the standard event plane method [1] for , we reconstruct a third-harmonic event plane from Time Projection Chamber (TPC ()) tracks and also from Forward TPC( ) tracks. The event plane vector and the event plane angle from the third harmonic of the particle azimuthal distribution are defined by the equations

(2) |

(3) |

(4) |

In event plane calculations, tracks have a weighting factor in units of GeV for GeV , and w = 2 for 2 GeV. This helps to maximize the event plane resolution and also helps to reduce the contribution from jets in event plane calculations. Although the STAR detector has good azimuthal symmetry, small acceptance effects in the calculation of the event plane azimuth were removed by the method of shifting [14]. When using the TPC event plane, we used the sub-event method with an additional gap of 0.05.

(5) |

This avoids self-correlations because the particles and the event plane are in opposite hemispheres. When using the FTPCs, we used the full event plane from both FTPCs [1].

(6) |

This introduces a large gap between the particles and the event plane. Since there is no overlap between the coverage of the TPC and FTPCs, there is no possibility of self-correlation when using the FTPC event plane.

### 2.2 Two-Particle Correlation Method

We study vs for two particles with indices and to understand the dependence of the triangular flow signal and to distinguish among different sources of non-flow correlations. This distribution of vs. can be well described by wide and narrow Gaussian peaks as shown in Fig. 1 for two centrality intervals. The narrow Gaussian is identified with short range non-flow correlations like Bose-Einstein correlations, resonance decay, Coulomb interactions, and effects from track merging. The narrow peak disappears above GeV/ so it is unlikely to be from jet correlations. The wide Gaussian represented by is the signal of interest and its fit parameters are used to calculate as a function of centrality and transverse momentum. For integrated over and we have

(7) |

where equals when weighted with the number of particle pairs, or is set to 1 for unit weight. The quantity can be obtained from the scalar product [15] relation

(8) |

where the particle is selected from the bin of interest.

## 3 Results and Discussion

Figure 2 shows the width and amplitude of the wide Gaussian as a function of centrality. We observe that the width peaks in mid-central collisions whereas the amplitude increases from central to peripheral collisions. Figure 3 shows the width and amplitude of the wide Gaussian as a function of for centralities 0–5% and 30–40%. Above 0.8 GeV/, the distribution can be described by a single wide Gaussian. The amplitude increases with and then saturates around 3 GeV/. The dependence of the width seems to depend on centrality, with the 0–5% most central data showing first an increase and then a gradual decrease, while at 30–40% centrality, the data appear to gradually decrease as a function of . Figure 4 shows the centrality dependence for -integrated from several different analyses: from Eq. 7 and Fig. 1 for the wide Gaussian, measured with third harmonic event plane reconstructed in the TPC, two-particle cumulants with a minimum pseudorapidity separation between particles of one unit, and where is measured relative to the third harmonic event plane reconstructed in the FTPCs.

Clearly the various analysis methods for differ greatly in Fig. 4. For the wide Gaussian and the TPC event plane, the results are similar, suggesting that the narrow Gaussian is not very important. When a large is required, the results decrease, especially for peripheral collisions. The variation between results shown in Fig. 4 is caused by the dependence as shown in Fig. 5. Non-flow effects, which are assumed to be largest when is small, are expected to contribute less to the methods sampling a larger . This hypothesis is consistent with previous studies of elliptic flow, based on two-particle correlations, where the corresponding wide Gaussian was ascribed to mini-jet correlations [27]. The decrease with of may arise from decreasing non-flow or from decreasing initial-state density fluctuations; it is argued that those fluctuations decrease with increasing separation [25]. Since depends so strongly on the separation of the particles, one must always quote for each measurement, and compare results to models with the same . If is related to the initial eccentricity fluctuations, then the reduction of at large would presumably require a decrease of the initial state fluctuations at large rapidity separations. Recent work has found such a decoherence effect with a hadron and parton cascade model [25].

The non-flow contributions due to short-range correlations are effectively suppressed by using the wide Gaussian or by the -gap method. We studied the influence of different pseudorapidity gaps. However, Fig. 5 shows that we did not find a pattern where stabilized at a constant value for large . We compared like and unlike charge-sign combinations, since they have different contributions from resonance decays, fluctuations, and final state interactions. The large discrepancy between the methods may have its origin in the dependence of the fluctuations, which seem to decrease with increasing separation. Thus it is not clear if one should extrapolate to large to avoid non-flow, or to small to measure all the fluctuations [25]. More theoretical input is necessary to completely understand the dependence of this signal. Figure 6 shows the dependence of using event plane methods. For particles in the TPC using the sub-event method, is somewhat peaked at mid-rapidity. With the event plane in the FTPCs, there is a large gap and is flat at all centralities. This implies that acceptance effects at the edges of the TPC are not important. Thus, even though a large means that one of the particles is probably at large , this apparently does not strongly influence the dependence.

In Fig. 7, STAR results with the event plane in the TPC are very similar to those of PHENIX [18]. This is surprising because the mean of their RXN detector is larger than in the case of the sub-events in the STAR TPC. However, the STAR FTPC results are lower than the PHENIX results. This is expected, because the mean is considerably larger in the FTPC than in the PHENIX RXN detector.

Third harmonic flow coefficient has been studied in event-by-event ideal hydrodynamics [12] with MC Glauber initial conditions. These authors concluded that instead of averaged initial conditions, event-by-event calculations are necessary to compare with experimental data. The first prediction of with viscous hydro was Ref. [10]. In Fig. 8, and as a function of transverse momentum obtained with TPC sub event plane method are compared with several models for 0–5% and 30–40% central collisions. The experimental results for the TPC sub-event plane method which may have small non flow contribution are shown because they eliminate the short-range correlations but yet have small like the theory calculation. The specific models are the viscous hydrodynamic model of Ref. [12], where the ratio of viscosity to entropy is 0.08 and 0.16, and the AMPT model [13]. AMPT model results are for and from Ref [4] compared with 0-5% and 30-40% central data . Predictions for from Parton-Hadron String Dynamics [20] at 30–40% centrality for have been made by the sub-event method with the event planes at , and are also plotted in Fig. 8 (lower right panel).

Elliptic flow() results are mostly described by ideal hydrodynamics in the case of the most central collisions, and by 0.08 in the case of mid-central collisions. We find that the third harmonic coefficient( ) results are also described by this model with similar viscosities. The PHSD model also agrees with data.

## 4 Summary

We present STAR’s measurements of third harmonic flow of charged particles from Au+Au collisions at GeV as a function of pseudorapidity, transverse momentum, and centrality. We report result rom a two-particle method for particle pairs with an gap or fit with a wide Gaussian in pseudorapidity separation, as well as from the standard event-plane method with an event plane near midrapidity or at forward rapidity. Short-range correlations are eliminated either by an gap or by eliminating the narrow Gaussian in pseudorapidity separation. The measured values of continuously decrease as the mean pseudorapidity separation() of the particles increase within the range observable by STAR. It is not known whether this decrease is due to a decrease in non-flow correlations or due to a decrease in fluctuations. We observe that increases with transverse momentum before it levels-off above 3 GeV, similar to the case of elliptic flow. Our results are mostly described by hydrodynamic models with small viscosity with MC Glauber initial conditions.

## 5 Acknowledgments

We thank the RHIC Operations Group and RCF at BNL, the NERSC Center at LBNL and the Open Science Grid consortium for providing resources and support. This work was supported in part by the Offices of NP and HEP within the U.S. DOE Office of Science, the U.S. NSF.

## 6 References

## References

- [1] A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, 1671 (1998).
- [2] G. Agakishiev et al. (STAR Collaboration), Phys. Rev. C 85, 014901(2012). Y. Pandit (for the STAR Collaboration) J. Phys. Conf. Ser.316, 012001 ( 2011). B. I. Abelev et al. (STAR Collaboration), Phys. Rev. Lett. 101, 252301 (2008). B. I. Abelev et al. (STAR Collaboration), Phys. Rev. C 81, 044902 (2010).
- [3] P. F. Kolb, Phys. Rev. C 68, 031902(R) (2003); J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 062301 (2004).
- [4] B. Alver and G. Roland, Phys. Rev. C 81, 054905 (2010).
- [5] D. Teaney and L. Yan, Phys. Rev. C83, 064904 (2011).
- [6] P. Sorensen, B. Bolliet, A. Mocsy, Y. Pandit and N. Pruthi, Phys. Lett. B 705 , 71( 2011).
- [7] A. P. Mishra, R. K. Moohapatra, P. S. Saumia and A. M. Shrivastava, Phys. Rev. C 81, 034903 (2010).
- [8] G.-Y. Qin, H. Petersen, S. A. Bass, and B. Müller, Phys. Rev. C 82 , 064903 (2010).
- [9] J. L. Nagle and M. P. McCumber, Phys. Rev. C83, 044908 (2011).
- [10] B. H. Alver, C. Gombeaud, M. Luzum and J. Y. Ollitrault, Phys. Rev. C 82, 034913 (2010).
- [11] H. Petersen, G.-Y. Qin, S. A. Bass, and B. Müller, Phys. Rev. C 82, 041901(R) (2010).
- [12] B. Schenke, S. Y. Jeon, and C. Gale, Phys. Rev. Lett. 106, 042301 (2011).
- [13] J. Xu and C. M. Ko, Phys. Rev. C 84 014903 (2011).
- [14] J. Barrette et al. (E877 Collaboration), Phys. Rev. C 56, 3254 (1997).
- [15] C. Alder et al. (STAR Collaboration), Phys. Rev. C 66, 034904 (2002); N. Borghini, P. M. Dinh, and J. Y. Ollitrault, Phys. Rev. C 63, 054906 (2001).
- [16] J. -Y. Ollitrault, A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 80, 014904 (2009).
- [17] P. Bozek and W. Broniowski, arXiv:1204.3580 [nucl-th].
- [18] A. Adare et al. (PHENIX Collaboration), Phys. Rev. Lett. 107, 252301 (2011).
- [19] J. Adams et al. (STAR Collaboration), Nucl. Phys. A 757, 102 (2005).
- [20] V. P. Konchakovski, E. L. Bratkovskaya, W. Cassing, V. D. Toneev, S. A. Voloshin, and V. Voronyuk, Phys. Rev. C 85, 044922 (2012).
- [21] A. Bilandzic, R. Snellings, and S. Voloshin, Phys. Rev. C 83, 044913 (2011).
- [22] J. Y. Ollitrault, Phys. Rev. D 46, 229 (1992).
- [23] J. Adams et al. (STAR Collaboration), Phys. Rev. C 72, 014904 (2005).
- [24] M. L. Miller, K. Reygers, S. J. Sanders and P. Steinberg, Ann. Rev. Nucl. Part. Sci. 57, 205 (2007); B. Alver et al., Phys. Rev. C 77, 014906 (2008).
- [25] K. Dusling, F. Gelis, T. Lappi and R. Venugopalan, Nucl. Phys. A 836, 159 (2010).
- [26] H. Petersen, C. Greiner, V. Bhattacharya and S. A. Bass, arXiv:1105.0340 [nucl-th].
- [27] M. Daugherity (STAR Collaboration), J. Phys. G 35, 104090 (2008); D. Kettler (STAR Collaboration), PoS C ERP2010, 011 (2010); T. A. Trainor and D. T. Kettler, Phys. Rev. C 83, 034903 (2011).