Analysis of Ward identities in supersymmetric Yang-Mills theory

In numerical investigations of supersymmetric Yang-Mills theory on a lattice, the supersymmetric Ward identities are valuable for finding the critical value of the hopping parameter and for examining the size of supersymmetry breaking by the lattice discretisation. In this article we present an improved method for the numerical analysis of supersymmetric Ward identities, which takes into account the correlations between the various observables involved. We present the first complete analysis of supersymmetric Ward identities in $\mathcal{N}=1$ supersymmetric Yang-Mills theory with gauge group SU(3). The results show that lattice artefacts scale to zero as $O(a^2)$ towards the continuum limit in agreement with theoretical expectations.


Abstract:
In numerical investigations of supersymmetric Yang-Mills theory on a lattice, the supersymmetric Ward identities are valuable for finding the critical value of the hopping parameter and for examining the size of supersymmetry breaking by the lattice discretisation. In this article we present an improved method for the numerical analysis of supersymmetric Ward identities, which takes into account the correlations between the various observables involved. We present the first complete analysis of supersymmetric Ward identities in N = 1 supersymmetric Yang-Mills theory with gauge group SU (3). The results indicate that lattice artefacts scale to zero as O(a 2 ) towards the continuum limit in agreement with theoretical expectations.

Introduction
Ward identities are the key instruments for studying symmetries in quantum field theory. They represent the quantum counterparts to Noether's theorem, expressing the realisation of a classical symmetry at the quantum level in terms of relations between Green's functions. They also allow to characterise sources of explicit symmetry breaking. In the case of theories that are regularised non-perturbatively by means of a space-time lattice, Ward identities are a useful tool for the investigation of lattice artefacts, which are related to the breaking of symmetries. In lattice QCD, for example, chiral Ward identities in the form of the PCAC relation are being used to quantify the breaking of chiral symmetry by the lattice discretisation, and thereby to control the approach to the continuum limit [1].
For supersymmetric (SUSY) theories the corresponding relations are the supersymmetric Ward identities. In the context of numerical investigations of supersymmetric Yang-Mills theory on a lattice, SUSY Ward identities are being employed for a two-fold purpose [2]. First, in numerical simulations using Wilson fermions a gluino mass is introduced, which breaks supersymmetry softly. With the help of SUSY Ward identities the parameters of the model can be tuned such that an extrapolation to vanishing gluino mass is possible. Second, the discretisation on a lattice generically breaks supersymmetry [3], leading to lattice artefacts of order a in the lattice spacing. By means of SUSY Ward identities it can be checked if lattice artefacts are small enough for an extrapolation to the continuum limit.
Our collaboration has employed SUSY Ward identities in previous investigations of N = 1 supersymmetric Yang-Mills theory with gauge group SU(2); for recent result see [4]. In the analysis of SUSY Ward identities, following the methods introduced in [2], the correlations between the various quantities entering the calculation are, however, not being taken into account. Therefore, for our present studies with gauge group SU(3) we developed a method, based on a generalised least squares fit, that incorporates these correlations. In this article we describe the method and present the results of the first complete analysis of SUSY Ward identities for supersymmetric Yang-Mills theory with gauge group SU(3).

