Continuum extrapolation of Ward identities in N=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathcal {N}}=1}$$\end{document} supersymmetric SU(3) Yang–Mills theory

In N=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {N}}=1$$\end{document} supersymmetric Yang–Mills theory, regularised on a space-time lattice, in addition to the breaking by the gluino mass term, supersymmetry is broken explicitly by the lattice regulator. In addition to the parameter tuning in the theory, the supersymmetric Ward identities can be used as a tool to investigate lattice artefacts as well as to check whether supersymmetry can be recovered in the chiral and continuum limits. In this paper we present the numerical results of an analysis of the supersymmetric Ward identities for our available gauge ensembles at different values of the inverse gauge coupling β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} and of the hopping parameter κ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}. The results clearly indicate that the lattice artefacts vanish in the continuum limit, confirming the restoration of supersymmetry.


Introduction
Supersymmetry (SUSY) is an elegant idea which relates fermions and bosons, whose spin differs by 1/2, through supercharges [1]. SUSY provides dark matter candidates, arising from the lightest supersymmetric particles [2]. In addition to that, supersymmetric extensions of the Standard Model would resolve the hierarchy problem [3]. N = 1 supersymmetric Yang-Mills (SYM) theory, which is being a e-mail: sajid.ali@uni-muenster.de (corresponding author) b e-mail: georg.bergner@uni-jena.de c e-mail: henning.gerber@posteo.de d e-mail: montvay@mail.desy.de e e-mail: munsteg@uni-muenster.de f e-mail: stefano.piemonte@ur.de g e-mail: scior@physik.uni-bielefeld.de considered in this article, provides an extension of the pure gluonic part of the Standard Model [4]. It describes the strong interactions between gluons and gluinos, the superpartners of the gluons. Gluinos are Majorana particles that transform under the adjoint representation of the gauge group. The onshell Lagrangian of N = 1 SYM theory, which consists of the gluon fields A a μ (x) and the gluino fields λ a (x), where a = 1, . . . , N 2 c − 1, can be written in Minkowski space as where the first term, containing the field strength tensor F a μν , is the gauge part, and D μ in the second term is the covariant derivative in the adjoint representation of the gauge group SU(N c ), N c being the number of colors. The last part of the above Lagrangian is a gluino mass term which breaks SUSY softly for mg = 0, which means that it does not affect the renormalisation properties of the theory and that the spectrum of the theory depends on the gluino mass in a continuous way. The physical spectrum of this theory is expected to consist of bound states of gluons and gluinos, arranged in mass degenerate supermultiplets if SUSY is not broken [5,6].
In order to perform Monte-Carlo simulations of the theory, we discretise the Euclidean action and put it onto a fourdimensional hypercubic lattice. We use the Curci-Veneziano version [7] of the lattice action S = S g + S f , where the gauge part S g is defined by the usual plaquette action with the inverse gauge coupling given by β = 2N c /g 2 , and the fermionic part implements the gluinos as Wilson fermions. Here the adjoint link variables are defined by V ab,xμ = 2 tr (U † xμ T a U xμ T b ), where T a are the generators of the gauge group, and the hopping parameter κ is related to the bare gluino mass mg by κ = 1/(2mg + 8). In order to approach the limit of vanishing gluino mass, the hopping parameter has to be tuned properly. In our numerical investigations the fermionic part is additionally O(a) improved by adding the clover term −(c sw /4)λ(x)σ μν F μν λ(x) [8].
In our previous investigations we have determined the low-lying mass spectrum of the theory with gauge group SU(2) and SU(3) non-perturbatively from first principles using Monte Carlo techniques [4,[9][10][11], and obtained mass degenerate supermultiplets [12].

SUSY Ward identities
In classical physics, Noether's theorem provides a relation between symmetries and conservation laws. In the case of quantum field theories, symmetries are translated to Ward identities, representing quantum versions of Noether's theorem. In N = 1 supersymmetric Yang-Mills theory a gluino mass term breaks SUSY softly. The soft breaking effects vanish in the chiral limit, a limit where theory is characterised by massless gluinos. In order to analyse this breaking of supersymmetry and to identify the chiral limit, we employ the Ward identities for supersymmetry. Moreover, on the lattice supersymmetry is broken explicitly due to the introduction of the discretisation of space-time lattice as a regulator of the theory. SUSY Ward identities can be used to check whether supersymmetry is restored in the continuum limit.
In the Euclidean continuum, on-shell supersymmetry transformations of the gauge and gluino fields are given by where the transformation parameter ε is an anticommuting Majorana spinor. From the variation of the action under a supersymmetry transformation with a space-time-dependent parameter ε(x) one derives the SUSY Ward identities. For any suitable gauge invariant local operator Q(y), they read where S μ (x) = (S α μ (x)) is the supercurrent of spin 3/2, and the term mg χ(x)Q(y) is due to the gluino mass in the action of the theory. In the continuum the supercurrent S μ (x) and the operator χ(x) are given by The last term of Eq. (5) is a contact term, which contributes only if x = y, and it can be avoided if Q(y) is not localised at x. Therefore the contact term is ignored in the following discussions.
The four-dimensional space-time lattice breaks SUSY explicitly. As a consequence, the lattice versions of the Ward identities differ from their continuum counter parts by an additional term X S (x)Q(y) . The explicit form of this term is known, but need not be displayed here. At tree level this term is proportional to the lattice spacing a and vanishes in the limit of zero lattice spacing. At higher orders in perturbation theory, nevertheless, the contribution of this term is finite in the continuum limit due to divergences proportional to 1/a that multiply the factor a. This plays a role for the renormalisation of the supercurrent and of the gluino mass [7,13]. In the renormalisation of SUSY Ward identities, operators of dimensions ≤ 11/2 have to be taken into account. They lead to a modification of the gluino mass, and in addition a current T μ , mixing with the supercurrent, appears, corresponding to an operator of dimension 9/2. Consequently, on the lattice the following Ward identities are obtained where Z S and Z T are renormalisation coefficients. The subtracted gluino mass is defined as m S = mg −m, wherem is the mass subtraction coming from the operators of dimension 7/2. The mixing current is defined as Regarding the local insertion operator Q(y), our choice is the spinor Q(y) = χ (sp) (y), with where the indices i, j ∈ {1, 2, 3}. The reason behind this choice is that it gives the best signal [13].

Numerical analysis of SUSY Ward identities
We have analysed the SUSY Ward identities numerically, employing the configurations produced in our project on N = 1 supersymmetric Yang-Mills theory with gauge group SU(3). Numerically it is convenient to use integrated Ward identities where integration or sum is performed over all three spatial coordinates. The resulting identities will then hold for every time-slice distance t. In the analysis the data from all time-slice distances in an interval t min ≤ t ≤ t max are included. The lower limit t min is always taken to be larger or equal than 3 in order to avoid contamination from contact terms. The choice of t min for the different ensembles of configurations is discussed below. Since the correlation functions are symmetric or antisymmetric in t, the upper limit t max is chosen to be half of the time extent of the lattice. Each term in Eq. (8) is a 4× 4 matrix in spin-space and can be expanded in the basis of 16 Dirac matrices, i. e. 1, γ 5 , γ μ , γ μ γ 5 , iσ μν . It can be shown, with the help of discrete symmetries, that only the following two contributions are non-zero [13]: where In In these equations thex iα = x iα are the expectation values of random variables x iα , which themselves are considered to be the results of a finite Markov chain. We compute the estimators x iα for the correlation functionsx iα numerically using high performance facilities. The Eq. (13), including all time-slice distances t from t min to t max , are solved simultaneously for A α by means of minimal chi-squared methods. Two methods, namely the so-called local method and global method, have been used in the past by our collaboration [4,13]. These methods, however, do not take properly into account correlations between the different quantities appearing in Eq. (13). For this purpose we have developed a new method based on a generalised least squares fit, the so-called GLS method [14], based on the maximum likelihood. For fixed A α (α = 1, 2, 3) and given numerical data x iα , the probability distribution P ∼ exp(−L) of the quantitiesx iα , subject to the constraints (13), has its maximum at a point where L = L min , with where Next, the desired coefficients A α have to be found such that L min as a function of A 2 and A 3 is minimised. This cannot be solved analytically, and we find A α numerically such that the global minimum of L min (A 2 , A 3 ) is reached; for details see Ref. [15]. In particular, owing to A 3 = −am S Z −1 S this provides us with the subtracted gluino mass m S up to the renormalisation factor. To estimate the statistical uncertainties we employ the standard Jackknife procedure.

Discretisation effects
All terms in the Ward identity (8), including the O(a) term X S (x)Q(y) , are correlation functions of gauge invariant operators. In the corresponding Eq. (11) they are correlation functions of operators localised on time slices or pairs of adjacent time slices at distance t. As for any gauge invariant correlation function of this type, they decay exponentially in t, with a decay rate given by the mass gap of the theory. For very small t the contributions of higher masses will affect the impact of the O(a) term on the Ward identities. Therefore we expect that the value of the obtained gluino mass will depend on the minimal time slice distance t min . This effect should become negligible at sufficiently large t min . On the other hand, if t min is chosen too large, noise in the data will dominate. The behaviour that can be observed in Fig. 1 is compatible with these expectations.
An adequate choice of t min is therefore important for the quality of the results. We cope with this in two ways.
In order to avoid perturbing effects at too small t min and a poor signal-to-noise ratio at too large t min , for each hopping parameter and inverse gauge coupling, the value of t min is selected by finding an optimal starting point where a plateau in the subtracted gluino mass begins. The results are presented in Table 1.
In the second approach, we consider that our simulations of the theory are done at different values of the lattice spacing a, which leads to different O(a) terms in the Ward identities. A fixed value of t min in lattice units would mean a lower limit on the time-slice distances in physical units, that is on the cutoff-scale and shrinks to zero in the continuum limit. Instead it would be more appropriate to consider t min at constant physical distance for all gauge ensembles. This is done in the following way. At the coarsest lattice spacing, at inverse gauge coupling β 0 , the value of t min is selected according to the plateau criterion explained above. For finer lattice spacings at inverse gauge couplings β i the corresponding t min are then obtained by scaling with a physical scale. In order to determine the physical scale we use the mass m gg of the gluino-glue parti-cle and the Wilson flow parameter w 0 . Correspondingly, t min is scaled according to where β 0 = 5.4, β 1 = 5.45, β 2 = 5.5, and β 3 = 5.6. The resulting t min is rounded to the nearest integer value. The values obtained by this method are collected in Table 2. In most points they are equal or almost equal to those in Table 1. 3.2 Adjoint pion and remnant gluino mass The chiral limit is defined by the vanishing of the subtracted gluino mass. Its measured values can therefore be employed for the tuning of the hopping parameter κ to the chiral limit.
On the other hand, we can also use the vanishing of the adjoint pion mass m a-π for the tuning [16]. The adjoint pion a-π is an unphysical particle in the SYM theory, that can be defined in partially quenched chiral perturbation theory [17]. In the numerical simulations its correlation function can be computed as the connected piece of the correlation function of the a-η particle. Similar to the Gell-Mann-Oakes-Renner relation of QCD [5], in the continuum limit there is a linear relation between the adjoint pion mass squared and the gluino mass: m 2 a-π ∝ mg. The numerical results for the subtracted gluino mass from the Ward identities and the adjoint pion mass squared in lattice units are shown for β = 5.6 in Fig. 2 together with their extrapolations towards the chiral limit.
In the continuum the subtracted gluino mass and the adjoint pion mass should vanish at the same point. On the lattice, however, this is not the case due to lattice artefacts. As an estimate for this discrepancy we determine the value of the subtracted gluino mass at vanishing adjoint pion mass. This quantity is called the remnant gluino mass (am S Z −1 S ), and it is expected to vanish in the continuum limit. The values of the remnant gluino mass, obtained by taking an average of the values calculated using the procedures explained above, are presented in Table 3.

Continuum limit
The remnant gluino mass is a lattice artefact and should vanish in the continuum limit a → 0. It is therefore a quantity to check on whether supersymmetry is recovered or not. Concerning the dependence of the remnant gluino mass on the lattice spacing, arguments based on partially quenched chiral perturbation theory suggest that the remnant gluino mass is of order a 2 at m 2 a-π = 0 [13]. In order to investigate this relation, the remnant gluino mass has to be expressed in physical units. Our choice for the scale is the Wilson flow parameter w 0 , which is defined through the gradient flow [10]. We use its values extrapolated to the chiral limit, w 0,χ . Similarly the lattice spacing is represented by a/w 0,χ . Our numerical results for the remnant gluino mass as a function of the lattice spacing and its extrapolation towards the continuum limit are shown in Fig. 3. The data points in Fig. 3a show the results from separate chiral extrapolations for each   Fig. 3 The remnant gluino mass (w 0 m S Z −1 S ) in physical units w 0 as a function of the lattice spacing squared, and its linear extrapolation towards the continuum limit lattice spacing and the corresponding extrapolation to the continuum limit. The extrapolation to the continuum and the error of this extrapolation are obtained by means of parametric bootstrap with linear fits. On the other hand, Fig. 3b is obtained by means of a simultaneous fit of the dependence on the hopping parameter and the lattice spacing [18].
The remnant gluino mass in the continuum limit is compatible with zero within one standard-deviation, confirming the preliminary results present in Ref. [15] with only two data points. Lattice artefacts vanish in the continuum limit as expected, and supersymmetry is recovered in the chiral and continuum limits, in agreement with our findings from the mass spectrum [12].

Conclusion
In this paper we have presented numerical results of an analysis of SUSY Ward identities in N = 1 supersymmetric Yang-Mills theory on the lattice with gauge group SU(3). Contact terms and O(a) lattice artefacts in the Ward identities have been controlled by suitable choices of time-slice distances. Ensembles of gauge configurations at four different values of the lattice spacing and various hopping parameters have been analysed, allowing us for the first time to perform an extrapolation to the continuum limit, where the lattice artefacts vanish. The remnant gluino mass has been extrapolated in two alternative ways, on the one hand by extrapolating to the chiral limit at each lattice spacing separately and then to the continuum limit, and on the other hand by means of a simultaneous extrapolation to the chiral and continuum limit. With both extrapolations the lattice artefacts in the subtracted gluino mass appear to scale to zero as of order a 2 in agreement with the theoretical expectations. Our findings support the validity of SUSY Ward identities and the restoration of supersymmetry in the continuum limit.