Fluid Dynamics Study of the $\Lambda$ Polarization for Au+Au Collisions at $\sqrt{s_{NN}}=200$ GeV

With a Yang-Mills field, stratified shear flow initial state and a high resolution (3+1)D Particle-in-Cell Relativistic (PICR) hydrodynamic model, we calculate the $\Lambda$ polarization for peripheral Au+Au collisions at RHIC energy of $\sqrt{S_{NN}}=200$ GeV. The obtained longitudinal polarization in our model agrees with the experimental signature and the quadrupole structure on transverse momentum plane. It is found that the relativistic correction (2nd term), arising from expansion and from the time component of the thermal vorticity, plays a crucial role in our results. This term is changing the signature and exceeds the first term, arising from the classical vorticity. Finally, the global polarization in our model shows no significant dependence on rapidity, which agrees with the experimental data. It is also found that the second term flattens the sharp peak arising from the classical vorticity (1st term).


Introduction
In non-central relativistic heavy ion collisions, after penetrating each other, the Lorentz contracted nuclei will break down into quarks and gluons, forming the so called quark gluon plasma (QGP), which carries substantial amount of initial orbital angular momentum [1,2]. The initial shear flow in the viscous QGP will lead to the rotation/vorticity and then the induced vorticity, via the spin-orbit interaction, will eventually give rise to the spin alignment of particles. The non-trivial local and global polarization, which is aligned with the initial angular momentum, was observed in many heavy ion collision experiments [3,4,5] and raising great interest among the physicists [6,7,8,9].
The recent measurement of Λ hyperon polarization in heavy ion collisions launched by the STAR collaboration at, RHIC has shown non-vanishing signals, even at 200 GeV [10,11,12], and in low energy range, the signal could be as significant as 8%. Presently, the experimental signals seem conform with the theoretical calculations and predictions, but still many puzzles remain in this field [13].
One of them is the sign problem: the longitudinal polarization on transverse momentum space shows flipped sign distribution with respect to (w.r.t) the sign distribution of vorticity induced by elliptic flow. Another problem is that theΛ polarization is significantly larger than Λ polarization, e.g. by a factor of 4 at √ s N N = 7.7 GeV, which heretofore has not been predicted and thereafter still has not been interpreted satisfactorily, by theory. There are proposals that the magnetic field induced by the spectators, plays a crucial role in splitting theΛ and Λ polarizations, but the magnitude and the duration of the magnetic field emerging in early stage is still unclear [14,15,16]. However, a recent work proposed that the magnetic field induced by the vortical baryonic QGP of participant system could last much longer time until freeze-out, although the magnitude is much smaller than that induced by charged spectators [17]. More precise measurements of Λ at the RHIC 200GeV Au+Au collisions reveal more disagreement between theory and experiment. E.g., the azimuthal distribution of longitudinal polarization shows flipped sign distribution compared to the hydro-simulated, z-directed polarization over transverse momentum space; the y-directed global polarization is larger in in-plane than in out-of-plane, which is opposite to the hydrodynamic simulations [18,19], as well as transport model results [20]; the global polarization shows no significant dependence on pseudo-rapidity and transverse momentum, while the simulations from a multiple phase transport (AMPT) model shows a normal distribution on pseudo-rapidity and decreasing with the transverse momentum.
In another words, the existing calculations and simulations with the approach developed by F. Becattini et al. [21], using the thermal vorticity at freeze-out to predict the final particle polarization, presently can agree with experiments on the energy dependence behavior of the global polarization and the quadrupole structure of lon-arXiv:1907.00773v1 [hep-ph] 1 Jul 2019 gitudinal polarization. When it comes to the differential measurements, there exists many discrepancies between them, as described above. This indicates some underlying misunderstandings among/inside the theory, simulation and experiments. E.g. the ref. [22] proposed that the theoretically deduced polarization vector should be averaged over the momentum direction of emitted Λ hyperons, to sample the experimental data. This results in a global polarization that is only proportional to the (spatial) thermal vorticity, differing from previous calculations where the temporal component of the relativistic vorticity was also included. As we can see later in this paper, this approach will have a global polarization that is larger in in-plane than out-of-plane, agreeing with the experimental data. Likewise, another ref. [23], also attempted to use the projected thermal vorticity as the source of spin polarization, and surprisingly they seem obtain a correct sign distribution of longitudinal polarization compared to experimental data.
However, up to now, there is still not a comprehensive calculation for the Λ polarization at 200 GeV Au-Au collisions, which we could directly compare to experimental data. Therefore, the main task of this paper is to simulate the 200 GeV Au+Au collisions, using a Yang-Mills fluxtube initial state and a high resolution (3+1)D Particle-in-Cell Relativistic (PICR) hydrodynamics model. We calculate comprehensively the local and global Λ polarization (as function of different variables), that could be compared to the experimental data. By doing this, we expect to see directly how many discrepancies and agreements there exist among the theory, simulations and experimental results.

