Evolution of global polarization in relativistic heavy-ion collisions within a perturbative approach

Extremely large angular orbital momentum can be produced in non-central heavy-ion collisions, leading to a strong transverse polarization of partons that scatter through the quark-gluon plasma (QGP) due to spin-orbital coupling. We develop a perturbative approach to describe the formation and spacetime evolution of quark polarization inside the QGP. Polarization from both the initial hard scatterings and interactions with the QGP have been consistently described using the quark-potential scattering approach, which has been coupled to realistic initial condition calculation and the subsequent (3+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(3+1)$$\end{document}-dimensional viscous hydrodynamic simulation of the QGP for the first time. Within this improved approach, we have found that different spacetime-rapidity-dependent initial energy density distributions generate different time evolution profiles of the longitudinal flow velocity gradient of the QGP, which further lead to an approximately 15% difference in the final polarization of quarks collected on the hadronization hypersurface of the QGP. Therefore, in addition to the collective flow coefficients, the hyperon polarization may serve as a novel tool to help constrain the initial condition of the hot nuclear matter created in high-energy nuclear collisions.


I. INTRODUCTION
High-energy nucleus-nuclues collisions at the Relativistic Heavy-Ion Collider (RHIC) and the Large Hadron Collider (LHC) create a color deconfined state of nuclear matter, known as the Quark-Gluon Plasma (QGP), whose strongly-coupled nature has been confirmed by both the collective flow of soft hadrons and the quenching phenomena of jets observed in these energetic collisions [1,2].In addition to the hottest and densest environment one may obtain in laboratory, non-central heavy-ion collisions also deposit huge amount of angular momentum (on the order of 10 5 ) into the nuclear matter [3], leading to the most vortical system (vorticity on the order of 10 21 s −1 ) one may create [4,5].
It was proposed in Ref. [6] that this large angular momentum can generate polarization of quarks via their spin-orbital interactions, and in the end be observed as the global polarization of hyperons, e.g.Λ, Σ and Ξ.Other possible observables include vector meson spin alignment [7] and the emission of circularly polarized photons [8].On the experimental side, search for the Λ polarization in heavy-ion collisions was initiated by the STAR Collaboration [9] and recently confirmed in Refs.[10,11].The non-trivial dependences of the Λ polarization on the collision energy, rapidity region and the transverse momentum of Λ have attracted tremendous theoretical efforts on exploring the detailed mechanisms that generate polarization [12].
A major category of theoretical approaches is based on the assumption that the spin degrees of freedom are at local equilibrium on the hadronization hypersurface of the QGP [13,14], allowing one to extract the polarization of the final-state hadrons from the extended Cooper-Frye formalism for particles with spin.The key quantity that drives the development of polarization is the thermal vorticity [13,15].The evolution of the QGP can be either simulated with a perturbative-based transport model (e.g.APMT) [16][17][18] or be considered as a strongly-coupled medium (e.g.hydrodynamics) [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33].While many of these studies provide a reasonable description of the experimental data, in order to address non-equilibrium effects on polarization, it is necessary to introduce the quantum kinetic theory.This have been recently developed in Refs.[34][35][36][37][38], where crucial questions like conservation of angular momentum in a transport model and evolution towards local equilibrium have been explored.Instead of a full transport description of interactions between quark pairs, an alternative microscopic approach was proposed in Ref. [39], where each quark is considered interacting with a mean field of the medium background.The perturbative scattering picture from the earliest study [6] has been extended to multiple scatterings of quarks through the QGP, from which one can observe the development of the quark polarization driven by the longitudinal velocity gradient of the fluid background.
In the present study, we will further develop this perturbative approach [39] in several aspects.We will consistently apply this perturbative scattering picture to both the initial hard collisions between nuclei and quark scatterings through the QGP for the first time, which allows us to study the initial production of polarization and its subsequent evolution within the same framework.A simplification was usually applied in earlier perturbative calculations [6,39] where a projectile quark is constrained in a half hemisphere relative to the target potential.Realistic spatial distributions of projectile and target will be taken into account for a more precise evaluation of the average quark polarization in this work.In addition, a (3+1)-D viscous hydrodynamic model CLVisc [40,41] will be introduced to provide a realistic evolution of the longitudinal flow velocity profile of the QGP.Impacts of using different initial conditions for the hydrodynamic expansion on the final-state global polarization will be investigated in detail.We will focus on the global polarization that results from the large angular momentum with a direction perpendicular to the reaction plane of noncentral heavy-ion collisions.The other crucial branch of polarization, the local (or longitudinal) polarization due to the QGP expansion in the transverse plane [42][43][44] is beyond the scope of this work.
The rest of this paper will be organized as follows.In Sec.II, we will discuss the quark polarization produced by the initial hard collisions, and investigate how it is affected by the nucleon density distribution.In Sec.III, we will study the evolution of the global polarization through the QGP phase, and explore how the final-state polarization depends on the initial energy density profile of the hydrodynamic evolution of the QGP.A summary will be presented in Sec.IV.

II. INITIAL POLARIZATION FROM HARD SCATTERINGS
We consider a non-central collision event between two nuclei, as demonstrated in Fig. 1.The momentum ( p) of each nucleon from the projectile nucleus (B) is assigned along the +ẑ direction, while the momentum of each nucleon from the target (A) is the opposite.The impact parameter ( b), pointing from the target to the projectile, is assigned along the +x direction.The reaction plane is then defined as the z-x plane, whose direction can be determined as n = +ŷ.This non-central collision would deposit huge angular momentum into the overlapping region between the two nuclei, leading to the polarization of quarks that constitute the QGP medium via the spinorbital coupling.In this section, we study the formation of the polarization from this initial hard nucleus-nucleus scattering, and discuss how the magnitude of polarization depends on the initial nucleus geometry.
Following Ref. [6], we investigate a scattering between a quark and a static potential A 0 (q T ) = g/(q 2 T + µ 2 ), where g 2 = α s /(4π) is the strong coupling constant, q T represents the momentum transfer and µ is Debye screening mass of the exchanged gluon.This potential scattering picture can be consistently applied to the initial hard scattering and the subsequent quark scattering through the QGP medium.For the initial hard scattering in vacuum, we take µ = 0.5 GeV as in Ref. [6]; while for the in-medium scattering, µ 2 = g 2 (N c + N f /2)T 2 /3 will be adopted in the next section, with T being the local temperature of the QGP and N c = N f = 3 in the present study.For a given initial momentum (E, p) and a final spin λ/2 of the outgoing quark along n, the differential cross section for an initially unpolarized quark can be obtained as [6] with In the above equations, C T = 2/9 is the color factor associated with the target, A µ = (A 0 , 0) is the scattering potential, q T ( k T ) is the momentum transfer between the quark and the potential in the real (complex conjugate) space, x T is the relative transverse distance pointing from the scattering center to the quark (as illustrated in Fig. 2), and m is the quark mass.After being scattered, the outgoing quark momentum is p q(k) = p + q T ( k T ).
Here, we have averaged over the initial spin states (λ i ) for the unpolarized projectile quark.For a high energy quark, one may assume a small scattering angle (q T , k T ∼ µ E), which simplifies the I λ part as The first part in the square bracket is spin-independent, which contributes to the total cross section as while the second term is spin-dependent, contributing to the difference between positive and negative spin configurations as Here, xT = x T /x T and K n 's are the modified Bessel functions.The polarization of quarks is then defined as the ratio between the integrated ∆σ and σ above, P = ∆σ/σ.Although this ratio does not explicitly depend on α s (or g), the strong coupling constant will affect the magnitude of polarization through the Debye screening mass, as will be discussed later in this work.When using the above equation to evaluate the initial polarization, we take the center-of-mass momentum p = 100 GeV along ±ẑ directions for Au-Au collisions at √ s NN = 200 GeV, and the quark mass m = 100 MeV for estimating their energy E. We have verified that varying the quark mass within 5 ∼ 300 MeV has no visible impact on the polarization result.Observed from Eq. ( 5), the sign of the cross section difference depends on the relative position between the projectile quark and the scattering center.As sketched in Fig. 2, if the projectile quark (e.g.quark 1) is in the x > 0 hemisphere with respective to the scattering center, d∆σ/d 2 x T is negative; contrarily, for the x < 0 hemisphere (e.g.quark 2), d∆σ/d 2 x T is positive.We note that in earlier studies [6,39], a simplified assumption is adopted for two-particle scattering, where one quark from B is always in the x > 0 hemisphere relative to the other quark from A. By integrating over the corresponding half hemisphere, a heuristic approximation for the quark polarization is obtained as 1 To improve the estimation of polarization from non-central nuclear collisions, we extend this two-particle scattering picture with relative position restricted in the half hemisphere to scatterings between all pairs of nucleons that are realistically distributed in the two colliding nuclei.We will show that the finite polarization in the end is not because one quark from B is always above (x > 0) the other quark from A, but because the nuclear matter contributed by B is on average above that by A.
To take into account the configurations of nuclei A and B, we need to integrate Eqs. ( 4) and ( 5) over the whole transverse plane, weighted by the thickness functions of A and B. Two different models of the nuclear matter distribution are used and compared here.The first one is the Hard Sphere Model, where the nuclear matter is assumed to be strictly constrained within a radius R. The corresponding thickness function reads where the centers of nuclei A and B are placed at (−b/2, 0) and (b/2, 0) in the transverse plane respectively, and A on the right hand side denotes the nucleon number (same for A and B in this work).Another model for the nuclear matter density is the Woods-Saxon distribution, written as where ρ 0 is the equilibrium density of nuclear matter, and a is the surface thickness parameter.Unlike the Hard Sphere Model, the nuclear matter density in the Woods-Saxon model decreases smoothly, though rapidly, when the position is over R + a away from the nucleus center.Since only nucleons that participate in the initial nucleus-nucleus collisions contribute to the QGP medium, we only include those participant nucleons in our calculation.For the above two models, the participant number density read [3] respectively, in which σ NN is the inelastic cross section of nucleon-nucleon collision, and the subscript 1(2) is used to represent participants from nucleus A(B).For Au-Au collisions at √ s NN = 200 GeV, parameters in Eqs. ( 6) -( 9) are set as A = 197, R = 6.38 fm, ρ 0 = 0.17 fm −3 , a = 0.535 fm and σ NN = 42 mb [45].With these participant number density distributions, the final-spin-summed probability for a quark from nucleus A at a given location (x A , y A ) to scatter with the entire nucleus B can be written as 7.5 5.0 2.5 0.0 2.5 5.0 7.5 x>IP@  where an integral weighted by the participant density in B is implemented.The same definition can also be applied for the spin difference in the scattering probability ∆P A (x A , y A ) from its corresponding cross section d∆σ.Similarly, one can also obtain the spin-summed and spindifferent probability -P B (x B , y B ) and ∆P B (x B , y B ) -for a quark in nucleus B at (x B , y B ) to scatter with the entire nucleus A. By combining contributions from A and B, we obtain the following polarization of nuclear matter at a given location: , y) .
(11) In Fig. 3, we first present the spatial distribution of the initial polarization calculated via Eq.( 11) above.A Au-Au collision system at √ s NN = 200 GeV is used here, with the impact parameter set as b = 6.7 fm, approximating the 5-40% centrality region.Comparing between the two subfigures, one can observe the dependence of polarization on the nuclear distribution function.In the upper panel, we find large polarization values mainly distribute around the boundary of the overlapping region between the two nuclei for the Hard Sphere distribution.This could be understood with the illustration in Fig. 4. For location 1 that is at the center of the overlapping region, a quark from A interacts with similar number of participants from B on its left (x < 0) and right (x > 0), cancelling polarization due to our earlier discussion for Eq. ( 5).On the other hand, for location 2 that resides on the boundary, a quark from A "sees" most participants from B on its right, contributing to a negative polarization.Although the opposite conclusion (positive polarization) is drawn for a quark from B at location 2, its density [T 2 in Eq. ( 11)] is much smaller than that from A [T 1 ] at this location, leading to a net negative value of polarization in the end.Since there is no interaction outside the overlapping region within the Hard Sphere model, polarization is also zero in the corresponding region.Different distribution of polarization can be observed in the lower panel for the Woods-Saxon model.While the polarization is also small at the center and maximized around the boundary of the overlapping region, its magnitude decreases smoothly outside the boundary instead of suddenly vanishes because of the different location (x, y) dependences of the nuclear density between Eq. ( 6) and Eq.(7).It is interesting to note that for the case of Woods-Saxon model, one can also observe positive polarization value far away from the center.As illustrated in Fig. 4, at location 3, a quark from A is polarized along −ŷ by interacting with nucleus B; while a quark from B is polarized along ŷ by interacting with A. Since the scattering cross section Eqs. ( 4) and ( 5) rapidly decreases with distance, the magnitude of the above polarization for the quark in A [∆P A in Eq. ( 11) from Eq. ( 10)] is much smaller than that of the quark in B [∆P B ], leading to a net positive polarization in the end.However, since the participant number densities [both T 1 and T 2 in Eq. (11)] are small at location 3, this positive value has little contribution to the global polarization after we integrate over the entire transverse plane, as will be shown later in this work.
By integrating over the transverse plane, we may define the global polarization as follows, b>IP@ P>@ Hard Sphere Woods-Saxon The corresponding value is presented in Fig. 5 as a function of the impact parameter for Au-Au collisions at √ s NN = 200 GeV.As expected, the global polarization is along the −ŷ direction, reflected by its negative value, because on average participants from B is on the right (x > 0) of participants from A. As the impact parameter increases, the magnitude of polarization first increases due to more asymmetric collisions, but then decreases due to vanishing participant nucleons in the initial hard scatterings.The maximum magnitude is obtained when b is around 2R for sideswipe between the two nuclei.Beyond that, the polarization suddenly disappears for the Hard Sphere model due to the hard cut-off of its nuclear overlap function, while smoothly decreases to zero for the Woods-Saxon model.Results from these two nuclear density distributions are almost identical to each other for small impact parameter (b R); noticeable difference is only observed at larger b.In this work, we initialize polarization according to smooth nuclear density distributions.Effect from event-by-event fluctuations on the global polarization was shown small in Ref. [46].On the other hand, the local polarization along the longitudinal direction due to the radial flow of the QGP might be more sensitive to the initial state fluctuation.This will be explored in our future work by implementing our current framework using the Monte-Carlo method.

