Symanzik improvement of the gradient flow in lattice gauge theories.

We apply the Symanzik improvement programme to the [Formula: see text]-dimensional local re-formulation of the gradient flow in pure SU(N) lattice gauge theories. We show that the classical nature of the flow equation allows one to eliminate all cutoff effects at [Formula: see text], which originate either from the discretised gradient flow equation or from the gradient flow observable. All the remaining [Formula: see text] effects can be understood in terms of local counterterms at the zero flow-time boundary. We classify these counterterms and provide a complete set as required for on-shell improvement. Compared to the 4-dimensional pure gauge theory only a single additional counterterm is required, which corresponds to a modified initial condition for the flow equation. A consistency test in perturbation theory is passed and allows one to determine all counterterm coefficients to lowest non-trivial order in the coupling.

Introduction

In recent years the Yang–Mills gradient flow has been established as a very promising new tool to study non-perturbative aspects of strongly coupled gauge theories [1-4]. The gradient flow defines a deterministic mapping from the original gauge field to a smoothed gauge field configuration, , at flow time t, which is obtained as the solution of the gradient flow equation (see Appendix A for a summary of our conventions),where denotes the gauge covariant derivative andis the associated field strength tensor. The name relates to the fact that the right hand side of (1) is equal to minus the gradient of the Yang–Mills gauge action. Hence, with increasing flow time t, the solution, , is driven towards a minimum of the action and thus approaches a smooth classical field configuration.

There is quite some freedom when translating the gradient flow equation to a Euclidean space-time lattice. A simple possibility is to choose Wilson’s plaquette action, , and to define the lattice gauge field at finite flow time, , as the solution of the Wilson flow equation,where denotes the Lie-algebra valued derivative with respect to . It should be noted that similar smoothing operations have long been successfully applied in lattice QCD. For example, the stout link smearing technique of Ref. [5] can be understood as a discretised flow time version of Eq. (3). The essential new element is a theoretical understanding of the renormalisation properties of the Yang–Mills gradient flow. In particular, in [3, 4] it was proved to all orders of perturbation theory that QCD at finite flow time t is renormalised once it is renormalised at flow time through the usual renormalisations of the gauge coupling and the quark mass parameters. Furthermore, gauge invariant fields at positive flow time are automatically renormalised and do not mix with other fields of the same or lower dimensions. These properties allow one to define a new class of renormalised gauge invariant observables which can be used to probe the theory in various ways. It also opens new ways to define renormalised composite operators at zero flow time; the study of Ward identities at positive flow times [4, 6, 7] and the applications of the so called “small flow-time expansion” have received much attention recently in this context [8-10].

Many current lattice QCD applications of the gradient flow only involve the simplest possible gauge invariant field, the action density,As initially proposed in [2], the expectation value can be used for a non-perturbative definition of either a reference scale or a coupling constant. This has proven very attractive: in large volume simulations it leads to the most precise determination of a reference scale (for a recent review cf. [11]). On the other hand, when considered in a finite space-time volume the scale evolution of the corresponding coupling [12-16] can be traced with high statistical precision (see [17] for a recent review).

Notwithstanding these nice properties a major practical problem consists in the relatively large cutoff effects which have been observed in several applications (cf. [17] and references therein). On general grounds, the leading effects are expected to be of order . Their size depends on the detailed choices made when translating the flow equation (1) to the lattice, but also on the discretisation of the observable and on the lattice action. Alternative flow equations have been tried, e.g. in Ref. [18] where the Wilson action was replaced by the tree-level improved Lüscher–Weisz action,  [19, 20]. For some attempts to reduce cutoff effects in the particular observable cf. Refs. [21, 22]. Here we would like to proceed more systematically by applying the Symanzik procedure [19, 23] to the -dimensional local formulation of the theory [3, 24]. This will lead us to a particular choice for the lattice flow equation, referred to as the “Zeuthen flow” and defined bywith the initial condition . Here and are the lattice forward and backward covariant derivatives, respectively. We will show that the integration of the Zeuthen flow equation does not generate any cutoff effects at ). If combined with classical improvement of the observable all effects are eliminated apart from those corresponding to local counterterms in the action at zero flow time. We will give a complete list of such counterterms and test our framework to lowest non-trivial order in perturbation theory.

The paper is organised as follows: in Sect. 2 we recall the definition of the -dimensional local theory, with flow time as the added dimension. In Sect. 3 we discuss the general Symanzik procedure and the simplifications due to the special properties of this theory. We present the classical a-expansion of both the flow action and the gradient flow observable E(t, x), as part of the simplified Symanzik procedure, and carry out the standard Symanzik analysis for the counterterms at the boundary. Section 4 presents a number of perturbative tests of the improved theory, and Sect. 5 our conclusions. We have included three appendices regarding our notations and conventions (Appendix A), some technical details pertaining to the classical a-expansion (Appendix B), and some explicit expressions used in Sect. 4 (Appendix C), respectively.

