Weyl fermions in a non-abelian gauge background and trace anomalies

We study the trace and chiral anomalies of Weyl fermions in a non-abelian gauge background in four dimensions. Using a Pauli-Villars regularization we identify the trace anomaly, proving that it can be cast in a gauge invariant form, even in the presence of the non-abelian chiral anomaly, that we rederive to check the consistency of our methods. In particular, we find that the trace anomaly does not contain any parity-odd topological contribution, whose presence has been debated in the recent literature.


Introduction
In this paper we calculate the trace (and chiral) anomalies of Weyl fermions coupled to non-abelian gauge fields in four dimensions. One of the motivations to study this problem arises from a debate on whether a topological, parity-odd term is present in the trace anomaly of the stress tensor of chiral fermions. We find that it does not.
We start by considering Bardeen's method [1] that embeds the Weyl theory into the theory of Dirac fermions coupled to vector and axial non-abelian gauge fields. Using a Pauli-Villars (PV) regularization we calculate its trace anomaly. As an aside we rederive the well-known non-abelian chiral anomaly to check the consistency of our methods. A chiral limit produces the searched for anomalies of the non-abelian Weyl fermions.

Bardeen's model
We consider the Bardeen's model of massless Dirac fermions ψ coupled to vector and axial non-abelian gauge fields, A a and B a . The lagrangian is given by where D / (A, B) = γ a D a (A, B), with D a (A) = ∂ a + A a + B a γ 5 the covariant derivative for the gauge group G × G. Taking an appropriate limit on the background (by setting A a = B a → Aa 2 ) one finds the theory of left-handed Weyl fermions. We expand the gauge fields on the generators of the gauge group as A a = −iA α a T α and B a = −iB α a T α . The components A α a and B α a are real, and T α denote the hermitian generators in the representation of G chosen for ψ (we allow for the presence of an abelian subgroup, for example one could consider the group U (N ) with the fermion ψ sitting in the fundamental representation) 1 .
This model is classically gauge invariant and conformally invariant. We wish to compute systematically the anomalies. The chiral anomaly is well-known, of course, and we recompute it to test our methods. The main aim is to obtain the trace anomaly.
Let us first review the classical symmetries. The lagrangian is invariant under the G × G gauge transformations. Using infinitesimal parameters α = −iα α a T α and where ψ c = C −1 ψ T is the charge conjugated spinor. The transformations of the gauge fields can be written more compactly as where A a = A a + B a γ 5 andα = α + βγ 5 . The corresponding field strength contains the Bardeen curvaturesF ab andĜ ab In the following we prefer to use the more explicit notation with γ 5 . One can use A α a and B α a as sources for the vector J aα = iψγ a T α ψ and axial J aα 5 = iψγ a γ 5 T α ψ currents, respectively. These currents are covariantly conserved on-shell, with the conservation law reading or, equivalently, as Indeed, under an infinitesimal gauge variation of the external sources A α a and B α a , the action S = d 4 x L varies as 8) and the gauge symmetries implies that J aα and J aα 5 are covariantly conserved onshell, as stated above.
Similarly, to study the stress tensor, it is useful to couple the theory to gravity by introducing the vierbein e µ a and related spin connection. One may verify that the model acquires a Weyl invariance, i.e. an invariance under arbitrary local scalings of the vierbein. This suffices to prove conformal invariance in flat space. The vierbein is used also as an external source for the stress tensor where e denotes the determinant of the vierbein. The Weyl symmetry implies that the stress tensor is traceless on-shell. Indeed, an infinitesimal Weyl transformation with local parameter σ is of the form Then, the full Weyl symmetry implies that the trace vanishes on-shell, T a a = 0.

PV regularization
To regulate the one-loop effective action we introduces massive PV fields. The mass term produces the anomalies, which we will compute with heat kernel methods. We denote by ψ the PV fields as well (for the moment this does not cause any confusion) and add a Dirac mass term to their massless lagrangian in (2.1) It preserves vector gauge invariance but breaks axial gauge invariance. Indeed under (2.2) the mass term varies as where β = −iβ α T α , which shows that the vector gauge symmetry is preserved, leaving room for an anomaly in the axial gauge symmetry. The mass term sources also a trace anomaly, as the curved space version of (3.1) varies under the infinitesimal Weyl transformation (2.10) as However, it preserves the general coordinate and local Lorentz symmetries. One concludes that only axial gauge and trace anomalies are to be expected. Casting the PV lagrangian L P V = L + ∆L in the form where φ = ψ ψ c , allows us to recognize the operators The latter identifies the regulators, as we shall see in the next section.

