Topological fermion condensates from anomalies

We show that a class of fermion theory formulated on a compact, curved manifold will generate a condensate whose magnitude is determined only by the volume and Euler characteristic of the space. The construction requires that the fermions be treated as K\"{a}hler-Dirac fields and the condensate arises from an anomaly associated with a $U(1)$ global symmetry which is subsequently broken to a discrete subgroup. Remarkably the anomaly survives under discretization of the space which allows us to compute the condensate on an arbitrary triangulation. The results, being topological in character, should hold in a wide range of gravitationally coupled fermion theories both classical and quantum


Introduction
Kähler-Dirac fermions have been proposed as an alternative to Dirac fermions in [1] and play an important role in certain supersymmetric field theories where they have been used to construct lattice theories with exact supersymmetry [2]. They also bear a close connection to staggered fermions on regular hypercubic lattices [1,3].
In this paper we will show that Kähler-Dirac fermions possess another important property; on compact even dimensional spaces with non-zero Euler character they exhibit an anomaly which breaks an exact U (1) symmetry and generates a non-zero fermion condensate.
Rather remarkably this anomaly can be reproduced exactly in a discrete version of the theory, a result which is tied to the structure of the Kähler-Dirac operator and homology theory. We use this fact to compute the anomaly and condensate directly on an arbitrary triangulation. where Ω = (ω 0 , ω 1 , . . . ω p , ψ D ) is a collection of p-forms (antisymmetric tensors). This is the Kähler-Dirac equation and was first proposed as an equation describing fermions in [4]. The action of d and d † on these component forms is given by [1]

Kähler-Dirac fermions in the continuum
An inner product of two such Kähler-Dirac fields is then given by .1 can be obtained by variation of the action where Ω is an independent (in Euclidean space) Kähler-Dirac field. From the Kähler-Dirac field Ω one can build a matrix Ψ using the Clifford algebra of gamma matrices In flat space it is straightforward to show that this fermionic matrix satisfies the usual Dirac equation and describes 2 D/2 degenerate Dirac spinors. Remarkably, eqn. 2.1 holds on any smooth manifold so that the coupling of Kähler-Dirac fields to gravity is quite different from that of Dirac fermions. This difference can be made manifest if one introduces vielbeins γ µ = e a µ γ a in eqn. 2.6 and replaces the derivative in eqn. 2.1 with an appropriate covariant derivative containing the spin connection ω ab µ . This yields a equation which couples the four spinors encountered in flat space through gravitational interactions. It is easily seen that this is not (four copies of) the usual Dirac equation in curved space.
In this paper will exploit the fact that Kähler-Dirac fermions do not require the introduction of this extra spinor structure and work directly with eqn. 2.1.
Consider the linear operator Γ which transforms the p-form fields according to It is straightforward to verify that this operator anticommutes with the Kähler-Dirac operator {K, Γ} = 0. This in turn allows us to define a U (1) symmetry of the massless Kähler-Dirac action Ω → e iθΓ Ω (2.9) This symmetry, which is the analog of chiral symmetry for the Dirac theory, ensures that any perturbative mass renormalization of the theory is multiplicative in the bare fermion mass on any background space. It also guarantees that non-zero eigenvalues of K come in pairs (iλ n , −iλ n ) with corresponding eigenvectors φ n and Γφ n where Kφ n = λ n φ n . (2.10) Of course the symmetry will be broken or anomalous if the fermion measure is not invariant under Γ. Performing such a transformation yields a Jacobian J = det e 2iΓθ . Thus the existence of the anomaly depends on a non-vanishing trace for Γ where the states with non zero eigenvalue are orthogonal and do not contribute while n + and n − denote the number of zero modes with eigenvalue ±1 respectively 1 . All of this discussion mirrors the usual arguments given for chiral fermions 2 and a derivation of the anomaly can be gotten in the continuum after the theory is suitably regulated. However, in this case a simpler method is available. The Kähler-Dirac operator can be discretized using a triangulation of the space in such a way as to preserve the anomaly exactly.

