Vortical effects in Dirac fluids with vector, chiral and helical charges

Helicity of free massless Dirac fermions is a conserved, Lorentz-invariant quantity at the level of the classical equations of motion. For a generic ensemble consisting of particles and antiparticles, the helical and chiral charges are different conserved quantities. The flow of helicity can be modelled by the helicity current, which is again conserved in the absence of interactions. Similar to the axial vortical effect which generates an axial (chiral) current, the helicity current is induced by vorticity in a finite temperature medium with vector (electrical) charge imbalance via the helical vortical effects, leading to new nondissipative transport phenomena. These phenomena lead to the appearance of a new hydrodynamic excitation, the helical vortical wave. Our results suggest the existence of a new type of triangle anomalies in QED which involve the helicity currents in addition to the standard vector and axial currents. Further exploiting the conservation of the helical current, we show that a finite helical chemical potential may be used to characterise thermodynamic ensembles of fermions similarly to, but independently of, the vector charge and chirality. We derive the pressure $P$ for fermions at finite vector, axial and helical chemical potentials and show that the quantities arising in anomalous transport, including various vortical and circular conductivities and the shear-stress coefficients, can be obtained by differentiation of $P$ with respect to the appropriate chemical potentials. Finally, we calculate the helicity relaxation time in the quark-gluon plasma above the crossover and show that it is similar to that for the axial charge.


Introduction
Massless or nearly-massless fermions appear in many areas of physics, including theories of fundamental interactions, cosmological models of the early Universe, ultra-hot relativistic plasmas, and superfluids [1].Many relativistic phenomena are now available for experimental verification in recently discovered crystals of Dirac and Weyl semimetals, where the massless fermions appear as quasiparticle excitations [2].The most important properties of these excitations are usually associated with their vector (gauge) and axial (chiral) symmetries that affect, in the case of semimetals, electromagnetic [3], thermal [4], and elastic [5] responses of these materials.Many of the unusual features of these semimetals are associated with the quantum anomaly that breaks the continuous axial symmetry of an underlying classical theory [6].Similar anomalies lead to exotic transport phenomena of quarks, mediated by the topology of evolving gluon fields in expanding quark-gluon plasma (QGP) of heavy-ion collisions (HIC) [7].
In addition to the vector and axial charges, massless fermions can be characterised via a third, well-known, and, simultaneously, oftenforgotten quantity: the helicity.The fermionic helicity is sometimes confused with the chirality, even though these quantities reflect different physical properties of fermions [8].To highlight the importance of helicity, we demonstrate the existence of a set of new transport phenomena emerging in a gas of rotating massless fermions, the Helical Vortical Effects (HVEs), that differ substantially from their chiral counterparts, the Chiral Vortical Effects (CVEs) [9][10][11][12][13][14].
Hydrodynamics of relativistic plasmas with nonzero vorticity has attracted significant attention [15][16][17][18][19].The HVEs may see their applications in noncentral ultrarelativistic HIC that create a nearly perfect fluid of QGP, the most vortical fluid ever known [20].Exploiting the conservation of the helicity current J µ H for free fermions, we model helicity imbalance by means of the helicity chemical potential µ H , alongside the more familiar vector and axial chemical potentials, µ V and µ A , respectively.At a classical (non-quantum) level, we consider an ensemble of fermions at finite temperature and chemical potentials µ distributed according to the Fermi-Dirac distribution.Within this thermodynamic framework, we take derivatives of the pressure P with respect to µ A and show that the resulting quantities can be related to the vortical transport coefficients (e.g., the vortical and circular conductivities σ ω and σ τ , respectively) appearing in anomalous transport.We then employ finite-temperature quantum field theory to compute the expectation values of the charge currents J µ of vector ( = V ), axial ( = A), and helical ( = H) charges and of the stressenergy tensor T µν at finite temperature T , finite chemical potentials (including finite µ H ) and finite rotation parameter Ω.We show that all "classical" quantities receive quantum corrections that are of higher order with respect to the angular frequency Ω.
This paper is organised as follows.In Sec. 2, the vector (J µ V ), axial (J µ A ) and helical (J µ H ) charge currents are introduced as a triad corresponding to three Abelian symmetry groups of the classical Dirac Lagrangian.Also here, their operator form after second quantisation and their properties under the C, P, and T transformations are derived and the chemical potentials µ V /A/H associated to the conserved charge operators Q V /A/H are introduced.Section 3 introduces the classical description of fermions at finite T and µ , based on the Fermi-Dirac distribution.The properties of thermal states under rigid rotation at finite chemical potentials are discussed in Sec. 4. The full details of the computation are presented in Appendix A. We comment on the possible implications for new types of triangle anomalies involving the helical vertex (Helical triangle anomalies), the helicity relaxation time in the strongly-interacting quark-gluon plasma, new wave-like excitations in the helical-vortical matter (Helical vortical waves) and new vortices in Dirac fluids in Section 5. Section 6 concludes this paper.

U (1) H group of helical transformations and helicity current
We consider one species of free massless Dirac fermions in a flat (3 + 1)d Minkowski spacetime, with the Lagrangian where ψ = ψ † γ 0 is the Dirac adjoint of the 4spinor ψ.We take the 4 × 4 gamma matrices γ µ (µ = 0, . . .3) in the Dirac representation.
The classical Dirac Lagrangian ( 1) is known to be invariant under the vector U (1) V global Abelian symmetry: which allows the vector current to be derived via Noether's theorem, satisfying ∂ µ J µ V = 0 when ψ obeys the Dirac equation, where / ∂ = γ µ ∂ µ is the Feynman slash notation.Another well-known symmetry is related to the U (1) A symmetry of the Lagrangian (1) under chiral transformations, where γ 5 = iγ 0 γ 1 γ 2 γ 3 is the fifth gamma matrix.The U (1) A symmetry is also global and it can be understood by splitting ψ = ψ R + ψ L into its right and left chiral parts, Substituting the above decomposition into Eq.( 1), it can be seen that where f ← → ∂ µ g = f ∂ µ g − (∂ µ f )g is the bilateral derivative.In analogy to the vector current J µ V introduced in Eq. ( 3), one can introduce the axial current J µ A via which corresponds to a conserved axial charge, ∂ µ J µ A = 0, provided ψ satisfies the Dirac equation, Eq. ( 4).
The free theory (1) is also invariant under a "helical" U (1) H symmetry: (9) where the helicity operator, is determined via the spin operator s i = 1 2 ε 0ijk Σ jk which is a part of the covariant spin tensor Σ µν = i 4 [γ µ , γ ν ].The Hamiltonian of free fermions (1) is: In order to understand the U (1) H transformations in Eq. ( 9), let us consider a decomposition of ψ = ψ ↑ + ψ ↓ into its positive and negative helicity components, in a manner equivalent to the one in Eq. ( 6): We now seek to obtain a decomposition of the Lagrangian (1) similar to that in Eq. (7), As opposed to the chiral case, the above "helical" Lagrangian differs from the Lagrangian L (1) via We will now show that the above contribution can be interpreted as a 3-divergence which cancels when considering the action S = d 4 x L.
For this purpose, it is convenient to introduce the Fourier transform of the spinor ψ with respect to the spatial coordinates, Using Eq. ( 10), the positive and negative helicity parts are Inserting now the above into ( 14) and integrating with respect to d 3 x, we arrive at where the bilateral derivative ← → ∂ t acts on the spinors ψ p and ψp .Writing the spin matrix as it can be seen that both γ 0 and γ • p commute with with s•p/p.Further noting that [(s•p)/p] 2 = 1/4, it is straightforward to show that the factor between ψp and ψ p appearing in Eq. ( 17) vanishes separately for each value of λ, leading to Coming back to the U (1) H transformations, these can be applied on the spinor ψ at the level of the decomposition in Eq. ( 12), Denoting by L = L H + δL H the form of the Lagrangian after the U (1) H transformation, it can be seen that L H = L H , while Under integration with respect to d 3 x, δL H makes vanishing contributions and thus We are thus lead to interpret the U (1) H transformations as non-local.Nevertheless, the symmetry implied by Eq. ( 22) allows us to introduce the helicity current J µ H , which satisfies ∂ µ J µ H = 0 when ψ is a solution of the Dirac equation ( 4).Importantly, we remind here that the helicity current J µ H remains conserved when the fermions have finite mass m, as discussed in Ref. [49].Our paper is devoted to the massless case only.

