Non-Hermitian Chiral Magnetic Effect in Equilibrium

We analyze the Chiral Magnetic Effect for non-Hermitian fermionic systems using the biorthogonal formulation of quantum mechanics. In contrast to the Hermitian chiral counterparts, we show that the Chiral Magnetic Effect may take place in thermal equilibrium of an open non-Hermitian system with, generally, massive fermions. The key observation is that for non-Hermitian charged systems, there is no strict charge conservation as understood in the Hermitian case, so the Bloch theorem preventing currents in the thermodynamic limit in equilibrium does not apply.


Introduction
The Chiral Magnetic Effect (CME) is the generation of an electric current J in presence of an external magnetic field B [1]: The current (1) appears naturally in a particular set of physical systems characterized by a broken invariance under the spatial reflection P. The broken P-invariance may be realized, for example, in condensed matter systems with massless fermionic quasiparticles of a Weyl or Dirac type.Such fermions are characterized by different (left and right) chiralities which are often said to be ascribed to different "two Weyl nodes".If the number of fermions with different chiralities is not equal to each other, then the system is P-broken.
The chiral imbalance is convenient to characterize by a difference denoted by µ 5 = µ L − µ R between the chemical potentials in the right (µ R ) and left (µ L ) Weyl nodes.The difference in the chemical potentials determines the magnitude of the CME current in Eq. ( 1).
The CME is a relevant transport phenomenon that has its roots in the physics of quantum anomalies.The theory is said to be anomalous if there exists a quantity which is conserved at the classical level and which fails to do so when going to the quantum realm.In particular, the CME stems from the axial anomaly which leads to non-conservation of the chiral current in Weyl systems described by the Weyl Hamiltonian of a massless particle with the wavefunction ψ s : Here the parameter s = ±1 denotes the chirality of the particle propagating with the velocity v F and momentum k and σ is the vector of the Pauli matrices.In relativistic systems, Lorentz invariance forces v F = c, while in condensed matter systems, the velocity v F is not constrained to any particular value.It is often said that the parameter s = ±1 labels two distinct "Weyl points".The Weyl systems (2) possess two quantities that are conserved at the level of classical equations of motion.These are the electric charge Q and the axial charge Q 5 that are described, respectively, by electric four-current J µ and the chiral current 1 j µ 5 : J µ ≡ (Q, J) = Ψγ µ Ψ, j µ 5 ≡ (q 5 , j 5 ) = Ψγ 5 γ µ Ψ.
Mathematically, the conservation implies that the four-divergence of the both currents is identically zero, ∂ µ J µ = ∂j µ 5 = 0, provided the wave function Ψ = (ψ R , ψ L ) T satisfies the classical equation of motion HΨ = 0 where the Hamiltonian H = diag (H R , H L ) incorporates both chiral modes (2).We use the conventional nomenclature of γ matrices in the chiral basis: In presence of conventional electromagnetic fields, the quantum fluctuations lead to nonconservation of the chiral current.Technically, the loss of conservation appears as a result of the so-called triangle diagrams of virtual fermions that lead to [2,3]: The triangle diagrams that give rise to the non-conserved current (5) are also responsible for the generation of the current in Eq. (1).The form of the anomaly can be partially fixed by some algebraic constraints on an effective action of the theory that leads the right hand side of Eq.( 5).These constraints are imposed by the algebraic structure of the symmetry that is anomalously broken [4].Technically, the relation of the chiral (triangular) anomaly to the generation CME current (1) requires a rigorous derivation which takes into account the Wess-Zumino consistency conditions [5].In the presence of the chiral gauge fields A 5  µ that are coupled to the chiral current j µ 5 , the anomalous effects become more subtle and the currents (3) have to be modified consistently.In our work A 5  µ ≡ 0 so that we will use the straightforward definition (3) for vector and axial currents.
In general, two ways have been proposed to create an environment that is able to generate the anomalous current (1).The first approach is to drive the system out of equilibrium in order to reach a stationary regime where µ 5 = 0.This regime may be achieved by applying simultaneously an electric field E parallel to B so that the chiral anomaly creates a charge imbalance via the chiral anomaly (5) and generates a nonzero chiral chemical potential µ 5 .Then the system generates a non-equilibrium electric current via the CME mechanism (1).Notice that the chiral imbalance µ 5 = 0 does not exist in an equilibrium regime as the population of the left-handed particles and the right-hand particles mix with each other due to interactions and then relaxes towards µ 5 = 0 equilibrium.Therefore, the anomalous electric current (1) is zero in thermal equilibrium.
The second approach could consist in moving the position of the Weyl nodes in energies without carrying the system out of equilibrium: We reserve the notation j for another current to be defined below (12).
< l a t e x i t s h a 1 _ b a s e 6 4 = " y w + 7 4 h G B 1 M M P g T Z D 6 9 8 k e a 9 F Y 1 j 1 A n k q + / A f y M 8 l g = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " y w + 7 4 h G B 1 M M P g T Z D 6 9 8 k e a 9 F Y 1 j 1 A n k q + / A f y M 8 l g = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " y w + 7 4 h G B 1 M M P g T Z D 6 9 U j e J 5 + S z 8 m X P 9 J k s K p 5 i H q R f P 0 N K J r p q Q = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " K E o 8 4 r h J O w 7 U j e J 5 + S z 8 m X P 9 J k s K p 5 i H q R f P 0 N K J r p q Q = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " K E o 8 4 r h J O w 7 where the energies L,R are the position of left and right Weyl points in energy.The equilibrium current ( 6) is, however, vanishing in Hermitian systems.Physically, the current is zero because the difference in energies of the right-and left-handed chiral fermions does not create the true chiral imbalance.Mathematically, the current ( 6) does not exist because the difference in energies ( R − L ) is sensitive to the chirality of the fermion and, therefore, is nothing but the zeroth component of a chiral gauge field, A 5 0 .Once the chiral gauge field appears, the definition of the physical electric current starts to differ from the naive covariant version (3) by an addition of an extra term coming from the so-called Bardeen polynomials.This term cancels out this energy difference precisely, and physical (so called, "consistent") version of the current, vanishes in thermal equilibrium (6).For further details and technicalities we refer an interested reader to Ref. [5].
Despite in thermodynamic equilibrium the axial chemical potential is zero, µ 5 = 0, the vector chemical µ = µ R + µ L for a generic fermionic system may still be nonzero.In the presence of the background magnetic field B, the system generates -via the same chiral anomaly -the chiral current: which is the direct analogue of the CME (1) but now the chiral sector.Equation ( 7) describes the Chiral Separation Effect (CSE) which will play a role in our derivations below along with the CME.While the non-equilibrium situation has been explored extensively in the literature leading, for instance, to the celebrated negative quadratic magnetoresistivity in Weyl metals, the equilibrium scenario appears to be not possible, and to date there is consensus that the CME is not possible in equilibrium [6][7][8].
The statement of absence of CME in equilibrium can be seen as an extension of a no-go theorem given by Bloch, concerning the existence of equilibrium currents in solids in the thermodynamic limit [9].This theorem has been extended to chiral matter in Ref. [7], and refined in Ref. [8] (the absence of CME in equilibrium using the chiral kinetic formalism has been obtained in Ref. [6]).There are three elements usually associated to this theorem in chiral matter: The existence of Weyl nodes that always come by pairs [10], [11], (local) gauge invariance and, of course, the assumption that the system is in the equilibrium state.As, we have mentioned, it is known how to break the second condition and drive the system out of equilibrium.Recently, it has been proposed that the first assumption of having pairs of Weyl nodes can be broken in Weyl superconductors, where an external magnetic field indices a gap in one of the Weyl nodes (and its particle-hole conjugate), leaving effectively a single Weyl node [12].However, we stress here that the presence of Weyl nodes is not a strict requirement to the absence (or presence) of CME [7,8].
Interestingly, non-Hermitian fermionic systems appear to be a promising physical environment that can be realized in real experiments.Nowadays, there is a surge of interest in non-Hermitian systems for many different reasons, ranging from very fundamental questions in the quantum (and statistical) theory of fields and the role of topology in non-Hermitian systems [13][14][15], to applied science.Among them, specially interesting are the non-Hermitian systems that display a real spectrum, as the P * T -symmetric systems or the quasi-Hermitian systems.Although non-Hermitian, they display a unitary evolution, and it is possible to define a consistent thermodynamics for them [16].