Lattice gauge theory in dimensions

The gradient flow equation can be viewed as a way to define a particular class of observables, i.e. fields which are functionals of the fundamental gauge field . The flow time thus appears as an additional parameter which measures the range in space-time over which the fundamental gauge field enters into an observable defined in terms of the flowed gauge field . The flow time t has dimension length squared and the “smearing radius” is usually taken as the corresponding length scale.1 Thus, gradient flow observables are non-local objects from the perspective of the 4-dimensional gauge theory and their properties under renormalisation are difficult to assess. Moreover, the non-locality prevents a straightforward application of the Symanzik expansion, which is our main theoretical tool for understanding the cutoff dependence of the theory. For this purpose, it is therefore highly beneficial to follow [4] and view the theory from a -dimensional perspective, with flow time as the added dimension. In this re-formulation locality is restored in the -dimensional sense, and dimensional counting can be applied to classify counterterms to the action and observables.

We start with the formulation of the lattice set-up, including the introduction of a flow-time lattice. The latter should be regarded as an intermediate regularisation which helps to resolve certain technical issues [4]. While none of this is original it serves for later reference and to fix our notation.

The 4-dimensional lattice action

On-shell improvement of the 4-dimensional gauge theory can be achieved by introducing, besides the 4-link plaquette action, further 6-link Wilson loops with appropriately chosen coefficients [19]. We will consider a general class of lattice gauge actions parameterised by the coefficients , defined by,where the second sum extends over all oriented Wilson loops of type . As illustrated in Fig. 1, these Wilson loops are the usual plaquettes, , the planar loops or “rectangles”, , the bent rectangles or “chairs”, , and finally the “parallelograms”, .

Fig. 1

Wilson loops of type and

It is conventional to normalise the action by requiringsuch that the standard continuum Yang–Mills action is obtained in the classical continuum limit, with any choice of the three free parameters.2 Popular choices are the Wilson plaquette (W) action (, ) and the tree-level improved Lüscher–Weisz (LW) action (, , ).

-dimensional set-up

Given the 4-dimensional action, the flow equation is now incorporated in the action as a constraint, by introducing the Lagrange multiplier field , which is hermitian and such that is Lie-algebra valued. The -dimensional action of this theory then takes the formwhere the boundary condition,is assumed andis a shorthand notation which allows one to write the lattice gradient flow equation in the form . The action is some 4-dimensional lattice gauge action for the flowed field , the Wilson action being the simplest choice [cf. Eq. (3)]. In any case it is unrelated to the gauge action in (8). How to best define is at the core of this work and will be discussed in the next section.

Given the action for the -dimensional half space , expectation values of composite fields O[V, L] are defined as usual,with the condition . A few remarks are in order: first, the integration over the gauge field includes the integration over its boundary values at i.e. the standard 4-dimensional gauge field . Hence, for observables which only depend on , the functional integrals over and cancel between numerator and denominator, reproducing the standard expectation value of the 4-dimensional theory. To see this more explicitly it is convenient to pass to a flow-time lattice with spacing and lattice points  [4],where we have assumed the discretisation,with the correct limit. Inserting this representation of the action into the functional integral, the integration over the fields produces a string of -functions3These can be eliminated one by one, by integrating over for strictly positive n, leaving the unconstrained integration over the fundamental gauge field intact, as expected.

Gauge symmetry

By construction, the -dimensional action is gauge invariant under t-independent gauge transformations,where is an SU(N)-valued gauge function. This leads to the transformation,so that gauge invariance of the action is guaranteed provided thati.e.  must be in the adjoint representation of the gauge group. The integration measure is invariant under such a change of variables, so that the gauge symmetry of the 4-dimensional boundary theory is inherited by the bulk theory.

It is occasionally useful to generalise the gauge symmetry to the flow-time coordinate t, i.e. admit flow-time dependent gauge functions . In the continuum theory this amounts to replacing , and  [1]. In the presence of the lattice cutoff (but continuous ) we define the covariant -derivative bywhere . This, together with the transformation under an -dependent gauge transformationleads to the left hand side of the covariant flow equation transforming asRendering the t-derivative covariant in the definition of (Eq. 10), one then obtains,so that -dimensional gauge invariance is established, provided that transforms just like . Discretising the flow-time coordinate is also straightforward, one just needs to elevate the fourth component of the gauge field to a link field , with corresponding changes in the covariant derivative and gauge transformation behaviour.

Finally we note that the Yang–Mills flow equation in the continuum can be written aswhich shows that the -dimensional theory, while exactly gauge invariant, does not enjoy any generalised Lorentz symmetry. This is of course already clear from the dimensions, in particular, and thus must have mass dimension 2, in contrast to the usual derivatives and gauge fields in 4 dimensions.

Symanzik improvement to

Generalities

