# Helical Separation Effect and helical heat transport for Dirac fermions

Victor E Ambruş, Maxim N Chernodub

► To cite this version:

Victor E Ambruş, Maxim N Chernodub. Helical Separation Effect and helical heat transport for Dirac fermions. European Physical Journal C: Particles and Fields, 2024, 84 (3), pp.282. 10.1140/epjc/s10052-024-12636-6 . hal-04179325

HAL Id: hal-04179325

https://hal.science/hal-04179325v1

Submitted on 10 Jun 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Helical Separation Effect and helical heat transport for Dirac fermions

Victor E. Ambrus¹ and Maxim N. Chernodub²,¹

¹Department of Physics, West University of Timișoara,
Bd. Vasile Pârvan 4, Timișoara 300223, Romania

²Institut Denis Poisson, CNRS UMR 7013, Université de Tours, 37200 France

(Dated: August 4, 2023)

An ensemble of massless fermions can be characterized by its total helicity charge given by the sum of axial charges of particles minus the sum of axial charges of antiparticles. We show that charged massless fermions develop a dissipationless flow of helicity along the background magnetic field. We dub this transport phenomenon as the Helical Separation Effect (HSE). Contrary to its chiral cousin, the Chiral Separation Effect, the HSE produces the helical current in a neutral plasma in which all chemical potentials vanish. In addition, we uncover the Helical Magnetic Heat Effect which generates a heat flux of Dirac fermions along the magnetic field in the presence of non-vanishing helical charge density. We also discuss possible hydrodynamic modes associated with the HSE in neutral plasma.

## I. INTRODUCTION

Massless or nearly massless fermions can be found in various areas of physics. They are present in theories of fundamental interactions and reveal themselves in cosmological models of the early Universe, ultra-hot relativistic plasmas, superfluids, and many other physical systems [1]. Discoveries of Dirac and Weyl semimetals, where the massless (chiral) fermions manifest as quasi-particle excitations, have opened up intriguing possibilities for experimental verification of relativistic phenomena associated with chiral fermions [2].

The pivotal properties of these excitations are typically linked to their vector (gauge), axial (chiral) and conformal (scale) symmetries. In the case of semimetals, these symmetries impact the electromagnetic [3, 4], thermal [5, 6], and elastic [7] responses of the materials. The peculiar characteristics of these semimetals are closely tied to the quantum anomaly that disrupts a continuous symmetry of an underlying classical theory due to quantum fluctuations [8]. In fundamental field theories, the vivid example is represented by the axial anomaly which gives rise to exotic transport phenomena of quarks, mediated by the evolving gluon fields' topology in the expanding quark-gluon plasma of heavy-ion collisions [9].

The continuous symmetries of massless fermions, leaving aside their conformal properties, are usually restricted to vector and axial symmetries, which lead to the conservation of vector (electric) and axial (chiral) charges. The excess of the number of particles over antiparticles, regardless of their spins, gives the vector charge. The sign of the spin projection on the momentum of a fermion or an anti-fermion, regardless of their vector charges, determines the chirality of the (anti-)particle in question. However, an ensemble of massless fermions can be characterized by a third, familiar, and, at the same time, intermittently disregarded quantity: helicity.

The similarity of fermionic helicity to fermion chirality often leads to the identification of these quantities, even though they represent distinct physical properties

of fermions [10]. The crucial difference – and the source of the persistent confusion – is rooted in the fact that the chirality of a single fermion identically coincides with its helicity while the chirality and helicity of a single anti-fermion are precisely opposite to each other. These properties lead to erroneous generic statements like “chirality is the same as helicity,” which is invalid for systems containing both particles and antiparticles. For example, a system made of a fermion and an anti-fermion can reside in a state with double total chirality and vanishing helicity (when the chiralities of the fermion and anti-fermion are the same), or in another state in which the total chirality equals zero, while helicity is non-vanishing (when the chiralities of the fermion and anti-fermion are opposite).

The helicity and chirality numbers coincide only for ensembles of massless fermions that consist only of particles with vector charges of the same sign. However, in various physically relevant environments – for example, in neutral finite-temperature plasmas or systems at high temperatures – this condition is not met, and the helicity should be treated on the same footing as chirality.

The vector, axial, and helical charges form a triad of quantities that characterize an ensemble of massless fermions and anti-fermions. The conservation of the vector $Q_V$, axial $Q_A$ and helical $Q_H$ charges in closed systems of free massless fermions allows us to introduce the thermodynamically conjugated chemical potentials that determine these charges in thermodynamic system. In particular, we introduce the helical chemical potential $\mu_H$ to encode the helicity imbalance, alongside the more familiar vector and axial chemical potentials, $\mu_V$ and $\mu_A$, respectively. As a result, the fermion ensemble is described by the Fermi-Dirac distribution at finite temperature and chemical potentials $\mu_\ell$ associated with the vector ($\ell = V$), axial ($\ell = A$), and helical ($\ell = H$) charges. The presence of the helical charge in interacting fermionic matter affects its thermodynamics, leading to a rich phase diagram [11, 12], in clear distinction with the thermodynamic effects of the axial charge density [13, 14]. Moreover, in addition to the standard vector
$j_V$ and axial $j_A$ currents, the helical degree of freedom allows us to construct the helical current $j_H$ which determines how the helicity is transferred in the medium by fermions.

In our paper, we concentrate on anomalous transport phenomena associated with helical properties of massless fermions in the presence of the background magnetic field following our earlier study to the fermion helicity in vortical backgrounds [15–17]. The main result of this paper is represented by the following formula for the helical current, which we dub the Helical Separation Effect (HSE):

$$j_H = \frac{T \ln 2}{\pi^2} q \mathbf{B} \quad \text{for } \mu_V = \mu_A = \mu_H = 0, \quad (1)$$

which shows that in the background of a time-independent uniform magnetic field $B$, the neutral matter develops a helical current along the direction of the magnetic field. The associated helical magnetic conductivity is proportional to the temperature $T$, implying that the higher the temperature, the stronger the helical current. Equation (1) is invariant under time reversal transformation, $t \to -t$, implying that the generated helical current is a dissipationless quantity.<fnref n="1" />

As the helicity of massless fermions is similar to the zilch quantum number for photons [21], the HSE (1) is a counterpart of the Zilch Vortical Effect which shows how the polarization of photons is transferred in the rotating photon gas [22]. We have shown in Ref. [17] that the helicity and axial vortical effects play a complementary role that elegantly accounts for the discrepancy between the polarization of particles and antiparticles. Similarly, one can expect that the new HSE will also play a complementary role to the influence of the CSE on the polarization of particles in heavy-ion collisions discussed in Ref. [23].

The HSE (1) has a simple generalization to the case of dense matter characterized by the triad of the chemical potentials $\mu_\ell$ ($\ell = V, A, H$):

$$j_H = \frac{T \ln 2}{\pi^2} \kappa_\beta(\boldsymbol{\mu}) q \mathbf{B}, \quad (2)$$

where

$$\kappa_{\beta}(\boldsymbol{\mu}) = \frac{1}{4 \ln 2} \sum_{\substack{a,b,c=\pm 1 \\ abc=1}} \ln \left( 1 + e^{\beta(a\mu_V + b\mu_A + c\mu_H)} \right), \quad (3)$$

is a thermodynamic factor normalized as $\kappa_\beta(\mathbf{0}) = 1$.

The name “Helical Separation Effect” (1) is inspired by its similarity with the Chiral Separation Effect (CSE) [24]:

$$j_A = \frac{\mu_V}{2\pi^2} q \mathbf{B}. \quad (4)$$

The CSE implies that the dense matter with $\mu_V \neq 0$ develops the current of axial (chiral) charge $j_A$ along a background magnetic field $B$. Contrary to the CSE (4), the HSE generates the helical current (1) in the absence of matter imbalance (with all chemical potentials vanishing). Notice that, while the strength of the CSE (4) is determined by the axial anomaly [24], the relevance of this anomaly to the helical transport (1) is not certain due to the $\ln 2$ factor in (1). A relevant discussion on the existence of possible helical anomalies in the context of the helical vortical conductivity is given in Ref. [15, 16]. Finally, we mention that the Helical Separation Effect (1) should not be confused with the Helical Magnetic Effects which appear in the context of Weyl semimetals [25] or in hydrodynamic transport [26].

Another important result of our article is that the magnetic field in fermionic matter with nonzero helical density generates a heat current. Namely, the background magnetic field $B$ in the presence of a nonvanishing helical density encoded in the helical chemical potential $\mu_H$, produces a non-zero heat flux $j_\epsilon^i \equiv T^{0i}$ along the axis of the magnetic field:

$$j_{\epsilon} = \frac{\mu_{H} T \ln 2}{\pi^{2}} q \mathbf{B}, \quad (5)$$

where $T^{0i}$ is an off-diagonal component of the energy momentum tensor for massless Dirac fermions. The Helical Magnetic Heat Effect (5) represents the leading term in the strong field limit, $|qB| \gg T^2$, with higher-order corrections in chemical potentials of the order of $O(\mu^3)$ being neglected. Notice that the heat current (5) is unexpectedly proportional to the first power of temperature $T$, contrary to the $T^2$ heat current produced by the mixed axial-gravitational anomaly [27]. A similar unusual feature was also noted for the helical vortical effects which produce a heat current in vortical helical matter [16].

This paper is organized as follows. In Sec. II, we introduce the conserved helicity operator for massive Dirac fermions in a magnetic field. A disadvantage of considering massive fermions is that, in addition to the explicit breaking of chirality conservation, they have a frame-dependence in the definition of helicity. Therefore, while discussing the applications of the Helical Separation Effect at the end of the paper, we will turn to the massless fermionic limit, similar to the generic tactics applied to the chiral degrees of freedom. We start our discussion from the discussion of the formal definition of the helicity operator in magnetic field, Sec. II, which leads us to the construction of the common mode solutions for the helicity operator and the Hamiltonian in Sec. III. Using the derived eigenmodes, the finite-temperature field theory at finite helical density is developed in Sec. IV which

<footnote n="1">¹ Notice that magnetic fields can also generate dissipationless (persistent) electric current in small metallic rings via the Aharonov-Bohm effect [18–20]. This effect, however, has a different nature than the anomalous transport effects of the kind (1) since the persistent current in a ring (i) emerges in the transverse direction to the magnetic field and (ii) disappears in the thermodynamic limit, when the radius of the ring becomes infinite.</footnote>
brings us to the anomalous transport in Sec. V. We derive not only the Helical Separation Effect (1), but also the corrections of the helical charges to the chiral magnetic and chiral vortical conductivities as well as to the heat transport. The interplay of the fluctuations of the helical charge density and the energy density – supported by the Helical Separation Effect (1) and the Helical Magnetic Heat Effect (5) – could, in principle, lead to the appearance of new hydrodynamic excitations. In Sec. VI we analyse this possibility for a neutral plasma and show that the helical degree of freedom does not generate a new hydrodynamic wave if all chemical potentials vanish. Our conclusions are summarized in the last section. Throughout this paper, Planck units ($\hbar = c = k_B = 1$) and the (+ – –) metric signature are adopted.

## II. HELICITY OPERATOR IN MAGNETIC FIELD

The symmetric Lagrangian for the Dirac field $\psi$ with mass $M$ coupled to an external electromagnetic potential reads as follows:

$$ \mathcal{L} = \frac{i}{2} \overleftarrow{\psi} \overleftarrow{\partial} \psi - q \overline{\psi} A \psi - M \overline{\psi} \psi, \quad (6) $$

where we use the standard convention that $q = -e < 0$ is the (negative) electron charge. The Dirac equation (6) can be obtained using the Euler-Lagrange formalism for both $\psi$ and $\overline{\psi}$:

$$ \begin{aligned} (i \overleftarrow{\partial} - q \overline{A} - M) \psi &= 0, \\ \overline{\psi} (i \overleftarrow{\partial} + q \overline{A} + M) &= 0. \end{aligned} \quad (7) $$

The equation for $\psi$ can be cast in the form:

$$ i \frac{\partial \psi}{\partial t} = H \psi, \quad H = \gamma^0 (-i \boldsymbol{\gamma} \cdot \boldsymbol{\nabla} + q \overline{A} + M), \quad (8) $$

where $H$ is the Hamiltonian operator.

For the case of an electron in a magnetic field, the four-potential $A^\mu$ can be taken in the Coulomb gauge,

$$ A^\mu = (0, \mathbf{A}), \quad (9) $$

where $\mathbf{A} \equiv \mathbf{A}(\mathbf{x})$ is a time-independent vector field. In this case, the Hamiltonian becomes:

$$ \begin{aligned} H &= M\gamma^0 + \gamma^0\boldsymbol{\gamma} \cdot \boldsymbol{\pi} \\ &= M\gamma^0 + 2\gamma^5\mathbf{S} \cdot \boldsymbol{\pi}, \end{aligned} \quad (10) $$

where

$$ \boldsymbol{\pi} = -i \boldsymbol{\nabla} - q \mathbf{A} \quad (11) $$

is the generalized momentum of a charged particle in an electromagnetic field and $\mathbf{S} = \frac{1}{2}\gamma^5\gamma^0\boldsymbol{\gamma}$ is the spin matrix. Due to the time-independence of the Hamiltonian, $\partial_t H =$

0, we can take particle normal-mode solutions $U_j$ which are eigenfunctions of the operator $H$:

$$ H U_j = E_j U_j \Rightarrow U_j = e^{-i E_j t} \tilde{U}_j, \quad (12) $$

where $\tilde{U}_j \equiv \tilde{U}_j(\mathbf{x})$ is a time-independent spinor field and the cumulative symbol $j$ labels distinct solutions to be discussed in detail later. Similarly, the antiparticle solutions of Eq. (7) can be obtained by charge conjugation,

$$ V_j = i\gamma^2 U_j^* \Big|_{q \to -q} = e^{iE_j t} \tilde{V}_j, \quad (13) $$

with $\tilde{V}_j \equiv \tilde{V}_j(\mathbf{x}) = i\gamma^2 \tilde{U}_j^*$, $q \to -q$. Clearly, $HV_j = -E_j V_j$.

In order to describe the polarization degree of freedom, we employ the helicity operator, defined as the projection of the spin along the generalized momentum:

$$ h = \frac{\mathbf{S} \cdot \boldsymbol{\pi}}{\sqrt{H^2 - M^2}}. \quad (14) $$

Taking into account that $[\gamma^5, \mathbf{S}] = [\gamma^0, \mathbf{S}] = 0$, as well as $[H, \boldsymbol{\pi}] = 0$, it is not difficult to show that $[H, \mathbf{S} \cdot \boldsymbol{\pi}] = 0$, from which we directly conclude that $H$ commutes with the helicity operator $h$:

