Covariant spin kinetic theory I: collisionless limit ^*

Kinetic theory is widely used to study transport phenomena in many-particle systems. The classical Boltzmann kinetic theory has been established as the framework to describe the evolution of the distribution function in phase space. To study the effects of spin, the quantum kinetic theory must be used, and a single distribution function is usually insufficient for such a purpose. For massless Dirac fermions, the leading spin effects appear at $O(\hbar)$; thus, two distribution functions are required, one for right-handed chirality and the other for left-handed chirality^①. This established framework is the chiral kinetic theory (CKT) [1-3], which has been intensively investigated recently [4-14]. The out-of-equilibrium dynamics of anomaly-induced phenomena, such as the chiral magnetic effect [15, 16] and chiral vortical effect [17-19], have also been thoroughly investigated in the framework of CKT.

Unlike the massless case, the spin of a massive fermion is independent of the kinetic momentum. As a result, the dynamic evolution of massive fermions is specified with more than two degrees of freedom. If $\theta^ \mu$ is the unit spacelike vector specifying the spin quantizing orientation, the dynamical variables are two parameters used to determine $\theta^ \mu$, and two distribution functions, $f_\pm$, for particles with spin parallel and anti-parallel to $\theta^ \mu$, respectively. Therefore, the kinetic theory of massive fermions is more complicated than the CKT, and extensive investigation of this framework is needed [20-23].

One of the motivations to develop the aforementioned kinetic theory is the spin-polarization phenomenon in heavy-ion collisions, which is an important probe of the hot and dense quark gluon matter [24-27]. The first signal of the global spin polarization of $\Lambda$ hyperons (hereafter, $\Lambda$ polarization) [28] indicates the existence of a very strong fluid vorticity [29-31]. The subsequent measurements exhibit nontrivial features that cannot be understood based on the simple vorticity interpretation of the spin polarization [32, 33]. For example, the measured longitudinal and transverse $\Lambda$ polarizations show the opposite azimuthal angle dependence compared with the so-called thermal vorticity [34-41]. This strongly indicates that the spin polarization has independent dynamic evolution, in a non-equilibrium state, rather than being chained to the fluid vorticity. A covariant kinetic theory for spin transport (hereafter, spin kinetic theory for short) would be a promising tool for capturing the dynamics of spin polarization.

In this paper, we derive the collisionless spin kinetic theory at $O(\hbar)$ in a curved background spacetime and an external electromagnetic field. As an application, we investigate spin polarization using the spin kinetic theory. We give the general expression for spin polarization in terms of $f_\pm$ and $\theta^ \mu$. We then specify the equilibrium conditions from the spin kinetic theory and derive the spin polarization at both local and global equilibria. We stress that the present study differs from earlier works [20-23] in the following aspects. First, we include curved geometric background spacetime as well as an electromagnetic field. Such a general formalism should be applicable to spin transport, not only in heavy-ion collisions and astrophysical systems, but also in deformed materials and thermal-gradient systems, which has attracted significant attention in condensed matter physics. Second, we show that the frame-choosing vector can always be eliminated in the covariant kinetic theory of massive fermions, unlike the massless case. Third, we discuss the underlying physics of the Clifford components and their constraint equations. Fourth, we provide the kinetic equations in a more transparent way, exhibiting clear physical contexts. In particular, we verify that, in the classical limit, these equations are correctly reduced to the Vlasov equation, Bargmann–Michel–Telegdi (BMT) equation [42], and Mathisson–Papapetrou–Dixon (MPD) equations [43-45]. Finally, we discuss the global equilibrium in terms of the spin vector $\theta^ \mu$. The validity of this equilibrium state is qualified by the resulting spin polarization, which is consistent with that in Refs. [46, 47].

This paper is organized as follows. In Sections 2 and 3, we introduce the Wigner function and discuss the physical meaning of the dynamic equation for each Clifford component of the Wigner function. In Section 4, we derive the semi-classical kinetic theory for massive fermions. In Section 5, we derive the kinetic representation of the spin polarization for both massive and massless fermions and investigate spin polarization at both local and global equilibria. In this paper, we adopt the same notations and conventions as those in Ref. [13]; for instance, $\nabla_ \mu$ denotes the covariant derivative in terms of diffeomorphism and the local Lorentz transformation, and $p_ \mu$ ( $y^ \mu$) is the momentum variable (its conjugate).

The Wigner operator covariant under the ${\rm{U}}(1)$ gauge, local Lorentz transformations, and diffeomorphism is defined as [13]

$$\hat{W}(x,p) = \int {\rm d}^4y\sqrt{-g(x)}\, {\rm e}^{-{\rm i}p\cdot y/\hbar} \hat{ \rho}(x,y) \,, \tag{ 1 }$$

$$\hat{ \rho}(x,y) = \bar\psi(x){\rm e}^{y\cdot \overleftarrow{D}/2} \otimes {\rm e}^{-y \cdot D/2} \psi(x) \,, \tag{ 2 }$$

where $\psi(x)$ is the Dirac spinor operator. Here, we introduce the following notations: $\bar\psi(x)\equiv \psi^\dagger(x)\gamma^ {\hat{0}}$ and $\bar\psi \overleftarrow{O}\equiv [O \psi]^\dagger \gamma^ {\hat{0}}$ for an operator O, and $[\bar\psi\otimes \psi]_{ab} = \bar\psi_b \psi_a$, with $a,b$ being the spinor indices. The derivative $D_\mu$ is called the horizontal lift of $\nabla_\mu$: $D_ \mu = \nabla_ \mu - \Gamma^{ \lambda}_{ \mu \nu}y^ \nu \partial_ \lambda^y$ in the tangent bundle [i.e., the $(x,y)$-space]. Similarly the horizontal lift in the cotangent bundle [i.e., the $(x,p)$-space] is given by

$$D_ \mu = \nabla_ \mu + \Gamma^{ \lambda}_{ \mu \nu}p_ \lambda \partial^ \nu_p \,. \tag{ 3 }$$

This $D_\mu$ gives us a great advantage regarding analysis because of the property $[D_\mu,y^\nu] = [D_\mu,p_\nu] = 0$. We note that the gauge field $A_ \mu$ should also be involved in $D_ \mu$ when acting on the Dirac spinor: $D_ \mu \psi(x,y) = ( \nabla_ \mu- \Gamma^{ \lambda}_{ \mu \nu}y^ \nu \partial_ \lambda^y+ {\rm i}A_ \mu/\hbar) \psi(x,y)$ with $\psi(x,y)\equiv {\rm e}^{y\cdot D}\psi(x)$.

The Wigner function is defined by replacing the operator $\hat{ \rho}(x,y)$ with the ensemble average $\rho(x,y)\equiv\langle\hat{ \rho}(x,y)\rangle$ in Eq. (1). In this paper, we focus on the collisionless limit; thus, we impose the spinor field to obey the Dirac equation $({\rm i}\hbar\gamma^\mu D_\mu-m)\psi(x) = \bar{\psi}(x)\,({\rm i}\hbar{ \overleftarrow D}_\mu\gamma^\mu +m) = 0$. In this case, we derive the kinetic theory of massive fermions in the same manner as that in Ref. [13] (in particular, see Section 3 and Appendices C, D, and E therein). After the semi-classical expansion^②, and the decomposition in terms of the Clifford algebra as $W = \dfrac{1}{4}[{\cal F}+{\rm i} \gamma^5{\cal P}+ \gamma^ \mu{\cal V}_ \mu + \gamma^5 \gamma^ \mu{\cal A}_ \mu+ \dfrac{1}{2} \sigma^{ \mu \nu}{\cal S}_{ \mu \nu}]$, we arrive at the following system of equations:

$$\Delta_ \mu {\cal V}^ \mu = \frac{\hbar^2}{24}( \nabla_\mu R_{\nu\rho}) \partial^\mu_p \partial^\nu_p {\cal V}^\rho \,, \tag{ 4 }$$

$$\hbar \Delta_ \mu {\cal A}^ \mu = - 2m {\cal P} \,, \tag{ 5 }$$

$$\hbar \Delta_{[ \mu} {\cal A}_{ \nu]} - \epsilon_{ \mu \nu \rho \sigma}{ \Pi}^ \rho {\cal V}^ \sigma = -\frac{\hbar^2}{8}\tilde{R}_{\rho\sigma\mu\nu}\partial_p^\rho {\cal V}^\sigma \,, \tag{ 6 }$$

$${ \Pi}_ \mu {\cal V}^ \mu = m {\cal F} + \frac{\hbar^2}{8}R_{ \mu \nu} \partial^ \mu_p {\cal V}^ \nu \,, \tag{ 7 }$$

$${ \Pi}_ \mu {\cal A}^ \mu = \frac{\hbar^2}{8}R_{ \mu \nu} \partial^ \mu_p {\cal A}^ \nu \,, \tag{ 8 }$$

$$\hbar \Delta_{[ \mu} {\cal V}_{ \nu]} - \epsilon_{ \mu \nu \rho \sigma}{ \Pi}^ \rho {\cal A}^ \sigma = m {\cal S}_{ \mu \nu} - \frac{\hbar^2}{8} \tilde{R}_{\rho\sigma\mu\nu} \partial^\rho_p {\cal A}^\sigma \,, \tag{ 9 }$$

$${ \Pi}_ \mu {\cal F}+\frac{\hbar}{2} \Delta^ \nu {\cal S}_{ \mu\nu} = m {\cal V}_ \mu \,, \tag{ 10 }$$

$$\frac{\hbar}{2} \Delta_ \mu {\cal P}- { \Pi}^ \nu \tilde{ {\cal S}}_{\mu\nu} = m {\cal A}_ \mu - \frac{\hbar^2}{8} \tilde{R}_{\nu\rho\sigma\mu}\partial_p^\nu {\cal S}^{\rho\sigma} \,, \tag{ 11 }$$

