Local spin polarization in high energy heavy ion collisions

We revisit the azimuthal angle dependence of the local spin polarization of hyperons in heavy-ion collisions at 200 GeV in the framework of the (3+1)D viscous hydrodynamic model CLVisc. Two different initial conditions are considered in our simulation: the optical Glauber initial condition without initial orbital angular momentum and the AMPT initial condition with an initial orbital angular momentum. We find that the azimuthal angle dependence of the hyperon polarization strongly depends on the choice of the so-called"spin chemical potential"$\Omega_{\mu\nu}$. With $\Omega_{\mu\nu}$ chosen to be proportional to the temperature vorticity, our simulation shows qualitatively coincidental results with the recent measurements at RHIC for both the longitudinal and transverse polarization. We argue that such a coincidence may be related to the fact that the temperature vorticity is approximately conserved in the hot quark-gluon matter.

We revisit the azimuthal angle dependence of the local spin polarization of hyperons in heavy-ion collisions at 200 GeV in the framework of the (3+1)D viscous hydrodynamic model CLVisc. Two different initial conditions are considered in our simulation: the optical Glauber initial condition without initial orbital angular momentum and the AMPT initial condition with an initial orbital angular momentum. We find that the azimuthal angle dependence of the hyperon polarization strongly depends on the choice of the so-called spin chemical potential Ωµν . With Ωµν chosen to be proportional to the temperature vorticity, our simulation shows qualitatively coincidental results with the recent measurements at RHIC for both the longitudinal and transverse polarization. We argue that such a coincidence may be related to the fact that the temperature vorticity is approximately conserved in the hot quark-gluon matter.