Kähler-Dirac fermions on the lattice
A triangulation of a D-dimensional manifold is constructed by gluing together D-simplices.
A p-simplex consists of (p + 1) vertices with every vertex connected to every other vertex. We will restrict our discussion to triangulations in which each D-simplex is equilateral with edge length a. Such triangulations appear in so-called dynamical triangulation models of quantum gravity [5,6]. It is straightforward to generalize this discussion to triangulations with variable edge lengths. A p-simplex contains (p + 1) boundary components which are (p − 1)-simplices. The simplest triangulations or simplicial complexes are constructed by gluing D-simplices together along their faces with each face common to two and only two D-simplices. The possible vertex labelings of a given D-simplex fall into two categories according to whether they consist of an odd or even permutation of the ordered list of vertices. These two categories are called orientations. To define the Kähler-Dirac operator requires that the triangulation be oriented which means that the orientation of each D-simplex be chosen in such a way that a given face is held with opposite orientations in the two D-simplices that contain it.
There is a natural mapping of differential forms on a continuum manifold to such an oriented discrete triangulation. The idea is to associate a continuum p-form with a p-cochain 1 Γ and K commute on the λ = 0 subspace and hence zero modes can be labeled with a given eigenvalue of Γ.
2 In flat space it is easy to see that Γ acts like an axial rotation on the spinors corresponding to columns of the matrix Ψ.
which is a real or complex linear function over p-simplices f (C p ) where an individual p-simplex is denoted C p [7,8]. The operators d and d † have corresponding discrete analogs in terms of so-called co-boundary ∂ and boundary ∂ operators that act on such p-simplices. For example the boundary operator ∂ acts on a p-simplex consisting of the (ordered) vertices [a 1 , . . . , a p ] as where the hat on vertexâ i indicates that the vertex is omitted. Thus the operation of ∂ on some p-simplex yields an oriented sum of (p − 1)-simplices in its boundary. Applied to a p-cochain f (C p ) this yields a mapping to a (p − 1)-cochain f (C p−1 ) via the relation where I(p, p − 1) is the incidence matrix for the triangulation which takes the value +1 if C p−1 is contained in the boundary of C p with the correct orientation, −1 if it is contained in the boundary with the opposite orientation, and is zero otherwise. Similarly the co-boundary operator ∂ is defined by From these definitions one can verify that ∂ 2 = ∂ 2 = 0. The discrete Kähler-Dirac operator is defined by which squares to (minus) the discrete Laplacian. Furthermore it should be clear that {Γ, K D } = 0 still so there is an exact U (1) symmetry present just as for the continuum theory acting now on p-cochains via ω p → e iθ(−1) p ω p . For any given triangulation the fermion measure is just Under the U (1) symmetry this transforms by the phase factor e i2θχ where If χ = 0 this breaks the U (1) down to Z 2χ , which depends only on the Euler characteristic and is otherwise independent of the triangulation. Since the non-zero eigenvalues of the lattice operator K D again come in complex conjugate pairs this phase must arise from the zero modes of the discrete Kähler-Dirac operator, which must again be either even or odd under Γ. Homology theory guarantees that such exact zero modes exist even in the discrete system and are simultaneously zero modes of the discrete Laplacian − = ∂∂ + ∂∂ = I(p, p − 1)I T (p − 1, p) + I T (p, p + 1)I(p + 1, p) (3.7) By the Hodge theorem (which applies equally to both discrete and continuum theories) we know that n + − n − = χ, which is consistent with the previous computation of the phase change of the measure. Notice that the anomaly vanishes in odd dimensions since χ = 0 there and n + = n − .

