Arxiv:1701.03368V2 [Hep-Lat] 12 May 2021 Hspofrle Ntetcnqeo Inrtasomto [30–33] Transformation Wigner of Technique Supercon the Without U Theory

Home , Lattice field theory

arXiv:1701.03368v2 [hep-lat] 12 May 2021 hspofrle ntetcnqeo inrtasomto [30–33] transformation Wigner of technique supercon the without u theory. systems on to Ho the relied seems (in [27]. proof systems CME theorem This theoretic equilibrium Bloch the field Wey go of the absence of - to the no model theorem the Bloch particular to the the contradict extend within may remain reported CME current nonze equilibrium was CME the the CME total system integrated of the the absence of systems the boundary given the theo particula the to field in In close lattice that but of exist. context exist, not the not does in does investigated CME was of current version CME the [20] recently reported equilibrium in was original observed it was the Although [25] conductivity 20–23]. ordinary [17, to r in contribution and example, [17–19] observe for example, to discussed, for was possibility (see, the collisions physics ion energy - high heavy relativistic the the of context the In olwn h omls eeoe n[8 ntepeetpprw consid t we paper that present indication, the m an in - [28] as 3 in this the developed t over formalism consider the We integration if Following the Namely, undefined. if con current. remains while the current CSE first, gives the current performed axial of is of frequencies calculation calculation the the regulariz in then ultraviolet ambiguity temperatures, the [8 an curre that in is the however, regularization there for argue, lattice expression we of the framework paper, present the because theory in regularization field discussed ultraviolet quantum was no nonregularized it naive Also the [3]. of framework the In potential 1]i n ftemmeso hsfml.I aiet tefi h equilib the in t itself (corresponding manifests field It magnetic family. external this of conside of presence members been the the in has where theor of semimetals matter one Weyl condensed is and of [16] Dirac context discovered the recently in the and physics energy high ∗ [email protected] h S tafis okhstesm rgn(h hrlaoay sthe as anomaly) chiral (the origin same the has look first a at CSE The ntepeetpprw rce hsln frsac n investiga and research of line this proceed we paper present the In h aiyo h o isptv rnpr ffcshsbe ieyd widely been has effects transport dissipative - non the of family The µ h xa urn per ie by given appears current axial the aetecnrbto oteCEo eo n oe fteGree the of poles and zeros become of doublers CSE fermion the the to this contribution for the price case a propos regularization as overlap disappea but the symmetry, of CSE dou version fermion the modified the symmetry the of is contributions chiral the exact when the conventi fermions, the lattice of of presence and fermions the ma Wilson in the of for - expression cases exact conventional particular without its two and by with described - is regularization effect of types two discuss ecnie hrlSprto ffc CE nteltiereg lattice the in (CSE) Effect Separation Chiral consider We hrlSprto ffc nltieregularization lattice in Effect Separation Chiral .Ceeuhisaa2,Mso,175,Russia 117259, Moscow, 25, Cheremushkinskaya B. nttt o hoeia n xeietlPhysics, Experimental and Theoretical for Institute ..hiuo n ..Zubkov M.A. and Z.V.Khaidukov .INTRODUCTION I. Dtd a 3 2021) 13, May (Dated: j 5 k = − 4 π 1 2 ǫ ijk 0 etoa euti h umto vrteMatsubara the over summation the if result ventional µF to siprat edmntaeta ihu it without that demonstrate We important. is ation oCEcretmyapa.I a demonstrated, was It appear. may current CME ro h edstrength field the o tde o oti lrvoe iegne.I the In divergences. ultraviolet contain not does nt drcnl,weeteei h xc chiral exact the is there where recently, ed 18.Tepsil perneo uheet in effects such of appearance possible The [1–8]. y icvrdDrcsmmtl 2] twssonthat shown was it [24], semimetals Dirac discovered lr aclec te.Aohrexample Another other. each cancel blers ij sesfrin.Ti silsrtdb the by illustrated is This fermions. ssless o ioos h ucetaayia ro of proof analytical sufficient The rigorous. - non s eo ntecneto odne atrtheory matter condensed of context the In zero. s y twsage,ta h qiiru ukCME bulk equilibrium the that argued, was It ry. nloelpfrin.A h aetime same the At fermions. overlap onal hrlsmer.I h atrcs this case latter the In symmetry. chiral httepsil xsec ftecia magnetic chiral the of existence possible the that therein). eferences mnai efre rt h xrsinfrthe for expression the first, performed is omenta h S a icse eety o xml,in example, for recently, discussed was CSE the S a icse,i atclr nrlto to relation in particular, in discussed, was CSE utvt)wspeetdb n fu n[8 29]. [28, in us of one by presented was ductivity) ple oteltierglrzdqatmfield quantum regularized lattice the to applied ergru atc euaiaini ob used. be to is regularization lattice rigorous he e 91] h hrlsprto ffc (CSE) effect separation chiral The [9–15]. red eteeulbimCEo h aegrounds. same the on CSE equilibrium the te eiea 2] eie,i a rud that argued, was it Besides, [26]. semimetal l eo fteGenfnto.I this In function. Green the of zeros ucincne ahother. each cancel function n rCEo h ai fWge transformation Wigner of basis the on CSE er s hsi lutae ytenaive the by illustrated is This rs. lrzdqatmfil hoy We theory. field quantum ularized ,weei a rud htteCEneeds CSE the that argued, was it where ], ,i 58 sn aiu ueia methods numerical various using [5–8] in r, susdrcnl ohi h otx fthe of context the in both recently iscussed emdli osdrda ml u finite but small at considered is model he imter ihmsls ia fermion, Dirac massless with theory rium ∗ ee,tewyteato f[7 re to tried [27] of author the way the wever, o-cle hrlMgei ffc,which Effect, Magnetic Chiral called - so F ij n h riaychemical ordinary the and ) (1) 2 technique [30, 31] applied to the Green functions. Unlike the case of the CME in general case the coefficient in the linear response of the axial current to the external magnetic field is not a topological invariant. However, it appears, that the coefficient in the CSE current standing at the product of magnetic field and chemical potential approaches the topological invariant when mass of the fermion tends to zero. This allows us to derive the conventional expression for the CSE current. Thus the link between the CSE effect in lattice regularization and momentum space topology is established. It is worth mentioning, that Momentum space topology is a powerful method, which was developed earlier mainly within condensed matter theory. It allows to describe in a simple way, for example, the stability of the Fermi points, the Anomalous Quantum Hall effect, and the fermion zero modes on vortexes (for the review see [34, 35]). Recently certain aspects of momentum space topology were discussed in the framework of the four - dimensional lattice gauge theory (see for example [36–38]). The paper is organized as follows. In Sect. II following [28, 29] we describe how the technique of the Wigner transformation may be applied to the lattice regularized quantum field theory. In Sect. III we consider the CSE in conventional lattice regularizations (Wilson fermions and overlap fermions). In Sect. IV we consider the CSE in the case of exact chiral symmetry - for the naive lattice regularization with 16 doublers and in the lattice regularization with deformed overlap fermions, where instead of the 15 doublers the zeros of the Green function appear. In Sect. V we demonstrate the ambiguity in the calculation of the CSE current that takes place in the naive continuum theory. In Sect. VI we end with the conclusions.