Supersymmetric Ward identities on the lattice
The N = 1 supersymmetric Yang-Mills (SYM) theory is the supersymmetric extension of Yang-Mills theory with gauge group SU(N c ). It represents the simplest field theory with supersymmetry and local gauge invariance. In the present investigations of our collaboration [5] we are focussing on gauge group SU (3). SYM theory describes the carriers of gauge interactions, the "gluons", together with their superpartners, the "gluinos", forming a massless vector supermultiplet. The gluons are represented by the non-Abelian gauge field A a µ (x), a = 1, . . . , N 2 c − 1. The gluinos are massless Majorana fermions, described by the gluino field λ a (x) obeying the Majorana conditionλ = λ T C with the charge conjugation matrix C, thus being their own antiparticles. Gluinos transform under the adjoint representation of the gauge group, so that the gauge covariant derivative is given by (D µ λ) a = ∂ µ λ a + g f abc A b µ λ c . In the Euclidean continuum the (on-shell) Lagrangian of the theory, where auxiliary fields have been integrated out, is where F a µν is the non-Abelian field strength. Adding a gluino mass term (m 0 /2)λ a λ a , which is necessary in view of the numerical simulations, breaks supersymmetry softly.
Infinitesimal supersymmetry transformations, that leave the action of the massless theory invariant, are given by where σ µν = (1/2)[γ µ , γ ν ], and the parameter ǫ is a Grassmann valued spinor. Noether's theorem, applied to the classical theory, yields a supercurrent [6] whose divergence is proportional to the gluino mass, where Both S µ (x) and χ(x) are spinorial quantities. The corresponding formal SUSY Ward identities in the quantised theory with a mass term are Here Q(y) is any suitable insertion operator, and the last term represents a contact term given by the SUSY variation of Q(y), which vanishes if Q(y) is localised at space-time points different from x. A quantised theory is, however, only properly defined once it is regularised. Regularisation on a lattice and renormalisation leads to significant modifications of the Ward identities [7,2]. For details we refer to the cited articles, and just report the main results. In addition to the soft breaking by the gluino mass term, supersymmetry is broken by the lattice regularisation. Analysis of the relevant operators indicates that a continuum limit should exist with the following characteristics. First, the gluino mass receives an additive renormalisation, leading to a subtracted gluino mass m S . Second, and more important, the supercurrent mixes with another dimension 7/2 current, namely Based on suitably defined SUSY transformations on the lattice [7,8], the resulting SUSY Ward identity, omitting contact terms, reads where Z S and Z T are renormalisation coefficients. A renormalised supercurrent can then be defined through S R µ = Z S S µ + Z T T µ . In our numerical simulations we use a lattice action proposed by Curci and Veneziano [7], which is built in analogy to the Wilson action of QCD for the gauge field and Wilson fermion action for the gluino. Both supersymmetry and chiral symmetry are broken on the lattice, but they are expected to be restored in the continuum limit if the gluino mass m S is tuned to zero. The Curci-Veneziano action for SYM theory on the lattice is given by is the gauge field action with inverse gauge coupling β = 2N c /g 2 , summed over the plaquettes p, and is the gauge field variable in the adjoint representation (T a are the generators of SU(N c )), and the hopping parameter κ is related to the bare gluino mass via κ = 1/(2m 0 + 8). In our numerical simulations the fermion action is O(a) improved by addition of the clover term with the one-loop coefficient specific for this model [9].
The supercurrent S µ (x) and the density χ(x) can be defined on the lattice in various ways, differing by O(a) terms. We choose the local transcriptions of the continuum forms, which have led to the best signals in previous numerical studies. For this choice, ∇ µ indicates the symmetric lattice derivative, and P (cl) ρν (x) is the clover plaquette. The supersymmetric continuum limit is obtained at vanishing gluino mass m S . The value of the critical hopping parameter κ c , where m S is zero, has to be determined numerically. With suitable choices of Q(y), this can be achieved with the lattice SUSY Ward identity. The expectation values appearing in Eq. (8) can be evaluated in the Monte Carlo calculations. This allows to obtain the coefficient m S /Z S , which in turn enables us to locate the point m S = 0. An alternative tuning is obtained from the signals of a restored chiral symmetry, see below. It is expected that both are consistent up to lattice artefacts. The investigation of the SUSY Ward identities allows to confirm this scenario and to estimate the relevant lattice artefacts.

Numerical analysis of SUSY Ward identities
In the numerical analysis it is convenient to project to zero momentum by summing the operators over the three spatial coordinates. As a result one obtains a Ward identity for each time slice separation t = x 4 − y 4 . Each term in Eq. (8) is a 4 × 4 matrix in Dirac space and can be expanded in the basis of 16 Dirac matrices. Using discrete symmetries one can show that only two non-trivial independent equations survive [2]: where O(a) terms are omitted, and Here, traces over spinorial indices are implied. Concerning the insertion operator, it turned out that gives the best signal. The signal-to-noise ratio is improved further by applying APE and Jacobi smearing to this operator.
The six different correlatorsx b,t,α are estimated numerically in our Monte Carlo simulations for gauge group SU(3). The usual estimators for these expectation values are the numerical averages of the corresponding observables over the Monte Carlo run. Let us call these averages x b,t,α . They are random variables with expectation valueŝ x b,t,α ≡ x b,t,α . It should be noted that only data at t ≥ 3 are being considered in order to avoid contamination by contact terms.
For each t the two equations (13) could be solved for Taking all t together, however, we have an overdetermined set of equations for these two coefficients. The aim is to find solutions for A and B numerically such that with the measured values x b,t,α the equations are satisfied approximately in an optimal way. In previous studies for gauge group SU(2) the coefficients A and B have been calculated by means of a minimal chi-squared method, as proposed in [2]. The correlators x b,t,α are, however, statistically correlated amongst each other, in particular for nearby values of t, and these correlations have not been taken into account.
In order to improve on this point, we have developed a method, which takes all correlations fully into account, so that more reliable results and error estimates can be obtained. The approach is based on the method of generalised least squares [10].
The equations (13) hold for the expectation values. With the notation and the double index i = (b, t), they can be written Let C iα,jβ = x iα x jβ − x iα x jβ be the covariance matrix of x iα . The probability distribution of the x iα is given by P ∼ exp(−L) with For estimating A α we employ the method of maximum likelihood in the following way.
1. For given x iα , consider A α to be fixed and determinex iα such that P is maximal under the constraint α A αxiα = 0. The value P max (A α ) at maximum depends on A α .

