On Torsion Contribution to Chiral Anomaly via Nieh-Yan Term

In this note we present a solution to the question of whether or not, in the presence of torsion, the topological Nieh-Yan term contributes to chiral anomaly. The integral of Nieh-Yan term is non-zero if topology is non-trivial; the manifold has a boundary or vierbeins have singularities. Noting that singular Nieh-Yan term could be written as a sum of delta functions, we argue that the heat kernel expansion cannot end at finite steps. This leads to a sinusoidal dependence on the Nieh-Yan term and the UV cut-off of the theory (or alternatively the minimum length of spacetime). We show this ill-behaved dependence can be removed if a quantization condition on length scales is applied. It is expected as the Nieh-Yan term can be derived as the difference of two Chern class integrals (i.e. Pontryagin terms). On the other hand, in the presence of a cosmological constant, we find that indeed the Nieh-Yan term contributes to the index with a dimensionful anomaly coefficient that depends on the de Sitter length or equivalently inverse Hubble rate. We find similar result in thermal field theory where the anomaly coefficient depends on temperature. In both examples, the anomaly coefficient depends on IR cut-off of the theory. Without singularities, the Nieh-Yan term can be smoothly rotated away, does not contribute to topological structure and consequently does not contribute to chiral anomaly.


I. INTRODUCTION
There has been a long debate with conflicting results on the contribution of torsion to chiral anomaly via the Nieh-Yan term [1][2][3][4][5][6][7][8][9][10].One suspects that in the presence of torsion, the parity odd Nieh-Yan four-form could potentially result in non-conservation of chiral current 1 .However, as the vierbein is chosen dimensionless and torsion has mass dimension one, the anomaly coefficient of the Nieh-Yan term is dimensionful.Naive computations imply that the mass parameter in the coefficient is the UV cut-off of the theory, namely the anomaly term is regulator dependent which is not normally accepted.It is the source of confusion that led some authors to completely abandon the Nieh-Yan part.
In this short paper, we revisit this problem and we show that depending on the underlying spacetime, the Nieh-Yan term may or may not contribute to chiral anomaly.In fact for smooth geometry, the Nieh-Yan contribution is completely rotated away.On the other hand, in geometries with singularities in vierbein structure, it needs more elaboration.In this case, in a locally flat gauge condition, the Nieh-Yan term can be written as a sum of delta functions about singular points.We will argue that, in the Fujikawa method with a regulator, the Heat kernel expansion cannot be truncated at finite steps.It gives a sinusoidal dependence on the Nieh-Yan term as well as on the UV cut-off.This ill contribution can be removed via applying a quantization condition on the parameters.We argue that it is natural as the Nieh-Yan term can be written as a difference between two Pontryagin terms.However, we find that in the presence of cosmological constant, there is a residual and indeed the Nieh-Yan term contributes to the index and chiral anomaly.The anomaly coefficient is proportional to the square of the Hubble rate or inverse de Sitter radius.Similarly, we find that in a finite temperature field theory there is a contribution to chiral anomaly through the Nieh-Yan term.In this case, the anomaly coefficient is proportional to temperature squared.In both scenarios, the anomaly is related to a mass scale, in fact the IR scale in the theory.
In the next section, after a brief review of chiral anomaly and torsional geometry we present the main results and finally conclude in the section III.

II. TORSION AND CHIRAL ANOMALY
We briefly review the calculation of chiral anomaly using the Fujikawa method.In this method one uses the invariance of partition function under a chiral rotation, which is merely a field redefinition, and uses the fact that the measure in field space does not remain invariant under this redefinition.Considering an infinitesimal chiral rotation given by parameter β(x) as the Dirac Lagrangian transforms as