I. INTRODUCTION
It is well known that the rotation and spin polarization are correlated and can be converted to each other in materials [1,2]. Recently, the polarization of electrons in a vortical fluid has been observed [3]. Similar phenomena also exist in high-energy heavy-ion collisions in which huge orbital angular momenta (OAM) are produced in peripheral collisions [4][5][6][7][8] (for a recent review, see, e.g. [8]). The huge OAM are distributed into the quark gluon plasma created in heavy-ion collisions in the form of local vorticity [9][10][11][12], which result in the local polarization of hadrons along the vorticity direction [13,14] due to the spin-orbit coupling [4,6]. The net effect of the local polarization at all space-time points on the freeze-out hyper-surface gives the global polarization in the direction of the reaction plane or the OAM of two colliding nuclei [4-6, 8, 15-17].
To describe the STAR data on the global polarization of hyperons which is along the direction of the reaction plane, the hydrodynamic and transport models have been used to calculate the vorticity field [11,12,[56][57][58][59][60][61][62]. In the hydrodynamic framework, the velocity and in turn the vorticity fields in the fireball can be obtained naturally. The transport models describe the phase space evolution of a particle system through collisions among particles, so the position and momentum of each particle in the system at any time is given. To obtain the fluid velocity and then the vorticity at one space-time point, the suitable coarse graining procedure has to be used. Once the vorticity field is obtained, the global polarization of hyperons can be calculated from an integral over the freeze-out hyper-surface which agrees well with the data [61][62][63][64][65].
The polarization of hyperons as a function of the azimuthal angle in the transverse plane has been recently measured in the STAR experiment [19,66]. However the data for the polarization along both the longitudinal and the transverse directions cannot be described by the hydrodynamic models (including A Multi-phase Transport (AMPT) model from which the vorticity field is extracted by the coarse graining method) [62,67,68] based on the coupling of the thermal vorticity and the spin at equilibrium. There is a sign difference between the data and these model calculations. Although one model based on the chiral kinetic theory can explain the sign of the data [69], it cannot reproduce the magnitude of the data. Recent studies showed that the feed-down effects cannot resolve the sign difference [70,71].
The disagreement between theories and experiments indicates that the spin degree of freedom may not reach equilibrium in the fireball and thus the spin polarization may not be determined by the thermal vorticity. The relation between spin and thermal vorticity is dictated by the condition of local thermodynamic equilibrium if the spin tensor does not play a physical role [72]. This calls for new approaches, for examples, the spin can be treated as an independent dynamical variable in the spin kinetic theory and spin hydrodynamics, or dissipative terms should be considered which are possibly larger than believed. Recently, the framework of spin hydrodynamics was developed [73][74][75]. The spin evolution based on particle collisions was derived [76]. The purpose of the present paper is not to make a numerical study based on these new approaches, instead, our purpose is not that ambitious: we will explore different choices of the so-called "spin chemical potential" Ω µν and calculate the corresponding local hyperon polarization. The underlying reason is that, beyond global equilibrium, the thermal vorticity is not guaranteed to be the spin chemical potential, and thus the latter becomes a free parameter [72,74,75]. In the (3+1)D hydrodynamic model CLVisc [77,78], we will assume that the spin chemical potential Ω µν is still determined by the fluid velocity and temperature (or equivalently the energy density). This means that Ω µν , being an anti-symmetric tensor, can be regarded as a type of vorticity (with appropriate normalization to make the dimension correct). We will thus explore four different definitions for Ω µν or vorticity and calculate the local hyperon polarization and compare with the data. In our hydrodynamic simulation, we will examine two different initial conditions: the optical Glauber initial condition without initial OAM and the AMPT initial condition with an initial OAM.
The paper is organized as follows. In Section II we give a brief discussion about our motivation. In Section III we introduce our hydrodynamic model which we use for the simulation. We present our numerical results in Section IV. We give some discussions in Section V. Finally, we give a summary of our results in Section VI.

II. SPIN POLARIZATION AND VORTICITY
The thermodynamic equilibrium in quantum field theory can be described by the density operatorρ. Its form at local equilibrium can be obtained by maximizing the entropy S = −Tr(ρ lnρ) with fixed densities of the energymomentum, the angular momentum, and the conserved charge current on a space-like hyper-surface Σ µ = n µ Σ pointing to a time-like direction n µ [79][80][81][82], n µ Tr(ρT µν ) = n µ T µν , n µ Tr(ρĴ µ,αβ ) = n µ J µ,αβ , whereT µν ,Ĵ µ,αβ andN µ are the density operators of the energy-momentum tensor, the angular momentum tensor, and the conserved charge current, respectively. Note thatT µν is not necessarily symmetric. The quantities T µν , J µ,αβ and N µ are their expectation values. For simplicity we will callT µν (T µν ) andĴ µ,αβ (J µ,αβ ) the energy-momentum and angular momentum tensor respectively though they are actually tensor densities. The angular momentum density operator includes the orbital and spin partŝ Thus, the second constraint in Eq. (1) can be equivalently expressed as n µ Tr(ρŜ µ,αβ ) = n µ S µ,αβ .
The form of the density operator under the constraints (1), or with the second constraint in Eq. (1) being replaced by the constraint (3), that maximizes the entropy readŝ where β ν , Ω αβ and ζ are Lagrangian multipliers which have physical meanings: β ν = u ν /T with u ν being the fourvelocity and T being the temperature, ζ = µ/T with µ being the chemical potential, and Ω αβ plays the role of the chemical potential for the angular momentum 1 . In the following, we will simply call Ω αβ the spin chemical potential as it determines the spin polarization at local equilibrium. The density operatorρ LE in (4) defines the local thermal equilibrium and in general depends on the time.
In relativistic hydrodynamics, in order to obtain the spin vector, we need to first obtain T , u µ , and Ω µν by solving the hydrodynamic equations in which the spin degree of freedom (or equivalently Ω µν ) is treated on the same footing as T and u µ . Such a framework is the spin hydrodynamics [73,75]. However, the numerical spin hydrodynamics has not been established yet. Therefore we will adopt an usual (3+1)D hydrodynamic model, CLVisc [77,78], which can give the space-time evolution of T and u µ . Since Ω µν is antisymmetric, we then assume that Ω µν can be constructed from T and u µ as Ω µν = −(1/2)λ(T )[∂ µ (g(T )u ν ) − ∂ ν (g(T )u µ )] ≡ λ(T )ω µν or its projections where λ and g are scalar functions of T and ω µν is the vorticity tensor 2 . In our numerical simulation, four types of vorticity will be considered, namely, the kinematic vorticity, the relativistic extension of the non-relativistic vorticity (NR vorticity), the thermal vorticity, and the temperature vorticity (T-vorticity). All these four types of vorticities have definite physical meaning and have been widely studied for quite some time in literature. For example, in Ref. [9] the properties of the kinematic, thermal and T-vorticity have been studied in heavy ion collisions, in Ref. [83] the NR vorticity is defined in a study of Lambda polarization in heavy ion collisions.
The kinematic vorticity is defined by where u µ = γ(1, v) is the four-velocity and γ is the Lorentz factor. The mass dimension of the kinematic vorticity is 1. The kinematic vorticity tensor can be decomposed into the part parallel and the part orthogonal to the fluid velocity 1 More precisely, it is T Ω αβ that plays the role of a chemical potential for the angular momentum. 2 In principle, it is also allowed to use the Hodge dual of the vorticity tensor to construct Ωµν . However, when the global equilibrium is approached, it is known that Ωµν should approach the thermal vorticity up to a constant (depending on the symmetry properties of T µν ). We therefore do not consider such a possibility here.
where ε µ = −(1/2)u ν ∂ ν u µ and ω µ = (1/2) µνρσ u ν (∂ ρ u σ ). In comparison with the decomposition of the electromagnetic field strength F µν , the vector ε µ is like an 'electric' field while the vorticity vector ω µ is like a 'magnetic' field. It is clear that ω µ is a direct extension of the vorticity defined in non-relativistic hydrodynamics, ω = (1/2)∇ × v. We thus define the last term in Eq. (6) as the NR vorticity tensor [83] ω (NR) Similar to ω (NR) µν , there has been an attempt to use the spatial components of the thermal vorticity as the spin chemical potential to study the longitudinal spin polarization of hyperons [84].
The temperature vorticity or T-vorticity is defined by where the temperature enters the space-time derivative. The mass dimension of the temperature vorticity is 2. We have decomposed ω µν into the part involving the space-time gradient of the temperature and the part without it which is proportional to the kinematic vorticity.
An important property of the T-vorticity is that it obeys a conservation law [9,11,85]. Suppose Ξ is a twodimensional hyper-surface and C is its boundary, thus the flux of the temperature vorticity on Ξ is equal to the corresponding circulation of T u µ along the boundary Ĉ Since the viscosity of the hot matter in the fireball is small, we can approximately apply the Euler equation for an ideal fluid where ε and P are the energy density and pressure respectively, d/dτ = u µ ∂ µ is co-moving time derivative, and ∇ µ = ∂ µ − u µ (d/dτ ). Rewriting the Euler equation in the following form one easily finds This is the relativistic Helmholtz-Kelvin theorem: the flux of the T-vorticity tensor is conserved with the fluid cell along u µ . We will see that this imposes a strong influence on the spin polarization. The thermal vorticity ω (th) µν is defined by where β = 1/T . The thermal vorticity is dimensionless. Similar to ω (T ) µν into the part involving the space-time gradient of the temperature and the part without it which is proportional to ω (K) µν . We see in Eq. (8) and (13) that ω µν provided that the energy-momentum tensorT µν has a non-vanishing anti-symmetric component [72,75,86]. This can be seen from the following procedure (the analysis based on the dissipative spin hydrodynamics or the kinetic theory gives the same conclusion). The global equilibrium is the state that the density operator (4) becomes independent of the choice of the hyper-surface Σ µ , so that where we used ∂ µŜ µ,αβ =T βα −T αβ . The above condition is fulfilled when 3 Note that ifT µν is symmetric, the spin tensor is conserved ∂ µŜ µ,αβ = 0 and the third condition in (15) does not hold which means that Ω µν remains an independent variable even at global equilibrium. Note that the first line of Eq. (15) is called the Killing equation [87] whose solution is µα are constants. For spin-1/2 fermions at local equilibrium, when ω (th) µν is small, the average spin vector (defined as the Pauli-Lubanski vector) over the hyper-surface Σ µ can be expressed as [13,14,86] where is the Fermi-Dirac distribution, normally Σ µ is chosen as the freeze-out hypersurface for hyperon polarization at the freeze-out. In the calculation we will set ζ = 0 as the net baryon density is almost zero in the hot fireball created in heavy ion collisions at high energies. In this paper, we assume that Eq. (16) can be generalized by replacing ω (th) µν with the spin chemical potential Ω µν as where we will consider four types of vorticities as the spin chemical potentials, namely, ρσ . Here we have chosen suitable factors, λ(T ) = 1/T, 1/T 2 , 1, 1/T , respectively, to make the spin chemical potential dimensionless. Note that Eq. (17) is the main assumption of this paper.
In the following, we will use the (3+1)D hydrodynamic model CLVisc to calculate the four types of vorticities and then use Eq. (17) to obtain the spin vector and then the corresponding spin polarization.

III. HYDRODYNAMIC MODEL
The space-time evolution of the hot quark gluon plasma and dense hadronic matter is described by second order relativistic hydrodynamic equations, where T µν = (ε + P )u µ u ν − P g µν + π µν is the energy-momentum tensor which is symmetric, ε is the local energy density in the co-moving frame of the fluid, P = P (ε) is the pressure determined by the QCD equation of state, u µ is the fluid four-velocity obeying u µ u µ = 1, g µν = diag(1, −1, −1, −1/τ 2 ) is the metric tensor, π µν is the shear-stress tensor whose evolution is solved using a separate group of equations, where η υ is the shear viscous coefficient, σ µν ≡ 2∇ <µ u ν> ≡ 2 µναβ ∇ α u β is the symmetric shear tensor, τ π = 5η υ /(T s) is the relaxation time for the shear viscosity, µν = g µν − u µ u ν is the projection operator that makes the resulting contracted vector orthogonal to u µ , µναβ = 1 2 ( µα νβ + µβ να ) − 1 3 µν αβ is the double projection operator that makes the resulting contracted tensor symmetric, traceless and orthogonal to u µ , θ ≡ ∇ µ u µ is the expansion rate. The operator ∇ µ is the covariant derivative operator defined as for vectors λ µ and tensors λ µν . The Γ's are Christoffel symbols solved as a function of g µν , The above (3+1)D viscous hydrodynamic equations are solved numerically using CLVisc [77,78] with s95p-pce lattice QCD equation-of-state [88], and two different initial conditions: optical Glauber initial condition without initial OAM and AMPT initial condition with initial OAM are applied to check the dependences of the results on initial conditions.

IV. NUMERICAL RESULTS FOR HYPERON POLARIZATION
In this section we will present our numerical results for the polarization of Λ hyperons through vorticity fields by Eq. (17). We choose the coordinate system for collisions of two gold nuclei at 200 GeV in 20-50% centrality, see Fig.  1. The spatial indices µ = 1, 2, 3 in S µ (p) in (17) correspond to the x, y and z direction respectively, so sometimes we write µ = 1, 2, 3 as µ = x, y, z. To test effect of different choices for the spin chemical potential coupled to the spin tensor, we choose four types of vorticities: the kinematic vorticity, the T-vorticity, the thermal vorticity and the NR vorticity. We use the hydrodynamic model CLVisc to compute the vorticity field on the freeze-out hyper-surface.
We use two types of the initial condition: the optical Glauber initial condition without initial OAM and AMPT initial condition with an initial OAM. In the optical Glauber initial condition, the initial energy density distribution is boost invariant at mid-rapidity and is symmetric about the y-axis. As a result, it does not provide any initial OAM. On the other hand, the AMPT initial condition uses HIJING strings. These strings are attached to forward and backward going participants whose distributions are not symmetric about the y-axis in non-central collisions. On one side of the y-z plane, the midpoints of those strings are shifted to forward rapidity in the projectile-going direction. On the other side of the y-z plane, the midpoints of those strings are shifted to backward rapidity in the target-going direction. This forward-backward asymmetry in the AMPT model introduces non-zero initial OAM along the negative y-axis. One should keep in mind that there is no forward-backward asymmetry if the length of strings is infinity (before string fragmentation). In that case, all strings cover mid-rapidity. The system would be perfectly boost invariant along the space-time rapidity and symmetric about the y-axis. This corresponds to extreme high energy collisions where the initial OAM disappears at mid-rapidity. It is consistent with experimental observation that the global polarization is stronger in lower energy collisions. For collisions at low beam energies, the strings from the AMPT model have finite fluctuating lengths [89]. The lengths of strings are determined by the longitudinal-momentum differences between their two end points which are quarks and diquarks from the projectile and the target, whose longitudinal momenta are sampled from parton distribution functions. In this way, the lengths of strings are finite and fluctuating. This helps to propagate the forward-backward asymmetry to left-right asymmetry at mid-rapidity, which is responsible for the initial OAM.
We calculate the transverse and longitudinal polarization of Λ hyperons in the rapidity range Y ∈ [−∆Y /2, ∆Y /2] where S µ (p) is given by Eq. (17) and p denotes the four-momentum of the Λ hyperon, with m T = m 2 Λ + p 2 x + p 2 y . We choose the rapidity range Y ∈ [−1, 1] or ∆Y = 2 in the calculation.

A. Results with optical Glauber initial condition
With the optical Glauber initial condition, we present the results for the longitudinal polarization for Λ. Figure  2 shows P (i) z (p) for four types of spin chemical potentials or vorticities i = K, T, th, NR. The transverse momentum p x and p y all range from −3 to 3 GeV. We see that the T-vorticity has the sign (+, −, +, −) from the first to fourth quadrants consistent with the data. The kinematic, thermal and NR vorticity have the sign (−, +, −, +) opposite in comparison with the data. Note that there is no contribution from the space-time gradient of the temperature in the kinematic vorticity, the temperature gradient part in the T-vorticity has the opposite sign to that in the thermal vorticity. In Fig. 2 we see that the magnitude of P   The longitudinal polarizations from four types of vorticities as functions of azimuthal angles in transverse momenta are shown in Fig. 3. The azimuthal angle relative to the reaction plane (zx-plane) is defined as tan φ p = p y /p x . The azimuthal angle distribution of the polarization − → P (φ p ) is obtained by taking an average over p T = p 2 where ∆p T = p max T − p min T denotes the range of the transverse momentum. In Fig. 3 we see that P (T ) z (φ p ) ∼ sin(2φ p ) which is consistent with the data, while all P µν (T ) have the opposite sign as shown in Eqs. (8,13).  Since there is no initial OAM in the optical Glauber initial condition, the polarizations in the y direction P y (p) are vanishing for all four types of vorticities.

B. Results with AMPT initial condition
In this subsection, we present the results for the AMPT initial condition which encodes the initial OAM of two nuclei in −y direction.
The results of P z (p) and P y (p) for four types of vorticities are shown in Fig. 4 and Fig. 5 respectively. We see in Fig. 4 that the signs of P z (p) with the AMPT initial condition are the same as those with the Glauber initial condition but the magnitudes of P z (p) with the AMPT initial condition are smaller than those with the Glauber initial condition except P z (p) for the NR vorticity that has almost the same magnitude with both initial conditions.
We can take an average over p T for P z (p) in a transverse momentum range to obtain P z (φ p ). The results for P z (φ p ) are shown in Fig. 6 for two transverse momentum ranges. We see that the magnitudes of P z (φ p ) become smaller in the p T range with smaller transverse momenta. For the range p T ∈ [0, 1.2] GeV, the magnitude of P z (φ p ) matches the data. But if we choose p T ∈ [0, 3] GeV, the magnitude of P z (φ p ) is one order of magnitude larger than the data.
In contrast to the vanishing P y (p) with the Glauber initial condition, we obtain finite values of P y (p) in the AMPT initial condition as shown in Fig. 5. The results for P y (φ p ) are displayed in Fig. 7. All four types of vorticities give the correct sign of the initial OAM in −y direction. Note that only P (T ) y (φ p ) for the T-vorticity gives the falling trend in φ p consistent with the data. Although we have the correct trend in φ p in P (T ) y (φ p ), our results fall slower than the data as φ p increases. Our results for P (T ) y (φ p ) match the data at φ p = 0, but at φ p = π/2 our results are P (T ) y (φ p ) ≈ 0.25 while the data approach zero.
We also calculated P y (Y ) in −y direction as functions of the rapidity Y , the results are shown in Fig. 8. We see that P y (Y ) is an even function of Y and increases slowly with |Y |. The values of P y (Y ) are very close for the kinematic and NR-vorticity, and P y (Y ) for the T-vorticity is the largest while that for thermal vorticity is the smallest.     Fig. 4.

C. With different average method in momentum
We can also choose a different method of the average over transverse momenta and rapidity to replace Eqs. (23,25). From Eq. (17) we can take an average of the denominator and numerator separately to obtain the i-th component of the polarization vector, Note that we have introduced an additional p T factor into the p T integrals in both the denominator and numerator since it corresponds to the Lorentz invariant integral d 3 p/E p . The numerical results for P z (φ p ) are presented in Fig.  9. We see that with the same cutoffs for p T , the results for P z (φ p ) from Eq. (26) are a little larger than from Eqs. (23,25). The same behavior also occurs in the results for P y (φ p ) with two different average methods.

V. DISCUSSIONS
We make some remarks about the results. We have checked the first two conditions in (15) and found that they are not fulfilled, so the use of the thermal vorticity as the spin chemical potential is not justified in the hydro-simulation.
For P y (φ p ), we find that only the T-vorticity gives the right trend in φ p comparing to the data, although it decreases slower than the data. Except the trend in φ p , all vorticities can give the global polarization consistent with the data. The reason why the T-vorticity can give the right trend in φ p may be understood as follows. The T-vorticity is conserved so that the T-vorticity flux is frozen in the fluid and move with the fluid cell. In this sense, we can regard the T-vorticity flux as a kind of conserved charge. At the early stage of a non-central collision, the T-vorticity in the out-of-plane direction may be induced by the global OAM, then as the pressure gradient is stronger in the in-plane direction than the out-of-plane direction, the T-vorticity will have a positive elliptic flow which results in the unique φ p dependence as shown in Fig. 7. This suggests that if the spin is (quasi-)conserved, after polarized in the early stage by the OAM, the pressure gradient would lead to a similar φ p dependence as that for the T-vorticity. This may be verified by the simulation using spin hydrodynamics [73,75].
We see very different and even opposite behaviors of P z (φ p ) from different vorticities. This might be related to the fact that P z (φ p ) is one order magnitude smaller than P y (φ p ) since there is no initial OAM in the z direction. Also P z (φ p ) is almost independent of P y (φ p ). This can be seen from the observation that the results of P z (φ p ) from all types of vorticities in the Glauber initial condition (without initial OAM) have the same behaviors as in the AMPT initial condition (with initial OAM). In the optical Glauber initial condition, we found that P y (φ p ) from all types of vorticities are vanishing since there is no orbital angular momentum encoded in the initial state.
Only the T-vorticity in our simulation can describe the data of P z (φ p ) which is the main finding of the paper. The temperature part ω  µν . The implication of the T-vorticity by the data may possibly indicate: (1) The time behavior of the temperature at the freeze-out is essential for the T-vorticity to reproduce the correct sign of P z (φ p ). (2) The T-vorticity might be coupled with the spin in a similar way that a magnetic moment is coupled to a magnetic field. Considering an ideal fluid without a conserved charge density (such as the baryon number density) which is the case in the current hydrosimulation for high energy heavy ion collisions, T u µ can be regarded as a vector potential and the T-vorticity tensor is then the corresponding field strength tensor, so the conservation of T-vorticity flux is similar to the conservation of the magnetic flux in an ideally conducting fluid, see Eqs. (11,12). However, such a picture is not yet rigorously established and it is also unclear how the roles of T-vorticity and thermal vorticity change when the system approaches global equilibrium. Nevertheless, for collisions at lower energies in which the baryon number density is finite, the conservation of the T-vorticity flux does not hold anymore [85]. Thus, the behavior of P z (φ p ) in low energy collisions might provide a test of this point of view. (3) The assumption that the spin chemical potential can be constructed using T and u µ might not be correct, so the fact that the T-vorticity can qualitatively reproduce the experimental data for P z (φ p ) and P y (φ p ) is just accidental. This may be tested by using the spin hydrodynamics which is, however, beyond the scope of this work and we leave it for future. (4) It is also possible that it is a coincidence from the main assumption that the spin vector is given by the T-vorticity in the same way as the thermal vorticity. The true relationship between the spin vector on the freeze-out hyper-surface and all these vorticities is unclear and has to be figured out.
All our results depend on a set of parameters and assumptions. For P z (φ p ) and P y (φ p ), one of the most sensitive parameter is the cutoffs in p T in Eq. (25). For example, as shown in Fig. 6, if we choose the range p T ∈ [0, 1.2] GeV, the theoretical results match the data of P z (φ p ). But if we choose a larger range p T ∈ [0, 3] GeV, our theoretical results are much larger the data of P z (φ p ). The aim of this paper is a qualitative study instead of a quantitative one. We will carry out a detailed and quantitative study of the effects in the future.

VI. SUMMARY
There is a disagreement between theoretical model calculations and recent experimental data about the azimuthal angle dependence of both the longitudinal and transverse polarization of hyperons. These theoretical models are mainly based on the hydrodynamic or kinetic descriptions of the fluid vorticity and express the spin polarization in terms of the thermal vorticity. However, away from global equilibrium, the linear relationship between the spin polarization and thermal vorticity may not be valid (higher order contribution might be relevant). Instead, the spin polarization (or equivalently the spin chemical potential) itself should be regarded as a dynamical variable. Recently there have been attempts in formulating the theory of relativistic hydrodynamics with the spin chemical potential as a (quasi-)hydrodynamic variable, but so far there has been no reliable numerical implementation of the spin hydrodynamics in the market yet.
In this paper, we assume that the spin vector is determined from the spin chemical potential Ω µν in the same way as from the thermal vorticity when the thermal vorticity is small, see Eq. (16) and (17). We also assume that the spin chemical potential Ω µν is still determined by the fluid velocity and temperature, which means that Ω µν can be regarded as being proportional to a type of vorticity. In relativistic hydrodynamics there are various types of vorticities such as the kinematic, temperature and thermal vorticity. There is also a relativistic extension of the non-relativistic vorticity. We thus explore the possibility that the spin chemical potential is proportional to these four vorticities and the spin vector is given by Eq. (17).
We use CLVisc, a (3+1)D viscous hydrodynamic model, to compute the vorticity field. We choose two different initial conditions for the hydro-simulation: the optical Glauber one without initial orbital angular momentum and the AMPT one with an initial orbital angular momentum. We calculated P z (φ p ) and P y (φ p ) as functions of φ p , the azimuthal angle in transverse momentum, for four types of vorticities: the kinematic, temperature, thermal and relativistic extension of the non-relativistic vorticity. Our results show: (1) All types of vorticities have the correct sign of P y for the AMPT initial condition. With the optical Glauber initial condition, they all give vanishing results for P y (φ p ) since there is no orbital angular momentum encoded in the initial state. For P y (φ p ) with the AMPT initial condition, only the temperature vorticity has the same trend as the data, although its magnitude does not agree with the data. (2) For the azimuthal angle distribution in the longitudinal polarization, P z (φ p ), only the temperature vorticity reproduces the sign of the oscillation in the azimuthal angle in data. Other three types of vorticities have a sign difference from the data. (3) The oscillation behavior of P z (φ p ) (not the magnitude) is insensitive to the initial conditions with or without the orbital angular momentum encoded.