Find
Minimising L with the help of Lagrange multipliers gives and where For given A α the matrix D ij is estimated, up to an irrelevant constant factor, from the measured values by whereC iα,jβ is the covariance matrix of the primary observables. Now the minimum of L min (A α ) as a function of the parameters A 2 and A 3 (A 1 = 1) has to be found. Because D ij depends on the A α , it is not possible to do this analytically, and we determine the global minimum numerically, thus obtaining A 2 and A 3 . To get the statistical errors we re-sample the data and apply the jackknife method, repeating the whole procedure for each jackknife sample. In this way we arrive at our final result for B = am S Z −1 S .

Results for SU(3) SYM
For SYM theory with gauge group SU(3) we have applied the method to our current simulation ensembles obtained with O(a) improved clover fermion action [11] at different inverse gauge couplings β and hopping parameters κ. At two lattice spacings, corresponding to β = 5.4 and 5.5, the available statistics has allowed to obtain reliable results for the Ward identities. From the results for the gluino mass parameter am S Z −1 S the value of κ c , where m S vanishes, can be estimated.
Comparing the results for am S Z −1 S with those from the earlier method, which does not properly take the correlations into account, we find that the values are compatible within errors, but this time we have a precise and reliable estimate of the errors. As examples, the results of both methods for β = 5.5 are shown in Tab Table 1: Results for am S Z −1 S from the previous method and from the generalised least squares (GLS) method for our ensembles at β = 5.5.
An alternative way to estimate κ c in the Monte Carlo calculations employs the mass of the adjoint pion a-π, see e. g. [12]. The a-π is an unphysical particle in SYM theory. However, by arguments based on the OZI-approximation [13], and in the framework of partially quenched chiral perturbation theory [14], the squared mass m 2 a-π is expected to vanish linearly with the gluino mass close to the chiral limit.
In Fig. 1 we show am S Z −1 S and (am a-π ) 2 as a function of 1/(2κ) for our two values of β. Both quantities depend linearly on κ −1 within errors, as expected, and yield independent estimates of the value of κ c .
The values of κ c obtained from the Ward identities and from m 2 a-π are very close to each other, but there is a small difference. This discrepancy should be due to lattice artefacts, and we expect it to disappear in the continuum limit.
In the case of lattice QCD, Wilson chiral perturbation theory to leading order shows a shift linear in a in the dependence of the squared pion mass on the quark mass: with certain low-energy constants B 0 and W 0 [15,16]. On the other hand, for the PCAC quark mass, defined by means of the chiral Ward identity, exactly the same shift is present in leading order [17], Consequently, at vanishing pion mass, the remnant m PCAC is of order a 2 , and this result is not changed in higher orders of chiral perturbation theory, In SYM the adjoint pion mass can be calculated in partially quenched chiral perturbation theory [14]. We haven't evaluated the contributions from the lattice terms explicitly, but the structure of terms is similar to those for QCD, and therefore we expect that in SYM the remnant gluino mass ∆m S at vanishing adjoint pion mass is of order a 2 , too. In order to check this numerically, the masses have to be expressed in a physical scale. We use the scale w 0 , defined through the gradient flow; for details see [11]. In Fig. 2 we show the remnant gluino mass as a function of the squared lattice spacing a 2 . The line through the points extrapolates to zero within errors. For an analogous plot linear in a this is by far not the case. Having only two points available, one has to be cautious drawing conclusions, but the result clearly indicates that the remnant gluino mass ∆m S vanishes proportional to a 2 in the continuum limit.

Conclusions
We have presented a method for the numerical analysis of SUSY Ward identities in supersymmetric Yang-Mills theory on a lattice, which employs the expectation values of the relevant operators on a range of time slices. The statistical correlations between all observables are taken into account by means of a generalised least squares procedure. Applied to SUSY Yang-Mills theory with gauge group SU(3), the value of the hopping parameter, where the renormalised gluino mass vanishes, can be estimated, and is in rough agreement with the estimation using the adjoint pion mass. The difference between the estimates appears to vanish in the continuum limit. Our results represent the first continuum extrapolation of SUSY Ward identities. The scaling of lattice artefacts as of O(a 2 ) is in agreement with theoretical expectations.