$${ \Pi}_ \mu {\cal P} + \frac{\hbar}{2}\Delta^\nu \tilde {\cal S}_{\mu\nu} = 0 \,, \tag{ 12 }$$

$$\frac{\hbar}{2} \Delta_ \mu {\cal F}-{ \Pi}^ \nu {\cal S}_{ \mu \nu} = -\frac{\hbar^2}{16} (R_{ \mu \nu \rho\sigma} \partial^ \nu_p {\cal S}^{ \rho\sigma} + 2 R^{\nu\rho} \partial^p_ \nu {\cal S}_{ \rho \mu}) \quad\ \ \tag{ 13 }$$

with $X_{[\mu}Y_{\nu]} = \frac{1}{2}(X_\mu Y_\nu - X_\nu Y_\mu)$. Here, ${R^\rho}_{\sigma\mu\nu} = 2\partial_{[\nu}\Gamma_{\mu]\sigma}^\rho + 2\Gamma_{\lambda[\nu}^\rho\Gamma^\lambda_{\mu]\sigma}$ is the Riemann tensor, $R_{\mu\nu} = {R^\rho}_{\mu\rho\nu}$ is the Ricci tensor, and we define $\tilde{ {\cal S}}^{\mu\nu} = \dfrac{1}{2}\epsilon^{\mu\nu\rho\sigma} {\cal S}_{\rho\sigma}$ and $\tilde{R}^{\mu\nu\alpha\beta} = \dfrac{1}{2}\epsilon^{\alpha\beta\rho\sigma}{R^{\mu\nu}}_{\rho\sigma}$. The operators $\Pi_\mu$ and $\Delta_\mu$ are given by

$$\begin{split} { \Pi}_\mu = & p_\mu - \frac{\hbar^2}{12}(\nabla_\rho F_{\mu\nu})\partial^ \nu_p\partial^\rho_p + \frac{\hbar^2}{24}{R^\rho}_{\sigma\mu\nu}\partial^\sigma_p\partial_p^\nu p_\rho \\ & +\frac{\hbar^2}{4}R_{\mu\nu}\partial_p^\nu, \\ \Delta_\mu = & D_\mu -F_{\mu\lambda}\partial_p^\lambda - \frac{\hbar^2}{12}(\nabla_ \rho R_{ \mu \nu})\partial_p^ \rho\partial_p^\nu \\ & - \frac{\hbar^2}{24}(\nabla_ \lambda{R^\rho}_{\sigma\mu\nu}) \partial_p^\nu\partial_p^ \sigma\partial_p^\lambda p_\rho + \frac{\hbar^2}{8}{R^\rho}_{\sigma\mu\nu}\partial_p^\nu\partial_p^\sigma D_\rho \\ & + \frac{\hbar^2}{24} ( \nabla_\alpha\nabla_\beta F_{\mu\nu} + 2{R^\rho}_{\alpha\mu\nu}F_{\beta\rho} ) \partial_p^\nu \partial_p^\alpha\partial_p^\beta \,. \end{split} \tag{ 14 }$$

In equations (4)-(13), the spacetime curvature enters at least at $O(\hbar^2)$. We note that the Clifford coefficients ${\cal F}$, ${\cal P}$, ${\cal V}_ \mu$, ${\cal A}_ \mu$, and ${\cal S}_{ \mu \nu}$ are not totally independent. To proceed, we choose ${\cal V}_ \mu$ and ${\cal A}_ \mu$ as the independent variables^③. Thus, ${\cal P}$, ${\cal F}$, and ${\cal S}_{ \mu \nu}$ can be expressed in terms of ${\cal V}_ \mu$ and ${\cal A}_ \mu$ using Eqs. (5), (7), and (9). In Minkowski spacetime, the same set of equations up to $O(\hbar)$ was first derived in Ref. [48].

In the kinetic description, various physical quantities are built from W, which is (the Wigner transformation of) a two-point correlator of Dirac fields. For instance, the vector and axial current are computed as

$$J^\mu = \int_p {\rm{tr}}\Bigl[\gamma^\mu W\Bigr] = \int_p {\cal V}^\mu \,, \tag{ 15 }$$

$$J_{5}^\mu = \int_p {\rm{tr}}\Bigl[\gamma^\mu \gamma^5 W\Bigr] = \int_p {\cal A}^\mu \tag{ 16 }$$

with $\int_p \equiv\int \dfrac{{\rm d}^4 p}{(2 \pi)^4 [-g(x)]^{1/2}}$. From these, the Clifford coefficients ${\cal V}_ \mu$ and ${\cal A}_ \mu$ can be identified as the corresponding current densities in phase space (see more discussion in Sec. 5). In a similar way, ${\cal F}$ is the scalar condensate density (which in the classical limit is also interpreted as the distribution function of the vector charge); ${\cal P}$ is the axial condensate density; and ${\cal S}_{ \mu \nu}$ is the electromagnetic dipole moment density, up to a factor of m [see Eqs. (50)-(54)]. For convenience, we further represent the canonical energy-momentum tensor, spin current, and total angular momentum current as

$$T^{\mu\nu} = \int_p {\rm{tr}}\biggl[ \frac{i\hbar}{2}\gamma^\mu \overset{\leftrightarrow\;\;}{D^\nu} W ] = \int_p {\cal V}^ \mu p^ \nu \,, \tag{ 17 }$$

$$S^{\mu,\nu\rho} = \int_p {\rm{tr}}\biggl[ \frac{\hbar}{4}\{\gamma^\mu,\sigma^{\nu\rho}\} W ] = -\frac{\hbar}{2}\int_p \epsilon^{\mu\nu\rho\lambda} {\cal A}_ \lambda \,, \tag{ 18 }$$

$$M^{\lambda,\mu\nu} = x^\mu T^{\lambda\nu} - x^\nu T^{\lambda\mu} + S^{\lambda,\mu\nu} \,, \tag{ 19 }$$

where we define $\overset{\leftrightarrow}{D}_\mu (\bar\psi_b \otimes \psi_a) \equiv \bar\psi_b \otimes D_\mu \psi_a - \bar\psi_b \overleftarrow{D}_\mu \otimes \psi_a$.

Note that $S^{\lambda,\mu\nu}$ is not anticipated to be an observable for spin because there exists the Belinfante-Rosenfeld type pseudo-gauge ambiguity [49]; $S^{\lambda,\mu\nu}$ in Eq. (19) can be absorbed into a redefinition of the energy momentum tensor once the Belinfante tensor has been introduced. Instead, an unambiguous way to represent particle spin is to employ the Pauli-Lubanski (PL) vector operator:

$$\hat{ {\cal W}}^{ \mu} \equiv -\frac{1}{\hbar} \epsilon^{ \mu \nu \rho \sigma}\hat{P}_ \nu \hat{M}_{ \rho \sigma} \,, \tag{ 20 }$$

where the hat symbols denote quantum mechanical operators, $\hat{P}_ \nu$ and $\hat{M}_{ \rho \sigma}$ are the canonical momentum and total angular momentum operators, respectively, and the prefactor $-1/\hbar$ is introduced as our convention. It is important to note that the orbital part of the canonical angular momentum does not contribute to the above equation.

Following Eq. (20), we may define the PL vector in our kinetic theory as [50]^④

$${\cal W}^\mu (x,p) \equiv -\frac{1}{\hbar(p\cdot \nu)} \epsilon^{ \mu \nu \rho \sigma} p_\nu M_{\rho\sigma} \,, \tag{ 21 }$$

where we define $M_{\rho\sigma}\equiv\nu^\lambda M_{\lambda,\rho\sigma}$ with $\nu^\mu$ being a unit timelike vector, and the factor $1/{p\cdot \nu}$ is introduced for normalization. It can be easily checked that the above definition of ${\cal W}^\mu$ excludes the orbital angular momentum part and can be reduced as

$${\cal W}^\mu (x,p) = {\cal A}^\mu (x,p) \,. \tag{ 22 }$$

Note that this relation is independent of the vector $\nu^\mu$. The coincidence between ${\cal W}^\mu$ and ${\cal A}^\mu$ is expected. As an example, for massless fermions, the magnetic-field-induced spin polarization can be considered as the axial current generation, which is the chiral separation effect [51, 52]. The spin polarization density defined with the PL vector is hence equivalent to the axial current:

$$\mathscr{W}^{ \mu}(x) \equiv\int_p {\cal W}^ \mu(x,p) = \int_p {\cal A}^\mu(x,p) = J_5^\mu(x) \,. \tag{ 23 }$$

In this section, we discuss the physical meanings of Eqs. (4)-(13). To show this, we perform integration over momentum space, which results in much simpler expressions; most of the total-derivative terms vanish as surface integrals [an exception in Eq. (5) is discussed later]. For simplicity, we hereafter only focus on $O(\hbar)$ terms, so that the Riemann curvature is neglected in Eqs. (4)-(13) and Eq. (14) is reduced to $\Pi_\mu = p_\mu$ and $\Delta_\mu = D_\mu - F_{\mu\lambda}\partial_p^\lambda$.

First, we demonstrate that Eqs. (4) and (6) lead to fundamental Ward identities. After integrating over momentum space, Eq. (4) yields the vector current conservation law:

$$\nabla_ \mu J^ \mu = 0 \,. \tag{ 24 }$$

It is obvious from this that Eq. (4) is the kinetic equation of the vector charged particle. Integrating Eq. (4) after multiplying by $p^ \nu$, we obtain the energy-momentum conservation law in the presence of an external field:

$$\nabla_ \mu \left( T^{ \mu \nu} + T_{\rm{ext}}^{ \mu \nu} \right) = 0 \,, \tag{ 25 }$$

where $T_{\rm{ext}}^{ \mu \nu} = -F^{ \mu \lambda}F^ \nu_{\;\; \lambda} + \frac{1}{4} g^{ \mu \nu} F^{ \rho \sigma}F_{ \rho \sigma}$ is the energy-momentum tensor of an electromagnetic field. Here, we have used Maxwell's equation $\nabla_ \mu F^{ \mu \nu} = J^ \nu$ and the Bianchi identity $\nabla_ \mu \tilde{F}^{ \mu \nu} = 0$ with $\tilde{F}^{ \mu \nu} = \frac{1}{2} \epsilon^{ \mu \nu \rho \sigma}F_{ \rho \sigma}$. On integrating Eq. (6) over momentum space, we obtain