The re-formulation of gradient flow observables in terms of a local -dimensional lattice gauge theory creates the standard situation to which Symanzik’s effective theory [23] can be applied in the usual way. We start with Symanzik’s effective action which is given as an expansion in powers of ,One might worry about odd powers of a arising in a -dimensional theory. However, as we will show in detail in Appendix B, gauge invariance, reflection symmetries and the fact that the flow-time parameter t has mass dimension imply that non-trivial counterterms to the action must be even dimensional. In Eq. (23) we have separated the effective action of the flow in the -dimensional volume, , from the action with support restricted to the 4-dimensional boundary at . Both parts will be discussed in turn below. Besides the effective action, also local observables are described by an effective continuum field, again expanded in powers of . For a generic local observable we writeTo the Symanzik expansion of lattice expectation values then takes the form,Here, the expectation values on the RHS are defined in the continuum theory with respect to the continuum action , and the notation serves as a reminder that only the connected part contributes to the correlation functions with counterterm insertions, for instanceAs the next step in the Symanzik procedure one determines a basis of counterterms both for the action and the observables of interest. In the case of the action these take the formwhere the fields are gauge invariant polynomials in the fundamental fields , and their (space-time and/or flow-time) derivatives, and the are similarly constructed, but evaluated at . Since must be dimensionless, the fields must have mass dimension 8 and otherwise share all the symmetries with the lattice theory. The fields are dimension 6 fields, localised at the boundary. One of the important outcomes of the Symanzik analysis are the numbers and of basis elements, where fields differing by total (space-time) derivative terms are considered equivalent. Furthermore, restricting to on-shell improvement the field equations for , and can be used to simplify the basis. Given a basis of counterterms the final step of Symanzik’s procedure consists in adding lattice representatives of these operators to the lattice action, such that, with appropriately chosen coefficients, the terms and are eliminated in Symanzik’s effective action for the improved lattice action.

A similar analysis then needs to be carried out for each observable of interest, i.e.  in Eq. (24) is given as a linear combination of local fields of mass dimension , which share all the lattice symmetries with . While this procedure applies to any observables, we will here focus on gradient flow observables, i.e. gauge invariant composite fields with support at strictly positive flow times.

If the full Symanzik procedure as outlined above were really necessary, improvement would probably remain an academic curiosity. In particular, a rather long list of dimension 8 counterterms for could be written down, with little hope for practical relevance, so that one might be tempted to give up on systematic improvement.

Before proceeding along these lines, however, it is advisable to have a closer look at this particular theory. As shown by Lüscher and Weisz, the theory is perturbatively renormalizable to all orders in the 4-dimensional gauge coupling g [3]. More precisely, if one restricts attention to gauge invariant observables, one just needs to renormalise the gauge coupling in the usual way, and also the quark masses if the boundary theory is generalised to QCD.4 Moreover, any composite fields defined at finite flow time are automatically renormalised and do not mix with any other fields of the same or lower canonical dimension. The action density (4) is a typical example: its renormalisation at flow time requires the subtraction of both a quartic and a logarithmic divergence. None of this is required at finite t. It is instructive to consider leading order perturbation theory to get a basic understanding of the mechanism at work. Effectively, at finite flow time t, integrals over the loop-momentum p are cut off by an exponential suppression factor in the integrand. This renders most momentum integrals finite, so that one is only left with those divergences which are cancelled by the standard counterterms in the boundary theory.

Hence the -dimensional theory enjoys rather special properties. In particular, the field plays the rôle of a Lagrange multiplier field which enforces the gradient flow equation as a constraint. The smoothening properties of this equation are related to the fact that perturbation theory only generates tree diagrams for the correlation functions of gradient flow observables [3]. The Symanzik expansion is then very much simplified as we expect the following to hold:In the following we first remind the reader of the classical a-expansion and then address these points in the subsequent subsections one at a time.

The absence of bulk loop diagrams in the perturbative expansion of gradient flow observables implies that classical improvement of the flow action cancels the effects exactly, i.e. without any corrections.

By the same argument, non-perturbative improvement of composite operators at positive flow time can be achieved by choosing discretisations that do not generate effects when expanded classically.

The only counterterms which receive genuine quantum corrections are the ones living in the 4-dimensional boundary at . The full Symanzik procedure outlined above thus needs to be applied only to the boundary part, , of the Symanzik action, and of course to any observable which is at least in part localised at the boundary.

The classical a-expansion

