Arxiv:1711.06008V2 [Nucl-Th] 16 Mar 2018 H Lblplrzto.Rcnl H Trclaoainh Momentu Collaboration STAR Angular the Orbital Recently the Polarization

Home , Vector meson

arXiv:1711.06008v2 [nucl-th] 16 Mar 2018 h lblplrzto.Rcnl h TRclaoainh momentu collaboration STAR angular the orbital Recently the polarization. of transfe global direction are the momenta the angular along the quarks of of Some review). recent a for h emeeg cnporm[,9.A l nrisblw62. below energies all At 9]. [8, program for scan energy beam the hs uigteerysaeo h vlto 1,2–6.Th 20–26]. So 28]. [10, [27, evolution energies. moments the high of magnetic at stage momentatheir other early angular each the orbital during through phase the pass as nuclei two direction after same the in average oaiainpeoeo ntefia tt fhdos otp Most hadrons. of state final Λ the in phenomenon polarization sWge ucinapoc 2,2,5–5.Tepolarizati The 51–55]. 27, qua [26, field. the approach from magnetic momenta function their Wigner s of coale and functions as this as octet With baryons spin-1/2 and straightforward. the mesons are of [ multiplets baryons baryon models and and recombination mesons or vector on 48] focus [47, qua models polarized coalescence of pola coalescence quark the through compute can decuplet we spin-3/2 which with and quark matrix, a density spin such the in on included not comput was to dependence used was momentum f model But polarization recombination hyperon [2]. quark the same of The estimate an model. us, [42–44] (Fermi-Dirac) of hydrodynamic one thermal by and the Liang obtained are from which arise fields to vorticity assumed are hadrons h otct n rmtemgei edsprtl.W also We appro separately. some field in magnetic discuss the we Sect ρ from VI In and Section vorticity decuplet. In the we spin-3/2 functions. and IV, Wigner Section octet from spin-1/2 In qua the a dependence. in formulate momentum baryons we explicit III, with Section mesons In polarization. particle’s eto egv oebsc bu h pndniymatrix. density spin the about basics some give we section 00 h ytmi o-eta ev-o olsosa ihener high at collisions heavy-ion non-central in system The togmgei ed r lopoue nnncnrlheav non-central in produced also are fields magnetic Strong o hs oaie ursfr n eaet oaie hadr polarized to relate and form quarks polarized these How nti ae ewl omlt o-eaiitcqakcoa quark non-relativistic a formulate will we paper this In h ae sognzda olw.I eto I eintroduc we II, section In follows. as organized is paper The ic emk s ftesi est arxtruhu h pa the throughout matrix density spin the of use make we Since and Λ o etrmsn.I h atscinw ocueb iiga giving by conclude we section last the In mesons. vector for 3 ula cec iiin S7R39 arneBree Na Berkeley Lawrence 70R0319, MS Division, Science Nuclear and 1 Λ ¯ eateto oenPyis nvriyo cec n Tec and Science of University Physics, Modern of Department yeos[93]aebsdo ttsia-yr oes[35 models statistical-hydro on based are [29–34] hyperons Λ ¯ ur olsec oe o oaie etrmsn n ba and mesons vector polarized for model coalescence quark A ihepii oetmdpnec.Terslsaecompare densi are f hadrons. spin results baryons for The the and models statistical mesons With dependence. vector momentum explicit of decuplet. with polarizations spin-3/2 the way and uniform octet spin-1/2 of aebe esrd[8]. measured been have o-eaiitcqakcaecnemdli omltdf formulated is model coalescence quark non-relativistic A 2 e aoaoyo ur n etnPyis(O)adInstitu and (MOE) Physics Lepton and Quark of Laboratory Key aggagYang, Yang-guang eta hn omlUiest,Whn 309 China 430079, Wuhan, University, Normal China Central I AISO PNDNIYMATRIX DENSITY SPIN OF BASICS II. 1 e-ogFang, Ren-hong .INTRODUCTION I. 2 kcaecnemdlfrtesi est arxo vector of matrix density spin the for model coalescence rk smaue h lblplrzto of polarization global the measured as k n niqak nasseai a.Teconventional The way. systematic a in anti-quarks and rks u Wang, Qun ecnemdlwt xlctmmnu eednebased dependence momentum explicit with model lescence cnemdlw a opt h oaiain fvector of polarizations the compute can we model scence ishv ag ria nua oet 16 seRf [7] Ref. (see [1–6] momenta angular orbital large have gies togmgei ed a loplrz ursthrough quarks polarize also can fields magnetic strong e rtasotmdl 3,4,4] na al ae 1 by [1] paper early an In 46]. 45, [32, models transport or o oaie urswsmd naqakrecombination quark a in made was quarks polarized rom rteeetpae hstp fplrzto scalled is polarization of type This plane. event the or m rdt h o n es atradla opolarization to lead and matter dense and hot the to rred 9 0 onticuetesi ere ffedm We freedom. of degrees spin the include not do 50] 49, o eso o ocmueteqakpolarization quark the compute to how show we V ion n nld hs nue ytevriiyadb the by and vorticity the by induced those include ons iain fvco eosadbroso pn12octet spin-1/2 of baryons and mesons vector of rizations o olsos[01] h antcfilsaeo the on are fields magnetic The [10–19]. collisions ion y hyaeepce ohv ffcsol ntepartonic the in only effects have to expected are they iesm etrso h pndniymti element matrix density spin the of features some give e,nnvnsigvle ftegoa polarization global the of values non-vanishing GeV, 4 u uhsrn antcfilsdcyvr quickly very decay fields magnetic strong such But . eiu acltoso h lblplrztoso the of polarizations global the of calculations revious itiuinwt pndgeso reo u othe to due freedom of degrees spin with distribution h est arxwihw iluet opt a compute to use will we which matrix density the e kplrztosotie rmohrapoce such approaches other from obtained polarizations rk i-/ eult h eeaiaint te meson other to generalization The decuplet. pin-3/2 inlLbrtr,Bree,Clfri 42,USA 94720, California Berkeley, Laboratory, tional iain h ur n arnplrztosfrom polarizations hadron and quark the ximations h pndniymti lmnso etrmesons vector of elements matrix density spin the e omlt h aemdlfrteplrztosof polarizations the for model same the formulate nlg fCia ee,Ahi202,China 230026, Anhui Hefei, China, of hnology umr ftemi results. main the of summary eobnto model. recombination n sa motn usint nesadthe understand to question important an is ons rplrzdvco eosadbaryons and mesons vector polarized or e ocmuetepril oaiain nthis in polarization, particle the compute to per 4] nteecluain h oaiain of polarizations the calculations these In –41]. o hs fqak n antiquarks and quarks of those rom ota bandfo iei and kinetic from obtained that to d 1 ymti,oecncmuei a in compute can one matrix, ty n i-inWang Xin-nian and eo atcePhysics, Particle of te ,3 2, Λ ryons and Λ ¯ yeosin hyperons 2

We consider an ensemble of particles with spin quantum number S. The normalized spin states are labeled by ψi from which the spin density operator can be deﬁned | i

ρ = P ψ ψ , (1) i| iih i| Xi where Pi is the probability of the spin state ψi . We list here the properties of the spin density operator [56]: (a) ρ is Hermitian, ρ = ρ†; (b) The trace of ρ is 1, Tr| ρ =i ψ ρ ψ = P =1; (c) ρ is a positive-semidefinite operator, j h j | | j i i i P 2 P i.e., for any state φ , we have φ ρ φ = P φ ψ > 0. | i h | | i i i h | ii P The dimension of the spin space for a spin-S particle is (2S +1), so the spin density operator is a (2S +1) (2S +1) matrix with 2(2S + 1)2 real parameters. The conditions that ρ is Hermitian and has unit trace reduce the number× of real parameters to 4S(S+1). The positive-semidefinite condition can also impose restrictions on these real parameters. If the particle is unpolarized, the spin density matrix ρ has a simple form 1 ρ = diag (1, 1, , 1). (2) 2S +1 ··· Any deviation from (2) indicates some degree of spin alignment or polarization. Let us look at a few examples. The first example is the spin-1/2 particle. The spin density matrix has three real parameters and can be written as 1 ρ = (1 + −→P σ), (3) 2 · where −→P is a unit three-vector serving as the polarization vector and σ = (σ1, σ2, σ3) are Pauli matrices. The polarization can be read out by Tr(ρσ) −→P = . (4) Tr(ρ) The octet hyperon Λ0 is a spin-1/2 particle whose polarization can be measured by its weak decay, Λ0 p+ π− [9, 57]. In Section IV we will use the quark coalescence model to compute the polarization of Λ0 through→ its spin density matrix. The second example is the spin-1 particle whose spin density matrix has eight real parameters. We look at the vector meson K∗0 and φ for illustration with following strong decay modes K∗0 K+ + π−, ( 100%), → ∼ φ K+ + K−, ( 49%), (5) → ∼ where the fractions are the branching ratios. The polarization of K∗0 and φ can be measured through these decays [58, 59]. We note that parity is conserved in strong decays so the final states must be in the p-wave with L =1. We take K∗0 for an example, all discussions apply to φ equally. The decay amplitude can be described by the S matrix K+, π− S K∗0; S = Y (θ, φ) (6) h | | zi 1,Sz where Sz = 1, 0, 1, (θ, φ) denote the polar and azimuthal angle of one decay daughter, and Y1,Sz (θ, φ) denote the − + − spherical harmonic functions of L =1. The angular distribution of K or π for a specific initial spin state 1,Sz of K∗0 is | i dN 2 = K+, π− S K∗0; S = Y (θ, φ) 2 . (7) dΩ h | | zi | 1,Sz |

∗0 We now consider an ensemble of K at rest having the probability Pi in the spin state ψi . Then the spin density operator ρ is given by Eq. (1). The angular distribution for the decay daughter can be written| i as

dN 2 = P K+, π− S ψ dΩ i h | | ii Xi

∗0 ∗0 Inserting a completeness relation K ; Sz K ; Sz =1 into (8), we obtain Eq. (10) of Ref. [60], Sz | ih | P dN = K+, π− S K∗0; S dΩ h | | z1i SzX1,Sz2 K∗0; S ρ K∗0; S K∗0; S S † K+, π− h z1| | z2ih z2| | i ∗ = ρSz1,Sz2 Y1,Sz1 (θ, φ)Y1,Sz2 (θ, φ) SzX1,Sz2 3 = (1 ρ )+(3ρ 1)cos2 θ 8π − 00 00 − 2 2 2Reρ− sin θ cos(2φ) 2Imρ− sin θ sin(2φ) − 1,1 − 1,1 +√2Re(ρ− ρ ) sin(2θ)cos φ 1,0 − 01 +√2Im(ρ−1,0 ρ01) sin(2θ) sin φ , (9) − i where ρ K∗0; S ρ K∗0; S denote the spin density matrix elements for K∗0. Note that dN/dΩ is auto- Sz1,Sz2 ≡ h z1| | z2i matically normalized to 1, dΩ(dN/dΩ) = 1. If the parent particle is unpolarized, ρSz1,Sz2 = δSz1,Sz2 /3, then we have dN/dΩ=1/(4π). ´ In principle, if one measures dN/dΩ in experiments, one can determine ﬁve real parameters out of eight ones in the spin density matrix for vector mesons according to Eq. (9),

ρ , Reρ− , Imρ− , Re(ρ− ρ ), Im(ρ− ρ ). (10) 00 1,1 1,1 1,0 − 01 1,0 − 01 The polarization vector −→P can be determined from the spin density matrix [30]

1 Tr(ρSˆ) −→P = , (11) S Tr(ρ) where S =1 and the spin operators for vector mesons, Sˆ = (Sˆ1, Sˆ2, Sˆ3), are deﬁned by

0 1 0 0 i 0 10 0 1 1 − Sˆ1 =  1 0 1  , Sˆ2 =  i 0 i  , Sˆ3 =  00 0  . (12) √2 0 1 0 √2 0 i −0 0 0 1      − 

Then we obtain the polarization −→P = (P1, P2, P3), 1 P = √2 Re(ρ + ρ ), 1 Tr(ρ) −1,0 01 1 P = √2 Im(ρ + ρ ), 2 Tr(ρ) −1,0 01 1 P = (ρ ρ− − ) . (13) 3 Tr(ρ) 11 − 1, 1

Note that the above do not match the ﬁve quantities in (10) that can be measured from dN/dΩ. Therefore the polarization of vector mesons cannot be measured from dN/dΩ in strong decays in which the parity is conserved. This is very diﬀerent from Λ0 in the weak decay where the broken parity can be used to determine the polarization. When the azimuthal angle φ in Eq. (9) is integrated out, we obtain the polar angle distribution

dN 2π dN 3 = dφ = (1 ρ )+(3ρ 1)cos2 θ . (14) d cos θ ˆ dΩ 4 − 00 00 − 0

So one can measure ρ00 from the polar angle distribution of the daughter particle. Any deviation from 1/3 for ρ00 may indicate some degree of polarization of vector mesons. This type of measurements have been performed in the STAR experiment [58]. The third example is the spin-3/2 decuplet baryon ∆++ which mostly decays into a proton and a π+ by strong interaction. The spin-parity of ∆++ is (3/2)+, while those of the proton and pion are (1/2)+ and 0− respectively. By the parity conservation in the strong decay the angular distribution of the proton must be in the P-wave state with 4

L =1. The spin density matrix of the spin-3/2 particle like ∆++ is a 4 4 complex matrix with ﬁfteen independent parameters. The angular distribution for the decay daughter can be written× as

dN p + S S † p + = pˆ,Sz ; π ρ pˆ,Sz ; π dΩ p h | | i XSz p ∗ p = ρSz1,Sz2 f(Sz ,Sz1)f (Sz ,Sz2), (15) p SzX1,Sz2 XSz

p + where the final state of the proton and pion is denoted as pˆ,Sz ; π with pˆ being the proton’s momentum direction p | i and Sz denoting the proton’s spin state along the z-direction. In deriving the second line of Eq. (15) we have inserted ++ ++ ++ into the first line of Eq. (15) a complete set of spin states of ∆ , Sz ∆ ,Sz ∆ ,Sz = 1, where Sz denotes p | i h | the spin state along the z-direction. In Eq. (15), f(Sz ,Sz) is the decayP amplitude given by Eq. (A3) and ρSz1,Sz2 are ++ the spin density matrix elements for ∆ given by Eq. (A6). The summation over Sz1 and Sz2 in Eq. (15) is made ++ 1 3 for all spin states of ∆ with Sz1,Sz2 = 2 , 2 . The final form of the angular distribution in Eq. (15) is shown in Eq. (A5). We see that the angular distribution± ± involves only five parameters out of fifteen independent ones, the rest ten parameters cannot be determined by the daughter particle’s mometum distribution alone. The polarization vector for ∆++ can be obtained from Eq. (11) and the result is given by Eq. (A8). We see that the parameters appearing in the polarization vector are independent of those in the angular distribution (A5). This is the feature of the parity conserved decay of ∆++, same as in the strong decay of K∗0 in (5).

III. SPIN DENSITY MATRIX FOR VECTOR MESONS

We consider an ensemble of free quarks and antiquarks in a deconfined system. These quarks and antiquarks will combine to form hadrons in the hadronization process. For mesons we define a quark-antiquark state with momenta and spins along a fixed direction (the z-direction)

q , ¯q ; s ,s ; p , p q , ¯q ; s ,s q; p , p , (16) | 1 2 1 2 1 2i ≡ | 1 2 1 2i | 1 2i where q1 = u, d, s and ¯q2 = ¯u, d¯,¯s denote the ﬂavors of the quark and anti-quark respectively, s1,s2 = 1/2 denote the spin of the quark and antiquark in the z-direction respectively, and ’q’ labels the quark-antiquark± momentum state. The momentum states of the quark and antiquark in coordinate representation are plane waves, 1 x , x q; p , p = exp(ip x + ip x ) , (17) h 1 2| 1 2i V 1 · 1 2 · 2 where V is the volume. Now we can write down the density operator for quarks and antiquarks,

d3p d3p ρ = V 2 1 2 w (p )w (p ) ˆ (2π)3 (2π)3 q1,s1 1 ¯q2,s2 2 sX1,s2 qX1,¯q2 q , ¯q ; s ,s ; p , p q , ¯q ; s ,s ; p , p , (18) × | 1 2 1 2 1 2i h 1 2 1 2 1 2| where wq1,s1(p1)/w¯q2,s2(p2) are the weight functions for the quark/antiquark which satisfy the normalization conditions

wq,s(p)=1, w¯q,s(p)=1. (19) Xs Xs The weight functions are related to the polarization of the quark and antiquark by 1 w (p) = [1 P (p)] , q,±1/2 2 ± q 1 w (p) = [1 P (p)] , (20) ¯q,±1/2 2 ± ¯q where Pq(p) and P¯q(p) are the polarization of the quark and antiquark in the z-direction respectively. We will use shorthand notations for the weight functions, wq(¯q),± wq(¯q),±1/2, in the rest of the paper. Now we look at the meson state deﬁned by ≡

M; S,S , p M; S,S M; p , (21) | z i ≡ | zi | i 5 where M; S,S denotes the meson’s spin-ﬂavor wave function and ’M’ labels the meson state. For a vector meson | zi we have S =1 and Sz = 1, 0, 1 for the spin state in the z-direction. In the quark model, the momentum state of the meson in coordinate representation− is given by 1 x , x M; p = exp(ip x)ϕ (y). (22) h 1 2| i V 1/2 · M

Here x = (x1 + x2)/2 is the center position of the quark and antiquark, y = x1 x2 is their distance, and ϕM(y) is the meson wave function satisfying the normailzation condition d3y ϕ (y) 2 =1−. We choose the Gaussian form for | M | ϕM(y) whose Fourier transform is given by [47, 48] ´

ϕ (q) = d3ye−iq·yϕ (y) M ˆ M 2√π 3/2 q2 = exp 2 , (23) aM −2aM where aM is the coefficient characterizing the average width of the momentum in the meson wave function. One can check that ϕ (q) satisfies the normailzation condition d3q ϕ (q) 2 = (2π)3. M | M | We can compute the spin density matrix element ρS´z1Sz2 on the vector meson state with, for example, the quark flavor q¯q, where Sz1,Sz2 = 1, 0, 1. The element ρ00 for Sz1 = Sz2 =0 is of special interest since it can be measured in experiments (as discussed− in Section II), ρS=1(p) = M; 1, 0, p ρ M; 1, 0, p 00 h | | i 1 d3p d3p = V 2 1 2 2 ˆ (2π)3 (2π)3 2 [w (p )w −(p )+ w −(p )w (p )] q; p , p M; p × q,+ 1 ¯q, 2 q, 1 ¯q,+ 2 |h 1 2| i| 1 d3q p p q q = 3 wq,+ + w¯q,− 2 ˆ (2π) h 2 2 − p p 2 +wq,− + q w¯q,+ q ϕM(q) , (24) 2 2 − i | | where we have used the Clebsch-Gordan coefficients in the quark spin wave function of the vector meson state S =1,S =0 and the amplitude of a vector meson state projected onto the quark-antiquark state | z i q; p , p M; p = d3x d3x q; p , p x , x x , x M; p h 1 2| i ˆ 1 2 h 1 2| 1 2i h 1 2| i (2π)3 = δ(3)(p p p )ϕ (q), (25) V 3/2 − 1 − 2 M where we have used Eqs. (17,22,23) with q = (p p )/2. From Eq. (25) we obtain 1 − 2 (2π)3 q; p , p M; p 2 = δ(3)(p p p ) ϕ (q) 2 , (26) |h 1 2| i| V 2 − 1 − 2 | M | where we have used δ(3)(0)= V/(2π)3. Eq. (26) has been used in Eq. (24). Similarly we can obtain other elements of the spin density matrix for vector mesons as well as for scalar (spin-0) mesons. We list all diagonal elements as follows, 1 d3q p p S=0 p q q ρ00 ( ) = 3 wq,+ + w¯q,− 2 ˆ (2π) h 2 2 − p p 2 +wq,− + q w¯q,+ q ϕM(q) , 2 2 − i | | 1 d3q p p S=1 p q q ρ00 ( ) = 3 wq,+ + w¯q,− 2 ˆ (2π) h 2 2 − p p 2 +wq,− + q w¯q,+ q ϕM(q) , 2 2 − i | | d3q p p S=1 p q q q 2 ρ11 ( ) = 3 wq,+ + w¯q,+ ϕM( ) , ˆ (2π) 2 2 − | | d3q p p S=1 p q q q 2 ρ−1,−1( ) = 3 wq,− + w¯q,− ϕM( ) . (27) ˆ (2π) 2 2 − | | 6

S=1 Note that all oﬀ-diagonal elements for the spin triplet are vanishing, i.e. ρSz1Sz2 = 0 for Sz1 = Sz2. Since the trace of ρ runs over all meson spin states, it fulﬁlls 6

S=0 S=1 S=1 S=1 ρ00 (p)+ ρ00 (p)+ ρ11 (p)+ ρ−1,−1(p)=1. (28)

S=1 p If we only consider vector mesons, we can obtain the normalized spin density matrix element ρ¯SzSz ( ) (Sz = 1, 0, 1) for the vector meson which is related to the quark and antiquark polarization functions. The normalized spin− density S=1 matrix element ρ¯00 (p) is given by d3q[1 P (p/2+ q)P (p/2 q)] ϕ (q) 2 S=1 p q ¯q M ρ¯00 ( ) = 3 − − | |2 . (29) ´ d q[3 + Pq(p/2+ q)P¯q(p/2 q)] ϕM(q) ´ − | | P P S=1 We see that if the quark and antiquark are not polarized ( q = ¯q = 0), we have ρ¯00 = 1/3. If the quark polarization is small, P P 1, Eq. (29) can be approximated as q ∼ ¯q ≪ 1 4 d3q p p S=1 p P q P q q 2 ρ¯00 ( ) 3 q + ¯q ϕM( ) . (30) ≈ 3 − 9 ˆ (2π) 2 2 − | | If the quark polarization is independent of momentum, Eq. (30) becomes 1 4 ρ¯S=1 P P . (31) 00 ≈ 3 − 9 q ¯q S=1 P P From Eqs. (30,31) we see that the deviation from non-polarized case ρ¯00 =1/3 is of second order in q or ¯q. For φ(1020) with the quark content s¯s, we obtain its polarization along z-direction from Eq. (13),

ρφ ρφ P 11 −1−1 φ = φ −φ φ , (32) ρ11 + ρ00 + ρ−1−1 where we have from Eq. (27),

1 d3q p p φ φ q 2 P q P q ρ11 ρ−1,−1 = 3 ϕM( ) s + + ¯s , − 2 ˆ (2π) | | h 2 2 − i 3 1 d3q p p φ φ φ q 2P q P q ρ11 + ρ00 + ρ−1,−1 = + 3 ϕM( ) s + ¯s . 4 4 ˆ (2π) | | 2 2 − φ P P P P For the φ(1020) meson, ρ00 is given by Eq. (30) or Eq. (31) with replacement q s and ¯q ¯s. If the quark polarization is independent of momentum, Eq. (32) becomes → →

2(Ps + P¯s) Pφ = , (33) 3+ PsP¯s

From Section II we know that one cannot measure the polarization Pφ in the strong decay (5) except the element ρ00, φ P P see Eq. (14). However, by measuring ρ00 one can determine the values of s and ¯s from which one can indirectly determine Pφ.

IV. SPIN DENSITY MATRIX FOR OCTET AND DECUPLET BARYONS

Similar to the quark-antiquark state defined in Eqs. (16,17) for a meson, we can define a three-quark state for a baryon in the deconfined quark system,

q , q , q ; s ,s ,s ; p , p , p q ,s , p q ,s , p q ,s , p | 1 2 3 1 2 3 1 2 3i ≡ | 1 1 1i | 2 2 2i | 3 3 3i = q , q , q ; s ,s ,s q; p , p , p , | 1 2 3 1 2 3i | 1 2 3i 1 x , x , x q; p , p , p = exp(ip x + ip x + ip x ) , (34) h 1 2 3| 1 2 3i V 3/2 1 · 1 2 · 2 3 · 3 where s1,2,3 = 1/2 and q1,2,3 = u, d, s denote the spin states in the z-direction and flavors of three quarks respectively. The first line implies± that the spin and flavor part of the three-quark state is independent of its momentum part. The 7 second line means that three quarks are assumed to be in plane waves. The spin density operator for three quarks has the form d3p d3p d3p ρ = V 3 1 2 3 w (p )w (p )w (p ) ˆ (2π)3 (2π)3 (2π)3 q1,s1 1 q2,s2 2 q3,s3 3 s1X,s2,s3 q1X,q2,q3 q , q , q ; s ,s ,s ; p , p , p q , q , q ; s ,s ,s ; p , p , p , (35) × | 1 2 3 1 2 3 1 2 3i h 1 2 3 1 2 3 1 2 3| where wq,s are the weight functions for quarks related to the quark polarization Pq by Eq. (20). Now let us look at a baryon state defined by

B; S,S , p = B; S,S B; p , (36) | z i | zi | i where B; S,S is the spin-ﬂavor wave function of the baryon and B; p is the baryon wave function in momentum | zi | i space. For a spin-1/2 baryon, we have S =1/2 and Sz = 1/2. The momentum state of the baryon has the following form in coordinate representation ± 1 x , x , x B; p = exp(ip x)ϕ (y, z), (37) h 1 2 3| i V 1/2 · B where x = (x1 + x2 + x3)/3, y = (x1 + x2)/2 x3, and z = x1 x2 whose conjugate momenta are p = p1 + p2 + p3, r = (p + p 2p )/3, and q = (p p )/2,− respectively. Note− that the Jacobian between (x , x , x ) and (x, y, z) 1 2 − 3 1 − 2 1 2 3 is 1, so is the Jacobian between (p1, p2, p3) and (p, r, q). The baryon wave function ϕB(y, z) is normalized as 3 3 2 d yd z ϕB(y, z) =1, whose Fourier transform is assumed to have the Gaussian form [47, 48] ´ | | ϕ (r, q) = d3yd3z exp( iq z ir y)ϕ (y, z) B ˆ − · − · B 3/2 2 2 3 1 r q = (2√π) exp 2 2 , (38) aB1aB2 −2aB1 − 2aB2 where aB1 and aB2 are two coeﬃcients in the Gaussian wave packet for the baryon. One can check the normalization 3 3 2 6 condition d rd q ϕB(r, q) = (2π) . The spin´ density| matrix| elements for the baryon are given by projecting the baryon states on the spin density operator,

ρ (p) = B; S,S , p ρ B; S,S , p Sz1Sz2 h z1 | | z2 i d3p d3p d3p = V 3 1 2 3 w (p )w (p )w (p ) ˆ (2π)3 (2π)3 (2π)3 q1,s1 1 q2,s2 2 q3,s3 3 s1X,s2,s3 q1X,q2,q3 B; S q , q , q ; s ,s ,s q , q , q ; s ,s ,s B; S × h z1| 1 2 3 1 2 3i h 1 2 3 1 2 3| z2i B; p q; p , p , p 2 . (39) × |h | 1 2 3i| Due to the factorization forms of the quark and baryon wave function in Eqs. (34,36), we can calculate the overlapping amplitudes of the spin-flavor parts and those of the momentum parts separately. We now look at the spin-flavor quark wave functions of spin-1/2 octet baryons and spin-3/2 decuplet baryons which belong to the 56-plet of SU(6) [61]. The spin-flavor quark wave functions of these baryons must be symmetric with repsect to any interchange of two quark labels since their color wave function is anti-symmetric and their orbital angular momentum is L =0. Therefore the SU(6) or spin-flavor quark wave functions are in the form [61]

1 1 B; S = ,S = (F χ + F χ ) 2 z √ MS MS MA MA 8 2 3 B; S = ,S = F χ (40) 2 z S S 10 where F denotes the flavor wave function and χ denotes the spin wave function with specific symmetry for interchange of two quark labels, the supscript 8 and 10 mean octet and decuplet respectively, and the subscript S, MS and MA mean symmetric, mixed symmetric and mixed anti-symmetric respectively. Then we compute the spin density matrix elements for baryons in Eq. (39). We take the Λ hyperon (B=Λ) as an example, whose flavor-spin parts of the spin density matrix elements can be evaluated by the SU(6) quark wave 8 function of Λ. The projection amplitude of a baryon momentum state on the three-quark one is given by

q; p , p , p B; p = d3x d3x d3x q; p , p , p x , x , x x , x , x B; p h 1 2 3| i ˆ 1 2 3 h 1 2 3| 1 2 3i h 1 2 3| i 1 = (2π)3δ(3)(p p p p )ϕ (r, q) , (41) V 2 − 1 − 2 − 3 B whose modulus square is (2π)3 q; p , p , p B; p 2 = δ(3)(p p p p ) ϕ (r, q) 2 . (42) |h 1 2 3| i| V 3 − 1 − 2 − 3 | B | Λ Λ Then we obtain the diagonal element ρ++ ( ρ 1 , 1 ) in Eq. (39), ≡ 2 2 1 1 1 1 ρΛ (p) = Λ; , , p ρ Λ; , , p ++ 2 2 2 2

1 d3r d 3q = ϕ (r, q) 2 24 ˆ (2π)3 (2π)3 | Λ | w (p ) [2 P (p )P (p ) P (p )P (p )] ×{ s,+ 1 − u 2 d 3 − u 3 d 2 +w (p ) [2 P (p )P (p ) P (p )P (p )] s,+ 2 − u 3 d 1 − u 1 d 3 +w (p ) [2 P (p )P (p ) P (p )P (p )] , (43) s,+ 3 − u 1 d 2 − u 2 d 1 } where we have applied p = p/3+ r/2+ q, p = p/3+ r/2 q and p = p/3 r. Note that the subscript ’s, +’ 1 2 − 3 − of ws,+ means the s quark with spin up (please do not confuse s with the quark spin). Similarly we obtain another Λ Λ diagonal element ρ−− ( ρ− 1 ,− 1 ) ≡ 2 2 1 1 1 1 ρΛ (p) = Λ; , , p ρ Λ; , , p −− 2 −2 2 −2

1 d3r d 3q = ϕ (r, q) 2 24 ˆ (2π)3 (2π)3 | Λ |

w −(p ) [2 P (p )P (p ) P (p )P (p )] ×{ s, 1 − u 2 d 3 − u 3 d 2 +w −(p ) [2 P (p )P (p ) P (p )P (p )] s, 2 − u 3 d 1 − u 1 d 3 +w −(p ) [2 P (p )P (p ) P (p )P (p )] . (44) s, 3 − u 1 d 2 − u 2 d 1 } All oﬀ-diagonal elements are vanishing, i.e. ρ+− = ρ−+ =0. The polarization formula for spin-1/2 baryons can be expressed as ρB (p) ρB (p) P ++ −− B,1/2(p) = B − B , (45) ρ++(p)+ ρ−−(p) where the superscript ’B’ denotes a type of spin-1/2 baryons. For the spin-3/2 baryons, the polarization is given by Eq. (A8), namely, 1 B p B p B p B p 3 ρ 1 , 1 ( ) ρ− 1 ,− 1 ( ) + ρ 3 , 3 ( ) ρ− 3 ,− 3 ( ) P 2 2 − 2 2 2 2 − 2 2 B,3/2(p)= h B B i B B , (46) ρ 1 1 (p)+ ρ 1 1 (p)+ ρ 3 3 (p)+ ρ 3 3 (p) 2 , 2 − 2 ,− 2 2 , 2 2 , 2 where the superscript ’B’ denotes a type of spin-3/2 baryons. For the Λ hyperon we have 1 d3r d3q ρΛ (p)+ ρΛ (p) = ϕ (r, q) 2 ++ −− 24 ˆ (2π)3 (2π)3 | Λ | [6 P (p )P (p ) P (p )P (p ) × − u 2 d 3 − u 3 d 2 P (p )P (p ) P (p )P (p ) − u 3 d 1 − u 1 d 3 P (p )P (p ) P (p )P (p )] , − u 1 d 2 − u 2 d 1 1 d3r d3q ρΛ (p) ρΛ (p) = ϕ (r, q) 2 ++ − −− 24 ˆ (2π)3 (2π)3 | Λ | P (p ) [2 P (p )P (p ) P (p )P (p )] ×{ s 1 − u 2 d 3 − u 3 d 2 +P (p ) [2 P (p )P (p ) P (p )P (p )] s 2 − u 3 d 1 − u 1 d 3 +P (p ) [2 P (p )P (p ) P (p )P (p )] . (47) s 3 − u 1 d 2 − u 2 d 1 } 9

Now we look at the limit that the quark polarization is much smaller than 1. In this case we can neglect the quardratic terms in the quark polarization, from Eq. (45) the polarization of Λ can be approximated as

1 d3r d3q P (p) ϕ (r, q) 2 [P (p )+ P (p )+ P (p )] , (48) Λ ≈ 3 ˆ (2π)3 (2π)3 | Λ | s 1 s 2 s 3 which shows that the polarization of Λ is mainly determined by that of the s-quark. For antihyperons, we can make the replacement Λ Λ¯ and u, d, s ¯u, d¯,¯s in Eqs. (43-48). Similarly with the→ SU(6) quark wave→ functions of the octet baryons Σ0 and Ξ− and the decuplet baryons ∆++ and Ω, we can obtain their polarizations from Eq. (45) and (46) respectively. If the quark polarization is independent of momentum, we obtain

PΛ0 = Ps, 2(Pu + Pd) Ps 3PuPdPs P 0 = − − , Σ 3 2P (P + P )+ P P − s u d u d P P P P2 4 s d 3 d s P − = − − , Ξ 3 4P P + P2 − d s s P (5 + P2) P = s s , Ω P2 3(1 + s ) P P2 u(5 + u) P ++ = . (49) ∆ P2 3(1 + u) If we further have P = P P , the above polarization results are consistent to Ref. [1]. u d ≡ q

V. FERMION POLARIZATIONS FROM WIGNER FUNCTIONS

The covariant Wigner function method [62–65] for spin-1/2 fermions is a useful tool to study the chiral magnetic effect, the chiral vortical effect and other related effects [26, 27, 51, 53–55]. The Wigner function is equivalent to the quantum field and contains all information that the quantum field does. One can also obtain the quark polarization density from the Wigner function. The Wigner function for spin-1/2 fermions in an external electromagnetic field is given by [62–65]

1 4 −ip·y ¯ y 1 1 y Wαβ(x, p)= 4 d ye ψβ x + U(A; x + y, x y)ψα x , (50) (2π) ˆ 2 2 − 2 − 2

x2 µ where ψ denotes the fermionic ﬁeld and α, β the spinor indices, U(A; x2, x1) = exp iQ x dx Aµ(x) is the gauge h 1 i link which renders the Wigner function gauge invariant with Aµ being the gauge potential´ of the electromagnetic ﬁeld, and Oˆ denotes the ensemble average of the operator Oˆ over thermal states. The spin tensor component of the WignerD functionE can be extracted as 1 M µαβ(x, p) Tr γµ,Sαβ W (x, p) ≡ 2 { } 1 = ǫµαβρA (x, p), (51) −2 ρ µ where Aρ(x, p) is the axial vector component of the Wigner function and can be obtained by A (x, p) = µ µ να µναλ 5 αβ i α β M µαβ Tr[γ γ5W (x, p)], here we have used γ ,S = ǫ γ γλ and S = 4 [γ ,γ ]. Note that (x, p) are real quantities. The spin angular momentum{ tensor} of a fermionic system is given by

M ρσ = d3x d3p dp u M λρσ(x, p) ˆ ˆ ˆ 0 λ 1 = u ǫλρσα d3x d3p dp A (x, p). (52) −2 λ ˆ ˆ ˆ 0 α

λρσ λ We see that dp0uλM (x, p) plays the role of the spin tensor density in phase space. Here u is the fluid velocity and p u ´p. A four-momentum can be decomposed as pλ =p ¯λ + (p u)uλ. So in a general frame the spatial 0 ≡ · · 10 momentum integral in Eq. (52) means d3p d3p¯. Then an on-shell momentum satisfies p2 = m2 which can be ≡ rewritten as p2 +¯p2 = m2. The on-shell energy is then given by E m2 p¯2. In the local rest frame of the fluid 0 p ≡ − we have uµ = (1, 0), p¯λ = (0, p), p¯2 = p 2, E = p 2 + m2. In thisp section we sometime use the covariant form −| | p | | and sometime boldface version of a four-vector. Byp boldface of the four-vector we refer to its the spatial component in the local rest frame of the fluid. The axial vector component at the zeroth or leading order in powers of the field strength tensor and space-time derivative is given by [27] A µ µ 5 (0) = Tr[γ γ W(0)] µ µ 2 2 = m [θ(p )n (p, n) θ( p )n ( p, n)] δ(p m )(f f−), (53) 0 − − 0 − − − + − where p u p and f± is given by 0 ≡ ·

2 0 0 f± θ(p )f (p µ±)+ θ( p )f ( p + µ±) . (54) ≡ (2π)3 FD 0 − − FD − 0 In Eqs. (53,54) upper/lower sign denotes the spin state along n (n is the spin quantization direction) in the ± βE particle’s rest frame, µ± denotes the corresponding chemical potentials, and fFD(E)=1/(e + 1) is the Fermi-Dirac distribution. In Eq. (53) nµ(p, n) is the spin quantization direction in the lab frame and given by

µ µ ν n p (n p)p n (p, n) = Λ ν ( vp)n (0, n)= · , n + · . (55) − m m(m + Ep)

µ ν Here Λ ν ( vp) is the Lorentz transformation for vp = p/Ep and n (0, n) = (0, n) is the four-vector of the spin quantization− direction in the particle’s rest frame. One can check that nµ(p, n) satisﬁes n2 = 1 and n p =0, so it behaves like a spin four-vector up to a factor of 1/2. − · The ﬁrst or next-to-leading order contribution for the axial vector component of the Wigner function for massive fermions can be obtained by generalizing the solution for massless fermions [26, 27, 51, 53],

2 2 α 1 ασ 2 2 d(f+ + f−) αλ δ(p m ) A (x, p) = ~βΩ˜ p δ(p m ) Q~F˜ p − (f + f−), (56) (1) −2 σ − d(βp ) − λ p2 m2 + 0 − ˜ρλ 1 ρλµν ˜ ξη 1 ξηνσ 1 µνσβ where F = 2 ǫ Fµν and Ω = 2 ǫ Ωνσ with Ωνσ = 2 (∂ν uσ ∂σuν ). Here ǫ and ǫµνσβ are anti- µνσβ − symmetric tensors with ǫ = 1( 1) and ǫµνσβ = 1(1) for even (odd) permutations of indices 0123, so we have 0123 − − ǫ = ǫ0123 =1. In Eq. (56) we have also assumed that β =1/T is a constant. − ρ We now compute dp0A (x, p) from Eqs. (53,56). At the leading order, we obtain ´ 1 m 1 1 α α p n α p n dp0A(0) = 3 n ( , ) β(E −µ ) n ( , ) β(E +µ ) ˆ (2π) Ep e p + +1 − − − e p + +1 1 1 nα(p, n) + nα( p, n) . (57) − eβ(Ep−µ−) +1 − − eβ(Ep+µ−) +1

α If µ± = µ does not depend on the spin state along n, we see from Eq. (57) that A(0) = 0. In this case the ± α non-vanishing polarization can only come from the ﬁrst-order contribution A(1)(x, p) in Eq. (56) whose vorticity part reads

α 1~ ˜ ασ d(f+ + f−) 2 2 dp0A(ω) = β dp0 Ω pσ δ(p m ) ˆ −2 ˆ d(βp0) − 1 1 eβ(Ep−µ) eβ(Ep+µ) = ~β Ω˜ ασ p p 3 σ p0=Ep β(E −µ) 2 σ p0=−Ep β(E +µ) 2 (2π) Ep | [e p + 1] − | [e p + 1] 1 1 eβ(Ep−µ) eβ(Ep+µ) ~β Ω˜ ασp + , (58) 3 σ β(E −µ) 2 β(E +µ) 2 → (2π) Ep [e p + 1] [e p + 1] where in the last line we have made the replacement p p in the antifermion term and implied the same on-shell σ →− momentum p = (Ep, p) for both fermions and antifermions. This makes no diﬀerence when carrying out the integral over momentum. We have assumed the chemical potentials are the same for both spin orientations, i.e. µ± = µ, so the sum over spin states gives a factor 2. Note that we can further simplify Eq. (58) by keeping the time-like 11

σ σ component of the momentum, p u Ep, in the last line if there is no dependence of fFD(Ep µ) on the direction → α ασ ∓ of p. This will result in that the polarization density is proportional to ω = Ω˜ uσ. The electromagnetic ﬁeld part α dp0A(EM) from Eq. (56) gives ´ 2 2 α αλ δ(p m ) dp A = Q~ dp F˜ p − (f + f−) ˆ 0 (EM) − ˆ 0 λ p2 m2 + − β(Ep−µ) β(Ep+µ) 1 1 αλ e e βQ~ F˜ pλ , (59) 3 2 β(Ep−µ) 2 β(Ep+µ) 2 → (2π) Ep [e + 1] − [e + 1]

λ where we have used the same on-shell momentum p = (Ep, p) for both fermions and antifermions. The derivation of (59) is given in Appendix B. λρσ We now obtain the spin tensor density in phase space, dp0uλM (x, p), in Eq. (52) as ´ 1 dp u M λρσ (x, p) = u ǫλρσα dp A (x, p), (60) ˆ 0 λ −2 λ ˆ 0 α where dp0Aα(x, p) are given by Eqs. (58,59). It is convenient to distinguish fermions from antifermions, for which we obtain´ their spin tensor densities in phase space

1 1 1 eβ(Ep∓µ) M ρσ(x, p) = ~β u ǫλρσα Ω˜ Q F˜ pσ , (61) ± 3 λ ασ ασ β(E ∓µ) 2 −2(2π) Ep ± Ep [e p + 1] where the upper/lower sign is for the fermion/antifermion. Using the number density for fermions and antifermions in phase space, 1 2 n±(x, p) = , (62) (2π)3 eβ(Ep∓µ) +1 we obtain the spin tensor per particle in phase space

ρσ ¯ ρσ M± (x, p) M± (x, p) = n±(x, p)

1 λρσα 1 σ = ~β uλǫ Ω˜ ασ Q F˜ασ p [1 fFD(Ep µ)] . (63) − 4Ep ± Ep − ∓

We then obtain the spin vector per particle in phase space for the massive fermion or antifermion according to the Pauli-Lubanski pseudovector, 1 Sµ (x, p) = ǫµρσν M¯ ± p ± −2m ρσ ν β µλ 1 µλ = ~ Ω˜ Q F˜ pλ [1 fFD(Ep µ)] , (64) 4m ± Ep − ∓ where p is an on-shell four-momentum for both fermion and antifermion and E m2 p¯2 = p u in a general p ≡ − · frame. p In the non-relativistic limit and the local rest frame of the ﬂuid, the spatial component of the spin per particle in momentum space takes the form β S±(p) ~ (E ω QB) [1 fFD(E µ)] , (65) ≈ 4m p ± − p ∓ where we have dropped the coordinate dependence of the spin vector for simplicity. The polarization per particle in momentum space can be obtained by

−→P±(p)=2S±(p). (66)

We can apply Eqs. (65,66) to quarks and antiquarks and use Eq. (30) to compute ρ00 for vector mesons and use Eq. (48) to compute the polarization of Λ (Λ¯). 12

VI. POLARIZATIONS: FROM QUARKS TO VECTOR MESONS AND BARYONS

In this section we look closely at the quark and hadron polarizations induced by the vorticity and magnetic field by applying Eqs. (65,66). We look at the contributions from the vorticity and from the magnetic field separately. For the simplicity of illustration we take a limit with three conditions: (a) small constant polarizations; (b) Boltzmann limit with 1 fFD(Ep µ) 1; (c) non-relativitic limit with Ep mq (mq is the quark mass). At these limits the magnitudes of− quark polarizations∓ ≈ have the simple form by Eqs. (65,66)≈ as 1 P (ω) βω, q ≈ 2 P (B) βµ B, (67) q ≈ mq where µmq = Qq/(2mq) is the magnetic moment of the quark with the electric charge Qq, and we have used ω = ω and B = B . Substituting Eq. (67) into Eq. (49), we obtain the polarizations of baryons to the leading order in the| | quark polarization.| | We can also obtain the spin density matrix elements and polarizations for the φ meson through Eqs. (31,33). All these results are listed in Table I. We make some remarks about these results. The first remark is about the spin density matrix element of the φ φ φ meson ρ00. The deviations from the non-polarized value, ρ00 1/3, have the opposite sign from the vorticity and φ 2 φ − 2 2 2 φ magnetic field: ρ00(ω) 1/3 (βω) /9 < 0 while ρ00(B) 1/3 4β µmsB /9 > 0. Our result for ρ00(ω) is consistent with the hadron− statitsical≈ − model [30] in which the− mesons≈ are treated as elementary particles, but our φ result for ρ00(B) is different from Eqs. (C5,C6) in the hadron statitsical model [30], see Appendix C for a brief introduction to the particle polarization in the hadron statitsical model from the magnetic field and vorticity. In general cases, ρ00 for other vector mesons is given by Eq. (31). Since the vorticity contribution ρ00(ω) is independent of quark flavors, so we always have 1 1 1 ρ (ω) (βω)2 < , (68) 00 ≈ 3 − 9 3 for all vector mesons. In contrast the magnetic field contribution ρ00(B) does depend on the electric charges of the quark and antiquark inside the vector meson, 1 4 ρ (B) β2µ µ B2 00 ≈ 3 − 9 mq1 m¯q2 1 1 Q Q = β2 1 2 B2, (69) 3 − 9 m1m2 where Q1 and Q2 are the electric charge of the quark and antiquark in the vector meson respectively, and m1 and 0 ∗0 m2 are the mass of the quark and antiquark respectively. So for electrically neutral vector mesons such as ρ , K , ¯ ∗0 S=1 + − ∗+ ∗− K and φ, we have ρ00 (B) > 1/3. But for electrically charged vector mesons such as ρ , ρ , K and K , we S=1 have ρ00 (B) < 1/3. So we may conclude that for electrically charged vector mesons we have ρ00 < 1/3, while the situation is inconclusive for electrically neutral vector mesons, i.e. the magnitude of ρ00 can be either ρ00 < 1/3 or ρ00 > 1/3 depending on the competition between the vorticity and magnetic field contribution. The second remark is about the polarizations of baryons in magnetic fields which we can express in terms of the baryon magnetic moments with following formula in the constituent quark model [61]:

µmΛ = µms, 1 µ = (2µ +2µ µ ), mΣ 3 mu md − ms 1 µ = (4µ µ ), mΞ 3 ms − md µm∆++ = 3µmu,

µmΩ = 3µms. (70) By these relations in the constituent quark model about baryon magnetic moments, our results for these baryons’ polarizations are consistent with the statistical model for hadrons [30] in which hadrons are treated as elementary particles, see Appendix C for details. Our results for baryon polarizations from the vorticity are also consistent with the hadron statistical model [30]. It is also interesting to look at the ratios among the polarizations of diﬀerent baryons: 1 1 1 5 5 P(ω) Λ:Σ0 : Ξ− : ∆++ :Ω= : : : : , → 2 2 2 6 6 0 − ++ 5 5 P(B) Λ:Σ : Ξ : w∆ :Ω= µ : µ 0 : µ − : µ ++ : µ , (71) → mΛ mΣ mΞ 9 m∆ 9 mΩ 13

Table I: The polarizations of the φ meson and baryons in the quark coalescence model. The spin density matrix elements φ ρ00 for the φ meson are also listed. In the fisrt/second line are listed the results for the φ meson and baryons from the vorticity/magnetic field. In the third line are listed the polarizations of baryons in terms of baryon magnetic moments given by the constituent quark model. These results are consistent with those from the hadron statistical model in which hadrons are treated as elementary particles. In the third line are also listed the results for the φ meson in the hadron statistical model in which φ is treated as an elementary particle. The results are all expanded to the leading order in the small polarization. φ Pφ ρ00 PΛ PΣ0 PΞ− P∆++ PΩ P φ 2 1 − 1 2 1 1 1 5 5 (ω) or ρ00(ω) 3 βω 3 9 (βω) 2 βω 2 βω 2 βω 6 βω 6 βω P φ 1 4 2 1 − 1 − 5 5 (B) or ρ00(B) 0 3 + 9 (βµmsB) βµmsB 3 βB(2µmu + 2µmd µms) 3 βB(4µms µmd) 3 βµmuB 3 βµmsB 2 1 1 2 5 5 P − − (B) 3 βµmφB 3 9 (βµmφB) βµmΛB βµmΣ0 B βµmΞ B 9 βµm∆++ B 9 βµmΩB where the first line is the ratios for the vorticity contributions while the second line is for the magnetic field contributions. The final remark is that the polarization of any hadron from the vorticity in the hadron statistical model is consistent with that in the quark coalescence model under the three conditions as listed in the beginning of this section.

VII. SUMMARY AND CONCLUSIONS

In this paper we formulate a non-relativistic quark coalesence model with explicit momentum dependence based on the spin density matrix, with which one can describe the spin alignments of vector mesons and polarizations of baryons in a uniform way. The building blocks of the quark coalescence model is the quark and antiquark polarizations as the functions of momenta. The quark and antiquark polarizations can be calculated from the spin-tensor component of the Wigner function. Then we obtain the quark and antiquark polarizations induced by the vorticity and the magnetic ﬁeld, with which the polarizations of vector mesons and baryons can be built up. For vector mesons we start from the density matrix of quark-antiquark states with propability functions related to the quark and antiquark polarizations. Then we project vector meson states onto the density matrix to obtain the spin density matrix element with momentum dependence. The overlapping amplitude between the vector meson state and the quark-antiquark state has to be evaluated. The ﬁnal result for the spin density matrix elements of vector mesons are obtained as a functional of the quark and antiquark polarization functions. In the same way, we can also describe the baryon polarizations in terms of the quark ones through the density matrix of three-quark states. By projection of the baryon state onto the density matrix of three-quark states we obtain the spin density matrix of baryons as a functional of the quark and antiquark polarization functions. The overlapping amplitude between the baryon state and the three-quark state has to be evaluated. From the spin density matrix of baryons we can compute the baryon polarizations. Note that the current quark coalescence model is non-relativistic and valid for mesons and baryons with small momenta compared to their masses. For hadrons with large momenta, one has to formulate a relativistic version of the model. The current model can be the basis for further numerical simulations with event generators to give realistic predictions about global polarizations of vector mesons and baryons, which we will investigate in the future.

Acknowledgments

QW thanks H. Li and X.-L. Xia for helpful discussions. QW is supported in part by the Major State Basic Research Development Program (973 program) in China under the Grant No. 2015CB856902 and 2014CB845402 and by the National Natural Science Foundation of China (NSFC) under the Grant No. 11535012.

Appendix A: Spin amplitudes and polarization vector of ∆++

In this appendix we give the spin amplitudes of the decay ∆++ p+ π+ and the polarization vector of ∆++. ++ ++ → 1 3 The initial spin states of ∆ are written as ∆ ,Sz , where Sz = 2 , 2 . The angular momentum states of the proton and pion are in the multiplet of total| angular momei ntum J =3± /2±which we denote as J, J . The state | zif J, J is a coupled state of the proton spin state 1 ,Sp and its orbital angular momentum state L,L with | zif 2 z p | ziL

L =1 and L =0, 1, z ± 3 3 1 1 , = 1, 1 , , 2 2 | iL 2 2 f p

3 3 1 1 , = 1, 1 , , 2 −2 | − iL 2 −2 f p

3 1 1 1 1 2 1 1 , = 1, 1 , + 1, 0 , , 2 2 r3 | iL 2 −2 r3 | iL 2 2 f p p

3 1 2 1 1 1 1 1 , = 1, 0 , + 1, 1 , . (A1) 2 −2 r3 | iL 2 −2 r3 | − iL 2 2 f p p

We can define the decay transition matrix S as 3 S ∆++,S = ,S , (A2) z 2 z f from which we can compute the decay transition amplitude bet ween the initial and final state, f(Sp,S ) = pˆ,Sp; π+ S ∆++,S z z h z | | zi 3 3 = pˆ,Sp; π+ ,S = pˆ,Sp ,S , (A3) h z 2 z h z 2 z f f p + p where the final state can be denoted as pˆ,Sz ; π θ,φ,S z with pˆ = (θ, φ ) being the proton’s momentum direction in the polar and azimuthal angle. With| pˆ L =1i,L ≡ | = Y i (θ, φ), we obtain f(Sp,S ) as h | ziL 1,Lz z z 1 3 1 3 f , = Y (θ, φ), f , =0, 2 2 1,1 −2 2 1 3 1 3 f , = Y − (θ, φ), f , =0, −2 −2 1, 1 2 −2 1 1 2 1 1 1 f , = Y (θ, φ), f , = Y (θ, φ), 2 2 r3 1,0 −2 2 r3 1,1 1 1 1 1 1 2 f , = Y − (θ, φ), f , = Y (θ, φ). (A4) 2 −2 r3 1, 1 −2 −2 r3 1,0 With the decay amplitudes (A4), we can evaluate Eq. (15) and obtain the angular distribution as

dN 3 2 2 = 1 ρ− 1 ,− 1 + ρ 1 1 1 2 ρ− 1 ,− 1 + ρ 1 , 1 cos θ dΩ 8π − 3 2 2 2 2 − − 2 2 2 2 2 + Reρ− 3 ,− 1 Reρ 1 , 3 sin(2θ)cos φ √3 2 2 − 2 2 2 + Imρ− 3 ,− 1 Imρ 1 , 3 sin(2θ) sin φ √3 2 2 − 2 2

2 2 Reρ− 3 , 1 + Reρ− 1 , 3 sin θ cos(2φ) −√3 2 2 2 2

2 2 Imρ− 3 , 1 + Imρ− 1 , 3 sin θ sin(2φ) , (A5) −√3 2 2 2 2 where the spin density matrix elements for ∆++ are deﬁned by ρ ∆++,S ρ ∆++,S . (A6) Sz1,Sz2 ≡ h z1| | z2i ++ The polarization vector −→P = (P1, P2, P3) for ∆ can be determined from Eq. (11) with S =3/2 and the spin operators for spin-3/2 particles being deﬁned by 0 √30 0 0 √30 0 3 00 0 − 2 1  √  i  √   1  ˆ 30 2 0 ˆ 3 0 2 0 ˆ 0 2 0 0 S1 = , S2 = − , S3 = 1 , (A7) 2  0 2 0 √3  2  0 2 0 √3   0 0 0     −   − 2   0 0 √3 0   0 0 √3 0   00 0 3       − 2  15 whose result is

P 1 2 4 1 = Reρ− 3 ,− 1 + Reρ 1 , 3 + Reρ− 1 , 1 , Tr(ρ) √3 2 2 2 2 3 2 2 P 1 2 4 2 = Imρ− 3 ,− 1 + Imρ 1 , 3 + Imρ− 1 , 1 , Tr(ρ) √3 2 2 2 2 3 2 2 P 1 1 3 = ρ 1 , 1 ρ− 1 ,− 1 + ρ 3 , 3 ρ− 3 ,− 3 . (A8) Tr(ρ) 3 2 2 − 2 2 2 2 − 2 2

Appendix B: Derivation of Eq. (59)

We give a detailed derivation of Eq. (59). To this end we look at the four-momentum integral,

4 µ 4 µν ′ 2 2 d pA = Q~ d pF˜ p δ (p m )(f + f−) ˆ (EM) ˆ ν − + 4 = Q~F˜µν d4pp δ′(p2 m2) (2π)3 ˆ ν − [θ(u p)f (u p µ)+ θ( u p)f ( u p + µ)] × · FD · − − · FD − · 4 = Q~F˜µν d4pp δ′(p2 m2)θ(u p)[f (u p µ) f (u p + µ)] (2π)3 ˆ ν − · FD · − − FD · 2 dδ(p2 m2) = Q~F˜µν d4p − θ(u p)[f (u p µ) f (u p + µ)] (2π)3 ˆ dpν · FD · − − FD · 2 d = Q~F˜µν d4pδ(p2 m2)θ(u p) [f (u p µ) f (u p + µ)] −(2π)3 ˆ − · dpν FD · − − FD · 2 = ~QβF˜µν u d4pδ(p2 m2)θ(u p) (2π)3 ν ˆ − · [f (u p µ)(1 f (u p µ)) f (u p + µ)(1 f (u p + µ))] × FD · − − FD · − − FD · − FD · 1 ~ ˜µν 3 pν = 3 QβF d p 2 (2π) ˆ Ep [f (E µ)(1 f (E µ)) f (E + µ)(1 f (E + µ))] . (B1) × FD p − − FD p − − FD p − FD p In the second line we have used µ± = µ, then the sum over spin states gives a factor 2. In the ﬁfth line we have used the integral by part and dropped the complete derivative term. In the last equality we have replaced u p /E ν → ν p and carried out the integral over p p u, where E = m2 p¯2. Then we arrive at Eq. (59). 0 ≡ · p − p

Appendix C: Particle polarization in statistical model for hadrons

In the statisitical model for hadrons, we treat a hadron as an elementary particle with inner structure. We take the nonrelativistic limit and consider a hadron with the spin S and the magnetic moment µm = µmS/S in a magnetic ﬁeld B = Bez along the z-direction. We know that the spin along the z-direction takes the value Sz = S, S + 1, ,S 1,S. Suppose that the ensemble of the particles is in an equilibrium state with the temperature−T , so− the ··· − probability for the state with a speciﬁc value of Sz is

eβµmBSz/S w(Sz)= . (C1) S βµ BS′ /S ′ e m z Sz=−S P Then the average polarization along the z-direction is

S 1 P = S w (S ) S S z z SzX=−S 1 1 1 βµ B = 1+ coth 1+ βµ B coth m . (C2) 2S 2S m − 2S 2S 16

If the magnetic ﬁeld is weak, in the leading order we have

(S + 1) P βµ B. (C3) S ≈ 3S m For spin-1/2 and 3/2 particles, the above becomes

5 P = βµ B, P = βµ B, (C4) S=1/2 m,1/2 S=3/2 9 m,3/2 which are consistent with the quark coalescence model. For spin-1 particles, we have 2 P = βµ B. (C5) S=1 3 m,1 We can also obtain the density matrix element in the weak magnetic ﬁeld

S=1 1 ρ00 = w(S =1,Sz =0)= 1 + 2cosh(βµm,1B) 1 1 1 (βµ B)2 , (C6) ≈ 3 − 9 m,1 ≤ 3 P 2 where the equality holds for particles without the magnetic moment. For the φ meson, we have φ = 3 βµmφB and ρ = 1 1 (βµ B)2. φ 3 − 9 mφ The polarization from the vorticity can be obatined by replacing µmB/S by ω.

[1] Z.-T. Liang and X.-N. Wang, Phys. Rev. Lett. 94, 102301 (2005), [Erratum: Phys. Rev. Lett.96,039901(2006)], nucl- th/0410079. [2] Z.-T. Liang and X.-N. Wang, Phys. Lett. B629, 20 (2005), nucl-th/0411101. [3] S. A. Voloshin (2004), nucl-th/0410089. [4] B. Betz, M. Gyulassy, and G. Torrieri, Phys. Rev. C76, 044901 (2007), 0708.0035. [5] F. Becattini, F. Piccinini, and J. Rizzo, Phys. Rev. C77, 024906 (2008), 0711.1253. [6] J.-H. Gao, S.-W. Chen, W.-t. Deng, Z.-T. Liang, Q. Wang, and X.-N. Wang, Phys. Rev. C77, 044902 (2008), 0710.2943. [7] Q. Wang, Nucl. Phys. A967, 225 (2017), 1704.04022. [8] L. Adamczyk et al. (STAR), Nature 548, 62 (2017), 1701.06657. [9] B. I. Abelev et al. (STAR), Phys. Rev. C76, 024915 (2007), [Erratum: Phys. Rev.C95,no.3,039906(2017)], 0705.1691. [10] D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, Nucl. Phys. A803, 227 (2008), 0711.0950. [11] V. Skokov, A. Yu. Illarionov, and V. Toneev, Int. J. Mod. Phys. A24, 5925 (2009), 0907.1396. [12] V. Voronyuk, V. D. Toneev, W. Cassing, E. L. Bratkovskaya, V. P. Konchakovski, and S. A. Voloshin, Phys. Rev. C83, 054911 (2011), 1103.4239. [13] W.-T. Deng and X.-G. Huang, Phys. Rev. C85, 044907 (2012), 1201.5108. [14] J. Bloczynski, X.-G. Huang, X. Zhang, and J. Liao, Phys. Lett. B718, 1529 (2013), 1209.6594. [15] L. McLerran and V. Skokov, Nucl. Phys. A929, 184 (2014), 1305.0774. [16] U. Gursoy, D. Kharzeev, and K. Rajagopal, Phys. Rev. C89, 054905 (2014), 1401.3805. [17] V. Roy and S. Pu, Phys. Rev. C92, 064902 (2015), 1508.03761. [18] K. Tuchin, Phys. Rev. C91, 064902 (2015), 1411.1363. [19] H. Li, X.-l. Sheng, and Q. Wang, Phys. Rev. C94, 044903 (2016), 1602.02223. [20] A. Vilenkin, Phys. Rev. D22, 3080 (1980). [21] K. Fukushima, D. E. Kharzeev, and H. J. Warringa, Phys. Rev. D78, 074033 (2008), 0808.3382. [22] D. T. Son and P. Surowka, Phys.Rev.Lett. 103, 191601 (2009), 0906.5044. [23] D. T. Son and B. Z. Spivak, Phys. Rev. B88, 104412 (2013), 1206.1627. [24] D. T. Son and N. Yamamoto, Phys. Rev. Lett. 109, 181602 (2012), 1203.2697. [25] M. A. Stephanov and Y. Yin, Phys. Rev. Lett. 109, 162001 (2012), 1207.0747. [26] J.-H. Gao, Z.-T. Liang, S. Pu, Q. Wang, and X.-N. Wang, Phys.Rev.Lett. 109, 232301 (2012), 1203.0725. [27] R.-h. Fang, L.-g. Pang, Q. Wang, and X.-n. Wang, Phys. Rev. C94, 024904 (2016), 1604.04036. [28] N. Mueller and R. Venugopalan (2017), 1702.01233. [29] F. Becattini, L. Csernai, and D. J. Wang, Phys. Rev. C88, 034905 (2013), 1304.4427. [30] F. Becattini, I. Karpenko, M. Lisa, I. Upsal, and S. Voloshin, Phys. Rev. C95, 054902 (2017), 1610.02506. [31] I. Karpenko and F. Becattini, Eur. Phys. J. p. 213 (2017), 1610.04717. 17

[32] H. Li, L.-G. Pang, Q. Wang, and X.-L. Xia (2017), 1704.01507. [33] Y. L. Xie, M. Bleicher, H. Stoecker, D. J. Wang, and L. P. Csernai, Phys. Rev. C94, 054907 (2016), 1610.08678. [34] Y. Xie, D. Wang, and L. P. Csernai, Phys. Rev. C95, 031901 (2017), 1703.03770. [35] C. van Weert, Ann. Phys. 140, 133 (1982). [36] D. Zubarev, A. Prozorkevich, and S. Smolyanskii, Teor. Mat. Fiz. 40, 394 (1979). [37] F. Becattini and L. Tinti, Annals Phys. 325, 1566 (2010), 0911.0864. [38] F. Becattini, Phys. Rev. Lett. 108, 244502 (2012), 1201.5278. [39] F. Becattini, V. Chandra, L. Del Zanna, and E. Grossi, Annals Phys. 338, 32 (2013), 1303.3431. [40] F. Becattini and E. Grossi, Phys. Rev. D92, 045037 (2015), 1505.07760. [41] T. Hayata, Y. Hidaka, T. Noumi, and M. Hongo, Phys. Rev. D92, 065008 (2015), 1503.04535. [42] L. P. Csernai, V. K. Magas, and D. J. Wang, Phys. Rev. C87, 034906 (2013), 1302.5310. [43] L. P. Csernai, D. J. Wang, M. Bleicher, and H. Stoecker, Phys. Rev. C90, 021904 (2014). [44] L.-G. Pang, H. Petersen, Q. Wang, and X.-N. Wang, Phys. Rev. Lett. 117, 192301 (2016), 1605.04024. [45] Y. Jiang, Z.-W. Lin, and J. Liao, Phys. Rev. C94, 044910 (2016), 1602.06580. [46] Y. Sun and C. M. Ko, Phys. Rev. C96, 024906 (2017), 1706.09467. [47] V. Greco, C. M. Ko, and P. Levai, Phys. Rev. Lett. 90, 202302 (2003), nucl-th/0301093. [48] R. J. Fries, B. Muller, C. Nonaka, and S. A. Bass, Phys. Rev. C68, 044902 (2003), nucl-th/0306027. [49] Q.-B. Xie and X.-M. Liu, Phys. Rev. D38, 2169 (1988). [50] R.-q. Wang, F.-l. Shao, J. Song, Q.-b. Xie, and Z.-t. Liang, Phys. Rev. C86, 054906 (2012), 1206.4373. [51] J.-W. Chen, S. Pu, Q. Wang, and X.-N. Wang, Phys. Rev. Lett. 110, 262301 (2013), 1210.8312. [52] D. T. Son and N. Yamamoto, Phys. Rev. D87, 085016 (2013), 1210.8158. [53] J.-h. Gao and Q. Wang, Phys. Lett. B749, 542 (2015), 1504.07334. [54] Y. Hidaka, S. Pu, and D.-L. Yang, Phys. Rev. D95, 091901 (2017), 1612.04630. [55] N. Mueller and R. Venugopalan (2017), 1701.03331. [56] P. Minnaert, Phys. Rev. 151, 1306 (1966). [57] G. Bunce et al., Phys. Rev. Lett. 36, 1113 (1976). [58] B. I. Abelev et al. (STAR), Phys. Rev. C77, 061902 (2008), 0801.1729. [59] B. Mohanty (ALICE), in 17th International Conference on Strangeness in Quark Mat- ter (SQM 2017) Utrecht, the Netherlands, July 10-15, 2017 (2017), 1711.02018, URL https://inspirehep.net/record/1634609/files/arXiv:1711.02018.pdf. [60] K. Schilling, P. Seyboth, and G. E. Wolf, Nucl. Phys. B15, 397 (1970), [Erratum: Nucl. Phys.B18,332(1970)]. [61] F. E. Close, An Introduction to Quarks and Partons (Academic Press, 1980), ISBN 012175152X. [62] U. W. Heinz, Phys. Rev. Lett. 51, 351 (1983). [63] H.-T. Elze, M. Gyulassy, and D. Vasak, Nucl.Phys. B276, 706 (1986). [64] D. Vasak, M. Gyulassy, and H.-T. Elze, Annals Phys. 173, 462 (1987). [65] A. Abada, M. C. Birse, P.-f. Zhuang, and U. W. Heinz, Phys. Rev. D54, 4175 (1996), hep-ph/9605249.