Regulators and consistent anomalies
Using the Pauli-Villars regularization, we relate the anomaly computation to a sum of heat kernel traces, following the scheme of refs. [3,4] which we briefly review. Starting with a lagrangian for ϕ invariant under a linear symmetry acting also on the backgroud fields contained in the operator T O, one constructs the one-loop effective action Γ by a path integral. The latter is regulated by subtracting a loop of a massive PV field φ with lagrangian (3.5) where it is understood that one should take the decoupling limit M → ∞, with all divergences canceled by renomalization. The anomalous response of the path integral under a symmetry is due to the PV mass term only, as one can define the measure of the PV field to make the whole path integral measure invariant. In a hypercondensed notation, where a term like φ T φ includes in the sum of the (suppressed) indices a spacetime integration as well, the lagrangian in (4.1) is equivalent to the action, and one computes the symmetry variation of the regulated path integral as follows where brackets ... indicate normalized correlators. It is convenient to manipulate this expression further, by using the identity and invariance of the massless action, to cast it in the equivalent form In the derivation we have considered a fermionic theory, used the PV propagator taken into account the opposite sign for the PV loop, and considered an invertible matrix T . In the limit M → ∞ the regulating factor . This substitution allows us to use well-known heat kernel formulae. Obviously, a symmetry remains anomaly free if one finds a symmetrical mass term.
Thus, denoting the anomaly is related to the trace of the heat kernel of the regulator R with insertion of the operator J iδΓ It has the same form of the regulated Fujikawa's trace producing the anomalies [5,6], where J is the infinitesimal part of the jacobian arising from a change of variables in the path integral under a symmetry transformation, and R is the regulator. The limit extracts only the mass independent term (negative powers of the mass vanish in the limit, while positive powers are renormalized away, usually by employing additional PV fields). The PV method guarantees that the regulator R together with the jacobian J produces consistent anomalies, which follows from the fact that one is computing directly the variation of the effective action.
Let us now go back to the specific case of the Bardeen's model, and extract the heat kernel traces that compute the anomalies. For each symmetry we must consider the transformation generated by K and obtain the corresponding form of J.
To start with, the vector current J aα remains covariantly conserved also at the quantum level, as the PV mass term is invariant under vector gauge transformations.
For the axial current, recalling the transformation laws in (2.2), one finds as δT vanishes, while the contribution from δO is also seen to vanish (all possible terms vanish under the Dirac trace). Here, T αT denotes the transposed of T α . Removing the spacetime integration and the local parameters β α from (4. where the remaining trace is the finite dimensional one on the gamma matrices and color space. Here, we find the so-called Seeley-DeWitt coefficients a 2 (R i ) corresponding to the regulators R i associated to the fields assembled into φ The a 2 coefficients are the only ones that survive renormalization and the limit M → ∞.
Similarly, for the Weyl symmetry one uses the transformations in (2.10) to find T a a = − 1 2(4π) 2 tr a 2 (R ψ ) + tr a 2 (R ψc ) (4.12) where now also δT contributes to (4.7), while δO drops out as before. Again, all remaining traces are in spinor and color spaces. Since the mass term is general coordinate and local Lorentz invariant, no anomalies arise in those symmetries.

Anomalies
We are left to compute the anomalies produced by the traces of the heat kernel coefficients a 2 in (4.10) and (4.12), with the regulators (4.11). The heat kernel formulae needed in the calculation are well-known, and for commodity we have reported them in appendix A. The vector symmetry is guaranteed to remain anomaly free by the invariance of the mass term. As a check one may verify, using the explicit traces given in appendix A, that the would-be anomaly vanishes (5.1)

Chiral anomaly
Evaluation of (4.10) produces the chiral anomaly where the remaining trace is only in color space (the trace on gamma matrices has been implemented). P ET s indicate the parity-even terms that take the form They are canceled by the chiral gauge variation of a local counterterm (5.4) and the remaining answer coincides with the famous result obtained by Bardeen [1].

Trace anomaly
Evaluation of (4.12) produces the trace anomaly where CT T s are the cohomologically trivial terms that are canceled by the Weyl variation of the following counterterm where µ, ν are curved indices, and R the Ricci scalar. Of course, one restricts to flat space after variation. The counterterms (5.4) and (5.7) merge consistently into a unique counterterm that in curved space reads (5.8) One may already notice that, on top of the complete gauge invariance of the trace anomaly, there is no parity-odd term present.

Chiral and trace anomalies of Weyl fermions
We are now ready to study the chiral limit of the Bardeen's model, and identify the chiral and trace anomalies of Weyl fermions. We take the limit A a = B a → 1 2 A a , which creates a chiral projector in the coupling to the gauge field, normalized as usual after the scaling. Then,F ab =Ĝ ab → 1 2 F ab (A) and J a = J a 5 → J a = 1 2 (J a + J a 5 ), so that from (5.2) and (5.5) (without the cohomologically trivial terms) we immediately derive the searched for anomalies for the left-handed Weyl fermions coupled to nonabelian gauge fields The chiral anomaly is the standard one, rederived as a check on the methods used here. The trace anomaly is our new result, that verifies the absence of parity-odd terms. It is just half the trace anomaly of non-chiral Dirac fermions.

Conclusions
We have calculated the chiral and trace anomaly in the Bardeen's model of Dirac fermions coupled to non-abelian vector and axial gauge fields, rederiving the famous result for the chiral anomaly and finding the trace anomaly. Then, by a chiral limit we have obtained the chiral and trace anomalies for left-handed Weyl fermions coupled to non-abelian gauge fields.
The main aim of this paper was to find the explicit form of the trace anomaly for Weyl fermions, verifying that it does not contain any parity-odd term proportional to the topological Chern-Pontryagin density. The latter was conjectured to be a possibility in [7], see also comments in [8,9]. However, it was found to be absent in the abelian gauge coupling of a single Weyl fermion [2]. Here we prove that it is absent also in the more general case of the coupling to non-abelian gauge fields. The analogous case of a Weyl fermion on a curved spacetime background has been debated more extensively in the literature, where a topological term proportional to the Pontryagin density had been reported in [10], and confirmed in [11,12], where the concept of a MAT background, that extends the Bardeen construction to curved space, has been developed. However, the topological term was found to be absent in [13], as confirmed also in [14]. We believe that the latter are the correct results. This conclusion indeed finds support from the analogous situation studied in this paper. Also, an analysis of a Dirac fermion on the MAT background, suitably regularized with PV fields, does not seem to produce parity-odd terms in the trace anomaly [15].

Acknowledgments
We wish to thank Loriano Bonora for stimulating discussions.

A The heat kernel
Let us consider a flat D-dimensional spacetime and an operator H of the form where V is a matrix potential and ∇ 2 = ∇ a ∇ a , with ∇ a = ∂ a + W a the gauge covariant derivative satisfying The trace of the corresponding heat kernel has a small time expansion given by (A.3) where the symbol "tr" is a trace on the remaining discrete matrix indices, J(x) is an arbitrary matrix function, and a n (x, H) are the heat kernel, or Seeley-DeWitt, coefficients. They are matrix valued, and the first ones are a 0 (x, H) = 1 where ∇ a V = ∂ a V + [W a , V ], and so on. More details on the heat kernel expansion can be found in [16,17]. They have been computed with quantum mechanical path integrals in [18,19], while a useful report is [20].
In the main text, the role of the hamiltonian H is played by the regulators R ψ and R ψc , and is ∼ 1 M 2 , see eq. (4.8) (here we use a minkowskian set-up). In D = 4 the s-independent term contains a 2 (x, H), which produces the anomalies.
Let us now specialize to the regulator R ψ = −D / 2 (A, B) which is expanded as and contains the Bardeen curvaturesF ab andĜ ab given in (2.5), the covariant derivative D a (A) = ∂ a + A a , and the covariant divergence of Comparing it with the heat kernel operator H in eq. (A.1) allows one to fixes Now the coefficient a 2 (R ψ ) can be made explicit using (A.4). We compute directly the relevant Dirac traces, and list some intermediate results for the reader interested in checking our calculations. Recalling the three contributions in the last line of (A.4), we find (with D a ≡ D a (A)): i) from a 2 = DBDB .
(A.12) The analogous results for a 2 (R ψc ) are obtained by replacing A → −A T and B → B T (and also T α → T αT for the explicit T α appearing in the traces). Their effect is just to double the contribution from a 2 (R ψ ) in the chiral and trace anomalies.