According to the preceding discussion the counterterms appearing in and in for gradient flow observables are completely determined by classically expanding the lattice action in the -dimensional volume and the observables under consideration to order . The classical expansion assumes that the lattice approximates an underlying continuum space-time manifold on which a smooth continuum gauge field, , is defined. The lattice gauge field, , is then related to the continuum gauge field by parallel transport along the lattice links. Parameterising the path along the lattice link from to x by (with parameter ), the precise relation is obtained by iteratively solving the differential equation,The solution, , can be concisely written in terms of a path-ordered exponential,While it is straightforward to carry out the expansion around , in practice, even a simple gauge invariant quantity like the trace of the plaquette contains four link variables which need to be expanded and combined to fourth order in a to obtain the leading non-trivial term. It is therefore highly advisable to perform the expansion efficiently (cf. e.g. [19, 25]). We here follow Lüscher and Weisz [19], who, for fixed indices and , proposed to work in the following gauge:As a result, the expansion around is very much simplified. For example, the plaquette field,is reduced to a single link,Recalling the definition of the path-ordered exponential (30) one needs the expansion of the B-field around ,where is a constant and neglected terms are of order . Following [19] the gauge covariant expressions can be unambiguously restored, with the resultInserting into the path-ordered exponential with appropriate replacements for , we thus obtain the gauge covariant expansion for the plaquette field,which holds for any argument (t, x). Similar expressions can be derived for the other three plaquettes in the plane:and the next few orders can be obtained with moderate additional effort.

Determination of

To find the bulk counterterm action we simply need to apply the classical expansion to the bulk action in Eq. (8). This essentially amounts to the a-expansion of the gradient flow equation, i.e.  in Eq. (10). For the first term we find, in the Lüscher–Weisz gauge (31),as all other terms are proportional to . The Taylor expansion can easily be performed to all orders in a with the resultWe therefore expect that the correct gauge covariant expression at any lattice point x must readAt this point one may wonder whether the gauge covariant expression really follows unambiguously from the gauge fixed expansion, in particular, whether the t-derivative always has to be to the right of the covariant -derivatives. That this is indeed correct can be established by using the -dimensional gauge symmetry (cf. Sect. 2), which implies that the a-expansion of this term must be given as covariant derivatives acting on .

Turning to the second term of (10), i.e. the gradient force term, we choose a quite general lattice gauge action parameterised by , which includes all four- and six-link Wilson loops (plaquettes, rectangles, chairs) except the parallelograms. We decompose the action as follows:We first express the gradient force in terms of plaquettes and their covariant derivatives. For the plaquette action we then findwhere we have introduced the projection on the traceless antihermitian part, i.e. for an matrix M in colour space we defineFor the rectangle action we findand a similar but slightly more complicated expression is obtained for the chairs. Expanding each term to order and recombining them we getwhere the arguments (t, x) on the RHS have been omitted. Collecting all results we define the expansion coefficientswhere the leading term defines the continuum limit,Hence the correct normalisation to reproduce the Yang–Mills gradient flow equation (1) is , which we use to eliminate in the higher order terms:Before proceeding we remark on the presence of odd powers of a in the expansion, which seems at odds with our expectation that only even powers of a occur in this theory. The resolution of this apparent contradiction lies in the fact that the lattice fields and should be defined on the lattice link connecting x and , rather than at the lattice site x. In Appendix B we demonstrate how the covariant re-expansion about the midpoint of the link, , eliminates such terms. While this problem will not affect our discussion of the counterterms, it clarifies that the correction terms are indeed of order .

We now proceed and work out the simplifications due to the field equations for and . Varying the continuum action with respect to one obtains the Yang–Mills flow equation (1), whereas the variation with respect to yieldsUsing the flow equation eliminates the O(a) term , and this is the reason why the terms remain unaffected by the symmetrisation about the midpoint , once the field equations are taken into account. From the continuum flow equation we deriveThis allows one to rewrite the term as follows:From the corresponding flow action,one may now directly read off the counterterm structures that correspond with a given choice of the coefficients . Unfortunately, there does not seem to be a choice such that vanishes. We also attempted to use Eq. (46) as follows: considering the termone may perform an integration by parts with respect to t. This generates a surface term at ,which re-defines a coefficient of the counterterms entering (cf. Sect. 3.5). Equation (46) then leads to space-time derivatives acting on , which can be integrated by parts (no surface terms are generated here) to redefine . Unfortunately, this does not yield a solution with either. We notice, however, that with the Lüscher–Weisz choice of coefficients and , has a rather simple structure,To cancel this term is relatively straightforward. Starting from the lattice gradient force defined with the Lüscher–Weisz action, , we simply act withon this gradient force, which yields the “Zeuthen flow” equation (5). The flow action for the Zeuthen flow does indeed vanish, i.e. we have successfully implemented improvement in the -dimensional bulk.

improvement of E(t, x)