The model
We will study is a non-Hermitian extension of the massive Dirac model in (3 + 1) dimensions, where, together with the usual mass term m, an anti-Hermitian mass m 5 is introduced [17][18][19]: Here we used the original Dirac notations α ≡ γ 0 γ and β ≡ γ 0 where the Dirac gamma matrices are given in Eq. ( 4), and k describes the momentum of the particle.The advantage of this model is that the two first terms of the right hand side of Eq.( 8) are Hermitian by themselves, so the only non-Hermitian (anti-Hermitian) term is m 5 γ 5 .When m 5 the model Eq.( 8) corresponds to the usual Dirac model for relativistic fermions.Also, it constitutes the low energy model for the bulk states of topological insulators, and when m = 0 or m = m(k) becomes a function with nodal points in momentum space, this model describes Weyl fermions [20].
To date, there is not known experimental realization of an electronic system with the non-Hermitian mass term m 5 βγ 5 .However, we are not aware of any no-go theorem that would forbid this term to appear in open systems.Therefore, we consider the model (8) as a generic system which captures the essential properties of the non-Hermitian mass.Our aim is to show, conceptually, that the equilibrium Chiral Magnetic Effect is, in principle, possible in a generic non-Hermitian system.
It is already stated in the literature that non-Hermitian hamiltonians are not gauge-invariant in general.This can be viewed as the fact that the Noether theorem relating continuous symmetries and conserved currents in field theories, does not hold in non-Hermitian systems [18,19,21].For this reason, there is some arbitrariness when defining a coupling to electromagnetic fields in the Hamiltonian (8).In the present work, we will interested to compare our results with the ones in Hermitian systems, so we will consider the coupling to electromagnetic fields to the model ( 8) with m 5 = 0, which is Hermitian (and the principle of local gauge invariance holds), and later we switch on the non-Hermitian term m 5 βγ 5 .
Also, that we cannot apply the Noether theorem for Hermitian systems to non-Hermitian ones does not imply that conserved currents exist for the latter.
In a conventional Hermitian quantum mechanics, the time dependence of any operator O can be determined in the Heisenberg picture: If we naively try to proceed further following the same steps for a NH system, we will get an unconventional expression for the time variation of the operator O(t) [22,23]: For Hermitian systems, H + = H and we recognize the commutation with H as the condition any operator must satisfy to be a conserved quantity.For non-Hermitian systems, we immediately see that an operator is a conserved quantity if, instead of commuting with the Hamiltonian, it fulfills the quasi-Hermiticity condition: H + O = OH so Ȯ = 0.In the case of the U(1) charge symmetry, it is clear that the generator of this symmetry in Hermitian systems commutes with H but does not satisfy the quasi-Hermiticity condition, so it is not a conserved quantity for non-Hermitian systems.As we will see in next paragraphs, it is possible to find operators that, not commuting with H, they satisfy the quasi-Hermiticity condition, thus defining conserved quantities.
Before computing the CME and CSE for the system at hands, it is convenient to understand what is the physical meaning of the conserved currents associated to the Hamiltonian (8).The most convenient way is the following [24]: The Hamiltonian ( 8) is a quasi-Hermitian Hamiltonian that satisfies the relation Hη = ηH + , where η is some positive definite operator called metric operator.The condition of the metric operator of being positive definite, allows us, among other things, to define a non-unitary similarity transformation S (η = S + S) that maps the NH Hamiltonian H in Eq.( 8) onto a hermitian Hamiltonian Ĥ (for further details we refer to the Appendix (B)).Then, we can find the conserved currents in the auxiliary hermitian Hamiltonian and use the mapping S on these conserved current to find the corresponding conserved quantities in the non-hermitian model.In certain systems, where the similarity matrix S cannot be constructed explicitly, it is still possible to identify certain conserved quantities.These quantities are associated with the operators that are symmetries of the system.Namely, given an operator O, we can construct another operator O = ηO, whose time evolution is described, according to Eqs.(9,10), as follows: We see that this new operator O now possesses a conventional time evolution of the quantum mechanics in the Heisenberg picture.Also, if the original operator O is a symmetry of the problem (that is [H, O] = 0), then the new operator O = ηO defines a conserved quantity as well.This discussion allows us to motivate the use of a bi-orthogonal formulation in our paper.It is clear, indeed, that the expression in Eq.( 11), constructed with the help of the operator O = ηO, follows the standard time evolution in terms of the conjugate wavefunctions ( ψ| , |ψ ).
Alternatively, instead of using the modified operator O , we could had perfectly maintained the operator O and defined a modified conjugate wavefunction ψ| η.The pair ( ψ| η, |ψ ) is called a bi-orthogonal set.We will make use of this formulation in the next sections.
Let us now apply the mentioned results to the model defined in Eq.( 8).The corresponding procedure, developed in Ref. [19], utilizes the metric operator η = 1 + m 5 m γ 5 .It turns our that the Hermitian model associated to the non-Hermitian Hamiltonian H corresponds to a massive Dirac spinor χ with mass M = m 2 − m 2 5 .We construct the U(1) conserved current ĵµ = χγ µ χ associated with the spinor χ.After using the inverted mapping S, we get the corresponding current in the non-Hermitian system in terms of the field ψ + : Since ∂ µ ĵµ = 0, can trivially show that the current j µ corresponds to a conserved quantity ∂ µ j µ = 0. We thus see that the current j µ = η J µ is a conserved current with J µ = γ µ being a symmetry of the Hamiltonian in Eq.( 8).As it will discussed in the next section, the most important consequence of having the conserved current j µ is that we can define a chemical potential µ associated to j 0 = η.
We immediately see that the current is made of a piece proportional to the identity, as it corresponds to an abelian current in the Hermitian case, together with a chiral current weighted by m 5 /m that implies a chiral imbalance.We will show in the next Section and in the Appendix (B) that a system defined by the Hamiltonian ( 8) that exchange particles in a manner defined by this precise chemical potential µ defines is in a trully equilibrium thermal state with non vanishing CME and CSE.

