Non-abelian lattice gauge theory with a topological action

SU(2) gauge theory is investigated with a lattice action which is insensitive to small perturbations of the lattice gauge fields. Bare perturbation theory can not be defined for such actions at all. We compare non-perturbative continuum results with that obtained by the usual Wilson plaquette action. The compared observables span a wide range of interesting phenomena: zero temperature large volume behavior (topological susceptibility), finite temperature phase transition (critical exponents and critical temperature) and also the small volume regime (discrete beta-function or step-scaling function). In the continuum limit perfect agreement is found indicating that universality holds for these topological lattice actions as well.


Introduction
Universality for lattice gauge theory is generally thought to mean that any lattice action can be used in simulations provided (1) it reproduces the desired continuum action in the classical continuum limit, (2) is local, (3) has the correct symmetries, (4) the continuum extrapolation is performed. Due to asymptotic freedom bare perturbation theory then correctly predicts the approach to the continuum.
The possibility that universality in field theory holds more broadly was first addressed in the non-linear O(3) model in 2 dimensions [1]. It was shown that a topological lattice action which is insensitive to small perturbations of the lattice fields, gives the correct results in the continuum. The topological nature of the lattice action means that the classical vacuum is infinitely degenerate and hence bare perturbation theory can not be set up and is inherently meaningless. The only result that may be obtained from a topological action is the fully non-perturbative one and it is apparently the same as the one with the usual lattice action of the O(3) model. Hence it seems that requirement (1) above can be dropped and universality still holds. Similar results were also shown [2,3] to hold for the 2d XY model as well. 1 From the point of view of the path integral a useful way of thinking about the type of topological actions we investigate is the following. Field space is divided into two sets: one, where the action is zero and two, where the action is infinite. Hence configurations from the former enter with equal weights and fluctuate freely while configurations from the latter are forbidden. The only non-trivial information about the field theory is then encoded in the boundary separating the two sets.
We investigate the same phenomenon for non-abelian gauge theories. A topological lattice action can easily be defined analogously to the O(3) model and hence bare perturbation theory is again meaningless. Nevertheless we show that for SU (2) pure gauge theory the non-perturbative continuum results obtained with the topological action agree with that of the usual Wilson plaquette action which we know is in the correct universality class of Yang-Mills theory. The compared observables are sensitive to a wide range of interesting physics. At T = 0 we calculate the topological susceptibility and set the scale t 0 by the gradient flow. At T = T c we compare the dimensionless ratio T c √ 8t 0 and also the critical exponent 1/ν from the Binder cumulant of the Polyakov loop. In small physical volumes we calculate the discrete β-function (or step-scaling function) at two values of the renormalized coupling in the gradient flow scheme. In all cases perfect agreement is found between continuum extrapolated results using the topological and Wilson plaquette gauge actions. It is worthwhile to point out that a smooth action that effectively prohibits large plaquettes and is very close to the Wilson plaquette action for small plaquettes was investigated recently in [6]. The organization of the paper is as follows. In section 2 the topological lattice action and the details of our simulations are given. Section 3 is dedicated to the topological susceptibility and scale setting, in section 4 the results related to the deconfinement phase transition are discussed and in section 5 the results for the discrete β-function are presented. Finally in section 6 we end with a conclusion and possible future aspects.