$$ [H, h] = 0. \quad (15) $$

Writing conversely:

$$ h = \frac{\gamma^5 H - M\gamma^0}{2\sqrt{H^2 - M^2}}, \quad (16) $$

it is easy to see that $h^2 = 1/4$.

We can thus take the Hamiltonian eigenfunctions $U_j$, $V_j$ to be simultaneously eigenfunctions of $h$, satisfying:

$$ hU_j = \lambda_j U_j, \quad hV_j = \lambda_j V_j, \quad (17) $$

with the helicity eigenvalues $\lambda_j = \pm 1/2$ such that $\lambda_j^2 = 1/4$. Note that, at the level of the classical theory, the helicity operator $h$ preserves its eigenvalue under charge conjugation, in contradistinction to the Hamiltonian $H$ or chirality $\gamma^5$ (in the case of massless fermions). Using Eq. (12), the eigenvalue relations (17) become as follows:

$$ h_j \tilde{U}_j = \lambda_j \tilde{U}_j, \quad h_j^c \tilde{V}_j = \lambda_j \tilde{V}_j, \quad (18) $$

where the helicity operator and its conjugated form reduce, respectively, to the following expressions:

$$ \begin{aligned} h_j &= \frac{1}{2\sqrt{E_j^2 - M^2}} \begin{pmatrix} 0 & E_j + M \\ E_j - M & 0 \end{pmatrix}, \\ h_j^c &= -\frac{1}{2\sqrt{E_j^2 - M^2}} \begin{pmatrix} 0 & E_j - M \\ E_j + M & 0 \end{pmatrix}. \end{aligned} \quad (19) $$

The above equations are automatically solved by taking:

$$ \tilde{U}_j = \frac{1}{\sqrt{2}} \begin{pmatrix} \mathcal{K}_j^+ \\ 2\lambda_j \mathcal{K}_j^- \end{pmatrix} \otimes \phi_j, \quad \tilde{V}_j = \frac{1}{\sqrt{2}} \begin{pmatrix} 2\lambda_j \mathcal{K}_j^- \\ -\mathcal{K}_j^+ \end{pmatrix} \otimes \phi_j^c, \quad (20) $$
where we introduced the notation

$$ \mathcal{K}_j^{\pm} = \sqrt{1 \pm \frac{M}{E_j}}, \qquad (21) $$

The Pauli two-spinors $\phi_j$ and $\phi_j^c$ satisfy

$$ \begin{aligned} \sigma \cdot \pi \phi_j &= 2\lambda_j \sqrt{E_j^2 - M^2} \phi_j, \\ \sigma \cdot \pi \phi_j^c &= 2\lambda_j \sqrt{E_j^2 - M^2} \phi_j^c, \end{aligned} \qquad (22) $$

being related through charge conjugation:

$$ \phi_j^c = i\sigma^2 \phi_j^* \Big|_{q \to -q} . \qquad (23) $$

The helicity operator (14) shares similarity with the zilch quantum number which – despite a seemingly non-local definition – is a well-defined quantity associated with the polarization of photons [21].

### III. MODE SOLUTIONS IN A CONSTANT MAGNETIC FIELD

In this section, we review the eigensystem of a charged spinor particle in the background of the uniform homogeneous magnetic field, with an emphasis on the role of the helical degrees of freedom.

#### A. Energy levels

We seek to describe a particle moving in a constant magnetic field, oriented along the z axis:

$$ B = B e_z, \qquad (24) $$

where $e_z$ is a unit vector pointing in the positive direction along the z axis and B is a constant real number. The above magnetic field can be implemented via the asymmetric four-potential

$$ A^{\mu} = B x \delta_{y}^{\mu}, \qquad (25) $$

giving rise to the Faraday tensor

$$ F_{\mu\nu} = \partial_{\mu} A_{\nu} - \partial_{\nu} A_{\mu} = -B(g_{\mu x}g_{\nu y} - g_{\mu y}g_{\nu x}). \qquad (26) $$

We note that the symmetric four-potential $A_{\text{sym}}^{\mu} = \frac{1}{2}B(x\delta_{y}^{\mu} - y\delta_{x}^{\mu})$ can be obtained from the four-potential $A^{\mu}$ in Eq. (25) via the gauge transformation $A_{\text{sym}}^{\mu} = A^{\mu} + \partial^{\mu} \Lambda$, with $\Lambda = -\frac{1}{2}xyB$. For simplicity, we will continue the discussion in the context of the asymmetric choice (25), bearing in mind that the choice of gauge has no impact on the expectation values that form the main objectives of our paper.

It is easy to check that when $A^{\mu}$ is given by Eq. (25), the Hamiltonian commutes with both $P^y = -i\partial_y$ and

$P^z = -i\partial_z$, such that $\phi_j$ and $\phi_j^c$ can be taken as eigenfunctions of these operators,

$$ \phi_j = \frac{e^{ip_j^y y + ip_j^z z}}{2\pi} \begin{pmatrix} f_j^+ \\ f_j^- \end{pmatrix}, \quad \phi_j^c = \frac{e^{-ip_j^y y - ip_j^z z}}{2\pi} \begin{pmatrix} f_j^{c;+} \\ f_j^{c;-} \end{pmatrix}, \qquad (27) $$

where $f_j^{\pm} \equiv f_j^{\pm}(x)$ are scalar functions depending only on the transverse spatial coordinate $x$. The charge conjugates $f_j^{c;\pm}$ are related to $f_j^{\pm}$ via

$$ f_j^{c;\pm} = \pm f_j^{\mp;*} \Big|_{q \to -q} . \qquad (28) $$

Focussing now on the particle modes only, substituting Eq. (27) into Eq. (22) leads to

$$ \begin{pmatrix} 2\lambda_j \sqrt{E_j^2 - M^2} - p_j^z & i\sqrt{2|qB|}(\partial_{\xi_j} - \frac{\sigma}{2}\xi_j) \\ i\sqrt{2|qB|}(\partial_{\xi_j} + \frac{\sigma}{2}\xi_j) & 2\lambda_j \sqrt{E_j^2 - M^2} + p_j^z \end{pmatrix} \begin{pmatrix} f_j^+ \\ f_j^- \end{pmatrix} = 0, \qquad (29) $$

where we introduced

$$ \xi_j = \sqrt{2|qB|} \left( x - \frac{\sigma p_j^y}{|qB|} \right), \quad \sigma = \text{sgn}(qB). \qquad (30) $$

With the above notation, the functions $f_j^{\pm}(\xi_j)$ satisfy the following differential equation:

$$ \begin{gathered} \left( \frac{\partial^2}{\partial \xi_j^2} - \frac{1}{4} \xi_j^2 + \nu_j^{\pm} + \frac{1}{2} \right) f_j^{\pm} = 0, \\ \nu_j^{\pm} = \frac{E_j^2 - M^2 - (p_j^z)^2}{2|qB|} - \frac{1 \mp \sigma}{2}. \end{gathered} \qquad (31) $$

As is well known, the above equation has regular solutions given in terms of the Hermite polynomials when $\nu_j^{\pm}$ is a non-vanishing integer:

$$ f_j^{\pm} = C_j^{\pm} e^{-\xi_j^2/4} H_{\nu_j^{\pm}}(\xi_j), \qquad (32) $$

where $C_j^{\pm}$ are integration constants. From the above discussion, it is easy to write the quantization relation for the energy:

$$ E_j^2 = M^2 + (p_j^z)^2 + 2n_j |qB|, \quad \nu_j^{\pm} = n_j - \frac{1 \mp \sigma}{2}. \qquad (33) $$

Using now the recurrence relations $H'_{\nu}(z) = \nu H_{\nu-1}(z)$ and $zH_{\nu}(z) = H_{\nu+1}(z)+\nu H_{\nu-1}(z)$, it is possible to derive the relation

$$ \begin{aligned} & \left( \frac{\partial}{\partial \xi_j} - \frac{1 \mp \sigma}{2} \xi_j \right) H_{\nu}(\xi_j) \\ & = \nu \frac{1 \pm \sigma}{2} H_{\nu-1}(\xi_j) - \frac{1 \mp \sigma}{2} H_{\nu+1}(\xi_j), \end{aligned} \qquad (34) $$
which implies that the coefficients $C_j^\pm$ satisfy the following matrix equation:

$$ \left( \begin{pmatrix} 2\lambda_j \sqrt{E_j^2 - M^2} - \sigma p_j^z & -i\sqrt{2|qB|} \\ in_j \sqrt{2|qB|} & 2\lambda_j \sqrt{E_j^2 - M^2} + \sigma p_j^z \end{pmatrix} \begin{pmatrix} C_j^\sigma \\ C_j^{-\sigma} \end{pmatrix} \right) = 0. \quad (35) $$

The above expression gives one relation between the two normalization constants $C_j^\pm$. The second relation allowing these constants to be determined is provided by the requirement of unit norm of the modes $U_j$ under the Dirac inner product,