Lorentz covariance of the helicity current J µ H
We now discuss the covariance of the helicity current under Lorentz transformations, using Wigner's method of induced representations, as discussed in Ref. [50] (see Chapter 2 therein).As is well-known, the helicity of a massive Dirac particle depends on the reference frame in which it is measured.For example, let us consider a positivehelicity particle travelling along the positive z axis, with four-momentum p µ = (E, 0, 0, p) µ .To an observer also travelling along the positive z axis with velocity V > p/E, the particle will appear to move towards the negative z axis.At the same time, in the rest frame of such an observer, the spin of the particle will remain aligned along the positive z axis, since the spin projection along the z axis does not change under boosts along the same axis.Therefore, the observer will measure a negative helicity in their reference frame.
In contrast, the helicity of a massless particle is Lorentz-invariant.This becomes manifest when considering the transformation properties of the Dirac spinor under a Lorentz transformation Λ: where D(Λ) = e −iS αβ ω αβ is the unitary representation of the Lorentz group, while S αβ = i 4 [γ α , γ β ] are the spin part of the generators of the Lorentz transformations for the Dirac field.For a one-particle momentum eigenstate U p,λ (x), satisfying one can show that [50] where θ is a real number defined by In the above, W (Λ, p) represents an element of the little group ISO(2) of the massless Dirac field, while ) is a pure boost transforming the null vector k µ = (k, 0, 0, k) µ into p µ characterized by the spherical coordinates (p, θ p , ϕ p ).The transformation S(α, β) and the rotation R 3 (θ) are given explicitly by where ζ = 1 2 (α 2 + β 2 ).Let us now consider the transformation properties of the helicity current.In a reference frame related to the laboratory frame by a Lorentz transformation Λ, the helicity current is given by (29) To see how J Λ; µ H (x) is connected to the laboratory-frame helicity current J µ H (x), we first replace ψ Λ (x) as shown in Eq. (24): (30) where h.c.denotes the Hermitian conjugate and we used the property D(Λ)γ µ D(Λ) = Λ µ ν γ ν .In order to evaluate D(Λ)h x D(Λ)ψ(Λ −1 x), we consider the expansion where V p,λ (x) = iγ 2 U * p,λ (x) represent the antiparticle states, satisfying When applying D(Λ)h x D(Λ) on ψ(Λ −1 x), Eq. ( 26) can be employed to replace where we used the property where the phase factor θ(Λ −1 , Λp) corresponds to the transformation (35) The above inverse can be computed as the matrix W (Λ −1 , Λp) can be put in the form (27) with parameters while θ(Λ −1 , Λp) = −θ(Λ, p).Using this latter result together with Eq. (34) into Eq.(33), we arrive at Noting that a similar relation holds also for the antiparticle states V p,λ (Λ −1 x), we find With the above relation, we can conclude that for massless particles, J µ H (x) transforms covariantly as a four-vector under the Lorentz group transformations:

Second quantisation and the Q V /A/H charge operators
Since ∂ µ J µ V /A/H = 0 for the classical solutions of the free-fermion theory (1), the vector-axialhelical charges, form a "triad" of the classically conserved U (1) quantities of the massless Dirac fermions [22].The Dirac equation (4) admits a complete set of solutions comprised of the particle and antiparticle modes, U j and V j = iγ 2 U * j , where j cumulatively indexes the eigenmodes.The commutation relations [γ 5 , h] = [H, γ 5 ] = [H, h] = 0 indicate that these modes are simultaneous eigenfunctions of the Hamiltonian, chirality and helicity operators: where E j = p j is the (positive) mode energy.The relation between the chirality χ j = ±1 and helicity λ j = ± 1 2 can be established using Eq. ( 10): The field operator is where canonical anticommutation rules for the particle bj and antiparticle dj operators are assumed (we use hats to denote operators acting on Fock space).Since we consider E j > 0, we are essentially working with the maximally-symmetric (Minkowski) vacuum, corresponding to an inertial Lorentz frame.Other choices are possible, for example the rotating vacuum introduced by Iyer [51], in which the modes with positive co-rotating energy E j = E j − Ωm j correspond to particle modes (m j is the total angular momentum along the axis rotation).The rotation parameter Ω will be introduced below.
The operators of the conserved charges (42) are: : : where the colons denote Wick (normal) ordering.The particle and antiparticle states contribute differently to the axial (46b) and helicity (46c) charges.By virtue of Eq. ( 44), Eqs. ( 43) and ( 46) also imply that the chiralities and helicities are indeed equal (opposite) to each other for particle (antiparticle) modes.

Discrete symmetries of the J µ V /A/H charge currents
In this section, the transformation properties of the helicity charge current (HCC) operator under the charge conjugation C, parity P, and time reversal T operations are derived, following the conventions of Ref. [23].In order to introduce the helicity charge current operator, J µ H (x), we first note that the classical quantity 2hψ(x) can be decomposed as After second quantisation, the equivalent relation can be obtained by taking the commutator of the quantum operator ψ with the helicity charge operator, Q H , defined in Eq. (46c): Based on the classical definition (23) of the helicity charge operator, J µ H , its quantum form can be introduced as follows: The above definition is exactly mirrored in the case of the vector charge current: where the commutators of the field operator and its adjoint with the vector charge operator are trivially The CPT symmetries of the HCC can be described in close analogy to those of the VCC.Specifically, we are interested in computing the following quantities: where the charge conjugation ( C), parity ( P) and time reversal ( T ) operators are assumed to be unitary.Furthermore, in the Dirac theory, these operators have the following action on the one-particle operators [23]: where η C , η P , η T , e −iθ b,j and e iθ d,j are irrelevant phases.For the set of eigenvalues j ≡ (E j , k j , m j , λ j ), the parity transformation gives j P ≡ (E j , −k j , m j , −λ j ), while the time reversal gives j T ≡ (E j , −k j , −m j , λ j ).
Let us now consider the action of a discrete symmetry, S, on the helicity and vector charge currents: The symmetry transformation can be applied to the commutator of two operators, A and B, using the following rule: Looking again at Eq. ( 54), it can be seen that the CPT properties of the HCC can be inferred from those of the VCC if the transformation properties of the corresponding charge operators, Q V and Q H , are known.It is not difficult to show that The properties of Q H differ from those of Q V due to the behaviour of the helicity, λ j , under the CPT transformations, as can be seen in the left half of Table 1.
Using Eq. ( 56), the CPT properties of the HCC can be inferred from those of the VCC, as follows: where x = (t, −x) when x = (t, x).When µ = i is a spatial index, J H i = − J i H , so that the temporal and spatial parts of the HCC are even under the T and P transformations, respectively, as can be seen in the right half of Table 1.

