Next Article in Journal
Explaining the ‘Outliers’ Track in Black Hole X-ray Binaries with a BZ-Jet and Inner-Disk Coupling
Next Article in Special Issue
Exploring Gauge Theories with Adjoint Matter on the Lattice
Previous Article in Journal
Numerical Simulation of Ionospheric Disturbances Due to Rocket Plume and Its Influence on HF Radio Waves Propagation
Previous Article in Special Issue
Lattice Computations for Beyond Standard Model Physics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Topological Gauge Actions on the Lattice as Overlap Fermion Determinants

by
Nikhil Karthik
1,2,* and
Rajamani Narayanan
3,*
1
Department of Physics, College of William & Mary, Williamsburg, VA 23185, USA
2
Thomas Jefferson National Accelerator Facility, Newport News, VA 23606, USA
3
Department of Physics, Florida International University, Miami, FL 33199, USA
*
Authors to whom correspondence should be addressed.
Submission received: 10 April 2022 / Revised: 31 May 2022 / Accepted: 10 June 2022 / Published: 17 June 2022

Abstract

:
Overlap fermion on the lattice has been shown to properly reproduce topological aspects of gauge fields. In this paper, we review the derivation of Overlap fermion formalism in a torus of three space-time dimensions. Using the formalism, we show how to use the Overlap fermion determinants in the massless and infinite mass limits to construct different continuum topological gauge actions, such as the level-k Chern–Simons action, “half-CS” term and the mixed Chern–Simons (BF) coupling, in a gauge-invariant lattice UV regulated manner. Taking special Abelian and non-Abelian background fields, we demonstrate numerically how the lattice formalism beautifully reproduces the continuum expectations, such as the flow of action under large gauge transformations.

1. Introduction

The gauge theories in three space-time dimensions admit a parity-odd Chern–Simons (CS) topological gauge action in addition to the parity-even Maxwell gauge action. The Maxwell theory can be nonperturbatively regulated via the lattice discretization of space-time and by using the local plaquette gauge action. The CS theories are not so straightforward to regulate on the lattice, mainly due to the fact that the CS action is only gauge-invariant up to integer winding under nontrivial gauge transformations (e.g., [1]) and it is not possible to realize such a term simply as a local Wilson loop gauge action. Vigorous research work is being conducted on CS theories coupled to matter content and certain infrared duality relations [2,3,4] have been conjectured to exist at critical points separating different topological phases. Therefore, the question of how to study such theories numerically on the lattice is important. The aim of this paper is to elucidate how to introduce topological gauge actions, such as the Chern–Simons action, on the lattice in a completely gauge-invariant manner by identifying such actions as the induced gauge actions of lattice fermions.
Let us first consider gauge theories in even dimensions to see how gauge field topology is realized using lattice fermions. The space of Euclidean continuum gauge fields, A = A μ ( x ) d x μ , in even dimensional space, D = 2 m , usually has infinitely many disconnected pieces and each piece has an associated topological number. This is well known and a chapter or more is attributed to this topic in all modern books on quantum field theory; we find it useful to refer to the lecture notes by Bilal [5] which has a complete self-contained description and has citations to other relevant lecture notes and books. The topological number is given by
Q = 1 m ! ( 2 π ) m P m ( F ) ; P m ( F ) = Tr F m ;
where F = d A + i A A is the Euclidean field strength associated with A μ ( x ) and F m = F F F . As such not all gauge fields can be connected to the trivial one, A μ ( x ) = 0 . One way to nonperturbatively regularize a gauge theory is using lattice, where one introduces gauge fields via gauge-links that connect neighboring lattice sites. Link variables belonging to the Lie group defined by the path ordered product of the Lie group elements,
U μ ( x ) = P x x + μ ^ e i A μ ( y ) d y e i A μ lat ( x )
along the path connecting x and x + μ ^ (we have set the lattice spacing to unity and x takes on integer values) are lattice gauge fields. Naively, U μ q ( x ) = e i q A μ lat ( x ) for some real valued parameter q, continuously connects any gauge field configuration on the lattice to the trivial one, U μ ( x ) = 1 , by sliding the value of q from 0 to 1 seemingly without encountering any singular behavior in gauge-links or the plaquettes at any x during the process. Notwithstanding the apparent lack of discontinuity on the lattice between any two gauge-fields that could otherwise be topologically distinct from each other in the continuum, an assignment of a topological integer to every gauge field configuration is still possible. A straightforward approach is to invoke the Atiyah–Singer index theorem [6] and use fermions to match Q with the index of a lattice Dirac operator. For every lattice gauge field background in even dimensions and the associated massive Hermitian Wilson–Dirac operator, H w ( U ; m w ) , the index is the difference between the total number of negative eigenvalues of H w ( U ; ± m w ) [7]. If the index associated with a particular, U μ ( x ) = e i A μ ( x ) is not zero, we will see an eigenvalue of H w ( U q ; m w ) cross zero as one smoothly changes q [ 0 , 1 ] in U μ q ( x ) = e i q A μ ( x ) . Therefore, there is one value of q where the ground state of the many body operator
H w ( U ; m w ) = a H w ( U ; m w ) a ,
for a D + 1 dimensional auxiliary fermionic system, with a , a being canonical fermion creation and annihilation operators, is doubly degenerate. As is also well known, chiral gauge anomalies in even dimensions are closely related to the topological index [5] and this can also be understood in terms of the ground state, | 0 ; U ; m w , of H w ( U ; m w ) as explained in [8]. Having defined the one form,
d | 0 ; U ; m w = μ , x A μ ( x ) | 0 ; U ; m w d A μ ( x )
it is shown in [8] that
d j cons j cov = Tr P d P d P ; P = | 0 ; U ; m w 0 ; U ; m w |
is a well defined function of the lattice gauge field background and j cons and j cov are the consistent and covariant currents. The problem of anomaly cancellation can be studied using Equation (5) and the need to fine tune the lattice Wilson–Dirac operator is discussed in [8]. The above discussion on the ability of massless overlap fermion to detect and classify topologically distinct gauge sectors on the lattice is well-known. In this paper, we review the aspects of overlap fermions in odd-dimensions, especially in 2 + 1 dimensions, and how the parity anomaly of overlap fermions can be used to introduce topological gauge actions that are characteristic of odd-dimensional gauge theories.
Chiral anomaly inducing topological index in even dimensions and parity anomaly inducing the Chern–Simons action in odd-dimensions are locally related as [5]
P m ( F ) = d Q 2 m 1 ( A , F )
where Q 2 m 1 ( A , F ) is the Chern–Simons form in one dimension lower, namely, D = 2 m 1 . Setting a one-parameter family of gauge fields equal to A t = t A , and noting that F t = t d A + i t 2 A A ,
Q 2 m 1 ( A , F ) = m 0 1 d t Tr ( A F t m 1 ) ; m > 1 .
Focusing on m = 2 , we have
Q 3 ( A , F ) = Tr A d A + i 2 3 A A A .
Similar to our discussion on the challenge in defining the topological index simply as a local operator constructed out of local Wilson-loop operators on the lattice, it is not simple to define the above Chern–Simons form as a local gauge-link-based operator and be able to satisfy invariance under large gauge transformations of the type we will discuss later in this paper. Solution to this problem again is to introduce the Chern–Simons action using the fermions on the lattice; concretely, through the parity-odd part of the induced gauge action from overlap fermions. An early study in Ref. [9] showed that the Abelian parity anomaly is reproduced using lattice perturbation theory with a single-flavor of two-component Wilson fermion with non-zero mass at lattice UV scales [9]. The important point we stress in this paper is that the massive two-component Wilson Dirac operator X on any background field A latt , immediately leads to a gauge covariant unitary operator [10], V,
V X X 1 / 2 X ,
and the gauge-invariant phase of det ( V ) is parity-odd and becomes the lattice realization of the Chern–Simons action for any gauge field background [10,11,12,13]. The unitary operator V is nothing but the overlap operator of a two-component fermion of mass of inverse lattice spacing. The phase within lattice regularization has been extensively analyzed in [14] for various Abelian backgrounds. In addition to the Chern–Simons action, the recent literature on fractional quantum Hall states rely heavily on parity-anomalous two-component exactly massless Dirac fermions that leads to the so-called “half-Chern–Simons” term. Subtleties arise when discussing half the Chern–Simons action while maintaining gauge invariance [2,15,16]. We also show how the construction of the unitary lattice operator V also immediately leads to the generalization of the Chern–Simons term to include the BF terms such as B F = B d A .
In order to keep this paper as self-contained as possible, we first review the derivation and the salient features of overlap fermions in three dimensions in Section 2. In Section 3, we focus on the variation of overlap fermion determinant as fermion mass is varied from 1 / a to massless limit; the point of this discussion is to show that the infinite mass and zero fermion mass limits indeed correctly reproduce the Chern–Simons and “half-Chern–Simons” terms correctly in the continuum limit and independent of any lattice UV regulator parameters, such as the mass term in the Wilson fermion kernel. More interestingly, in Section 4, we take specific Abelian backgrounds with non-trivial topology on 2d spatial planes and show how the flow from infinite mass to zero mass limit preserves gauge invariance. For this, we follow the discussion in [10]. In Section 5, we take a non-Abelian background to discuss how the A A A part of CS term present for non-Abelian case is correctly reproduced. After the discussion of the Chern–Simons terms, in Section 6, we focus on straight-forward extensions of overlap formalism to implement mixed Chern–Simons terms that couple two different gauge field backgrounds, and as a consequence, provide dictionary between some of the recently proposed fermion-boson dualities in the continuum to those on the lattice.