$$ \begin{aligned} \langle U_j, U_{j'} \rangle &= \int d^3 x \overline{U}_j(x) \gamma^0 U_{j'}(x) = \delta(j, j') \\ &\equiv \delta_{n_j n_{j'}} \delta_{\lambda_j \lambda_{j'}} \delta(p_j^y - p_{j'}^y) \delta(p_j^z - p_{j'}^z). \end{aligned} \quad (36) $$

The integration with respect to $y$ and $z$ leads to the $\delta(p_j^y - p_{j'}^y)$ and $\delta(p_j^z - p_{j'}^z)$ factors, while the integration with respect to $x$ can be performed using the orthogonality relation for the Hermite polynomials,

$$ \int_{-\infty}^{\infty} d\xi_j e^{-\xi_j^2/2} H_{n_j}(\xi_j) H_{n_j'}(\xi_j) = \sqrt{2\pi n_j!} \delta_{n_j, n_j'}. \quad (37) $$

In the following subsections, we will discuss independently the cases $n_j = 0$ (the so-called lowest Landau level, LLL) and $n_j > 0$ (the higher Landau levels, HLL).

## B. Lowest Landau level (LLL)

On the lowest Landau level (LLL) corresponding to $n_j = 0$, we avoid the divergent solution with $\nu_j^{-\sigma} = -1$ by setting

$$ C_j^{-\sigma} \Big|_{n_j=0} = 0. \quad (38) $$

The only non-trivial solution corresponds to the polarization satisfying

$$ 2\sigma\lambda_j = \operatorname{sgn}(p_j^z), \quad (39) $$

where we remind the reader that $\sigma = \operatorname{sgn}(qB)$, cf. Eq. (30).

Normalizing now the modes using the Dirac inner product in Eq. (36), we find the normalization $C_j^\sigma = (|qB|/\pi)^{1/4}$ corresponding to the LLL mode:

$$ \begin{aligned} U_j(x) \Big|_{n_j=0} &= \theta(\sigma \lambda_j p_j^z) \frac{e^{-iE_j t}}{\sqrt{2}} \begin{pmatrix} K_j^+ \\ 2\lambda_j K_j^- \end{pmatrix} \otimes \phi_j^\sigma(\mathbf{x}) \Big|_{n_j=0}, \\ \phi_j^\sigma \Big|_{n_j=0}(\mathbf{x}) &= \frac{e^{ip_j^y y + ip_j^z z}}{2\pi} e^{-\frac{\xi_j^2}{4}} \left(\frac{|qB|}{\pi}\right)^{1/4} \begin{pmatrix} \frac{1+\sigma}{2} \\ \frac{1-\sigma}{2} \end{pmatrix}. \end{aligned} \quad (40) $$

Similarly, the antiparticle modes can be found as

$$ \begin{aligned} V_j(x) \Big|_{n_j=0} &= \theta(-\sigma \lambda_j p_j^z) \frac{e^{iE_j t}}{\sqrt{2}} \begin{pmatrix} 2\lambda_j K_j^- \\ -K_j^+ \end{pmatrix} \otimes \phi_j^{c;\sigma}(\mathbf{x}) \Big|_{n_j=0}, \\ \phi_j^{c;\sigma}(\mathbf{x}) \Big|_{n_j=0} &= \frac{e^{-ip_j^y y - ip_j^z z}}{2\pi} e^{-\frac{\xi_j^2}{4}} \left( \frac{|qB|}{\pi} \right)^{1/4} \begin{pmatrix} \frac{1+\sigma}{2} \\ -\frac{1-\sigma}{2} \end{pmatrix}, \end{aligned} \quad (41) $$

where $\xi_{j;c}$ is the charge conjugate of $\xi_j$,

$$ \xi_{j;c} = \xi_j \Big|_{q \to -q} = \sqrt{2|qB|} \left( x + \frac{\sigma p_j^y}{|qB|} \right). \quad (42) $$

These expressions allow us to establish a link between $\phi_j^{c;\sigma}$ and its charge conjugate as follows:

$$ \phi_{0,py,pz,\lambda}^{c;\sigma}(\mathbf{x}) = \sigma \phi_{0,-py,-pz,\lambda}^{\sigma}(\mathbf{x}). \quad (43) $$

## C. Higher energy levels

For $n_j > 0$, we have

$$ C_j^{-\sigma} = - \frac{in_j \sqrt{2|qB|}}{2\lambda_j \sqrt{E_j^2 - M^2} + \sigma p_j^z} C_j^{\sigma}, \quad (44) $$

which leads to

$$ \begin{aligned} C_j^\sigma &= \frac{\mathcal{N}_j}{\sqrt{n_j!}} \left( \frac{|qB|}{4\pi} \right)^{1/4} \mathfrak{p}_j^{2\sigma\lambda_j}, \\ C_j^{-\sigma} &= \frac{-i\mathcal{N}_j \sqrt{n_j}}{2\lambda_j \sqrt{n_j!}} \left( \frac{|qB|}{4\pi} \right)^{1/4} \mathfrak{p}_j^{-2\sigma\lambda_j}, \end{aligned} \quad (45) $$

where $\mathcal{N}_j$ is a phase factor, $|\mathcal{N}_j| = 1$, while

$$ \mathfrak{p}_j^{\pm} = \sqrt{1 \pm \frac{p_j^z}{\sqrt{E_j^2 - M^2}}}. \quad (46) $$

Specifically, we have

$$ U_j(x) = \frac{e^{-iE_j t}}{\sqrt{2}} \begin{pmatrix} K_j^+ \\ 2\lambda_j K_j^- \end{pmatrix} \otimes \phi_j^\sigma(\mathbf{x}), \quad (47a) $$

where the Pauli spinors are given by

$$ \begin{aligned} \phi_j^+(\mathbf{x}) &= \frac{e^{ip_j^y y + ip_j^z z}}{2\pi\sqrt{n_j!}} e^{-\xi_j^2/4} \left(\frac{|qB|}{4\pi}\right)^{1/4} \\ &\quad \times \begin{pmatrix} \mathfrak{p}_j^{2\lambda_j} H_{n_j}(\xi_j) \\ -2i\lambda_j \mathfrak{p}_j^{-2\lambda_j} \sqrt{n_j} H_{n_j-1}(\xi_j) \end{pmatrix}, \\ \phi_j^-(\mathbf{x}) &= \frac{e^{ip_j^y y + ip_j^z z}}{2\pi\sqrt{n_j!}} e^{-\xi_j^2/4} \left(\frac{|qB|}{4\pi}\right)^{1/4} \\ &\quad \times \begin{pmatrix} \mathfrak{p}_j^{2\lambda_j} \sqrt{n_j} H_{n_j-1}(\xi_j) \\ 2i\lambda_j \mathfrak{p}_j^{-2\lambda_j} H_{n_j}(\xi_j) \end{pmatrix}. \end{aligned} \quad (47b) $$
In the above, we used $N_j = 1$ for $\sigma = 1$ and $N_j = 2i\lambda_j$ for $\sigma = -1$. The antiparticle modes $V_j = i\gamma^2 U_j^*$ can be obtained through the charge conjugation operation, being given explicitly by

$$ 
\begin{aligned}
V_j(x) &= \frac{e^{iE_j t}}{\sqrt{2}} \begin{pmatrix} 2\lambda_j \mathcal{K}_j^- \\ -\mathcal{K}_j^+ \end{pmatrix} \otimes \phi_j^{c;\sigma}(\mathbf{x}), \\
\phi_j^{c;+}(\mathbf{x}) &= \frac{e^{-ip_j^y y - ip_j^z z}}{2\pi\sqrt{n_j!}} e^{-\xi_{j;c}^2/4} \left(\frac{|qB|}{4\pi}\right)^{1/4} \times \begin{pmatrix} -2i\lambda_j \mathfrak{p}_j^{-2\lambda_j} H_{n_j}(\xi_{j;c}) \\ -\mathfrak{p}_j^{2\lambda_j} \sqrt{n_j} H_{n_j-1}(\xi_{j;c}) \end{pmatrix}, \\
\phi_j^{c;-}(\mathbf{x}) &= \frac{e^{-ip_j^y y - ip_j^z z}}{2\pi\sqrt{n_j!}} e^{-\xi_{j;c}^2/4} \left(\frac{|qB|}{4\pi}\right)^{1/4} \times \begin{pmatrix} 2i\lambda_j \mathfrak{p}_j^{-2\lambda_j} \sqrt{n_j} H_{n_j-1}(\xi_{j;c}) \\ -\mathfrak{p}_j^{2\lambda_j} H_{n_j}(\xi_{j;c}) \end{pmatrix}.
\end{aligned} 
\quad (48a) $$

It can be seen that the charge-conjugated Pauli spinors $\phi_j^{c;\sigma}(\mathbf{x})$ corresponding to higher levels satisfy the relation

$$ \phi_{n,p^y,p^z,\lambda}^{c;\sigma}(\mathbf{x}) = -2i\lambda\sigma\phi_{n,-p^y,-p^z,\lambda}^{\sigma}(\mathbf{x}), \quad (48b) $$

which shares certain similarity with Eq. (43).

# IV. FINITE-TEMPERATURE FIELD THEORY AT FINITE HELICAL DENSITY

## A. Second quantization

The full solution of the Dirac equation reads

$$ \psi(x) = \sum_j [U_j(x)a_j + V_j(x)b_j^\dagger], \quad (49) $$

where the symbol $j$, which labels the eigenmodes, collectively denotes the eigenvalues $p_j^y, p_j^z, n_j$ and $\lambda_j$. The sum over the eigenmodes is shorthanded as follows:

$$ \sum_j = \sum_{\lambda_j = \pm \frac{1}{2}} \int_{-\infty}^{\infty} dp_j^y \int_{-\infty}^{\infty} dp_j^z \sum_{n_j=0}^{\infty}, \quad (50) $$

where the $n_j = 0$ term corresponds to the contribution coming from the lowest Landau level (LLL), which exists only when $2\sigma\lambda_j p_j^z > 0$.

Demanding that the equality (49) is exact and noting that $a_j = \langle U_j, \psi \rangle$ and $b_j^\dagger = \langle V_j, \psi \rangle$, we arrive to the completeness relation for the set of modes $U_j$ and $V_j$:

$$ \sum_j [U_j(t, \mathbf{x}) U_j^\dagger(t, \mathbf{x}') + V_j(t, \mathbf{x}) V_j^\dagger(t, \mathbf{x}')] = \delta^3(\mathbf{x} - \mathbf{x}'). \quad (51) $$

Promoting now $\psi(x)$ to a Fock-space operator $\hat{\psi}$, the canonical anti-commutation relation

$$ \{\hat{\psi}(t, \mathbf{x}), \hat{\psi}^\dagger(t, \mathbf{x}')\} = \delta^3(\mathbf{x} - \mathbf{x}'), \quad (52) $$

can be achieved by imposing

$$ \{\hat{a}_j^\dagger, \hat{a}_{j'}\} = \{\hat{b}_j^\dagger, \hat{b}_{j'}\} = \delta(j, j'), \quad (53) $$

where $\delta(j, j') = \delta(p_j^y - p_{j'}^y)\delta(p_j^z - p_{j'}^z)\delta_{n_j, n_{j'}}\delta_{\lambda_j, \lambda_{j'}}$, as in Eq. (36).

## B. Conserved charges

In this paper, we describe fermionic ensembles at finite temperature in the background of a constant homogeneous magnetic field. We will investigate states which exhibit imbalance with respect to the electric and helical charges, as well as with respect to the axial charge (in the case of massless fermions). These conserved charges can be obtained by considering the vector, axial and helical charge currents, defined, respectively, as follows:

$$ 
\begin{aligned}
J_V^\mu &= \bar{\psi}\gamma^\mu\psi, \\
J_A^\mu &= \bar{\psi}\gamma^\mu\gamma^5\psi, \\
J_H^\mu &= \bar{\psi}\gamma^\mu h\psi + \bar{h}\psi\gamma^\mu\psi.
\end{aligned} 
\quad (54) $$

We now compute the divergences of the above currents. In the case of the vector current, we have

$$ \partial_\mu J_V^\mu = \overline{\partial} \psi \psi + \overline{\psi} \partial \psi. \quad (55) $$

Using the Dirac equation (7), the derivatives $\overline{\partial}\psi$ and $\overline{\partial}\psi$ can be replaced via

$$ \overline{\partial}\psi = -i(q\mathcal{A} + M)\psi, \quad \overline{\partial}\psi = i\overline{\psi}(q\mathcal{A} + M). \quad (56) $$

Substituting the above into Eq. (55) gives the conservation of the vector current, $\partial_\mu J_V^\mu = 0$.

Employing the same steps as above for the derivative of the axial four-current $\partial_\mu J_A^\mu$ reveals the partial conservation of the axial current (PCAC),

$$ \partial_\mu J_A^\mu = 2M(i\overline{\psi}\gamma^5\psi), \quad (57) $$

where the expression between the parentheses on the right-hand side is the pseudoscalar condensate.

Finally, in the case of the helicity current, we have

$$ \partial_\mu J_H^\mu = \overline{\partial} \psi h \psi + \overline{\psi} \partial \psi h \psi + \text{h.c.}, \quad (58) $$

where h.c. denotes the Hermitian conjugate. Taking advantage of the commutation relation (15), it is easy to show that if $\psi$ is a solution of the Dirac equation, then $h\psi$ also satisfies the Dirac equation. This leads to the relation

$$ \overline{\psi} \partial \psi h \psi = -i \overline{\psi} (q \mathcal{A} + M) h \psi = - \overline{\partial} \psi h \psi, \quad (59) $$

implying that the helical current is identically conserved even for massive fermions:

$$ \partial_\mu J_H^\mu = 0. \quad (60) $$
Note however that the conservation of the helicity current is broken perturbatively by helicity-violating pair annihilation processes, as discussed in Ref. [16]. In a setting of real physical systems, such as the quark-gluon plasma, the non-conservation of helicity due to perturbative scattering has the same timescale as the non-conservation of chirality due to dynamical mass generation [28], which allows us to treat the helicity and chirality on equal basis [16].

The above discussions prompt us to introduce the total vector, axial and helical charges,

$$Q_{\ell} = \int d^3 x J_{\ell}^0, \quad \ell = V, A, H, \qquad (61)$$

which are identically conserved for the vector and helical degrees of freedom:

$$\partial_t Q_V = \partial_t Q_H = 0 \quad [\text{for fermions with any mass}], \qquad (62)$$

for both chiral and massive fermions. However, due to the PCAC property (57), the axial charge does not represent a conserved quantity for a massive particle:

$$\partial_t Q_A = 2M\eta, \quad \eta = \int d^3x i\bar{\psi}\gamma^5\psi. \qquad (63)$$

Thus, while at finite mass $M \neq 0$, the axial charge $Q_A$ is not a good conserved number, the vector and helical charges are identically conserved.

Thanks to the helicity eigenmode solutions discussed in Sec. III, both the vector and the helicity charge operators $\hat{Q}_{V/H}$ are diagonal and, together with the Hamiltonian $\hat{H}$, they admit the following Fock space representation:

$$\begin{aligned}
\hat{H} &:= \sum_j E_j (\hat{a}_j^\dagger \hat{a}_j + \hat{b}_j^\dagger \hat{b}_j), \\
\hat{Q}_V &:= \sum_j (\hat{a}_j^\dagger \hat{a}_j - \hat{b}_j^\dagger \hat{b}_j), \\
\hat{Q}_H &:= \sum_j 2\lambda_j (\hat{a}_j^\dagger \hat{a}_j - \hat{b}_j^\dagger \hat{b}_j),
\end{aligned} \qquad (64)$$

where the colons: $\hat{\mathcal{O}} := \hat{\mathcal{O}} - \langle 0 | \hat{\mathcal{O}} | 0 \rangle$ denote Wick (normal) ordering of an operator $\hat{\mathcal{O}}$. It is worth pointing out that, while $\hat{H}$ is even under charge conjugation (both particles and antiparticles contribute with the same sign to $\hat{H}$), both the vector and the helical charge operators are odd and therefore discriminate between particles and antiparticles. The CPT properties of $\hat{Q}_H$ are discussed in Ref. [15].

As a side remark, the Hamiltonian $H$ can be obtained from the Dirac energy-momentum tensor, $\Theta^{\mu\nu}$, defined via Noether's theorem as

$$\Theta^{\mu\nu} = \frac{\partial \mathcal{L}}{\partial(\partial_\mu \psi)} \partial^\nu \psi + \partial^\nu \bar{\psi} \frac{\partial \mathcal{L}}{\partial(\partial_\mu \bar{\psi})} - g^{\mu\nu} \mathcal{L}, \qquad (65)$$

where the Dirac Lagrangian $\mathcal{L}$ is given in Eq. (6). Taking into account that $\mathcal{L}$ vanishes when $\psi$ satisfies the

Dirac equation, we arrive at the following relation for the energy-momentum tensor:

$$\Theta^{\mu\nu} = \frac{i}{2} \bar{\psi} \gamma^{\mu} \overleftrightarrow{\partial}^{\nu} \psi. \qquad (66)$$

Due to the presence of interactions, $\Theta^{\mu\nu}$ is, in general, not conserved:

$$\partial_{\mu} \Theta^{\mu\nu} = q J_{V}^{\mu} \partial^{\nu} A_{\mu}. \qquad (67)$$

The right-hand side of the above expression is not gauge invariant. To restore gauge invariance, it is customary to alter $\Theta^{\mu\nu}$ by adding a gauge-fixing total divergence term, $\Delta\Theta^{\mu\nu} = \partial_\rho(F^{\mu\rho}A^\nu) = -qJ_V^\mu A^\nu + F^{\mu\rho}\partial_\rho A^\nu$, where we used Maxwell's equations, $\partial_\mu F^{\mu\nu} = qJ_V^\nu$. Adding the first term, $-qJ^\mu A^\nu$, to $\Theta^{\mu\nu}$ yields $\tilde{\Theta}^{\mu\nu} = \frac{i}{2}\gamma^\mu \overleftrightarrow{\partial}^\nu \psi - qJ^\mu A^\nu$. Its divergence, $\partial_\mu \tilde{\Theta}^{\mu\nu} = -qJ_{V;\mu}F^{\mu\nu}$, represents the relativistic Lorentz force acting on the electric charge and is a gauge-invariant quantity. The second term, $F^{\mu\rho}\partial_\rho A^\nu$, serves to make the electromagnetic contribution to the energy-momentum tensor, $\Theta^{\mu\nu}_{e.m.} = -F^{\mu\lambda}\partial^\nu A_\lambda$, as well as its divergence, gauge-invariant. However, the gauge-fixing term becomes position-dependent, since $A^\nu = Bx\delta_y^\nu$, as shown in Eq. (25). For this reason, we will not consider the gauge-fixing term and we will focus instead on the form of the energy-momentum tensor $\Theta^{\mu\nu}$ given in Eq. (66).

The static homogeneous magnetic field considered in this paper, $A_\mu = Bxg_{\mu y}$, does not produce work as the Lorentz force is always normal to the direction of the current. Therefore, for the energy component of the conservation relation (67), $\nu = 0$, we have $\partial_\mu \Theta^{\mu 0} = 0$, so that the energy is conserved, as expected. This property allows the Hamiltonian to be defined as the total conserved energy,

$$H = \int d^3x \Theta^{00}. \qquad (68)$$

Upgrading $H$ to the Fock space operator $\hat{H}$ and employing the explicit mode decomposition (49) gives the expression in Eq. (64).

Finally, in the case of chiral (massless) fermions, the axial charge is conserved, $\partial_t Q_A = 0$, and the corresponding Fock-space operator $\hat{Q}_A$ admits the decomposition

$$\hat{Q}_A := \sum_j 2\lambda_j (\hat{a}_j^\dagger \hat{a}_j + \hat{b}_j^\dagger \hat{b}_j) \quad [\text{for } M = 0]. \qquad (69)$$

Note that, contrary to the vector and helical charge operators $\hat{Q}_V$ and $\hat{Q}_H$, the axial charge operator $\hat{Q}_A$ is even with respect to charge conjugation (see Table 1 in Ref. [16]).

### C. Thermal expectation values

We now construct thermal states at finite vector and helical chemical potentials, using the statistical opera-
tor [29–31]

$$ \hat{\rho} = e^{-\beta(:\hat{H}: - :\mu\cdot\hat{Q}:)}, \quad (70) $$

where $\hat{H}$ is the Hamiltonian, while the conserved charges sector contains the vector and helical charges for massive fermions. In the case of massless fermions, we allow also for a chiral charge, so that we introduce a generic shift of the Fermi energy levels controlled by the whole triad of the corresponding chemical potentials:

$$ \mu \cdot \hat{Q} = \mu_V \hat{Q}_V + \mu_A \hat{Q}_A + \mu_H \hat{Q}_H. \quad (71) $$

One should bear in mind that the axial chemical potential $\mu_A$ vanishes for a massive particle, since the axial charge is not conserved if the particle mass is nonzero, cf. (63).

The statistical operator (70) describes a system at rest, having four-velocity

$$ u^{\mu} \partial_{\mu} = \partial_{t}. \quad (72) $$

With the above four-velocity vector, we can construct the covariant magnetic field four-vector,

$$ B^{\mu} = \frac{1}{2} \epsilon^{\mu\nu\alpha\beta} u_{\nu} F_{\alpha\beta} = B \delta_{z}^{\mu}, \quad (73) $$

where $\epsilon^{\mu\nu\alpha\beta}$ is the Levi-Civitta tensor and we took the sign convention such that $\epsilon^{0123} = 1$.

The thermal expectation value of an operator $\hat{A}$ is defined as

$$ \langle \hat{A} \rangle = \mathcal{Z}^{-1} \mathrm{tr}(\hat{\rho} \hat{A}), \quad \mathcal{Z} = \mathrm{tr}(\hat{\rho}). \quad (74) $$

In the following, we focus only on operators which are quadratic with respect to the field operator $\hat{\psi}$, thus allowing for the following decomposition:

$$ \hat{A} = \sum_{j,j'} \left[ \begin{aligned} & \mathcal{A}(U_j, U_{j'}) \hat{a}_j^\dagger \hat{a}_{j'} + \mathcal{A}(V_j, V_{j'}) \hat{b}_j \hat{b}_{j'}^\dagger \\ & + \mathcal{A}(U_j, V_{j'}) \hat{a}_j^\dagger \hat{b}_{j'}^\dagger + \mathcal{A}(V_j, U_{j'}) \hat{b}_j \hat{a}_{j'} \end{aligned} \right], \quad (75) $$

where $\mathcal{A}(\psi, \chi)$ is a sesquilinear form in the sense that

$$ \begin{aligned} \mathcal{A}(\alpha\psi + \beta\phi, \chi) &= \alpha^* \mathcal{A}(\psi, \chi) + \beta^* \mathcal{A}(\phi, \chi), \\ \mathcal{A}(\psi, \alpha\chi + \beta\phi) &= \alpha \mathcal{A}(\psi, \chi) + \beta \mathcal{A}(\psi, \phi). \end{aligned} \quad (76) $$

Evaluating $\langle \hat{A} \rangle$ essentially boils down to evaluating the thermal expectation values of quadratic products of the one-particle operators, which we discuss below.

First, we note that the decompositions (64) imply the following commutation relations:

$$ \begin{aligned} [\hat{H}, \hat{a}_j^\dagger] &= E_j \hat{a}_j^\dagger, & [\hat{H}, \hat{b}_j^\dagger] &= E_j \hat{b}_j^\dagger, \\ [\hat{Q}_V, \hat{a}_j^\dagger] &= \hat{a}_j^\dagger, & [\hat{Q}_V, \hat{b}_j^\dagger] &= -\hat{b}_j^\dagger, \\ [\hat{Q}_H, \hat{a}_j^\dagger] &= 2\lambda_j \hat{a}_j^\dagger, & [\hat{Q}_H, \hat{b}_j^\dagger] &= -2\lambda_j \hat{b}_j^\dagger. \end{aligned} \quad (77) $$

A fourth line can be added for massless fermions, namely

$$ M = 0: [\hat{Q}_A, \hat{a}_j^\dagger] = 2\lambda_j \hat{a}_j^\dagger, \quad [\hat{Q}_A, \hat{b}_j^\dagger] = 2\lambda_j \hat{b}_j^\dagger. \quad (78) $$

With the above relations, it can be established that

$$ \hat{\rho} \hat{a}_j^\dagger \hat{\rho}^{-1} = e^{-\beta \epsilon_j^+} \hat{a}_j^\dagger, \quad \hat{\rho} \hat{b}_j^\dagger \hat{\rho}^{-1} = e^{-\beta \epsilon_j^-} \hat{b}_j^\dagger, \quad (79) $$

where we introduced the shifted energy levels

$$ \mathcal{E}_j^{\zeta_j} = E_j - \zeta_j \mu_V - 2\lambda_j \mu_A - 2\zeta_j \lambda_j \mu_H, \quad (80) $$

where the quantum number $\zeta_j = \pm 1$ distinguishes particles and antiparticles:

$$ \zeta_j = \begin{cases} +1, & \text{for particles,} \\ -1, & \text{for antiparticles.} \end{cases} \quad (81) $$

Equation (80) takes into account the difference in the positions of the Fermi levels for particles with vector, axial, and helical charges which are anchored in the corresponding chemical potentials $\mu_V$, $\mu_A$, and $\mu_H$, respectively. Despite the seemingly complicated form of this equation, it has a simple meaning, consistent with Eq. (71): the vector chemical potential $\mu_V$ controls the excess of particles over antiparticles which, in turn, are distinguished by the quantum number $\zeta_j$. The axial chemical potential $\mu_A$ applies equally to particles and antiparticles but distinguishes the polarization states described by $\lambda_j$. Finally, the helical chemical potential $\mu_H$ is sensitive to the helicity of the particle which is opposite for particles and antiparticles, hence the product $\zeta_j \lambda_j$ in the last term of Eq. (80).

The thermal expectation value of the quadratic form $\hat{a}_j^\dagger \hat{a}_{j'}$ thus evaluates to

$$ \langle \hat{a}_j^\dagger \hat{a}_{j'} \rangle = e^{-\beta \epsilon_j^+} \langle \hat{a}_{j'} \hat{a}_j^\dagger \rangle = \frac{\delta(j, j')}{e^{\beta \epsilon_j^+} + 1}, \quad (82) $$

where on the last line we used $\langle \hat{a}_{j'} \hat{a}_j^\dagger \rangle = \delta(j', j) - \langle \hat{a}_j^\dagger \hat{a}_{j'} \rangle$. Similarly,

$$ \langle \hat{b}_j^\dagger \hat{b}_{j'} \rangle = \frac{\delta(j, j')}{e^{\beta \epsilon_j^-} + 1}. \quad (83) $$

Finally, it is not difficult to see that

$$ \langle \hat{a}_j^\dagger \hat{b}_{j'}^\dagger \rangle = \langle \hat{b}_j \hat{a}_{j'} \rangle = 0. \quad (84) $$

We are now in a position to evaluate the thermal expectation value of the operator $\hat{A}$ in Eq. (75). Focussing on thermal effects, we consider the expectation value of the normal-ordered operator : $\hat{A} := \hat{A} - \langle 0|\hat{A}|0 \rangle$, which evaluates to

$$ A \equiv \langle: \hat{A} :\rangle = \sum_j \left[ \frac{\mathcal{A}(U_j, U_j)}{e^{\beta \epsilon_j^+} + 1} - \frac{\mathcal{A}(V_j, V_j)}{e^{\beta \epsilon_j^-} + 1} \right]. \quad (85) $$
For simplicity, we will often consider separately the contributions due to the LLL ($n_j = 0$) and to the higher-energy modes ($n_j > 0$) by writing

$$A = A_0 + \sum_{n_j=1}^{\infty} A_{n_j}, \qquad (86)$$

where $A_0$ corresponds to the lowest Landau level (with $n_j = 0$). The terms $A_{n_j}$ are obtained after performing the $p_j^y$ and $p_j^z$ integrations, as well as the sum over polarization $\lambda_j$:

$$A_{n_j} = \sum_{\lambda_j} \int dp_j^y \int dp_j^z \left[ \frac{\mathcal{A}(U_j, U_j)}{e^{\beta \mathcal{E}_j^+} + 1} - \frac{\mathcal{A}(V_j, V_j)}{e^{\beta \mathcal{E}_j^-} + 1} \right]. \quad (87)$$

### D. Scalar and pseudoscalar condensates

Let us apply the previously-discussed formalism to the computation of the scalar and pseudoscalar condensates, $\bar{\psi}\psi$ and $i\bar{\psi}\gamma^5\psi$. In the latter case, it is not difficult to see that the mode solutions in Eqs. (40) and (47) give vanishing contributions to the pseudoscalar condensate,

$$i\overline{U}_j\gamma^5 U_j = 0, \qquad (88)$$

and similarly $i\overline{V}_j\gamma^5 V_j = 0$.

In the case of the scalar condensate, the antiparticle contribution can be obtained via charge conjugation, as follows:

$$\overline{V}_j V_j = -(\overline{U}_j U_j)^*_{q \to -q}. \qquad (89)$$

On the LLL, we have

$$\overline{U}_j U_j|_{n_j=0} = \frac{Me^{-\xi_j^2/2}}{4\pi^2 E_j} \sqrt{\frac{|qB|}{\pi}} \theta(\sigma \lambda_j p_j^z), \quad (90)$$

while for $n_j > 0$,

$$\overline{U}_j U_j|_{n_j>0} = \frac{Me^{-\xi_j^2/2}}{8\pi^2 E_j n_j!} \sqrt{\frac{|qB|}{\pi}} \times (\mathfrak{p}_{j;2\sigma\lambda_j}^2 H_{n_j}^2 + n_j \mathfrak{p}_{j;-2\sigma\lambda_j}^2 H_{n_j-1}^2), \quad (91)$$

where it is understood that the Hermite polynomials take the argument $\xi_j$. In order to evaluate the chiral condensate $\langle \bar{\psi}\psi \rangle$, the integration with respect to $p_j^y$ can be performed using the orthogonality relation (37) for the Hermite polynomials,

$$\int_{-\infty}^{\infty} dp^y e^{-\xi^2/2} H_n(\xi) H_m(\xi) = \sqrt{\pi |qB|} n! \delta_{nm}. \quad (92)$$

Specifically, we have

$$\overline{U}_j U_j = \frac{M |qB|}{4\pi^2 E_j} \begin{cases} \theta(\sigma \lambda_j p_j^z), & n_j = 0, \\ 1, & n_j > 0. \end{cases} \quad (93)$$

The final result reads

$$\bar{\psi}\psi = \frac{M|qB|}{4\pi^2} \sum_{n_j=0}^{\infty} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} \frac{dp_j^z g_j}{E_j(e^{\beta\mathcal{E}_j} + 1)}, \quad (94)$$

where the summation runs over the helicity quantum number $\lambda_j = \pm 1/2$, particle-antiparticle number $\zeta_j = \pm 1$ and $\mathcal{E}_j \equiv \mathcal{E}_j^\zeta$ was introduced in Eq. (80). The degeneracy factor evaluates to $g_j = 1$ for $n_j = 0$ and $g_j = 2$ for $n_j > 0$.

### E. Charge currents

The sesquilinear forms corresponding to the charge current operators read

$$\begin{aligned} \mathcal{J}_V^\mu(\psi, \chi) &= \bar{\psi}\gamma^\mu\chi, \\ \mathcal{J}_A^\mu(\psi, \chi) &= \bar{\psi}\gamma^\mu\gamma^5\chi, \\ \mathcal{J}_H^\mu(\psi, \chi) &= \bar{\psi}\gamma^\mu h\chi + \bar{h}\psi\gamma^\mu\chi. \end{aligned} \quad (95)$$

Replacing now $\psi = \chi = U_j$ or $\psi = \chi = V_j$ in terms of the helicity eigenmodes discussed in Sec. III, it can be seen that the sesquilinear forms for the helical charge current can be obtained from the ones for the vector charge current via

$$\begin{aligned} \mathcal{J}_H^\mu(U_j, U_j) &= 2\lambda_j \mathcal{J}_V^\mu(U_j, U_j), \\ \mathcal{J}_H^\mu(V_j, V_j) &= 2\lambda_j \mathcal{J}_V^\mu(V_j, V_j). \end{aligned} \quad (96)$$

The expectation values $\mathcal{J}_\ell^\mu(x) \equiv \langle : \hat{J}_\ell^\mu(x) : \rangle$ with $\ell \in \{V, A, H\}$ can be evaluated using Eq. (85), where the sesquilinear forms $\mathcal{J}_\ell^\mu(V_j, V_j)$ corresponding to the antiparticle modes can be obtained from those corresponding to the particle modes by using the following properties:

$$\begin{aligned} \overline{V}_j \gamma^\mu V_j &= (\overline{U}_j \gamma^\mu U_j)^*_{q \to -q}, \\ \overline{V}_j \gamma^\mu \gamma^5 V_j &= -(\overline{U}_j \gamma^\mu \gamma^5 U_j)^*_{q \to -q}. \end{aligned} \quad (97)$$

In the following, we will focus on the explicit computation of the sesquilinear forms corresponding to the particle modes, which we will denote in shorthand notation as

$$\mathcal{J}_{\ell; j}^{\mu} \equiv \mathcal{J}_{\ell}^{\mu}(U_j, U_j). \quad (98)$$

On the LLL, an explicit computation shows that

$$\phi_j^{\sigma; \dagger} \sigma^x \phi_j^{\sigma} \Big|_{n_j=0} = \phi_j^{\sigma; \dagger} \sigma^y \phi_j^{\sigma} \Big|_{n_j=0} = 0, \quad (99)$$

such that

$$\mathcal{J}_{\ell; j}^{x} \Big|_{n_j=0} = \mathcal{J}_{\ell; j}^{y} \Big|_{n_j=0} = 0. \quad (100)$$
The temporal and vertical components can be computed as follows:

$$
\begin{align}
\mathcal{J}_{V; j}^{t} |_{n_j=0} &= \frac{e^{-\xi_j^2/2}}{4\pi^2} \sqrt{\frac{|qB|}{\pi}} \theta(\sigma\lambda_j p_j^z), \nonumber \\
\mathcal{J}_{V; j}^{z} |_{n_j=0} &= \frac{2\lambda_j \sigma |p_j^z| e^{-\xi_j^2/2}}{4\pi^2 E_j} \sqrt{\frac{|qB|}{\pi}} \theta(\sigma\lambda_j p_j^z), \nonumber \\
\mathcal{J}_{A; j}^{t} |_{n_j=0} &= \frac{2\lambda_j |p_j^z| e^{-\xi_j^2/2}}{4\pi^2 E_j} \sqrt{\frac{|qB|}{\pi}} \theta(\sigma\lambda_j p_j^z), \nonumber \\
\mathcal{J}_{A; j}^{z} |_{n_j=0} &= \frac{\sigma e^{-\xi_j^2/2}}{4\pi^2} \sqrt{\frac{|qB|}{\pi}} \theta(\sigma\lambda_j p_j^z). \tag{101}
\end{align}
$$

Using Eq. (92), the $p_j^y$ integral can be evaluated, leading to

$$
\begin{align}
\begin{pmatrix} J_{V;0}^t \\ J_{A;0}^t \\ J_{H;0}^t \end{pmatrix} &= \frac{|qB|}{4\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z}{e^{\beta\mathcal{E}_j} + 1} \begin{pmatrix} \zeta_j \\ 2\lambda_j p_j^z/E_j \\ 2\lambda_j \zeta_j \end{pmatrix}, \nonumber \\
\begin{pmatrix} J_{V;0}^z \\ J_{A;0}^z \\ J_{H;0}^z \end{pmatrix} &= \frac{qB}{4\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z}{e^{\beta\mathcal{E}_j} + 1} \begin{pmatrix} 2\lambda_j p_j^z/E_j \\ \zeta_j \\ p_j^z/E_j \end{pmatrix}, \tag{102}
\end{align}
$$

where $E_j = \sqrt{M^2 + p_{z; j}^2}$ on the LLL.

At higher Landau levels with $n_j > 0$, we compute $\mathcal{J}_{\ell; j}^{\mu}$ as follows:

$$
\begin{align}
\mathcal{J}_{V; j}^{t} &= \frac{e^{-\xi_j^2/2}}{8\pi^2 n_j!} \sqrt{\frac{|qB|}{\pi}} (\mathfrak{p}_{2\sigma\lambda_j}^{j;2} H_{n_j}^2 + \mathfrak{p}_{-2\sigma\lambda_j}^{j;2} n_j H_{n_j-1}^2), \nonumber \\
\mathcal{J}_{V; j}^{y} &= -\frac{qB e^{-\xi_j^2/2}}{4\pi^2 E_j (n_j - 1)!} \sqrt{\frac{2}{\pi}} H_{n_j} H_{n_j-1}, \nonumber \\
\mathcal{J}_{A; j}^{z} &= \frac{\sigma e^{-\xi_j^2/2}}{8\pi^2 n_j!} \sqrt{\frac{|qB|}{\pi}} (\mathfrak{p}_{2\sigma\lambda_j}^{j;2} H_{n_j}^2 - \mathfrak{p}_{-2\sigma\lambda_j}^{j;2} n_j H_{n_j-1}^2), \nonumber \\
\begin{pmatrix} \mathcal{J}_{A; j}^{t} \\ \mathcal{J}_{A; j}^{y} \\ \mathcal{J}_{V; j}^{z} \end{pmatrix} &= 2\lambda_j \sqrt{1 - \frac{M^2}{E_j^2}} \begin{pmatrix} \mathcal{J}_{V; j}^{t} \\ \mathcal{J}_{V; j}^{y} \\ \mathcal{J}_{A; j}^{z} \end{pmatrix}, \tag{103}
\end{align}
$$

while $\mathcal{J}_{V; j}^{x} = \mathcal{J}_{A; j}^{x} = 0$. As before, the Hermite polynomials take the argument $\xi_j$. Performing the $p_j^y$ integral using Eq. (92) leads to

$$
\begin{align}
\int dp_j^y \begin{pmatrix} \mathcal{J}_{V; j}^t \\ \mathcal{J}_{A; j}^t \\ \mathcal{J}_{H; j}^t \end{pmatrix} &= \frac{|qB|}{4\pi^2} \begin{pmatrix} 1 \\ 2\lambda_j \mathcal{K}_j^+ \mathcal{K}_j^- \\ 2\lambda_j \end{pmatrix}, \nonumber \\
\int dp_j^y \mathcal{J}_{\ell; j}^x &= \int dp_j^y \mathcal{J}_{\ell; j}^y = 0, \nonumber \\
\int dp_j^y \begin{pmatrix} \mathcal{J}_{V; j}^z \\ \mathcal{J}_{A; j}^z \\ \mathcal{J}_{H; j}^z \end{pmatrix} &= \frac{p_j^z |qB|}{4\pi^2 E_j} \begin{pmatrix} 1 \\ 2\lambda_j / \mathcal{K}_j^+ \mathcal{K}_j^- \\ 2\lambda_j \end{pmatrix}, \tag{104}
\end{align}
$$

where the quantities $\mathcal{K}_j^{\pm}$ are given in Eq. (21). It can be seen that $\mathcal{J}_{\ell; j}^{z}$ are odd with respect to $p_j^z \to -p_j^z$ when $n_j > 0$, thus making vanishing contributions to the cor-

responding terms $J_{\ell; n_j}^z$:

$$
\begin{align}
\begin{pmatrix} J_{V; n_j > 0}^t \\ J_{A; n_j > 0}^t \\ J_{H; n_j > 0}^t \end{pmatrix} &= \frac{|qB|}{2\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z}{e^{\beta\mathcal{E}_j} + 1} \begin{pmatrix} \zeta_j \\ 2\lambda_j \mathcal{K}_j^+ \mathcal{K}_j^- \\ 2\lambda_j \zeta_j \end{pmatrix}, \nonumber \\
J_{\ell; n_j > 0}^z &= 0. \tag{105}
\end{align}
$$

F. Energy-momentum tensor

In the case of the energy-momentum tensor $\Theta^{\mu\nu}$, the sesquilinear forms $\mathcal{T}^{\mu\nu}(V_j, V_j)$ corresponding to the antiparticle modes can be related to those pertinent to the particle modes via

$$
\mathcal{T}^{\mu\nu}(V_j, V_j) = -[\mathcal{T}^{\mu\nu}(U_j, U_j)]_{q \to -q}^* \quad (106)
$$

The sesquilinear forms $\mathcal{T}_j^{\mu\nu} \equiv \mathcal{T}^{\mu\nu}(U_j, U_j)$ corresponding to the non-vanishing entries of the canonical energy-momentum tensor (66) can be written in terms of those encountered for the vector current,

$$
\begin{align}
\mathcal{T}_j^{tt} &= E_j \mathcal{J}_{V; j}^t, \quad \mathcal{T}_j^{yy} = p_j^y \mathcal{J}_{V; j}^y, \quad \mathcal{T}_j^{zz} = p_j^z \mathcal{J}_{V; j}^z, \nonumber \\
\mathcal{T}_j^{tz} &= p_j^z \mathcal{J}_{V; j}^t, \quad \mathcal{T}_j^{zt} = E_j \mathcal{J}_{V; j}^z. \tag{107}
\end{align}
$$

The other components vanish, except for $\mathcal{T}_j^{xx}$, which receives non-vanishing contributions only for $n_j > 0$:

$$
\begin{align}
\mathcal{T}_j^{xx} &= -\frac{i}{2} \sqrt{2|qB|} \overline{U}_j \gamma^x \overleftrightarrow{\partial}_{\xi_j} U_j \nonumber \\
&= \frac{|qB|^{3/2} e^{-\xi_j^2/2}}{4\pi^2 E_j (n_j - 1)! \sqrt{\pi}} W_{\xi_j} (H_{n_j-1}, H_{n_j}), \tag{108}
\end{align}
$$

where $W_z(f,g) \equiv f(z)g'(z) - f'(z)g(z)$ is the Wronskian. Due to symmetry constraints, we expect that $\Theta^{xx} = \Theta^{yy}$.

On the LLL, when $n_j = 0$, we have $\mathcal{T}_j^{xx}|_{n_j=0} = \mathcal{T}_j^{yy}|_{n_j=0} = 0$, while the non-vanishing components of $\Theta_0^{\mu\nu}$ can be obtained as

$$
\begin{align}
\Theta_0^{tt} &= \frac{|qB|}{4\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z E_j}{e^{\beta\mathcal{E}_j} + 1}, \nonumber \\
\Theta_0^{tz} &= \Theta_0^{zt} = \frac{qB}{4\pi^2} \sum_{\lambda_j, \zeta_j} 2\zeta_j \lambda_j \int_0^\infty \frac{dp_j^z p_j^z}{e^{\beta\mathcal{E}_j} + 1}, \nonumber \\
\Theta_0^{zz} &= \frac{|qB|}{4\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z p_{z; j}^2}{E_j (e^{\beta\mathcal{E}_j} + 1)}. \tag{109}
\end{align}
$$

In order to compute the contributions due to the higher energy levels ($n_j > 0$), we replace $p_j^y = xqB - \sigma\xi_j \sqrt{|qB|/2}$ and use the relation $H'_{n_j}(\xi_j) = n_j H_{n_j-1}(\xi_j)$,
such that

$$
\begin{align}
\int_{-\infty}^{\infty} dp_j^y \mathcal{T}_j^{tt} \Big|_{n_j>0} &= \frac{E_j |qB|}{4\pi^2}, \nonumber \\
\int_{-\infty}^{\infty} dp_j^y \mathcal{T}_j^{tz} \Big|_{n_j>0} &= \int_{-\infty}^{\infty} dp_j^y \mathcal{T}_j^{zt} \Big|_{n_j>0} = \frac{p_j^z |qB|}{4\pi^2}, \nonumber \\
\int_{-\infty}^{\infty} dp_j^y \mathcal{T}_j^{xx} \Big|_{n_j>0} &= \int_{-\infty}^{\infty} dp_j^y \mathcal{T}_j^{yy} \Big|_{n_j>0} = \frac{n_j (qB)^2}{4\pi^2 E_j}, \nonumber \\
\int_{-\infty}^{\infty} dp_j^y \mathcal{T}_j^{zz} \Big|_{n_j>0} &= \frac{(p_j^z)^2 |qB|}{4\pi^2 E_j}. \tag{110}
\end{align}
$$

It can be seen that $\mathcal{T}_j^{tz}|_{n_j>0} = \mathcal{T}_j^{zt}|_{n_j>0}$ are odd with respect to $p_j^z \rightarrow -p_j^z$, thus $\Theta_{n_j}^{tz} = \Theta_{n_j}^{zt} = 0$ for $n_j > 0$. The nonvanishing components of the contribution of the higher Landau levels to the expectation value of the energy-momentum tensor $\Theta_{n_j}^{\mu\nu}$ are given by

$$
\begin{align}
\Theta_{n_j>0}^{tt} &= \frac{|qB|}{2\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z E_j}{e^{\beta\mathcal{E}_j} + 1}, \\
\Theta_{n_j>0}^{xx} &= \Theta_{n_j>0}^{yy} = \frac{n_j |qB|^2}{2\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z}{E_j(e^{\beta\mathcal{E}_j} + 1)}, \\
\Theta_{n_j>0}^{zz} &= \frac{|qB|}{2\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^\infty \frac{dp_j^z p_{z,j}^2}{E_j(e^{\beta\mathcal{E}_j} + 1)}. \tag{111}
\end{align}
$$

V. ANOMALOUS TRANSPORT

For a classical (non-quantum) fluid in thermodynamic equilibrium, the charge currents $J_\ell^\mu$ and energy-momentum tensor $\Theta^{\mu\nu}$ take the ideal fluid form,

$$
J_{\ell;\text{ideal}}^{\mu} = Q_{\ell} u^{\mu}, \quad \Theta_{\text{ideal}}^{\mu\nu} = \epsilon u^{\mu} u^{\nu} - P \Delta^{\mu\nu}, \quad (112)
$$

where $u^\mu$ is the local fluid four-velocity, given by Eq. (72) for the present case, while $\Delta^{\mu\nu} = g^{\mu\nu} - u^\mu u^\nu$ is the pro- jector on the hypersurface orthogonal to $u^\mu$. The ther- modynamic pressure $P \equiv P(T, \mu_\ell) = -\frac{1}{3}\Delta_{\mu\nu}\Theta^{\mu\nu}_{\text{ideal}}$ is re- lated to the charge densities $Q_\ell = u_\mu J_\ell^\mu$ and the entropy density $s$ via

$$
Q_{\ell} = \frac{\partial P}{\partial \mu_{\ell}}, \qquad s = \frac{\partial P}{\partial T}, \qquad (113)
$$

while the energy density $\epsilon = u_\mu T^{\mu\nu} u_\nu$ is given via the Euler relation by

$$
\epsilon = sT - P + \boldsymbol{\mu} \cdot \mathbf{Q}. \qquad (114)
$$

The purpose of this section is to establish the expres-
sions for the quantities appearing above in the presence
of a constant magnetic field, as well as to highlight devia-
tions from the perfect fluid form, which we will interpret
as a signature of an anomalous transport phenomenon
associated with helical, rather than axial degrees of free-
dom. It is worth mentioning that the term “anomalous

transport” denotes dissipationless – often off-equilibrium – transport effects that appear due to anomalous breaking of a continuous symmetry. The CSE (4) results from the anomalous breaking of the axial symmetry [24]. The presence of an anomaly for helical degrees of freedom has been speculated in Refs. [15, 16] and below we will show that the magnetic field produces the Helical Separation Effect which appears to be very similar to its famous chiral magnetic counterpart, the CME.

A. Hydrodynamic decomposition

In general, the charge currents $J_\ell^\mu$ ($\ell \in \{V, A, H\}$) can be decomposed with respect to the fluid four-velocity $u^\mu$ as follows:

$$
J_{\ell}^{\mu} = Q_{\ell} u^{\mu} + j_{\ell}^{\mu}, \quad Q_{\ell} = J_{\ell}^{t}, \qquad (115)
$$

where $j_\ell^\mu = \Delta_\alpha^\mu J_\ell^\alpha$ represents the charge flow in the fluid rest frame. Since the only vector orthogonal to $u^\mu$ available in this problem is $B^\mu$, we have the following charge flows:

$$
j_{\ell}^{\mu} = \sigma_{\ell}^{B} B^{\mu}, \qquad (116)
$$

with

$$
\sigma_{\ell}^{B} = \frac{1}{B} J_{\ell}^{z}, \qquad (117)
$$

being the magnetic conductivity of the *l*-th charge.

We now consider the fluid described by a symmetric
energy-momentum tensor $\Theta^{\mu\nu}$. Once the velocity $u^\mu$ is
fixed (we take $u^\mu = \delta_\ell^\mu$), $\Theta^{\mu\nu}$ can be decomposed as

$$
\Theta^{\mu\nu} = \epsilon u^{\mu} u^{\nu} - (P + \Pi) \Delta^{\mu\nu} + \pi^{\mu\nu} + j_{\epsilon}^{\mu} u^{\nu} + u^{\mu} j_{\epsilon}^{\nu}. \quad (118)
$$

Besides the energy density $\epsilon$ and thermodynamic pres-
sure $P$, the above decomposition introduces as devi-
ations from the ideal fluid form the dynamic (bulk)
pressure $\Pi = -P - \frac{1}{3}\Delta_{\mu\nu}\Theta^{\mu\nu}$, the shear-stress tensor
$\pi^{\mu\nu} = \Delta_{\alpha\beta}^{\mu\nu}\Theta^{\alpha\beta}$ and the heat flux in the fluid rest frame
$j_{\epsilon}^{\mu} \equiv \Delta_{\alpha}^{\mu} u_{\beta} \Theta^{\alpha\beta}$. In the above, we employed the projector
$\Delta_{\alpha\beta}^{\mu\nu}$ defined as

$$
\Delta_{\alpha\beta}^{\mu\nu} = \frac{1}{2}(\Delta_{\alpha}^{\mu}\Delta_{\beta}^{\nu} + \Delta_{\beta}^{\mu}\Delta_{\alpha}^{\nu}) - \frac{1}{3}\Delta^{\mu\nu}\Delta_{\alpha\beta}. \quad (119)
$$

Specializing the above decomposition to the present case with $u^\mu \partial_\mu = \partial_\ell$, the scalar quantities $\epsilon$, $P$ and $\Pi$ can be obtained as

$$
\epsilon = T^{tt}, \quad P + \Pi = \frac{1}{3}(\Theta^{xx} + \Theta^{yy} + \Theta^{zz}). \quad (120)
$$

In classical fluids, the dynamic pressure $\Pi$ is usually re-
lated to the fluid expansion rate $\theta \equiv \partial_{\mu} u^{\mu}$, which van-
ishes in the present case. However, as we will see later,
the Dirac fluid in a magnetic field exhibits non-vanishing
dynamic pressure.
As in the case of the charge flow $j_{\epsilon}^{\mu}$, the requirement of orthogonality to the fluid velocity $u^{\mu}$ restricts the heat flux $j_{\epsilon}^{\mu}$ to the form

$$j_{\epsilon}^{\mu} = \sigma_{\epsilon}^{B} B^{\mu} \quad \Rightarrow \quad \sigma_{\epsilon}^{B} = \frac{1}{B} \Theta^{tz}, \qquad (121)$$

which we dub the Helical Magnetic Heat Effect [cf. Eq. (5)]. The shear-stress tensor $\pi^{\mu\nu}$ can be written in terms of the available tensors of the problem,

$$\pi^{\mu\nu} = au^{\mu}u^{\nu} + bB^{\mu}B^{\nu} + c(u^{\mu}B^{\nu} + u^{\nu}B^{\mu}) + d\Delta^{\mu\nu}, \quad (122)$$

where a, b, c, and d are certain coefficients. Imposing $\pi^{\mu\nu}u_{\nu} = au^{\mu} + cB^{\mu} = 0$ reveals that $a = c = 0$. Furthermore, the tracelessness condition $\pi^{\mu}_{\mu} = a - bB^{2} + 3d = 0$ gives $b = 3d/B^{2}$, such that

$$\begin{aligned} \pi^{\mu\nu} &= d \left( \Delta^{\mu\nu} + 3 \frac{B^{\mu}B^{\nu}}{B^2} \right) \\ &= \pi_B \times \operatorname{diag} \left( 0, -\frac{1}{2}, -\frac{1}{2}, 1 \right), \end{aligned} \qquad (123)$$

where we defined $\pi_B = 2d$ for notational convenience.

Before ending this section, we remind the reader that in non-equilibrium relativistic hydrodynamics, the mass-energy equivalence gives rise to an ambiguity in defining the energy density, in the sense that the heat flux can be in part or totally absorbed into a new timelike vector that plays the role of the fluid velocity. In the above discussion, we have implicitly employed the so-called beta (or thermometer) velocity frame [32–35], which gives the velocity in terms of the temperature four-vector defining the density operator $\hat{\rho}$. Another frequently used hydrodynamic frame is the so-called Landau (or energy) frame, defined by the eigenvalue equation

$$\Theta^{\mu}_{\nu} u^{\nu}_{L} = \epsilon_{L} u^{\mu}_{L}. \qquad (124)$$

In the present case, it is not difficult to establish that

$$u_{L}^{\mu} \partial_{\mu} = \Gamma_{L}(\partial_{t} + \beta_{L} \partial_{z}), \quad \beta_{L} = \frac{\Theta^{tz}}{\epsilon_{L} + \Theta^{zz}}, \qquad (125)$$

where $\Gamma_L = (1 - \beta_L)^{-1/2}$ is the Lorentz factor of the Landau velocity and the Landau energy density is

$$\epsilon_L = \frac{1}{2} \left[ \Theta^{tt} - \Theta^{zz} + \sqrt{(\Theta^{tt} + \Theta^{zz})^2 - 4\Theta_{tz}^2} \right]. \qquad (126)$$

It is easy to see that if $\Theta^{tz} \neq 0$ then, in general, $\epsilon_L \neq \epsilon$.

## B. Thermometer frame analysis

Let us now write down the full expressions for the charge densities $Q_{\ell} = J_{\ell}^{t}$, energy density $\epsilon = \Theta^{tt}$ and

isotropic pressure computed in the thermometer frame:

$$\begin{aligned} \left( \begin{matrix} Q_V \\ Q_A \\ Q_H \end{matrix} \right) &= \frac{|qB|}{4\pi^2} \sum_{n_j=0}^{\infty} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} \frac{dp_j^z}{e^{\beta\varepsilon_j} + 1} \left( \begin{matrix} \zeta_j \\ 2\lambda_j K_j^{+} K_j^{-} \\ 2\lambda_j \zeta_j \end{matrix} \right), \\ \epsilon &= \frac{|qB|}{4\pi^2} \sum_{n_j=0}^{\infty} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} \frac{dp_j^z}{e^{\beta\varepsilon_j} + 1}, \\ P + \Pi &= \frac{|qB|}{12\pi^2} \sum_{n_j=0}^{\infty} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} \frac{dp_j^z}{E_j(E_j^2 - M^2)} \frac{E_j(E_j^2 - M^2)}{E_j(e^{\beta\varepsilon_j} + 1)}, \end{aligned} \qquad (127)$$

where $g_j = 1$ for $n_j = 0$ and $g_j = 2$ for $n_j > 0$, $K_j^{\pm} = \sqrt{1 \pm M^2/E_j^2}$ [cf. Eq. (21)], and $E_j\sqrt{p_{z; j}^2 + M^2 + 2n_j|qB|}$ as in Eq. (33), while $\mathcal{E}_j = E_j - \zeta_j\mu_V - 2\lambda_j\mu_A - 2\lambda_j\zeta_j\mu_H$ according to Eq. (80). Similarly, $\pi_B = \Theta^{zz} - (P + \Pi)$ can be obtained as

$$\pi_B = \frac{|qB|}{6\pi^2} \sum_{n_j=0}^{\infty} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} \frac{dp_j^z g_j (p_{z; j}^2 - n_j |qB|)}{E_j (e^{\beta \varepsilon_j} + 1)}. \qquad (128)$$

The above expressions do not explicitly distinguish between the LLL and the higher energy-level contributions. On the contrary, the charge conductivities $\sigma_{\epsilon}^{B}$ and the heat conductivity $\sigma_{\epsilon}^{B}$ are fully determined by the contribution from the LLL, receiving vanishing contributions from the $n_j > 0$ energy levels:

$$\begin{aligned} \left( \begin{matrix} \sigma_{V}^{B} \\ \sigma_{A}^{B} \\ \sigma_{H}^{B} \end{matrix} \right) &= \frac{q}{4\pi^2} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} \frac{dp_j^z}{E_j(e^{\beta\varepsilon_j} + 1)} \left( \begin{matrix} 2\lambda_j p_j^z \\ \zeta_j E_j \\ p_j^z \end{matrix} \right), \\ \sigma_B^\epsilon &= \frac{q}{4\pi^2} \sum_{\lambda_j, \zeta_j} 2\zeta_j \lambda_j \int_0^\infty \frac{dp_j^z p_j^z}{e^{\beta\varepsilon_j} + 1}, \end{aligned} \qquad (129)$$

where it is understood that $n_j = 0$ in the above relations.

It is worth stressing that the property that only the single lowest Landau level determines the value of a current, while the degenerate higher level contributions mutually cancel each other, is the characteristic feature of the anomalous currents. At the level of the charged currents, it is an expected feature for the chiral magnetic conductivity $\sigma_{V}^{B}$, responsible for the Chiral Magnetic Effect [36], and the axial magnetic conductivity $\sigma_{A}^{B}$, which gives the Chiral Separation Effect [24]. These conductivities are determined by the axial anomaly (133). What is less expected is that a similar statement is also true for the helical conductivity $\sigma_{H}^{B}$, which, apparently, is not connected to the axial anomaly but may be related to another form of an anomaly associated with the helical degrees of freedom [15, 16].

In the limit of massless fermions, $E_j = |p_j^z|$ on the LLL and the $p_j^z$ integral can be performed analytically, when

$$\begin{aligned} \int_0^\infty \frac{dp_j^z}{e^{\beta\varepsilon_j} + 1} &= \frac{1}{\beta} \ln(1 + e^{\beta q_j \cdot \mu}), \\ \int_0^\infty \frac{dp_j^z p_j^z}{e^{\beta\varepsilon_j} + 1} &= -\frac{1}{\beta^2} \mathrm{Li}_2(-e^{\beta q_j \cdot \mu}), \end{aligned} \qquad (130)$$
with $\mathbf{q}_j \cdot \boldsymbol{\mu} = \zeta_j \mu_V + 2\lambda_j \zeta_j \mu_H + 2\lambda_j \mu_A$ and $\mathrm{Li}_n(z) = \sum_{k=1}^{\infty} z^k/k^n$ being the polylogarithm function [37]. Thus, in the high-temperature limit we obtain the following conductivities for massless fermions:

$$ \sigma_V^B = \frac{q\mu_A}{2\pi^2} + \frac{q\beta\mu_V\mu_H}{4\pi^2} [1 + O(\beta^2)], \quad (131a) $$

$$ \sigma_A^B = \frac{q\mu_V}{2\pi^2} + \frac{q\beta\mu_A\mu_H}{4\pi^2} [1 + O(\beta^2)], \quad (131b) $$

$$ \sigma_H^B = \frac{q}{\pi^2 \beta} \ln 2 + \frac{q\beta\boldsymbol{\mu}^2}{8\pi^2} + O(\beta^3), \quad (131c) $$

where the higher-order terms appearing in $\sigma_V^B$ and $\sigma_A^B$ vanish when $\mu_V\mu_H = 0$ and $\mu_A\mu_H = 0$, respectively. It is remarkable that the helicity current appears even in a neutral plasma.

Let us now discuss the conductivities (131) in more detail. The meaning of these formulas is straightforward: in the background of a magnetic field $\mathbf{B}$, the fermionic system produces electric ($\ell = V$), axial ($\ell = A$), and helical ($\ell = H$) currents [cf. Eq. (116)],

$$ j_{\ell} = \sigma_{\ell}^{B} \mathbf{B}, \quad \ell = V, A, H, \quad (132) $$

with the appropriate conductivities (131).

The first term in Eq. (131a) gives us the well-known expression for the chiral magnetic conductivity for the chiral separation effect (4) in the absence of the helical chemical potential, $\mu_H = 0$. In particular, the coefficient $1/(2\pi^2)$ in this equation is determined by the axial anomaly which dictates, in turn, the non-conservation of the axial charge in the electromagnetic background:

$$ \partial_{\mu} j_{A}^{\mu} = \frac{q^{2}}{16\pi^{2}} \epsilon^{\mu\nu\alpha\beta} F_{\mu\nu} F_{\alpha\beta}, \quad (133) $$

where $\epsilon^{\mu\nu\alpha\beta}$ is the Levi-Civita tensor with $\epsilon^{0123} = +1$. The same statement is true for the first term in the axial magnetic conductivity (131b) which determines the chiral separation effect (4). The next-to-leading terms in both these expressions, Eqs. (131a) and (131b), give the new high-temperature contributions from the helical imbalance to the appropriate conductivities which are suppressed by the first power of temperature, $O(1/T)$. These conductivities vanish in the electrically neutral and chirally neutral plasmas, at $\mu_V = \mu_A = 0$.

Surprisingly, the helical magnetic conductivity (131c) differs substantially from its known CME (131a) and CSE (131b) counterparts. First of all, the coefficient of proportionality does not show an obvious association with the axial anomaly (133) due to the presence of the $\ln 2$ factor (see, however, the discussion on suspected role of a new helical anomaly in Refs. [15, 16]). Secondly, the generation of the helical current is possible in the absence of any (electric, axial or helical) imbalance with all chemical potential set to zero. This property makes it somewhat similar to the mixed axial-gravitational contribution to the heat conductivity which also operates in neutral matter being proportional to the (second power

of) temperature [27]. The helical conductivity is given by the first term in Eq. (131c) proportional to temperature $T = 1/\beta$ and we recover our principal result (1). Finally, similarly to the CME and the CSE effects, the Helical Separation Effect (1) does not change sign under time reversal, $t \to -t$, which is a characteristic feature of anomalous dissipationless transport.

The magnetic heat conductivity evaluates to

$$ \sigma_{\epsilon}^{B} = \frac{q\mu_{H}}{\pi^{2}\beta} \ln 2 + \frac{q\mu_{V}\mu_{A}}{2\pi^{2}} + \frac{q\beta\mu_{H}}{24\pi^{2}} [3\boldsymbol{\mu}^{2} - 2\mu_{H}^{2} + O(\beta^{2})], \quad (134) $$

where the higher-order terms vanish when $\mu_H = 0$. In deriving the above expression, we employed the following expansion of the polylogarithm:

$$ \mathrm{Li}_2(-e^x) = -\frac{\pi^2}{12} - x \ln 2 - \frac{x^2}{4} - \frac{x^3}{24} + \frac{x^5}{960} + O(x^7). \quad (135) $$

As mentioned already, the magnetic conductivities related to anomalous transport receive contributions only from the LLL. The other quantities appearing in Eqs. (127)–(128) receive contributions also from the HLL. However, in the limit of a strong magnetic field, $|qB| \gg T^2$, the HLL contributions are exponentially suppressed. Moving as before to the case of massless fermions, we find

$$ \begin{aligned} Q_{\ell} &= \frac{|qB|\mu_{\ell}}{2\pi^2} + \frac{|qB|\beta}{4\pi^2} \frac{\partial \mu_{\times}^3}{\partial \mu_{\ell}} [1 + O(\beta^2)], \\ \epsilon &= 3(P + \Pi) = \frac{3}{2}\pi_B = \frac{|qB|}{12\beta^2} \\ & \quad + \frac{|qB|\mu^2}{4\pi^2} + \frac{|qB|\beta}{4\pi^2} \mu_{\times}^3 [1 + O(\beta^2)], \end{aligned} \quad (136) $$

where we denoted $\boldsymbol{\mu}^2 = \mu_V^2 + \mu_A^2 + \mu_H^2$ and $\mu_\times^3 = \mu_V\mu_A\mu_H$.

### C. Thermodynamic analysis

Let us consider the grand canonical potential of the system,

$$ \Phi = - \frac{|qB|V}{4\pi^2\beta} \sum_{n_j=0}^{\infty} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} dp_j^z g_j \ln(1 + e^{-\beta\mathcal{E}_j}), \quad (137) $$