III. EVOLUTION OF POLARIZATION THROUGH THE QGP PHASE A. Polarization of initially polarized quarks
In the previous section, we discussed the production of global polarization from initial hard collisions between nuclei, where Eqs. ( 1)-( 5) were derived for initially unpolarized quarks.To investigate the further evolution of polarization through the QGP phase, one needs to ex-tend these equations to the scenario where the projectile quark already possesses non-zero polarization.This has been developed in Ref. [39], where Eq. ( 1) is re-written as in which R λi = (1 + λ i P i )/2, with P i being the initial polarization of the projectile quark.Using the small scattering angle approximation as before, the I λ f λi part can be approximated with which returns to Eq. ( 3) after the initial spin states (λ i ) are averaged over.The differential cross section is then simplified to The final-spin-independent sum of the cross section reads while the final-spin-dependent difference reads Instead of integrating over the half x-y plane for scattering centers always on one side of the projectile quark, as assumed in Ref. [39], we obtain the total cross section by integrating over the entire transverse plane, weighted by the local entropy density (or particle number density) of the QGP.Therefore, the final polarization of a quark residing at a given location (x, y) is given by where the local entropy density s can be taken from a hydrodynamic simulation of the QGP medium, as will be discussed in the next subsection.
The change of polarization after the scattering is then obtained as For multiple scatterings inside the QGP medium, the relative momentum p in Eq. ( 19) is taken as Ev z for estimating the polarization along the ŷ direction, where E = /ρ represents the energy of a quark with and ρ being the local energy and particle number densities inside the QGP, and v z is the local fluid velocity along the ẑ direction.In addition, we assume the mean free path of a quark is τ q , which can be roughly related to the shear viscosity of the QGP via η v ≈ (1/3)ρ p th (4/9)τ q ≈ (4/9)T ρτ q for a thermal ensemble of gluons [39,47], with p th = 3T being the average thermal momentum at temperature T .This leads to the following equation for the time evolution of polarization: We will use this equation to evolve the quark polarization at each location (x, y) inside the QGP.The local temperature (T ), number density (ρ) and entropy density (s) can be provided by the hydrodynamic model, and the shear-viscosity-to-entropy-density ratio is taken as η v /s = 0.08.