We here consider only the simplest observable, namely the action density E(t, x) of Eq. (4) The two most popular lattice discretisations of E(t, x) are referred to as plaquette (pl) and clover (cl) definitions, respectively. They are either obtained from the Wilson plaquette action or based on the so called clover leaf definition of the field strength tensor,which uses the four plaquettes (32), (35)–(37) in the plane. The plaquette and clover lattice versions of E(t, x) are now given byandPushing the classical a-expansion of the plaquette (32) to one obtainswith the continuum limit given by Eq. (4). Proceeding in this way for all four plaquettes of the clover leaf we obtain the classical expansionSeveral remarks are in order. First, the a-expansion of the plaquette yields contributions at every order in a, whereas the symmetries of the clover definition imply only even powers of a. The odd powers of a could be eliminated by re-defining as an average over the four plaquettes of the clover leaf, which, due to the trace operation, coincide with for appropriately displaced arguments x. Second, note the total derivative terms which may appear at any order in a. Such terms do not contribute to the expectation value , provided that the chosen set-up is translation invariant. This would e.g. be the case in a finite volume with periodic or twisted periodic boundary conditions, and thus in the limit of infinite volume. However, translation invariance no longer holds with either Dirichlet or Neumann conditions5 as required for the Schrödinger functional [13] or with open boundary conditions [14]. Similarly, when considering higher correlation functions such as the 2-point correlator of two fields E(t, x) total derivative terms cannot be ignored. We will here focus on the translation invariant case and from now on consider such total derivative terms negligible. This eliminates all the odd powers of a in the expansion of . Hence, both discretisations are on equal footing and counterterms for and are now easily identified as the coefficients in the classical expansion. Given both a-expansions we observe that the terms have the same structure, with the coefficients in the clover definition being larger by a factor of 4. In any case we observe that the linear combinationdefines an improved observable for which vanishes. An alternative improved definition of E(t, x) can be obtained from the action density of a tree-level improved lattice action such as the Lüscher–Weisz action (Eq. (6) with , and ). Here again, any ambiguity in the definition of a density from the action amounts to total derivative terms, which we consider negligible in the present context.

In this subsection we list the gauge invariant local fields of dimension 6 which may appear in the boundary action of Symanzik’s effective action. Disregarding total derivative terms with respect to the space-time coordinates x, we find the following list of seven candidate counterterms:where denotes the field strength tensor of the fundamental gauge field.

Again we apply the field equations. The Yang–Mills flow equation impliesso that, after taking into account the boundary condition , we haveThis eliminates . The field equation (46) is not useful here. However, a third field equation can be derived by varying the action at with respect to the fundamental gauge field . Technically this is best done by discretising only the flow time in the -dimensional continuum action and taking the limit of continuous flow time in the end. The resulting field equation is6This equation leads to the relationsAt this point it is useful to recall the situation in the standard 4-dimensional theory [19]. In fact there is a 1-parameter family of improved actions, which, to tree-level, are parameterised by as follows:Expanding the action classically, the free parameter is seen to multiply the counterterm . The counterterm is thus redundant for the improvement of standard observables. In principle one may thus tune the coefficients (66) to achieve improvement of both standard and gradient flow observables. In practice, however, these coefficients define the gauge action used in the Monte-Carlo simulation and the corresponding effective coefficient of should be regarded as fixed. Therefore we choose as our basis of counterterms (cf. Sect. 3.6).

Finally, we remark that the use of the field equation (64) in the counterterm basis holds for counterterm insertions only up to contact terms, namely whenever the counterterm argument coincides with the location of some field in the correlation function under study. Such contact terms are thus absent for gradient flow observables localised at strictly positive flow times. However, we expect these relations to hold more generally, i.e. even if some fields in the correlation functions are defined at zero flow time. In this case we expect that the contact terms which make the difference are of the same form as the counterterms to the fields in the correlation function and therefore just redefine these counterterm coefficients. This parallels the discussion in Ref. [26] of on-shell improvement in lattice QCD with Wilson quarks.

Summary of Section 3 and some practical considerations

Section 3 contains the main results of this paper and may appear rather technical. We therefore provide a short summary and comment on the practical implementation of the lattice counterterm .

There is a natural way of interpreting the gradient flow as a -dimensional local quantum field theory. The flow time t plays the role of the coordinate in the fifth dimension, which only takes on non-negative values (). The dynamics of the theory in the bulk () is completely fixed by the deterministic flow equation. The classical nature of the theory for allows one to implement the Symanzik improvement programme in a rather simple way: all cutoff effects produced by integrating the flow equation can be eliminated via a suitable discretisation of the flow equation, which can be determined by the classical expansion to . Similar considerations allow one to define discretised flow observables that are free from lattice artefacts. The only remaining effects are generated by the action at the boundary , and are genuine quantum effects. They correspond to the usual counterterms (55)–(57) in the 4-dimensional action affecting all lattice observables.

To implement an improved lattice action one first has to choose an improved 4-dimensional lattice gauge action which amounts to choosing coefficients in Eq. (6) appropriately. It is well known how improvement can be implemented at tree-level, and also to order in the case of the pure gauge theory [20]. In addition one needs to incorporate a lattice version of such as to cancel the insertion of on observables without changing the coefficients .