$$\nabla_ \mu S^{ \mu, \rho \sigma} = T^{ \sigma \rho}-T^{ \rho \sigma} \,. \tag{ 26 }$$

This, combined with Eq. (25), gives the conservation law of the canonical total angular momentum:

$$\nabla_ \mu \left( M^{ \mu, \rho \sigma} + M^{ \mu, \rho \sigma}_{\rm{ext}} \right) = 0 \tag{ 27 }$$

with $M^{\lambda,\mu\nu}_{\rm{ext}} = x^\mu T^{\lambda\nu}_{\rm{ext}} - x^ \nu T^{\lambda\mu}_{\rm{ext}}$ being the angular momentum of electromagnetic field. This reflects the absence of the Lorentz anomaly [53].

Next, we consider Eqs. (8)-(13), the physical contents of which are less transparent than those of Eqs. (4) and (6). Equation (8) involves only ${\cal A}^ \mu$; thus, it is a subsidiary condition for ${\cal A}^ \mu$. Up to $O(\hbar)$, it reduces to

$$p_ \mu {\cal A}^ \mu = 0 \,. \tag{ 28 }$$

Based on the identification (22), ${\cal W}^\mu = {\cal A}^\mu$, the above equation implies the following facts: spin must be either perpendicular to the momentum (i.e., for massive fermions) or parallel to the momentum (i.e., for massless fermions so that $p^2 = 0$ on-shell classically). In Section 4, we discuss the details with quantum corrections. The electromagnetic dipole moment is derived from Eq. (9):

$$m\int_p {\cal S}_{ \mu \nu} = -\int_p \epsilon_{ \mu \nu \rho \sigma}p^ \rho {\cal A}^ \sigma+\hbar \nabla_{[ \mu} J_{ \nu]} \,, \tag{ 29 }$$

where the first (second) term on the right-hand side represents the spin (orbital) contribution. Equations (10) and (11) are Gordon decompositions for the vector and axial currents. Upon integration over momentum, they separate the convection and gradient currents:

$$mJ^ \mu =\int_p p^ \mu {\cal F}+\frac{\hbar}{2} \nabla_ \nu \int_p {\cal S}^{ \mu \nu} \,, \tag{ 30 }$$

$$m J_5^{ \mu} = - \int_p p_ \nu \tilde{ {\cal S}}^{ \mu \nu}+\frac{\hbar}{2} \nabla^ \mu \int_p {\cal P} \,. \tag{ 31 }$$

We note that the second term in Eq. (30) is the covariant form of the well-known magnetization current. Similarly, Eqs. (12) and (13) give

$$0 =\int_p p^ \mu {\cal P}+\frac{\hbar}{2} \nabla_ \nu \int_p \tilde{ {\cal S}}^{ \mu \nu} \,, \tag{ 32 }$$

$$0 = -\int_p p^ \nu {\cal S}_{ \mu \nu}+ \frac{\hbar}{2} \nabla^ \mu \int_p {\cal F} \,, \tag{ 33 }$$

where the right-hand sides are dual to those of Eqs. (30) and (31). We note that ${\cal P}$ is regarded as the source of spin [see Eq. (5)]. Thus, Eqs. (32) and (33) imply that there do not exist vector and axial currents carrying 'magnetic charges' in Dirac theory.

Finally, we consider quantum anomalies related to Eqs. (5) and (7). The momentum integral of Eq. (5) generates a nonvanishing surface term. After a technical evaluation of such a term, we derive the anomalous axial Ward identity:

$$\nabla_ \mu J_{5}^ \mu = \mathscr{A} -\frac{2m}{\hbar}\int_p {\cal P}, \tag{ 34 }$$

where $\mathscr{A}$ is the chiral anomaly originating from the surface integral (see Appendix A). Similarly, from Eq. (7), we obtain

$${T^\mu}_\mu = m \int_p {\cal F} \,, \tag{ 35 }$$

which represents the Ward identity in terms of the dilatation. We emphasize that up to $O(\hbar)$, no surface integral contributes to Eq. (35). As a result, the trace anomaly does not emerge here, while the chiral anomaly does, as given in Eq. (34). Indeed, one can confirm from dimensional analysis that the trace anomaly is $O(\hbar)$ higher than the chiral anomaly^④. For the same reason, the chiral anomaly in Eq. (34) is not involved in the gravitational contribution, which is $O(\hbar^2)$ higher than the electromagnetic one [13]. In the kinetic theory involving $O(\hbar^2)$ or $O(\hbar^3)$ terms, these additional contributions should enter the right-hand sides of Eqs. (34) and (35). We will leave discussion of the higher-order kinetic theory to a future publication.

In the kinetic theory up to $O(\hbar)$, the general solutions for ${\cal V}_ \mu$ and ${\cal A}_ \mu$ take the following forms (see Appendix B):