2. Overlap Formalism in Three Dimensions

This section follows [17] very closely and we repeat the derivation while keeping a phase ambiguity intact till the very end. Despite this paper being about nonperturbative regularization of topological field theories, the lattice formalism is strictly presented on toroidal S × S × S manifold tessellated into uniform cubes of volume a 3 , with a being the lattice spacing. The naïve massless Dirac operator on a three dimensional lattice (we will set the lattice spacing to unity) is given by
D = 1 2 μ = 1 3 σ μ T μ T μ ; D = D ;
where σ μ are Pauli matrices satisfying σ μ σ ν = δ μ ν + i ϵ μ ν λ σ λ , and the action of translation operator ( T μ ψ ) ( x ) = U μ ( x ) ψ ( x + μ ^ ) up to lattice periodicity. Under parity ( x μ x μ ),
T μ T μ D D .
and under a gauge transformation G ,
T μ G T μ G ; ( G ϕ ) ( x ) = g ( x ) ϕ ( x ) ; GG = 1 ,
which implies
D G D G .
The naïve massless Dirac operator has a two fold degeneracy in all gauge field backgrounds. Furthermore, for every eigenvalue there is one with the opposite sign. To see these two features, we observe that the anti-Hermitian operator only couples odd lattice sites with even lattice sites. The eigenvalues come in ± i λ pairs and the fermion determinant is real and positive in all gauge backgrounds and there is no parity anomaly. In order to realize a single flavor two-component massive Dirac fermion without any doublers in the overlap formalism [7], we define two Hamiltonians that act on four component spinors:
H = 1 0 0 1 ; H + = B D D B ,
where 1 denotes an identity matrix of the same size as D. We have added the Wilson term,
B = 1 2 μ = 1 3 2 T μ T μ m w ; B = B ,
with a Wilson mass parameter 0 < m w < 2 and B B under parity. Under a gauge transformation
B G B G ; H + G H + G .
Define the many body Hamiltonians by
H ± = a b H ± a b
with a , b and a , b being canonical creation and annihilation operators for fermions. With | 0 ± denoting the ground states of H ± , the generating functional for a single two-component overlap fermion with a mass, m, is
Z ( η , η ¯ ) = 0 | exp η ¯ b + a η + m a b | 0 + = d ξ ¯ d ξ e ξ ¯ ξ 0 | exp χ ¯ b + a χ | 0 + ; χ ¯ = η ¯ + m ξ ¯ ; χ = η + ξ ,
where η , η ¯ , ξ , ξ ¯ are Grassmann variables.
The problem of diagonalizing H + in three dimensions is simplified by going to a new basis. Let
Σ = 1 2 1 1 1 1 ; Σ = Σ ; Σ 2 = 1 .
The rotated Hamiltonian is
H + = Σ H + Σ = 0 B D B + D 0 0 X X 0 .
We can write
X = L Λ R ; Λ i j = λ i δ i j ; λ i > 0 ; R R = L L = 1 .
We define the unitary operator as
V = L R = 1 X X X
and does not suffer from the phase ambiguity present in R (L is fixed once R is fixed). Under parity,
X X ; V V
and under a gauge transformation
X G X G ; V G V G .
Let us make the dependence of V on U explicit and derive the relation under charge conjugation ( U U * ):
V ( U * ) = σ 2 V t ( U ) σ 2 .
We first note that
T μ * ( U ) = T μ ( U * ) B * ( U ) = σ 2 B ( U * ) σ 2 ; D * ( U ) = σ 2 D ( U * ) σ 2 .
From this we obtain
X ( U * ) = σ 2 X t ( U ) σ 2 , X ( U * ) X ( U * ) = σ 2 X * ( U ) X t ( U ) σ 2 ,
and our relation, Equation (25), follows.
We can diagonalize H + as
H + = U Λ 0 0 Λ U ; U = 1 2 R + L R L R L R + L ; UU = 1 .
We define new sets of canonical creation and annihilation operators by
c = a R + L 2 + b R L 2 ; d = a R L 2 + b R + L 2 , c = R + L 2 a + R L 2 b ; d = R L 2 a + R + L 2 b ,
and we can write
H + = c Λ c + d Λ d .
The ground states, | 0 ± , are obtained by filling all the states corresponding to c and a , respectively. Therefore, we have
c | 0 + = 0 ; d | 0 + = 0 ; 0 | a = 0 ; 0 | b = 0 .
Using Equation (29), we can write
a = c 2 R + L b A ; b = 2 R + L d + G a ,
where
G = 1 V 1 + V .
Using the above equations, we can write
χ ¯ b + a χ = Q + + Q ,
where
Q + = c 2 R + L χ + χ ¯ 2 R + L d ; Q = χ ¯ G a b G χ .
Since
Q + = a χ + b G χ χ ¯ G a + χ ¯ b
it follows that
[ Q + , Q ] = 2 χ ¯ G χ .
Therefore, we have
Z ( η , η ¯ ) = d ξ ¯ d ξ e ξ ¯ ξ 0 | exp [ Q + + Q ] | 0 + = d ξ ¯ d ξ e ξ ¯ ξ exp 1 2 [ Q + , Q ] 0 | e Q e Q + | 0 + = d ξ ¯ d ξ e ξ ¯ ξ exp 1 2 [ Q + , Q ] 0 | 0 + = d ξ ¯ d ξ exp [ ξ ¯ ξ ( η ¯ + m ξ ¯ ) G ( η + ξ ) ] det R + L 2 = exp η ¯ G 1 + m G η det ( 1 + m G ) det 1 + V 2 det R = exp η ¯ 1 V 1 + m + ( 1 m ) V η det ( 1 + m ) + ( 1 m ) V 2 det R .
The fermion mass is in the range m [ 1 , 1 ] and the fermion determinant is gauge invariant. There is a phase ambiguity present in the fermion determinant due to det R and the fermion determinant at m = 1 is det R . The choice of fixing this phase is tied to the choice of preserving parity symmetry at m = 0 at the cost of introducing gauge anomaly, and the choice of preserving gauge invariance at the cost of losing parity symmetry at m = 0 . For the latter option, the choice of det ( R ) = 1 fixes the phase of infinite mass, m = 1 , fermion and preserves gauge invariance for all values of m. In this paper, we will set det R = 1 from here on.