arXiv:2308.00578v2 [hep-th] 30 Dec 2023
The partition function can be written as Under the chiral transformation (1) the measure dψd ψ is not left invariant and transforms as where we have and ψ n are the eigenmodes of the Dirac operator.
Invariance of partition function under this transfromation further implies where J µ5 = ψγ µ γ 5 ψ is the chiral current.
In general this sum is not well defined as it stands since this is equal to trγ 5 .δ 4 (0), therefore it is common to use a regulator to make it well defined.The typical regulator that one can use is exp − λ 2 n M 2 where λ n are the eigenvalues of the Dirac operator and M is some UV cut-off that is sent to infinity in the calculation.
On the other hand we can see that the integral of nonconservation of chiral current is determined by a topological quantity, namely the index of the Dirac operator i / D on a manifold M. Considering the integral of (7) we see that only the zero modes contribute to the integral and we can write where n + and n − are respectively the number of positive and negative chirality zero modes.
Using this background we can write the index of the Dirac operator as where ψ n are the eigenstates of the Dirac operator and we have used the Fujikawa method with a Gaussian cutoff fixed by M to regularize the sum.Anomaly computations imply that the chiral current non-conservation is proportional to parity odd Pontryagin and Nieh-Yan [11][12][13][14] four forms where the field strengths (curvature and torsion respectively) are which satisfy Bianchi identities DF ab = 0 and DT a = F ab ∧ e b .Here e a = e a µ dx µ is the vierbein and A is a connection 1-form.
We note the appearance of a dimensionful coefficient M in the second term which is also expected on dimensional reasons as the vierbein is dimensionless.If we decompose the Lorentz connection into spin-connection plus contorsion i.e.A ab = ω ab + K ab , then the curvature is written as b is Riemann curvature and torsion would be T a = K ab ∧ e b .Then, the Pontryagin term is computed as where the axial vector field S is the fully antisymmetric part of torsion (a.k.a.H-torsion) or in component form S µ = ϵ µνρσ T νρσ .This is the part of torsion that couples to fermions through the Dirac operator [15] / D = / D 0 − i 8 / Sγ 5 ( / D 0 is the torsion-free operator).Moreover, the axial vector S can be decomposed into a transverse pseudo-vector S ⊥ µ (namely ∂ µ S ⊥µ = 0) and a pseudo-scalar σ Finally, the Nieh-Yah term which only receives a contribution from the pseudo-scalar part of torsion is where 1 is the unit (volume) four-form.
On the other hand, the Nieh-Yan term can be written as the difference of two Chern class integrals [14][3].We consider two SO (1,4) connections where l is a length scale.Then, the field strengths are We find that the difference between corresponding Chern class integrals is proportional to the Nieh-Yan term Since these are topological terms with integer integrals, this leads to a quantization condition on the Nieh-Yan term that has a dependence on the length scale l.We think a natural choice for this length scale where we want to have dimensionless geometry [16,17] is the Planck's length considered to be the minimum length of spacetime.