B. Hydrodynamic simulation of the QGP
In this work, we use the (3+1)-dimensional viscous hydrodynamic model CLVisc [40,41] to simulate the spacetime evolution of the QGP medium.Following our earlier work [48], three different model setups are used to generate the initial energy density distribution of the medium, which is fed as the initial condition into the hydrodynamic evolution.By comparing between these three model setups, one may explore how the final-state global polarization depends on the initial geometry of the medium.
Our first parametrization of the initial condition is based on Refs.[49,50], which is one of the earliest and most widely applied initialization scheme that takes into account the deformation of nuclear matter in the reaction plane due to non-central collisions, as illustrated in Fig. 6.In this model, the transverse (x, y) and longitudinal (spacetime rapidity η s ) distribution of the energy density is parametrized with a wounded nucleon weight function as follows, (21) in which T 2,1 is the participant nucleon density distributions of the ±ẑ-going nucleus, as discussed in Eq. ( 9).The Woods-Saxon model of the nucleon distribution is used for initializing the hydrodynamic evolution.The f ± function is designed to introduce the geometric asymmetry along the longitudinal direction as where η m defines the range of rapidity correlations, which affects the relative contribution from forward and backward participating nucleons.This parametrization introduces not only a longitudinal shift but also a counter-clockwise tilt in the x-η s plane to the 3-dimensional geometry of the medium, as illustrated in the left panel of Fig. 6.The full expression of the initial energy density distribution is then given by (x, y, η s ) = K • W (x, y, η s ) • H(η s ), (22) where K is an overall normalization factor that is fitted to the final charged particle yield (dN ch /dη with η being the pseudorapidity) observed in relativistic heavy-ion collisions.A function H(η s ) is introduced to describe the plateau pattern of dN ch /dη with respect to η, with η w and σ η as two model parameters.The total weight function combines contributions from wounded nucleons and binary collisions -the latter is given by The parameter α determines the relative contribution from wounded nucleons and binary collisions, which can be extracted from the impact parameter dependence of dN ch /dη.Related model parameters will be listed later in this subsection when we calculate for a given collision system.
Case (B) CCNU parametrization.An alternative parameterization of the deformed initial energy density was developed in Ref. [51], which is similar to the above Bozėk-Wyskiel ansatz -Eqs.( 22) - (25), except that the longitudinal dependence of the wounded nucleon weight function is parametrized as Two parameters -H t and η t -are introduced to describe the unbalanced energy deposition, between the projectile and target nuclei at different transverse locations, into the medium at different spacetime rapidities.Similar to the Bozėk-Wyskiel setup, this CCNU parameterization will also cause both shift and tilt of the initial energy density distribution as illustrated in the left panel of Fig. 6.Case (C) Shen-Alzhrani parametrization.The third η s -dependent initial condition model was adopted from Refs.[32,52], which ensures the local energy-momentum conservation when convert the two colliding nuclei into the energy density profile of the hot nuclear medium.One first defines the local invariant mass M (x, y) and the center-of-mass rapidity y CM as, where y beam = arccosh( √ s NN /2m N ) is the rapidity of each nucleon inside the colliding beams, with m N being its mass.The initial energy density profile is then constructed as [52], Same as the previous two models, K is the overall normalization factor, and η w and σ η are the width parameters for the plateau width of dN ch /dη distribution with respect to η.A new parameter -y L = f y CM with f ∈ [0, 1] -is introduced in this model to describe the deformation of the medium along the longitudinal direction.
The transverse density distribution N e is determined by the local invariant mass M (x, y) as where erfc(x) is the complementary error function.Different from Bozėk-Wyskiel and CCNU parameterizations, this Shen-Alzhrani parametrization only generates a shift deformation along the longitudinal direction, as illustrated in the right panel of Fig. 6; tilt of the medium profile has not been included yet.A detailed comparison of this deformation can be found in our earlier work [48].In many hydrodynamic calculations, including our previous study [48], local flow velocities are initialized with zero.While this simplification has minor effect on observables that are mainly driven by the QGP expansion, e.g. the yield of charged particles and their harmonic flow coefficients, it ignores the initial orbital angular momentum deposited into the system and thus would fail in describing observables related to the global polarization.As revealed in Ref. [39], the velocity gradient ∂v z /∂x is the main origin of the global polarization generated inside the QGP.Therefore, for all cases (A, B and C above), we follow Ref. [52] to initialize the off-diagonal components of the energy-momentum tensor as The initial flow velocity in the longitudinal direction is then given by their ratio as From the above equations, one notices that the y L (or f ) parameter designed in Case (C) determines the amount of the longitudinal momentum from the beam nucleons that is deposited into the QGP medium as its initial longitudinal velocity.In the present study, we only introduce the non-zero initialization of the longitudinal velocity, the transverse components of the energy-momentum tensor (T τ x and T τ y ), or the corresponding flow velocities (v x and v y ) are still initialized as zero, since they are not expected to affect the global polarization which aligns with the ŷ direction.Based on the above initial energy density and flow velocity, we then apply the CLVisc hydrodynamic model [40,41] to simulate the subsequent spacetime evolution of the QGP profiles, starting from an initial proper time τ 0 .The hydrodynamic equations read: where the energy-momentum tensor is given by Here, is the local energy density, u µ is the fluid fourvelocity, P is the pressure, π µν is the shear stress tensor, and Π is the bulk pressure (taken as zero in our current calculations); is the metric tensor and ∆ µν = g µν − u µ u ν is the projection tensor.These hydrodynamic equations are solved together with the lattice QCD Equation of State (EoS) from the Wuppertal-Budapest work (2014) [53].
In this study, we apply the isothermal freeze-out conditions [40], in which the freeze-out hypersurface is determined by a constant temperature T frz .Interaction between a projectile quark and its surrounding medium, and thus the evolution of its polarization, ceases when the local temperature drops below T frz .
In the rest of this paper, we will use the 5-40% Au-Au collisions at √ s NN = 200 GeV, if not otherwise specified, as an example to study the evolution of quark polarization inside the QGP.In Tab.I, we summarize all our model parameters mentioned earlier in this subsection, which were constrained in our previous study [51] by the soft hadron yield and their directed flow coefficient.As previously discussed for Eq. ( 19), the evolution of quark polarization inside the QGP depends on the t>IP@ < T > >*H9@ Bozek-Wyskiel CCNU Shen-Alzhrani medium temperature (T ) and longitudinal velocity (v z ) profiles.Therefore, we first investigate how different initial conditions affect these quantities during the hydrodynamic expansion.In Fig. 7, we study the time evolution of the average temperature of the QGP fireball.The average is conducted on the η s = 0 (or z = 0) plane over the hydrodynamic cells with local temperature above T frz .The local entropy density s is applied as the weight for the average.As shown in Fig. 7, the average temperature decreases with time towards T frz .No visible difference can be observed for the temperature evolution between the three initial condition setups, because they are all adjusted to describe the soft hadron yield data in our earlier work [51].
Shown in Fig. 8 is the v z distribution along the x-axis (y = z = 0) at different times.Due to the asymmetric collision, one expects to see v z > 0 in the x > 0 half plane while v z < 0 in the x < 0 half plane at the initial time.As time evolves, the magnitude of v z becomes smaller.It is interesting to note that while the magnitude of v z starting from the Shen-Alzhrani initial condition (bottom panel) decreases with time towards zero, the sign of v z from the other two initial conditions (top and middle panels) can flip at later time (e.g.t = 6 fm), leading to v z < 0 at x > 0 and v z > 0 at x < 0. This could be understood with the tilted geometry in the Bozėk-Wyskiel and CCNU initial conditions as illustrated in the left panel of Fig. 6, whose further expansion produces negative v z component in the x > 0 plane, while positive v z component in the x < 0 plane.This drives a quicker decay of v z from these two initial conditions (compared to the Shen-Alzhrani initial condition without such tilt) and in the end can also reverse the sign of v z .Since the velocity gradient is the key origin of the development of quark polarization inside the QGP, we expect to obtain different magnitudes of the final-state global polarization from these initial conditions, as will be presented in the coming subsection.We also note that the magnitude of v z obtained from the realistic hydrodynamic simulation here is much smaller than the relativistic laminar flow model applied in the earlier study [39], which will affect the magnitude of global polarization for the final state.