3. Introducing Chern–Simons and Half-Chern–Simons Terms on the Lattice

With the choice of phase as explained in the last section, the fermion determinant becomes
Z ( m ) = det ( 1 + m ) + ( 1 m ) V 2 ,
and it satisfies
Z ( m ) Z * ( m ) = det V .
This is the parity anomaly which has the built in feature that if we write
Z ( m ) = | Z ( m ) | e i Φ ( m )
then
e i [ Φ ( m ) + Φ ( m ) ] = e 2 i Φ ( 0 ) = det V ,
and Φ ( 0 ) is usually written using the η -invariant as π η 2 (Note, the propagator satisfies G ( m ) = G ( m ) and preserve parity. Thus the anomaly is in the fermion induced gauge measure.) With this lattice formalism, we have all the required ingredients for constructing Chern–Simons theories on lattice by the identification of the parity-odd phase of det ( V ) with level-1 Chern–Simons action. As the simplest case, we can introduce a level-k Chern–Simons action as
e i S CS ( k ) = det ( V ) k .
In the massless limit, it is easy to see that the phase of det ( 1 + V ) is half of det ( V ) up to ± 1 . We can introduce the so-called U k + N f / 2 ( 1 ) “half-Chern–Simons” theories on the lattice as
e i S CS ( k + N f / 2 ) = det 1 + V 2 N f det ( V ) k .
First, how do we know that det ( V ) is the same as Chern–Simons term? In the study in [9], it is analytically shown that the phase of det ( X ) in the massive Wilson fermion case is the same as Chern–Simons term. The pure phase det ( V ) in the case of overlap fermions is the same as the phase of det ( X ) , and, hence, we can borrow their results for overlap fermions. In the subsequent two sections, we will also take an empirical approach and show that for cases of Abelian and non-Abelian background fields where Chern–Simons term can be exactly be worked out, the phase det ( V ) indeed approaches the expectations in the continuum limit. Second, how did we manage to introduce “half-Chern–Simons” term in an evidently gauge-invariant manner? Using a non-trivial Abelian background in the next section, we demonstrate this through the flow of the phase of massless overlap fermion determinant as a function of the Wilson loop, e i 2 π h 3 for h 3 [ 0 , 1 ] , and show that at specific h 3 where there is a discontinuity in the phase at m = 0 , the determinant also vanishes.