Helical chemical potential
The chirality and helicity are different quantities.The chirality of a particle is equal to its helicity (a right-chiral particle has a right-handed helicity) while the chirality of an antiparticle is opposite to its helicity (a right-chiral antiparticle has a left-handed helicity) [33].For a single fermion the vector (particle/antiparticle) and the axial (right-/left-chiral) charges determine uniquely its helicity.This is no longer the case for a fermion ensemble.
In a generic ensemble of Dirac fermions, all three classically conserved charges, Q = Q , can be expressed via linear combinations of four independent numbers that count particles (N R ↑ and N L ↓ ) and antiparticles ( N R ↓ and N L ↑ ).Let us attribute, for a moment, a "chemical potential" to each of the mentioned individual numbers (µ R ↑ is associated with N R ↑ , etc).The system can be described both via the V , A, H charges and alternatively, via the particle numbers: The individual particles contribute to each charge differently: Then Eq. (58) gives us the relations: that are readily understood.Usually, we focus on chirality and neglect helicity.Then the right-and left-handed chiral chemical potentials are given by the sum of particle and antiparticle contributions (60b , shown in the chiral-cone representation of the energy dispersion in Fig. 1(a).The axial chemical potential takes the familiar form: If we care about the helicity while being uninterested in the chirality, Eq. (60c) gives the right-and left-handed helical chemical potentials, respectively: These potentials can be used to label the occupation numbers of the "helical Weyl cones", thus characterizing the energy branches in a chirality-independent (and non-equivalent) manner, as shown in Fig. 1(b).The helical chemical potential takes the suggestive form: Physically, the property of helicity stands in between the vector and the axial characteristics of an ensemble of Dirac fermions: the helical charge (current) density is a difference between the axial charge (current) density of positive and negative energy eigenstates in the ensemble.
In practice, the similarities and differences between helicity and chirality [33] imply that an ensemble which contains particles only should possess a total chirality which exactly equals the total helicity of this particle-only ensemble.However, an ensemble of particles and antiparticles does not have the same chirality and helicity.Therefore, as both quantities are conserved in the massless theory, one needs to introduce separate chemical potentials for helicity and chirality to describe the grand canonical ensemble for this system.

Relativistic kinetic theory model
In this section, we consider an analysis of fluids possessing vector, axial and helical charge imbalance under the assumption of global thermodynamic equilibrium.In Subsec.3.1, we introduce Fermi-Dirac distribution functions to model the four types of populations, corresponding to particles and antiparticles of positive and negative helicity, serving as starting point for the computation of the pressure P .We further introduce derivatives of the pressure of the form ∂ n+1 P/∂ n µ A ∂µ , where the case n = 0 corresponds to the familiar charge densities Q , while at n > 0 we uncover the vortical, circular and shear conductivities σ ω , σ τ and σ Π appearing in anomalous transport.In Subsec.3.2, we discuss the kinematics of rigid motion, which will prove useful when discussing the vortical effects and corrections revealed by means of a full quantum field theoretical treatment in Sec. 4.

Classical V/A/H fluids
Let us consider an effective kinetic theory description of a fluid comprised of spin 1/2 particles and their corresponding antiparticles, described by a set of four distribution functions, where p is the particle four-momentum (p 2 = 0 for massless constituents), σ = ±1 distinguishes between particles (+1) and antipartices (−1), λ = ±1/2 denotes the polarization and u µ is the local velocity of the fluid.The V /A/H chemical potentials appearing in Eq. ( 61) are contained in the three-vector µ = (µ V , µ A , µ H ) and q σ,λ represent the eigenvalues of the : Q V /A/H : charge operators, namely Under the assumption of global thermodynamic equilibrium, the local temperature and chemical potentials are related to the macroscopic fluid motion via where Γ = u 0 is the local Lorentz factor of the fluid, such that the constants T 0 and µ ; 0 represent the values of the temperature and chemical potentials where the fluid is at rest.Using the distribution in Eq. ( 61), the stressenergy tensor (SET) and charge currents can be computed via where dP = d 3 p/[p 0 (2π) 3 ] and the subscript "cl" indicates classical quantities.Due to the form of f eq σ,λ , it is easy to see that both T µν cl and J µ cl take the perfect fluid form, where the energy density for massless particles satisfies E cl = 3P cl , while the pressure P cl and charge densities Q ; cl are Taking into account the expression (61) of f eq; σ p,λ , it can be seen that the charge densities can be obtained from the pressure via As we will show later in this paper, the higher order derivatives of P with respect to the axial chemical potential µ A can be related to the conductivities arising in anomalous transport.In particular, we identify the vortical and circular conductivities, σ ω ; cl and σ τ ; cl , as well as the shear conductivities σ Π ; cl as follows: We now compute the pressure by noting that since both chirality and helicity are frameindependent properties for massless particles, the Lorentz boost to the local rest frame can be performed as usual: (69) The integration with respect to p can be performed in terms of the polylogarithm function Li [53] and the results for ; cl , σ τ ; cl and Π ; cl can be written compactly as Li 4 (−e α σ,λ ), where Li 1 (−e α ) = − ln(1 + e α ) and Li 0 (−e α ) = −(1 + e −α ) −1 .The notation α σ,λ = q σ,λ • µ/T was introduced for brevity.In the above expressions, we have affixed the label "cl" to quantities which are inherently quantum, in order to distinguish the above "thermodynamic" contributions from quantum corrections which will be discussed in the following section (see also the Appendix).It is worthwhile remarking that the computation of σ ω ; cl and σ Π ; cl involves the product q A σ,λ q σ,λ , which is equal to 1 when = A. When ∈ {V, H}, this product becomes In the large temperature limit when µ /T 1, the pressure can be expanded as where The O(T −1 ) terms vanish when either one of µ V , µ A and µ H vanishes.The charge densities can be written explicitly as The above relations neglect O(T −1 ) contributions to Q ; cl which are proportional to µ ).The vortical conductivities can be expanded at high temperatures as The circular conductivities satisfy Finally, the shear stress conductivities σ Π ; cl satisfy Within the classical relativistic kinetic theory, the relations (65) giving the stress-energy tensor and the charge currents remain unchanged regardless of the kinematic state of the fluid, such that T µν cl and J µ cl; always retain their perfect fluid form.In Sec. 4, we will obtain in the frame of quantum field theory the corrections to the perfect fluid forms that lead to anomalous transport.The classical equilibrium quantities P cl and Q ; cl , as well as the coefficients σ ω ; cl , σ τ ; cl and σ Π ; cl , are related to thermodynamic integrals involving the Fermi-Dirac equilibrium distribution for fermions with the V , A and H charges and will thus serve as the basis to express all terms appearing as quantum corrections.

