Investigation of the elliptic flow fluctuations of the identified particles using the A Multi-Phase Transport model

A Multi-Phase Transport (AMPT) model is used to study the elliptic flow fluctuations of identified particles using participant and spectator event planes. The elliptic flow measured using the first order spectator event plane is expected to give the elliptic flow relative to the true reaction plane which suppresses the flow fluctuations. However, the elliptic flow measured using the second-order participant plane is expected to capture the elliptic flow fluctuations. Our study shows that the first order spectator event plane could be used to study the elliptic flow fluctuations of the identified particles in the AMPT model. The elliptic flow fluctuations magnitude shows weak particle species dependence and transverse momentum dependence. Such observation will have important implications for understanding the source of the elliptic flow fluctuations.

Many studies of the ultra-relativistic heavy-ion collisions at the Relativistic Heavy Ion Collider and the Large Hadron Collider show that an exotic state of matter named Quark-Gluon Plasma (QGP) is created in these collisions. A large number of studies are focused on identifying the dynamical evolution and the transport properties of the QGP.
In heavy-ion collisions, the produced particle azimuthal anisotropy measurements have been used in various studies to show the viscous hydrodynamic response of the QGP to the initial energy density spatial distribution produced in the early stages of the collisions [1][2][3][4][5][6][7][8][9][10][11][12][13][14]. The azimuthal anisotropy of the particles emitted relative to the reaction plane Ψ R can be described by the Fourier expansion [15,16] of the final-state azimuthal angle φ distribution, The first Fourier harmonic, v 1 , is the directed flow; v 2 is called the elliptic flow, and v 3 is the triangular flow, etc. A wealth of information on the characteristics of the QGP has been gained via the anisotropic flow studies of directed and elliptic flow. Hydrodynamic studies suggest that anisotropic flow stems from the evolution of the medium in the presence of initial-state anisotropies, determined by the eccentricities ε n . The v 2 and v 3 flow harmonics are recognized to be linearly correlated to ε 2 and ε 3 , respectively [7, 28, 32-38]. Therefore for these flow harmonics, v n = κ n ε n , where κ n encodes knowledge about the medium properties such as the specific shear viscosity (η/s) of the QGP. Accurate extraction of η/s requires certain restrictions on the initial-state models employed in such extractions. Such constraints can be achieved via measurements of the flow harmonics and the event-by-event flow fluctuations [39]. Flow fluctuations could be arising from several sources: one of which has attracted considerable attention is the initial eccentricity fluctuations [40][41][42]. Recent theoretical studies have begun to take into account initial conditions that include energy density fluctuations, initial flow [13,37,43], and the full shear stress tensor [44] at µ B = 0 and at µ B > 0 [45][46][47][48]. Also, the partonic structure inside the nucleons has been considered in Reference [49].
Recently, Reference [50] presented more realistic eventby-event fluctuating initial conditions, Initial Conserved Charges in Nuclear Geometry (ICCING), of not only the initial energy density profile but also the initial conserved charges of baryon number (B), strangeness (S), and electric charge (Q) density distributions. This work pointed out that while baryon number and electric charge have almost the same geometries to the energy density profile, the initial strangeness distribution is considerably more eccentric. Such an effect predicts that the elliptic flow fluctuations will be larger for the strange and multi-strange hadrons. This effect can be detected experimentally via studying the elliptic flow fluctuations of the identified hadrons.
The ratio between four-particles elliptic flow, v 2 {4}, and the two-particles elliptic flow, v 2 {2}, is often used to estimate the strength of the elliptic flow fluctuations as a fraction of the measured flow harmonic strength [51,52]. However, important caveats to studying the elliptic flow fluctuations using (v 2 {4}/v 2 {2}) for the identified hadrons are, first, the demand for high statistical power, and second, the multi-strange hadron identification process [53]. Consequently, the ratio of v 2 {4}/v 2 {2} is of limited experimental use for carrying out these investi-gations for the multi-strange hadrons.
In this work, we investigate an alternative validation scheme, which employs the use of the first-order spectator event plane ,Ψ SP 1 , along with the second-order event plane Ψ EP 2 to study the elliptic flow fluctuations of the identified hadrons. Here, the underlying notion is that v SP 2 (with respect to the spectator first-order event plane) will reduce the elliptic flow fluctuations due to the strong correlations between the Ψ SP 1 and the true reaction plane. Therefore, the ratio v SP 2 /v EP 2 is expected to reflect the elliptic flow fluctuations.
For RHIC highest energy and using the STAR detector, we propose a similar investigation to be performed using the first-order spectator event plane from spectator neutrons, measured by the zero-degree calorimeters (ZDC) [54] and the second-order event plane using the new installed Event-Plane-Detector (EPD) [55]. Consequently, we think that conducting a similar experimental study will reveal important information about the elliptic flow fluctuations and will shed light on the ICCING scenario suggested in Reference [50].