4. Fermion Determinant in an Abelian Background with Uniform Magnetic Flux and Non-Trivial Temporal Wilson Loop

We now analyze the complex fermion determinant of a two-component three dimensional fermion in a well known Abelian background of interest both from the view point of showing subtle properties under gauge invariance and also from its relevance in condensed matter physics [2]. The gauge field background on a continuum 3 torus is
A 1 = 2 π Q x 2 2 ; A 2 = 0 ; A 3 = 2 π h 3 .
Since A 1 ( x 1 , , x 3 ) has to be gauge equivalent to A 1 ( x 1 , 0 , x 3 ) , Q has to be an integer. In addition, gauge invariance sets all h 3 + n to be equivalent for any integer n. The evaluation of the Chern–Simons action for this background in Equation (45) is tricky [16] and yields
S cs = 2 π h 3 Q .
Since F = 2 π Q 2 d x 1 d x 2 , for this background, Q is the topological charge in all two-dimensional slices at a fixed x 3 and the deformation of A to t A has to connect two dimensional gauge fields in disconnected spaces. With a lattice regularization, t Q , as t goes from 0 1 , will result in Q levels of the two dimensional Wilson–Dirac operator crossing zero [14] and the phase within lattice regularization properly reproduces the gauge invariant Chern–Simons action [14]. Overlap fermions can be used to study the complex fermion determinant strictly in the massless limit with the lattice regularization in place and we will show that the massless fermion determinant has a zero in the path connecting h 3 and h 3 + 1 for a fixed Q enabling it to correctly reproduce (1) a smooth function of h 3 , (2) that is gauge invariant under h 3 h 3 + 1 , (3) equal to half of S cs in Equation (46) at all values of h 3 , and (4) has a jump in the phase at the location of the zero of the fermion determinant.
We can implement the above Abelian background on the lattice by using the gauge-links as
U 1 ( x ) = 1 x 1 L 1 e i 2 π Q L x 2 x 1 = L 1 ; U 2 ( x ) = e i 2 π Q L 2 x 1 ; U 3 ( x ) = e i 2 π h 3 L ;
on a three dimensional periodic lattice defined by the points x 1 , x 2 , x 3 [ 0 , L 1 ] and
U μ ( x + L ν ^ ) = U μ ( x ) ; μ , ν = 1 , 2 , 3 .
Only the plaquettes in the ( 1 , 2 ) plane have a non-zero flux and they are given by
U 12 ( x ) = U 1 ( x ) U 2 ( x + 1 ^ ) U 1 * ( x + 2 ^ ) U 2 * ( x ) = exp i 2 π Q L 2 x ( L 1 , L 1 , x 3 ) exp i 2 π Q L 2 i 2 π Q x ( L 1 , L 1 , x 3 ) .
We note that the flux is not uniform and singular in the continuum limit if Q is not an integer. Therefore, we will set Q to be integers.
Since the gauge field background does not depend on x 3 , one can go to momentum space in this direction. We will assume fermions obey antiperiodic boundary conditions in this direction. Setting these momenta to be 2 π k 3 L π L , k 3 [ 0 , L 1 ] , the operators B and D reduce to
B ( k 3 ) = 1 2 μ = 1 2 ( 2 T μ T μ ) + 2 sin 2 π h 3 1 2 + k 3 L m w ; D ( k 3 ) = 1 2 μ = 1 2 σ μ ( T μ T μ ) + i σ 3 sin 2 π h 3 1 2 + k 3 L .
with the gauge fields in the (1-2) plane being U 1 ( x ) and U 2 ( x ) . Let us denote the fermion determinant by Z L ( h 3 , Q , m ; m w ) on the L 3 periodic lattice in this background and note that
Z L 0 ( h 3 , Q , m ; m w ) = Z L 0 ( h 3 + 1 , Q , m ; m w ) .
We define
Z L 0 ( h 3 , Q , m ; m w ) = Z L ( h 3 , Q , m ; m w ) Z L ( 0 , Q , m ; m w ) = Z L 0 ( h 3 , Q , m ; m w ) exp i Φ L 0 ( h 3 , Q , m ; m w ) .
and
Z L ( 0 , Q , m ; m w ) Z L ( 0 , 0 , m ; m w ) = exp F L ( Q , m ; m w )
as the determinant with reference to h 3 = 0 and the determinant at h 3 = 0 with respect to the free determinant, respectively and F L is a real function. The key properties of the overlap fermion determinant are shown in the figures from Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5. Let us start with the top panel of Figure 1 which focuses on the Chern–Simons action, namely, Φ L 0 ( h 3 , Q , 1 ; m w ) . We have shown the results only for m w = 1 but the L limit is independent of m w and we should find
lim L Φ L 0 ( h 3 , Q , 1 ; m w ) = 2 π h 3 Q .
The top panel clearly shows that the correct limit is approached for Q = 1 as L ( L = 8 and L = 16 fall on top of each other) and the dependence on Q is also as expected and the overlap fermion correctly reproduces the first subtle properly and this is an obvious consequence of the same result with Wilson fermions seen in [14]. We move on to behavior of the phase for the massless fermions in the bottom panel of Figure 1. We should find
lim L 2 Φ L 0 ( h 3 , Q , 0 ; m w ) = lim L Φ L 0 ( h 3 , Q , 1 ; m w )
and
lim L Φ L 0 ( h 3 , Q , 0 ; m w ) = lim L Φ L 0 ( h 3 + 1 , Q , 0 ; m w ) .
This necessitates a jump in the phase when the flux quantum, Q, takes on odd values. First of all, we see that the phase has a limit when L as seen by comparing the behavior for L = 8 , m w = 1 and L = 16 , m w = 1 . Furthermore, the results for L = 16 , m w = 0.5 and L = 16 , m w = 1.5 are indistinguishable from L = 16 , m w = 1 showing the independence on the regulator parameter, m w , as L . Finally, we see that the phase shows a jump of π at h 3 = 1 2 for Q = 1 and Q = 3 .
The plot of the full determinant, Z L 0 ( h 3 , Q , m ; 1 ) , is shown for Q = 1 , Q = 2 and Q = 3 in Figure 2, Figure 3 and Figure 4 respectively. In these plots, h 3 [ 0 , 1 ] , and the motion along the closed curve is clockwise starting from the normalized value of Z L 0 ( 0 , Q , m ; 1 ) = 1 . When m = 1 , the closed curves are unit circles that wind Q times and this is shown for reference in all three plots. We set m L to be a constant when m ( 0 , 1 ) to maintain a constant physical mass. On the one hand, we see that Z L 0 ( h 3 , Q , m ; 1 ) winds around Q times for m < 0 and its magnitude changes with h 3 . On the other hand, we see that the phase of Z L 0 ( h 3 , Q , m ; 1 ) , reaches a maximum and minimum value in the range π 2 , π 2 for m > 0 and its magnitude changes with h 3 . With the behavior in place for m < 0 and m > 0 , we see that Z L 0 1 2 , Q , 0 ; 1 is zero and enables a jump in the phase for odd values of Q with it being a smooth function of h 3 . Finally, we show the results for F L ( Q , m ; m w ) in Figure 5. It remains finite as L , which is the continuum limit of the background field, and independent of the regulator parameter, m w .