Kinematics of rigid motion
To reveal the role of the helical potential, we consider a gas of Dirac fermions at a finite temperature, uniformly rotating about the axis z.It Behaviour of the Levi-Civita symbol ε µναβ and of the spatial parts of the elements of the kinematic tetrad (velocity u, acceleration a, vorticity ω and fourth vector τ ) under the C-, P -, and T -inversions.The signs +/− indicate the even/odd nature of these quantities under the corresponding discrete transformations.
is convenient to introduce a particular basis of kinematic vectors for the rigid motion with the four-velocity where Ω is the angular velocity and Γ is the Lorentz factor.We use cylindrical coordinates (t, ρ, ϕ, z).The acceleration a µ = u ν ∇ ν u µ and vorticity ω µ = 1 2 ε µνλσ u ν ∇ λ u σ four-vectors are: The fourth vector, which is orthogonal to u, a and ω, is The properties of the elements of the kinematic tetrad under the C−, P − and T − inversions are summarised in Table 2.The squares of the vectors comprising the kinematic tetrad are while u 2 = 1.

Quantum field theoretical analysis
In this section, we discuss the properties of the expectation values of the stress-energy tensor T µν (SET) and of the vector, axial and helical charge currents J µ V /A/H , computed in a global equilibrium at finite temperature T , finite chemical potentials µ V /A/H and finite rotation Ω.The method employed for the computation relies on thermal mode sums and is summarised in Subsec.4.1.The results leading to quantum corrections to the classical results obtained in Sec. 3, as well as the new terms deviating from the perfect fluid form in Eq. ( 65) that give rise to vortical effects are discussed in Subsec.4.2.The computational details are given in Appendix A.

Thermal mode sums
The basis for the construction of quantum thermal states under rotation is the thermal average A of an operator A, defined as: where Z = Tr(ˆ ) is the partition function and the trace runs over Fock space [23][24][25][26].We employ the density operator ˆ given by where H is the Hamiltonian (11), M z = −i∂ ϕ +S z is the z component of the total angular momentum operator, while β 0 ≡ 1/T 0 and µ ,0 are, respectively, the values of the inverse temperature and the chemical potentials at the rotation axis ρ = 0 [27].
We now seek to evaluate the trace in Eq. (77) using the mode sum decomposition of the field operator ψ given in Eq. (45).Using the fact that [ H, M z ] = [ H, P z ] = 0, the eigenmodes of H, γ 5 and h introduced in Eq. ( 43) can be taken to be, simultaneously, the eigenmodes of M z and P z : where m j = ± 1 2 , ± 3 2 , . . ., takes only half-oddinteger values.These mode solutions have been previously derived in Ref. [28] and are reproduced, for convenience, in Eq. (123) of Appendix A.
The decompositions of the charge operators in Eq. ( 46), together with equivalent ones for H and M z given below, allow the following relations to be derived: where E j = E j − Ωm j is the co-rotating energy.
Let us apply this formalism for the computation of the thermal expectation values (t.e.v.s) of the charge current and stress-energy tensor operators, Using the convention A = : A : and the expansion (45) of the field operator with respect to the particle and antiparticle modes U j and V j , Eq. ( 84) can be used to obtain The sesquilinear forms J µ (ψ, χ) and T µν (ψ, χ) are and their explicit expressions are given in Eqs. ( 126) and (127).
The sum over j appearing in Eq. ( 85) is in principle sensitive to the choice of vacuum state.For the construction of rigidly-rotating thermal states, Iyer [51] argued that the modes with positive corotating energy E j = E j − Ωm j > 0 should be interpreted as particle modes.With this convention, the modes with E j E j < 0 will have opposite interpretation as compared to the stationary vacuum case.For example, a mode with E j < 0 and E j > 0 will represent a particle with respect to the rotating vacuum and an anti-particle with respect to the stationary vacuum.The difference in results obtained with respect to the rotating and stationary vacua is independent of the temperature T and chemical potentials µ of the medium and can thus be easily identifiable at the end of the calculation.We therefore present in this section the discussion in the simpler case of the computation with respect to the stationary vacuum, when the summation over j can be written as The results for the non-vanishing components of the charge current can be summarised as: where the summation runs over m = ± 1 2 , ± 3 2 , . . ., σ = ±1 and λ = ± 1 2 .The functions J * m ≡ J * m (qρ) ( * ∈ {+, −, ×}) are quadratic with respect to the Bessel functions, The components of the stress-energy tensor are while T ρρ = T zz .

Vortical effects in V/A/H fluids
The rigidly rotating gas of Dirac fermions generates the vector, axial, and helical 4-currents ( = V, A, H) according to Eq. ( 88), which can be rearranged as follows: along the four-vectors u µ , ω µ and τ µ .The radial components along a µ are absent for all currents (91).Henceforth, we work in the β (thermometer) frame, by fixing the four-velocity u µ to be equal to the one given in Eq. ( 73) [29][30][31].Combining Tables 1 and 2, the CPT parities of σ ω and σ τ can be obtained, as summarised in Table 3.
Looking now at the thermodynamic relations proposed in Eq. ( 68), it can be seen that the CPT properties of σ ω and σ τ should be the same as those of Q Q A and Q , respectively, which is easily confirmed by comparing Tables 1 and 3. Table 3 highlights, in particular, that helicity and chirality are different quantities that should not be confused with each other [33].
Springer Nature 2021 L A T E X template tor, chiral and helical charges Table 3 CPT properties of the vortical and circular conductivities σ ω and σ τ appearing in Eq. ( 91).The shear charge conductivities σ Π are also shown for completeness.The last four columns show the vortical and circular heat conductivities σ ω ε and σ τ ε , together with the shear coefficients Π 1 and Π 2 appearing in Eq. ( 93).The T -even nature of all transport coefficients in this table highlights the dissipationless nature of the transport phenomena, including, in particular, the new transport related to the helical degrees of freedom.
The thermal expectation value of the stress energy tensor T µν admits the general decomposition with respect to u µ as follows: where P is the isotropic pressure and E = 3P is the energy density.Both the anisotropic stress Π µν and the heat flux W µ are orthogonal to u µ and in addition, Π µν is required to be traceless.The most general decomposition of Π µν and W µ compatible with the tensor structure of T µν exhibited in Eq. ( 90) is The properties of the transport coefficients σ ω ε , σ τ ε , Π 1 and Π 2 under the C−, P − and T − inversions, summarised in Table 3, can be obtained based on those of the kinematic vectors u µ , a µ , ω µ and τ µ , given in Table 2.
We now report the results for the terms appearing in the decompositions of the charge currents and of the stress-energy current: where the quantities bearing the "cl" subscript are defined in Eq. ( 70) [their high temperature expansion can be found in Eq. ( 72)].The terms with the "(h)" superscript represent higher-order, subleading corrections, summarised in Appendix A, along with further computational details.The above expressions are fully consistent (and indeed, constrained) by the CPT properties revealed in Tables 1 and 3.These restrictions are even more important in determining the coefficients appearing on the SET sector.For example, the pressure has even parity (+ + +) with respect to all discrete transformations and can therefore receive corrections proportional only to σ ω A; cl and σ Π A; cl .Since σ ω ε has units of (energy) 3 and CPT parity equal to (+−+), its leading order contribution must be proportional (in fact, equal) to Q A; cl , while its first correction of order (energy) 1 is necessarily proportional to σ τ A; cl .Conversely, since σ τ ε has dimensions (energy) 2 and parity (+ + +), its leading order contribution is proportional to σ ω A; cl , while its dimensionless correction is proportional to σ Π A; cl , just like the leading order contribution to Π 1 .Finally, the dimensionality and (+−+) parity of Π 2 indicate that its leading order term must be related to σ τ A; cl .The vortical transport effects (91) and (93) are consistent with the C-, P -, and T -symmetries of the vector, axial and helical currents and charges, as shown in Table 1.These are dissipationless effects because the laws (91)-( 93) are even under the T -inversion.
On the rotation axis (ρ = 0), the circular currents vanish and the currents (91) and heat flux (93) point exactly along the vorticity Ω.The rotating dense (charged) Dirac matter generates on the rotation axis the helical current J H = σ ω H ω that is linearly proportional to the vector chemical potential µ V and temperature T (72c).On the other hand, the neutral Dirac matter with nonzero helicity (µ H = 0) generates the vector (charge) current J V = σ ω V ω.Remarkably, the mentioned helical terms, linearly proportional to a chemical potential and temperature, are allowed for the helical effects and, at the same time, are forbidden for the chiral phenomena by virtue of the C-, P -, and T -symmetries, as pointed out above.

