Measurement of the flow harmonic correlations in pp, p+Pb and low multiplicity Pb+Pb collisions with the ATLAS detector at the LHC

XXVIIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2018) Measurement of the flow harmonic correlations in pp,p+Pb and low multiplicity Pb+Pb collisions with the ATLAS detector at the LHC Dominik Derendarz on behalf of the ATLAS Collaboration Institute of Nuclear Physics PAS, ul. Radzikowskiego 152, 31-342 Krak´ow, Poland Abstract Recent measurements of the correlations between flow harmonics obtained using four-particle symmetric cumulants and three-particle asymmetric cumulants with the ATLAS detector at the LHC are described. The data sets of pp,p+Pb and peripheral Pb+Pb collisions at various energies are analyzed, aiming to probe the long-range collective nature of multi-particle production in small systems. The sensitivity of the standard cumulant method to non-flow correlations is investigated by introducing the subevents method. A systematic reduction of non-flow effects is observed when using the two-subevent method. Further reduction is observed with the three-subevent method that is consistent with the results obtained with the four-subevent one. A negative correlation between v2and v3and a positive correlation between v2 and v4, for all studied collision systems and over full multiplicity range, is observed. The correlation strength computed as symmetric cumulants normalized by the v2 nis similar for all collision systems and weakly depends on multiplicity. These measurements provide new evidence for long-range multi-particle collectivity in small collision systems and quantify the nature of its event-by-event fluctuations. Keywords: QGP, heavy ion collisions, small systems, collectivity, symmetric cumulants, subevent cumulants 1. Introduction Collective phenomena have been key observables in studies of the properties of the quark-gluon plasma (QGP). Typically the azimuthal distribution of charged particles is characterized by the Fourier decomposition of dN/dφ∝1+2Σ∞ n=1vncos n(φ−Φn), where vnand Φnrepresent the magnitude and the event-plane angle of the n-th flow harmonic. It is understood that in high energy A+Acollisions vncoefficients carry information about the hydrodynamic medium response to the asymmetric initial conditions [1]. The presence of multi-particle correlations in small collision systems such as p+Aand pp collisions [2, 3] open up the discussion on whether the physical processes responsible for the observed long range rapidity correlations Available online at www.sciencedirect.com Nuclear Physics A 982 (2019) 479–482 0375-9474/© 2018 Published by Elsevier B.V. www.elsevier.com/locate/nuclphysa https://doi.org/10.1016/j.nuclphysa.2018.09.076 This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). and their azimuthal structure are of the same nature in small systems as in heavy ion collisions. The ”symmetric cumulants”, scn,m{4}=v2 nv2 m−v2 nv2 m, are novel observables that measure the correlation between two flow harmonics of different orders, vnand vm, using four-particle correlations, and were shown in A+A collisions to have better sensitivity to initial conditions and QGP transport properties than the magnitudes of the single harmonics vn[4]. The measurement of symmetric cumulants in small systems could shed more light on the properties of the medium created in those collisions. However, experimentally the measurement of flow effects in small systems is sensitive to non-flow arising from correlations among limited number of particles due to, for example, resonance decays or jets production. It was clearly shown that the subevent cumulant method effectively suppresses non-flow in small systems [3] and thus would be essential in the unbiased measurement of the correlation between flow harmonics. The goal of the analysis is to study the influence of non-flow on the scn,m{4}measurement and to compare the vn−vmharmonics correlation across small and large collision systems. Additionally, the observable that involves correlations among three-particles, termed as ”asymmetric cumulants”, acn{3}=v2 nv2ncos 2n(Φn− Φ2n), is proposed to study the correlation of vn−v2nharmonics. In contrast to scn,m{4}, which measures only the correlation between the magnitude of vnand vm,ac n{3}is sensitive to both the flow magnitude vn and the flow phase Φn. 2. Data analysis This analysis is based on ATLAS [5] data corresponding to integrated luminosities of 0.9 pb−1of pp data recorded at √s=13 TeV, 28 nb−1of p+Pb data recorded at √sNN =5.02 TeV and 7 μb−1of Pb+Pb data at √sNN =2.76 TeV. The mathematical framework for cumulants is based on the Q-cumulant [6] and has been extended to the case of symmetric cumulants [7]. The formalism of the subevent method used for both standard cumulants and symmetric cumulants has been recently presented in [8]. Symmetric and asymmetric cumulants are constructed from the two-, threeand four-particle correlations in the following way: acn{3}={3}n ,scn,m{4}={4}n,m − {2}n{2}m .(1) The double angular brackets () indicate that the averaging procedure has been performed in two steps. In the first step, the average is taken over all distinct particle multiplets in an event. In the second step, the average is calculated over an ensemble of events with similar Nch,efficiency-corrected average number of charged particles. In the standard cumulant method all k-particle multiplets involved in {k}n and {k}n,m are selected using tracks from the entire detector pseudorapidity acceptance (|η|<2.5). While four-particle correlations are sufficient to suppress non-flow effects in A+Acollisions, this is not the case in small systems. To suppress further the non-flow correlations that typically involve a few particles within a localized region in η, the sample of charged particle tracks is divided into subevents, each covering a unique ηinterval. The multi-particle correlations are then constructed by only correlating tracks between different subevents. In the case of using two subevents, labeled aand b, the detector ηacceptance is split in half and the per-event k-particle azimuthal correlations are constructed by correlating tracks from subevent awith tracks from subevent b. Here the acn{3}and scn,m{4}are defined as [8]: ac2a|b n{3}={3}n2a|b,sc2a|2b n,m{4}={4}n,m2a|2b−{2}na|b{2}ma|b.(2) Similarly, in the three-subevent method, the ηacceptance is divided in three equal subevents a,band c, which gives: aca,b|c n{3}={3}na,b|c,sca,b|2c n,m{4}={4}n,ma,b|2c−{2}na|c{2}mb|c.(3) In the case of the four-subevent method, each of the subevents covers one fourth of the ηacceptance, and the corresponding cumulant is defined as: sca,b|c,d n,m{4}={4}n,ma,b|c,d−{2}na|c{2}mb|d.(4) In order to improve statistical precision, results obtained from all combinations of subevents a,b,c(and din the case of scn,m{4}) are averaged to obtain the final values. The complete analysis description can be found in the ATLAS publication [9]. D. Derendarz / Nuclear Physics A 982 (2019) 479–482480 Fig. 1. The symmetric cumulants sc2,3{4}(left column), sc2,4{4}(middle column) and asymmetric cumulant ac2{3}(right column) as a function of Nchobtained for pp collisions (top row), p+Pb collisions (middle row) and low-multiplicity Pb+Pb collisions (bottom row) [9]. In each panel, the results are obtained from the standard method (filled symbols), the two-subevent method (open circles), three-subevent method (open squares) and in the case of symmetric cumulants four-subevent method (open diamonds). The error bars and shaded boxes represent the statistical and systematic uncertainties, respectively. 3. Results Figure 1 compares the sc2,3{4},sc 2,4{4}and ac2{3}(columns) obtained from the standard, two-, threeand four-subevent methods in pp,p+Pb and Pb+Pb collisions (rows). Values of the symmetric and asymmetric cumulants from the standard method in pp collisions (top row) are significantly larger than the values from the subevent methods that are in good agreement with each other. This is due to the strong influence of non-flow effects. In the most extreme case, the sc2,3{4}obtained with the standard method shows a positive correlation of v2and v3that is reduced to an anti-correlation after suppressing the non-flow with the subevent technique. In p+Pb collisions (middle row) the most pronounced difference between the standard and subevent cumulants is seen for the events with Nch<140. Again, the sc2,3{4}at lower multiplicities shows a positive correlation between v2and v3that changes sign when the subevent method is applied. On the other hand, the v2−v4correlation is positive over the entire measured Nchrange, as measured by sc2,4{4} and ac2{3}. Interestingly, for Nch>200 values from the standard method are systematically larger than those obtained from the subevent method. Similar behaviour is observed in Pb+Pb collisions (bottom row) in the high Nchregion. This difference may be caused by non-flow contributions or could be attributed to flow decorrelation effects that are more significant in the case of the v4harmonic [10]. No differences are observed between the threeand four-subevent methods, indicating that with three subevents most of the non-flow effects are suppressed. As can be seen in Figure 1, ac2{3}reproduces all features observed in sc2,4{4}, while it is measured with much better statistical precision thanks to the larger magnitude and larger number of three particle multiplets. A direct comparison of the cumulants in the three collision systems using the three-subevent method is shown in Figure 2. The top row shows the results for sc2,3{4},sc 2,4{4}and ac2{3}. An anti-correlation D. Derendarz / Nuclear Physics A 982 (2019) 479–482 481 Fig. 2. The top row shows the Nchdependence of sc2,3{4}(left column), sc2,4{4}(middle column) and ac2{3}(right column) obtained for pp collisions (solid circles), p+Pb collisions (open circles) and low-multiplicity Pb+Pb collisions (open squares) [9]. The bottom row shows the normalized versions of symmetric and asymmetric cumulants. The error bars and shaded boxes represent the statistical and systematic uncertainties, respectively. between v2and v3and a correlation between v2and v4is observed in all systems. The strength of the correlation in the Nchrange covered by pp collisions is approximately the same across all systems and, for higher Nch, increases faster for Pb+Pb than for p+Pb collisions. To further investigate the strength of the correlation across different colliding systems and separate the contribution from the magnitudes of the flow harmonics, the sc2,3{4},sc 2,4{4}and ac2{3}are normalized by the values of the corresponding vncoefficients obtained from the two-particle correlation method [9] (Figure 2 bottom row). The normalized cumulants show much weaker multiplicity dependence than the unnormalized cumulants and are similar for different collision systems even at high Nch. The only exception is the normalized sc2,3{4}in pp interactions, which is significantly different from the sc2,3{4}measured in p+Pb and Pb+Pb collisions. As was already noticed in [11], the template fit with perhipheral subtraction tends to underestimate the odd harmonics, in this case v3. The results presented here provide a further evidence that azimuthal correlations in small systems behave similarly to heavy ion collisions. 4. Acknowledgment This work was supported in part by the National Science Centre, Poland grant 2016/23/B/ST2/00702 and by PL-Grid Infrastructure. References [1] U. Heinz, R. Snellings, Ann. Rev. Nucl. Part. Sci. 63 (2013) 123–151. arXiv:1301.2826 [nucl-th], doi:10.1146/annurev-nucl102212-170540. [2] ATLAS Collaboration, Eur. Phys. J. C77 (6) (2017) 428. arXiv:1705.04176 [hep-ex], doi:10.1140/epjc/s10052-017-4988-1. [3] ATLAS Collaboration, Phys. Rev. C97 (2) (2018) 024904. arXiv:1708.03559 [hep-ex], doi:10.1103/PhysRevC.97.024904. [4] ALICE Collaboration, Phys. Rev. Lett. 117 (2016) 182301. arXiv:1604.07663 [nucl-ex], doi:10.1103/PhysRevLett.117.182301. [5] ATLAS Collaboration, JINST 3 (2008) S08003. doi:10.1088/1748-0221/3/08/S08003. [6] A. Bilandzic, R. Snellings, S. Voloshin, Phys. Rev. C83 (2011) 044913. arXiv:1010.0233 [nucl-ex], doi:10.1103/PhysRevC.83.044913. [7] A. Bilandzic, C. H. Christensen, K. Gulbrandsen, A. Hansen, Y. Zhou, Phys. Rev. C89 (6) (2014) 064904. arXiv:1312.3572 [nucl-ex], doi:10.1103/PhysRevC.89.064904. [8] J. Jia, M. Zhou, A. Trzupek, Phys. Rev. C96 (3) (2017) 034906. arXiv:1701.03830 [nucl-th], doi:10.1103/PhysRevC.96.034906. [9] ATLAS Collaboration, arXiv:1807.02012 [nucl-ex], submitted to Phys. Lett. B. [10] ATLAS Collaboration, Eur. Phys. J. C78 (2) (2018) 142. arXiv:1709.02301 [nucl-ex], doi:10.1140/epjc/s10052-018-5605-7. [11] ATLAS Collaboration, Phys. Rev. C96 (2) (2017) 024908. arXiv:1609.06213 [nucl-ex], doi:10.1103/PhysRevC.96.024908. D. Derendarz / Nuclear Physics A 982 (2019) 479–482482