5. Fermion Determinant in a Non-Abelian SU ( 2 ) Background with Non-Zero Tr ( A A A )

The second background we will consider is a constant s u ( 2 ) background on a 3 torus given by
A 1 = 2 π q 1 t 1 ; A 2 = 2 π q 2 t 2 ; A 3 = 2 π q 3 t 3 ;
where t μ are the s u ( 2 ) generators in color space given by Pauli matrices, normalized such that t μ t ν = δ μ ν + i ϵ μ ν ρ t ρ . In this case q i 0 , 2 are all gauge inequivalent and the Chern–Simons action reduces to
S cs = 16 π 2 q 1 q 2 q 3 .
Contrary to the Abelian background the phase of the massless fermion determinant is simply given by 8 π 2 q 1 q 2 q 3 and we will show this to be the case. Defining Z L ( q 1 , q 2 , q 3 , m ; m w ) as the lattice regulated overlap fermion determinant on a L 3 periodic lattice with m being the fermion mass and m w being another regulator parameter, we will show that both
lim L Z L ( q 1 , q 2 , q 3 , m ; m w ) Z L ( q 2 , q 2 , 0 , m ; m w ) ; lim L Z L ( q 1 , q 2 , 0 , m ; m w ) Z L ( 0 , 0 , 0 , m ; m w )
are both finite and independent of the regulator m w . This constant SU(2) background can be introduced on the lattice as the link variables
U μ = e i 2 π q μ t μ L .
We will consider this background on a three dimensional periodic lattice defined by the points n 1 , n 2 , n 3 [ 0 , L 1 ] and
U μ ( n + L ν ^ ) = U μ ( n ) ; μ , ν = 1 , 2 , 3 .
All values of q μ that remain finite as L are gauge inequivalent.
One can go to momentum space in all three directions and write
B i a , j b ( k ) = b δ i j δ a b + μ s μ ( q ) s μ a ( k ) δ i j t μ a b ; D i a , j b ( k ) = i μ c μ ( q ) s μ a ( k ) σ μ i j δ a b + i μ s μ ( q ) c μ a ( k ) σ μ i j t μ a b ,
where
b = 3 m w μ c μ ( q ) c μ a ( k ) ; c μ ( p ) = cos 2 π p μ L s μ ( p ) = sin 2 π p μ L ; c μ a ( k ) = cos 2 π k μ L μ = 1 , 2 cos 2 π k μ + π L μ = 3 ; s μ a ( k ) = sin 2 π k μ L μ = 1 , 2 sin 2 π k μ + π L μ = 3 .
We have assumed anti-periodic boundary conditions for fermions in the μ = 3 direction. The matrix X ( k ) is given by
X ( k ) = α 1 α 3 α 4 α 5 α 3 * α 2 α 6 α 4 α 4 * α 6 α 1 * α 3 α 5 α 4 * α 3 * α 2 *
where
α 1 = b + s 3 ( q ) s 3 a ( k ) + i c 3 ( q ) s 3 a ( k ) + i s 3 ( q ) c 3 a ( k ) α 2 = b s 3 ( q ) s 3 a ( k ) + i c 3 ( q ) s 3 a ( k ) i s 3 ( q ) c 3 a ( k ) α 3 = s 1 ( q ) s 1 a ( k ) i s 2 ( q ) s 2 a ( k ) α 4 = i c 1 ( q ) s 1 a ( k ) + c 2 ( q ) s 2 a ( k ) α 5 = i s 1 ( q ) c 1 a ( k ) i s 2 ( q ) c 2 a ( k ) α 6 = i s 1 ( q ) c 1 a ( k ) + i s 2 ( q ) c 2 a ( k ) .
Let us denote the fermion determinant by Z L ( q 3 , m ; m w ) on the L 3 periodic lattice in this background and define
Z L 0 ( q 1 , q 2 , q 3 , m ; m w ) = Z L ( q 1 , q 2 , q 3 , m ; m w ) Z L ( q 1 , q 2 , 0 , m ; m w ) = Z L 0 ( q 1 , q 2 , q 3 , m ; m w ) exp i Φ L 0 ( q 1 , q 2 , q 3 , m ; m w ) .
and
Z L ( q 1 , q 2 , 0 , m ; m w ) Z L ( 0 , 0 , 0 , m ; m w ) = exp F L ( q 1 , q 2 , m , ; m w )
as the determinant with reference to q 3 = 0 and the determinant at q 3 = 0 with respect to the free determinant, respectively and F L is a real function. We will set q 1 = 1 4 and q 2 = 1 2 π and vary q 3 . The nonabelian Chern–Simons action given in Equation (58) reduces to S cs = 2 π q 3 and we show the phase of the overlap fermion correctly reproduces this result as L in the top panel of Figure 6. Since all q 3 are gauge inequivalent, we should find
lim L Φ L 0 1 4 , 1 2 π , q 3 , 0 ; m w = π q 3
and we should also find
lim L Φ L 0 1 4 , 1 2 π , q 3 , m L ; m w + Φ L 0 1 4 , 1 2 π , q 3 , m L ; m w = 2 π q 3
Both these features are correctly reproduced in the top panel of Figure 6. Since all q 3 are gauge inequivalent, we see that the phase at q 3 = 1 and q 3 = 2 only approaches 2 π q 3 as L . Note that unlike the Abelian case, the determinant winds around the origin for all values of fermion mass and the fermion determinant remains non-zero for all values of q 3 . This is made clear through a plot of Z L 0 1 4 , 1 2 π , q 3 , m L ; 1 in Figure 7.
We note a curious observation in this particular background. The fermion determinant for massless fermions becomes very small for certain values of q 3 and it has zeros even at finite L that remains stable as L as seen in Figure 7. For our choice of q 1 and q 2 , we find zeros a pair of zeros at q 3 = 0.508 and q 3 = 0.646 and another pair at q 3 = 1.502 and q 3 = 1.556 that remain stable across L. In spite of the fact that all q 3 are gauge inequivalent, we see non-trivial behavior seen in the complex determinant for massless fermions in this particular background. Finally, similar to the Abelian background, we found the results for F L 1 4 , 1 2 π , m ; m w to be finite as L and independent of the regulator parameter, m w .