where $V$ is the (infinite) system size. The grand potential $\Phi$ is related to the partition function $\mathcal{Z} = \mathrm{tr}(\hat{\rho})$ via

$$ \mathcal{Z} = e^{-\beta \Phi}. \quad (138) $$

Given the thermodynamic relation

$$ d\Phi = \beta^{-2} S d\beta - P dV - \mathcal{Q} \cdot d\boldsymbol{\mu}, \quad (139) $$
we can identify the total entropy $S$, thermodynamic pressure $P$ and total charge $Q$ as

$$S = \beta^2 \frac{\partial \Phi}{\partial \beta}, \quad P = -\frac{\partial \Phi}{\partial V}, \quad Q_{\ell} = -\frac{\partial \Phi}{\partial \mu_{\ell}}. \qquad (140)$$

Since the grand potential is linearly proportional to the system volume, $\Phi = \phi V$ with $\phi$ being a volume-independent quantity, it is clear that

$$P = -\phi = \Theta^{zz}. \qquad (141)$$

It is easy to check that $Q_{V/H} \equiv Q_{V/H}/V = \partial\Theta^{zz}/\partial\mu_{V/H}$. A similar relation holds for $Q_A$ only in the case of massless fermions ($M = 0$ and $K_j^+ = K_j^- = 1$), confirming that the chiral charge is conserved only for massless fermions.