To achieve this we remind the reader that the -dimensional set-up is used only for the theoretical analysis, whereas in practice one integrates the gradient flow equation numerically and evaluates any observable such as E(t, x) along the flow. It turns out that the insertion of can be realised by a change in the initial condition at for the gradient flow equation. Since in this case and are not the same we need to fix the integration variables in the -dimensional field theory. We choose to integrate over the fundamental gauge field and the flow field for . Therefore on the lattice we choose to integrate over and for . A shift in the initial condition can be implemented viawhere is the free improvement coefficient, and any 4-dimensional lattice action. In the -dimensional formulation with -discretised flow time, the fields and only enter in the termswhere is, up to terms of , the RHS of the flow equation. Now we can trade all references to into , that is our path integral variable. Using Eq. (67) we can writewhere the dots represent higher order terms in the lattice spacing. Therefore the shift in the initial condition is equivalent (up to higher order corrections in a) to the insertion of the countertermThis can be reinterpreted as the previous situation with standard boundary conditions, , except for the extra term in the action. Hence we have successfully traded the modified boundary conditions for the flow equation for the term in the lattice action. In the next section we will determine its coefficient at tree-level of perturbation theory.

Perturbative analysis

In this section we will study the Symanzik improvement of the gradient flow in perturbation theory. This will allow us first to determine the improvement coefficient to tree-level. Second, the study of the Zeuthen flow both in small volumes and for different observables will allow us to check explicitly that the use of a tree-level improved action for the simulation together with the tree-level value7 of , the Zeuthen flow and a classically improved definition of the observable yields expectation values that are free from effects at tree level. As observable we choose first E(t, x). The contributions by the action, flow and observable to the cutoff effects of at tree level have been computed recently [22, 27]. Here we will show that the tree-level cutoff effects are absent not only in infinite volume, but also in a finite volume with twisted periodic boundary conditions, where the additional scale L leads to more stringent tests. Second we will consider the connected correlation function for E(t, x)E(s, y) and show that improvement by the Zeuthen flow is also obtained in this case.

Gauge fixing

In perturbation theory one parametrises the links in a neighbourhood of a classical configuration as follows:Note that this standard convention implies a re-scaling of the fields,compared to the preceding sections. In perturbation theory it is convenient to use gauge symmetry to simplify explicit computations. In the context of the gradient flow, gauge fixing is performed by studying the generalised flow equationThe superscript serves as a reminder that covariant derivatives and field strength are made of the modified flow field , i.e. the solution of the above equation. Note that the original flow equation is recovered by setting . The key observation is that gauge invariant observables are independent of  [1, 3, 24]. In order to see this, one only has to check that the gauge transformationwheretransforms a solution of the flow equation with arbitrary into one with .

On the lattice the procedure is completely analogous. We consider the generalised flow equationor, for the case of the Zeuthen flow,with initial condition . One then easily verifies that the gauge transformationtransforms a solution with an arbitrary function into one with . A natural choice for the function then is given as the solution of the equation,Note that this is a particular application of the -dimensional gauge transformations described in Sect. 2.3 and it is thus clear that gauge invariant observables remain unaffected by the choice of . This can be turned around to provide checks on the correctness of a given calculation. In the following we drop the indices (or ) from the fields and we will quote any intermediate results in Feynman gauge (). Some elements used for our checks of gauge parameter independence are given in Appendix C.

Determination of to tree level

We first assume that the lattice is infinitely extended and expand the general class of actions, Eq. (6), to leading order in the coupling,8where is a gauge fixing parameter and explicit expressions for the lattice kernels, , are given in Appendix A.

Similarly, the flow equation contains the gradient of a lattice action which, to leading order in the coupling, is parameterised by another action kernel, . The flow equation to this order then takes the form of the heat equation,The initial condition for the flow equation Eq. (67) reads to leading order in the fields,9where is yet another action kernel. No gauge fixing term is required here, so that the gauge parameter is set to zero. The linearised flow equation (77) can now be solved easilywhere is the heat kernel given byNote that we have used here as a matrix with respect to the Lorentz indices and the exponential has to be taken of that matrix. In the following we will often make use of such a matrix notation, in order to avoid an abundance of Lorentz indices.

Finally, the observable E(t, x), being an action density, can be parameterised by a further lattice action kernel, , with gauge fixing parameter set to zero. To this order we then obtain for the expectation valuewhere the trace is over Lorentz indices only and the gauge field propagator at positive flow time is defined byDue to the relation (78), this propagator depends implicitly on both gauge parameters, and , of the flow equation and of the action, respectively. Introducing the standard 4-dimensional gauge field propagatorthis propagator is the matrix inverse of the action kernel,and the gauge fixing parameter must be non-zero for the inverse to exist. Using these ingredients, the gauge field propagator at positive flow time can now be written as follows:where we have denoted the matrix transpose by the superscript T.

In summary, the choices of action, flow and observable discretisation correspond to the choice of three action kernels. Finally the shift in the initial condition is encoded in a fourth choice of kernel. Explicit expressions for some popular choices of kernels are given in Appendix C.