$$\begin{split} {\cal V}^ \mu = & 4 \pi \biggl[ \delta(\xi) \biggl( p^ \mu f + \frac{\hbar}{2p\cdot n}\epsilon^{ \mu \nu \rho \sigma} n_ \nu \Delta_ \rho \bar{ {\cal A}}_{ \sigma} ) \\ & + \delta'(\xi) \hbar \tilde{F}^{ \mu \nu} \biggl( \bar{ {\cal A}}_{ \nu} {-\frac{p\cdot\bar{ {\cal A}}}{p\cdot n} n_ \nu} ) ]\,, \end{split} \tag{ 36 }$$

$${\cal A}^ \mu = 4 \pi \biggl[ \delta(\xi) \bar{ {\cal A}}^{ \mu} + \delta'(\xi)\hbar \tilde{F}^{ \mu \nu}p_ \nu f ],\ \tag{ 37 }$$

where we utilize $x\delta'(x) = -\delta(x)$, and denote $\Delta_ \mu = D_ \mu - F_{ \mu \lambda}\partial_p^\lambda$ and

$$\xi = p^2 - m^2 \,. \tag{ 38 }$$

The scalar function $f = f(x,p)$ is to be considered as the distribution function of vector charge, and $n^ \mu$ is a unit timelike vector that satisfies $p\cdot n\neq0$. According to Eq. (8), the vector $\bar{ {\cal A}}^{ \mu}$ must satisfy the condition

$$p_ \mu \bar{ {\cal A}}^{ \mu} \delta(\xi) = 0 \,. \tag{ 39 }$$

Here, $\bar{ {\cal A}}^{ \mu}$ is not necessarily perpendicular to the momentum because of the presence of the delta function. To proceed, we decompose $\bar{ {\cal A}}^{ \mu}$ as

$$\bar{ {\cal A}}^ \mu = p_ \mu f_5 + \bar{ {\cal A}}_{\perp}^ \mu \,, \tag{ 40 }$$

where $\bar{ {\cal A}}_{\perp}^ \mu$ is perpendicular to the momentum: $p\cdot\bar{ {\cal A}}_\perp = 0$.

4.1. Massless case

In the massless limit, plugging Eqs. (36) and (37) into Eq. (9), we identify

$$\bar{ {\cal A}}_{\perp}^ \mu = \hbar \Sigma_n^{ \mu \nu} \Delta_ \nu f \,, \quad \Sigma_n^{ \mu \nu}\equiv\frac{ \epsilon^{ \mu \nu \rho \sigma}p_ \rho n_ \sigma}{2p\cdot n} \,. \tag{ 41 }$$

Then, the solutions (36) and (37) are reduced to

$$\begin{split} ( {\cal V}, {\cal A})^ \mu = & 4 \pi \biggl[ \delta(p^2) \Bigl\{ p^ \mu \left( f,f_5 \right) +\hbar \Sigma_n^{ \mu \nu} \Delta_ \nu \left( f_5,f \right) \Bigr\} \\ & + \hbar \delta'(p^2) \tilde{F}^{ \mu \nu}p_ \nu \left( f_5,f \right) ] \,. \end{split} \tag{ 42 }$$

This indicates that $f_5$ plays the role of the axial charge distribution function. The second term in the above equation is called the side-jump term (or the magnetization current), and $\Sigma_n^{ \mu \nu}$ is known as the spin tensor at the spin-defining vector $n^\mu$ [6, 7]; e.g., $\Sigma^{ij}_n = \epsilon^{ijk0}p_k/2p_0$ in the rest frame $n^\mu = (1, {\bf 0})$. Additionally, it is important to note from the above ${\cal A}^\mu$ that $\Sigma_n^{\mu\nu}$ is connected to the canonical spin current (18) through

$$S^{ \lambda, \mu \nu}n_ \lambda = \hbar \int_p 4 \pi \delta(p^2)\, p\cdot n\, f_5 \Sigma_n^{ \mu \nu}\,. \tag{ 43 }$$

This relation more transparently accounts for why $\Sigma^{\mu\nu}_n$ characterizes the particle spin, and $n^ \mu$ represents the frame of the spin^⑤.

Note that Eq. (42) correctly reproduces the solution in the CKT, with the replacement as ${\cal R}^\mu/ {\cal L}^\mu = \frac{1}{2}( {\cal V}\pm {\cal A})^\mu$ and $f_{R/L} = \frac{1}{2}(f\pm f_5)$ [13]. Accordingly, the chiral kinetic equations are also obtained from Eqs. (4) and (5) with the above solutions (42):

$$\begin{split} 0 & = \delta(p^2\mp\hbar F_{ \alpha \beta} \Sigma_n^{ \alpha \beta}) \biggl[ p_ \mu \Delta^ \mu f_{R/L} \\ & \quad \pm\frac{\hbar}{p\cdot n}\tilde{F}_{ \mu \nu} n^ \mu \Delta^ \nu f_{R/L} \pm \hbar \Delta^ \mu \left( \Sigma^n_{ \mu \nu} \Delta^ \nu f_{R/L} \right) ] \,. \end{split} \tag{ 44 }$$

More discussion about the CKT can be found, e.g., in Refs. [8, 11, 13].

Now, we re-consider the chiral anomaly in the CKT. Using the $O(\hbar)$ solution ${\cal A}^ \mu$ in Eq. (42), we derive the anomalous Ward identity

$$\begin{split} & \qquad\quad \nabla_ \mu J^ \mu_{5} = \mathscr{A} \,,\\ \mathscr{A} & = -\frac{\hbar}{8} F^{\mu\nu} \tilde{F}_{\mu\nu} \int_ { p} f( { p}) \partial_p^i \frac{p_i}{| { p}|^3} \\ & = -\frac{\hbar}{16\pi^2} F^{\mu\nu} \tilde{F}_{\mu\nu} f( {\bf 0}) \,, \end{split} \tag{ 45 }$$

with $\int_ { p} \equiv\int \dfrac{{\rm d}^3 p}{(2 \pi)^3}$. This reproduces the usual chiral anomaly when f is (twice that of) the Fermi-Dirac distribution (see details in Appendix A). The important fact is that $\mathscr{A}$ receives the contribution only from the singular term at $p^2 = 0$, which generates the Berry monopole $\partial_p^i \dfrac{p_i}{| { p}|^3} = 4\pi\delta^3( { p})$. The above covariant expression hence manifests that the chiral anomaly is a topological nature of massless fermions in an electromagnetic field [1, 2].

4.2. Massive case

We now focus on the massive case, in which we can perform two reductions for the solutions (36) and (37). First, Eq. (39) for $m\neq 0$ implies

$$f_5 \delta(\xi) = 0 \,, \tag{ 46 }$$

with which one can remove the parallel part in $\bar{ {\cal A}}^ \mu$ from the solutions (36) and (37). Second, the frame vector $n^\mu$ in Eqs. (36) and (37) can be absorbed into the distribution function f, redefined as (see Appendix C):

$$f \rightarrow f + \frac{\hbar}{m^2}\frac{\epsilon_{ \mu \nu \rho \sigma} p^ \mu n^ \nu}{2 p\cdot n} \Delta^ \rho \bar{ {\cal A}}_{\perp}^\sigma \,. \tag{ 47 }$$

We emphasize that this redefinition is equivalent to identifying the frame $n^ \mu$ as the particle's rest frame $n_{\rm{rest}}^ \mu = p^ \mu/m$. The frame vector $n^ \mu$ can be removed because, in the massive case, there is a special choice of $n^\mu$, i.e., the rest frame $n_{\rm{rest}}^ \mu$. Thus, we can always redefine the scalar distribution function f from $n^ \mu$ to $n_{\rm{rest}}^ \mu$ through a local Lorentz boost. This procedure does not work for massless fermions because of the lack of such a special frame, making it inevitable to introduce $n^ \mu$.

Because of the constraint $p\cdot\bar{ {\cal A}}_\perp = 0$, there are three degrees of freedom in $\bar{ {\cal A}}_\perp^\mu$. One of them is interpreted as the axial charge distribution, which specifies the norm of $\bar{ {\cal A}}_\perp^\mu$. Because the axial current density is identified as the particle spin, as shown in Eq. (22), the other two correspond to the parameters of the spin direction. In the massive case, we hence parametrize $\bar {\cal A}^\mu_\perp$ as

$$\bar{ {\cal A}}_{\perp}^{ \mu} \equiv m \theta^ \mu f_A \,, \tag{ 48 }$$

where $f_A$ is the axial distribution function and $\theta^\mu$ is a unit vector (with two degrees of freedom). Note that $\theta^\mu$ is normalized with the spacelike condition $\theta^ \mu \theta_ \mu = -1$ and $p_ \mu \theta^ \mu = 0$. In addition, it is useful for later discussion to introduce the following tensor:

$$\Sigma_S^{ \mu \nu}\equiv\frac{1}{2m} \epsilon^{ \mu \nu \rho \sigma} \theta_ \rho p_ \sigma \,, \tag{ 49 }$$

which may be regarded as the spin tensor of massive fermions, as $\Sigma_n^{\mu\nu}$ is for massless fermions (41). Indeed, it is easily verified that $\Sigma^{ij}_S = \epsilon^{ijk0}\theta_k/2$ for the rest particle with $p_\mu = (m, {\bf 0})$.

Collecting the discussions above, we present the solutions of the Clifford coefficients ${\cal F}$, ${\cal P}$, ${\cal V}_ \mu$, ${\cal A}_ \mu$, and ${\cal S}_{ \mu \nu}$, as follows:

$${\cal F} = 4 \pi \biggl[ \delta(\xi) m f - \delta'(\xi)\hbar m F^{ \mu \nu} \Sigma^S_{ \mu \nu} f_A ], \tag{ 50 }$$

$${\cal P} ={-2\pi\hbar \Delta_ \mu \left[ \theta^ \mu f_A \delta(\xi) \right]}, \tag{ 51 }$$

$$\begin{split} \quad\quad {\cal V}^ \mu = & 4 \pi \biggl[ \delta(\xi) \biggl( p^ \mu f +\hbar\frac{\epsilon^{ \mu \nu \rho \sigma}}{2m}p_ \nu \Delta_ \rho \left( \theta_ \sigma f_A \right) ) \\ & + \delta'(\xi)\hbar\,m\tilde{F}^{ \mu \nu} \theta_ \nu f_A ], \end{split} \tag{ 52 }$$

$${\cal A}^ \mu = 4 \pi \biggl[ \delta(\xi) m \theta^ \mu f_A + \delta'(\xi)\hbar \tilde{F}^{ \mu \nu}p_ \nu f ], \tag{ 53 }$$

$$\begin{split} \quad\quad {\cal S}^{ \mu \nu} = & 4 \pi \biggl[ \delta(\xi) \biggl( 2 m f_A \Sigma_S^{ \mu \nu} {-\frac{\hbar}{m}p^{[ \mu} \Delta^{ \nu]}f} ) \\ & - \delta'(\xi) \hbar m F^{ \mu \nu} f ] \end{split} \tag{ 54 }$$

with $\xi = p^2 -m^2$. In Eqs. (50)-(54), there are four independent variables: two for the distribution functions f and $f_A$, and the other two for the spatial orientation of the spin vector $\theta^ \mu$. Therefore, the covariant spin kinetic theory up to $O(\hbar)$ is described by the following four independent evolution equations:

$$\begin{split} 0 = & \delta(\xi\mp \hbar \Sigma_S^{ \alpha \beta} F_{ \alpha \beta}) \\ & \times \biggl[ \biggl( p\cdot \Delta \pm\frac{\hbar}{2} \Sigma_S^{ \mu \nu} \left( \nabla_ \rho F_{ \mu \nu} -p_ \lambda{R^ \lambda}_{ \rho \mu \nu} \right) \partial^ \rho_p ) f_{\pm} \\ & +\frac{\hbar}{2}(f_+ - f_-)\biggl( \left( \nabla_ \rho F_{ \mu \nu} -p_ \lambda{R^ \lambda}_{ \rho \mu \nu} \right) \partial^ \rho_p \Sigma_S^{ \mu \nu}\\ & - \frac{1}{2m} \tilde{F}^{ \nu \sigma} \partial^p_ \nu \left( p\cdot \Delta \theta_ \sigma -F_{ \sigma \lambda} \theta^ \lambda \right) ) ] \,, \end{split} \tag{ 55 }$$

$$\begin{split} 0 = & \delta(\xi) \biggl[ f_Ap\cdot \Delta \theta^ \mu - f_AF^{ \mu \nu} \theta_ \nu + \theta^ \mu p\cdot \Delta f_A \\ & -\frac{\hbar}{4m }\epsilon^{ \mu \nu \rho \alpha} p_ \alpha \left( \nabla_ \sigma F_{ \nu \rho} -p_ \lambda{R^ \lambda}_{ \sigma \nu \rho} \right) \partial^ \sigma_p f \\ & -\frac{\hbar}{2m} \tilde{F}^{ \mu \nu} \partial^p_ \nu \left( p\cdot \Delta f \right) ] \,, \end{split} \tag{ 56 }$$

with $f_{\pm}\equiv \dfrac{1}{2} \left( f \pm f_A \right)$. In Appendix D, we present the derivation of the above kinetic equations. With given initial conditions, Eqs. (55) and (56) determine the time evolutions of $f_\pm$ and $\theta^ \mu$ for massive fermions at the collisonless limit. The flat-spacetime counterparts of Eqs. (55) and (56) were discussed recently in Refs. [20-22].

We provide some comments about Eqs. (55) and (56):

(1) $\Sigma_S^{ \mu \nu}$ is related to $\theta^ \mu$ through its definition (49). Thus, in Eq. (55), it is sufficient to keep only the $O(1)$ order contribution in $\Sigma_S^{ \mu \nu}$, which is always accompanied by an additional $\hbar$ factor.

(2) The delta function in Eq. (55) shows that the on-shell condition is shifted by $\mp \hbar \Sigma_S^{ \alpha \beta} F_{ \alpha \beta}$. This term should be regarded as the magnetization coupling, similar to $\mp \hbar \Sigma_n^{ \alpha \beta} F_{ \alpha \beta}$ in the massless kinetic equation (44).

(3) Note that $f_+$ ( $f_-$) represents the distribution for fermions that have spin parallel (antiparallel) to $\theta^ \mu$. Indeed, the particle number of such spin-aligned fermions can be written with the Wigner function, as follows:

$$N_\pm \equiv \int_p {\rm{tr}} [ \mathscr{P}_\pm W] = \int_p 4 \pi \delta(\xi\mp \hbar \Sigma_S^{ \alpha \beta} F_{ \alpha \beta})\,m f_{\pm} \,, \tag{ 57 }$$

where $\mathscr{P}_\pm \equiv \dfrac{1}{2}(1\pm \gamma^5 \gamma^ \mu \theta_ \mu)$ is the spin projection operator in terms of $\theta^\mu$ [54]. Moreover, this observation of $f_\pm$ is consistent with Eq. (55); the two kinetic equations of $f_\pm$ degenerate to the same Vlasov equation $\delta(\xi) p^\mu(\partial_\mu + \Gamma_{\mu\nu}^\rho p_\rho\partial_p^\nu - F_{\mu\nu}\partial_p^\nu) f_\pm = 0$ in the classical limit, where spin-up/-down particles are indistinguishable.

(4) The third term in Eq. (56) is of $O(\hbar)$ order, as we check by substituting Eq. (55). Therefore, in the classical limit, Eq. (56) is reduced to $p\cdot \Delta \theta^ \mu = F^{ \mu \nu} \theta_ \nu$ with the on-shell condition $p^2 = m^2$. This is the Bargmann-Michel-Telegdi (BMT) equation, which describes the Larmor-Thomas precession of the spin [42]; in Minkowski spacetime, the BMT equation for a rest particle under a magnetic field ${{B}}$ is written as the well-known form of the usual Larmor precession: $m\dot{ { \theta}} = {{B\times { \theta}}}$.

(5) From Eq. (55), we extract the following single-particle equations of motion:

$$\frac{D x^ \mu}{D \tau} = \frac{p^ \mu}{m}\,, \tag{ 58 }$$

$$\frac{D p^ \mu}{D \tau} = F^{ \mu \lambda}\frac{p_ \lambda}{m}\pm \frac{\hbar}{2m} \Sigma_S^{ \alpha \beta} \left(\nabla^ \mu F_{ \alpha \beta}-p_ \lambda{R^{ \lambda\mu}}_{ \alpha \beta} \right) \,.\quad \tag{ 59 }$$

Here, $D/D \tau$ is the covariant derivative in terms of $\tau$, which is the proper time along the trajectory of the particle, and the on-shell condition $\xi\mp \hbar \Sigma_S^{ \alpha \beta} F_{ \alpha \beta} = 0$ is implicitly applied. Equation (59) is known as the first Mathisson–Papapetrou–Dixon (MPD) equation [43-45]. The first, second, and third terms in Eq. (59) represent the Coulomb-Lorentz force, the Zeeman force, and the spin curvature coupling, respectively.

(6) Multiplying by $\epsilon_{ \alpha \beta \eta \mu}p^ \eta$, Eq. (56) becomes $p\cdot \Delta \Sigma_S^{ \mu \nu} = 2F_ \sigma^{\;\;[ \mu} \Sigma_S^{ \nu] \sigma}+O(\hbar)$. Combining this with Eq. (58), the following equation of motion is derived:

$$\frac{D\hbar \Sigma_S^{ \mu \nu}}{D \tau} = 2\frac{1}{m}F_ \sigma^{\;\;[ \mu}\hbar \Sigma_S^{ \nu] \sigma}+2p^{[ \mu}\frac{Dx^{ \nu]}}{D \tau} \,. \tag{ 60 }$$

This is the second MPD equation, which determines the spin motion in electromagnetic and gravitational backgrounds [43-45]. Note that the Tulczyjew-Dixon condition [55, 56] is automatically satisfied: $p_ \mu \Sigma_S^{ \mu \nu} = 0$.

5.1. General state

As an application of our spin kinetic theory, we calculate the spin polarization of Dirac fermions, which is an intensively studied topic in heavy-ion collisions. As already mentioned, an unambiguous definition of the spin polarization is the PL vector ${\cal W}^\mu = {\cal A}^\mu$ in Eq. (22) and $\mathscr{W}^\mu(x) = \int_p {\cal W}^\mu = \int_p {\cal A}^\mu$ in Eq. (23). Combined with Eqs. (42) and (53), this polarization vector is expressed by

$$\mathscr{W}^{ \mu} (x) = \begin{cases} \displaystyle \int_p 4\pi \delta(p^2) \biggl[ p^ \mu f_5 + \hbar \Sigma_n^{ \mu \nu} \Delta_ \nu f - \dfrac{\hbar}{2} \tilde{F}^{ \mu \nu} \partial^p_ \nu f ]\quad\quad {({\rm{massless}}) ,} \\ \displaystyle \int_p 4\pi \delta(\xi) \biggl[ m \theta^ \mu f_A - \dfrac{\hbar}{2} \tilde{F}^{ \mu \nu} \partial^p_ \nu f ] \quad\quad {({\rm{massive}}) .} \end{cases} \tag{ 61 }$$

For later use, we also define the polarization per particle in the phase space:

$$w^\mu (x,p) = \frac{ {\cal W}^\mu(x,p)}{4\pi f(x,p)} = \frac{ {\cal A}^\mu(x,p)}{4\pi f(x,p)} \,. \tag{ 62 }$$

These expressions are available in the nonequilibrium state. The last terms in each case stem from the Zeeman coupling, which gives an additional $O(\hbar)$ contribution. In addition to the magnetic field, other sources, such as the fluid vorticity (or rotation), also induce spin polarization. In Eq. (61), such contributions are found, only after the concrete forms of the distribution functions are determined. For this analysis, the collision terms are needed, which we will discuss in a subsequent paper. In the global equilibrium state, however, we can identify the vorticity-dependence of the distribution functions without knowing the collision terms, as shown below.

5.2. Equilibrium state

In this section, we study spin polarization in the equilibrium state. In kinetic theory, the local equilibrium state is specified by the distribution functions that eliminate the collision kernel. This implies that the distribution functions must depend only on the linear combination of the collisional conserved quantities: the particle number, the energy and momentum, and the angular momentum. Therefore, we consider the following ansatz for the local equilibrium distributions, $f^{\rm{LE}}_{\pm} = n_F(g_{\pm})$ with $g_{\pm} = p \cdot \beta + \alpha_\pm \pm \hbar \Sigma_S^{ \mu \nu} \omega_{ \mu \nu}$ for massive fermions (where we have absorbed the orbital angular momentum into a redefinition of the $\beta$ field), and $f^{\rm{LE}}_{R/L} = n_F(g_{R/L})$ with $g_{R/ L} = p \cdot \beta + \alpha_{R/ L} \pm \hbar \Sigma_n^{ \mu \nu} \omega_{ \mu \nu}$ for massless fermions. The coefficients $\beta_ \mu, \alpha$'s, $\omega_{ \mu \nu}$ (called the spin chemical potential) depend only on x, where $\beta^ \mu$ is assumed to be time-like. Although the actual functional form of $n_F$ is not essential, we assume it to be the Fermi-Dirac function for demonstration.

5.2.1. Massive case

In the massive case, at local equilibrium, the spin polarization vectors are readily computed from Eqs. (61) and (62), as follows:

$$w^\mu_ {\rm{LE}} (x,p) = - \delta(\xi) m\theta^\mu(\alpha_A + \hbar\theta\cdot\Omega)\bar{n}_F + \hbar\delta'(\xi)\tilde{F}^{\mu\nu}p_\nu \,, \\ \tag{ 63 }$$

$$\mathscr{W}_ {\rm{LE}}^{ \mu}(x) = 4 \pi\int_p \delta(\xi) \Bigl[ 2m \theta^ \mu( \alpha_A + \hbar\Omega\cdot\theta) -\hbar \tilde{F}^{ \mu \nu} \beta_ \nu\Bigr] n_{F}' \,, \\ \tag{ 64 }$$

with $n_F = n_F(p\cdot \beta+ \alpha)$, $\bar{n}_F = 1- n_F$, $\alpha = ( \alpha_++ \alpha_-)/2$, $\alpha_A = ( \alpha_+- \alpha_-)/2$, and $\Omega^ \mu = \epsilon^{ \mu \nu \rho \sigma}p_ \nu \omega_{ \rho \sigma}/(2m)$. Note that $\alpha_A$ is assumed to be of $O(\hbar)$; otherwise, a finite spin polarization would be generated, even in the classical limit.

It is more important to discuss the polarization at global equilibrium. For this purpose, we determine the necessary constraints imposed by the kinetic equations (55) and (56). Substituting $f^{\rm{LE}}_{\pm}$ into Eq. (55), one can show that the following conditions can fulfil Eq. (55) up to $O(\hbar)$ for an arbitrary spin vector $\theta^ \mu$:

$$\nabla_ \mu \beta_ \nu + \nabla_ \nu \beta_ \mu = 0, \tag{ 65 }$$

$$\nabla_{[ \mu} \beta_{ \nu]} - 2 \omega_{ \mu \nu} = 0, \tag{ 66 }$$

$$\nabla_ \mu \alpha_{\pm} = F_{ \mu \nu} \beta^ \nu \,. \tag{ 67 }$$

Furthermore, we verify that under the conditions (65)-(67), the following choice of $\alpha_A$ and $\theta^ \mu$ fulfills Eq. (56) (see Appendix E):

$$\alpha_A = 0\,, \quad \theta^ \mu = -\frac{1}{2m \Gamma} \epsilon^{ \mu \nu \rho \sigma}p_ \nu\nabla_{ \rho} \beta_{ \sigma}, \tag{ 68 }$$

where $\Gamma = (\dfrac{1}{2} \nabla_{[ \mu} \beta_{ \nu]} \Lambda^{ \mu \rho} \Lambda^{ \nu \sigma} \nabla_{[ \rho} \beta_{ \sigma]})^{1/2}$ with $\Lambda^{ \mu \nu} = g^{ \mu \nu}- p^ \mu p^ \nu/m^2$. We call the state specified by the conditions (65)-(68) the global equilibrium state and denote $f^ {\rm{GE}}$ as the corresponding distribution function. At global equilibrium, the thermal vorticity $\nabla_{[ \mu} \beta_{ \nu]}$ determines both the spin chemical potential $\omega_{ \mu \nu}$ and the spin vector $\theta^ \mu$. We emphasize that finite Riemann curvature or an external electromagnetic field is necessary to derive Eq. (66). Without the external electromagnetic field and the curved background geometry, the spin degree of freedom is inactive in the collisionless kinetic theory and we cannot link $\omega_{ \mu \nu}$ to $\nabla_{[ \mu} \beta_{ \nu]}$. Additionally, in Appendix F, we re-derive the conditions (65)-(67) for massive fermions [and (72)-(74) for massless fermions] based on the density operator.

At global equilibrium, the spin polarization vectors read

$$w^\mu_ {\rm{GE}}(x,p) = \hbar\frac{\delta(\xi) }{2} \tilde\omega^{\mu\nu} p_\nu \bar{n}_F + \hbar\delta'(\xi) \tilde{F}^{\mu\nu} p_\nu \,, \tag{ 69 }$$

$$\mathscr{W}_ {\rm{GE}}^{ \mu}(x) = 4 \pi\hbar \int_p \delta(\xi) \left[ - \tilde{\omega}^{\mu\nu}p_\nu - \tilde{F}^{ \mu \nu} \beta_ \nu \right] n_{F}' \tag{ 70 }$$

with $\tilde{\omega}^{\mu\nu} = \dfrac{1}{2}\epsilon^{\mu\nu\rho\sigma}\omega_{\rho\sigma}$. In addition to ${\cal W}^\mu$ and $w^\mu$, at global equilibrium, it is also practically useful to compute the space-integrated polarization. Suppose that the fermions are frozen out on a space-like hypersurface $\Sigma^ \mu(x)$. The average spin polarization per particle after freeze-out is given by ^⑥

$$\begin{split} \bar{ {\cal W}}_ {\rm{GE}}^{ \mu}(p) & \equiv \frac{\displaystyle\int {\rm d} \Sigma^ \lambda p_ \lambda \displaystyle\int_0^\infty {\rm d}(p\cdot u) {\cal W}^ \mu_ {\rm{GE}}(x,p)}{4 \pi \displaystyle\int {\rm d} \Sigma^ \lambda p_ \lambda f_ {\rm{GE}}(x,p)}\\ & = \frac{\displaystyle\int {\rm d} \Sigma^ \lambda p_ \lambda \dfrac{\hbar}{4E_p} \left[ -\tilde{\omega}^{\mu\nu}p_\nu - \tilde{F}^{ \mu \nu} \beta_ \nu \right] n_{F}'}{\displaystyle\int {\rm d} \Sigma^ \lambda p_ \lambda n_F} \,. \end{split} \tag{ 71 }$$

If we set $F_{\mu\nu} = 0$, the above equation is consistent with the result derived in Refs. [46, 47]^⑥, which has been widely used for the calculation of the hadron spin polarization. In the above, $u^ \mu = T \beta^ \mu$ is the fluid velocity^⑥ and the momentum in the second line is on-shell; in Minkowski spacetime and in the local rest frame of the fluid, $p^ \mu = (E_p = \sqrt{{{p}}^2+m^2},{{p}})$ where ${{p}}$ is the three momentum.

5.2.2. Massless case

In the same manner, Eq. (44) with $f^{\rm{LE}}_{R/L}$ yields the following global equilibrium conditions [13]:

$$\nabla_ \mu \beta_ \nu + \nabla_ \nu \beta_ \mu = \phi (x) g_{ \mu \nu}, \tag{ 72 }$$

$$\nabla_{[ \mu} \beta_{ \nu]} - 2 \omega_{ \mu \nu} = 0, \tag{ 73 }$$

$$\nabla_ \mu \alpha_{R/L} = F_{ \mu \nu} \beta^ \nu \,. \tag{ 74 }$$

Unlike the massive case, the first condition has an arbitrary function $\phi (x)$, which arises as a result of the conformal invariance in the massless case; thus, $\beta^ \mu$ is a conformal Killing vector. At global equilibrium, analogously to Eqs. (69)-(71), we calculate

$$w^\mu_ {\rm{GE}}(x,p)\!\! =\!\! \frac{\hbar\delta(p^2)}{2}\bigl(- \!2p^ \mu \alpha_5/\hbar\! +\! \tilde{\omega}^{\mu\nu}p_\nu \bigr)\bar{n}_F \! +\! \hbar\delta'(p^2) \tilde{F}^{\mu\nu} p_\nu \,, \tag{ 75 }$$

$$\mathscr{W}_ {\rm{GE}}^ \mu (x) = 4 \pi\hbar \int_p \delta(p^2) \Big[ 2p^ \mu \alpha_5/\hbar - \tilde{\omega}^{\mu\nu}p_\nu - \tilde{F}^{ \mu \nu} \beta_ \nu \Big]n_F'\,, \tag{ 76 }$$

$$\bar{ {\cal W}}_ {\rm{GE}}^{ \mu}(p) = \frac{\displaystyle\int {\rm d} \Sigma^ \lambda p_ \lambda \dfrac{\hbar}{4E_p} \Big[ 2p^ \mu \alpha_5/\hbar - \tilde{\omega}^{\mu\nu}p_\nu - \tilde{F}^{ \mu \nu} \beta_ \nu \Big] n_{F}'}{\displaystyle\int {\rm d} \Sigma^ \lambda p_ \lambda n_F} \tag{ 77 }$$

with $\alpha_5 = (\alpha_R - \alpha_L)/2$, which is of $O(\hbar)$ as well as $\alpha_A$^⑥. In the second equation, the on-shell condition is implicitly applied and we define $E_p = u\cdot p$; in Minkowski spacetime and the rest frame of the fluid, $E_p = |{{p}}|$. Note that spin polarization induced by the thermal vorticity and the electromagnetic field takes the same form for both massless and massive cases at global equilibrium, up to the difference in the on-shell conditions. Moreover, the results are independent of the choice of the frame vector $n^ \mu$, as they should be.

In this paper, we derive the collisionless covariant spin kinetic theory at $O(\hbar)$ for Dirac fermions in curved spacetime and an external electromagnetic field. We start by deriving the dynamic equation for each Clifford component of the Winger function up to $O(\hbar^2)$. We discuss the physical meaning of each such dynamic equation. We then take ${\cal V}^ \mu$ and ${\cal A}^ \mu$ as independent dynamic variables and derive two evolution equations for massless fermions (44) and four evolution equations for massive fermions (55) and (56), respectively. We introduce a time-like unit frame-choosing vector $n^ \mu$ to solve the Wigner function. In the massless case, $n^ \mu$ is necessary because it represents the frame in which the spin for the massless particle is defined. In the massive case, we show that the vector $n^ \mu$ can be removed by redefining the vector distribution function through a boost from the frame $n^ \mu$ to the rest frame of the particle.

As an application, we analyze spin polarization using the approach of the kinetic theory. We derive the global equilibrium conditions from the kinetic equations and find that the finite Riemann curvature or an external electromagnetic field is necessary to determine the spin-thermal vorticity coupling. We derive the expression of spin polarization induced by the electromagnetic field and the thermal vorticity at global equilibrium, which is consistent with the results in previous literature. We also derive expressions for spin polarization at local equilibrium and out of equilibrium. They may be used to study the local $\Lambda$ polarization puzzle found in heavy-ion collisions, which cannot be understood in the calculations based on global equilibrium assumption.

We expect the spin kinetic theory to be useful for the study of both the electromagnetic plasma and quark-gluon plasma in heavy-ion collisions. Furthermore, formulating the kinetic theory in curved spacetime may find fundamental applications in astrophysics and condensed matter physics. For example, our present theory may be used to study the deformed crystal or a material with a temperature gradient, which is described as an electron system in a fictitious gravity [57, 58]. Potentially, we could study the mass correction to the chiral magnetic effect and the generation and transport of spin currents in such systems. Numerical works to solve the kinetic theory and to simulate the evolution of spin polarization in heavy-ion collisions are also important tasks. Once the collision term is included, it would be interesting to derive the covariant spin hydrodynamics [59-61] from the covariant spin kinetic theory.

We are grateful to Francesco Becattini, Gaoqing Cao, Ren-Hong Fang, Lan-Lan Gao, Xingyu Guo, Koichi Hattori, Yoshimasa Hidaka, An-Ping Huang, Jin-Feng Liao, Xin-Li Sheng, Qun Wang, Xiao-Liang Xia, Di-Lun Yang, and Pengfei Zhuang for useful discussions. We also thank the Yukawa Institute for Theoretical Physics at Kyoto University, where this work was developed during the course 'Quantum kinetic theories in magnetic and vortical fields'.

In this Appendix, we derive the chiral electromagnetic anomaly from the solutions of the Wigner function. We consider the massless case for demonstration. Plugging ${\cal A}^\mu$ in Eq. (42) into the kinetic equation (5), and integrating it, we get $\nabla_ \mu J^ \mu_{5} = \mathscr{A}$ with

$$\begin{split} \mathscr{A} = & \int_p 4\pi F^{ \mu\lambda} \partial^p_\lambda \biggl[ \hbar \delta'(p^2) \tilde{F}_{ \mu\nu}p^\nu f \\ & + \delta(p^2) \biggl( p_ \mu f_5 + \hbar\frac{ \epsilon_{ \mu \nu \rho \sigma}}{2p\cdot n}n^ \nu p^ \sigma \Delta^\rho f ) ] \\ = & \frac{\hbar}{4} F^{ \mu \nu}\tilde{F}_{ \mu \nu} \int_p 4\pi \, \delta'(p^2) p^\lambda \partial^p_\lambda f \,, \end{split} \tag{ A1 }$$

where we employ the Schouten identity and $x\delta''(x) = - 2\delta'(x)$ and drop the surface terms without the singularity at $p^2 = 0$. We have chosen the local Lorentz coordinate to perform the computation, as $\mathscr{A}$ is a scalar. The roots of $p^2 = 0$ are $p_0 = \pm | { p}|$, with which the delta function is reduced to

$$\delta(p^2) = \frac{1}{2| { p}|} \Bigl[ \delta(p_0 - | { p}|) + \delta(p_0 + | { p}|) \Bigr] \,. \tag{ A2 }$$

Furthermore, when we carry out the $p_0$-integration, we need the replacement of the momentum derivatives, as follows:

$$\begin{split} \partial_p^i f(\pm| { p}|,p_i) & = \biggl( \partial_p^i + \frac{\partial p_0}{\partial p_i}\partial_p^0 ) f(p_0,p_i) \Big |_{p_0 = \pm| { p}|} \\ & = \biggl( \partial_p^i - \frac{p^i}{p^0}\partial_p^0 ) f(p_0,p_i) \Big |_{p_0 = \pm| { p}|} \\ & \equiv \tilde{\partial_p^i} f(p_0,p_i) \Big |_{p_0 = \pm| { p}|} \,. \end{split} \tag{ A3 }$$

Subsequently, the integral in Eq. (A1) is cast into

$$\begin{split} \int_p 4\pi\delta'(p^2) p^\lambda \partial^p_\lambda f = & \int_p 4\pi \frac{1}{2}\Bigl[\partial^\lambda_p \delta(p^2)\Bigr] \partial^p_\lambda f = - \frac{1}{2} \int_p 4\pi \delta(p^2) \partial^\lambda_p \partial^p_\lambda f \\ = & - \frac{1}{2} \int_p 4\pi\delta(p^2) \biggl[ \tilde{\partial}_p^i \tilde{\partial}^p_i + \frac{2}{p_0} \partial^p_0 +\frac{2}{p_0} p^i\tilde{\partial^p_i} \partial^p_0 ] f \,. \end{split} \tag{ A4 }$$

In the last line, the second and third terms in the integrand cancel out: performing the integration by parts, we can rewrite the third term as

$$\begin{split} \int_p \delta(p^2) \frac{p^i}{p_0} \tilde{\partial}_i^p \partial^p_0 f & = \int_p \delta(p^2) \biggl[ \frac{p^i}{p_0} \partial^p_i - \frac{p^i p_i}{p_0 p^0} \partial_p^0 ] \partial^p_0 f \\ & = - \int_p \delta(p^2) \frac{1}{p_0} \partial^p_0 f \,. \end{split} \tag{ A5 }$$

Finally, $\mathscr{A}$ in Eq. (A1) is calculated as

$$\begin{split} \mathscr{A} & = - \frac{\hbar}{8} F^{ \mu \nu}\tilde{F}_{ \mu \nu} \int_ { p} \frac{1}{| { p}|} \partial_p^i \partial^p_i f \\ & = -\frac{\hbar}{16\pi^2} F^{\mu\nu} \tilde{F}_{\mu\nu} f( { p} = {\bf 0}) \,, \end{split} \tag{ A6 }$$

where we utilize

$$\partial_p^i \frac{p_i}{| { p}|^3} = 4\pi\delta^3 ( { p}) \,. \tag{ A7 }$$

The usual chiral anomaly relation is recovered; hence, we take $f( { p} = {\bf 0}) = 2$, where the factor $2$ accounts for the spin degeneracy of Dirac fermions.

We parametrize the perturbative solutions as

$${\cal V}^\mu = {\cal V}^\mu_ {(0)} +\hbar {\cal V}^\mu_ {(1)} + O(\hbar^2)\,,\quad {\cal A}^\mu = {\cal A}^\mu_ {(0)} +\hbar {\cal A}^\mu_ {(1)} + O(\hbar^2) \,. \tag{ B1 }$$

According to Eqs. (7), (8), and (10), the general solutions in the classical limit are given by

$${\cal V}_{ {(0)}}^ \mu = 4 \pi p^ \mu f^{ {(0)}} \delta(\xi) \,, \tag{ B2 }$$

$${\cal A}_{ {(0)}}^ \mu = 4 \pi \bar{ {\cal A}}_{ {(0)}}^ \mu \delta(\xi) \tag{ B3 }$$

with $\xi = p^2 - m^2$. Here, $f^{ {(0)}} = f^{ {(0)}}(x,p)$ is the classical vector charge distribution function and the vector $\bar{ {\cal A}}_{ {(0)}}^{ \mu}$ satisfies the condition $p_ \mu \bar{ {\cal A}}_{ {(0)}}^{ \mu} \delta(\xi) = 0$. Substituting Eqs. (B2) and (B3) into Eqs. (6)-(13), we obtain the solutions at $O(\hbar)$:

$$\begin{split} {\cal V}_ {(1)}^ \mu = & 4 \pi \bigg[ \left( p^ \mu f^ {(1)} +\frac{1}{2p\cdot n}\epsilon^{ \mu \nu \rho \sigma} n_ \nu \Delta_ \rho \bar{ {\cal A}}^{ {(0)}}_ \sigma \right) \delta(\xi) \\ & +\tilde{F}^{ \mu \nu} \left( \bar{ {\cal A}}^{ {(0)}}_ \nu-\frac{p\cdot\bar{ {\cal A}}^{ {(0)}}}{p\cdot n} n_ \nu \right) \delta'(\xi)\bigg] \,,\\ {\cal A}_ {(1)}^ \mu = & 4 \pi \Big[ \bar{ {\cal A}}_ {(1)}^{ \mu} \delta(\xi) + \tilde{F}^{ \mu \nu}p_ \nu f^{ {(0)}} \delta'(\xi)\Big]\,,\\ \end{split} \tag{ B4 }$$

where $f^ {(1)} = f^ {(1)}(x,p)$ is the first-order quantum correction to the vector distribution function and $\bar{ {\cal A}}_ {(1)}^{ \mu}$ satisfies the same condition as $\bar{ {\cal A}}_{ {(0)}}^{ \mu}$ : $p_ \mu \bar{ {\cal A}}_ {(1)}^{ \mu} \delta(\xi) = 0$. Defining $f\equiv f^{ {(0)}}+\hbar f^ {(1)}$ and $\bar{ {\cal A}}^{ \mu}\equiv\bar{ {\cal A}}_{ {(0)}}^{ \mu}+\hbar\bar{ {\cal A}}_ {(1)}^{ \mu}$, we obtain Eqs. (36) and (37).

Here, we show that, with the redefinition of the distribution function f in Eq. (47), we can remove the frame vector $n^ \mu$ from the spin kinetic theory for massive fermions. The discussion is kept at $O(\hbar)$. Acting on Eq. (9) with $\Delta_ \alpha$, we derive

$$\begin{split} p\cdot \Delta {\cal A}^ \beta + F^{ \alpha \beta} {\cal A}_ \alpha-p^ \beta \Delta_ \alpha {\cal A}^ \alpha = \frac{m}{2} \epsilon^{ \alpha \beta \rho \sigma} \Delta_ \alpha {\cal S}_{ \rho \sigma}-\frac{\hbar}{2} \epsilon^{ \alpha \beta \rho \sigma} \Delta_ \alpha \Delta_ \rho {\cal V}_ \sigma \,. \end{split} \tag{ C1 }$$

Using Eqs. (5) and (12), we obtain

$$p^ \beta \Delta_ \alpha {\cal A}^ \alpha = \frac{m}{2} \epsilon^{ \beta \alpha \rho \sigma} \Delta_ \alpha {\cal S}_{ \rho \sigma}. \tag{ C2 }$$

Combining the above two equations, we find

$$p\cdot \Delta {\cal A}_ \mu = F_{ \mu \nu} {\cal A}^ \nu+\frac{\hbar}{2} \epsilon_{ \mu \nu \rho \sigma} \Delta^ \nu \Delta^ \rho {\cal V}^ \sigma \,. \tag{ C3 }$$

Next, we substitute the redefined distribution function in Eq. (47) into the solution of ${\cal V}^ \mu$ in Eq. (36) and obtain

$$\begin{split} {\cal V}^ \mu = & 4 \pi \delta(\xi) \bigg[ p^ \mu \left( f - \frac{\hbar \epsilon^{ \alpha \beta \rho \sigma} p_ \alpha n_ \beta}{2m^2 p\cdot n} \Delta_ \rho \bar{ {\cal A}}_{\perp \sigma} \right) \\ & +\frac{\hbar\epsilon^{ \mu \nu \rho \sigma} n_ \nu}{2p\cdot n} \Delta_ \rho \bar{ {\cal A}}^\bot_ \sigma \bigg] +4 \pi \delta'(\xi)\hbar\tilde{F}^{ \mu \nu} \bar{ {\cal A}}^\bot_ \nu \,. \end{split} \tag{ C4 }$$

To reduce the above equation, we utilize the Schouten identity:

$$\begin{split} p^ \mu \epsilon^{ \alpha \beta \rho \sigma} p_ \alpha n_ \beta \Delta_ \rho \bar{ {\cal A}}_{\perp\sigma} = & - \Bigl( p^2 \epsilon^{ \beta \rho \sigma \mu} n_ \beta \Delta_ \rho + p\cdot n \epsilon^{ \rho \sigma \mu \alpha} p_ \alpha \Delta_ \rho \\ & + \epsilon^{ \sigma \mu \alpha \beta} p_ \alpha n_ \beta p\cdot \Delta + \epsilon^{ \mu \alpha \beta \rho} p_ \alpha n_ \beta p^ \sigma \Delta_ \rho \Bigr) \bar{ {\cal A}}_{\perp\sigma} \,, \end{split} \tag{ C5 }$$

where the last two terms cancel, according to Eq. (C3). We then rewrite Eq. (C4) as the $n^ \mu$ independent form:

$${\cal V}^ \mu = 4 \pi\delta(\xi) \bigg[ p^ \mu f +\frac{\hbar\epsilon^{ \mu \nu \rho \sigma} p_ \nu}{2m^2} \Delta_ \rho \bar{ {\cal A}}^\bot_ \sigma -\frac{\hbar}{\xi}\tilde{F}^{ \mu \nu} \bar{ {\cal A}}^\bot_ \nu \bigg] \,. \tag{ C6 }$$

The frame vector $n^ \mu$ is also removed from the solution of ${\cal A}^ \mu$. Comparing the above equation with Eq. (36), we find that the redefinition of f is equivalent to replacing $n^ \mu$ with $p^ \mu/m$ in Eq. (36).

Here, we derive the kinetic equations for $f_{\pm}$. Substituting Eq. (52) into Eq. (4), one obtains the following equation:

$$\begin{split} 0 & = p^ \mu \Delta_ \mu f \delta(\xi) - \hbar \Sigma_S^{ \alpha \beta} F_{ \alpha \beta} p^ \mu \Delta_ \mu f_A \delta'(\xi)\\ & +\frac{\hbar}{2} \left( \nabla_ \rho F_{ \mu \nu} \partial^ \rho_p+ \left[ D_ \mu, D_ \nu \right] \right) \left( \Sigma_S^{ \mu \nu} f_A \right) \delta(\xi)\\ & {- \hbar\frac{1}{2m} f_A \tilde{F}^{ \nu \sigma} \partial^p_ \nu \left( p^ \rho \Delta_ \rho \theta_ \sigma -F_{ \sigma \lambda} \theta^ \lambda \right) \delta(\xi)}. \end{split} \tag{ D1 }$$

In addition, contracting Eq. (C3) with $\theta^ \mu$ and inserting Eqs. (52) and (53), we obtain

$$\begin{split} 0 = & p^ \mu \Delta_ \mu f_A \delta(\xi) - \hbar \Sigma_S^{ \alpha \beta} F_{ \alpha \beta} p^ \mu \Delta_ \mu f \delta'(\xi)\\ & +\frac{\hbar}{2} \Sigma_S^{ \mu \nu} \left( \nabla_ \rho F_{ \mu \nu} \partial^ \rho_p+ \left[ D_ \mu, D_ \nu \right] \right) f \delta(\xi). \end{split} \tag{ D2 }$$

The addition and subtraction of the above two equations result in Eq. (55). Moreover, the kinetic equation to determine $\theta^\mu$ is obtained from Eq. (C3), with the solutions (52) and (53).

In the massless case, the discussion of the global equilibrium conditions (65)-(67) was given in Ref. [13]. Following a similar strategy, one can show that, for the massive case, the conditions (65)-(67) can fulfill Eq. (55) for arbitrary $\theta^ \mu$ and for $\alpha_A = O(\hbar)$. Also, it is easy to see that the condition (68) fulfills Eq. (55) up to $O(\hbar)$. We verify that condition (68) also fulfills Eq. (56) under the conditions (65)-(67), as follows. The leading order of $f_A^{\rm{LE}}$ can be written as $f^ {\rm{LE}}_A = 2( \alpha_A+\hbar \Sigma^S_{ \alpha \beta} \omega^{ \alpha \beta}) n_F'( \beta\cdot p+ \alpha) + O(\hbar^2)$. Using Eqs. (65)-(68) and inserting $f^ {\rm{LE}}$ and $f_A^ {\rm{LE}}$, we obtain

$$\begin{split} & {\rm{RHS}}\;{\rm{of}}\;{\rm{Eq}}.\;({\rm{56}}) \\ & = 2 \delta(\xi) n'_F \biggl[ \hbar p\cdot \Delta \left( \theta^ \mu \Sigma^S_{ \alpha \beta} \omega^{ \alpha \beta} \right) - \left( \alpha_A+\hbar \Sigma^S_{ \alpha \beta} \omega^{ \alpha \beta} \right) F^{ \mu \nu} \theta_ \nu \\ & - \frac{\hbar}{4m}\epsilon^{ \mu \nu \rho \alpha} p_ \alpha \left( \nabla_ \sigma F_{ \nu \rho} -p_ \lambda{R^ \lambda}_{ \sigma \nu \rho} \right) \beta^\sigma ] \\ & = - \hbar\frac{ \delta(\xi)}{2m} n'_F \Bigl[ p\cdot \Delta \left( \epsilon^{ \mu \nu \rho \sigma}p_ \nu\nabla_{ \rho} \beta_{ \sigma} \right) - F^{ \mu \nu} \epsilon_{ \nu \lambda \rho \sigma}p^ \lambda\nabla^{ \rho} \beta^{ \sigma} \\ & + \epsilon^{ \mu \nu \rho \alpha} p_ \alpha \left( \nabla_ \sigma F_{ \nu \rho} -p_ \lambda{R^ \lambda}_{ \sigma \nu \rho} \right) \beta^ \sigma \Bigr] \,. \end{split} \tag{ E1 }$$

In the above equation, the second equality follows from $\alpha_A = 0$ and

$$\theta^\mu\Sigma^{\alpha\beta}_S\omega_{\alpha\beta} = \frac{1}{2} \theta^\mu \Gamma = -\frac{1}{4m} \epsilon^{\mu\nu\rho\sigma}p_\nu\nabla_\rho\beta_\sigma \,, \tag{ E2 }$$

which stems from Eq. (68). The above three terms in Eq. (E1) totally vanish, as follows:

$$\begin{split} p\cdot \Delta \left( \epsilon^{ \mu \nu \rho \sigma}p_ \nu\nabla_{ \rho} \beta_{ \sigma} \right) = & p_ \lambda \left( \epsilon^{ \nu \rho \sigma \lambda}F^{ \mu}_{\;\; \nu} + 2 \epsilon^{ \sigma \lambda \mu \nu}F^{ \rho}_{\;\; \nu} \right) \nabla_{ \rho} \beta_{ \sigma} \\ & + \epsilon^{ \mu \nu \rho \sigma}p_ \nu p_ \lambda R^ \lambda_{\;\; \alpha \rho \sigma} \beta^ \alpha \\ = & F^{ \mu}_{\;\; \nu} \epsilon^{ \nu \rho \sigma \lambda}p_ \lambda \nabla_{ \rho} \beta_{ \sigma} -p_ \lambda \epsilon^{ \mu \nu \sigma \lambda} \beta\cdot \nabla F_{ \nu \sigma} \\ & + \epsilon^{ \mu \nu \rho \sigma}p_ \nu p_ \lambda R^ \lambda_{\;\; \alpha \rho \sigma} \beta^ \alpha \,. \end{split} \tag{ E3 }$$

In the above equation, we have used the Schouten identity and the equilibrium conditions (65) and (67) with $\nabla_{\mu} \nabla_{[\nu}\beta_{\rho]} = - \beta^\lambda R_{\lambda\mu\nu\rho}$.

We discuss the global equilibrium using the maximum entropy principle, following Refs. [62-64]. The density operator for the local equilibrium state can be written as

$$\hat{ \rho}_ {\rm{LE}} \equiv \frac{1}{z} {\rm e}^{-\int {\rm d} \Sigma_ \mu\left(\hat{T}^{ \mu \nu} \beta_ \nu + \hat{S}^{ \mu, \lambda \nu} \omega_{ \lambda \nu} + \alpha \hat{J}^ \mu\right)}\,, \tag{ F1 }$$

where $z \equiv {{\rm{Tr}}} \left[ {\rm e}^{-\int {\rm d} \Sigma_ \mu\left(\hat{T}^{ \mu \nu} \beta_ \nu+ \hat{S}^{ \mu, \lambda \nu} \omega_{ \lambda \nu} + \alpha \hat{J}^ \mu\right)} \right]$ where $\hat{T}^{ \mu \nu}$, $\hat{S}^{ \mu, \lambda \nu}$, and $\hat{J}^ \mu$ are the canonical energy-momentum, spin, and charge current operators, respectively. Here, $\Sigma_ \mu$ is a space-like hypersurface, $\beta^ \mu$, $\alpha$, and $\omega_{ \mu \nu}$ have the same meanings as in the main text. The entropy is defined as

$$S \equiv -\langle\ln \hat{ \rho}_ {\rm{LE}}\rangle = - {{\rm{Tr}}} \left( \hat{ \rho}_ {\rm{LE}} \ln \hat{ \rho}_ {\rm{LE}} \right) \,. \tag{ F2 }$$

We denote $\ln z = \int {\rm d} \Sigma_ \mu \phi^ \mu$, where $\phi^ \mu$ is (the negative of) the thermodynamic potential density current. Then, the entropy is represented as $S = \int {\rm d} \Sigma_ \mu s^ \mu$, with

$$s^ \mu = \phi^ \mu + T^{ \mu \nu} \beta_ \nu + S^{ \mu, \lambda \nu} \omega_{ \lambda \nu} + \alpha J^ \mu \,. \tag{ F3 }$$

The global equilibrium condition is such that the local thermodynamic potential and entropy are maximized, so that $\nabla_ \mu \phi^ \mu = 0$ and $\nabla_ \mu s^ \mu = 0$. After some straightforward calculations, we arrive at

$$\begin{split} 0 = & T^{ \mu \nu}_{\rm sy} \nabla_{ \mu} \beta_{ \nu} + T^{ \mu \nu}_{\rm as} \left( \nabla_{ \mu} \beta_{ \nu} - 2 \omega_{ \mu \nu} \right) \\ & + S^{ \mu, \lambda \nu} \nabla_ \mu \omega_{ \lambda \nu} + J^ \mu ( \nabla_ \mu \alpha-F_{ \mu \nu} \beta^ \nu) \,, \end{split} \tag{ F4 }$$

where $T^{ \mu \nu}_{\rm sy/as}$ is the symmetric/antisymmetric part of $T^{ \mu \nu}$. In the massless case, ${T^ \mu}_ \mu = 0$, we obtain Eqs. (72)-(74) (with $\alpha_R = \alpha_L = \alpha$); in the massive case, we obtain Eqs. (65)-(67). Note that one further constraint from Eq. (F4), $\nabla_{[ \mu} \omega_{ \lambda \nu]} = 0$, is automatically fulfilled.