Importance of correlations and fluctuations on the initial source eccentricity in high-energy nucleus-nucleus collisions

In relativistic heavy-ion collisions, anisotropic collective flow is driven, event by event, by the initial eccentricity of the matter created in the nuclear overlap zone. Interpretation of the anisotropic flow data thus requires a detailed understanding of the effective initial source eccentricity of the event sample. In this paper, we investigate various ways of defining this effective eccentricity using the Monte Carlo Glauber (MCG) approach. In particular, we examine the participant eccentricity, which quantifies the eccentricity of the initial source shape by the major axes of the ellipse formed by the interaction points of the participating nucleons. We show that reasonable variation of the density parameters in the Glauber calculation, as well as variations in how matter production is modeled, do not significantly modify the already established behavior of the participant eccentricity as a function of collision centrality. Focusing on event-by-event fluctuations and correlations of the distributions of participating nucleons, we demonstrate that, depending on the achieved event-plane resolution, fluctuations in the elliptic flow magnitude ${v}_{2}$ lead to most measurements being sensitive to the root-mean-square rather than the mean of the ${v}_{2}$ distribution. Neglecting correlations among participants, we derive analytical expressions for the participant eccentricity cumulants as a function of the number of participating nucleons, ${N}_{\mathrm{part}}$, keeping nonnegligible contributions up to $\mathcal{O}(1/{N}_{\mathrm{part}}^{3})$. We find that the derived expressions yield the same results as obtained from mixed-event MCG calculations which remove the correlations stemming from the nuclear collision process. Most importantly, we conclude from the comparison with MCG calculations that the fourth-order participant eccentricity cumulant does not approach the spatial anisotropy obtained assuming a smooth nuclear matter distribution. In particular, for the Cu+Cu system, these quantities deviate from each other by almost a factor of 2 over a wide range in centrality. This deviation reflects the essential role of participant spatial correlations in the interaction of two nuclei.