Applications and discussion
Let us now consider how the helical degrees of freedom reveal themselves in an interacting field theory.In the previous section, we pointed out that finite vorticity gives rise to transport related to the helical current and the corresponding helical chemical potential.This allows us to speculate that in the underlying (fundamental) interacting theory (such as QED), triangle diagrams involving the helical vertex may be anomalous similarly to the well-known anomalous triangular diagrams which incorporate the axial current.We address this issue in Subsec.5.1.In Subsec.5.3, we point out that wave-like excitations involving both vector and helical currents appear in neutral Dirac materials in the presence of vorticity.Finally, we discuss in Subsec.5.4 potential implications of dynamical helical degrees of freedom for condensed matter systems by pointing out the possibility of whirlpool structures emerging in Dirac fluids.

Helical anomalies in QED
Fermions with electric charge e couple to electromagnetism through a source term eJ µ V A µ added to the Lagrangian (1): 2) is unbroken at the quantum level, thus reflecting the fundamental property of the electric charge conservation in quantum theory.The axial symmetry (5), however, is broken at the quantum level, leading to the nonconservation of J µ A in the presence of the electric E and magnetic B background fields [6]: where While elastic processes due to the electromagnetic interaction conserve both chirality and helicity [23,33,34], the inelastic processes of the form pair annihilation → pair creation may lead to helicity non-conservation.For example, in the process e + R e − L → e + L e − R , the axial charge of the initial and final states is 0, thus being conserved.The helicity charge of the initial state is −2 (here e − L is a particle with negative helicity), while in the final state, it is +2.We will discuss the efficiency of such processes in destroying helicity imbalance in the next subsection.
It is known that the vector σ ω V and axial σ ω A vortical conductivities (72c) at µ H = 0 are determined by the axial quantum anomalies [35].For example, the µ V µ A term in σ ω V is generated by the axial-vector-vector (AV V ) vertex of the axial anomaly (95), which is also responsible for the µ 2 V term in σ ω A .Both these terms share similar prefactors with the axial anomaly (95).The axial conductivity σ ω A contains also the µ 2 A term due to the axial-axial-axial (AAA) triangular anomaly, as well as the T 2 term which originates from the axial-graviton-graviton (AT T ) vertex of the mixed axial-gravitational anomaly [36].
The presence of the helical component in the vortical conductivities (72c) strongly suggests the existence of new types of triangle anomalies which involve the helicity vertex (10).The new helical anomalies must reveal themselves in the background of a "helical vector field" A H µ which couples with the Dirac fermions via the source term A H µ J µ H added to the Lagrangian (1).For instance, the leading µ H T (µ V T ) terms in σ ω V (σ ω H ), as well as the µ V µ H T term in σ ω ε , could have their origin in the new triangle HT V anomaly involving vector "V " (γ µ ), helical "H" (γ µ h), and graviton "T " ( T µν ) vertices.The quadratic µ H dependence (72c) of σ ω A implies the existence of a particular form of the mixed axial-helical anomaly responsible for the nonconservation of the axial current via a AHH vertex that shares similarity with the standard AV V vertex of the axial anomaly (95), to be discussed in Ref. [37].

Helicity relaxation time
In QED, helicity conservation is violated in the pair annihilation processes of the form e + R e − L → e + L e − R .The differential cross section corresponding to such helicity-violating pair-annihilation (HVPA) processes is [33] where α is the fine-structure constant.The centerof-momentum energy E cm and scattering angle θ cm are related to the Mandelstam variables s and t via where (p, k) and (p , k ) are the incoming and outgoing momenta, respectively.In the QGP, the relevant fermionic degrees of freedom are the (massless) quarks and the HVPA processes take place via gluon exchange.Eq. ( 96) is modified by replacing α by α QCD and by including the SU (3) generators from the gluon vertices [33]: where t a ij (1 ≤ a ≤ 8) are the SU(3) generators in the fundamental representation, normalized such that tr(t a t b ) = 1 2 δ ab , while (i, j) and (i , j ) are the colour indices of the initial and final quarks.Note that in Eq. (98), we are neglecting higherorder diagrams such as gluon emission or gluon exchange between the outgoing quarks.
Let us estimate the helicity relaxation time corresponding to an infinitesimal but small helical density.The Boltzmann equation for the distribution function of fermions (σ = 1) / anti-fermions (σ = −1) with on-shell four-momentum p µ = (E p , p) and helicity λ reads where C[f ] collects all processes in which such fermions are involved.The helical charge can be computed via Eq.( 66), reproduced below for convenience: where dP = d 3 p/[(2π) 3 E p ], u µ is the local plasma four-velocity, while σ and λ take the values ±1 and ±1/2, respectively.The degeneracy factor g = N c N f accounts for N c = 3 QCD colours and N f quark flavours.
The time derivative of Q H in a spatiallyhomogeneous system at rest will thus be given by the helicity-violating pair-annihilation (HVPA) processes contribution to C[f ] [39,60], where where the f σ p,λ = 1 − f σ p,λ factors originate from the Pauli blocking for fermions.The first (second) term in the square bracket is the gain (loss) term.The variables (p µ , k µ ) and (p µ , k µ ) denote the momenta of the outgoing (incoming) and incoming (outgoing) particles in the gain (loss) term.Focusing on the loss term, a quark of type σ with momentum p and helicity λ, having colour i and flavour f , annihilates with a quark of type −σ, with momentum k, colour j and opposite helicity −λ, but with the same flavour f .Under the HVPA process, the emerging quark with type σ has flipped helicity, −λ, colour i and flavour f .Its partner of type −σ will have helicity λ, colour j and the same flavour f .The arbitrariness of f is accounted for by the N f factor on the last line of Eq. (102), while the colours i , j and j can be summed starting from: Taking into account that tr(t b t a ) = 1 2 δ ab and that a t a t a = 4  3 I 3×3 , the above sums evaluate to 2/3, such that We consider now a near-equilibrium neutral plasma with a small helical imbalance, such that with f 0p = [e βEp +1] −1 being the Fermi-Dirac distribution for a neutral, unpolarized fluid.Ignoring quadratic terms in µ H , the sum over λ and σ can be performed in Eq. ( 101), leading to Writing with Q H = gµ H /3β 2 , the relaxation time τ H can be evaluated via The above integral can be evaluated in the CM frame (the details can be found in Appendix B), giving For two light fermionic flavors, N f = 2, the QGP in the strongly coupled (α QCD 1) regime above the crossover (T = 250 MeV) gives us, theoretically: On the other hand, estimations of the relaxation rate of the axial charge τ A from firstprinciple simulations [59] give at temperatures above the QCD crossover.An estimate based on the NJL model [60] gives the axial relaxation time in the range τ A (0.1 . . . 1) fm/c.The relaxation times τ A and τ H are lower than (or comparable with) the typical lifetime for the QGP, which amounts to several fm/c, indicating that non-elastic processes are rather effective in altering the helicity imbalance in the context of heavy ion collisions.Therefore, our crude estimations show that the helical charge decays is comparable or slower than that of the axial charge, τ H τ A , implying that both of them should be treated at the same level.