Condensates
From this point we confine our discussion to the sphere S D with χ = 2 although it is trivial to generalize to arbitrary Euler characteristic. In this case there are precisely 2 zero modes of the Laplacian associated with the 0 and D-forms (co-chains) both of which have Γ = +1. Let us now consider the partition function on a discrete triangulation T of the sphere with N D D-simplices and N 0 vertices in the presence of fermionic sources η, η of Kähler-Dirac type, paying careful attention to the zero modes that arise in the massless limit.
where the prime denotes that zero modes satisfying Φ i 0 = 0 have been omitted and corresponds to the i th zero mode component of the source. Notice that the first non-zero contribution to Z as m → 0 occurs by expanding the zero mode sources to quadratic order and yields The structure of this expression resembles the usual 't Hooft vertex seen in QCD [9] and generates a non-local four fermion interaction in the theory. Differentiating with respect to the sources will yield a non-zero four fermion correlation function that survives the m = 0 limit and breaks the U (1) symmetry down to Z 4 .
where x, y, z, w range over the discrete set of coordinates on the triangulation. To compute this correlator we consider a triangulation in which all simplices are equilateral with fixed side length. The discrete Laplacian is block diagonal in the p-forms. The 0-simplex block matrix follows from eq. 3.7 and takes the form where C ij = 1 if the vertex j is neighbor to vertex i and q i is the total number of neighbor vertices. The first zero mode solution then corresponds to a vector with constant non-zero 0-form components only Φ p 0 = Aδ p0 (1, 1, . . .) . (4.6) Normalizing this vector to unity yields A = 1/ √ N 0 . The second zero mode is gotten from the D-simplex block of the Laplacian and has a structure similar to that given in eqn 4.5 with the matrix element C ij non-zero for each face separating simplices and q i = (D + 1) since each D-simplex possesses (D + 1) neighbors. The corresponding zero mode is (1, 1, . . .) .
Inserting these results into eqn. 4.4 shows that the correlator in eqn 4.4 becomes Next set x = y and z = w and sum over x and z to see that Using cluster decomposition the four point function factorizes into the product of two separate 2 point functions for large |x − z|. But since there is no dependence on these coordinates this implies the factorization should hold for all values of the arguments. Thus the Z 4 should spontaneously break to Z 2 corresponding to the non-zero bilinear condensates These lattice expressions can be expressed in physical units by writing N 0 and N D in terms of the physical volume V and the cut-off a. For example the D-simplex condensate is given by 1 a D−1 V /a D −1 . Note that the magnitude of this condensate depends only on the volume of the space and the Euler number and is independent of the triangulation or continuum metric.

Numerical results
We can illustrate this effect using results obtained from an ensemble of random triangulations of the two sphere corresponding to the partition function where the action This action contains both a bare cosmological constant κ 2 and a coupling β that controls the local connectivity of the triangulation. Again we should emphasize that both the anomaly and  condensates are completely independent of the triangulation depending only on the number of 2-simplices which is held fixed in each ensemble (by tuning κ 2 ).
We show in fig. 1 both the 0 and 2-form bilinear condensates as a function of the bare mass m (the 1-form condensate vanishes as expected for m → 0). We extract the anomaly driven contribution by computing m < ω p ω p > and see that the observed value matches very well with theoretical expectations. The volume dependence of these condensates at β = 0 and m = 0.001 is shown in fig. 2 and shows the expected 1/V behavior.
We can further explore the topological character of these condensates by varying β, which changes the types of triangulation that dominate Z but does not change the global topology. Table. 1

Conclusions
In this paper we have discussed an anomaly that arises in a theory of Kähler-Dirac fermions formulated on a compact manifold. While anomalies are usually thought of as arising only in continuum theories, in this case the anomaly survives intact in a discretized version of the theory. The reason is simple; the anomaly is sensitive only to the topology of the background space which can be captured exactly in the lattice theory as a result of the precise correspondence between the theory of differential forms and the theory of co-chains -the subject of homology theory. As in Fujikawa's derivation of the chiral anomaly, this topological anomaly arises from the non-invariance of the fermion measure under a particular U (1) symmetry specific to Kähler-Dirac fermions. We have explained in some detail this correspondence and how the discrete Kähler-Dirac operator is constructed on a triangulation. In the context of dynamical triangulation models of quantum gravity the latter operator can potentially play a crucial role in determining both the phase structure of the system and the properties of fermions propagating on such general discrete spaces. We hope to examine some of these issues in upcoming work [10]. Here, however, we would like to stress that the anomaly, being topological in character, is completely agnostic about the properties of background space or the nature of any gravitational fluctuations -the effect we describe will be true for a wide range of gravitationally coupled theories involving Kähler-Dirac fermions in regimes described by both classical and quantum gravity both on and off the lattice. It is interesting to ask whether these results could have any implication for cosmology. Since Kähler-Dirac fermions behave locally like multiple copies of Dirac fermions in regions where the curvature is small they are not necessarily ruled out by experiment. Of course there are also global differences between Kähler-Dirac fermions and Dirac fermions. Unlike Dirac fermions, Kähler-Dirac fermions typically possess zero modes even on positive curvature manifolds which is the origin of the anomaly discussed in this paper 3 The physical manifestation of this anomaly is the existence of a fermion condensate. The magnitude of the latter is set by the Planck scale times the volume of the Universe in Planckian units. Presumably such a condensate would feed into the dynamics of the Universe at early times and might have implications for cosmology.