The total energy $\mathcal{E} = \langle \hat{H} \rangle = \mathcal{Z}^{-1}\text{tr}(\hat{\rho}\hat{H})$ can be obtained from $\mathcal{Z}$ via

$$\mathcal{E} = \mu \cdot \mathcal{Q} - \frac{\partial \ln \mathcal{Z}}{\partial \beta}. \qquad (142)$$

Since $\ln \mathcal{Z} = -\beta\Phi$, it is easy to uncover the Euler relation

$$s = \beta(\epsilon + P - \mu \cdot Q), \qquad (143)$$

where $s = S/V$, $\epsilon = \mathcal{E}/V$ and $Q = \mathcal{Q}/V$.

The above discussion shows that the entropy density can be obtained as

$$s = \beta \left( \Theta^{tt} + \Theta^{zz} - \mu \cdot \frac{\partial \Theta^{zz}}{\partial \mu} \right) = \frac{\partial \Theta^{zz}}{\partial T}. \qquad (144)$$

Having identified $P = \Theta^{zz}$, we are in a position to evaluate the dynamical pressure $\Pi$ using Eq. (120):

$$\Pi = \frac{1}{3}(\Theta^{xx} + \Theta^{yy} - 2\Theta^{zz}) = -\pi_B \\ = \frac{|qB|}{6\pi^2} \sum_{n_j=0}^{\infty} \sum_{\lambda_j, \zeta_j} \int_0^{\infty} \frac{dp_j^z g_j(n_j|qB| - p_{z;j}^2)}{E_j(e^{\beta\mathcal{E}_j} + 1)}. \qquad (145)$$

