Covariant Spin Kinetic Theory I: Collisionless Limit

  • We develop a covariant kinetic theory for massive fermions in curved spacetime and external electromagnetic field based on quantum field theory. We derive four coupled semi-classical kinetic equations accurate at $O(\hbar)$, which describe the transports of particle number and spin degrees of freedom. The relation with the chiral kinetic theory is discussed. As an application, we study the spin polarization in the presence of finite Riemann curvature and electromagnetic field in both local and global equilibrium states.
  • 加载中
  • [1] D. T. Son and N. Yamamoto, Phys. Rev. Lett., 109: 181602 (2012), arXiv:1203.2697[cond-mat.mes-hall doi: 10.1103/PhysRevLett.109.181602
    [2] M. A. Stephanov and Y. Yin, Phys. Rev. Lett., 109: 162001 (2012), arXiv:1207.0747[hep-th doi: 10.1103/PhysRevLett.109.162001
    [3] J.-H. Gao, Z.-T. Liang, S. Pu, Q. Wang, and X.-N. Wang, Phys. Rev. Lett., 109: 232301 (2012), arXiv:1203.0725[hep-ph doi: 10.1103/PhysRevLett.109.232301
    [4] D. T. Son and N. Yamamoto, Phys. Rev., D87: 085016 (2013), arXiv:1210.8158[hep-th
    [5] J.-W. Chen, S. Pu, Q. Wang, and X.-N. Wang, Phys. Rev. Lett., 110: 262301 (2013), arXiv:1210.8312[hep-th doi: 10.1103/PhysRevLett.110.262301
    [6] J.-Y. Chen, D. T. Son, M. A. Stephanov, H.-U. Yee, and Y. Yin, Phys. Rev. Lett., 113: 182302 (2014), arXiv:1404.5963[hep-th doi: 10.1103/PhysRevLett.113.182302
    [7] J.-Y. Chen, D. T. Son, and M. A. Stephanov, Phys. Rev. Lett., 115: 021601 (2015), arXiv:1502.06966[hep-th doi: 10.1103/PhysRevLett.115.021601
    [8] Y. Hidaka, S. Pu, and D.-L. Yang, Phys. Rev., D95: 091901 (2017), arXiv:1612.04630[hep-th
    [9] N. Mueller and R. Venugopalan, Phys. Rev., D96: 016023 (2017), arXiv:1702.01233[hep-ph
    [10] S. Carignano, C. Manuel, and J. M. Torres-Rincon, Phys. Rev., D98: 076005 (2018), arXiv:1806.01684[hep-ph
    [11] A. Huang, S. Shi, Y. Jiang, J. Liao, and P. Zhuang, Phys. Rev., D98: 036010 (2018), arXiv:1801.03640[hep-th
    [12] X.-G. Huang and A. V. Sadofyev, JHEP, 03: 084 (2019), arXiv:1805.08779[hep-th
    [13] Y.-C. Liu, L.-L. Gao, K. Mameda, and X.-G. Huang, Phys. Rev., D99: 085014 (2019), arXiv:1812.10127[hep-th
    [14] S. Lin and A. Shukla, JHEP, 06: 060 (2019), arXiv:1901.01528[hep-ph
    [15] D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, Nucl. Phys., A803: 227 (2008), arXiv:0711.0950[hep-ph
    [16] K. Fukushima, D. E. Kharzeev, and H. J. Warringa, Phys. Rev., D78: 074033 (2008), arXiv:0808.3382[hep-ph
    [17] A. Vilenkin, Phys. Rev., D20: 1807 (1979)
    [18] J. Erdmenger, M. Haack, M. Kaminski, and A. Yarom, JHEP, 01: 055 (2009), arXiv:0809.2488[hep-th
    [19] N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Dutta, R. Loganayagam, and P. Surowka, JHEP, 01: 094 (2011), arXiv:0809.2596[hep-th
    [20] J.-H. Gao and Z.-T. Liang, Phys. Rev., D100: 056021 (2019), arXiv:1902.06510[hep-ph
    [21] N. Weickgenannt, X.-L. Sheng, E. Speranza, Q. Wang, and D. H. Rischke, Phys. Rev., D100: 056018 (2019), arXiv:1902.06513[hep-ph
    [22] K. Hattori, Y. Hidaka, and D.-L. Yang, Phys. Rev., D100: 096011 (2019), arXiv:1903.01653[hep-ph
    [23] Z. Wang, X. Guo, S. Shi, and P. Zhuang, Phys. Rev., D100: 014015 (2019), arXiv:1903.03461[hep-ph
    [24] Z.-T. Liang and X.-N. Wang, Phys. Rev. Lett. 94, 102301 (2005), [Erratum: Phys. Rev. Lett.96, 039901(2006)], arXiv: nucl-th/0410079[nucl-th].
    [25] J.-H. Gao, S.-W. Chen, W.-T. Deng, Z.-T. Liang, Q. Wang, and X.-N. Wang, Phys. Rev., C77: 044902 (2008), arXiv:0710.2943[nucl-th
    [26] F. Becattini, F. Piccinini, and J. Rizzo, Phys. Rev., C77: 024906 (2008), arXiv:0711.1253[nucl-th
    [27] X.-G. Huang, P. Huovinen, and X.-N. Wang, Phys. Rev., C84: 054910 (2011), arXiv:1108.5649[nucl-th
    [28] L. Adamczyk EM et al. (STAR), Nature, 548: 62 (2017), arXiv:1701.06657[nucl-ex doi: 10.1038/nature23004
    [29] W.-T. Deng and X.-G. Huang, Phys. Rev., C93: 064907 (2016), arXiv:1603.06117[nucl-th
    [30] Y. Jiang, Z.-W. Lin, and J. Liao, Phys. Rev. C94, 044910 (2016), [Erratum: Phys. Rev.C95, no.4, 049904(2017)], arXiv: 1602.06580[hep-ph].
    [31] X.-G. Deng, X.-G. Huang, Y.-G. Ma, and S. Zhang, (2020), arXiv: 2001.01371[nucl-th].
    [32] J. Adam EM et al. (STAR), Phys. Rev., C98: 014910 (2018), arXiv:1805.04400[nucl-ex
    [33] J. Adam EM et al. (STAR), Phys. Rev. Lett., 123: 132301 (2019), arXiv:1905.11917[nucl-ex doi: 10.1103/PhysRevLett.123.132301
    [34] I. Karpenko and F. Becattini, Eur. Phys. J., C77: 213 (2017), arXiv:1610.04717[nucl-th
    [35] F. Becattini and I. Karpenko, Phys. Rev. Lett., 120: 012302 (2018), arXiv:1707.07984[nucl-th doi: 10.1103/PhysRevLett.120.012302
    [36] X.-L. Xia, H. Li, Z.-B. Tang, and Q. Wang, Phys. Rev., C98: 024905 (2018), arXiv:1803.00867[nucl-th
    [37] Y. Sun and C. M. Ko, Phys. Rev., C99: 011903 (2019), arXiv:1810.10359[nucl-th
    [38] D.-X. Wei, W.-T. Deng, and X.-G. Huang, Phys. Rev., C99: 014905 (2019), arXiv:1810.00151[nucl-th
    [39] X.-L. Xia, H. Li, X.-G. Huang, and H. Z. Huang, Phys. Rev., C100: 014913 (2019), arXiv:1905.03120[nucl-th
    [40] F. Becattini, G. Cao, and E. Speranza, Eur. Phys. J., C79: 741 (2019), arXiv:1905.03123[nucl-th
    [41] H.-Z. Wu, L.-G. Pang, X.-G. Huang, and Q. Wang, Phys. Rev. Research., 1: 033058 (2019), arXiv:1906.09385[nucl-th doi: 10.1103/PhysRevResearch.1.033058
    [42] V. Bargmann, L. Michel, and V. L. Telegdi, Phys. Rev. Lett., 2: 435 (1959) doi: 10.1103/PhysRevLett.2.435
    [43] M. Mathisson, Acta Phys. Polon., 6: 163 (1937)
    [44] W. G. Dixon, Proc. Roy. Soc. Lond., A314: 499 (1970)
    [45] A. Papapetrou, Proc. Roy. Soc. Lond., A209: 248 (1951)
    [46] F. Becattini, V. Chandra, L. Del Zanna, and E. Grossi, Annals Phys., 338: 32 (2013), arXiv:1303.3431[nucl-th doi: 10.1016/j.aop.2013.07.004
    [47] R.-H. Fang, L.-G. Pang, Q. Wang, and X.-N. Wang, Phys. Rev., C94: 024904 (2016), arXiv:1604.04036[nucl-th
    [48] D. Vasak, M. Gyulassy, and H. T. Elze, Annals Phys., 173: 462 (1987) doi: 10.1016/0003-4916(87)90169-2
    [49] F. Becattini, W. Florkowski, and E. Speranza, Phys. Lett., B789: 419 (2019), arXiv:1807.10994[hep-th
    [50] B. Maybee, D. O’Connell, and J. Vines, JHEP, 12: 156 (2019), arXiv:1906.09260[hep-th
    [51] D. T. Son and A. R. Zhitnitsky, Phys. Rev., D70: 074018 (2004), arXiv:hep-ph/0405216[hep-ph
    [52] M. A. Metlitski and A. R. Zhitnitsky, Phys. Rev., D72: 045011 (2005), arXiv:hep-ph/0505072[hep-ph
    [53] L. N. Chang and H. T. Nieh, Phys. Rev. Lett., 53: 21 (1984) doi: 10.1103/PhysRevLett.53.21
    [54] M. Kaku, EMQuantum field theory: A Modern introduction (Oxford Univ. Pr., New York, 1993).
    [55] W. Tulczyjew, Acta Phys. Pol., 18: 393 (1959)
    [56] W. G. Dixon, Phil. Trans. R. Soc. Lond. Ser. A, 277: 59 (1974) doi: 10.1098/rsta.1974.0046
    [57] J. M. Luttinger, Phys. Rev., 135: A1505 (1964) doi: 10.1103/PhysRev.135.A1505
    [58] L. Dong and Q. Niu, Phys. Rev., B98: 115162 (2018), arXiv:1802.02887[cond-mat.mtrl-sci
    [59] W. Florkowski, B. Friman, A. Jaiswal, and E. Speranza, Phys. Rev., C97: 041901 (2018), arXiv:1705.00587[nucl-th
    [60] W. Florkowski, R. Ryblewski, and A. Kumar, Prog. Part. Nucl. Phys., 108: 103709 (2019), arXiv:1811.04409[nucl-th doi: 10.1016/j.ppnp.2019.07.001
    [61] K. Hattori, M. Hongo, X.-G. Huang, M. Matsuo, and H. Taya, Phys. Lett., B795: 100 (2019), arXiv:1901.06615[hep-th
    [62] D. N. Zubarev, EMNonequilibrium statistical thermodynamics (Springer, 1974).
    [63] F. Becattini, Phys. Rev. Lett., 108: 244502 (2012), arXiv:1201.5278[gr-qc doi: 10.1103/PhysRevLett.108.244502
    [64] F. Becattini and D. Rindori, Phys. Rev., D99: 125011 (2019), arXiv:1903.05422[hep-th
  • 加载中

Get Citation
Yu-Chen Liu, Kazuya Mameda and Xu-Guang Huang. Covariant Spin Kinetic Theory I: Collisionless Limit[J]. Chinese Physics C. doi: 10.1088/1674-1137/44/9/094101
Yu-Chen Liu, Kazuya Mameda and Xu-Guang Huang. Covariant Spin Kinetic Theory I: Collisionless Limit[J]. Chinese Physics C.  doi: 10.1088/1674-1137/44/9/094101 shu
Article Metric

Article Views(289)
PDF Downloads(7)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌,
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article


Covariant Spin Kinetic Theory I: Collisionless Limit

  • 1. Physics Department and Center for Particle Physics and Field Theory, Fudan University, Shanghai 200433, China
  • 2. RIKEN Nishina Center for Accelerator-Based Science, Wako 351-0198, Japan
  • 3. Key Laboratory of Nuclear Physics and Ion-beam Application (MOE), Fudan University, Shanghai 200433, China

Abstract: We develop a covariant kinetic theory for massive fermions in curved spacetime and external electromagnetic field based on quantum field theory. We derive four coupled semi-classical kinetic equations accurate at $O(\hbar)$, which describe the transports of particle number and spin degrees of freedom. The relation with the chiral kinetic theory is discussed. As an application, we study the spin polarization in the presence of finite Riemann curvature and electromagnetic field in both local and global equilibrium states.


    1.   Introduction
    • Kinetic theory is widely used to study the transport phenomena in many-particle systems. The classical Boltzmann kinetic theory is established as the framework to describe the evolution of the distribution function in phase space. To study the effects of spin, one has to use the quantum kinetic theory and single distribution function is usually insufficient for such a purpose. For massless Dirac fermions, the leading spin effects appear at $ O(\hbar) $, and we need two distribution functions, i.e., one for right-handed chirality and the other for left-handed chirality. Such established framework is the so-called chiral kinetic theory (CKT) [1-3], which has been under intensive investigations recently [4-14]. The out-of-equilibrium dynamics of anomaly-induced phenomena, e.g., the chiral magnetic effect [15, 16] and the chiral vortical effect [17-19], have been commendably 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 dynamical evolution of massive fermions is specified with more than two degrees of freedom. Suppose that $ \theta^ \mu $ is the unit spacelike vector specifying the spin quantizing orientation. Then the dynamical variables are two parameters 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 we need extensive concerns of this framework [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 very 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 $ polarization show the opposite azimuthal angle dependence compared with the so-called thermal vorticity [34, 41]. This strongly indicates that the spin polarization has independent dynamical evolution, in a non-equilibrium state, rather than being enslaved to the fluid vorticity. The covariant kinetic theory for spin transport (hereafter, spin kinetic theory for short) would be one of the promising tools to capture the dynamics of the spin polarization.

      In this paper, we derive the collisionless spin kinetic theory at $ O(\hbar) $ in curved background spacetime and external electromagnetic field. As an application, we investigate the spin polarization from the spin kinetic theory. We give the general expression for the 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 equilibrium. We stress that the present study differs from the earlier works [20-23] in the following aspects: (i) We include curved geometric background spacetime as well as electromagnetic field. Such a general formalism should be applicable to the spin transport not only in heavy-ion collisions and astrophysical systems, but also in deformed materials and thermal-gradient systems, which attracts huge attention in condensed matter physics. (ii) We show that the frame-choosing vector can always be eliminated in the covariant kinetic theory of the massive fermions, unlike the massless case. (iii) We discuss the underlying physics of the Clifford components and their constraint equations. (iv) We write down the kinetic equations in a more transparent way, exhibiting clear physical contexts. Especially, we verify that in the classical limit, these equations are correctly reduced to the Vlasov equation, the Bargmann–Michel–Telegdi (BMT) equation [42], and Mathisson–Papapetrou–Dixon (MPD) equations [43-45]. (v) 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 II and III, we introduce the Wigner function and discuss the physical meaning of the dynamical equation for each Clifford component of the Wigner function. In Section IV, we derive the semi-classical kinetic theory for massive fermions. In Section V, we derive the kinetic representation of the spin polarization for both massive and massless fermions and investigate the spin polarization at both local and global equilibrium. In this paper, we adopt the same notations and conventions as those of 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 one).

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

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


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


      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 \,. $


      This $ D_\mu $ gives us a great advantage of the analysis due to 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+ iA_ \mu/\hbar) \psi(x,y) $ with $ \psi(x,y)\equiv 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, and thus we impose the spinor field to obey the Dirac equation $ (i\hbar\gamma^\mu D_\mu-m)\psi(x) = \bar{\psi}(x)\,(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] (especially, see Section III and Appendices C, D, and E therein). After the semiclassical expansion, and the decomposition in terms of the Clifford algebra as $ W = \frac{1}{4}[{\cal F}+i \gamma^5{\cal P}+ \gamma^ \mu{\cal V}_ \mu + \gamma^5 \gamma^ \mu{\cal A}_ \mu+ \frac{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 \,,$


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


      $ \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 \,, $


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


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


      $ \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 \,, $


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


      $ \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} \,, $


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


      $ \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\ \ $


      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} = \frac{1}{2}\epsilon^{\mu\nu\rho\sigma} {\cal S}_{\rho\sigma} $ and $ \tilde{R}^{\mu\nu\alpha\beta} = \frac{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} $


      In the above 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. Then $ {\cal P} $, $ {\cal F} $, and $ {\cal S}_{ \mu \nu} $ are expressed in terms of $ {\cal V}_ \mu $ and $ {\cal A}_ \mu $ through 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 \,, $


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


      with $ \int_p \equiv\int \dfrac{d^4 p}{(2 \pi)^4 [-g(x)]^{1/2}} $. From these, one identifies the Clifford coefficients $ {\cal V}_ \mu $ and $ {\cal A}_ \mu $ as the corresponding current densities in phase space (see more discussions in Sec. V). In a similar way, one finds that $ {\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 latter 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 \,, $


      $ 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 \,, $


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


      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 $.

      It is worthwhile to 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 one introduces the Belinfante tensor. Instead, an unambiguous way to represent particles 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} \,, $


      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 notice that the orbital part of the canonical angular momentum does not contribute the above equation.

      Following Eq. (20), we may define the PL vector in our kinetic theory as

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


      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. One can readily check that the above definition of $ {\cal W}^\mu $ excludes the orbital angular momentum part, and it is reduced as

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


      Note that this relation is independent of the vector $ \nu^\mu $. The coincidence between $ {\cal W}^\mu $ and $ {\cal A}^\mu $ is expectable; e.g., 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) \,. $

    3.   Physical Interpretation
    • In the following we discuss the physical meanings of Eqs. (4)-(13). To show this, we perform the 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, hereafter we 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) leads to fundamental Ward identities. After integrating over momentum, Eq. (4) gives the vector current conservation law:

      $ \nabla_ \mu J^ \mu = 0 \,. $


      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 the external field:

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


      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 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} $. Also integrating Eq. (6) over momentum, we obtain

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


      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 $


      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), whose physical contents are less transparent than Eqs. (4) and (6). Equation (8) involves $ {\cal A}^ \mu $ only, so it is a subsidiary condition for $ {\cal A}^ \mu $. Up to $ O(\hbar) $, it reduces to

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


      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 IV we discuss the details even 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]} \,, $


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

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


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


      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} \,,$


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


      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)]. In this sense, we novelly understand that Eqs. (32) and (33) imply as the following fact: 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}, $


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

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


      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 the discussion of the higher order kinetic theory in future publication.

    4.   Kinetic equations at $ O(\hbar) $
    • 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} $


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


      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 \,. $


      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 \,. $


      Here $ \bar{ {\cal A}}^{ \mu} $ is not necessarily perpendicular to the momentum due to 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 \,, $


      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} \,. $


      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} $


      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}) $. Besides, it is important to notice 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}\,. $


      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} $


      More discussions 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} $


      with $ \int_ { p} \equiv\int \dfrac{d^3 p}{(2 \pi)^3} $. This reproduces the usual chiral anomaly when f is (the twice 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 electromagnetic field [1, 2].

    • 4.2.   Massive case

    • Let us 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 \,, $


      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{1}{m^2}\frac{\epsilon_{ \mu \nu \rho \sigma} p^ \mu n^ \nu}{2 p\cdot n} \Delta^ \rho \bar{ {\cal A}}_{\perp}^\sigma \,. $


      We emphasize that this redefinition is equivalent to identify 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 since 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 due to the lack of such a special frame, and it is inevitable to introduce $ n^ \mu $.

      Due to 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 $. Since 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 \,, $


      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 \,, $


      which may be regarded as the spin tensor of massive fermions, as so is $ \Sigma_n^{\mu\nu} $ for massless fermions (41). Indeed it is readily checked 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 ], $


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


      $ \begin{split} {\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} $


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


      $ \begin{split} {\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}$


      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} $


      $ \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} $


      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 collisonless limit. The flat-spacetime counterparts of Eqs. (55) and (56) were discussed recently in Refs. [20-22].

      We give 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 onshell condition is shifted by $ \mp \hbar \Sigma_S^{ \alpha \beta} F_{ \alpha \beta} $. This term should be regarded as the magnetization coupling, similarly 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 whose spin is parallel (antiparallel) to $ \theta^ \mu $. Indeed, the particle number of such spin-aligned fermions are 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} \,, $


      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 on-shell condition $ p^2 = m^2 $. This is nothing but the the Bargmann–Michel–Telegdi (BMT) equation, which describes the Larmor-Thomas procession 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 procession: $ 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}\,, $


      $ \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 $


      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 term in Eq. (59) represent the Coulomb-Lorentz force, the Zeeman force, and the spin curvature coupling, respectively.

      (6) By multiplying $ \epsilon_{ \alpha \beta \eta \mu}p^ \eta $, one recasts Eq. (56) into $ p\cdot \Delta \Sigma_S^{ \mu \nu} = 2F_ \sigma^{\;\;[ \mu} \Sigma_S^{ \nu] \sigma}+O(\hbar) $. Combing this with Eq. (58), one derives the following equation of motion:

      $ \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} \,. $


      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.   Application: spin polarization

      5.1.   General state

    • As an application of our spin kinetic theory, we calculate the spin polarization of Dirac fermions, which is a topic under intensive study in heavy-ion collisions. As we have 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} \int_p 4\pi \delta(p^2) \biggl[ p^ \mu f_5 + \hbar \Sigma_n^{ \mu \nu} \Delta_ \nu f - \frac{\hbar}{2} \tilde{F}^{ \mu \nu} \partial^p_ \nu f ]\quad\quad {({\rm{massless}}) ,} \\ \int_p 4\pi \delta(\xi) \biggl[ m \theta^ \mu f_A - \frac{\hbar}{2} \tilde{F}^{ \mu \nu} \partial^p_ \nu f ] \quad\quad {({\rm{massive}}) .} \end{cases} $


      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)} \,. $


      These expressions are available in nonequilibrium state. The last terms in each case stem from the Zeeman coupling, which gives an additional $ O(\hbar) $ contribution. On top of magnetic field, other sources, e.g., the fluid vorticity (or rotation), also induces the spin polarization. In Eq. (61), such contributions are found, only after the concrete forms of the distribution functions are determined. For this analysis, one needs the collision terms, which we will leave in a subsequent paper. In 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 the following, we study the 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 the massive fermions (we have absorbed the orbital angular momentum into a redefinition of $ \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 the massless fermions. The coefficients $ \beta_ \mu, \alpha $'s, $ \omega_{ \mu \nu} $ (called spin chemical potential) depend only on x, where $ \beta^ \mu $ is assumed to be timelike. Although the actual functional form of $ n_F $ is not essential, let us take 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 \,, \\ $


      $ \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}' \,, \\ $


      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, let us find 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 arbitrary spin vector $ \theta^ \mu $:

      $\nabla_ \mu \beta_ \nu + \nabla_ \nu \beta_ \mu = 0, $


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


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


      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}, $


      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 the finite Riemann curvature or the 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]} $. Besides, 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 \,, $


      $ \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}' $


      with $ \tilde{\omega}^{\mu\nu} = \dfrac{1}{2}\epsilon^{\mu\nu\rho\sigma}\omega_{\rho\sigma} $. On top of $ {\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) $. Then the average spin polarization per particle after freeze-out is given by

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


      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}}) $ with $ {{p}} $ 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}, $


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


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


      Unlike the massive case, the first conditions has an arbitrary function $ \phi (x) $, which arises due to the conformal invariance in the massless case, so $ \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 \,, $


      $ \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'\,, $


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


      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 at the rest frame of the fluid, $ E_p = |{{p}}| $. Note that spin polarization induced by the thermal vorticity and the electromagnetic field take the same form for both the massless and massive cases at global equilibrium, up to the difference in the on-shell conditions. Also the results are independent of choice of the frame vector $ n^ \mu $, as they should be.

    6.   Summary and outlook
    • In this paper, we derive the collisionless covariant spin kinetic theory at $ O(\hbar) $ for Dirac fermions in curved spacetime and external electromagnetic field. We start with deriving the dynamical equation for each Clifford component of the Winger function up to $ O(\hbar^2) $. We discuss the physical meaning of each such dynamical equation. We then pick up $ {\cal V}^ \mu $ and $ {\cal A}^ \mu $ as independent dynamical variables and derive two evolution equations for massless fermions (44), and four evolution equations for massive fermions (55) and (56), respectively. We introduce a timelike 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 the spin polarization from the approach of the kinetic theory. We derive the global equilibrium conditions from the kinetic equations, and find that the finite Riemann curvature or the external electromagnetic field is necessary to determine the spin-thermal vorticity coupling. At global equilibrium, we given the expression of the spin polarization induced by the electromagnetic field and the thermal vorticity, which is consistent with the results in the early literatures. We also give expressions for the 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 an material with the temperature gradient, which is described as the 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 in Kyoto University, where this work was developed during 'Quantum kinetic theories in magnetic and vortical fields'.

    Appendix A: Derivation of chiral anomaly
    • 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

      $\tag{A1} \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} $


      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 $ is given by $ p_0 = \pm | { p}| $, with which the delta function is reduced to

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


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

      $\tag{A3} \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} $


      Then the integral in Eq. (A1) is cast into

      $\tag{A4} \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} $


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

      $\tag{A5} \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} $


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

      $ \tag{A6} \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} $


      where we utilize

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


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

    Appendix B: General solutions at $ O(\hbar) $
    • We parametrize the perturbative solutions as

      $\tag{B1} {\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) \,. $


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

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


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


      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) $:

      $ \tag{B4} \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} $


      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 get Eqs. (36) and (37).

    Appendix C: Elimination of $ n^\mu $
    • In this Appendix, 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 $ \Delta_ \alpha $ to Eq. (9), we derive

      $\tag{C1} \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} $


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

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


      Combining the above two equations, we find

      $\tag{C3} 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 \,. $


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

      $\tag{C4} \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} $


      In order to reduce the above equation, we utilize the Schouten identity:

      $\tag{C5} \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} $


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

      $ \tag{C6} {\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] \,. $


      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).

    Appendix D: Derivation of Eqs. (55) and (56)
    • Here, we derive the kinetic equations for $ f_{\pm} $. Substituting Eq. (52) into Eq. (4), one obtains the following equation:

      $ \tag{D1} \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} $


      Besides, contracting Eq. (C3) with $ \theta^ \mu $ and inserting Eqs. (52) and (53), we get

      $\tag{D2} \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} $


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

    Appendix E: Global equilibrium condition from kinetic theory
    • In the massless case, the discussion of the global equilibrium conditions (65)-(67) was already given in Ref. [13]. Following the 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) $. In the following we verify that the condition (68) also fulfills Eq. (56) under the conditions (65)-(67), as follows. The leading order of $ f_A^{\rm{LE}} $ is 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

      $\tag{E1} \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} $


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

      $ \tag{E2} \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 \,, $


      which stems from Eq. (68). One finds that the above three terms in Eq. (E1) totally vanishes, as follows:

      $ \tag{E3} \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} $


      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} $.

    Appendix F: Global equilibrium condition from density operator
    • We discuss the global equilibrium using the maximum entropy principle, following Refs. [62-64]. The density operator for local equilibrium state is written as

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


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

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


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

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


      The global equilibrium is given by the condition 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

      $ \tag{F4} \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} $


      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.

Reference (64)



DownLoad:  Full-Size Img  PowerPoint