II. LATTICE FERMIONS IN THE PRESENCE OF EXTERNAL U(1) GAUGE FIELD

A. Lattice models in momentum space

In this section we consider briefly the lattice models in momentum space following the methodology of [28, 29]. For the more detailed description of the method and for the references see [28, 29]. In the absence of the external gauge field the partition function of the theory defined on the infinite lattice may be written as D d p T 1 Z = DψDψ¯ exp ψ¯ (p) − (p)ψ(p) (2) − G Z ZM |M| where is the volume of momentum space , D is the dimensionality of space - time, while ψ¯ and ψ are the Grassmann|M| - valued fields defined in momentumM space . is specific for the given system. For example, the model with 3 + 1 D Wilson fermions corresponds to that hasM theG form G 1 (p)= i γkg (p) im(p) − (3) G − k − Xk where γk are Euclidean Dirac matrices defined in chiral representation 0 1 0 iσi γ4 = , γi = , i =1, 2, 3 1 0 iσi 0 − 1 0 γ5 = 0 1 − i where σ is the Pauli matrix. gk(p) and m(p) are the real - valued functions (k =1, 2, 3, 4) given by g (p) = sin p , m(p)= m(0) + (1 cos p ) (4) k k − a a=1,2,3,4 X The fields in coordinate space are related to the fields in momentum space as follows dDp ψ(r)= eiprψ(p) (5) ZM |M| At the discrete values of r corresponding to the points of the lattice this expression gives the values of the fermionic field at these points, i.e. the dynamical variables of the original lattice model. However, Eq. (5) allows to define formally the values of fields at any other values of r. The partition function may be rewritten in the form

T 1 Z = DψDψ¯ exp ψ¯ (rn) − ( i∂r)ψ(r) (6) − G − r=rn r Z Xn h i 3

Here the sum in the exponent is over the discrete coordinates rn. However, the operator i∂r acts on the function ψ(r) deﬁned using Eq. (5). In order to derive Eq. (6) we use identity −

eipr = δ(p) (7) |M| r X In the particular case of Wilson fermions we may rewrite the partition function in the conventional way as

Z = DψDψ¯ exp ψ¯T (r )( )ψ(r ) (8) − m Drn,rm n r ,r Z Xn m with

1 i i (0) x,y = [(1 + γ )δx+ei,y + (1 γ )δx ei,y] + (m + 4)δxy (9) D −2 − − i X

Here ei is the unity vector in the i - th direction.

B. Introduction of the gauge ﬁeld

Gauge transformation of the lattice ﬁeld takes the form

ψ(r ) eiα(rn)ψ(r ) (10) n → n In case of Wilson fermions the U(1) gauge ﬁeld is typically introduced as the following modiﬁcation of operator D:

1 i iA = [(1 + γ )δ e x+ei,y Dx,y −2 x+ei,y i X − i iAx ei,y (0) +(1 γ )δx ei,ye ] + (m + 4)δxy (11) − − Here Ax,y = Ay,x is the gauge ﬁeld attached to the links of the lattice. In the same way the gauge ﬁeld is typically incorporated− into the models of solid state physics. One can easily check, that Eq. (8) may be rewritten as

dDp Z = DψDψ¯ exp ψ¯T (p) ˆ(i∂ ,p)ψ(p) (12) − Q p Z ZM |M| Here