On the LLL and for massless fermions, the relation $\Pi = -\pi_B$ can be used in conjunction with Eq. (136) to obtain:

$$\epsilon = P, \quad \pi_B = -\Pi = \frac{2}{3}P. \qquad (146)$$

## VI. SEARCH FOR HELICAL EXCITATIONS IN NEUTRAL PLASMA

Let us consider a globally-neutral plasma, characterized by vanishing chemical potentials $\bar{\mu}_{\ell} = 0$ ($\ell = V, A, H$) and temperature $\bar{T}$, where the bar over the symbol indicates that we consider a mean global coordinate-independent quantity. We work in the hydrodynamic approximation, which assumes small departures from local thermodynamic equilibrium. We also assume the limit of a strong magnetic field, $|qB| \gg T^2$.

In a strong magnetic field, the excess in vector or axial densities generates a coherent propagation of these charges along the direction of the magnetic field in the form of a linear hydrodynamic excitation, known as the Chiral Magnetic Wave [38]. In this wave, an excess in vector (axial) charge generates an axial (vector) current along the magnetic field, which leads to a build-up of an axial (vector) charge and a cyclic repetition of the process.

We look for a similar excitation in terms of helical degrees of freedom. To this end we notice that the helical current $j_H = \sigma_H^B B$ and the heat current $j_{\epsilon} = \sigma_{\epsilon}^B B$, can be expressed, respectively, via Eq. (131c) for the helical conductivity $\sigma_H^B$ and Eq. (134) for the heat conductivity $\sigma_{\epsilon}^B$. These relations show that helical fluctuations are coupled to the fluctuations in temperature similarly to the coupling of vector and axial charges (currents) in the Chiral Magnetic Wave. Therefore, one could expect that the helical degrees of freedom may propagate in terms of a hypothetical “helical heat wave,” which would corresponds to a coherent hydrodynamic excitation that combines energy density and helical fluctuations.