Simulation and Polarization Vector
The nucleus-nucleus impact in our initial state is divided into many slab-slab collisions, and Yang-Mills flux-tubes. These are assumed to form streaks [24,25]. In this scenario, the initial state naturally generates longitudinal velocity shear flow, which when put into the subsequent high resolution (3+1)D Particle-in-Cell Relativistic (PICR) hydrodynamics model, will develop into substantial vorticity. Since our initial state+hydrodynamic model characterize the shear and vorticity in heavy ion collisions fairly well, its simulations to the Λ polarization has achieved many successes.
As the first hydrodynamic model applied to the polarization study [26], this model was then widely used in our other previous works, e.g. refs. [18,27], and exhibited good descriptions to the vorticity and Λ polarization. In the present work, we chose the modeling parameters as follows: the cell size is 0.343 3 fm 3 , the time increment is 0.0423 fm/c, and the freeze-out (FZ) time is 4+4.91 fm/c (4 fm/c for the initial state's stopping time and 4.91 fm/c corresponds to the hydro-evolution time, which is similar to that in ref. [26]). Based on the simulations to the RHIC's Au+Au collisions at the energy √ s N N = 200 GeV, we calculate the Λ polarization with the widely adopted formula from [26]. where is the inverse temperature four-vector field, and n F (x, p) is the Fermi-Jüttner distribution of the Λ, that is 1/(e β(x)·p−ξ(x) + 1), being ξ(x) = µ(x)/T (x) with µ being the Λ's chemical potential and p its four-momentum. dΣ λ is the freeze out hyper-surface element, for t =const. freeze-out, dΣ λ p λ → dV ε, where ε = p 0 being the Λ's energy. The formula indicates that the polarization originates from the relativistic thermal vorticity defined as: and is proportional to the mean spin vector: where ω 0 = ω 0j = 1 2 (∂ t β + ∇β 0 ) is the temporal component of relativistic thermal vorticity, and ω = ω ij = 1 2 (∇ × β) (i, j =x, y, z) is the spatial component. Therefore the polarization vector is just the normalized spin vector S by weighting over the Λ's number density on the freeze-out surface. E.g. the mean spin vector projected on y direction is, Besides, the polarization vector defined in eq. (1) is divided into 2 terms: the first term arises from the spatial component of relativistic thermal vorticity, and the second term is relativistic modification from the temporal component.
The Λ polarization is determined by measuring the angular distribution of the decay protons in the Λ's rest frame. In this frame the Λ polarization is Π 0 (p), which can be obtained by boosting the polarization Π(p) from the participant frame to the Λ's rest frame, [26], Finally, since the experimental results for Λ polarization are averaged polarizations over the Λ momentum, we evaluated the average of the y component of the polarization Π 0y p . We integrated the y component of the obtained polarization, Π 0y , over the momentum space as follows: Π 0y p = dp dx Π 0y (p, x) n F (x, p) dp dx n F (x, p) = dp Π 0y (p) n F (p) dp n F (p) .

Results and Discussion
We calculate the Λ polarization at 3 impact parameter ratios: b 0 = b/b max = 0.5, 0.6, 0.68, (where b is the impact parameter and b max is the maximum impact parameter), corresponding to 3 centrality points: 25%, 36%, 46%. The freeze-out time is 4+4.91 fm/c, except for the 46% case the FZ time is shorter, i.e. 3.5+4.75fm/c (since smaller system is usually assumed to has shorter evolution time). Fig. 1 shows the centrality dependence of the y-directed polarization boosted into Λ's rest frame, Π 0y . As expected, the y-directed polarization, Π 0y , increases with increasing centrality, since it is already known that the polarization arises from the initial angular momentum, which are related to the impact parameter. This nearly linear dependence on the impact parameter/ centrality were already shown in our previous work [27]. Fig. (1) is to show that our results (black symbols) of the global polarization for Au+Au 200 GeV, are within the boundary of the STAR data (red and blue symbols), exhibiting a good agreement between them. Fig. 2 (upper panel) shows the y component of polarization in the Λ's rest frame, for Au-Au 200 GeV collisions with impact parameter ratio b 0 = 0.6 at the rapidity bin |y| < 1. One can see that the y-directed global polarization in Λ's rest frame, increases in magnitude from in-plane to out-of-plane, contradicting experimental data. This increasing tendency is the same as our previous results for FAIR's U+U 8GeV collisions [18], as well as the other model results [19,26]. However, as suggested in refs. [22,23], we also calculate the global polarization by only considering the first term of the polarization vector in eq. (1) and boosting it into Λ's rest frame, as shown in down panel of Fig. 2. One can see that the polarization decreases from in-plane to out-of-plane, agreeing with the experimental data, whose magnitude, however, is about ten times smaller. Actually this is nothing new, but had been shown in our previous work [26]. Fig. 3 shows the transverse momentum distribution of longitudinal polarization, Π 0z , with impact parameter ratio b 0 = 0.68 for the 200GeV Au+Au collisions. It has the correct sign distribution compared to the experimental data, which is (+, -, +, -) counting from the first coordinate quadrant to fourth quadrant. The peak value at p T = 1.4 GeV is about 0.5%, the same as the global polarization at b = 0.5 fm/c. This is in agreement with the experimental data, i.e. the peak value of longitudinal polarization at p T = 1.4 GeV has similar magnitude with the global polarization.
Furthermore, as shown in Fig. 4 we find that the first term has a sign distribution (-, + , -, +), but the second term has the opposite signature and a larger magnitude, resulting in a (+, -, +, -) sign distribution that agrees with experimental data. Comparing to our previous results for FAIR's U+U 8GeV collisions [18], the first term   keeps the same sign distribution, i.e. (-, + , -, +), but with magnitude growing from about 2% to 8% at large transverse momentum. Meanwhile, the second term flips it sign distribution, from (-,+ ,-,+) to (+, -, +, -), and grows faster to a magnitude of 12%, which is larger than the first term. Two points are worthy to be noticed here: (1) The magnitude, of either the first/second term or the total of longitudinal polarization, increases from low energy (8 GeV) to high energy (200 GeV). This seems contradicts with a previous work [28], where the second harmonic coefficient of the longitudinal polarization decreases with energy increasing from 7.7 GeV to 2.76 TeV; (2) The second term, in our model, plays crucial role to obtain the experimentally observed sign structure and magnitude of the longitudinal polarization: it has a sign structure of (+, -, +, -), and a larger magnitude, covering the first term's opposite signature and amending the polarization value into a smaller but correct magnitude. This is similar to ref. [23], where the total longitudinal polarization flips its sign distribution with respect to that of the first term, although the signatures therein are just opposite to our results.
Then we explore also the global polarization as a function of rapidity, as shown in the lower panel of Fig. 5. The red dashed line in the lower panel figure is a rough approximation of the experimental data, which shows no significant dependence on the rapidity and fluctuates around the averaged value 3%. One can see that the global polarization from our model also shows no significant dependence on the rapidity. The global polarization, Π 0y , for b 0 = 0.5, 0.6, 0.68, fluctuates around the average value of 2.8%, 3.8% and 6% respectively, which are magnitudes similar to the global polarization. For more peripheral collisions, the fluctuations are relatively larger, e.g. at the case of b 0 = 0.68, there exists a dip in rapidity bin |y| < 0.4. Beyond the rapidity range |y| > 1 the global polarization goes down rapidly to zero. The first term of the y-directed polarization, as shown in the upper panel of Fig. 5, exhibits a normal distribution with respect to the rapidity, with peak value at center rapidity y = 0, which is similar to the vorticity distribution on pseudo-rapidity from AMPT model [29]. This similarity of structure simply demonstrates the definition of polarization vector's first term, i.e. Π 1y arises purely from the spatial component of relativistic vorticity, ω = 1 2 ∇ × β. For more peripheral collisions with larger impact parameter, the global polarization distribution peaks higher at center rapidity y = 0 and goes down faster to zero with a narrower width. Finally, the two figures together indicate that the second term related to the system expansion, flatten the peak of the first term induced by classical vorticity, resulting in an even distribution of global polarization on the rapidity.

Summary and Conclusions
With a Yang-Mills field stratified shear flow initial state and a high resolution (3+1)D Particle-in-Cell Relativistic (PICR) hydrodynamic model, we calculate the Λ polarization for Au-Au collisions at RHIC energy of √ S N N = 200 GeV. The transverse momentum distribution of global polarization shows a magnitude increasing from in-plane to out-of-plane, contradicting to the experimental data. However, the longitudinal polarization in our model shows correct signature of the experimentally observed quadrupole structure on transverse momentum plane. Besides, the peak value of the longitudinal polarization at p t = 1.4 GeV is similar to the global polarization, which is in good agreement with the experimental data. When delving into the two terms of the polarization vector, it is found that the second term arising from system expansion or the temporal component of the relativistic thermal vorticity, plays a crucial role to obtain our results, by changing the signature and magnitude of the first term which is induced by classic vorticity. Furthermore, we also plot the global polarization as a function of rapidity. Our results show no significant dependence on rapidity, which again conform with the experimental data, and again it is found that it is the second term that flattens the sharp peak of the first term induced by classical vorticity.