Non-Bessel-Gaussianity and Flow Harmonic Fine-Splitting Hadi Mehrabpour1, 2, ∗ and Seyed Farid Taghavi2, 3, y 1Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran 2School of Particles and Accelerators, Institute for Research in Fundamental Sciences (IPM), P.O. Box 19395-5531, Tehran, Iran 3Physik Department, Technische Universit¨atM¨unchen,Munich, Germany Both collision geometry and event-by-event fluctuations are encoded in the experimentally ob- served flow harmonic distribution p(vn) and 2k-particle cumulants cnf2kg. In the present study, we systematically connect these observables to each other by employing Gram-Charlier A series. We quantify the deviation of p(vn) from Bessel-Gaussianity in terms of flow harmonic fine-splitting. Sub- sequently, we show that the corrected Bessel-Gaussian distribution can fit the simulated data better than the Bessel-Gaussian distribution in the more peripheral collisions. Inspired by Gram-Charlier A series, we introduce a new set of cumulants qnf2kg that are more natural to study distributions near Bessel-Gaussian. These new cumulants are obtained from cnf2kg where the collision geometry effect is extracted from it. By exploiting q2f2kg, we introduce a new set of estimators for averaged ellipticityv ¯2 which are more accurate compared to v2f2kg for k > 1. As another application of q2f2kg, we show we are able to restrict the phase space of v2f4g, v2f6g and v2f8g by demanding the consistency ofv ¯2 and v2f2kg with q2f2kg equation. The allowed phase space is a region such that v2f4g − v2f6g & 0 and 12v2f6g − 11v2f8g − v2f4g & 0, which is compatible with the experimental observations. I. INTRODUCTION tum fluctuations at the initial state. As a result, the observed flow harmonics fluctuate event by event even if It is a well-established picture that the produced mat- we fix the initial geometry of the collision. In fact, the ter in the heavy ion experiment has collective behavior. event-by-event fluctuations are encoded in p(vn) and ex- Based on this picture, the initial energy density mani- perimentally observed flow harmonics as well. It is worth- fests itself in the final particle momentum distribution. while to mention that the observed event-by-event fluc- Accordingly, as the main collectivity consequence, the tuations are a reflection of the initial state fluctuations final particle momentum distribution is extensively stud- entangled with the modifications during different stages ied by different experimental groups in the past years. As of the matter evolution, namely collective expansion and a matter of fact, the experimental groups at Relativistic the hadronization. For that reason, exploring the exact Heavy Ion Collider (RHIC) and Large Hadron Collider interpretation of the flow harmonics is crucial to under- (LHC) measure the flow harmonics vn [1{8], the coeffi- stand the contribution of each stage of the evolution on cients of the momentum distribution Fourier expansion the fluctuations. Moreover, there is no well-established in the azimuthal direction [9{11]. All these observations picture for the initial state of the heavy ion collision so can be explained by several models based on the collec- far. The interpretations of the observed quantities con- tive picture. tain information about the initial state. This information Finding the flow harmonics is not straightforward, be- can shed light upon the heavy ion initial state models too. cause the reaction plane angle (the angle between the ori- According to the theoretical studies, flow harmonic entation of the impact parameter and a reference frame) vn 2k obtained from 2k-particle correlation functions is not under control experimentally. Additionally, the aref different,g and their difference is sourced by non- Fourier coefficients cannot be found reliably due to the flow effects [13, 14] and event-by-event fluctuations as low statistic in a single event. These enforce us to use well [19]. We should say that experimental observations arXiv:1805.04695v2 [nucl-th] 7 Mar 2019 an analysis more sophisticated than a Fourier analysis. show that v2 2 is considerably larger than v2 4 , v2 6 f g f g f g There are several methods to find the flow harmonics and v2 8 . Additionally, all the ratios v2 6 =v2 4 , f g f g f g experimentally, namely Event-Plane method [12], multi- v2 8 =v2 4 and v2 8 =v2 6 are different from unity f g f g f g f g particle correlation functions [13, 14] and Lee-Yang zeros [20, 21]. Alternatively, the distribution p(vn) is approx- [15, 16]. The most recent technique to find the flow har- imated by Bessel-Gaussian distribution (corresponds to monics is using the distribution of flow harmonic p(vn). a Gaussian distribution for vn fluctuations) as a simple This distribution has been obtained experimentally by model [10, 22]. Based on this model, the difference be- employing the unfolding technique [17, 18]. tween v2 2 and v2 4 is related to the width of the v2 It is known that the initial shape of the energy density fluctuations.f g However,f g this model cannot explain the dif- depends on the geometry of the collision and the quan- ference between other v2 2k . f g In the past years, several interesting studies about the non-Gaussian vn fluctuations have been done [23{ ∗ [email protected] 28]. Specifically, it has been shown in Ref. [27] that the y [email protected] skewness of v2 fluctuations is related to the difference 2 v2 4 v2 6 . Also, the connection between kurtosis of initial anisotropies and flow harmonics are used exten- f g − f g v3 fluctuations and the ratio v3 4 =v3 2 has been stud- sively to quantify the initial energy density and final mo- ied in Ref. [28]. It is worth notingf g thatf g the deviation mentum distribution. of vn fluctuations from Gaussian distribution immedi- The initial energy (or entropy) density of a single event ately leads to the deviation of p(vn) from Bessel-Gaussian can be written in terms of the initial anisotropies, namely distribution. In [29], the quantities vn 4 vn 6 and the ellipticity and triangularity. Specifically, ellipticity f g − f g vn 6 vn 8 for a generic narrow distribution is com- and triangularity (shown by "2 and "3, respectively) are puted.f g − f g cumulants of the distribution indicating how much it is In the present work, we will introduce a systematic deviated from a two dimensional rotationally symmetric method to connect vn 2k to distribution p(vn). In Gaussian distribution [30]. f g Sec. II, we have an overview of the known concepts of The final momentum distribution is studied by its cumulants, flow harmonic distributions and their rela- Fourier expansion in the azimuthal direction, tion with the averaged flow harmonicsv ¯n. Sec. III is " 1 # dedicated to the Gram-Charlier A series in which we find 1 dN 1 X = 1 + 2vn cos n(φ n) : (1) an approximated flow harmonic distribution in terms of N dφ 2π n=1 − cn 2k . Specifically for the second harmonics, we show f g that the deviation of p(v2) from Bessel-Gaussianity is In the above, vn and n are unknown parameters that quantified by the fine-splitting v2 2k v2 2` where in n inφ f g − f g can be found easily viav ^n vne = e . Here, the k; ` 2 and k = `. These studies guide us to define a averaging is obtained by using≡ the distributionh i 1 dN in ≥ 6 N dφ new set of cumulants qn 2k where they depend on the f g a given event. The parameterv ^n is called flow harmonic. event-by-event fluctuations only. In Sec. IV, we use the Instead of complex form, we occasionally use the flow new cumulants to introduce more accurate estimations harmonics in Cartesian coordinates, for the average ellipticity. As another application of new cumulants, we use q2 2k to constrain the v2 4 , v2 6 v = v cos n ; v = v sin n : (2) f g f g f g n;x n n n;y n n and v2 8 phase space in Sec. V. We show that the phase f g 1 space is restricted to a domain that v2 4 v2 6 0 Low multiplicity in a single event, randomness of the f g − f g & and 12v2 6 11v2 8 v2 4 & 0. We present the con- reaction plane angle and non-flow effects are the main clusion inf Sec.g − VI.f Theg − supplementaryf g materials can be challenges in extracting the experimental values of flow found in the appendices. We would like to emphasize that harmonics. By defeating these experimental challenges, in the Appendix C, we found a one-dimensional distribu- one would be able to find the distribution p(vn;x; vn;y) tion for p(vn) which is different from that mentioned in for each centrality class which is a manifestation of the Sec. III. Additionally, an interesting connection between event-by-event fluctuations. In finding p(vn;x; vn;y), we p(vn) expansion in terms of cumulants cn 2k and a rel- rotate the system at each single event such that the re- f g atively new concept of the multiple orthogonal polyno- action plane angle ΦRP is set to be zero. In this case, the mials in mathematics is presented in the Appendix D. averaged ellipticity v¯2 v2;x is a manifestation of the geometrical initial ellipticity≡ h fori events in a given cen- trality class irrespective of the fluctuations. In general, II. FLOW HARMONIC DISTRIBUTIONS AND we are able to define averaged flow harmonic v¯n vn;x 2k-PARTICLE CUMULANTS too.