To assess the existence of the helical heat wave as a linear hydrodynamic excitation, we consider the equations of motion derived from the charge (for all $\ell = V, A, H$), energy, and momentum conservation, namely:

$$\dot{Q}_{\ell} + Q_{\ell}\theta + \partial_{\mu}j_{\ell}^{\mu} = 0, \qquad (147a)$$

$$\dot{\epsilon} + (\epsilon + P + \Pi)\theta - \pi^{\mu\nu}\sigma_{\mu\nu} + \partial_{\mu}j_{\epsilon}^{\mu} - j_{\epsilon}^{\nu}\dot{u}_{\nu} = 0, \qquad (147b)$$

$$(\epsilon + P + \Pi)\dot{u}^{\mu} - \nabla^{\mu}(P + \Pi) + \Delta_{\lambda}^{\mu}\partial_{\nu}\pi^{\nu\lambda} + j_{\epsilon}^{\lambda}\nabla_{\lambda}u^{\mu} + Dj_{\epsilon}^{(\mu)} + j_{\epsilon}^{\mu}\theta = 0, \qquad (147c)$$

where $D$ and the overhead dot denote the comoving derivative, $Da \equiv \dot{a} \equiv u^{\mu}\partial_{\mu}a$, $\sigma_{\mu\nu} = \Delta_{\mu\nu}^{\alpha\beta}\partial_{\alpha}u_{\beta}$ is the shear tensor with $\Delta_{\mu\nu}^{\alpha\beta}$ introduced in Eq. (118) and $\theta = \partial_{\mu}u^{\mu}$ is the expansion scalar. In the following, we focus on the dissipationless transport phenomena supported by the above set of equations and ignore dissipative corrections to the non-ideal terms $\Pi$, $\pi^{\mu\nu}$, $j_{\epsilon}^{\mu}$, and $j_{\ell}^{\mu}$. Also, we consider that the electric field is negligible. In reality, one should allow for fluctuations in $B^{\mu}$ and in the electric field $E^{\mu}$, coming from the Maxwell equations. The backreaction from fluctuating electromagnetic field usually leads to a damping effect [38] which we neglect below.

We now consider infinitesimal fluctuations in a quiescent, neutral fluid of massless particles, at background temperature $\bar{T}$, in a strong magnetic field $|qB| \gg \bar{T}^2$. We consider the generic split of a hydrodynamic quantity $a$ into its global average $\bar{a}$ and fluctuation $\delta a$ as follows,

$$a = \bar{a} + \delta a, \quad \delta a = \int d^3k \tilde{\delta}a_k(\omega) e^{-i\omega t + i\mathbf{k}\cdot\mathbf{x}}, \qquad (148)$$

where $\tilde{\delta}a_k(\omega)$ represents the amplitude of its Fourier mode of wave vector $\mathbf{k}$. For brevity, we henceforth suppress the explicit $\mathbf{k}$ and $\omega$ dependencies of fluctuations.
We take as independent parameters the fluid temperature $T$, chemical potentials $\mu_\ell$, and macroscopic velocity $u^\mu$, with $\bar{\mu}_\ell = 0$ and $\bar{u}^\mu = \delta_0^\mu$. The velocity fluctuations lead to a modification $u^\mu = \bar{u}^\mu + \delta u^\mu$ that affects also the covariant magnetic field $B^\mu = \bar{B}^\mu + \delta B^\mu$, such that:

$$ \overline{B}^\mu = B\delta_z^\mu, \quad \delta B^\mu = B\delta u^z \delta_z^\mu. \quad (149) $$

Also, the comoving derivative $\dot{a}$ becomes

$$ \dot{a} \approx \partial_t \delta a \rightarrow -i\omega \delta \tilde{a}, \quad (150) $$

where the arrow indicates that the Fourier transform was taken as in Eq. (148). More generally, the derivative leads simply to the replacement

$$ \partial_{\mu} a^{\nu} \rightarrow -ik_{\mu} \delta \tilde{a}^{\nu}, \quad (151) $$

and in particular, $\theta \rightarrow -ik_\mu \delta \tilde{u}^\mu = i\mathbf{k} \cdot \delta \tilde{u}$.

In the case of the charge conservation equation, Eq. (147a), we have

$$ \begin{aligned} \bar{Q}_\ell &= 0, & \delta Q_\ell &= \frac{|qB|}{2\pi^2} \delta \mu_\ell, \\ \bar{\sigma}_{V/A}^B &= 0, & \delta \sigma_{V/A}^B &= \frac{q}{2\pi^2} \delta \mu_{A/V}, \\ \bar{\sigma}_H^B &= \frac{q\bar{T} \ln 2}{\pi^2}, & \delta \sigma_H^B &= \frac{q \ln 2}{\pi^2} \delta T. \end{aligned} \quad (152) $$

Taking into account that $\bar{\theta} = 0$, the second term $Q_\ell \theta$ in Eq. (147a) becomes negligible. The third term evaluates to

$$ \begin{aligned} \partial_{\mu} j_{\ell}^{\mu} ={}& B \partial_{z} \delta \sigma_{\ell}^{B} + B \bar{\sigma}_{\ell}^{B} \partial_{t} \delta u^{z} \\ & \rightarrow ik^{z} B \delta \tilde{\sigma}_{\ell}^{B} - i \omega B \bar{\sigma}_{\ell}^{B} \delta \tilde{u}^{z}. \end{aligned} \quad (153) $$

Since $\bar{\sigma}_{V/A}^B = 0$, the equations giving the conservation of $Q_V$ and $Q_A$ reduce to the closed set

$$ \left( \begin{array}{cc} \omega & -k^z \sigma \\ -k^z \sigma & \omega \end{array} \right) \left( \begin{array}{c} \tilde{\delta \mu}_V \\ \tilde{\delta \mu}_A \end{array} \right) = 0. \quad (154) $$

The above equations are solved when either $\omega = \pm k^z$ or when both $\tilde{\delta \mu}_V$ and $\tilde{\delta \mu}_A$ cancel. Hence, we have uncovered the chiral magnetic wave which comprises coherent oscillations of vector and axial charges propagating, in strong magnetic field background, with the speed of light [38].

Since $\bar{\sigma}_H^B \neq 0$, the conservation equation for the helicity charge $Q_H$ leads to a coupling between the helicity chemical potential $\tilde{\delta \mu}_H$ and the fluctuations in the hydrodynamic sector,

$$ \omega \delta \tilde{u}^z - k^z \frac{\delta \tilde{T}}{\bar{T}} + \frac{\sigma \omega}{2 \ln 2} \frac{\delta \tilde{\mu}_H}{\bar{T}} = 0. \quad (155) $$

Moving now to the $T^{\mu\nu}$ sector, Eq. (136) shows that $\epsilon = P$ and $\Pi = -\pi_B = -\frac{2}{3}P$. Considering that these equalities hold also in the perturbed system, we have

$$ \begin{aligned} \bar{\epsilon} &= \frac{|qB|\bar{T}^2}{12}, & \delta\epsilon &= \frac{|qB|\bar{T}}{6}\delta T, \\ \bar{\sigma}_{\epsilon}^{B} &= 0, & \delta\sigma_{\epsilon}^{B} &= \frac{q\bar{T}\ln 2}{\pi^2}\delta\mu_H. \end{aligned} \quad (156) $$

In the case of the shear-stress tensor $\pi^{\mu\nu}$, we ignore dissipative corrections and employ the form in Eq. (123). Writing $\pi^{\mu\nu} = \bar{\pi}^{\mu\nu} + \delta\pi^{\mu\nu}$, we have $\bar{\pi}^{\mu\nu} = \frac{1}{2}\bar{\pi}_B \times \text{diag}(0, -1, -1, 2)$ with $\bar{\pi}_B = \frac{2}{3}\bar{P}$ and

$$ \begin{aligned} \delta\pi^{\mu\nu} = \frac{1}{2}\delta\pi_B \operatorname{diag}(0, -1, -1, 2) &+ \frac{\bar{\pi}_B}{2} \begin{pmatrix} 0 & -\delta u^x & -\delta u^y & 2\delta u^z \\ -\delta u^x & 0 & 0 & 0 \\ -\delta u^y & 0 & 0 & 0 \\ 2\delta u^z & 0 & 0 & 0 \end{pmatrix}. \end{aligned} \quad (157) $$

The shear tensor $\sigma^{\mu\nu}$ has vanishing average, $\bar{\sigma}^{\mu\nu} = 0$, while $\delta\sigma^{\mu\nu} = \nabla_\mu(\delta u_\nu) - \frac{1}{3}\Delta_{\mu\nu}\theta$ has vanishing components when either $\mu = t$ or $\nu = t$, i.e. $\delta\sigma^{00} = \delta\sigma^{0i} = \delta\sigma^{i0} = 0$. On the spatial part, we have $\delta\sigma_{ij} = \partial_{(i}\delta u_{j)} - \frac{1}{3}g_{ij}\theta$, or in Fourier space,

$$ \tilde{\delta\sigma}^{ij} = -\frac{i}{2}(k^i \tilde{\delta u}^j + k^j \tilde{\delta u}^i) - \frac{i}{3}g^{ij} \mathbf{k} \cdot \tilde{\delta u}. \quad (158) $$