Topological lattice action
We seek an action which is gauge invariant. The simplest possibility is to use the usual plaquette P as the only building block, where the sum is over all plaquettes on the lattice. On a given lattice volume the only parameter is δ which will play the role of a bare coupling. Clearly, as δ → 0 only links close to unity are allowed hence δ → 0 will correspond to the continuum limit. Note that the action (2.1) has only two values, 0 or ∞, hence divides the space of links into two subsets. On one, which contain the unit links, the action is zero and hence the links fluctuate freely without any weight and contribute equally to the path integral. In particular, the vacuum is infinitely degenerate. On the second set of links the action is ∞ meaning that links are forbidden there and contribute nothing to the path integral. The only dynamical information is carried by the boundary of these two sets defined by δ. As the continuum is approached, δ → 0, links have less and less phase space to fluctuate but still always have equal weight. Gauge invariance is encoded in a gauge invariant definition of the boundary between allowed and not allowed links.
It is well-known [7] that if the plaquettes on a lattice are all restricted to be very small, 1− 1 2 Tr P < 0.015, then a geometric integer definition of the topological charge Q can be given. A slightly more permissive bound was later derived in [8], 1 − 1 2 Tr P < 1 12(2+ √ 2) ≃ 0.0244. Hence for very small bare couplings δ a local algorithm can not change topology. In practice though the values of δ we use in this work are much larger, δ ∼ 0.6 . . . 0.8, and we do encounter topology change frequently enough in all large volume runs. In section 5 the simulations are done in very small physical volumes where of course Q = 0 but this is not an algorithmic artifact but rather the consequence of being in the femto world. Even in this case δ > 0. 3. In all runs with the topological action we use a simple Metropolis algorithm whereas with the Wilson plaquette action a heat bath algorithm. Both Metropolis and heat bath sweeps are accompanied by two to five overrelaxation steps [9]. An allowed configuration by the topological action may turn into a forbidden configuration by an overrelaxation step, in this case the step is rejected and the original configuration is kept.

Topological susceptibility and t 0 scale
The first observable we would like to compare in the continuum is the topological susceptibility. The gradient flow t 0 scale [10][11][12][13][14][15] is used to make it dimensionless and set the scale, hence we will compare χt 2 0 = Q 2 t 2 0 /V . For both the topological action and the Wilson plaquette action the measurement of the gradient flow is done in the same way, using the plaquette discretization along the flow and the symmetric clover discretization for the observable E(t), where on the right hand side the choice 0.3 could in principle be different. Choosing a smaller value would lead to smaller finite volume effects, smaller errors but larger cut-off effects. The topological charge Q is also measured along the flow at t = t 0 requiring no renormalization for the topological susceptibility. We confirmed that the continuum results are insensitive to the choice of t as long as it is kept fixed in physical units. For example, using t = 0.75 t 0 . . . 1.25 t 0 for the susceptibility leads to identical results.
Finite volume effects are ensured to be negligible within our statistical errors by always using symmetric L 4 lattices such that √ 8t 0 /L < 0.25. Another indicator that finite volume effects are small is that we always have T c L ≥ 4 (see next section).
The measurement of t 0 /a 2 in our simulations is very precise, its relative error is at least an order of magnitude smaller than the relative error on the susceptibility. Hence the final errors are completely dominated by the latter. The data is shown in table 1, the generated number of configurations at each point is O(10 5 ) with O(100) configurations separating the measurements.
The continuum extrapolation of t 2 0 χ in a 2 /t 0 is straightforward and shown in figure 1. The continuum results for the two discretizations agree, 0.00159(4) and 0.00165(4) for the topological and Wilson plaquette actions, respectively.