C. Evolution of the global polarization
Using the temperature and longitudinal velocity profiles provided by the hydrodynamic simulation, we are able to calculate the evolution of the quark polarization at a given position via Eq.(20).
Shown in Fig. 9 are the snapshots of the spatial distribution of the quark polarization at different times.The strong coupling constant is taken as g = 2 (α s = 0.3) for these calculations.The upper four panels evolve from the initial polarization evaluated with the Hard Sphere nucleon density distribution (the upper panel of Fig. 3), while the lower four panels are from the Woods-Saxon distribution (the lower panel of Fig. 3).The hydrodynamic expansion is initialized with the CCNU model.From the figure, one can see that as time evolves, the quark polarization generally increases in magnitude (or becomes more negative).Compared between the two nucleon density functions, we find a smoother distribution for the quark polarization across the transverse plane from the Woods-Saxon than the Hard Sphere distribution.It is interesting to note that positive values of polarization exist at locations far away from the QGP center not only for the Woods-Saxon distribution, but also for the Hard Sphere distribution.This is driven by the negative ∂v z /∂x gradient at large |x| at the initial time (as shown in Fig. 8).As time evolves, these positive values can disappear at locations far away from the QGP center, but start to appear near the center, which can be understood with the sign flip of ∂v z /∂x both at large |x| and around x = 0 during the hydrodynamic expansion that starts with the CCNU initial condition (middle panel of Fig. 8).At each time step, the polarization of each quark is only affected by fluid cells with local temperature above the freeze-out temperature T frz here.We have verified that no visible difference can be observed between the CCNU and the Bozėk-Wyskiel initial condition.In Fig. 10, we present the similar snapshots of quark polarization to Fig. 9, except that the hydrodynamic simulation starts with the initial energy density from the Shen-Alzhrani model.Compared between Fig. 10 and Fig. 9, we observe a quicker increase (in magnitude) of the quark polarization within the QGP regime, which could be understood with the slower decay of the longitudinal flow velocity gradient with the Shen-Alzhrani initialization than the CCNU initialization, as we previously discussed in Fig. 8.No inversion is observed for either the positive value of polarization away from the center or the negative value near the center here, because the sign of ∂v z /∂x remains during the hydrodynamic evolution if it is initialized with the Shen-Alzhrani model.
In the end, we calculate the average polarization over the entire transverse plane.We assume that when the local temperature drops below T frz , the fluid cell hadronizes and stops participating in subsequent scatterings.Therefore, the corresponding polarization also freezes.For each time step of the hydrodynamic evolution, we first calculate the change of polarization at each location as ∆P (x, y, t).The change of the average polarization within this time step is then contributed by fluid cells above T frz and is given by ∆P (t) = T >T frz dxdys(x, y, t)∆P (x, y, t) dxdys(x, y, t) .
This ∆P (t) is then applied to calculate the time evolution of the transverse-plane-averaged polarization as which starts from an initial value produced by the primordial hard scatterings (shown in Fig. 5).
In Fig. 11 we present the time evolution of the transverse-plane-averaged polarization at mid-spacetimerapidity.Although positive values of polarization can exist at certain spacetime, due to the relatively small entropy (energy) densities at these locations, the average polarization is still negative as expected.In the upper panel, we use g = 2 as for Figs. 9 and 10, and compare between hydrodynamic expansion with different initial conditions.One observes that the fluid velocity profiles of the QGP have a non-negligible impact on the final average polarization.The tilted geometry of the initial energy density, as implemented in the Bozėk-Wyskiel and CCNU models, leads to a slower increase (in magnitude) of polarization compared to that from the Shen-Alzhrani initialization.This can be understood with the faster decay of the longitudinal velocity gradient in the former case than the latter, as discussed earlier.The dependence of the global polarization on the initial condition of the QGP has also been discussed in Refs.[11,54].We also notice that the global polarization mainly develops during the first 2 fm of the QGP evolution due to the large v z gradient and energy density at the early time.The global polarization we obtain for the final state is around 0.095%∼0.109%(about 15% uncertainty), depending on the selected initial condition.This is qualitatively consistent with (though quantitatively smaller than) the Λ polarization -0.277 ± 0.040 (stat) ± 0.039 0.049 (sys) [%]observed at the top RHIC energy [11].The magnitude of polarization we obtain here using a (3+1)-D relativistic hydrodynamic model is much smaller than that in the earlier study [39] using the relativistic laminar flow model, because of the weaker longitudinal flow velocity from our simulation.
In the middle panel of Fig. 11, we investigate effects of the initial nucleon density distribution and the strong coupling constant on the final-state global polarization.Although difference between the Hard Sphere and Woods-Saxon distributions could be observed in the t>IP@ P>@ spatial distribution of polarization previously in Figs. 9 and 10, it can be hardly seen here after the polarization has been averaged over the transverse plane.As discussed in Sec.II, the strong coupling constant can affect the quark polarization through the Debye screening mass.Since µ appears on both the numerator and the denominator of Eq. ( 19), and the value of polarization can be either positive or negative across the transverse plane, the g-dependence of the quark polarization differs from location to location.After averaging over the entire transverse plane, we find a larger magnitude of the quark polarization when g is reduced from 2 to 1.
In the lower panel of Fig. 11, we study the impact parameter dependence of this time evolution of polarization.As the impact parameter increases, a significantly larger magnitude of global polarization is obtained.This is consistent with the stronger orbital angular momentum deposited into the nuclear matter in more peripheral collisions.Note that as b increases, not only the initial polarization from the primordial hard scatterings becomes larger, as was also shown in Fig. 5, the increase of polarization during the QGP evolution becomes stronger as well because of the larger fluid velocity gradient ∂v z /∂x formed in more peripheral collisions.