In the case of the energy conservation equation, we employ $\pi^{\mu\nu}\sigma_{\mu\nu} = -\frac{1}{2}\pi_B(\delta\sigma^{xx} + \delta\sigma^{yy} - 2\delta\sigma^{zz}) \rightarrow \frac{i}{3}\bar{P}(\mathbf{k} \cdot \tilde{\delta}\mathbf{u} - 3k^z\tilde{\delta}u^z)$, $j_\epsilon^\nu \dot{u}_\nu \simeq 0$, and

$$ \partial_{\mu} j_{\epsilon}^{\mu} = B \partial_{z} \delta \sigma_{\epsilon}^{B} \rightarrow \frac{ik \bar{P}}{\bar{T}} \frac{12 \ln 2}{\pi^2} \delta \bar{\mu}_{H} \quad (159) $$

to arrive at

$$ 2\omega \frac{\tilde{\delta T}}{\bar{T}} - \mathbf{k} \cdot \tilde{\delta \mathbf{u}} - k^z \tilde{\delta u}^z - \frac{12\sigma \ln 2}{\pi^2} k^z \frac{\tilde{\delta \mu}_H}{\bar{T}} = 0. \quad (160) $$

Now we move on to the momentum conservation equation, i.e. Eq. (147c). The terms $j_\epsilon^\lambda \nabla_\lambda u^\mu$ and $j_\epsilon^\mu \theta$ are of higher order and can thus be ignored. Furthermore, the projector $\Delta_\nu^\mu$ together with $\dot{u}^0 \simeq 0$ ensures that the $\mu = 0$ relation reduces to a tautological $0 = 0$. For $\mu = i \in \{x, y\}$ and $\mu = z$, we have

$$ \begin{aligned} \Delta_{\lambda}^{i} \partial_{\nu} \pi^{\nu \lambda} ={}& \frac{1}{2} \partial^{i} \delta \pi_{B} + \frac{\bar{\pi}_{B}}{2} \partial_{t} \delta u^{i} \rightarrow -\frac{i}{3} \left( 2k^{i} \frac{\delta \tilde{T}}{\bar{T}} - \omega \delta \tilde{u}^{i} \right), \\ \Delta_{\lambda}^{z} \partial_{\nu} \pi^{\nu \lambda} ={}& \partial_{z} \delta \pi_{B} + \bar{\pi}_{B} \partial_{t} \delta u^{z} \rightarrow \frac{2i}{3} \left( 2k^{z} \frac{\delta \tilde{T}}{\bar{T}} - \omega \delta \tilde{u}^{z} \right). \end{aligned} \quad (161) $$

We thus arrive at:

$$ \omega \delta \tilde{u}^i = 0, \quad \omega \delta \tilde{u}^z - k^z \frac{\delta \tilde{T}}{\bar{T}} + \omega \frac{6\sigma \ln 2}{\pi^2} \delta \tilde{\mu}_H = 0. \quad (162) $$
The first relation above shows that the plasma does not support oscillations in the plane orthogonal to the magnetic field. The second relation in Eq. (162) can be used in conjunction with Eqs. (155) and (160) to derive

$$ \begin{pmatrix} \omega & -k^z & \frac{\sigma\omega}{2\ln 2} \\ -k^z & \omega & -\frac{6\sigma\ln 2}{\pi^2}k^z \\ \omega & -k^z & \frac{6\sigma\ln 2}{\pi^2}\omega \end{pmatrix} \begin{pmatrix} \tilde{\delta u}^z \\ \tilde{\delta T} \\ \tilde{\delta \mu}_H \end{pmatrix} = 0. \quad (163) $$

In order to have non-trivial solutions of the above equation, the determinant of the matrix on the left-hand side must vanish. This leads to the equation

$$ \sigma\omega \left( \frac{1}{2\ln 2} - \frac{6\ln 2}{\pi^2} \right) (\omega^2 - k_z^2) = 0, \quad (164) $$

where the solution $\omega = 0$ corresponds to the mode that supports helicity chemical potential fluctuations. The gapless modes $\omega = \pm k^z$ have identically $\tilde{\delta\mu}_H = 0$ and correspond to the standard sound modes present in the ideal fluid, which in the present case propagate at the speed of light.

## VII. CONCLUSION

In our paper, we have found the Helical Separation Effect (HSE) which implies that charged massless fermions develop a dissipationless flow of helicity along the background magnetic field (1). This dissipationless effect is similar to the Chiral Separation Effect (CSE) [24], which generates the flow of the chiral charge along magnetic field lines in dense fermionic matter. Contrary to the CSE, the HSE produces the helical current in a neutral plasma in which all chemical potentials vanish. Moreover, the expression for the HSE current (1) contains a “ln 2” factor thus making it difficult to associate this effect to a known anomaly (see, however, the discussions in Ref. [15, 16] on possible existence of helical anomalies). The associated helical conductivity is determined

by temperature which makes it distantly similar to the axial-gravitational effects [27]. The presence of matter plays an auxiliary role enhancing the helical current (2), (3) generated by the HSE. The same remarks are also valid for the Helical Magnetic Heat Effect (5) which generates, in the presence of non-vanishing helical charge density, the fermionic heat flux along the magnetic field.

Since the helical current generated by the HSE is linearly proportional to the background temperature and magnetic field, the helical current (1) and the heat current (5) become coupled together via the HSE in a strikingly similar way with the coupling mediated by the CME and CSE of the vector and chirality currents that result in the appearance of the chiral magnetic wave [38]. However, a detailed analysis shows that the fluctuations of the helical charge in a neutral plasma do not lead to the appearance of a new hydrodynamic mode. The helical degree of freedom may, however, lead to the appearance of additional hydrodynamic waves in dense plasma, which are left beyond the scope of this paper.

The helicity plays a complementary role to chirality, as the total helicity of a particle-antiparticle ensemble is given by the sum of axial charges carried by particles minus the sum of the axial charges possessed by antiparticles. As the helicity is different from chirality, its role in phenomenological applications in Quark-Gluon plasma may lead to new observational effects which remain to be explored.

## ACKNOWLEDGMENTS

The authors are grateful to Yago Ferreiros for collaborating on an earlier related project. VEA gratefully acknowledges the support through a grant of the Ministry of Research, Innovation and Digitization, CNCS - UE-FISCDI, project number PN-III-P1-1.1-TE-2021-1707, within PNCDI III. The authors gratefully acknowledge support by the European Union - NextGenerationEU through the grant No. 760079/23.05.2023, funded by the Romanian ministry of research, innovation and digitalization through Romania's National Recovery and Resilience Plan, call no. PNRR-III-C9-2022-18.

[1] Grigory E Volovik, The universe in a helium droplet, Vol. 117 (OUP Oxford, 2003).

[2] NP Armitage, EJ Mele, and Ashvin Vishwanath, “Weyl and dirac semimetals in three-dimensional solids,” Reviews of Modern Physics 90, 015001 (2018).

[3] MM Vazifeh and M Franz, “Electromagnetic response of weyl semimetals,” Physical review letters 111, 027201 (2013).

[4] M. N. Chernodub, Alberto Cortijo, and Maria A. H. Vozmediano, “Generation of a Nernst Current from the Conformal Anomaly in Dirac and Weyl Semimetals,” Phys. Rev. Lett. 120, 206601 (2018), arXiv:1712.05386 [cond-mat.str-el].

[5] Johannes Gooth, Anna C Niemann, Tobias Meng, Adolfo G Grushin, Karl Landsteiner, Bernd Gotsmann, Fabian Menges, Marcus Schmidt, Chandra Shekhar, Vicky Süß, et al., “Experimental signatures of the mixed axial-gravitational anomaly in the weyl semimetal nbp,” Nature 547, 324–327 (2017).

[6] Maxim N. Chernodub, Yago Ferreiros, Adolfo G. Grushin, Karl Landsteiner, and María A. H. Vozmediano, “Thermal transport, geometry, and anomalies,” Phys. Rept. 977, 1–58 (2022), arXiv:2110.05471 [cond-mat.mes-hall].
[7] Alberto Cortijo, Yago Ferreirós, Karl Landsteiner, and María AH Vozmediano, “Elastic gauge fields in weyl semimetals,” Physical review letters **115**, 177202 (2015).

[8] K. Fujikawa and H. Suzuki, *Path integrals and quantum anomalies* (2004).

[9] Dmitri E. Kharzeev, “The Chiral Magnetic Effect and Anomaly-Induced Transport,” *Prog. Part. Nucl. Phys.* **75**, 133–151 (2014), arXiv:1312.3348 [hep-ph].

[10] Palash B. Pal, “Dirac, Majorana and Weyl fermions,” *Am. J. Phys.* **79**, 485–498 (2011), arXiv:1006.1718 [hep-ph].

[11] M. N. Chernodub and Victor E. Ambrus, “Phase diagram of helically imbalanced QCD matter,” *Phys. Rev. D* **103**, 094015 (2021), arXiv:2005.03575 [hep-th].

[12] Shen-Song Wan and Marco Ruggieri, “Chiral phase transitions of helical matter,” *Phys. Rev. D* **103**, 094030 (2021), arXiv:2012.12577 [hep-ph].

[13] M. N. Chernodub and A. S. Nedelin, “Phase diagram of chirally imbalanced QCD matter,” *Phys. Rev. D* **83**, 105008 (2011), arXiv:1102.0188 [hep-ph].

[14] Marco Ruggieri, “The Critical End Point of Quantum Chromodynamics Detected by Chirally Imbalanced Quark Matter,” *Phys. Rev. D* **84**, 014011 (2011), arXiv:1103.6186 [hep-ph].

[15] Victor E. Ambrus, “Helical massive fermions under rotation,” JHEP **08**, 016 (2020), arXiv:1912.09977 [nucl-th].

[16] Victor E. Ambrus and M. N. Chernodub, “Vortical effects in Dirac fluids with vector, chiral and helical charges,” Eur. Phys. J. C **83**, 111 (2023), arXiv:1912.11034 [hep-th].

[17] Victor E. Ambrus and M. N. Chernodub, “Hyperon-antihyperon polarization asymmetry in relativistic heavy-ion collisions as an interplay between chiral and helical vortical effects,” Eur. Phys. J. C **82**, 61 (2022), arXiv:2010.05831 [hep-ph].

[18] N. Byers and C. N. Yang, “Theoretical Considerations Concerning Quantized Magnetic Flux in Superconducting Cylinders,” Phys. Rev. Lett. **7**, 46–49 (1961).

[19] Bascom S. Deaver and William M. Fairbank, “Experimental Evidence for Quantized Flux in Superconducting Cylinders,” Phys. Rev. Lett. **7**, 43–46 (1961).

[20] Ion I. Cotăescu, Doru-Marcel Băltățeanu, and Ion Cotăescu, “Relativistic persistent currents in ideal Aharonov-Bohm rings,” *International Journal of Modern Physics B* **30**, 1550245 (2016)

[21] Daniel M. Lipkin, “Existence of a new conservation law in electromagnetic theory,” *Journal of Mathematical Physics* **5**, 696–700 (1964).

[22] M. N. Chernodub, Alberto Cortijo, and Karl Landsteiner, “Zilch vortical effect,” *Phys. Rev. D* **98**, 065016 (2018), arXiv:1807.10705 [hep-th].

[23] Matteo Buzzegoli, “Spin polarization induced by magnetic field and the relativistic Bar-

nett effect,” Nucl. Phys. A **1036**, 122674 (2023), arXiv:2211.04549 [nucl-th].

[24] D. T. Son and Ariel R. Zhitnitsky, “Quantum anomalies in dense matter,” Phys. Rev. D **70**, 074018 (2004), arXiv:hep-ph/0405216.

[25] Dmitri E Kharzeev, Yuta Kikuchi, René Meyer, and Yuya Tanizaki, “Giant photocurrent in asymmetric weyl semimetals from the helical magnetic effect,” Physical Review B **98**, 014305 (2018).

[26] Naoki Yamamoto and Di-Lun Yang, “Helical magnetic effect and the chiral anomaly,” Phys. Rev. D **103**, 125003 (2021), arXiv:2103.13208 [hep-th].

[27] Karl Landsteiner, Eugenio Megias, and Francisco Pena-Benitez, “Gravitational Anomaly and Transport,” Phys. Rev. Lett. **107**, 021601 (2011), arXiv:1103.5006 [hep-ph].

[28] Nikita Astrakhantsev, V. V. Braguta, Massimo D'Elia, A. Yu. Kotov, A. A. Nikolaev, and Francesco Sanfilippo, “Lattice study of the electromagnetic conductivity of the quark-gluon plasma in an external magnetic field,” Phys. Rev. D **102**, 054516 (2020), arXiv:1910.08516 [hep-lat].

[29] V. Canuto and H. Y. Chiu, “Quantum theory of an electron gas in intense magnetic fields,” Phys. Rev. **173**, 1210–1219 (1968).

[30] V. Canuto and H. Y. Chiu, “Thermodynamic properties of a magnetized Fermi gas,” Phys. Rev. **173**, 1220–1228 (1968).

[31] V. Canuto and H. Y. Chiu, “Magnetic moment of a magnetized Fermi gas,” Phys. Rev. **173**, 1229–1235 (1968).

[32] P. Ván and T. S. Biró, “First order and stable relativistic dissipative hydrodynamics,” Phys. Lett. B **709**, 106–110 (2012).

[33] P. Ván and T. S. Biró, “Dissipation flow-frames: particle, energy, thermometer,” in *Proceedings of the 12th Joint European Thermodynamics Conference*, edited by M. Pilotelli and G. P. Beretta (2013) pp. 546–551.

[34] K. Landsteiner, E. Megías, and F. Pena-Benitez, “Anomalous transport from kubo formulae,” in *Strongly interacting matter in magnetic fields* (Lecture notes in physics **871**), edited by D. Kharzeev, K. Landsteiner, A. Schmitt, and H.-U. Yee (Springer-Verlag, Berlin Heidelberg, 2013) pp. 433–468.

[35] F. Becattini, L. Bucciantini, E. Grossi, and L. Tinti, “Local thermodynamical equilibrium and the β-frame for a quantum relativistic fluid,” Eur. Phys. J C **75**, 191 (2015).

[36] Kenji Fukushima, Dmitri E. Kharzeev, and Harmen J. Warringa, “The Chiral Magnetic Effect,” Phys. Rev. D **78**, 074033 (2008), arXiv:0808.3382 [hep-ph].

[37] F. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, *NIST handbook of mathematical functions hardback and CD-ROM* (Cambridge university press, 2010).

[38] Dmitri E. Kharzeev and Ho-Ung Yee, “Chiral Magnetic Wave,” Phys. Rev. D **83**, 085007 (2011), arXiv:1012.6026 [hep-th].