6. Mixed Chern–Simons (BF) Action and Dualities

Let V A denote the dependence of the unitary operator in Equation (22) on the Abelian gauge field background, A. A mixed Chern–Simons (BF) term can be written as
det V A V B V A + B = det V A V B V A B e i 2 π d 3 x ϵ μ ν λ A μ ν B λ .
One can formally verify the identity by inserting the naïve expressions for CS that are only valid for perturbative fields. The path integrals are defined over all gauge fields and a suitable measure such as a standard Maxwell action for gauge fields is needed to verify the integrals non-perturbatively. Therefore, the last step is essentially a mnemonic and it suggests relations of the form
[ d A ] det V A B V A C = δ ( B C ) , [ d A ] det V A B V A C V A D V A E = V B V C V D V E δ ( B + C D E ) ,
using a gauge action for A that is implicit and allows one to take a continuum limit (as pointed out explicitly in [3]).
Dualities among various three dimensional theories start with the conjecture [2] that a theory with one massless two component fermion coupled to a dynamical gauge field A and a classical background C defined by
Z ( C ) = [ d A ] e S g ( A ) det 1 + V A 2 det V A V A C = [ d A ] e S g ( A ) det 1 + V A 2 det V A C ,
is parity even and dual to a theory at the Wilson–Fisher fixed point. An explicit computation shows that
Z * ( C ) = [ d A ] e S g ( A ) det 1 + V A 2 det V A C = [ d A ] e S g ( A ) det 1 + V A 2 det V A C Z * ( C ) = [ d A ] e S g ( A ) det 1 + V A 2 det V A 2 V C 2 V A C .
If Z * ( C ) = Z ( C ) , we arrive at a non-trivial relation
det V A = det V C 2
where the expectation value is with respect to the measure in Equation (71) and the lattice regularization can be used to verify this relation. In fact, if we use
V 2 A = V A 4
which has been verified in the continuum limit when a measure for the gauge field is included [18], we see that if we assume that the dynamical fermion has a charge of 2 units,
Z 2 ( C ) = [ d A ] e S g ( A ) det 1 + V 2 A 2 det V A 2 det V A V C V A C
then Z 2 * ( C ) = Z 2 ( C ) is trivially satisfied.
Regularized versions of the various duality relations discussed in [2] can be obtained by following the steps found there. We multiply both sides of Equation (71) by det V B + C , promote C to the dynamical field with B being a background field and arrive at a regularized version of a fermion-boson duality,
e S g ( B ) det 1 + V B 2 = [ d C ] Z ( C ) det V B + C ,
after using Equation (70). If we assume Z ( C ) is real we arrive at a regularized version of a boson–boson duality
[ d C ] Z ( C ) det V B V B + C = [ d C ] Z ( C ) det V B + C .
We can multiply both sides of Equation (71) by det V B + C , promote C to a dynamical field with B being a background field and arrive at a regularized version of a fermion-fermion duality
[ d A ] [ d C ] e S g ( A ) det 1 + V A 2 det V A V A C V B + C = e S g ( B ) det 1 + V B 2
and we have used Equations (76) and (77).
A regularized version of a duality involving a fermion with charge of 2 units discussed in [19] can be obtained by setting B = 2 X in Equation (78). In this case, we can multiply both sides by det 2 V X to make the right-hand side even under parity. We also multiply both sides by det V X V Y V X Y to couple it to an external flux and promote X to a dynamical field. Then we have
Z ( Y ) = [ d A ] [ d C ] [ d X ] e S g ( A ) det 1 + V A 2 det V A det V A C V 2 X + C det 2 V X det V X V Y V X Y = [ d X ] e S g ( 2 X ) det 1 + V 2 X 2 det 2 V X det V X V Y V X Y
Defining a change of variable, X = Z + C , in the first integral, we obtain
Z ( Y ) = [ d A ] e S g ( A ) det 1 + V A 2 det V A det V Y [ d C ] [ d Z ] det V A C V 2 Z + C V Z + C V Z + C Y .
The integral over C can be performed using Equation (70) This forces A = Y and we arrive at the regularized version of a fermion-fermion duality
e S g ( Y ) det 1 + V Y 2 [ d Z ] det V Z V Y V Z + Y = [ d X ] e S g ( 2 X ) det 1 + V 2 X 2 det V X V Y V X Y
that connects a fermion with 2 units of charge to a fermion with 1 unit of charge.
We should remark that for the sake of simplicity and to a first degree of approximation, we assumed that the massless fermion limits of the odd-flavored theories considered above occurs at the “bare” fermion mass m = 0 . Unlike the parity-invariant theories with SU N f flavor symmetry with N f being even, where the mass term is protected by the symmetry, there is no such symmetry consideration in odd flavored theories. Thus, it could be possible that one needs to tune the overlap fermion mass m = m c in order to reach criticality, provided there is one. In that case, the above set of equations might have to be modified accordingly with such mass terms, but it is a straightforward exercise.