1 ˆ = − (p A(i∂ )) (13) Q G − p while the pseudo - differential operator A(i∂p) is defined as follows. First, we represent the original gauge field A(r) as a series in powers of coordinates r. Next, variable r is substituted in this expansion by the operator i∂p. Besides, in Eq. (13) each product of the components of p A(i∂p) is subsitituted by the symmetric combination (for the details see [28]). As it was mentioned above, for the− case of Wilson fermions the formulations of Eq. (12) and Eq. (8) are exactly equivalent. For the other regularizations there may be the difference, but it manifests itself in the terms that are proportional to the field strength times a2 (here a is the lattice spacing). Those extra terms may be neglected in continuum limit. Therefore, for any regularization we accept Eq. (13) as the definition of the model in the presence of external U(1) gauge field.

C. Electric current

Electric current is defined as the response of effective action log Z to the variation of external Electromagnetic field. This gives [28] −

D 1 k d p ˜ ∂ ˜(0) − j (R)= D Tr G(R,p) G (R,p) (14) − (2π) ∂pk ZM h i 4 where the Wigner transformation of the Green function is expressed as:

ipr G˜(R,p)= e− G(R + r/2, R r/2) (15) − r=r Xn while the Green function itself is 1 G(r , r )= DΨ¯ DΨ Ψ(¯ r )Ψ(r ) (16) 1 2 −Z 2 1 Z 1 exp Ψ(¯ r ) − ( i∂ − n G − r r Xn h h A(r))Ψ(r) − r=rn i i At the same time

G˜(0)(R,p)= (p A(R)) (17) G − The method of Wigner transformation in the application to the lattice models was developed in [28, 29] following its original formulation speciﬁc for the theory in continuous space - time [30–33]. In [28] the following expression was derived for the linear response of the electric current to external electromagnetic ﬁeld: 1 j(1)k(R)= ǫijkl A (R), (18) −4π2 Ml ij = Tr ν d4p (19) Ml l Z 1 1 i ∂ − ∂ ∂ − νl = 2 ǫijkl G G G (20) −3!8π G ∂pi ∂pj ∂pk h i

III. LATTICE REGULARIZATION WITH BROKEN CHIRAL SYMMETRY

A. Linear response of chiral current to external magnetic ﬁeld

In this section we consider the linear response of the chiral current to external Electromagnetic field. For the field system in continuous coordinate space this response may easily be calculated using Feynman diagrams. For the field system in lattice regularization this response may be calculated following the approach of [28, 29, 38] that was briefly described above. In continuum theory the naive expression for the chiral current is ψγ¯ µγ5ψ . Several different definitions for the particular lattice regularization may give this expression in the naive continuumh i limit. The evident choice of the definition of current in lattice regularization is the functional derivative over axial gauge field of the effective action. The latter field, in turn, may be defined through the covariant derivative, which acts on the left - handed and the right - handed fermions via opposite charges. For the particular choice of the lattice model with exact chiral symmetry (i.e. when commutes or anti - commutes with γ5) this definition gives the expression similar to that of [28, 29, 38]: G

D 1 5k d p 5 ˜ ∂ ˜(0) − j (R)= D Tr γ G(R,p) G (R,p) (21) (2π) ∂pk ZM h i where

G˜(0)(R,p)= (p A(R)) (22) G − Actually, we are able to adopt this definition to any lattice theory even without the exact chiral symmetry. One can easily check that in the naive continuum limit this definition gives ψγ¯ kγ5ψ . In order to regularize our expressions for the case of truly massleh ss fermionsi let us use the finite temperature version of the lattice theory. With the periodic boundary conditions in the spatial directions and anti-periodic in the imaginary time direction, the lattice momenta will be 2π pi (0, 2π); p4 = (n4 +1/2) (23) ∈ Nt 5 where i =1, 2, 3 while n =0, ..., N 1. Temperature is equal to T =1/N , in lattice units 1/a, where a is the lattice 4 t − t spacing. Thus the imaginary frequencies are discrete p4 = ωn = 2πT (n +1/2), where n = 0, 1, ...Nt 1, while the axial current (also in lattice units) is expressed via the Green’s functions as follows: −

Nt 1 − 3 5k i d p 5 1 j = T Tr γ ( (ω , p)∂ − (ω , p) −2 (2π)3 G n pi G n n=0 Z X 1 ∂ (ω , p)∂ − (ω , p))F (24) pj G n pk G n ij

B. The linear response of the chiral current to chemical potential and to external magnetic ﬁeld

Now let us consider the system without exact chiral symmetry. Recall, that the exact chiral symmetry is to be broken if we want to describe one Dirac fermion, which is related to Nielsen - Ninomiya theorem. For deﬁniteness, we may discuss ﬁrst Wilson fermions. But this is not necessary, and the results of this subsection are valid for any lattice models. We introduce the chemical potential in the standard way ωn ωn iµ. The derivative of the current with respect to µ gives → −

ijk j5k = N F µ (25) 4π2 ij with

Nt 1 − 3 ijk 1 d p 5 1 = T ∂ Tr γ (ω , p)∂ − (ω , p) N − 2 (2π) ωn G n pi G n n=0 Z X 1 ∂ (ω , p)∂ − (ω , p) (26) pj G n pk G n We assume that the singularities of the Green function may appear only at the ﬁnite sequence of values of ω = ω(0),ω(1), ... that do not coincide with the Matsubara frequencies. In the limit T 0 the sum over Matsubara frequencies becomes the integral that is regularized as follows →