A. Manifolds with Trivial Topology
In this section we consider the contribution of nonsingular Nieh-Yan term to chiral anomaly.Pseudo-scalar σ has a non-singular behavior and using the redefinition ψ(x) = e i 8 σ(x)γ 5 ϕ(x) we can see that the index does not depend on the total derivative part of axial torsion.In particular, we find By considering the Ward identity we note that the field redefinition is like a chiral transformation the measure in numerator and denominator is multiplied by the same factor thus we find We conclude the total derivative part of axial torsion doesn't contribute to chiral anomaly.
Moreover, if we consider the Gaussian regulator in Fujikawa method with where / D only has the S ⊥ µ part, we find that the contribution from torsion to chiral anomaly is Using ∇ µ S ⊥µ = 0 and the Bianchi identity we find Apparently, the first two terms depend on our choice of σ and S ⊥ µ .There is a residual symmetry transformation thus we can choose a gauge such that the sum of first two terms vanishes.Therefore, the chiral anomaly would be We provide an explanation for why calculating anomaly with the typical regulator e − / D 2 M 2 does not lead to the expected result.If we define a transformation using the partition function we can show that the anomaly of this transformation is So we see that using the regulator is equivalent to calculating the anomalous current ⟨∂ µ ( ψγ We can show that this current is not the same as the original current by looking at a symmetry that is present in the original calculation which is missing here.The transformation leaves ⟨∂ µ ( ψγ µ γ 5 ψ)⟩ invariant but the above current is not left invariant under this transformation.We note that even without the above explanation, we can redefine the chiral current such that the total derivative contributions like ∇ µ K µ = ∇ µ (S 2 S µ ) + ... are removed from the anomaly without changing its topological structure.This is also true for the contribution of the non-singular part of the Nieh Yan term as this part behaves like a total derivative and is not a topological term.

B. Manifolds with Non-Trivial Topology
We begin this section by an observation that on a manifold with singularities in its vierbeins, the Nieh-Yan term can be written as a sum of delta functions [9] T where l i are some length scales and x i are singular points.In the following, for the sake of simplicity, we ignore spacetime curvature and only consider the contribution of torsion to anomaly via the Nieh-Yan term.In the Fujikawa method, the regulator has an extra piece due to the Nieh-Yan term Then, we find the torsion contribution to the index as Given the delta function form of Nieh-Yan term, we cannot stop the sum at the first order term in □σ in the expansion and we must consider all terms.Then, the torsion contribution to the index is computed as In order to calculate this term we use a simple method.
Using the expansion of sine function we write Using the definition of delta function d 4 xδ 4 (x)f (x) = f (0) and the approximation δ 4 (0) = 1 ϵ 4 with ϵ related to the minimum length of spacetime (and therefore related to the UV cut-off) we can write and considering M = k ϵ with k some order 1 number we find ind ⊃ k 4 4π 2 sin( Thus, in order to remove this term in the index, we demand a quantization condition as where n i are integer numbers.We note that, if we had incorrectly kept only the first order term in the expansion of e i □σ 8M 2 , then we would have found a contribution to the index as which would be cut-off dependent and is not accepted.However, if we keep all terms in the expansion we find which is finite given that ϵ ∼ 1 M and vanishing if the quantization condition (38) is met.
Same result can be obtained via another approach.We recall that the zero modes of the Dirac operator with torsion ψ 0 n and the zero modes of the Dirac operator without torsion ϕ 0 n are related as With a finite torsion term, this simple reasoning states that the total derivative part of pseudo-vector torsion does not contribute to the index.However with singular point, the outcome is different.Given (31) we write as we have singular behavior at some points.Taking care of the divergent behavior close to one of the singularities we find 43) the index is computed as Given a minimum length scale ϵ we may write The first term in ( 45) is proportional to the torsion-free Index.The second part is divergent though.Again, if we demand a quantization condition for l i as where n i are integers, then the divergent piece will be absent.
In passing, we emphasis that this quantization condition is expected as the Nieh-Yan term could be seen as the difference between two Chern classes [2,14] where here we replace the length scale in (18) with the minimum length of spacetime.One might therefore conclude that adding torsion to gravity has a close relationship with spacetime discretization.In particular, if we choose l to be the minimum length of spacetime, we find the quantization condition the sinusoidal dependence of index to the Nieh-Yan term will be absent.
In the following we consider two cases both with a natural IR scale and show that conterary to what previously expected the coefficient of the Nieh-Yan term will be related to the IR scale of the theory.In this context we use the IR cut-off as the lower limit for the physical energy levels.In the case of a positive cosmological constant this natural scale is given by the cosmological constant.This can be expected if we consider de Sitter spacetime as a compact Euclidean space with effective radius l.In the presence of temperature we also have a natural IR scale which is expected given the fermionic Matsubara frequencies ω n = (2n + 1)πT .

Adding a Cosmological Constant
In the absence of the cosmological constant, as the Nieh-Yan term could be seen as the difference of two Pontryagin terms, we considered two Lorentz connections in (16) and derive where n i is an integer and we have chosen ϵ = l.
In the presence of the cosmological constant Λ ∼ l −2 0 with de Sitter length l 0 , we propose to change the second Lorentz connection in (16) as Then, the field strengths are We find that the difference between corresponding Chern class integrals is proportional to the Nieh-Yan term as (53) Then, following the same line of reasoning we find the quantization condition In the previous section, we computed the torsion contribution to anomaly Upon using (54) we find which is similar to the result suggested by [3].Therefore, we see from (40) that the coefficient of Nieh-Yan term in index should be the cosmological constant as an IR cut-off of the theory and we can write Note that the appearance of a potentially non-integer index is a result of spacetime quantization and is evident from ( 9).When we take the limit x ′ → x there is always a mismatch due to discreteness of spacetime.On the other hand due to the quantization condition (54) the final result for the index is equivalent to replacing the UV cut-off in the calculation with the IR scale 1 l0 .This is equivalent to summing over the eigenmodes of the Dirac operator with eigenvalues less than the IR scale in the index i.e. these modes effectively play the role of zeromodes in a theory with an IR scale.

Finite Temperature Effect
Finally, we study the effect of finite temperature on chiral anomaly in a background with torsion.It is in particular interesting in condensed matter experiments where torsion is realized as dislocations or defects in a sample [18][19][20].In the presence of temperature we use the concept of periodic Euclidean time in the definition of delta function.We have where β = 1 T is inverse temperature.We can approximate the delta function as where ϵ is related to the minimum length of spacetime.Then, after calculating the sum in (58) we find which is plotted in Fig. 1.Then, we find The torsion contribution to the index is computed as ). ( 62) Using (61) and M 2 = 2π ϵ 2 , we find It is consistent with the results presented in [21,22] if we require n i to be odd.