7. Conclusions

Overlap formalism was developed three decades ago [7] to properly reproduce all salient features of massless fermions in even dimensions. This was extended to odd dimensions in [10] and we showcase the salient features of massless fermions in odd dimensions; particularly, we extended the formalism and spell out the lattice constructions of topological gauge actions that are being investigated currently in the context of TQFTs coupled to fermions, and in the context of infrared dualities. We focused on the overlap fermion determinant and used two examples, one Abelian background and one non-Abelian background. We showed that the overlap fermion determinant correctly reproduces all known properties of the phase of the fermion determinant; especially, we discussed how the lattice regularization manages to implement the half-Chern–Simons term (or half-the-eta-invariant) in a gauge-invariant manner. While it is satisfying that we can nonperturbatively formulate the topological gauge theories on the lattice, an actual numerical study of such theories is not yet practical due to the sign problem and we did not address such issues in this paper.
An interesting possibility of having a lattice regularized Chern–Simons theory is the following. As we noted in this paper, it is important to realize that the identification of overlap det ( V ) with the continuum Chern–Simons action e i S CS is possible only in the continuum limit (as usual, the continuum limit taken at the trivial UV fixed point of the lattice gauge theory). However, as a lattice gauge theory that is away from any critical points, the overlap fermion determinants offer a great way to introduce new parity-odd gauge-invariant gauge-actions. Thus, one could now ask about the phase diagrams of such well defined lattice gauge theories as a function of different lattice couplings. This is an exciting direction to think about in the future.

Author Contributions

Conceptualization, N.K. and R.N.; methodology, N.K. and R.N.; software, N.K. and R.N.; investigation, N.K. and R.N.; writing—original draft preparation, R.N.; writing—review and editing, N.K.; visualization, N.K. and R.N.; funding acquisition, N.K. and R.N. All authors have read and agreed to the published version of the manuscript.

Funding

R.N. acknowledges partial support by the NSF under grant number PHY-1913010. N.K. is supported by Jefferson Science Associates, LLC under U.S. DOE Contract #DE-AC05-06OR23177 and in part by U.S. DOE grant #DE-FG02-04ER41302.

Data Availability Statement

Data available up on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dunne, G.V. Aspects of Chern–Simons theory. In Aspects Topologiques de la Physique en Basse Dimension. Topological Aspects of Low Dimensional Systems; Springer: Berlin/Heidelberg, Germany, 1999. [Google Scholar]
  2. Seiberg, N.; Senthil, T.; Wang, C.; Witten, E. A Duality Web in 2+1 Dimensions and Condensed Matter Physics. Ann. Phys. 2016, 374, 395–433. [Google Scholar] [CrossRef] [Green Version]
  3. Karch, A.; Tong, D. Particle-Vortex Duality from 3d Bosonization. Phys. Rev. 2016, X6, 031043. [Google Scholar] [CrossRef] [Green Version]
  4. Wang, C.; Nahum, A.; Metlitski, M.A.; Xu, C.; Senthil, T. Deconfined quantum critical points: Symmetries and dualities. Phys. Rev. 2017, X7, 031051. [Google Scholar] [CrossRef] [Green Version]
  5. Bilal, A. Lectures on Anomalies. arXiv 2008, arXiv:0802.0634. [Google Scholar]
  6. Atiyah, M.F.; Singer, I.M. The index of elliptic operators on compact manifolds. Bull. Am. Math. Soc. 1969, 69, 422–433. [Google Scholar] [CrossRef] [Green Version]
  7. Narayanan, R.; Neuberger, H. A Construction of lattice chiral gauge theories. Nucl. Phys. 1995, B443, 305–385. [Google Scholar] [CrossRef] [Green Version]
  8. Neuberger, H. Geometrical aspects of chiral anomalies in the overlap. Phys. Rev. 1999, D59, 085006. [Google Scholar] [CrossRef] [Green Version]
  9. Coste, A.; Luscher, M. Parity Anomaly and Fermion Boson Transmutation in Three-dimensional Lattice QED. Nucl. Phys. 1989, B323, 631. [Google Scholar] [CrossRef]
  10. Kikukawa, Y.; Neuberger, H. Overlap in odd dimensions. Nucl. Phys. 1998, B513, 735–757. [Google Scholar] [CrossRef] [Green Version]
  11. Narayanan, R.; Nishimura, J. Parity invariant lattice regularization of three-dimensional gauge fermion system. Nucl. Phys. 1997, B508, 371–387. [Google Scholar] [CrossRef]
  12. Bietenholz, W.; Nishimura, J. Ginsparg-Wilson fermions in odd dimensions. JHEP 2001, 7, 015. [Google Scholar] [CrossRef] [Green Version]
  13. Bietenholz, W.; Nishimura, J.; Sodano, P. Chern–Simons theory on the lattice. Nucl. Phys. B Proc. Suppl. 2003, 119, 935–937. [Google Scholar] [CrossRef] [Green Version]
  14. Karthik, N.; Narayanan, R. Phase of the fermion determinant in QED3 using a gauge invariant lattice regularization. Phys. Rev. 2015, D92, 025003. [Google Scholar]
  15. Alvarez-Gaume, L.; Della Pietra, S.; Moore, G.W. Anomalies and Odd Dimensions. Ann. Phys. 1985, 163, 288. [Google Scholar] [CrossRef]
  16. Witten, E. Three lectures on topological phases of matter. Riv. Nuovo Cim. 2016, 39, 313–370. [Google Scholar]
  17. Karthik, N.; Narayanan, R. Scale-invariance of parity-invariant three-dimensional QED. Phys. Rev. 2016, D94, 065026. [Google Scholar] [CrossRef] [Green Version]
  18. Karthik, N.; Narayanan, R. Parity Anomaly Cancellation in Three-Dimensional QED with a Single Massless Dirac Fermion. Phys. Rev. Lett. 2018, 121, 041602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Cordova, C.; Hsin, P.S.; Seiberg, N. Time-Reversal Symmetry, Anomalies, and Dualities in (2 + 1)d. SciPost Phys. 2018, 5, 6. [Google Scholar] [CrossRef]