Helical vortical waves
The emergence of the helicity degree of freedom allows us to uncover new hydrodynamic excitations in the helical sector, that are similar to the chiral magnetic [40] and chiral vortical [15] waves.To illustrate this fact, we neglect background electromagnetic fields and impose the conservation of the charge currents J µ in Eq. ( 91) for the case of a rigidly-rotating fluid, when u µ , ω µ and τ µ are those given in Eqs. ( 73), ( 74) and (75), respectively.Noting that ∂ µ u µ = ∂ µ ω µ = ∂ µ τ µ = 0, the divergences ∂ µ J µ of the vector, axial and helical currents can be evaluated on the z axis (when ρ = 0), reducing to tor, chiral and helical charges where we took into account that neither the axial nor the helical charge remain conserved in a strongly-interacting plasma.We take a globally neutral plasma (Q = 0, = V, A, H) at finite temperature (T = 0), and consider the simplest excitation that propagates along the vorticity vector on the rotational axis, where the mean fluid velocity vanishes, v = 0.The bar over a symbol means a local thermodynamic average.We consider linear modes in hydrodynamic fluctuations.Neglecting the quantum corrections such that (Q , σ ω ) (Q ; cl , σ ω ; cl ), computed in Eqs.(72b)-(72c), we have Plugging the above into Eq.( 112), it can be seen that the fluctuations in the axial chemical potential, δµ A , are related to the fluctuations in temperature, δT , which are in principle governed by the hydrodynamic equations following from ∇ µ T µν = 0.For this reason, the chiral vortical wave (CVW) -that involves the vector and axial sectors -does not propagate in the neutral plasma [41][42][43][44].Unexpectedly, the vector and helical chemical potentials are linearly cross-coupled and therefore the helical vortical wave (HVW) does propagate.This can be seen by substituting (δQ V /H , σ ω V /H ) given in Eq. ( 113) into Eq.( 112), leading to where τ H is the helicity relaxation time discussed in Sec.5.2.Combining the above relations leads to the wave equation where the HVW velocity v HVW satisfies (we restore the fundamental constants , k B , and c) Let us critically assess the possibility of emergence of the HVW in the realistic environment of heavy-ion collisions.To estimate the velocity of the HVW in ultra-relativistic heavy-ion collisions, we take the temperature T 250 MeV above the pseudocritical QCD value [45,46], and the vorticity Ω 6.6 MeV 10 22 s −1 revealed in a RHIC experiment [20,48].We find that at these parameters, the HVW propagates with the velocity v HVW 1 × 10 −2 c, such that Looking now at the full dispersion relation for Fourier modes with angular velocity ω and wavenumber k, We see that the HVW propagates only for sufficiently high wavenumbers, kv HVW > 1/(2τ H ). For wavelength λ ≡ 2π/k, this condition transforms to the requirement λ < 4πl H implying the wavelength of the HVW should be smaller than 4π times the distance l H = v HVW τ H the wave propagates during the helicity-relaxation time τ H .In our case, l H 2.5 × 10 −2 fm, so that the propagating waves should have rather short wavelength λ 0.3 fm, or as a function of temperature, At T = 250 MeV, the lowest wavenumber bound is k 20 fm −1 .However, despite the shortness of the wavelength of the propagating helical vortical wave, its energy scale, given by the real part of the frequency (118), is a rather modest quantity, ε HVW ≡ kv HVW 0.2 c•fm −1 ≈ 40 MeV, owing to the slow propagation velocity v HVW .With respect to the local temperature, the HVW energy must satisfy so that the propagating HVW may be thermally excited in the QGP.The mode is not gapless in the sense that the lowest wavenumber is bounded from below.We note that the lower bound for the energy of the propagating HVW was computed based on our perturbative estimation of the helicity relaxation time.A more accurate estimate requires lattice QCD techniques [59] and lies beyond the scope of our work.
It is instructive to compare the parameters of the HVW (116) with its chiral analog [15,41]: To this end, we set for the chemical potential µ V = µ q 30 MeV [47], and obtain the velocity of the CVW, v CVW 1 × 10 −3 c, which falls into a range of the original estimation of Ref. [15].The HVW propagates much faster than the CVW: since T µ V in the low-density QGP being created in the HIC experiments.Given the mentioned equivalence of relaxation times for axial (τ A ) and helical (τ H ) charges, τ A ∼ τ H , and the large difference in their velocities (122), one can estimate that the wavelength of the propagating CVW will be much lower then the wavelength of the HVW.Therefore, the prospects of the HVW to exist in the QGP are substantially higher than the ones of the CVW since the HVW propagates faster.
We leave a more thorough analysis of the spectrum of waves allowed by the transport laws uncovered in this paper for a future analysis [37].

Vortices in Dirac fluids
In the condensed-matter context, the nontopological vortical structures may be created in finite samples of graphene [54].These arguments are supported by a strong experimental evidence [55].However, the Dirac excitations in two-dimensional graphene do not posses helicity, while similar vortex (or, whirlpool) structures in three-dimensional systems are yet to be realized experimentally.The two-and three-dimensional Dirac materials are somewhat comparable, in particular, due to similarities in typical Fermi velocities, as well as in the electronic viscosities in the suspected hydrodynamic regimes.Taking the vortex size l ∼ 1 µm [55] and the Fermi velocity v F c/300 10 6 m/s one gets the vorticity in the terahertz region, ω ∼ v F /l ∼ 10 12 s −1 .Free Dirac excitations, on top of the rotating viscous electronic fluid, should thus generate the helical vortex excitation which propagates, at temperature T = 100 K, with the noticeable velocity v HVW ∼ 10 4 m/s.The existence of the helical vortical wave in this environment depends, as in the case of the quark-gluon plasma, on the helical relaxation time, the estimation of which is beyond the scope of this paper.