III. CONCLUSIONS
In this short note, we argued that there could exist an extra torsional contribution to chiral anomaly through Nieh-Yan term if the manifold has singular points and the field theory has a natural IR cut-off.Moreover for that to be a sound result, we found that there must exist a quantization condition on the minimum length scale in the present of torsion.We studied two scenarios which come with an IR cut-off; one with a cosmological constant and one in a finite temperature.In both cases we found a finite contribution to chiral anomaly proportional to the IR cut-off.These scenarios are in particular interesting as they are respectively relevant in early Universe cosmology and condensed matter experiments.The latter has been already studied in literature (see for instance [21,22]).In a future work, we will study the effect of this new torsional contribution to chiral current anomaly during inflation.
Acknowledgements This work is supported by the research deputy of Sharif University of Technology.We are thankful to Hessamaddin Arfaei for fruitful discussions.

APPENDIX A
In this appendix, we argue that we can subtract the effect of non-topological terms like ∇ µ (S 2 S µ ), ∇ µ (□S µ ) and the non-delta function part of the Nieh-Yan term ∇ µ S µ by redefining the current or equivalently adding counter terms to the Lagrangian.These counter terms include S 4 , S µ □S µ and M 2 S 2 which correspond to redefining the chiral current by adding terms proportional to S 2 S µ , □S µ and M 2 S µ .This subtraction is allowed as far as the topological behaviour of chiral anomaly is not altered.
We explain how we can subtract the effect of the extra terms using these counter terms.Considering the Lagrangian where a, b, c are some numbers and variating with respect to S µ we find Therefore we see that by an appropriate choice of a, b, c we can subtract the effect of the mentioned finite terms from chiral anomaly.
In passing we note that the presence of a coupling S µ C µ in the effective Lagrangian in [15], where C µ is related to Pontryagin as P = 2∇ µ C µ , is expected since we can write J µ5 = δL ef f δSµ and the C µ term in this variation is a confirmation of the presence of Pontryagin term in chiral anomaly.
In the case of delta-function type Nieh-Yan term we can write As an example we consider the case of a sphere as a compact manifold with two such terms one at the north pole and one at the south pole say x 1 and x 2 .We need to cover the sphere with two patches each including one of these singularities.We can write where S + and S − are respectively defined in patch 1 and patch 2. So the integral depends on the difference of S defined in the intersection of the two patches.Now given that S µ ≈ and considering the behaviour of da µ as we take the limit x → ∞ (that we assume to be the boundary of each patch) we find as expected that this integral receives a non-zero contribution from deltafunction type Nieh Yan term.
Using the same argument we can show that terms like S 2 S µ da µ and □S µ da µ receive no contributions from delta-function type Nieh-Yan term and therefore can be subtracted from anomaly by redefining the chiral current.
We show that this reasoning is in fact a result of our previous argument that a total derivative axial torsion should not contribute to anomaly as far as the topological behaviour is not altered.Our original Lagrangian is Here the S 2 term appears from torsional Ricci scalar.We have seen that when we use the regulator e − /