IV. SUMMARY AND OUTLOOK
We have investigated the production and evolution of the global polarization of quarks in relativistic heavyion collisions within a perturbative approach.The spinindependent and dependent parts of the quark-potential scattering cross sections have been consistently applied to both the initial hard scatterings between colliding nuclei and the subsequent quark scatterings through the QGP.Compared to earlier studies, we have improved this perturbative approach in two aspects: (1) the two-body scattering model where the projectile quark is constrained in a half hemisphere relative to the target potential has been extended to realistic spatial distributions of targets and projectiles; and (2) a (3+1)-D viscous hydrodynamic model is adopted for simulating the QGP expansion.Effects of the nucleon density function and the initial geometry of the QGP medium on the final-state quark polarization have been explored in detail.
Within this improved approach, we have found that the spatial distribution of quark polarization depends on the nucleon density function inside the colliding nuclei.While polarization sharply centers around the edge of the overlapping region between the colliding nuclei with the Hard Sphere model, a smoother distribution across the transverse plane with possible positive values far away from the overlapping region can be seen with the Woods-Saxon model.However, after averaging over the transverse plane, these two density functions provide consistent magnitudes of the quark polarization, except for very peripheral collisions.Three different setups of the initial energy density distributions of the QGP have been compared in this work.It has been found that with a counter-clockwise tilted initial geometry in the reaction plane, the QGP expansion leads to opposite longitudinal flow velocity (v z ) to its initial direction.This can accelerate the decay of v z and even reverse its direction at late time, thus resulting in a smaller magnitude of quark polarization compared to calculation without using the tilted initial condition.The transverse-plane-averaged global polarization we obtain is about 0.095%∼0.109%at mid-spacetime-rapidity when the strong coupling constant is taken as g = 2 inside the QGP.The approximately 15% uncertainty quantifies the sensitivity of the global polarization to the longitudinal flow velocity profiles of the QGP, and may serve as a novel tool to help constrain the initial energy density distribution of the QGP in the future, when both theoretical calculations and experimental measurements become more precise.
While our study constitutes a step forward in a more quantitative understanding of the production and evolution of the global polarization using perturbative calculations, it should be further extended in several directions.For instance, our current calculation is limited at the quark level.A sophisticated hadronization scheme [34,55] and decay contributions to polarized resonant states [22] should be introduced to connect our current result to the realistic polarization of Λ hyperons measured by experiments.In addition, it is also important to study the beam energy and rapidity dependences of the global polarization, as have already been measured by the STAR experiments [10,11].To achieve this, one may also need to release the small angle approximation applied in our current calculation [Eqs.( 3) and ( 14)], which might not be valid when the center-of-mass energy of a colliding quark pair is small [3].We will address these aspects in our follow-up efforts.

