Global polarization of Lambda and (Lambda)over-bar hyperons in Pb-Pb collisions at root s(NN)=2.76 and 5.02 TeV

The global polarization of the (cid:2) and (cid:2) hyperons is measured for Pb-Pb collisions at √ s NN = 2 . 76 and 5.02 TeV recorded with the ALICE at the Large Hadron Collider (LHC). The results are reported differentially as a function of collision centrality and hyperon’s transverse momentum ( p T ) for the range of centrality 5–50%, 0 . 5 < p T < 5 GeV/ c , and rapidity | y | < 0 . 5. The hyperon global polarization averaged for Pb-Pb collisions at √ s NN = 2 . 76 and 5.02 TeV is found to be consistent with zero, (cid:4) P H (cid:5) (%) ≈ 0 . 01 ± 0 . 06 (stat.) ± 0 . 03 (syst.) in the collision centrality range 15–50%, where the largest signal is expected. The results are compatible with expectations based on an extrapolation from measurements at lower collision energies at the Relativistic Heavy Ion Collider, hydrodynamical model calculations, and empirical estimates based on collision energy dependence of directed ﬂow, all of which predict the global polarization values at LHC energies of the order of 0 . 01%.


I. INTRODUCTION
The system created in a noncentral nucleus-nucleus collision might retain a significant fraction of the large orbital angular momentum of the colliding nuclei. Due to the spinorbit coupling, particles produced in such a collision can become globally polarized [1][2][3][4]. The global polarization is a phenomena when spins of all emitted final-state particles in a noncentral nucleus-nucleus collision are aligned along one preferential direction defined by the initial angular orbital momentum. Its measurement provides important information about the initial conditions and dynamics of the quark-gluon plasma (QGP), as well as the hadronization process [5][6][7]. The global polarization for a specific particle and corresponding antiparticle is expected to be very similar, or even identical, in a system with small or zero baryon chemical potential.
High-energy heavy-ion collisions are also characterized by ultrastrong magnetic fields [8,9], which on average are aligned with the direction of the angular momentum. While the peak values of the magnetic fields, which reach up to 10 18 Gauss, can be estimated rather accurately [8,9], the time evolution, which depends on the QGP electric conductivity, is practically unknown. These fields can also contribute to the global polarization, but their action on particles and antiparticles is expected to be in an opposite direction. Thus the measurement of the splitting between particle and antiparticle global polarizations provides very valuable information about the QGP properties.
The measurements of the global polarization P H for spin one half ( ) strange hyperons are experimentally favorable because their spin direction can be reconstructed via their weak decay topology into proton (antiproton) and charged pion. Recently, the STAR Collaboration observed nonzero global polarization of and hyperons in Au-Au collisions, first for the Relativistic Heavy Ion Collider (RHIC) beam energy scan (BES) energies of √ s NN = 7-39 GeV and, later, with more data, at √ s NN = 200 GeV [10,11]. The magnitude of the observed polarization varies from a few to a fraction of a percent. While the P H values for and agree within the experimental uncertainties, they are systematically higher for than for . Assuming that this difference originates due to the magnetic fields one estimates the field strength in units of elementary charge e to be eB ∼ 0.01m 2 π [12]. The exact nature of the spin-orbit interaction leading to the global polarization is not known. It is unclear at what stage of the system evolution (the QGP, hadronization, or the hadronic rescattering) the polarization is acquired, neither the corresponding relaxation times are known. Most of the recent calculations of the global polarization assume complete thermal equilibrium and validity of the hydrodynamical description of the system [6,[12][13][14][15]. They relate the particle polarization to the system's thermal vorticity at the hadronization time. In a nonrelativistic limit, assuming complete thermal equilibrium, the polarization of the particles can be evaluated as ζ = s /s = (s + 1)ω/(3T ), where s is the particle's spin, T is the system temperature, ω = (1/2)(∇ × v) is a nonrelativistic vorticity, and v is the local fluid velocity [12].
The global polarization is determined by the average vorticity component perpendicular to the collision reaction plane, which is spanned by the beam direction and the impact parameter vector. The global polarization measurements with the produced particle provide important information on both the nature of the spin-orbit interaction and the profile of velocity fields of the expanding system. Both the magnitude and the direction of the vorticity can strongly vary within the system [16]. In particular, a significant component along the beam direction can be acquired due to the transverse anisotropic flow [16,17].
The vorticity of the system, especially its component along the system's orbital momentum, is directly related to the asymmetries in the initial velocity fields. This links the vorticity with the directed flow v 1 , which is also strongly dependent on those asymmetries. The v 1 is defined by the first Fourier moment v 1 = cos(ϕ − RP ) of the produced particle's azimuthal asymmetry relative to the collision reaction plane angle RP . Hydrodynamic simulations show that the orbital angular momentum stored in the system and the directed flow of charged particles are almost directly proportional to each other [5]. This allows for an empirical estimate of the collision energy dependence of the global polarization [16]. The STAR results for the directed flow [18,19] and the hyperon global polarization [10,11] from the BES program show that the slopes of v 1 at midrapidity (dv 1 /dη) for charged hadrons (pions) and the hyperon polarization are indeed strongly correlated. The charged-particle directed flow in Pb-Pb collisions at √ s NN = 2.76 TeV is about three times smaller [20] than at the top RHIC energy of 200 GeV. This suggest that the global polarization at the Large Haron Collider (LHC) energies should be also about three times smaller than at RHIC (around ∼0.08%) and decreasing from √ s NN = 2.76 to 5.02 TeV by about ∼30% [21]. Even smaller polarization values at the LHC are expected when the directed flow is seen as a combination of the two effects: the tilt of the source in the longitudinal direction and the dipole flow originating from the asymmetry in the initial energy density distributions [22]. Taking into account that only the contribution to the directed flow from the tilted source is related to the vorticity and that its contribution relative to the dipole flow decreases with the collision energy [22], one arrives to an estimate for the global polarization at the LHC energies of the order of ∼0.04%. In this paper, the measurements of the global polarization of the and hyperons in Pb-Pb collisions at √ s NN = 2.76 and 5.02 TeV recorded with ALICE at the LHC are reported. The paper is organized as follows. In Sec. II the analysis details are presented, and the global polarization observable is introduced, as well as the measurement technique. The various sources of systematic uncertainties are discussed in Sec. III. The results for the and global polarization at two collision energies as a function of the hyperon transverse momentum and collision centrality are presented in Sec. IV.

A. The observable
In this measurement, and hyperons are reconstructed through their weak decay topologies → p + π − and → p + π + (64% branching ratio). The global polarization of the hyperons is determined from the angular distributions of their decay (anti-)protons. In the hyperon rest frame, the (anti-)proton angular probability distribution, dw/dn * p , is given by where n * p is the unit vector of the (anti-)proton direction. The same absolute value from [23] for the hyperon decay parameter α H = 0.642 ± 0.013, positive for and negative for , is used. Recent measurement by the BESIII Collaboration [24] extracted a new values for α = 0.750 ± 0.010 and α¯ = −0.758 ± 0.012, which are about 17% larger than in Ref. [23]. Given the statistical and systematical uncertainties of the results reported below for the hyperon global polarization the new values for α H will not change the conclusions of this paper but they should be considered for future hight precision measurements. The polarization vector ζ H in Eq. (1), which is satisfying condition |ζ H | 1, can be measured experimentally as [25] where the brackets . . . denote an event-by-event averaging over all hyperon decays. The polarization vector ζ H generally depends on the hyperon kinematics, namely the transverse momentum p T , rapidity y, and its azimuthal angle with respect to the reaction plane, ϕ − RP , as well as the collision centrality. The global polarization reported in this paper is determined by the component of the polarization vector perpendicular to the reaction plane. The magnitude P H of the global polarization can be measured by averaging a corresponding projection of the n * p vector, which in the laboratory coordinate system is given by n is the polar angle with respect to the collision axis (azimuthal angle) of the (anti-)proton direction in the hyperon rest frame. Substituting n * p⊥RP into Eq. (2) and assuming an ideal detector acceptance, an average over the θ * p yields Here the brackets . . . imply averaging over individual hyperons in all events. The polarization is defined to be positive if the hyperon spin has a positive component along the system's angular momentum-the same convention as employed in Ref. [10]. The detector acceptance effects are treated in this work as systematic uncertainty and discussed in Sec. III. A significant fraction of and hyperons originates from decay of heavier particles. Existing estimates [6,12,13] of the feed-down effect on the hyperon global polarization measurements, based on the assumption of thermal equilibrium and the particle yields from the statistical model [26], suggest that the primary and polarization should be by about 15-20% larger than what is measured at RHIC. Taking into account that feed-down effects are expected to be small compared to other uncertainties in the current analysis, no correction to the data has been applied.

B. Measurement technique
The main components of the ALICE detector system [27,28] used for this measurement are the Time Projection Chamber (TPC) [29], the silicon detectors of the Inner Tracking System (ITS) [30] and the two neutron Zero-Degree Calorimeters (ZDC) [31]. The analyzed data samples were During the data taking, the trigger required a hit in a pair of V0 detectors [32]. In 2010 and 2011, events with only one V0 hit and at least two hits in the outer layer of the Silicon Pixel Detector (SPD) were accepted as well. Contamination from beam-induced background was removed offline, as discussed in Refs. [33,34]. The events with poor correlation between multiplicities in V0, ITS, and TPC detectors were rejected. The analysis was restricted to the events with the primary vertex along the beam direction, V z , within ±10 cm from the nominal center of the TPC. This yielded for all collision centralities approximately 49 (75) million of Pb-Pb collisions at √ s NN = 2.76 (5.02) TeV. The collision centrality was determined using the energy deposition in the V0 detectors [35]. The reaction plane angle RP was estimated using the spectator plane angle SP , characterizing the deflection direction of the spectator neutrons. The spectator deflections at positive and negative rapidity were reconstructed with a pair of neutron ZDC detectors located 114 m away from the interaction point. The SP angle was evaluated separately for the two detectors using the transverse profile of the spectator energy distribution provided by the 2 × 2 segmentation of the ZDCs. For this, a pair of two-dimensional vectors Q t,p were constructed for two ZDCs following the procedure described in Ref. [20]: where indices p and t denote the ZDC on the projectile (η > 0) and target (η < 0) side of the interaction point, respectively; E i is the measured signal and n i = (x i , y i ) are the coordinates of the ith ZDC segment. The event averaged Q t,p revealed a strong dependence on the collision centrality as well as on the three collision vertex coordinates, which is imposed by the offset of the LHC beam transverse spot positions relative to the nominal center of the ZDCs. Moreover, the Q t,p demonstrated a strong time dependence in some of the 2015 runs. To compensate for these variations an event-by-event recentering correction [36] was applied as a function of collision centrality, three components of the collision vertex position, and beam-time variation, Additionally to the procedure described in Ref. [20], the width of the Q distributions was equalized as a function of centrality, which together with Eq. (5) resulted in an overall Q-vector correction, The width equalization, up to 20%, turned out to be particularly useful for the 2015 data sample, where one of the two neutron ZDC detectors lost signal from one of its four channels. SP for finite resolution of the spectator plane angle SP as a function of collision centrality for three data sets used in the analysis. Only statistical uncertainties, which are smaller than a symbol size, are shown.
The SP angle estimated from each ZDC is then given by the direction of the corrected Q vector, To account for a finite resolution of the spectator plane angle SP , a correction R (1) SP is introduce in the Eq. (3) following the method described in Ref. [25]: The polarization values obtained with Eq. (8) for each of the two ZDC detectors were found to be similar within the statistical uncertainties and were combined. The correction R (1) SP was extracted from correlations between Q-vector angles from different ALICE detectors following the technique described in Ref. [37]. Figure 1 presents the R (1) SP as a function of collision centrality for different data sets. The most central (0-5%) collisions are excluded from the analysis, because the small number of spectators does not allow for reliable estimation of the spectator deflection with the ZDCs. During the Pb-Pb data taking in 2011 ( √ s NN = 2.76 TeV) and 2015 ( √ s NN = 5.02 TeV), the beams were collided at a nonzero vertical crossing angle [28]. This resulted in a partial screening of the spectator neutrons by the LHC tertiary collimators [38] and in a degradation of the spectator plane resolution. The V0s, TPC, and two forward multiplicity detectors (FMD) [39] were used to estimate a possible uncertainties in R (1) SP extraction for the two ZDC detectors. In all data samples, these uncertainties turned out to be at a level of several percent.
The and hyperons were reconstructed via their weak decay topology following the method and selection criteria described in Refs. [28,40]. Charged daughter tracks from hyperon decay were required to have the pseudorapidity |η| < 0.9 and at least 70 space points in the TPC. The pion and (anti-)proton particle type assignments were based on the track charge and specific energy loss (dE/dx) measured in the TPC. The daughter tracks were paired to form and 044611-3 Bottom panels: Global polarization extraction via fit to sin(ϕ * p − SP ) as a function of the invariant mass, with SP reconstructed with the ZDC on the target (η < 0) side. Dashed and solid lines show L BG (M inv ) and the combined fit with Eq. (9). See text for more details.
candidates. The candidates were required to have transverse momentum p T > 0.5 GeV/c, rapidity |y| < 0.5, and the momentum vector pointing back to the primary collision vertex within a cone of opening angle less than 0.1 (0.08) radians for Pb-Pb collisions at √ s NN = 2.76 (5.02) TeV.
The hyperon global polarization P H was extracted with the fit to the measured correlation sin(ϕ * p − SP ) as a function of the or candidate invariant mass, M inv : (9) Here the constant S H gives the numerator of Eq. (8) and L(M inv ) is a linear parametrization of the background correlation as a function of M inv . The background fraction f BG (M inv ) was evaluated as Here N tot (M inv ) is the total measured yield of the or candidates. The N BG (M inv ) is given by the fourth-order polynomial fit to the N tot (M inv ) outside the and invariant mass peak range, M inv < 1.107 GeV/c 2 and M inv > 1.125 GeV/c 2 .
The fitting procedure is illustrated in Fig. 2 for 20-30% centrality range in Pb-Pb collisions at √ s NN = 2.76 TeV using data collected during the LHC operation in 2011.

III. SYSTEMATIC UNCERTAINTIES
The considered sources of systematic uncertainties are summarized in Table I. The significance of a given systematical variation is determined following the procedure presented in Refs. [41].
The selection on the primary vertex position V z was varied to ±7 cm and ±8 cm instead of the default ±10 cm. The analysis was repeated with the collision centrality estimated with either ITS or TPC detectors (instead of the default V0 detector). Variations of the results when changing the event selection criteria described above and centrality estimators are found to be negligible and were not included into the total systematic uncertainty.  The following and candidate selection criteria were varied: the distance of closest approach to the primary vertex, decay daughter track selection via specific energy loss in the TPC, the distance of closest approach between the pair of the daughter tracks, and the criteria on the hyperon candidates momentum pointing angle to the primary vertex. The corresponding contribution from these variations to the total systematic uncertainty is about 20-30% of the statistical un-certainty. The main contribution to the systematic uncertainty from the fitting procedure in Eq. (9) comes from a variation of the fit region.
In Eq. (8)  Results are compared with the STAR data at lower energies [10,11]. The insert shows zoomed-in comparison with the data at the top RHIC energy. The systematic uncertainties are shown as shaded boxes. Points are slightly shifted along the horizontal axis for better visibility.
distribution. Using the method described in Refs. [25,42] the corresponding relative uncertainty was estimated to be about 10% independent of centrality. This uncertainty comes primarily from the admixture of the higher-order harmonics when hyperon p T 2 GeV/c. A detailed study was performed for the evaluation of the systematic uncertainties of R (1) SP as well as for the evaluation of the difference between the two neutron ZDC detectors. The correlations between the flow vectors from different detectors, including TPC and pairs of ZDC, V0, FMD detectors were studied. The contribution to the total systematic uncertainties due to the R (1) SP extraction were found to be at a level of a few percentages.

Figures 3 and 4 show the measured hyperon global polar-
ization P H as a function of centrality and hyperon transverse momentum, p T , in Pb-Pb collisions for two collision energies. The results from 2010 and 2011 data samples were combined accounting for the corresponding statistical and systematic uncertainty. At RHIC energies, the global polarization ex-hibited a clear centrality dependence with three times large magnitude in peripheral collisions compared to that in central, while no significant p T dependence within the accessible p T range was observed [11]. The P H at the LHC is found to be consistent with zero within the experimental uncertainties for all studied centrality classes (Fig. 3) and p T ranges (Fig. 4). By repeating a similar analysis, no signal signal was observed as a function of rapidity either.
The average global polarization for two centrality ranges, 5-15% and 15-50%, are presented in Fig. 5, while numerical values are reported in Table II. Figure 5 also presents the comparison with the STAR data [10,11] for lower collision energies. Despite large uncertainties, the ALICE measurements confirm the trend of the global polarization decreasing with increasing collision energy.
Assuming the same values of the global polarization for and and neglecting the possible difference of about 30% (according to the empirical estimates discussed above) between the two LHC energies, one can average all four ALICE data points for 15-50% centrality, where the largest signal is expected. This yields a value P H (%) ≈ 0.01 ± 0.06 (stat.) ± 0.03 (syst.) for 15-50% centrality, which is consistent with the empirical estimates of P H (%) ≈ 0.04-0.08 based on directed flow measurements.

V. SUMMARY
The first measurements of and hyperons global polarization are reported for Pb-Pb collisions at √ s NN = 2.76 and 5.02 TeV recorded with ALICE at the LHC. The hyperon global polarization has been measured differentially as a function of centrality and transverse momentum (p T ) for the range of collision centrality 5-50%, 0.5 < p T < 5 GeV/c, and rapidity |y| < 0.5. No significant dependences neither splitting between the global polarization values for and has been observed. The average and polarization for 15-50% centrality range at two collision energies is found to be consistent with zero, P H (%) ≈ 0.01 ± 0.06 (stat.) ± 0.03 (syst.). This confirms the observed earlier trend of the global polarization decrease with increasing collision energy. The results are compatible with the extrapolation of the RHIC results and empirical estimates of P H (%) ≈ 0.04-0.08 based on the similarity of collision energy dependence of the global polarization and the slope of the directed flow in the midrapidity region. The high luminosity LHC run after the 2019-2021 long shutdown with the upgraded ALICE detector will bring more than 100 times more data, which should allow tests of the previously mentioned prediction with much better accuracy.

ACKNOWLEDGMENTS
The ALICE Collaboration thanks all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centers and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A. I. Alikhanyan National Science Laboratory (