ijk = ijk(ω(k) +0)+ ijk(ω(k) 0) (27) N − N3 N3 − Xk where (ω ) is given by N3 n 3 ijk 1 d p 5 1 (ω )= Tr γ (ω , p)∂ − (ω , p) N3 n −2 (2π)2 G n pi G n Z 1 ∂ (ω , p)∂ − (ω , p) (28) pj G n pk G n It is clear, that in the model, in which there are no poles or zeros of the Green function in the presence of exact chiral symmetry the linear response of the axial current to magnetic ﬁeld is the sum of topological invariants, i.e. it cannot be changed under the continuous deformations of the model. However, in general case when γ5 does not (anti) commute with the Green function, the terms in this expansion are not topological invariants. We may also rewrite

3 ijk 1 d p 5 1 = Tr γ (ω, p)∂ − (ω, p) N 2 (2π)2 G pi G ZΣ 1 ∂ (ω, p)∂ − (ω, p) (29) pj G pk G where Σ is the 3D hypersurface of inﬁnitely small volume that embraces the singularities of the Green function concentrated at the Fermi surfaces (or Fermi points). It has the form of the two inﬁnitely close pieces situated at ω = 0. The advantage of this representation is that Eq. (29) becomes the topological invariant if γ5 anti - commutes with± the Green function in a small vicinity of its poles. Namely, in this case we may represent

ijk [ij]k ǫ 1 5 1 = Tr γ (ω, p)d − (ω, p) N 3! (2π)2 G G ZΣ 1 d (ω, p) d − (ω, p) ∧ G ∧ G Here Σ is a small closed surface in 4D space. It may have an arbitrary form surrounding the singularities. 6

C. Regularization with Wilson fermions

For the case of Wilson fermions the singularities of the Green function may appear at ω =0, π (for m(0) > 0 they appear at ω = 0 only). Then the limit T 0 gives → ijk = ijk(+0) + ijk ( 0)+ ijk (π 0) ijk(π +0) (30) N −N3 N3 − N3 − − N3 The interesting particular case is when parameter m(0) vanishes. In this case at µ = 0 the only Fermi point appears at p = 0 and on Σ we have γ5, 0. In this particular case { G} ≈ ijk = ǫijk N which gives the regular expression for the Chiral Separation Eﬀect of Eq. (1). In order to conﬁrm this prediction we use also the numerical methods. Namely, we take Eq. (26) and calculate numerically for the component 123 the N integral over 3 - momenta and the sum over ωn using MAPLE package. It is seen, that at Nt the answer tends to 1 as it should. → ∞ In the presence of nonzero mass m(0) > 0 the situation is changed, and the poles of the Green function do not appear while µ

D. Overlap fermions

Let us discuss the regularization using overlap fermions [39]. The massless overlap Dirac operator is deﬁned as

+ 1/2 = m(1+ˆ ( m) ( ( m) ( m))− ) (31) Do D − D − D − where (m(0)) is the dimensionless Wilson-Dirac operator given by Eq. (9), where we substitute negative value of the D (0) mass parameter m = m. Operator o obeys Ginsparg - Wilson relation, which may be written in the following form − D 5 1 5 γ − ,γ = (32) {Do } m It is called sometimes ” exact ” chiral symmetry on the lattice. However, this statement is not precise, and actually 1 5 the conventional overlap propagator does not obey the exact chiral symmetry, which is o− ,γ = 0. In momentum space we have {D }

1 1 − = iγ C + , (33) Do − µ µ 2m with ˆ2 1 kµ kµ Cµ(p)= , A = m (34) 2m 2 2 2 − kµ + A + A q and p k = sin(p ), kˆ = 2sin( µ ) (35) µ µ 2 In this model the fermion Green function

1 = − G Do has the only pole at ω =0, p = 0. At p = (n1π,n2π,n3π,n4π) with integer ni =0, 1 such that n1 + n2 + n3 + n4 =0 1 6 we have the value (p)= 2m . Again, the aboveG consideration may be applied to the model in this regularization, and we have the expression for the linear response to the magnetic ﬁeld given by Eq. (25) with of Eq. (29). In particular, in the continuum limit when m 0, and µ =0 we get ijk = ǫijk that results in the usualN expression for the CSE current of Eq. (1). → N 7

IV. LATTICE REGULARIZATION WITH EXACT CHIRAL SYMMETRY

A. Naive lattice fermions

In this section we consider the case of the lattice model with exact chiral symmetry, when

γ5, =0 { G} The simplest example of such a system is given by the naive lattice fermions with the Green function in momentum space of the form

1 (p)= i γkg (p) im(0) − (36) G − k − Xk k where γ are Euclidean Dirac matrices while gk(p) are the real - valued functions (k =1, 2, 3, 4) given by

gk(p) = sin pk (37)

In this model at m(0) = 0 instead of one massless Dirac particle in continuum limit there are 16 massless particles. In this case the linear response of the chiral current to external magnetic field and chemical potential is given by Eq. (25) with of Eq. (29). The contributions of the doublers differ due to the orientation of effective vierbein, i.e. the correspondingN low energy effective theory for massless particle has the one - particle Euclidean lagrangian

= e eµγai L | | a ∇µ where

( 1)n1 0 0 0 −0 ( 1)n2 0 0 e eµ = − | | a  0 0 ( 1)n3 0  − n4  0 0 0( 1)   −  with ni =0, 1. As a result the contributions to CSE current of those 16 doublers cancel each other. Thus, unlike the case of the previous section, in the continuum limit when m 0, and µ = 0 all doublers contribute the sum thus giving = 0, and cancelling the overall CSE current. → N