Deconfinement phase transition
Next we compare quantities which are sensitive to the deconfinement phase transition which is second order for SU (2). These can again be compared to the results obtained with the Wilson plaquette gauge action or with the corresponding quantities in the 3-dimensional Ising model. First, using the Binder cumulant of the Polyakov loop we will determine the critical exponent 1/ν with the topological action and find that it agrees with 1/ν = 1.5878(4) from the 3-dimensional Ising model [16]. Then the critical couplings are determined on N t = 4, 6, 8, 10 lattices and the dimensionless ratio T c √ 8t 0 is obtained in the continuum. Again perfect agreement is found between the two actions.
We will use standard scaling theory for much of this section; see [17] for more details. The bare parameters, δ or β, are collectively denoted by α. Since we investigate the deconfinement phase transition at fixed temporal extent, the dependence on N t is often suppressed while the spatial volume is denoted by N s . The Binder cumulant of the Polyakov loop L is defined by The reduced temperature will be denoted by x = (α − α c )/α c . Scaling theory says that close to the critical point a scaling function f can be defined such that g(α, N s ) = f xN 1/ν s with some critical exponent ν. Expanding around x = 0 we obtain which means that by considering two spatial volumes, N s and bN s , for some b > 1, we may estimate the exponent as where S(N s ) is the slope of g(α, N s ) with respect to α close to the critical point.    Next, we turn to the determination of the critical couplings. Let us denote the intersection of the Binder cumulants corresponding to N s and bN s by α (N s , b). The dependence on N s and b is again fixed by scaling theory and for large enough volumes we have, with some constant A, where y 1 = −1 for the 3D Ising model. Having established that our result for 1/ν is compatible with the 3D Ising model, we will simply use y 1 = −1 and 1/ν = 1.5878 in the above infinite volume extrapolations. The infinite volume extrapolations, using (4.4) is shown in figure 3. Estimating the statistical error    on α c is non-trivial since the pair-wise intersections for various N s are correlated; we use a jackknife procedure starting from the independently measured Binder cumulants. The critical couplings with the Wilson plaquette action have been determined in [17] for N t = 4, 6, 8. Our results for the Binder cumulants and infinite volume extrapolation for N t = 10 is shown in figure 4.
The final results for the critical couplings are listed in table 3 which we subsequently use to determine T c in physical units, T c √ 8t 0 . In order to determine the dimensionless combination T c √ 8t 0 in the continuum only t 0 /a 2 needs to -  be measured at the critical couplings. Even though the statistical error of t 0 /a 2 is very small in a given simulation, the uncertainty of the critical coupling itself needs to be taken into account. By a simple interpolation of t 0 /a 2 in δ and β we estimate this uncertainty originating from the uncertainty on the critical couplings. It turns out that this is the dominant source of final uncertainty, the statistical error is negligible. In table 3 we list the results and the reason for the unusually large error on t 0 /a 2 , relative to table 1, is the aforementioned effect. Another way of saying it is that if only the central values of β c and δ c are taken, then there is an uncertainty on N t leading to an additional uncertainty on 8t 0 /a 2 /N t beside t 0 /a 2 . The continuum extrapolation is shown in figure 5, the results are 0.851(8) and 0.840(6) with the topological and Wilson plaquette action, respectively. Again, complete agreement is found.

