A study of Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document} and Λ¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\varLambda }$$\end{document} polarization splitting by meson field in PICR hydrodynamic model

With the PICR hydrodynamic model, we study the polarization splitting between Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document} and Λ¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\varLambda }$$\end{document} at RHIC BES energy range, based on the meson field mechanism. Our results fit to the experimental data fairly well. Besides, two unexpected effect emerges: (1) the baryon density gradient has non-trivial and negative contribution to the polarization splitting; (2) for 7.7 GeV Au+Au collisions within the centrality range of 20–50%, the polarization splitting surprisingly increases with the centrality decreases. The second effect might help to explain the significant signal of polarization splitting measured in STAR’s Au+Au 7.7 Gev collisions.


Introduction
Non-central heavy ion collisions create a participant system of extremely hot and dense matter, carrying substantial angular momentum that is perpendicular to the reaction plane [1][2][3]. Through the spin-orbital coupling, just as the Einstein-de-Hass effect [4] and Barnet effect [5] had revealed, the initial fireball angular momentum will eventually give rise to the spin alignment of final particles, such as Λ hyperons [6,7]. The Λ hyperon reveals its polarization by emitting preferentially the weak decay products along its spin direction, and thus is a fairly good choice of polarization measurement in experiments [8][9][10]. Many theories and simulations were also addressing this topic [6,[11][12][13].
Recently, the STAR collaboration measured the nonvanishing Λ polarization for Au-Au collisions at different energies √ s NN = 7.7-200 GeV [14][15][16], and as far as we know, the results conform with the theoretical predictions and simulations in two significant aspects: the global polarization of both Λ andΛ aligns with the initial angular momentum, and decreases with the energy; the local polarization along a e-mail: xieyl@cug.edu.cn (corresponding author) the beam direction shows quadrupolar structure on transverse momentum plane.
However, there still exist some puzzles in this field [17]. Locally, the longitudinal polarization on transverse momentum plane, from model simulations of both a multiple phase transport (AMPT) model [18] and the hydrodynamic model [19,20], exhibits opposite signature to the experimentally observed quadrupolar structure [21]. Many recent works have been devoted to this problem [22][23][24], and the feed-down effect from hyperon decay was proved to be too trivial to explain [25,26].Our recent work [27] using the high resolution (3+1)D Particle-In-Cell Relativistic (PICR) hydrodynamic model to calculate the polarization at 200 GeV Au-Au, shows a fairly good agreement to the experimentally observed longitudinal polarization.
Globally, the magnitude ofΛ polarization is larger than that of Λ polarization. Some might argue that due to the large errors in measurements, it is not sure that whether this polarization splitting really exists, but at least for collision energy of √ s NN = 7.7 GeV, the splitting effect can be identified with high confidence level (see Fig. 4). This splitting effect has raised great interests. It was proposed that the magnetic field induced by the charged spectators can give rise to the polarization splitting between Λ andΛ, but this will require a magnetic field that is long lasting and has a large magnitude, which are not very realistic. Besides, as indicated by Relativistic Magneto-hydrodynamics, the magnetic field can also be induced by charged particles in vortical Quark-Gluon-Plasma (QGP), and in this scenario the magnetic field could last long enough until freeze-out, but problem still exists: the charge density might not be large enough to produce a magnetic field that is strong enough. e.g., the upper limits of the estimated polarization difference at 7.7 GeV is below 1%, which is far away from the lower boundary of experimentally observed 3% difference [28].
Another novel mechanism was proposed by Ref. [29], that the vector meson's "magnetic" field, induced by the baryon vorticity at freeze-out, can split the polarization. However, the polarization splitting formula therein is driven mainly by the directed flow coefficient (c 1 ) and the shear flow coefficient (c 3 ) [29]. The coefficient C, which is proportional to Δc = c 1 − c 3 , is actually a free parameter. Therefore, in this paper, we are going to revisit the theory in Ref. [29] and modify the polarization splitting formula therein, by removing the free parameter C and explicitly bringing out the vorticity, which is essential in Λ polarization study. Then based on this meson field mechanism, we apply the PICR hydrodynamic model that has been previously used in polarization studies, to simulate and calculate the polarization splitting effect.
The paper is organized as follows. In Sect. 2, the theory in Ref. [29] is revisited and modified to explicitly include the vorticity and baryon density gradient. In Sect. 3, we use the PICR hydrodynamic model to simulate and calculate the polarization splitting for Au+Au collisions at RHIC BES energies √ s N N = 7.7-200 GeV. The global polarization and average freeze-out vorticity are also shown, and a special discussion is devoted to the polarization splitting in Au+Au 7.7 GeV collisions. Finally, a summary is drawn. Throughout this paper, we use the natural units:h = c = k B = 1.

Meson field in rotating system
Considering the strong interaction of any fermions mediated by any bosonic fields, one could always write down a general equation of Lagrangian density where L f denotes the Lagrangian density for the fermions, L b represents the Lagrangian density for the bosons, and L int is the interaction Lagrangian density between them. In a simplest case, this equation can be written as: where the first line corresponds to (L f + L int ), denoting the Lagrangian density of Dirac field for fermions with a Yuwaka interaction coupling. The second line corresponds to L b , being the Lagrangian density for the scalar boson σ and vector boson V μ . Here, g σ is the coupling constant between fermion ψ i (of species i) and the scalar boson σ , and g V is the coupling constant between the fermion ψ i and the vector boson V μ . m i , m σ and m V are respectively the mass of baryon, scalar meson and vector meson. The vector meson tensor is: The two constants, f σ and f V in the Yuwaka interaction term are parameters that should be determined case by case.
In relativistic heavy ion collisions, the hyperons are created at the chemical freeze out and then interact with other particles during the hadronic scattering phase. Given that the strong interaction of baryons (including hyperons) with other particles is mediated by a scalar meson field σ and a vector meson field V μ , then with the constants f σ = f V = 1, and following from the Euler-Lagrange equations, one finds the equations of motion for these fields: where n si = ψ ψ is the scalar density of species i, and J μ i = ψ γ μ ψ is the baryon current of species i. These equations are actually the Dirac field equations with scalar and vector field coupling, the Klein-Gordon equation and the Proca equations. The detailed treatments of the above three equations has been demonstrated in Ref. [29].
For the Proca equation (5), analogous to Maxwell equations of massless photon field, it could be decomposed into Maxwell-Proca equations for vector mesons whereḡ V is the mean coupling constant of vector meson, the baryon density is ρ B = i ψ + i ψ i and the baryon (three-)current is J B . These are components of the baryon (four-) Here the E V & B V are the 'electric' and 'magnetic' components of the vector meson field, defined as: where i, j, k = 1, 2, 3. Let us take the curl of Maxwell-Proca equations (7), and we obtain A simple solution was obtained: by neglecting the derivatives in Eqs. (10,11) due to large meson mass, m ω = 783 MeV and m σ = 550 MeV. Assuming global equilibrium of the system, so that ∇ρ = 0, then for the current where ω = ∇ × v is the vorticity of baryon current. Therefore, we could see that the vortical baryon current will induce a vector meson's 'magnetic' field, which, together with its 'electric' component, follow from the Maxwell-Proca equations (6), (7) and definition equations (8), (9). Then the non-relativistic Zeeman energy term in the Foldy-Wouthuysen (FW) Hamiltonian for the hyperon particle's spin (with effective mass M H ) and the vector meson's magnetic fields was written as [29]: where it was argued that the constant matrix β acting on the spin vector S, will result into the opposite signs for Λ andΛ, thus it might be the source of the polarization splitting. Supposing that the spin-1/2 hyperons are in a globally equilibrated system, one could add into the density matrix of the system, ρ, an extra term The ensemble average of the spin vector of spin-1/2 particles are given as S = tr(ρŜ) whereŜ is the spin operator. Then the ensemble averaged polarization vector in Boltzmann statistic limit can be obtained as [30] whereΩ is the unit vector along Ω direction. Taking Eqs. (12) and (13) into the above equation, the polarization splitting would be is a coefficient determined by strong coupling constants, hyperon and meson mass. Hence, if the baryons in high energy collisions have collectively vortical flow motion, the meson interaction with baryons can provide a mechanism for hyperon polarization splitting.
However, the equilibrium reached in high energy collision system is always assumed to be not global, but local. Thus the Eq. (13) is actually local, and should be modified: and the average polarization splitting Eq. (16) becomes: where ... denotes the average over the space. Here ΔP J is the average polarization splitting induced by rotating baryon current J B , ΔP ω = C(ρ B ω)/T originates from the vorticity ω only, and ΔP ρ = C(∇ρ B × v)/T results from the baryon density gradient. In this work, the values of coefficients in Eq. (15) are kept the same as in Ref. [29]: M Λ = 1115.6 MeV, M V = 780 MeV,ḡ V = 5, and g V Λ ≈ 0.55g VN ≈ 4.76. Besides, noting that the Foldy-Wouthuysen transformation used to deduce Eq. (14) is non-relativistic and so are the ensuing Eqs. (15), (16), (18), thus we assume that the post-freeze-out system is near to the Boltzmann limit, and the Λ particles are nonrelativistic. Then to compared with the experimental results, ΔP in Eqs. (15), (16), (18) should be Lorentz-boosted from pre-freeze-out center-of-mass frame into the Λ's rest frame, just like the polarization 3-vector Π( p) is Lorentz-boosted into the particle's rest frame via: However, according to our calculations, the polarization in Λ frame is only ∼ 0.3% smaller than that in QGP frame, or corresponding to only 5-10% correction (especially at the low energies 7.7-30 GeV). This is because most of the particles dwell in the low transverse momentum space, and thus the boost effect is also small. Thus we believe that the present calculations based on Eqs. (15) or (18) are satisfactory quantitative estimates.

7.7-200 GeV
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 [31,32]. In this scenario, the initial state naturally generates longitudinal velocity shear flow, which when placed into the subsequent high resolution (3+1)D PICR hydrodynamic model, will develop into substantial vorticity. Since our initial state+hydrodynamic model describes the shear and vorticity in heavy ion collisions fairly well, its simulations to the Λ polarization also achieved success.
Therefore, we use the PICR hydrodynamic model to simulate the Au+Au collisions at RHIC BES energy region √ s NN = 7.7-200 GeV, and calculate the global polarization with approaches developed in [7], as well as its difference between theΛ and Λ based on Eq. (18).
For the purpose of continuity, we do not perform a new simulation, but just use the same data in our previous Rapid Communication [33], which was then the first work to show the energy dependence of global polarization P H , and seemed to exhibit fairly good agreement with the experimental data. In that work, the simulation parameters were set as follows: the impact parameter ratio was: b 0 = b/b max = 0.7, (where b is the impact parameter and b max is the maximum impact parameter); the cell size was 0.343 3 fm 3 , the time increment is 0.0423 fm/c; the freeze-out time was fixed to be 7.24 fm/c = 2.5 + 4.74 fm/c for all collisions energies (2.5 fm/c for the initial state's stopping time and 4.74 fm/c corresponds to the hydro-evolution time). However, a fixed freeze-out time for different energies is actually not very physical, thus in this work, we are going to vary the freeze-out time for varied collision energies. More specifically, the freeze-out time increases from 5.9 to 7.9 fm/c with the collision energy increasing from 7.7 to 200 GeV, so that the average temperature of the system at freeze-out, as shown in Fig. 1, agrees with the theoretical calculations and experimental results [34,35]. The average baryon densities at the chosen freeze-out time are also shown in Fig. 2, whose values are at the same scale of the freeze-out charge densities in AMPT model as shown in Fig. 2 of Ref. [28].
The red squares in Fig. 3 show the global Λ polarization in our model for Au+Au collisions at √ s N N = 7.7 − 200 GeV, with varied freeze-out time t FO = 5.9-7.9 fm/c. As comparison, we also show in Fig. 3 our previous results with fixed t FO = 7.24 fm/c [33] , by pink stars. One can see that the new values of global polarization is larger than old ones, showing The red stars are extracted from our previous work [33], showing the global polarization as a function of collision energy, with the freeze-out time being fixed to t FO = 7.24 fm/c. The red squares correspond to the case of varied freeze-out time, t FO = 5.9-7.9 fm/c, with varied collision energy √ s = 7.7-200 GeV. The experimental data denoted by up or down triangles are extracted from Ref. [17] the sensitivity of global polarization to the freeze-out time, while the the energy dependence behavior is still kept. The large magnitude of new results are reasonable, since the collision is peripheral with centrality of c = b 2 0 = 49%. The STAR experiment results show that the average global polarization linearly dependent on the centrality, and its value at 20-50% centrality bin is about half of that at 50% centrality [17]. Assuming that the linear dependency of global polarization on centrality is similar for different RHIC BES energies, we estimate the global polarization at 20-50% centrality bin by scaling down the global polarization at b 0 = 0.7 with factor of 0.5, and show them with the dashed line in Fig.  3. One could see the estimated values are very close to the experimental results which are denoted by the up or down triangles. Besides, the correction from the feed-down effect The red squares represent the polarization splitting ΔP J in Eq. (15), with varied freeze-out time t FO = 7.24 fm/c. The green circles correspond to polarization splitting P ω which is induced by vorticity only. As comparison, ΔP J with fixed t FO = 7.24 fm/c are shown by purple stars. The experimental data denoted by cross symbols with error bars are extracted from Ref. [17], as an average for centralities c = 20-50%, corresponding to impact parameter ratio b 0 = 0.45-0.7 [27] of resonance decay would turn down the values by about 15-20% [25,26]. The red squares in Fig. 4 show the average polarization splitting ΔP J of Eq. (15) in our PICR hydrodynamic model for different RHIC BES collision energies, with varied freeze-out time. One can see that our calculation results has the same tendency with the STAR data, which are denoted by cross symbols and error bars, and the magnitudes are also in line with each other if we ignore the different centrality settings (b 0 = 0.7 or c ≈ 50% for our calculation and c = 20-50% for STAR). Besides, our results are similar to the Fig.  1(b) in Ref. [29], which was obtained with the free parameter C assumed to be dependent on collision energy. Actually the parameter C in Ref. [29] should be dependent on collision energy, since C is related to the system's vorticity at freezeout, and as shown by Fig. 5, the vorticity decreases with increasing energy. As comparison, we also show, by purple stars in Fig. 4, the ΔP J with fixed freeze-out time, and they are usually a bit smaller than that with varied freeze-out time.
The green circles in Fig. 4 denotes the polarization splitting ΔP ω induced by vorticity only. We can see that ΔP ω > ΔP J , which means the second term in Eq. (15) induced by baryon density gradient, ΔP ρ , is actually negative and nontrivial. For different collision energies, the ΔP ρ would downplay the final splitting effect for about 1/3 ∼ 1/4. Now we can compare the splitting effect denoted by red squares in Fig. 4, to the calculated polarization denoted by red squares in Fig. 3. One can see that the splitting effect induced by the meson field leads to a limited correction to the global polarization. Taking the case of 11.5 GeV for example: the splitting effect ΔP J with varied FZ time is about 1.33%, and The average y-directed vorticity −ω y at freeze out as function of collision energy with impact parameter ratio b 0 = 0.7. The black and red symbols respectively represent the vorticity at fixed freeze-out time t FO = 7.24 fm/c, and at varied freeze-out time t FO = 5.9-7.9 fm/c the corresponding polarization is about 3.65%. Thus the correction for the Λ polarization is 1.33%/2/3.65% ≈ 18%. Typically, the corrections for the global polarization for different collision energies are about 5-20%, except for 7.7 GeV (about 30%). Figure 5 shows the averaged vorticity along the y direction over the overlapping region as a function of collision energy. It is not surprising that the ydirected vorticity decreases with the collision energy, and the magnitude of the vorticity in our model is the same scale of that from the AMPT simulation. We take the case of √ s NN = 200 GeV for example. In our model the vorticity ω y at freeze-out for √ s NN = 200 GeV with b 0 =0.7 (or b=9.6fm), is about 0.25-0. 35. Meanwhile the AMPT model shows, with b = 9 fm, the vorticity value at late time ( Fig. 11 in ref. [36]) is around 0.2-0.3, and the thermal vorticity ( Fig. 1 in Ref. [37]) is about 0.15. 1 Finally we want to have a little discussion on the polarization splitting at 7.7 GeV. For the case of 7.7 GeV, the polarization difference from our model could be as significant as ΔP J ≈ 3.5%, which is already larger than the lower boundary of experimental measurement of 3%. Up to now several mechanisms were proposed, and quantitative calculations were performed to explain the Λ andΛ polarization splitting [28,29,38,39], but none of them can achieve 3% difference at 7.7 GeV. More specifically, our result for 7.7 GeV case is about 3 times larger than the upper boundary estimate in Ref. [28]. As discussed before that the values of density quantity and the vorticity between our model and the 1 To compare vorticity herein with the thermal vorticity = 1 2 ∇ × (γ v/T ) defined in Ref. [37], one could estimate the freeze-out temperature as around 170 MeV at √ s =39-200 GeV [34], and thus the factor 1 2h T ≈ 0.51 fm/c. AMPT model are very close, the reason why we have much larger polarization splitting effect than that in Ref. [28] lies on the coefficient C = 2(g VHḡV )/(M H m 2 V ), which contains strong coupling constants that are much 'stronger' than the weak coupling constants in Ref. [28].
One might argue that in the more central collisions of b 0 < 0.7, the vorticity will decrease and then the overall polarization splitting at 7.7 GeV would be suppressed to lower than 3%. To deal with this issue, we show the Table 1, where we calculate the polarization splitting ΔP J at different centralities, for Au+Au 7.7 GeV collisions. Surprisingly, within the centrality range of 20-50%, the polarization splitting ΔP J is actually larger for more central collisions. Two factors lead to this unexpected effect: (1) Note that the freeze-out condition herein is a constant temperature for the whole centrality range 20-50%, and as one can see in Table 1, this leads to a smaller freeze-out time for more central collisions than that for peripheral collisions. Meanwhile, Fig. 6 shows that the average vorticity at Au+Au 7.7 GeV collisions in our model has a very mild decreasing tendency with the evolution time (similar behavior was seen by UrQMD model at 2 GeV [40]). Therefore, the vorticity at freeze out for different centralities are rather close, which results into a similar value of ΔP ω in Table 1 for different centralities. (2) Then, the larger fluctuation of baryon density in peripheral collisions means a larger |ΔP ρ |, which of course results into a smaller polarization splitting ΔP J in peripheral collisions.
Therefore it indicates that for Au+Au 7.7 GeV collisions at the centrality range of 20-50%, the polarization splitting will be larger than 3.5%. Furthermore, we have checked that at the energy of √ s N N ≥ 11.5 GeV the above effect no longer exists, because the freeze-out time is larger for more central collisions, and then the vorticity as well as the polarization splitting is smaller for more central collisions. Thus for the energy of √ s N N ≥ 11.5 GeV, the polarization splitting at 20-50% centrality bin will be suppressed compared to that at b 0 = 0.7 (or c=49%) shown in Fig. 4. One can check this from Table 2 for the case of 11.5 GeV. The key point here is that the dependency of the freeze-out time on the centrality at 7.7 GeV is opposite to that at 11.5 GeV and beyond, and we are presently not sure the mechanism behind this. It is possible that our initial state + PICR model loses its validity for low energy collisions of 7.7 GeV, or it might imply the fluid dynamics is different in the collision system of 7.7 GeV and 11.5 GeV. This effect might be exactly the reason why the signal of polarization splitting observed at STAR Au+Au 7.7 GeV collisions seems so strong compared to that at other collision energies, but after all it needs more investigations and confirmations from other models.  Fig. 6 The time evolution of average y-directed vorticity −ω y with different impact parameters in our model for Au+Au collisions at 7.7 GeV Table 2 The average freeze-out temperature T FO , freeze-out time t FO , average vorticity −ω y , and ΔP ω , ΔP J defined in Eq. −ω y (fm −1 ) 0.095 0.120 ΔP ω 1.43% 2.23% ΔP J 0.83% 1.33%

Summary and conclusion
With the PICR hydrodynamic model, we study the polarization splitting between Λ andΛ , based on the meson field mechanism. Our results fit to the experimental data fairly well. Two unexpected effect emerges: (1) the baryon density gradient has non-trivial and negative contribution to the polarization splitting; (2) for 7.7 GeV Au+Au collisions within the centrality range of 20-50%, the polarization splitting surprisingly increases with the centrality decreases. The second effect might hint different hydrodynamics for collision system with energy below 7.7 GeV, and help to explain the significant signal of polarization splitting measured in STAR's Au+Au 7.7 Gev collisions. The reason why the Au+Au collisions at 7.7 GeV in our model has a longer evolution time for peripheral collisions, is still unclear and needs more investigations.