B. Modiﬁed overlap fermions

It was proposed (see, for example, [39]) to redeﬁne the overlap fermions as follows

1 1 = − G Do − 2m In this case the chiral symmetry is exact

,γ5 =0 (38) {G }

But the price for this is that at p = (n1π,n2π,n3π,n4π) with n1 + n2 + n3 + n4 = 0 we have the vanishing value of the Green function (p) = 0. 6 The zeros of the GreenG function are in many aspects similar to poles. In particular, they contribute to CSE in such a way, that the total current vanishes. Let us deﬁne

2m(√k2 + A2 + A) f(k2)= k2 Then k γµ = i µ (39) G − f(k2)k2 8

At ﬁnite temperatures we have

Nt 1 iT − j5k = − (ω )ǫijkF (40) 2π N3 n ij n=0 X

At any value of n the functional 3(ωn) is the topological invariant, i.e, it is not changed under any variation + δ if during such modificationN the poles or zeros of do not appear. It is given by G → G G G 3 1 ijk d p 5 1 (ω )= ǫ Tr γ (ω , p)∂ − (ω , p) N3 n 2 3! (2π)2 G n pi G n Z × 1 ∂ (ω , p)∂ − (ω , p) (41) pj G n pk G n Therefore, in the model, in which there are no poles or zeros of the Green function in the presence of exact chiral symmetry the linear response of the axial current to magnetic field is the sum of topological invariants, i.e. it cannot be changed under the continuous deformations of the model. The pole of the Green function at finite temperature may appear if there exists such integer n that ω = 2π (n +1/2) = π. This gives equation n Nt

2n +1= Nt which has a solution for odd values of Nt. Therefore, for simplicity in the following we assume that Nt is even. The ordinary chemical potential cannot cause the appearance of poles or zeros of the Green function, which is seen from the following consideration. Again, let us assume that the chemical potential appears as the imaginary contribution to Matsubara frequency. In this case the poles or zeros of the Green function may appear if:

3 sin2(ω iµ)+ sin2(p )=0 (42) n − i Xl=1 We obtain the following system