Conclusions
In this paper, we exploited the invariance of the free Dirac Lagrangian under the U (1) H helical transformations to introduce the helicity charge current, J µ H , and the corresponding helicity charge operator, Q H .The associated helical chemical potential, µ H , can be used to describe helical imbalance in thermal states.Physically, the helical charge (current) density can be considered as the difference between the axial charge (current) density of positive and negative energy eigenmodes.
We found that the helical degree of freedom participates in new transport phenomena, the helical vortical effects (HVE), which appear in a rigidly rotating ensemble of Dirac fermions at finite temperature and density.The HVE's generate currents of vector, axial, and helical charges and lead to new excitations propagating through the helical vortical matter (the helical vortical waves, HVW).Our estimations show that, while the effects of the helical vortical waves are unlikely to be observable in quark-gluon plasma in highenergy physics, they could be probed in new Dirac materials in solid-state context.We suggested that the HVE may possibly be related to new triangle anomalies in quantum electrodynamics.In a thermodynamic context, the helical imbalance may lead to profound consequences for the chiral symmetry breaking in dense quark-gluon plasma [56].
Based on a classical kinetic theory description of Fermi-Dirac particles with vector, axial and helical imbalance, we showed that the transport coefficients related to anomalous transport can be obtained thermodynamically as derivatives of the charge densities Q = ∂P/∂µ with respect to the axial chemical potential µ A .Following the constraints imposed by the parities under the discrete C-, P -and T -inversions, we were able to express the higher-order corrections to the anomalous transport coefficients in terms of such thermodynamic functions.The exact proportionality coefficients are computed on the basis of thermal field theory for matter under rotation.
We close this paper by remarking that the new effects uncovered herein do not rely on the persistence of helicity imbalance, which is expected to be suppressed due to quantum interactions.Despite the fact that helicity is not conserved in realistic interacting gauge theories, we demonstrate that its relaxation time is of the same order as the chirality relaxation time, τ H τ A .Thus, both the helicity and the chirality should be treated on an equal footing.A relevant example of the application of the helicity current to polarization of hyperons in noncentral heavy-ion collisions is discussed in Ref. [57].

A Computational details
This Appendix summarises the details regarding the computation of the thermal expectation values (t.e.v.s) discussed in Sec. 4. The computational technique was introduced in Ref. [28] and is presented in more detail in the book chapter [58].
This appendix is structured as follows.Subsection A.1 presents details required for the derivation of the expressions for the non-vanishing components of the charge currents and stressenergy tensor, given in Eqs. ( 88) and (90) in the main text.A short mathematical interlude in Subsec.A.2 summarises the procedure employed to derive analytical expressions in a power series with respect to the rotation parameter Ω.The results relevant for the charge currents and stressenergy tensor are derived in the Subsections A.3 and A.4, respectively.All calculations presented in this section are performed with respect to the non-rotating (stationary) vacuum.

A.1 Modes and sesquilinear forms
For systems undergoing rigid rotation, it is convenient to consider a set of modes with respect to which the statistical operator ˆ in Eq. ( 78) is diagonal.For this purpose, we employ the helicity eigenmodes U j derived in cylindrical coordinates (t, ρ, ϕ, z) in Ref. [28], which are summarised for the case of massless fermions below: In the above, m = ± 1 2 , ± 3 2 , . . . is the eigenvalue of the total angular momentum operator M z , while λ = ± 1 2 is the eigenvalue of the helicity operator h.The eigenvalue of the linear momentum operator along the z axis is k, while q = p 2 − k 2 is the magnitude of the transverse momentum.The anti-particle modes are obtained via charge conjugation, V j = iγ 2 U * j .For the construction of rigidly-rotating thermal states, Iyer [51] argued that the modes with positive co-rotating energy E j = E j − Ωm j > 0 should be interpreted as particle modes, while E j = ±|p j | is allowed to take both positive and negative values.For brevity, in this paper we discuss only t.e.v.s computed with respect to the non-rotating (Minkowski) vacuum, keeping in mind that the difference between the two approaches manifests itself as terms which depend solely on the rotation parameter Ω (via the vorticity ω or acceleration a), being independent of the medium properties T and µ .With respect to the Minkowski vacuum, the Minkowski energy E j = p j is always positive.In this case, ψ can be decomposed as We now summarise the explicit expressions for the sesquilinear forms appearing in Eq. ( 86).Due to charge conjugation, the sesquilinear forms for the anti-particle modes V j = iγ 2 U * j are related to those for the particle modes via Furthermore, it can be seen that J µ A (U j , U j ) = J µ H (U j , U j ) = 2λ j J µ V (U j , U j ), since U j is a simultaneous eigenvector of both h and γ 5 .The sesquilinear forms J µ V (U j , U j ) read [49]: q j 8π 2 p j J × mj (q j ρ), while J ρ V /A/H (U j , U j ) = 0.The terms proportional to k j appearing in J t (U j , U j ) and J z (U j , U j ) were dropped in Eq. (88) since they are odd with respect to k j → −k j .
On the stress-energy tensor (SET) sector, we have [49] T tt (U j , U j ) = p j 8π 2 J + mj (q j ρ) + 2λ j k j p j J − mj (q j ρ) , T ρρ (U j , U j ) = q 2 j 8π 2 p j J + mj (q j ρ) − m j q j ρ J × mj (q j ρ) , T ϕϕ (U j , U j ) = q j m j 8π 2 ρ 3 p j J × mj (q j ρ), Upon substitution of the above results in Eq. ( 85), the terms which are odd with respect to k j → −k j can be dropped and Eq. ( 90) is reproduced.