In order to obtain the leading order cutoff effects we now expand the kernels as follows:where the continuum kernel is given byUsing the continuum kernel only and neglecting cutoff effects we thus obtain the well-known continuum result in infinite volume,Explicit expressions for the correction terms are given in Appendix C. In order to compute the leading correction to the propagator and to the heat kernel it is convenient to work in Feynman gauge (), since in this case is proportional to . Working in a general gauge is, however, not much more difficult and serves as a check that the gauge dependence actually cancels in the final evaluation of the observable. A few technical details pertaining to such a check are given in Appendix C.

In the following we will use Feynman gauge and remove the gauge parameters as arguments of the action and flow kernels. We will also omit them in the kernels for the observable and initial conditions however, with the understanding that they must be set to zero in these cases. In Feynman gauge () it is straightforward to check that and finally, putting all the pieces together and after some algebra, we getwhere the constants are defined byandIn fact it is straightforward to evaluate the integrals, with the resultThe coefficients must be independent of the gauge parameters and and we have checked this explicitly. Their values depend on the choices made for the various kernels. For example, for a general action of the form Eq. (6) we have

Table 1 summarises the values of the coefficients for the most common choices. It is easy to see that the use of the Zeuthen flow together with the tree-level improved Lüscher–Weisz action and any classically improved discretisation for the observable (see Sect. 3.4) has no tree-level cutoff effects as long as . Therefore, to tree-level, the Lüscher–Weisz action () produces tree-level improved results for gradient flow observables. For the case of a generalised tree-level improved action Eq. (66) we have to choosein order to obtain tree-level improvement.

Table 1

Values of the coefficients in the terms of in infinite volume. The one-parameter family of tree-level improved actions corresponds to the choice of coefficients Eq. (66), the Lüscher–Weisz tree-level improved action being the particular choice with

Discretisation	d_1	d_2
Plaquette	-1/12	0
Lüscher–Weisz	-1/36	1 / 36
43 Plaquette -13 Clover	-1/36	1 / 36
One-parameter tree-level improved	-1/36	1/36-x_p
Clover	-1/4	-1/12
Zeuthen	0	0

As the reader can see, besides the Zeuthen flow there seem to be many ways to cancel the tree-level effects (see also [22]), as these are encoded in a single term, once the numerical values of (90) and for (cf. Table 1) are inserted into Eq. (87). We are thus led to look for more stringent tests of improvement by looking at a variety of observables and/or kinematics. After all, rather than improving a particular observable in a specific situation (e.g. in infinite volume), Symanzik improvement is designed to work for any observable in both finite and infinite volume.

Twisted periodic boundary conditions

A stringent test of our computations can be made when studying in a finite volume. Due to the presence of a new scale L, the cutoff effects will in general depend on the dimensionless ratio . Improvement requires that the tree-level cutoff effects vanish for all values of.

As a finite volume renormalisation scheme, we will use twisted boundary conditions for our gauge field. In this set-up, the gauge field changes by a gauge transformation when displaced by a period. Gauge invariant quantities are still periodic, and the absence of zero-modes turns out to be very convenient for perturbative computations. The gradient flow has already been studied in this set-up, and we will not give many details here but refer the interested reader to Ref. [15] and the references cited therein.

We will only need the perturbative expression of to leading order, given by with

Note that the expression is almost identical to the infinite volume one, except that the momentum integral has been substituted by a sum (hardly a surprise). The particularities of the twisted boundary conditions are hidden in the sum and momentum symbols. First notice that the momentum (with capital letters ) can be uniquely decomposed aswith andi.e. there is the usual space-momentum, but in the directions of the twisted plane the momentum lives in an apparently larger lattice of physical extent NL. Finally the sum symbol means sum both over and , but without the terms with . In particular the sum has no term with a zero total momentum. Notice that the colour factor is produced by the sum over .

The algebra is very similar to the one of the previous section, with the important difference that now the sums actually depend on the dimensionless ratio . In fact fixing the flow time in units of the volume in this way we getwhereand the third Jacobi theta function readsFinally the functions are given by10

In the limit , we recover the expressions of the infinite volume, in particularbut for non-zero c the functions are in general linearly independent (see Fig. 2). The coefficients are still the same, and the reader can check that the tree-level cutoff effects given by expression Eq. (96) vanish for all values of c when one uses our improved set-up (i.e. Lüscher–Weisz action, Zeuthen flow and Lüscher–Weisz observable). Any other choice of improved action together with the appropriate choice of also does the work. For this to happen it is crucial that the flow coefficients are both zero, since the functions and are linearly independent. In particular it is easy now to check that the so called Symanzik flow in the literature [28] or any set of coefficients in [22] does not remove the tree-level cutoff effects in finite volume. For the Zeuthen flow both coefficients do identically vanish, so that effects are indeed removed as expected on theoretical grounds.

Fig. 2