METHOD
The current study is conducted with simulated events for Au+Au collisions at √ s N N = 200 GeV, collected using the AMPT [56] model with the string-melting mechanism and hadronic cascade on. The AMPT model, which has been widely employed to study relativistic heavy-ion collisions [56-60, 60-62], includes four main dynamical components: initial condition, parton cascade, hadronization, and hadronic rescatterings. The initial conditions take into account soft string excitations and the phase space distributions of minijet partons, which are produced by the Heavy-Ion Jet Interaction Generator model (HIJING) [63] in which the Glauber model with multiple nucleon scatterings are used to define the heavy-ion collisions initial state.
The partons scatterings are handled according to the Zhang's Parton Cascade (ZPC) model [64], which contain only two-body elastic scatterings with a cross-section defined as: where α s = 0.47 is the strong coupling constant, µ is the screening mass and s and t are the Mandelstam variables.
In the AMPT with the string-melting mechanism, the excited strings and minijet partons are melted into partons.
The partons scatterings will lead to local energy density fluctuations, which are equivalent to the local transverse density of participant nucleons.
In the string-melting version and when partons stop interacting with each other, a quark coalescence model is used to couple partons into hadrons. Consequently, the partonic matter is then converted into hadronic matter and the hadronic interactions are given by the A Relativistic Transport (ART) model [65], which incorporates both elastic and inelastic scatterings for baryon-baryon, baryon-meson, and meson-meson interactions.
In this work, the centrality intervals are defined by selecting the impact parameter distribution, then the AMPT events are analyzed using (i) the event plane method and (ii) the multi-particle cumulant technique [66][67][68][69]. Using both methods, particle of interest (POI) comes from pseudorapidities |η| < 1, which matches the STAR experiment pseudorapidity acceptance, and with transverse momentum 0.1 < p T < 4.0 GeV/c.
The second-order event plane (Ψ EP 2 ), is estimated from the azimuthal distribution of final-state particles. The elliptic flow that will be obtained using this method will then be corrected with the corresponding event plane resolution (Res( Ψ EP 2 )) [16]. The Ψ EP 2 is reconstructed in a pseudorapidity range of 2.5 < |η| < 4.5, which matches the STAR experiment EPD acceptance, and 0.1 < p T < 2.0 GeV/c: where φ i is the final-state azimuthal angle of particle i, and ω i is its weight. The weight is chosen to be equal to p T . Also, the first order spectator plane Ψ SP 1 is constructed using the AMPT spectator x and y position information. Using the spectator or the event planes we can give the elliptic flow as: where Res( Ψ EP 2 ) and Res( Ψ SP 1 ) represent the resolution of the event planes. The event planes resolution is calculated using the two-subevent method [16].
On the other hand, the standard (subevents) cumulant methods framework is discussed in References [66][67][68][69]. In the standard cumulant method, the n-particle cumulants are constructed using particles from the |η| < 1.0 acceptance. Thus the constructed two-and four-particle correlations can be written as: where, represents the average over all particles in a single event, and then in average over all events, n is the harmonic number and ϕ i expresses the azimuthal angle Ratio of the i th particle. Then the four-particle elliptic flow harmonic can be given as: In general, when the flow fluctuation σ is smaller than the true reaction plan elliptic flow v 2 one can write [70,71]: The reliability of this elliptic flow fluctuations extraction will depend on the strength of the correlations between the spectator plane and the reaction plane. Ratio and v 2 {4}. Also a good agreement (within the errors) has been observed between the v 2 {4} and v SP 2 . The ratio v SP 2 /v EP 2 , presented in panel (b) presents the strength of the elliptic flow fluctuations which shows no p T dependence, consistent with the preliminary STAR measurements [77].
The centrality dependence of the identified particles v EP