Computation of CSE and CME with biorthogonal quantum mechanics
Here we will tackle the problem using the biorthogonal quantum mechanics formalism [25].Within this formalism, we distinguish between the eigenstates of H: Hψ s = ε s k ψ s , their complex conjugates: ψ + s H + = ψ + s ε s k , the bi-orthogonal states φ s : H + φ s = ε s k φ s , and their complex conjugates, φ + s H = ε s k φ + s .The point is that, because H is not Hermitian, ψ s = φ s , and ψ + s = φ + s .Also, for the same lack of Hermititivity, the states are not orthogonal ψ + s |ψ s = δ ss , where •|• is the standard scalar product in the corresponding Hilbert space.However, the state sets ψ s and φ s form a bi-orthogonal basis: For the model ( 8) we can define a metric operator η, that not only fulfills the quasi-Hermiticity condition, ηH = H + η but it is positive definite.The existence of such operator simplifies the construction of the bi-orthogonal basis sets, since these two bases are related to each other through η: With this particular normalization, we have ψ + s |φ s = φ + s |ψ s = δ ss .For the Hamiltonian at hands (8), such metric operator η is η = 1 + m 5 m γ 5 [19].The existence of a metric operator allowing us to define a well-defined inner product in the corresponding Hilbert space, defines an unitary time evolution of the states, as long as the spectrum is real, so a consistent description of quantum mechanics is allowed for the system, although being non-Hermitian.Also, it is now easy to see that any time operator will evolve using the conventional Heisenberg picture within the bi-orthogonal formalism.
Another relevant consequence of the existence of the metric operator is that η is a conserved quantity, since, as we mentioned, the matrix η fulfills the pseudo-Hermiticity condition (remember Eq.( 10)).Although η does not commute with the non-Hermitian Hamiltonian H [26], it allows for a construction of an unitary evolution.The existence of a conserved quantity makes it possible to define a Lagrange multiplier µ associated to the operator η.Since η is a conserved quantity, that Lagrange multiplier µ plays the role of the chemical potential.Consequently, we can define a new Hamiltonian as it is done in the standard Hermitian statistical mechanics.Of course, due to the non-Hermitian nature of the problem, the conserved quantity does not need to commute with H. Instead, to be conserved, the corresponding operator should satisfies the aforementioned pseudo-Hermiticity condition.However, the existence a common basis between H and any operator O is possible if and only if the operator O commutes with the Hamiltonian H, irrespective of the Hermiticity of H.This fact means that, we will not be able to find a common basis for η and H in terms of the eigenstates of the number operator, as it happens in conventional Hermitian Quantum Mechanics.This problem may be circumvented by building the bi-orthogonal basis which is, in turn, is constructed by diagonalizing the new Hamiltonian (15) H instead of the original Hamiltonian H.
In order to compute the non-Hermitian version of the Chiral Magnetic Effect, we consider the model Eq.( 15) in a classical background of an external constant magnetic field B that points along the third dimension.As it is a trivial exercise to obtain the Landau levels for this model, we are not presenting the details of the derivation.However, we are willing to highlight two properties of these Landau levels.
First, the algebraic wavefunction structure in the non-Hermitian system does not differ from the Hermitian case: The system is translationally invariant along the direction of the magnetic field B, so that the momentum along this direction is a conserved integral of motion.Thus, we may use a standard Fourier transformation of the wavefunction along this direction.The dispersion relation of H for a non-zero chemical potential is presented in Fig. 1(a).
Second, the wavefunctions are highly degenerate as in the Hermitian case, with the same Landau degeneracy.The current along the magnetic field becomes diagonal.Therefore, following a standard procedure used in the Hermitian case, we can integrate over the transverse spatial directions when computing quantum averages.This approach greatly simplifies the calculations and allows us to consider the problem as quasi-one dimensional.
Although from the perspective of constructing a thermal equilibrium ensemble, the lack of Hermiticity in the system might suppose a problem [27], here we will argue that this is not actually the case in our model system under general grounds, since the requirement the one-particle correlation function built from the bi-orthogonal basis satisfies the Kubo-Martin-Schwinger (KMS) periodicity condition: [28][29][30].It is known in the context of quantum statistical mechanics that states satisfying the KMS boundary condition extremize the Von Neumann entropy S = −Tr[ρ log ρ], being ρ the density matrix operator , so they describe equilibrium states [30].The point is to notice that, for non-zero µ, the time evolution of any field operator is done through the exponential of H (and not of H) so then the one-particle correlation function will satisfy the KMS boundary condition and we will be able to built an equilibrium ensemble [31].In the Appendix (B) we provide an explicit proof that the correlation functions of the non-Hermitian system considered in the present work can be mapped to the correlation functions of an equilibrium thermal state, thus satisfying the KMS condition.Also, this fact has been pointed out in the existing literature of non-Hermitian systems [32,33].
In our particular case, using the effective one-dimensional model, we will focus only in the lowest Landau level (LLL) after integrating over the perpendicular coordinates and explicitly writing the Landau level denegeracy ρ = 2πeB 3 the equilibrium thermal average of any observable O will be where G 0 (iω n , k 3 ) is the single-particle propagator in imaginary time: and (ψ s , φ s ) are the bi-orthogonal sets of single-particle eigenstates of the model (15) in presence of an external magnetic field B = B 3 ẑ: The generic label s comprises band labelling, the spin index τ, and the Landau level index N. = 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > = 0.4 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > = 1 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > 2⇡ 2 ⌦ J 3 ↵ e 2 mB 3 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > µ/m < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > = 1 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > = 0.6 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > = 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > µ/m < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > m = 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > m = 0.6 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > m = 0.3 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > For the operators defined as O = ∂H ∂λ , we can generalize the Feynman-Hellmann theorem to the basis (See AppendixA), if the eigenstates are real: obtaining, after performing the Matsubara summation, where n F (x) is the Fermi distribution function in absence of the chemical potential.The chemical potential is part of the spectrum.
For the case of CME, J 3 = ∂ε s µ (k 3 ) ∂k 3 , so The dispersion relation for the LLL (N = 0) sector is (see Fig.
where δ = m 5 m and s = ±1, while for N > 0, we have For the N > 0 Landau levels, the spin degree of freedom τ = ±1 appears explicitly.In Fig.
(1(b)) we have plotted the Landau level spectrum for N = 0 and N > 0. The all important difference between the eigen-energies for N = 0 and N > 0 is that, while ε s,τ,n µ (k 3 ) with N > 0 is an even function of k 3 for any value of m, δ = m 5 /m, and µ, the energy ε s,0 µ (k 3 ) with N = 0 is not.That means that, when taking the derivative with respect to k 3 and integrating over a symmetric interval, the N > 0 Landau levels will not contribute to the integral in (20), but the N = 0 will do.
The result turns out to be (see Fig.( 2)): This is the principal result of this Letter.For non-zero values of the mass m 5 , which is the parameter that controls the non-Hermiticity of H, there is a non-vanishing CME in equilibrium.
The chiral separation effect (CSE) is obtained by computing the average value of the chiral current, represented by the operator J i 5 = eα i γ 5 .We can follow the same route as in the case of the CME.We will add a term b 3 α 3 γ 5 to the Hamiltonian (15) and compute the spectrum in presence of the parameter b 3 .Then, we apply the Hellman-Feynman theorem to it, taking the derivative with respect to b 3 and constructing the expectation value for each Landau level.We send the parameter b 3 to zero after the calculation.It is a lengthy but straightforward calculation to check that for the n > 0 sector, ∂ε s,n µ (k 3 , b 3 )/∂b 3 is an odd function of k 3 in the limit b 3 → 0 for all values of m, m 5 , and µ.This implies that the integral over k 3 is zero and they do not contribute to the CSE.In contrast, for the n = 0 sector, we simply have ∂ε s µ (k 3 , b 3 )/∂b 3 = 1, so We have plotted J 3 5 in Figs.(3a,3b) as a function of δ = m 5 /m for fixed m, and as function of m for fixed δ.Performing the integral we finally have with 1.We note that there is a divergent contribution in the CSE.It is a particular feature of (1 + 1) dimensions that there is a duality between the chiral and charge currents.The charge current operator representing the CME is the chiral density, while the chiral current operator j 5  1 that is relevant for the study of the CSE is the same operator as the charge density.Having the same origin as the standard charge density, we regularize it in the same way.

Conclusions
In the present Letter we have demonstrated that CME in equilibrium is possible when non-Hermitian systems are considered.The key ingredient is to realize that the CME is zero if charge conservation is imposed in the system.However, charge conservation, associated to the U(1) symmetry, is not fulfilled in non-Hermitian systems as it is done in conventional Hermitian ones.
Another fact to pay attention is that there not an unique metric operator associated to the non-Hermitian Hamiltonian fulfilling the pseudo-Hermiticity condition.While there is no practical consequence of this regarding the construction of a bi-orthogonal basis (the average of observables do not depend on any particular choice of the metric operator), this observation is relevant as we can associate different chemical potentials to different metric operators understood as conserved quantities in the non-Hermitian sense.Interestingly, all the metric operators are related to each other by a similarity transformation [24], so we can generalize the results obtained here to other chemical potentials by modifying the spectrum correspondingly.
Finally, to the best of our knowledge, there are not experimental realizations of fermionic non-Hermitian systems with real spectrum to test our predictions.However, there are impressive experimental advances in the area of non-Hermitian P T −symmetric photonic systems and other condensed matter analogs [34,35].In fact, it has been recently proposed the experimental observation of the CME employing superconducting quantum circuit technology, and synthetic magnetic fields [36].We suggest the same experimental setup to test our theory, by extending the experimental setup with equal gain-loss [37].Besides, other topological equilibrium effects similar to the CME have been proposed to occur in electromagnetism [38][39][40][41], being the optical helicity and the optical chirality the electromagnetic symmetries that play the role of the chiral symmetry in ultrarelativistic fermionic systems.There, the biorthogonal formalism have probed to be useful to handle the effect of dissipation and loss in electromagnetism [42,43].The natural question is then to see how the topologically-related responses associated to these symmetries are modified by the presence of non-Hermitian effects.
Let us consider a Hamiltonian H depending on some parameter λ.To ease notation, we will keep the dependence with the generic parameter λ implicit in the eigenstates and eigenvalues.We  As |s and |n are eigenstates of H and H + with the same eigenvalue ε n .Simplifying a little, we finally have: which is the result we wanted to prove.
As demonstrated in [32], the existence of a metric operator η allows us to define the non-unitary mapping S of some quasi-Hermitian Hamiltonian H to an Hermitian partner Ĥ, with Ĥ = Ĥ+ (we will denote the Hermitian partners of operators with the hat symbolˆ): (A6) Also, we can define the Hermitian partner of any operator associated to the quasi-Hermitian system in the same way: This includes the field operators Ψ and Ψ + in the second quantization formalism.As discussed in the main text, the existence of the metric operator η allows us to construct a well behaved scalar product in the Hilbert space and to construct bi-orthogonal basis sets, |n and |n .In this way, we can define the following statistical average (here we will use the suffix bi to denote the statistical average with the bi-orthogonal basis): In the second line we have used η = SS (S + = S in our particular case), and that the eigenstates of the non-Hermitian H are related to the eigenstates of the Hermitian partner Ĥ through ψn = Sψ n .Also, we consider that the states ψ n of the Hermitian partner are conveniently normalized: ψ+ n ψn = 1 .To guarantee the proper normalization of (A8), we need to relate the partition function in the quasi-Hermitian system and its Hermitian partner.This is a particular case of the previous identity, as we can choose O = 1 and obtain the equality of the corresponding partition functions: so we conclude that the averages performed with the bi-orthogonal basis and with the density matrix ρ = e −βH satisfy a KMS boundary condition and thus this state defines a thermal state in equilibrium, since it is trivial to modify the previous reasoning by including the Dyson time ordering operation.Also, this reasoning justifies the definition of the discrete-frequency Green function G 0 (iω n ) in Eq.( 17) of the main text.
y T 4 K / w b v C G I b M J c / D z v z Y w 9 r 9 S c W Z d l v 3 r J t e s 3 b t 7 a u d 2 / c / f e / Q e D 3 Y e n V j W G w o g q r s x 5 S S x w J m H k m O N w r g 0 Q U X I 4 K + d v W v 7 s E o x l S r 5 3 C w 2 F I F P J a k a J i 6 n J Y A 9 f E g P a M q 7 k h b f P s C N N m H g s m r C P S + H n 4 e l k M M w O s m W k 2 y B f g S F a x c l k t / c d V 4 o 2 A q S j n F g 7 z j P t C k + M Y 5 R D 6 O P G g i Z 0 T q Y w j l A S A b b w y 8 + E 9 E m t T O p m k C 7 v 6 1 p P h L U L U U a N I G 5 m N 7 k 2 y T 4 K / w b v C G I b M J c / D z v z Y w 9 r 9 S c W Z d l v 3 r J t e s 3 b t 7 a u d 2 / c / f e / Q e D 3 Y e n V j W G w o g q r s x 5 S S x w J m H k m O N w r g 0 Q U X I 4 K + d v W v 7 s E o x l S r 5 3 C w 2 F I F P J a k a J i 6 n J Y A 9 f E g P a M q 7 k h b f P s C N N m H g s m r C P S + H n 4 e l k M M w O s m W k 2 y B f g S F a x c l k t / c d V 4 o 2 A q S j n F g 7 z j P t C k + M Y 5 R D 6 O P G g i Z 0 T q Y w j l A S A b b w y 8 + E 9 E m t T O p m k C 7 v 6 1 p P h L U L U U a N I G 5 m N 7 k 2 y T 4 K / w b v C G I b M J c / D z v z Y w 9 r 9 S c W Z d l v 3 r J t e s 3 b t 7 a u d 2 / c / f e / Q e D 3 Y e n V j W G w o g q r s x 5 S S x w J m H k m O N w r g 0 Q U X I 4 K + d v W v 7 s E o x l S r 5 3 C w 2 F I F P J a k a J i 6 n J Y A 9 f E g P a M q 7 k h b f P s C N N m H g s m r C P S + H n 4 e l k M M w O s m W k 2 y B f g S F a x c l k t / c d V 4 o 2 A q S j n F g 7 z j P t C k + M Y 5 R D 6 O P G g i Z 0 T q Y w j l A S A b b w y 8 + E 9 E m t T O p m k C 7 v 6 1 p P h L U L U U a N I G 5 m N 7 k 2 H h 7 k E b 9 7 P j x + v X J 8 B z 1 G e 2 g f 5 e g F O k Z v 0 Q k a I Y o + o m / o B / q Z 6 O R T 8 j n 5 8 k e a 9 F Y 1 j 1 A n k q + / A f y M 8 l g = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " y w + 7 4 h G B 1 M M P g T Z D 6 91 p 3 O l M j u Q = " > A A A C 7 n i c b V F N b x M x E H W W r x K + U j h y W T V C K h K q d i s k O F S o g g v H I p G 2 U r y N v N 7 Z x I q / Z H u L I t d 3 f g E 3 x J U T J y T 4 K / w b v C G I b M J c / D z v z Yw 9 r 9 S c W Z d l v 3 r J t e s 3 b t 7 a u d 2 / c / f e / Q e D 3 Y e n V j W G w o g q r s x 5 S S x w J m H k m O N w r g 0 Q U X I 4 K + d v W v 7 s E o x l S r 5 3 C w 2 F I F P J a k a J i 6 n J Y A 9 f E g P a M q 7 k h b f P s C N N m H g s m r C P S + H n 4 e l k M M w O s m W k 2 y B f g S F a x c l k t / c d V 4 o 2 A q S j n F g 7 z j P t C k + M Y 5 R D 6 O P G g i Z 0 T q Y w j l A S A b b w y 8 + E 9 E m t T O p m k C 7 v 6 1 p P h L U L U U a N I G 5 m N 7 k 2 + V + u F D F d K l 4 t q b T D 1 b C Q Q r d j W 0 B k W j E y N b F b 9 5 W u f l l 4 J n X j Q N K o j l z d 8 N S p t N 1 q L D J A H V + 0 1 d S w + M + U z o g h 1 M X d 9 7 G E D 1 S J 2 L 3 y e A 4 u j P P C Y w 6 1 u x r m 2 L D p

3 <
7 0 O d 9 0 d R u c H h 7 k E b 9 7 P j x + v X J 8 B z 1 G e 2 g f 5 e g F O k Z v 0 Q k a I Y o + o m / o B / q Z 6 O R T 8 j n 5 8 k e a 9 F Y 1 j 1 A n k q + / A f y M 8 l g = < / l a t e x i t > k l a t e x i t s h a 1 _ b a s e 6 4 = " + r h h U u 6 a 6 v 3 Y w a r 8 f 2 K v t h 7 Y 6 L I s e 9 c u I K f 4 Q / w 7 9 h H Y K I E + a y b + e 9 m d m d l 1 Z S O A z D X 4 N g 5 9 b t O 3 d 3 7 w 3 v P 3 j 4 6 P H e / s G F M 7 X l M O Z G G n u V M g d S a B i j Q A l X l Q W m U g m X 6 e y s 4 y + v w T p h 9 C e c V 5 A o V m i R C 8 7 Q p z 7 O p q + m e 6 P w M F w G 3 Q b R C o z I K s 6 n + 4 O f c W Z 4 r U A j l 8 y 5 S R R W m D T M o u A S 2 m F c O 6 g Y n 7 E C J h 5 q p s A l z f K t L X 2 e G 0 u x B L q 8 r 2 s b p p y b q 9 R r F M P S b X J d 8 r 9 c q n w 6 N T J b U r T H 5 T D X q u r G d o B p m g l W W N + t / 0 r M 3 y S N 0 F W N o L l X e y 6 v J U V D u 6 X 5 I g s c 5 b y r 5 l b 4 f 1 J e M s s 4 + t U O Y w 1 f u F G + e 9 b E M 8 B 2 E i V N L C H H x S i K r S h K P G n 7 q t S y f 6 r j v 6 r F h q q y b f O 5 O 4 S C d r t D b 9 T x + q i 1 B T j v U A l Z O / Q + R 5 u u b o O L l 4 e R x x + O R q d v V 4 7 v k q f k G X l B I v K a n J J 3 5 J y M C S c F + U a + k x / B + 6 A O F s H N H 2 k w W N U 8 I b 0 I v v 4 G a + P k Q A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " + r h h U u 6 a 6 v 3 Y w a r 8 f 2 K v t h 7 Y 6 L I s e 9 c u I K f 4 Q / w 7 9 h H Y K I E + a y b + e 9 m d m d l 1 Z S O A z D X 4 N g 5 9 b t O 3 d 3 7 w 3 v P 3 j 4 6 P H e / s G F M 7 X l M O Z G G n u V M g d S a B i j Q A l X l Q W m U g m X 6 e y s 4 y + v w T p h 9 C e c V 5 A o V m i R C 8 7 Q p z 7 O p q + m e 6 P w M F w G 3 Q b R C o z I K s 6 n + 4 O f c W Z 4 r U A j l 8 y 5 S R R W m D T M o u A S 2 m F c O 6 g Y n 7 E C J h 5 q p s A l z f K t L X 2 e G 0 u x B L q 8 r 2 s b p p y b q 9 R r F M P S b X J d 8 r 9 c q n w 6 N T J b U r T H 5 T D X q u r G d o B p m g l W W N + t / 0 r M 3 y S N 0 F W N o L l X e y 6 v J U V D u 6 X 5 I g s c 5 b y r 5 l b 4 f 1 J e M s s 4 + t U O Y w 1 f u F G + e 9 b E M 8 B 2 E i V N L C H H x S i K r S h K P G n 7 q t S y f 6 r j v 6 r F h q q y b f O 5 O 4 S C d r t D b 9 T x + q i 1 B T j v U A l Z O / Q + R 5 u u b o O L l 4 e R x x + O R q d v V 4 7 v k q f k G X l B I v K a n J J 3 5 J y M C S c F + U a + k x / B + 6 A O F s H N H 2 k w W N U 8 I b 0 I v v 4 G a + P k Q A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " + r h h U u 6 a 6 v 3 Y w a r 8 f 2 K v t h 7 Y 6 L I s e 9 c u I K f 4 Q / w 7 9 h H Y K I E + a y b + e 9 m d m d l 1 Z S O A z D X 4 N g 5 9 b t O 3 d 3 7 w 3 v P 3 j 4 6 P H e / s G F M 7 X l M O Z G G n u V M g d S a B i j Q A l X l Q W m U g m X 6 e y s 4 y + v w T p h 9 C e c V 5 A o V m i R C 8 7 Q p z 7 O p q + m e 6 P w M F w G 3 Q b R C o z I K s 6 n + 4 O f c W Z 4 r U A j l 8 y 5 S R R W m D T M o u A S 2 m F c O 6 g Y n 7 E C J h 5 q p s A l z f K t L X 2 e G 0 u x B L q 8 r 2 s b p p y b q 9 R r F M P S b X J d 8 r 9 c q n w 6 N T J b U r T H 5 T D X q u r G d o B p m g l W W N + t / 0 r M 3 y S N 0 F W N o L l X e y 6 v J U V D u 6 X 5 I g s c 5 b y r 5 l b 4 f 1 J e M s s 4 7 Y 6 L I s e 9 c u I K f 4 Q / w 7 9 h H Y K I E + a y b + e 9 m d m d l 1 Z S O A z D X 4 N g 5 9 b t O 3 d 3 7 w 3 v P 3 j 4 6 P H e / s G F M 7 X l M O Z G G n u V M g d S a B i j Q A l X l Q W m U g m X 6 e y s 4 y + v w T p h 9 C e c V 5 A o V m i R C 8 7 Q p z 7 O p q + m e 6 P w M F w G 3 Q b R C o z I K s 6 n + 4 O f c W Z 4 r U A j l 8 y 5 S R R W m D T M o u A S 2 m F c O 6 g Y n 7 E C J h 5 q p s A l z f K t L X 2 e G 0 u x B L q 8 r 2 s b p p y b q 9 R r F M P S b X J d 8 r 9 c q n w 6 N T J b U r T H 5 T D X q u r G d o B p m g l W W N + t / 0 r M 3 y S N 0 F W N o L l X e y 6 v J U V D u 6 X 5 I g s c 5 b y r 5 l b 4 f 1 J e M s s 4 W b O E b + B n + B k 8 I I p N w N z 6 + 5 9 x 7 7 X u o E d z 5 N P 0 1 S H Z u 3 L x 1 e / f O 8 O 6 9 + w / 2 9 g 8 O z 5 1 u L I M x 0 0 L b S 0 o c C K 5 g 7 L k X c G k s E E k F X N D 5 6 4 6 / + A D W c a 3 e + Y W B W b O E b + B n + B k 8 I I p N w N z 6 + 5 9 x 7 7 X u o E d z 5 N P 0 1 S H Z u 3 L x 1 e / f O 8 O 6 9 + w / 2 9 g 8 O z 5 1 u L I M x 0 0 L b S 0 o c C K 5 g 7 L k X c G k s E E k F X N D 5 6 4 6 / + A D W c a 3 e + Y W B W b O E b + B n + B k 8 I I p N w N z 6 + 5 9 x 7 7 X u o E d z 5 N P 0 1 S H Z u 3 L x 1 e / f O 8 O 6 9 + w / 2 9 g 8 O z 5 1 u L I M x 0 0 L b S 0 o c C K 5 g 7 L k X c G k s E E k F X N D 5 6 4 6 / + A D W c a 3 e + Y W B Q p K Z 4 h V n x M f U V U 5 l m L f T k B u w p p 3 u j 9 K j d B l 4 G 2 Q r M E K r O J s e D H 7 m p W a N B O W Z I M 5 N s t T 4 I h D r O R P Q D v P G g S F s T m Y w i V A R C a 4 I y 2 e 3 + E m l L f Y 1 4 O V 9 X R u I d G 4 h a d R I 4 m u 3 y X X J / 3 J U x j T V o l x S u M d V s F D S d G M 7 Q B Q u O Z n Z 2 K 3 / S l + 9 L A J X p v G g W F R H r m o E 9 h p 3 + 4 t F F p g X i 6 6 a W R 7 / i V l N L G E + b n m Y K / j I t I z d y 5 D P w b e T r A i 5 g M p f j 7 L c 8 l n t T 9 q + i l r y T 3 X 8 V 3 W 9 o T K 2 D V f d w S W 0 2 x 1 6 o 4 7 X R 6 0 t w E W H a i j b Y f Q 5 2 3 R 1 G 5 w / O 8 o i f v t 8 d P p q 5 f g u e o Q e o 6 c o

Figure 1 .
Figure 1.(color online) (a) Band structure of H at zero magnetic field but finite chemical potential.Contrary to Hermitian systems, the presence of chemical potentials might modify strongly the spectrum.(b) Landau level spectrum for the non-Hermitian model for finite chemical potential µ.Finite values of µ shifts the LLL spectrum (red) not only upwards or downwards, but it is laterally shifted.The lateral shift makes the contribution from LLL to be non-zero for the CME.

Figure 2 .
Figure 2. (color online) CME as a function of µ/m for three values of δ = m 5 /m.The vanishing CME for the Hermitian case, δ = 0 is recovered.

Figure 3 .
Figure 3. (color online) (a): Regularized CSE as a of µ for three values of δ = m 5 /m.We fix the mass parameter to be m = 0.3.(b): CSE as a function of µ for three values of m and fixed δ = 0.6. (1(b))):