Ratio of the finite volume functions (Eq. 90) over the infinite volume predictions (Eq. 99). When there are significant differences between them. Moreover, the different functions are in general linearly independent

The connected 2-point function of E(t, x)

Further tests of the Zeuthen flow can be obtained by considering different operators at positive flow time. In particular, we now consider the 2-point function of E(t, x) in a periodic box of size L with twisted periodic boundary conditions (i.e. the same set-up as above),Introducing the two dimensionless parameterswe write the leading order result in the formwithThe generalised kernel K(P, Q) encodes the discretisation of the observable. Up to terms of it is given bywith the continuum kernel given byThe finite volume calculation for has taught us that the contributions of the flow have to cancel by themselves, i.e. a cancellation with other contributions by the action or the observable is not possible, due to the linear independence of the momentum sums. In order to assess the improvement of the Zeuthen flow it is therefore enough to focus on these contributions. Using again the Feynman gauge for flow and action, we obtain terms from the flow of the formand a second term with a similar structure. In both cases it is useful to note the property of the kernel thatwhere T(P) is the transverse projector,The correction to the Zeuthen flow kernel, , has the nice property thatHence we can conclude that the Zeuthen flow does not contribute any effects to this 2-point function either. Due to the different Lorentz index structure of this case compared to the simpler case of E(t, x), and to the fact that now, in general, the cutoff effects are functions of two variables (c, d), this test imposes further constraints on the possible improvement solutions. In particular, the so called chair flow in [27], which happens to also cancel the effects of in a finite volume, can be shown to produce contributions to the 2-point function considered here.

Conclusions and outlook

We have systematically investigated the structure of effects in flow quantities using Symanzik’s approach applied to the -dimensional local formulation of the theory. Improvement to for gradient flow quantities appears to be easier than one might have thought, mainly due to the classical nature of the gradient flow equation. In particular the classical a-expansion is sufficient to obtain the counterterms for both local composite operators at positive flow time and the action in the -dimensional bulk (i.e. due to the absence of loops in the bulk, no new counterterms are generated).

Our main results are summarised in the Zeuthen flow equation (5) and the improved lattice definitions of the observable E(t, x), either as linear combination of clover and plaquette definitions (54) or as the action density of the tree-level improved Lüscher–Weisz action. We have shown that the integration of this Zeuthen flow equation and the evaluation of classically improved observables do not produce any effects to any order in the coupling or, indeed, non-perturbatively. At this point it is important to remark that although the analysis has been performed in the context of pure gauge theories, due to the classical nature of the flow equation, the aforementioned results are still valid in QCD or if any number of fermions in any representations are coupled to our gauge field. In the particular case of the pure gauge theory the only effects originate either from the 4-dimensional lattice action or from the additional counterterm parameterised by in the modified initial condition (67). Tree-level improvement is achieved with the Lüscher–Weisz gauge action and .

We have explicitly checked that the proposed Zeuthen flow equation does not generate any contribution to tree level for a variety of gluonic observables (different observables in arbitrary volumes). In doing so, we have shown that other proposals of the literature to improve the gradient flow (i.e. the -shift in [21], the coefficients in [22] or the chair flow in [27]) in fact do produce effects in some of the considered observables. In this sense, these proposals only produce vanishing cutoff effects in some particular situations (i.e. in infinite volume), and this cancellation should be regarded as accidental, and not as improvement.

Our results can be extended in various directions. First, it appears straightforward to extend the classical a-expansion to further observables, for example the energy-momentum tensor. When considering n-point correlation functions of such observables with or if boundary conditions do not respect translation invariance in some directions (as is the case with SF and open boundary conditions), some additional work is required to also eliminate total derivative terms which may contribute at any order in a. We also note that the improvement of observables and the flow equation are conceptually separate from the effects at . It is therefore conceivable to push the expansion further, in order to also cancel terms at . It is not clear how complicated this would be for the flow equation, but it is certainly an option for observables. However, one should be aware that higher order improvement would typically render these observables less local in lattice units. Another natural generalisation would be the extension of our work to include fermions and the fermionic flow equation, introduced in Ref. [4].

Table 2

Kernels corresponding to different choices of discretisation, and discretisation effect corrections for some of the most popular choices. See Appendix A for any unexplained notation

Discretisation	K_μν(p;λ)
Plaquette	p^^2δ_μν-(1-λ)p^_μp^_ν
Lüscher–Weisz	p^^2δ_μν-(1-λ)p^_μp^_ν+a^212[(p^^4+p^^2p^μ2)δ_μν-p^_μp^_ν(p^μ2+p^ν2)]
Clover	p˚^2c^μ2δ_μν-p˚_μc^_μp˚_νc^_ν
Zeuthen	(1-a^2p^μ2/12)p^^2δ_μν-p^_μp^_ν+a^212[(p^^4+p^^2p^μ2)δ_μν-p^_μp^_ν(p^μ2+p^ν2)]+λp^_μp^_ν