FIG 1 .
FIG 1. (Color online) Side view (left panel) and top view (right panel) on a non-central heavy-ion collision event.

FIG. 3 .
FIG. 3. (Color online) Initial spatial distribution of polarization in Au-Au collisions at √ sNN = 200 GeV with impact parameter b = 6.7 fm, compared between Hard Sphere (upper panel) and Woods-Saxon (lower panel) distributions of nuclear density.

FIG. 5 .
FIG. 5. (Color online) The initial global polarization as a function of the impact parameter, compared between Hard Sphere and Woods-Saxon nuclear density distributions.

FIG. 6 .
FIG. 6. (Color online) Illustration of the initial energy density distribution on the x-ηs plane, left for our Case (A) and Case (B) that include a counter-clockwise tilt, and right for our Case (C) that only includes a longitudinal shift.

FIG 7 .
FIG 7. (Color online) Time evolution of the average temperature T , compared between three different initial condition setups.

FIG. 8 .
FIG. 8. (Color online)The longitudinal velocity distribution along the x-axis at different times.Results from three different initial condition setups are presented in the three panels.

FIG 9 .
FIG 9. (Color online) Spatial distribution of quark polarization in the transverse plane at different times, compared between Hard Sphere (HS) and Woods-Saxon (WS) distributions of the initial nucleon density.The hydrodynamic evolution is initialized with the CCNU model.

FIG 10 .
FIG 10. (Color online) Spatial distribution of quark polarization in the transverse plane at different times, compared between Hard Sphere (HS) and Woods-Saxon (WS) distributions of the initial nucleon density.The hydrodynamic evolution is initialized with the Shen-Alzhrani model.

FIG. 11 .
FIG. 11. (Color online) Time evolution of the transverseplane-averaged polarization, compared between using different initial conditions for the hydrodynamic evolution (upper panel), different density functions of the initial nucleon distribution and different values of g (middle panel), and different impact parameters (lower panel).
Ht and ηt are for Case (B), and the other parameters are commonly applied to all models here.