A.2 Slow rotation limit
In order to derive analytical expressions for the t.e.v.s J µ ≡ : J µ : and T µν ≡ : T µν : , we consider the limit of slow rotation, in which the Fermi-Dirac factor n σ,λ (p) introduced in Eq. ( 83) can be expanded as follows: Note that the tilde no longer appears over p in the argument of n σ,λ on the right hand side.In order to illustrate the procedure, in what follows we take as an example the t component of the charge current J t .Substituting Eq. (128) in Eq. (88) gives The sum over m typically involves terms which are quadratic in Bessel functions [see Eqs.(88) and (90)], which can be performed analytically using the following formulas [28] where n is a non-negative integer and it is understood that m m 2n+1 J + m = m m 2n J − m = m m 2n J × m = 0.The coefficient s + n,j appearing above can be obtained from (131) and vanishes when j > n.For small values of n − j ≥ 0, the first few coefficients are Performing the sum over m using Eq.(130) in Eq. (129) gives The integration with respect to k of polynomials in q can be performed using which in the case of Eq. (133) gives Taking into account that the summation range of j is between 0 and n, the sums with respect to n and j can be exchanged, such that n now takes values from j to ∞. Shifting down the summation end for n from j to 0 gives in the case of Eq. ( 135) where integration by parts was employed 2j times to get from d 2n+2j n σ,λ (p)/dp 2n+2j to d 2n n σ,λ (p)/dp 2n .The integration with respect to p can be performed by changing the derivatives of n σ,λ (p) with respect to p to derivatives with respect to α Applying the above to Eq. (136) gives The summation over j can be performed in terms of the Lorentz factor Γ, vorticity ω 2 = Ω 2 Γ 4 and acceleration a 2 = Ω 2 Γ 2 (Γ 2 − 1).For example, in the case n = 0 in Eq. ( 138), the summation yields while when n = 1, one obtains The summation over σ and λ appearing in Eq. ( 138) can be performed in terms of the thermodynamic quantities introduced in Eq. (70) when r+1−2n ≥ 0. This corresponds to the n = 0 and n = 1 terms in Eq. ( 138), which give where J t,(h) is a high-order correction given by In what follows, we will label the contributions to thermal expectation values involving polylogarithms of negative index using the (h) superscript, as above.To understand the properties of such terms, one may look at the large temperature behaviour of the polylogarithm, (141) Summing now over σ and λ, we obtain where S s = σ,λ α s σ,λ is given explicitly by When s = 1, the sum vanishes, while for odd values of s higher than 1, it is not difficult to see that the result must be of order s with respect to µ V , µ A , µ H and their products.Explicitly, we have Charge prefactors q σ,λ can be obtained, e.g., by taking derivatives of Eq. ( 142) with respect to µ : Li n+r (−e α σ,λ ).(145) In the case of J t,(h) , the leading order contribution is proportional to There are no lower order contributions since ∂S 0 /∂µ = 0, S 1 = 0 and ζ(−2) = 0. Thus, we conclude that J t,(h) = O(µ V µ A µ H Ω 4 /µ T 3 ).

A.3 V/A/H charge currents
We now apply the strategy outlined in Subsec.A.2 for the case of the charge current J µ .We present as an intermediate step the relation similar to that in Eq. ( 136) in the case of the non-vanishing components of J µ : Springer Nature 2021 L A T E X template tor, chiral and helical charges Starting from the decomposition in Eq. ( 91), the charge density Q = u µ J µ , circular conductivity σ τ = −τ µ J µ /τ 2 and vortical conductivity σ ω = −ω µ J µ /ω 2 can be calculated via Using Eqs.(147), we find where q A σ,λ = 2λ was used in the expression for σ ω .After performing the p integral using Eq. ( 137) and the summation with respect to j as indicated in Eq. ( 139), it can be seen that the n = 0 contribution to Q coincides with Q ; cl introduced in Eq. (70).For n > 0, quantum corrections are uncovered: where Q (h) is subleading with respect to the temperature, taking the form The leading order contribution to Q (h) comes from the n = 0 term, while terms at higher n are penalised due to the T −2n factors.
As can be seen in Eq. ( 149), the circular conductivity receives its leading order contribution from the n = 1 term and is in fact equal to the classical term σ τ ; cl introduced in Eq. (70), where the quantum correction σ τ,(h) , exhibits the same structure and same leading order contribution as the correction Q (h) .
The computation of the vortical conductivity involves the product q A σ,λ q σ,λ , discussed in Eq. (71).Using the expressions for σ ω ; cl and σ Π ; cl introduced in Eq. ( 70), it can be seen that The quantum correction proportional to σ Π ; cl makes a T, µ -independent contribution to σ ω A , when σ Π A; cl = 1/π 2 + O(T −3 ), as indicated in Eq. (72e).This contribution can be regarded as a vacuum term which appears due to the fact that the thermal expectation values discussed in this appendix are computed with respect to the static Minkowski vacuum.If the calculation is performed with respect to the rotating vacuum proposed by Iyer [51], this term no longer appears, i.e. (σ ω A ) rot = σ ω A − lim T,µ →0 (σ ω A ).No such contributions can be found for σ ω V /H .The correction σ ω,(h) behaves differently for = A than when = V /H.Specifically, we find σ The results reported in this section are consistent with results previously reported in the literature.Our results correspond to the case of free massless fermions at finite µ V , µ A and µ H .In Sections 6 and 7 of Ref. [49], similar results are obtained in the case of vanishing axial chemical potential (µ A = 0), but at finite fermion mass.The case of finite µ A but vanishing µ H was considered in Ref. [32] (see Tables 1 and 2

A.4 Stress-energy tensor
We now discuss the properties of the stress-energy tensor.The pressure P , circular and vortical heat conductivities σ τ ε and σ ω ε and shear stress coefficients Π 1 and Π 2 can be obtained via Based on the SET components displayed in Eq. (90), the formalism from the previous subsection can be applied to derive the pressure, The heat vortical conductivity σ ω ε is given by where Q A; cl and σ τ A; cl can be found in Eq. ( 70).The heat circular conductivity σ τ ε is given by The shear stress coefficient Π 1 is given by

(158b)
As in the case of the charge currents, the results reported here for the stress-energy tensor are in perfect agreement with the massless limit of those found in Section 8 of Ref. [49] for helical fermions at vanishing µ A , as well as with those reported in Tables 1 and 2 of Ref. [32] for chiral fermions at vanishing µ H .

B Evaluation of the helicity relaxation time
In this appendix, we provide the details for the evaluation of the helicity relaxation time, based on the integral in Eq. ( 108), which we rewrite as The integral I Ω introduced above can be evaluated in the center of mass (CM) frame, where (161) where E 2 cm = (p + k) 2 .The above delta function allows the dK integral to be performed, leading to the replacement (k 0 , k ) = (p 0 , −p ) inside the integrand.We further write dP = p 0 dp 0 dΩ p /(2π) 3 and perform the p 0 integral, which amounts to replacing k 0 , p 0 → E cm /2 everywhere.
The boost to the CM frame changes the arguments of the FD distributions that depend on p and k .The exponential e βp 0 ≡ e βp •u l involves the four-product between the momentum p µ and the four-velocity u µ l = (1, 0, 0, 0) of the plasma in the fluid rest frame.In the center of mass frame, we have Taking into account also the effect of the delta function discussed above, the arguments of the equilibrium distributions will change according to where x = 2p •(p+k)/(E cm |p+k|) represents the cosine of the angle between p (in the CM frame) and p + k (in the laboratory frame).Finally, the factor (1 − cos θ cm ) 2 related to the angle between p and p measured in the CM frame can be computed in a Lorentz-covariant way via In order to evaluate the four-product p•p , we must boost p into the center of mass frame: where p cm is given by Going now back to Eq. (159), we have We now parametrize k in terms of its length k, cos γ = k • p/(kp) and the azimuthal angle ϕ k .With this parametrization, the integrand in Eq. (159) becomes independent of the angular parametrization of p: We now isolate the momentum magnitude by introducing the coordinates such that 0 ≤ z < ∞ and −1 ≤ δ ≤ 1.With this parametrization, the expressions appearing in the distribution functions are modified to where we introduced the variable ξ via such that dξ = 1−δ 2 4ξ d cos γ.We then have

Fig. 1
Fig. 1 The dispersion relations Ep = ±|p| for Dirac fermions in (a) chiral and (b) helical basis.The chemical potentials µ R/L determine the occupation numbers for the right-/left-handed chiralities, while the chemical potentials µ ↑/↓ dictate the densities of the right-/left-handed helicities.The spin orientations s with respect to the fermion's momentum p are shown by the horizontal arrows for each energy branch.