Small volume, perturbative regime
In this section we compare quantities in the perturbative regime. The simplest example is given by the discrete β-function or step-scaling function [18] in small physical volume or femto world [19][20][21] We will use the finite volume gradient flow scheme with periodic boundary conditions [22,23]. In this scheme the renormalized coupling is defined by where c = √ 8t/L is a constant so g 2 R only depends on one scale, L, the linear size of the periodic box. The factor is such that at tree level the above scheme agrees with MS, ϑ is the 3rd Jacobi elliptic function [22,23]. We set c = 3/10. The boundary condition is periodic in all directions and it is well-known that there are degenerate perturbative vacua within this setup. However for small renormalized couplings the system fluctuates around one of the vacua and tunnelling events are suppressed. In our simulations we work in this regime and indeed do not detect tunnelling at all as expected.
The expansion of g 2 R (L) in terms of g MS is unusual in a periodic 4-torus in the sense that both even and odd powers can appear, as well as non-analytic logarithms [21]. Particularly for SU (2) purely logarithmically suppressed terms appear right after tree level; see [22] for an extended discussion. Hence only the first coefficient of the β-function is the same in the gradient flow scheme and MS. It is important to note that the scheme as such is completely well-defined non-perturbatively, merely the perturbative expansion behaves in a somewhat unusual way. The β-function at fixed g 2 R is a perfectly well-defined and universal quantity. The subscript R will be dropped in what follows.
On the lattice the discrete β-function, (g 2 (sL) − g 2 (L))/ log(s 2 ), or step-scaling function is the most easily accessible quantity with a well-defined continuum limit for fixed g 2 (L). We set s = 3/2. It has a perturbative expansion where, as mentioned above, the terms in . . . contain both even and odd powers of g(L) as well as logarithms. The 1-loop term b 0 = 22/3 is nevertheless universal.
We computed the discrete β-function at two values of the renormalized coupling, g 2 (L) = 1.5 and g 2 (L) = 2.5. As we will see the first is small enough such that the continuum result is compatible with the 1-loop approximation (5.3). The continuum limit is approached by 3 sets of lattice volumes, 16 4 → 24 4 , 18 4 → 27 4 and 24 4 → 36 4 corresponding to the scale change s = 3/2. The desired values of g 2 (L) were tuned on the smaller lattices to high accuracy by tuning the bare couplings of the two actions, β and δ, respectively. Then at the same values of the bare couplings the renormalized couplings were measured on the larger lattices. Finally the results were extrapolated to the continuum linearly in a 2 /L 2 .
The tuned couplings and the measured values on the larger lattices are shown in table 4. The continuum extrapolations are shown in figure 6, all χ 2 /dof are less than unity. As emphasized above the 2-loop result is only shown for orientation, it is not universal in our scheme. There is perfect agreement for the continuum results between the two actions. The smaller coupling, g 2 (L) = 1.5  (7) L/a β g 2 (L) g 2 (sL) β g 2 (L) g 2 (sL) 16 4.68515 1.500 (1)   is expected to be small enough such that renormalized perturbation theory at 1-loop is trustworthy. This seems to be the case despite the only logarithmically suppressed terms and the results with the Wilson plaquette and topological actions are compatible with the 1-loop approximation (5.3) within 1 σ and 1.3 σ, respectively. Even though bare perturbation theory is not possible to set up, renormalized perturbation theory behaves as expected.

Conclusion and outlook
In this paper we investigated non-abelian lattice gauge theory with an action that is insensitive to small perturbations of the lattice fields. In particular, the classical vacuum is infinitely degenerate. By comparing continuum extrapolated results to those obtained with the Wilson plaquette action (which is known to be in the correct universality class) we conclude that even though the topological action has no classical continuum limit, the quantum continuum limit correctly reproduces the theory. There are of course only a finite number of comparisons that one can do in simulations but we have chosen quantities that span a wide range of interesting phenomena. The topological susceptibility, the critical temperature, critical exponents, and the β-function in small physical volumes were compared and perfect agreement was found. This suggests that universality on the lattice is very robust; one is free to modify the lattice action not only with higher dimensional irrelevant operators but also the relevant operators do not have to be reached in a smooth way dictated by the classical continuum Lagrangian. It may be worthwhile to remember that a positive transfer matrix can not be associated with a topological action at finite lattice spacing [24]. However, as duly pointed out in [24], this is not necessarily a problem if the effect of non-positivity disappears in the continuum limit. In other words if the positivity violation is merely a cut-off effect. Our results indicate that this is indeed the case since in the continuum the correct Euclidean Yang-Mills theory is obtained.
A less obvious question concerns the form of scaling violations, O(a 2 ) specifically, confirmed numerically in our work. Since bare perturbation theory is not applicable it is not immediately clear that the same scaling is expected as with the Wilson plaquette action. However, following the same line of reasoning as for the 2-dimensional O(3) model in [1], Symanzik's effective theory, formulated in the continuum, still applies and leads to the same O(a 2 ) scaling violations as with the Wilson plaquette action.
An intriguing application of topological actions might be with dynamical fermions. The Nielsen-Ninomiya no-go theorem [25] heavily relies on the classical action, namely that in momentum space the continuum Dirac operator is reproduced in the classical continuum limit. With a topological fermionic action the classical continuum limit is meaningless. Hence perhaps the no-go theorem may be circumvented by such actions although it is of course not at all clear how they would be first defined and then implemented.