1 cos(2ω )ch(2µ)+2 3 sin2(p )=0 − n l=1 i (43) sh(µ) sin(2ω )=0 ( n P

The second equation has a solution in real variables. We obtain ωn = π/2 or ωn = π. In the former case the ﬁrst equation does not have solutions. The latter case is realized for odd values of Nt only, and then the poles of the Green function appear as the solution of equation

3 2 1+2 sin (pi) = ch(2µ) Xl=1

However, as above we can always choose the even value of Nt, and therefore the poles of the Green function do not appear if we modify the value of µ. Therefore, we are able to calculate 3(ωn) for vanishing µ, and the result gives the answer for ﬁnite µ. This calculation is represented in Appendix A,N and as expected it gives the vanishing CSE current.

V. NAIVE CONTINUUM EXPRESSIONS FOR THE CSE CURRENT

A. The integration over 3 - momenta before the summation over Matsubara frequencies

Above we considered the CSE effect using rigorous lattice regularizations. Also we feel this instructive to present here the discussion of chiral separation effect in the framework of naive continuous field theory. We will see, that there is an ambiguity in this consideration, which is reflected by the lattice constructions with and without exact chiral symmetry. Let us consider the propagator of massless non-interacting Dirac fermions 1 p G(ωn, )= µ (44) γ pµ 9

Here pµ = (ω, p). We will consider the case when the magnetic ﬁeld is directed along the z axis (i.e F12 = B). The expression for the chiral current has the form: −

∞ d3p ω j5z =4Ti n B (45) (2π)3 ((ω )2 + p2)2 n= n X−∞ Z and after integration over 3 - momenta we arrive at

5z 2 ∞ j =4T iπ sign(ωn)B (46) n= X−∞ Thus, formally, in the case of massless fermions the chiral current is equal to the sum of the integer numbers. If, as a result of the interaction in medium, G is changed as ω f(ω ),p g(p ), then the result depends on sign f. n → n i → i We introduce the chemical potential in the standard way ωn ωn iµ. In this case we need an analytical continuation of the sign function. We may try to use, for example,→ the rule−

sign(ω iµ) = sign(Re(ω iµ)) (47) n − n − And then from this naive consideration the conclusion may be drawn, that the chemical potential does not inﬂuence the CSE current. Below we will see, that the formal expressions in the continuum theory will lead to the diﬀerent answer if the summation over the Matsubara frequencies is performed before the integration over momenta.

B. The integration over 3 - momenta after the summation over Matsubara frequencies

In the previous section we have shown that the axial current in an external magnetic field is expressed as the sum of the topological invariants multiplied by the field strength. This expression would become the exact result, but only if the theory does not contain divergences. We may extract another result from the above expressions. Namely, let us first perform the initial summation over the frequencies and only after that - the integration over momenta

∞ d3p ω iµ j5z =4Ti n − B = (48) (2π)3 ((ω iµ)2 + p2)2 n= n X−∞ Z −

dz d3p z z µ =2 th( − )B (49) 2iπ (2π)3 (z2 p2)2 T ZC Z − Where C is the contour that surrounds poles of the hyperbolic tangent function. We use the relation

m 1 1 d − m resz=z0f(z)= lim m 1 f(z)(z z0) (50) (m 1)! z z0 dz − − → − to calculate the value of the integral using the theory of residues. We deform contour C in such a way that it surrounds the points z = z . After this deformation we have ± 0 d3p z d z µ j5z = 2 th( − ) − (2π)3 (z + p)2 dz 2T |z=p Z z dh z µ + th( − ) z= p (51) (z p)2 dz 2T | − − i We can write the equation in this form because: d z =0 (52) dz (z p)2 ∓ Thus we can rewrite the integral as

d3p 1 d p µ 1 d p + µ j5z = 2B ( th( − ) th( )) (53) − (2π)3 4p dp 2T − 4p dp 2T Z 10 and after substitution d3p =4πp2dp we find that B Bµ j5z = dp(n (p µ) n (p + µ)) = (54) −2π2 f − − f −2π2 Z This expression coincides with the conventional expression Eq. (1) (see, for example, [16]) and also it coincides with the result obtained above using Wilson fermions and conventional overlap fermions. It is still protected from the renormalization of 3 - momentum (pi g(pi)). However, although the result of the→ rigorously regularized theory (using, say, lattice Wilson fermions) is reproduced by the approach of the present subsection, we would like to emphasize once again, that this approach itself is not self - consistent, and its application to the other problems may be limited. In particular, let us consider the modification of the system, which leads to the replacement of iωn by a function f(iωn) such that it tends to iωn at large n. Looking at Eq. (45) we may come to the conclusion, that such a modification cannot change the value of the axial current. However, Eq. (49) is not invariant under the substitution iωn f(iωn). This demonstrates once again, that the rigorous ultraviolet regularization is needed in order to calculate →the response of axial current to external field strength.

VI. CONCLUSIONS AND DISCUSSIONS

In the present paper we discuss the Chiral Separation Effect both in the framework of the naive continuum nonregularized quantum field theory and of the lattice regularized theory. In both cases we also regularize the theory using finite temperatures. We demonstrate, that the naive continuum formulation suffers from the ambiguities related to the order of taking the integral over the 3 - momenta and the sum over Matsubara frequencies. If the Matsubara frequencies are summed first, then the divergencies are not encountered and the conventional expression for the chiral current in the presence of external magnetic field is reproduced. At the same time if the 3 - momenta are integrated first, then the resulting expression is given by the sum of Eq. (46), where each term is equal to either 1 or 1. This sum is not well defined, but each term in this sum is indepedent of the chemical potential. − This ambiguity points out that although certain computational schemes of the CSE current do not encounter the ultraviolet divergences, the theory should be considered in the ultraviolet regularization in order to obtain rigorous results. Therefore, we consider several types of lattice regularization. First of all, we considered the naive lattice regularization, where 16 doublers represent the independent physical excitations. Those excitations differ by the orientation of effective vierbein, and as a result their contributions to the CSE current cancel each other. Recently the modification of the regularization using overlap fermions was proposed (see, for example, [39]), in which massless physical excitation appears at ω = p = 0 while at the positions of the other 15 doublers (of naive lattice fermions), the zeros of the Green function appear. As a result the exact lattice chiral symmetry is obeyed just like in the case of naive lattice fermions. The physical meaning of the zeros of the Green function remains unclear, but it is discussed in certain publications (mostly in the framework of condensed matter theory). We demonstrate, that the contribution of those zeros of the Green function to the CSE cancels the contribution of the physical massless excitation. Thus in both considered cases of the lattice theory with exact chiral symmetry the CSE does not appear. Typically, in the lattice models the exact chiral symmetry is broken, which is the price for the elimination of the fermion doublers. We consider the two particular cases of such conventional regularization - the case of lattice Wilson fermions and the case of conventional overlap fermions. In both cases the massless excitation appears at ω = p = 0 only, the other doublers disappear, and also there are no zeros of the Green function. The price for this is the absence of the exact chiral symmetry. However, in the case of overlap fermions there is the Ginsparg - Wilson relation instead. In both these regularizations we observe the emergence of the Chiral Separation Effect. The corresponding current tends to its conventional expression Eq. (1) in continuum limit of the model with massless fermions. In the case, when the theory describes massive fermions with mass m, the CSE current is absent at µ

We conclude, that in the physical regularizations with Wilson and overlap fermions the conventional CSE emerges. At the same time we suppose, that the model with the modiﬁed overlap fermions [39] with exact chiral symmetry Eq. (38) is unphysical. Although the zeros of the Green function do not contribute to the ordinary perturbation expansion on the same grounds as the physical excitations, they contribute the topological quantities responsible for the CSE in the same way as the fermion doublers. In this respect the modiﬁed overlap fermions with Eq. (38) are similar to the naive lattice fermions with 16 doublers, and they do not possess the CSE.

ACKNOWLEDGEMENTS

MAZ kindly acknowledges numerous discussions with G.E.Volovik. K.V.Z. is greatful to A.Yu.Kotov for discussions. Both authors are greatful to M.N.Chernodub for useful discussions, and to LMPT, the University of Tours, where this work was initiated, for kind hospitality. The present work was supported by Russian Science Foundation Grant No 16-12-10059.

APPENDIX A. AXIAL CURRENT FOR THE MODIFIED OVERLAP FERMIONS

Here we calculate the chiral current for the version of overlap fermions with exact chiral symmetry ,γ5 = 0. Let us use the following expression for the axial current {G }

Nt 3 5k i d p 5 1 j = T r(γ (ω , p)∂ − (ω , p) −2 (2π)3 G n pi G n n=1 Z X 1 ∂ (ω , p)∂ − (ω , p))F (55) pj G n pk G n ij We substitute Eq. (39) to this expression:

N i t d3p γµk j5k = T r(γ5 µ ∂i(γν k f(k2)) −2 (2π)3 k2f(k2) ν n=1 X Z γλk ∂j ( λ )(∂kγρk f(k2)))F (56) k2f(k2) ρ ij It may be written as

N t d3p f 2(k2)k ∂ik ∂j k ∂kk j5k = 2i ǫρµνλ µ ν λ ρ F (57) − (2π)3 f 2(k2)k4 ij n=1 M X Z µ We introduce the notation gµ = k . And the expression for the axial current is given by √k2

Nt 2i j5k = ǫµνλρ g dg dg dg ǫijkF (58) −3!(2π)3 µ ν ∧ λ ∧ ρ ij n=1 M X Z To calculate the topological invariant 1 (ω )= ǫµνλρ g dg dg dg (59) N3 n 12π2 µ ν ∧ λ ∧ ρ ZM we will use the method of [28], and the following parametrization:

g4 = sin(α),gi = kicos(α) (60)

2 π π where i =1, 2, 3 and k = 1, while α [ − , ]. Thus i i ∈ 2 2 P dg = cos(α)dα, dg = dk cos(α) k sin(α)dα (61) 4 i i − i And 3 = ǫijk cos2(α)k dα dk dk N3 12π2 i ∧ j ∧ k ZM 12

3 1 cos(2α) = ǫijk k ( − )dα dk dk = 12π2 i 2 ∧ j ∧ k ZM

3 α 1 = ǫijk k d( + sin(α)) dk dk 12π2 i 2 4 ∧ j ∧ k ZM

3 ijk α 1 = 2 ǫ ki( + sin(α))dkj dkk (62) − 12π ∂Ω 2 4 ∧ Xl Z In this expression ∂Ω is the small vicinity of point yl of momentum space where vector kl is undeﬁned. The absence π of the singularities of gk implies that α 2 at such points. Thus we see that the expression under→± the integral is the total derivative. We can rewrite it in the form

1 = sign(g (y ))Res(y ) (63) N3 −2 4 l l Xl Where we have used the notation [29] 1 Res(y )= ǫijk g dg dg (64) l 8π i j ∧ k Z∂Ω

It is worth mentioning, that this symbol obeys l Res(yl) =0. At each n the value of sign g4 is constant. Therefore, 3(ωn)=0 for any n. N P

[1] K. Landsteiner, E. Megias and F. Pena-Benitez, “Anomalous Transport from Kubo Formulae,” Lect. Notes Phys. 871 (2013) 433 [arXiv:1207.5808 [hep-th]]. [2] M. N. Chernodub, A. Cortijo, A. G. Grushin, K. Landsteiner, and M. A. Vozmediano, “A condensed matter realization of the axial magnetic effect”, Phys. Rev. B 89, 081407(R) (2014) [arXiv:1311.0878]. [3] E. V. Gorbar, V. A. Miransky, I. A. Shovkovy and P. O. Sukhachov, “Chiral separation and chiral magnetic effects in a slab: The role of boundaries,” Phys. Rev. B 92 (2015) 24, 245440 doi:10.1103/PhysRevB.92.245440 [arXiv:1509.06769 [cond-mat.mes-hall]]. [4] V. A. Miransky and I. A. Shovkovy, “Quantum field theory in a magnetic field: From quantum chromodynamics to graphene and Dirac semimetals,” Phys. Rept. 576 (2015) 1 [arXiv:1503.00732 [hep-ph]]. [5] S. N. Valgushev, M. Puhr and P. V. Buividovich, “Chiral Magnetic Effect in finite-size samples of parity-breaking Weyl semimetals,” arXiv:1512.01405 [cond-mat.str-el]. [6] P. V. Buividovich, M. Puhr and S. N. Valgushev, “Chiral magnetic conductivity in an interacting lattice model of parity- breaking Weyl semimetal,” Phys. Rev. B 92 (2015) 20, 205122 doi:10.1103/PhysRevB.92.205122 [arXiv:1505.04582 [cond- mat.str-el]]. [7] P. V. Buividovich, “Spontaneous chiral symmetry breaking and the Chiral Magnetic Effect for interacting Dirac fermions with chiral imbalance,” Phys. Rev. D 90 (2014) 125025 doi:10.1103/PhysRevD.90.125025 [arXiv:1408.4573 [hep-th]]. [8] P. V. Buividovich, “Anomalous transport with overlap fermions,” Nucl. Phys. A 925 (2014) 218 doi:10.1016/j.nuclphysa.2014.02.022 [arXiv:1312.1843 [hep-lat]]. [9] S. Parameswaran, T. Grover, D. Abanin, D. Pesin, and A. Vishwanath, “Probing the chiral anomaly with nonlocal transport in Weyl semimetals, Phys. Rev. X 4, 031035 (2014) [arXiv:1306.1234]. [10] M. Vazifeh and M. Franz, “Electromagnetic response of weyl semimetals”, Phys. Rev. Lett. 111, 027201 (2013) [arXiv:1303.5784]. [11] Y. Chen, S. Wu, and A. Burkov, “Axion response in Weyl semimetals”, Phys. Rev. B 88, 125105 (2013) [arXiv:1306.5344]. [12] Y. Chen, D. Bergman, and A. Burkov, “Weyl fermions and the anomalous Hall effect in metallic ferromagnets”, Phys. Rev. B 88, 125110 (2013) [arXiv:1305.0183]; David Vanderbilt, Ivo Souza, and F. D. M. Haldane Phys. Rev. B 89, 117101 (2014) [arXiv:1312.4200]. [13] S. T. Ramamurthy and T. L. Hughes, “Patterns of electro-magnetic response in topological semi-metals”, arXiv:1405.7377. [14] A. A. Zyuzin and A. A. Burkov, “Topological response in Weyl semimetals and the chiral anomaly,” Phys. Rev. B 86 (2012) 115133 [arXiv:1206.1868 [cond-mat.mes-hall]]. [15] Pallab Goswami, Sumanta Tewari, Axionic field theory of (3+1)-dimensional Weyl semi-metals, Phys. Rev. B 88, 245107 (2013), arXiv:1210.6352 [16] “Anomalous Axion Interactions and Topological Currents in Dense Matter”,Max A. Metlitski and Ariel R. Zhitnitsky,Phys. Rev. D 72, 045011 13

[17] D. E. Kharzeev, “The Chiral Magnetic Effect and Anomaly-Induced Transport,” Prog. Part. Nucl. Phys. 75 (2014) 133 doi:10.1016/j.ppnp.2014.01.002 [arXiv:1312.3348 [hep-ph]]. [18] D. E. Kharzeev, J. Liao, S. A. Voloshin and G. Wang, “Chiral Magnetic Effect in High-Energy Nuclear Collisions — A Status Report,” arXiv:1511.04050 [hep-ph]. [19] D. E. Kharzeev, “Chern-Simons current and local parity violation in hot QCD matter,” Nucl. Phys. A 830 (2009) 543C doi:10.1016/j.nuclphysa.2009.10.049 [arXiv:0908.0314 [hep-ph]]. [20] A. Vilenkin, Phys. Rev. D 22, 3080 (1980) [21] K. Fukushima, D.E. Kharzeev, H.J. Warringa, Phys.Rev.D 78:074033,2008 [22] D.T. Son and N. Yamamoto, ”Berry curvature, triangle anomalies, and chiral magnetic effect in Fermi liquids”, Phys.Rev.Lett.109:181602,2012, arXiv:1203.2697. [23] D. E. Kharzeev and H. J. Warringa, “Chiral Magnetic conductivity,” Phys. Rev. D 80 (2009) 034028 doi:10.1103/PhysRevD.80.034028 [arXiv:0907.5007 [hep-ph]]. [24] Q. Li, D. E. Kharzeev, C. Zhang, Y. Huang, I. Pletikosic, A. V. Fedorov, R. D. Zhong, J. A. Schneeloch, G. D. Gu, and T. Valla, arXiv:1412.6543. [25] H. B. Nielsen and M. Ninomiya, “Adler-bell-jackiw Anomaly And Weyl Fermions In Crystal,” Phys. Lett. B 130 (1983) 389. doi:10.1016/0370-2693(83)91529-0 [26] M. M. Vazifeh and M. Franz, Phys. Rev. Lett. 111, 027201 (2013) [27] N. Yamamoto, Phys. Rev. D 92, 085011 (2015). [28] M. A. Zubkov, “Wigner transformation, momentum space topology, and anomalous transport,” Annals Phys. 373, 298 (2016) [arXiv:1603.03665 [cond-mat.mes-hall]]. [29] M. A. Zubkov, “Absence of equilibrium chiral magnetic effect,” Phys. Rev. D 93, no. 10, 105036 (2016) [arXiv:1605.08724 [hep-ph]]. [30] E.P. Wigner, ”On the quantum correction for thermodynamic equilibrium”, Phys. Rev. 40 (June 1932) 749–759. doi:10.1103/PhysRev.40.749 [31] C. Zachos, D. Fairlie, and T. Curtright, ”Quantum Mechanics in Phase Space” ( World Scientific, Singapore, 2005) ISBN 978-981-238-384-6 [32] Robert G Littlejohn, ”The semiclassical evolution of wave packets”, Physics Reports Volume 138, Issues 4–5, May 1986, Pages 193-291 [33] Berezin, F.A. and M.A. Shubin, 1972, in: Colloquia Mathematica Societatis Janos Bolyai (North-Holland, Amsterdam) p. 21. [34] G.E. Volovik, The Universe in a Helium Droplet, Clarendon Press, Oxford (2003). [35] G. E. Volovik, “Topology of quantum vacuum,” Lecture Notes in Physics 870 (2013) 343 [arXiv:1111.4627 [hep-ph]]. [36] M.A.Zubkov, G.E.Volovik, Momentum space topological invariants for the 4D relativistic vacua with mass gap, Nuclear Physics B (2012) doi:10.1016/j.nuclphysb.2012.03.002, arXiv:1201.4185 [37] M.A.Zubkov, Generalized unparticles, zeros of the Green function, and momentum space topology of the lattice model with overlap fermions, Phys. Rev. D 86, 034505 (2012) [38] M. A. Zubkov, “Momentum space topology of QCD,” arXiv:1610.08041 [hep-th] [39] “Overlap Quark Propagator in Coulomb Gauge QCD and the Interrelation of Confinement and Chiral Symmetry Breaking”, M. Pak, M. Schröck , Phys. Rev. D 91 (2015) 074515. [40] M. Puhr and P. V. Buividovich, “A numerical study of non-perturbative corrections to the Chiral Separation Effect in quenched finite-density QCD,” arXiv:1611.07263 [hep-lat].