Figure 1. The top panel shows the flow of the phase Φ L 0 ( Q , h 3 ) in the infinite mass case, m = 1 , as a function of Wilson-loop variable h 3 [ 0 , 1 ] at Q = 1 , 2 , 3 . The Wilson mass entering the kernel of overlap operator is fixed at m w = 1 . For fixed Q = 1 , the variation with reduction in lattice spacing by increasing L from 4 to 8 is also shown. The bottom panel shows similar flow of the phase of the determinant in the massless case. The variability with respect to the regulator parameter m w and lattice spacing are shown.
Figure 1. The top panel shows the flow of the phase Φ L 0 ( Q , h 3 ) in the infinite mass case, m = 1 , as a function of Wilson-loop variable h 3 [ 0 , 1 ] at Q = 1 , 2 , 3 . The Wilson mass entering the kernel of overlap operator is fixed at m w = 1 . For fixed Q = 1 , the variation with reduction in lattice spacing by increasing L from 4 to 8 is also shown. The bottom panel shows similar flow of the phase of the determinant in the massless case. The variability with respect to the regulator parameter m w and lattice spacing are shown.
Universe 08 00332 g001
Figure 2. The flow of the overlap fermion determinant in the complex plane as a function of h 3 at a fixed Q = 1 on L = 16 lattice. The flows are shown for m = 1 (black), 1 / L (red), 0 (green), 1 / L (blue). The flow starts at ( 1 , 0 ) for h 3 = 0 , goes clockwise and returns back to ( 1 , 0 ) for h 3 = 1 .
Figure 2. The flow of the overlap fermion determinant in the complex plane as a function of h 3 at a fixed Q = 1 on L = 16 lattice. The flows are shown for m = 1 (black), 1 / L (red), 0 (green), 1 / L (blue). The flow starts at ( 1 , 0 ) for h 3 = 0 , goes clockwise and returns back to ( 1 , 0 ) for h 3 = 1 .
Universe 08 00332 g002
Figure 3. The flow of the overlap fermion determinant in the complex plane as a function of h 3 at a fixed Q = 2 on L = 16 lattice. The description is the same as in Figure 2.
Figure 3. The flow of the overlap fermion determinant in the complex plane as a function of h 3 at a fixed Q = 2 on L = 16 lattice. The description is the same as in Figure 2.
Universe 08 00332 g003
Figure 4. The flow of the overlap fermion determinant in the complex plane as a function of h 3 at a fixed Q = 3 on L = 16 lattice. The description is the same as in Figure 2.
Figure 4. The flow of the overlap fermion determinant in the complex plane as a function of h 3 at a fixed Q = 3 on L = 16 lattice. The description is the same as in Figure 2.
Universe 08 00332 g004
Figure 5. The plot demonstrates the existence of the continuum limit of the overlap fermion action in constant flux background at zero h 3 . The continuum extrapolations ( L ) are shown using an expansion in lattice spacing 1 / L . The consistency in the extrapolated values using different regulator parameter m w is seen.
Figure 5. The plot demonstrates the existence of the continuum limit of the overlap fermion action in constant flux background at zero h 3 . The continuum extrapolations ( L ) are shown using an expansion in lattice spacing 1 / L . The consistency in the extrapolated values using different regulator parameter m w is seen.
Universe 08 00332 g005
Figure 6. The figure is similar to Figure 1 showing the flow of the phase of the fermion determinant as a function of SU(2) gauge field magnitude q 3 . The top panel shows the result for infinitely massive fermion, m = 1 with regulator parameter m w = 1 . The convergence of the results at different L towards a continuum result is shown. The bottom panel shows the flow with q 3 at different fermion masses m.
Figure 6. The figure is similar to Figure 1 showing the flow of the phase of the fermion determinant as a function of SU(2) gauge field magnitude q 3 . The top panel shows the result for infinitely massive fermion, m = 1 with regulator parameter m w = 1 . The convergence of the results at different L towards a continuum result is shown. The bottom panel shows the flow with q 3 at different fermion masses m.
Universe 08 00332 g006
Figure 7. The dependence of the magnitude of the fermion determinant on q 3 . The result at zero and non-zero masses are shown. The determinant at zero mass vanishes are certain values of q 3 , though not for any reasoning from invariance under large gauge transformation as seen in the case of Abelian background field studied in this paper.
Figure 7. The dependence of the magnitude of the fermion determinant on q 3 . The result at zero and non-zero masses are shown. The determinant at zero mass vanishes are certain values of q 3 , though not for any reasoning from invariance under large gauge transformation as seen in the case of Abelian background field studied in this paper.
Universe 08 00332 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Karthik, N.; Narayanan, R. Topological Gauge Actions on the Lattice as Overlap Fermion Determinants. Universe 2022, 8, 332. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8060332

AMA Style

Karthik N, Narayanan R. Topological Gauge Actions on the Lattice as Overlap Fermion Determinants. Universe. 2022; 8(6):332. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8060332

Chicago/Turabian Style

Karthik, Nikhil, and Rajamani Narayanan. 2022. "Topological Gauge Actions on the Lattice as Overlap Fermion Determinants" Universe 8, no. 6: 332. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8060332

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop