Local and global Λ polarization in a vortical fluid

We compute the fermion spin distribution in the vortical ﬂuid created in o ﬀ -central high energy heavy-ion collisions. We employ the event-by-event (3 + 1)D viscous hydrodynamic model. The spin polarization density is proportional to the local ﬂuid vorticity in quantum kinetic theory. As a result of strong collectivity, the spatial distribution of the local vorticity on the freeze-out hyper-surface strongly correlates to the rapidity and azimuthal angle distribution of fermion spins. We investigate the sensitivity of the local polarization to the initial ﬂuid velocity in the hydrodynamic model and compute the global polarization of Λ hyperons by the AMPT model. The energy dependence of the global polarization agrees with the STAR data.


Introduction
Recently the STAR collaboration has measured the polarization of Λ and Λ hyperons [1].They observed that (a) the global polarization decreases with collision energies; (b) the polarization for Λ is always bigger than Λ at the same beam energy.These two features are very important to quantitatively constrain the fluid vorticity of the strongly coupled quark-gluon plasma (sQGP) and the magnitude of the magnetic field through the spin-vorticity and spin-magnetic coupling [2][3][4].The measured beam energy dependence is consistent with the predictions of the hydrodynamic or transport model [2,[5][6][7][8][9][10].The difference between Λ and Λ is caused by (1) pauli-blocking -it is more difficult to polarize Λs than Λs when there are more fermions than anti-fermions [11]; (2) the spin-magnetic coupling -it generates opposite contributions to Λ Email address: lgpang@lbl.gov(Long-Gang Pang* ) arXiv:1704.03569v1 [nucl-th] 11 Apr 2017 and Λ [2][3][4].This difference provides a unique opportunity to determine the magnetic field created in heavyion collisions.In this note, we investigate the effect of the initial flow on the local Λ polarization within a (3+1)D viscous hydrodynamic model and the beam energy dependence of the global Λ polarization within A Multi-Phase Transport (AMPT) [12] model.This study deepens our understanding of the most vortical fluid ever produced [1].

Λ polarization from (3+1)D viscous fluid dynamics
We employ the CLVisc model [10,13] to solve the energy-momentum conservation equation together with extended second-order Israel-Stewart equations for the shear viscosity, with T µν ≡ εu µ u ν − P∆ µν + π µν being the energy-momentum tensor, ε the energy density, u µ the fluid four-velocity normalized as u µ u µ = 1, P the pressure, ∆ µν ≡ g µν − u µ u ν the transverse projector obeying u µ ∆ µν = 0.The η is the shear viscosity, τ π is the relaxation time, is the anti-symmetric vorticity tensor and θ ≡ ∇ µ u µ is the expansion rate.
is the double projector that renders the contracted tensor traceless and transverse to the fluid four-velocity.We observed that the coupling between π µν and ω µν only introduces a small difference.The self-coupling term of ω µν has an undetermined coefficient λ 3 , which leads to a breakdown of the program for fluctuating initial conditions with λ 3 = 1.So these two terms are not included in the present work.The local polarization density on the freeze-out hyper-surface reads [11,14], where m is the fermion mass, dΣ α is the hyper-surface determined by the freeze-out temperature T f = 0.137 GeV, and Ω µν = 1 2 µνρσ ∇ ρ β σ is the thermal vorticity with β σ = u σ /T .Because of strong collectivity, the momentum distribution of the local polarization density is directly related to the local vorticity distribution in space-time [4].The AMPT model is employed to generate the initial energy density with the global angular momentum given by the asymmetry between forward and backward going participants separated by the impact parameter in the transverse plane.We start with the first assumption that at high beam energies there is only Bjorken flow at initial stage which leads to v x = v y = v η = 0.In this case the deposited angular momenta are caused by the asymmetric distribution of the matter for x > 0 (projectile side) and x < 0 (target side), for one specific space-time rapidity.This asymmetry is shown in Fig. 1.The local vorticity is 0 in the initial stage and generated during the fluid evolution.As shown in Fig. 2, the polarization are shifted to the −y direction which is the direction of the global angular momentum.Locally, the azimuthal angle distribution for P y has a cosine structure which indicates a vortex ring [10].The helicity of the vortex ring is opposite for forward and backward rapidity.The polarization at mid-rapidity shows a weaker azimuthal angle dependence.The maximal magnitude of the polarization at 62.4 GeV is about 5 times that at 2.76 TeV.
With initial flow For comparison, we compute the local polarization of Λ by the CLVisc model with the initial transverse flow v x and v y given by the energy-momentum tensor T µν of initial partons by the AMPT.With this configuration, the deposited initial angular momenta are not only given by the matter asymmetry (between x > 0 and x < 0) but also by the v x gradients along η.As a result, the local polarization at Φ = 0 now has the similar magnitude to that at Φ = π as shown in Fig. 3.This dramatic change indicates that the local polarization of Λ hyperons around mid-rapidity may provide rigorous constraints on the initial transverse flow.

Global Λ polarization from AMPT
We also compute the energy dependence of the global Λ polarization within the AMPT model.The fluid velocity and vorticity field are computed from the event average of the four-momentum of all particles in each space-time cell.We consider two impact parameters b = 7 fm and 9 fm at each specific collision energy.The product of the vorticity field Ω zx and the Λ distribution f Λ is integrated over on the hyper-surface to obtain the global Λ polarization.The numerical result with both primary Λ and feed-down contributions agrees with experimental data semi-quantitatively.The energy dependence is investigated in details using the distribution of f Λ and Ω zx .A visualization of these two quantities in the reaction plane indicates that the angular momentum deposition at mid-rapidity is quite small for high beam energies, which is consistent with the Bjorken scaling scenario.This scaling is broken for low beam energies with large asymmetry between the forward and backward going participants, which gives rise to a tilted shape at the mid-rapidity, see Ref. [15] for details.

Fig. 1 .
Fig. 1.The in-plane v x and v η distribution at τ = 0.7 fm after fluid dynamic evolution from fluctuating initial energy density distributions at τ 0 = 0.4 fm from AMPT model.

Fig. 3 .
Fig. 3.The rapidity and azimuthal angle distribution of Λ spin projected to y direction P y for (1) with initial transverse flow, (2) without initial transverse flow.

Fig. 4 .
Fig. 4. (Left) The global polarization of Λ from AMPT model for 2 different impact parameters b = 7 fm and 9 fm (approximating 20 − 50% collisions).(Right) The model data comparison for with and without feed-down contribution from Σ decay.