A study of center and chiral symmetry realization in thermal super Yang-Mills theory using the gradient flow
The realization of center and chiral symmetries in super Yang-Mills theory (SYM) is investigated on a four-dimensional Euclidean lattice by means of Monte Carlo methods. At zero temperature this theory is expected to confine external fundamental charges and to have a non-vanishing gaugino condensate, which breaks the non-anomalous Z chiral symmetry. In previous studies at finite temperatures, the phase transitions corresponding to deconfinement and chiral restoration were observed to occur at roughly the same critical temperature for SU(2) gauge group. We find further evidences for this observation from new measurements at smaller lattice spacings using the fermion gradient flow, and we discuss the agreement of our findings with conjectures from superstring theory. The implementation of the gradient flow technique allows us also to estimate, for the first time, the condensate at zero temperatures and zero gaugino mass with Wilson fermions.
Supersymmetric Yang-Mills theories (SYM) have been a useful laboratory for the understanding of the fundamental properties of quantum field theories. For instance, the study of supersymmetry combined with gauge symmetry has led to the discovery that duality is a feature of certain extended SYM. Besides the well-known gauge-gravity duality between the conformal SYM and string theory on the curved-space , the electro-magnetic duality has been conjectured for SYM and supersymmetric QCD [1, 2, 3]. The most interesting open challenge is the search for a dual theory of QCD that would open the possibility for the analytical understanding of confinement and strong interactions. SYM is the supersymmetric extension of the pure gauge sector of QCD describing the strong interactions of gluons and gauginos, and it shares many interesting aspects with QCD while retaining the minimal supersymmetry. There have been several theoretical predictions that this theory should exhibit confinement, chiral symmetry breaking and mass-gap [4, 5, 6, 7]. SYM is therefore a perfect candidate for the analytical understanding of non-pertubative phenomena of non-Abelian gauge theories, like confinement and chiral symmetry breaking.
Lattice simulations are the natural tool to explore non-perturbatively the phase diagram of strongly interacting gauge theories, such as SYM. The study of chiral symmetry breaking using the lattice discretization requires however a certain care. The naive discretization of the fermion action represents effectively sixteen fermion-species, the so-called doublers, instead of only one fermion. As a consequence, the fermion-boson state degeneracy and supersymmetry are broken. The Wilson term removes the doublers from the physical spectrum, at the cost of breaking chiral symmetry explicitly. The chiral condensate, the order parameter for chiral symmetry breaking, acquires a divergent additive renormalization term. The study of chiral symmetry breaking of SYM in the Wilson formulation is therefore challenging. Recently the gradient flow has been proposed as a regularization-scheme independent smoothing technique that is able to simplify drastically the renormalization of lattice bare composite operators [8, 9, 10]. In particular, the flowed chiral condensate is free from the additive renormalization of Wilson fermions .
In recent years the gradient flow has found a growing spectrum of applications. In particular the investigations of finite temperature QCD including the chiral transition have shown the benefits of this method [11, 12]. Furthermore, a novel approach to compute operator dimensions in conformal field theories was recently proposed in , which exploits the relation between the gradient flow and renormalization group transformations. The application for supersymmetric theories has been suggested in [14, 15], where it can help to renormalize the supercurrent. In several works also a supersymmetric version of the method has been developed [16, 17], which might even avoid the necessity of a multiplicative renormalization of the fermions.
In this contribution we present an extended study of the phase diagram of SYM with the gauge group SU(2) at zero and non-zero temperature. We measure the chiral condensate expectation value at positive gradient flow time. We show strong evidences that chiral symmetry is spontaneously broken at zero temperature by a non-vanishing expectation value of the gaugino condensate. We show that at high temperature chiral symmetry is restored and that a phase transition occurs in the massless limit. We also show that there are only two phases in the mass-less limit, characterized by both chiral symmetry breaking and confinement at low temperature and chiral symmetry restoration and deconfinement at high temperature. We have not found any evidences of mixed phases where deconfinement occurs while chiral symmetry is broken. A possible interpretation of our results from the point of view of superstring/M-theory is presented in the last section.
2 Super Yang-Mills theory and its lattice discretization
Supersymmetry is, according to the Coleman-Mandula theorem , the only possible non-trivial extension of the space-time symmetries. It extends the Poincaré algebra by irreducible spinorial generators. As the spinor charges act on the irreducible representations of the Poincare group by changing the spin by , the supermultiplets contain fermionic and bosonic fields. The simplest four-dimensional supersymmetric Yang-Mills theory corresponds to a theory with one super-multiplet, consisting on-shell and in Wess-Zumino gauge of one vector gauge field and one Majorana spinor, the so-called super Yang-Mills theory. In the present study we focus on the gauge group SU(2). For the Euclidean metric, the four-dimensional on-shell action can be written as
where is the usual Yang-Mills (YM) field strength, a Majorana spinor, the gaugino, and the gauge covariant derivative. The mass term breaks supersymmetry softly. Here all fields are su(N)-algebra valued and thus transform in its adjoint representation, i.e. for SU(2). The third term is a topological term, whose space-time integration yields the winding number of the gauge field. An important role is played by in the realization of chiral symmetry in the quantum theory, as will be discussed in the next section.
SYM is evidently very similar to QCD, as it both lacks of scalars and includes spinor fields interacting with a gauge field. This characteristic makes the model very attractive, if one intends to get some deeper understanding of non-perturbative phenomena of QCD as chiral symmetry breaking and confinement. In supersymmetric theories these phenomena are usually more accessible due to the constrains that the supersymmetry imposes.
Observables of the theory can be numerically computed by Monte-Carlo methods and based on a discretization of the Euclidean space-time, i.e. a lattice. There, the spinors are assigned to the sites and the Yang-Mills field to the links through the parallel transporter in the fundamental representation, being the lattice spacing. On the lattice, the SYM action takes on the form 
where are the parallel transporters in the adjoint representation, defined as 111 are the generators in the fundamental representation.. Here is proportional to the inverse squared gauge coupling, is the plaquette, i.e. the minimal closed loop of parallel transporters and is the tree-level Symanzik improvement term. One-level stout smearing with parameter has been employed for the links in the fermion action. Finally, , the hopping parameter, is equal to . A finite gaugino mass breaks supersymmetry softly and the bare parameter must be tuned to the point of a vanishing physical mass. A numerical accurate and inexpensive way to achieve that is by means of the adjoint-pion a-. Although this is not a physical degree of freedom of the theory, its mass is related to the gaugino mass as , as was shown in  within the frame of partially quenched perturbation theory.
In general, the lattice discretization breaks super Poincaré symmetry down to some subgroup of the hypercube isometries. In the massless and continuum limit, however, supersymmetry is recovered, as confirmed through analysis of the supersymmetry Ward identities in . In thermal quantum field theory, temperature corresponds to the inverse radius of a compactificatied direction of the path integral with thermal boundary conditions, i. e. periodic for bosons and antiperiodic for fermions.
3 The phase diagram of Super Yang-Mills theory
3.1 Confinement of static fundamental charges
At zero temperature, the non-supersymmetric YM vacuum is expected to behave as a confining medium for external static color-electric charges. The vacuum can be probed through the Polyakov loop in the fundamental representation of the gauge group, which is the path ordered product of the links in the fundamental representation along a line which wraps in the compact direction
where denotes the three-dimensional lattice volume of the non-compactified directions. It is not invariant under center symmetry transformations. Thus, a vanishing expectation value of the fundamental Polyakov loop implies unbroken center symmetry. The Polyakov loop can be related to the exponential of the free energy of an isolated static fundamental quark. Hence, a vanishing vacuum expectation value of the loop also indicates that isolated fundamental quarks are states with infinite free energy. As the temperature is increased, the vacuum expectation value of should become non-vanishing at some critical temperature, where center symmetry would be broken spontaneously. In the phase of broken center symmetry the energy of an isolated quark would be finite and, in that sense, one would speak of quark deconfinement.
In the case of QCD, dynamical fundamental quarks break center symmetry explicitly and there is no real confinement-deconfinement phase transition but a crossover. In the confined phase, the quark-antiquark potential grows linearly until it is screened by another quark-anti-quark pair popping up from the vacuum. In the deconfined phase, asymptotic freedom leads to a so-called quark-gluon plasma. The identification of the phase transition is much clearer for SYM. The N-ality of the fermion representation implies that the fermion fields are singlet under center transformations. Thus, the action of the theory remains invariant with respect to center transformations even for massless gauginos and connection of confinement and center symmetry breaking is still valid. A critical behavior is expected at some finite temperature and the Polyakov loop represents a good order parameter for the deconfinement phase transition for SYM.
3.2 Chiral symmetry
supersymmetric Yang-Mills theory is invariant under certain chiral transformations when no soft SUSY-breaking mass term is included in the action. Chiral symmetry coincides with the U(1) R-symmetry of the theory, which leaves the SUSY algebra invariant222In this section is a Weyl spinor, as this notation is more natural for the massless theory.
This symmetry is, however, broken at the quantum level by instanton contributions. In the presence of a non-vanishing -term in the action (1), the U(1) chiral rotation is equivalent to
where is the number of colors. The path integral is therefore invariant only for , i.e. the chiral symmetry of the quantum theory is actually the discrete subgroup ZU(1). An interesting question is whether the quantum chiral symmetry is spontaneously broken by a non-vanishing fermion condensate. Indeed, it can be shown analytically that the answer is affirmative [22, 19, 23, 24, 25, 26, 27, 4]. Thanks to remarkable properties of supersymmetric theories like non-renormalization and holomorphicity of the effective superpotential with respect to fields and couplings, it was proven [28, 29] that if the effective theory at long distances is massive with color-singlets as degrees of freedom, the gaugino condensate takes on the form
As the gaugino condensate is not invariant with respect to the full Z group, i.e. only if , the vacuum of the theory does not exhibit the same symmetry as the partition function and the global Z symmetry is spontaneously broken down to its discrete subgroup Z, i.e. to the sign flip . Hence, chiral symmetry is spontaneously broken at zero temperature with degenerate vacua, which are connected by Z transformations. As the broken symmetry is a discrete one, the existence of domain walls interpolating between the different vacua is expected. In  it was found that domain walls in SYM are BPS-saturated states and an expression of their energy density was computed.
3.3 The phase diagram
A non-trivial task is to determine what happens to the vacuum of the theory and specifically to the condensate and Polyakov loop distributions at finite temperature. At some critical temperature it is expected that the condensate vanishes and thus, that the vacuum becomes Z-symmetric. In addition, center symmetry breaks spontaneously at the critical deconfinement temperature and the Polyakov loop should acquire a non-vanishing expectation value. An interesting question arising from this analysis is whether both critical temperatures coincide, i.e. if confinement and chiral symmetry breaking share the same non-perturbative origin. This question is far from trivial. Although QCD, for example, is also characterized by confinement and broken chiral symmetry at low energies, some theories have been found to exhibit only one of those . In some numerical studies of QCD it has been found that chiral restoration and deconfinement phase transition occur at nearly the same critical temperature. However, taking QCD as a starting point implies the difficulty that these phase transitions are actually crossovers due to the finite-massive fundamental color charges. This means that no exact order parameter can be found, which unambiguously signals the pseudo-critical temperature. Hence, this observation strongly depends on the order parameters chosen. Fortunately, both the Polyakov loop and the gaugino condensate are exact order parameters for the respective phase transitions in SYM. Consequently, the question of the relation between the chiral and deconfinement transition can be answered by studying these quantities in a systematic way, the only remaining complications being related to the regularization and the renormalization of the these observables. The lattice realization with Wilson fermions on the lattice breaks chiral symmetry, which implies besides the multiplicative renormalization factor () an additional additive renormalization () of the gaugino condensate
In previous investigations , the additive constant was removed by subtracting the bare condensate at zero temperature,
One downside of this approach is the fact that the subtracted condensate is fixed to vanish at , or some very low temperature. Although this subtraction should preserve the behavior of the order parameter near the critical temperature, it is not possible to determine the renormalized condensate at zero temperature. Hence, the picture of the realization of chiral symmetry in thermal SYM is, following this method, incomplete. As it will be clear in the next session, the additive renormalization is not necessary when the gradient flow is used, allowing the computation of the condensate at zero temperature and, moreover, making the results comparable to fermion discretizations which satisfy the Ginsparg-Wilson relation.
4 The gradient flow
4.1 Flow equations
Motivated in the context of trivialising maps , Lüscher studied the correlation functions of fields flowed through the equations  333The term is a gauge parameter, which is included for mere technical reasons.
where (resp. ) and are the Yang-Mills gauge field and field strength, respectively. and are spinor fields and is the gauge-covariant derivative in the adjoint representation.
The parameter describes a flow on the vector space of gauge fields. These equations have a smoothening effect on the fields, which are Gaussian-like smeared over an effective radius , as it can be easily seen by integrating equation (3) in the non-interacting limit ,
where is the number of space-time dimensions and the heat-kernel. Moreover, the term regularizes the integral in momentum space when , removing the UV divergences at large momenta. Some years ago it was shown that through this kind of gradient flow the smearing property remains at all orders in perturbation theory . In this context, Lüscher developed a D+1 dimensional quantum field theory, where is taken as an spurious extra Euclidean dimension.444 The unflowed theory lives on the D-dimensional boundary of a D+1 dimensional manifold, where fields propagate through the kernel of the flow equations. Fields would propagate along the new dimension through the heat kernel, which is a retarded propagator (there are no flow loops) and the flow equations are imposed through Lagrange-multiplier-fields. It was found that in this theory the BRS-Ward identities hold only if there are no counter-terms beyond those of the unflowed theory. This led to the observation that correlation functions of monomials of flowed gauge-fields are renormalized without the necessity of extra counter-terms and that spinor operators renormalize multiplicatively . 555The monomials renormalize according to the field content since the flowed fields are effectively non-local. Furthermore, a specific advantage of the gradient flow relies on the fact that it is regularization-scheme independent and all these results are expected to hold also on a space-time lattice. Therefore, currents and densities which are explicitly broken by the lattice discretization should be more easily accessible within this method. As a special case, the additive renormalization constant, necessary for the computation of the gaugino condensate with Wilson fermions, is not required if these are flowed up to some finite flow-value. Consequently, it is possible to study if at even with Wilson fermions. The value of is irrelevant if the bare lattice gauge coupling is fixed. Since usual lattices are Euclidean, non-zero temperatures are achieved just by compactifying and changing the number of lattice sites in temporal direction, where fermion fields fulfill thermal (anti-periodic) boundary conditions.
4.2 Computation of the chiral condensate from the gradient flow
The flowed bare gaugino condensate can be defined as 
Hence, computing the flowed condensate amounts to the action of the heat-kernel on the fermion propagator . Analogously, its discrete version can be straightforwardly written:
where is the Wilson-Dirac operator and the trace runs over space-time (spinor) and gauge group indices. Following , the trace is estimated stochastically by inserting a complete set of random complex vectors with and :
Here denotes the average with respect to to both the Monte-Carlo time and any internal group-symmetry representations. Finally, to compute the new vectors , the so-called adjoint flow equation
must be integrated from to , i.e. backwards in comparison with the flow equations presented above. Here, the gauge connection in the covariant four-dimensional Laplacian is flowed up to the same as .
up to some by means of the following Runge-Kutta integrator with step-size
where S is the Wilson plaquette action. The intermediate fields are then kept on the computer memory in order to access them during the integration of the adjoint fermion equation. As a second step one random source vector is generated on the lattice and then integrated by means of the Runge-Kutta integrator down to
Further, the vector is computed, e.g. through conjugate gradient. Subsequently it is contracted with and the result averaged over the lattice sites, i.e. one calculates .
5.1 Gaugino condensate at zero temperature
As anticipated, the remarkable properties of the gradient flow method allow for the computation of the gaugino condensate at zero temperature, since no ambiguity is introduced by the additive renormalization. We have taken for ensembles from our previous investigations of SU(2) SYM at on a lattice with Hopping-parameter values
The data of these four points were then used to make an extrapolation to the chiral point, i.e. to vanishing renormalized gaugino mass. As shown in Fig. 1, we obtain a finite value of the condensate in the chiral limit. For the numerical simulations with Wilson fermions, this is the first clear indication that the discrete chiral symmetry is spontaneously broken at zero temperature and zero gaugino mass. In previous studies without the method of the gradient flow, only an indirect observation of the non-zero condensate has been possible from a double peak of the histogram that appears in rather unstable simulations close to the chiral point [34, 35]. Our result is remarkably compatible with the results obtained through Domain-Wall and overlap fermions in [36, 37, 38, 39, 40]. However, for a direct comparison of the results a consistent renormalization of the condensate obtained with the Domain-Wall fermions would be required. Although the existence of a non-vanishing condensate is a well-known fact, being able to compute it from Wilson fermions without the complication of a residual mass, shows the goodness of the gradient flow method. Moreover it opens up the possibility to study the different phases at zero temperature without the computational more expensive Ginsparg-Wilson fermions.
5.2 Discrete chiral symmetry restoration and quark deconfinement
The second main purpose of this study is to investigate the realization of the chiral and center symmetries in SU(2) SYM as temperature is turned on. Some first results in this direction were obtained in . This problem now revisited using the gradient flow on four different new lattice ensembles at and Hopping parameter . To further cross-check the validity of our results, we have also analyzed a set of ensembles generated at and using the tree-level clover improved fermion action with unsmeared links. Finite temperatures were achieved by fixing the lattice parameters and compactifying one dimension on a circle, imposing thermal boundary conditions on the fields, i.e. anti(periodic) for fermions (bosons). Thus, the lattice size was set to for , with the upper bound parametrizing the zero-temperature limit. Setting the scale, i.e. going from lattice to physical units, was achieved through the parameter, which is defined through the gradient flow as 666Here a mass-dependent scheme was chosen, i.e. a different for each value.
with the field energy density and the flowed Yang-Mills field strength.
The results are summarized in figure 2 and 3. The condensate is in general considerably reduced as the temperature increases, but a clear phase transition is difficult to identify for the larger gaugino masses corresponding to smaller values of . This smooth behavior of the chiral condensate at larger masses is expected since the jump of the order parameter around the pseudo-critical temperature is less pronounced for a larger explicit symmetry breaking. For the smallest gaugino mass () the signal is considerably better and a jump at the critical temperature can be identified. This suggests, as expected, that the chiral restoration becomes indeed a true phase transition and the gaugino condensate an adequate order parameter in the chiral (and supersymmetric) limit. This can be more clearly seen from the disconnected chiral susceptibility in Figure 3.777The disconnected part of the susceptibility is expected to represent the largest contribution to the phase transition’s peak  It is remarkable that the phase transition appears to happen at also with the improved lattice action, see Figure 3.
This temperature approximately coincides with the deconfinement phase transition, which was found in  to be second order and to have the critical behavior of the Z Ising model. The deconfinement transition occurs thus at the point where the Binder cumulant of the Polyakov loop
reaches the critical value 0.46548(5) . This point is shown in the plots on the right side of Figure 2. As the vanishing of the condensate in the chiral and continuum limit also signals a true phase transition at this point, our results show that, up to numerical uncertainties, in SU(2) supersymmetric Yang-Mills theory, the deconfinement phase transition and restoration of the discrete Z chiral symmetry occur simultaneously. This fact is far from trivial. Indeed, some non-supersymmetric QCD-like theories, like YM with several adjoint spinors, have been observed to exhibit phases with mixed deconfining and broken chiral symmetry phases .
An interesting final quantity determined by our study is the deconfinement temperature in the chiral/supersymmetric limit. We compare our result to the deconfinement critical temperature of SU(2) YM found in  through the ratio
where the two temperatures are compared in dimensionless units as , and the scale for pure gauge has been computed in . This value is compatible with our previous investigations in  at a smaller . The ratio furthermore roughly agrees with an analytical prediction found in , where it is claimed that .
5.3 Prediction from string theory
The result that deconfinement and restoration of the quantum chiral symmetry are occurring at the same critical point is far from trivial and can only be confirmed through non-perturbative methods. Lattice simulations can compute an estimation of the order parameters, but they do not provide immediately a qualitative physical interpretation of the mechanisms responsible for such results. A deep understanding of the dynamics of non-Abelian gauge theories in general and QCD in particular is, indeed, still missing and concepts like supersymmetry and string theory arose in the attempt to achieve it. An interpretation developed in the context of string theory might indeed provide such a deeper understanding of the mechanism behind our results.
In  Witten considered certain brane configuration consisting of two differently oriented NS five-branes and D-four-branes stretching between them in weakly coupled IIA superstring theory.888The framework is actually M-theory. The brane model is not equivalent to SYM but is in the same universality class. SYM is obtained when taking the IIA, i.e. ten dimensional limit. The effective theory on the world-volume of the four-brane is a 3+1 dimensional SU(N) gauge theory with supersymmetry. This theory has vacua, which can be identified with the pattern of chiral symmetry breaking of the supersymmetric QFT. Witten succeeded to show that the confining string emanating from an external quark is topologically equivalent to the IIA fundamental string. Furthermore he found that the BPS-saturated domain wall can actually be identified with a D-brane, on which the confining strings can end. This fact directly relates confinement and chiral symmetry breaking. Indeed, the domain wall is expected to have restored chiral symmetry in its core, while a color-electric source sufficiently near the wall behaves as a free quark, and the Polyakov loop expectation value does not vanish. This was investigated further in  by means of an effective field theory for the gaugino condensate and the Polyakov loop. The authors found that, for , Z restoration implies Z breaking and that the results of Witten, i.e. that the confining strings end on the domain wall, can only hold if both phase transitions occur simultaneously.
6 Conclusions and Outlook
The gradient flow has enabled us to explore the confining and the chiral properties of the SYM both at zero and finite temperature. We have been able to extrapolate the chiral condensate at zero temperature to the mass-less limit and to prove that chiral symmetry is broken. Our findings with the Wilson fermion action are in agreement with previous studies with Ginsparg-Wilson fermions, avoiding, however, the huge numerical cost of preserving chiral symmetry on the lattice. A precise quantitative comparison of the results obtained with the different fermion actions would be possible once a common renormalization scheme would be chosen to fix the multiplicative renormalization constant.
We have also explored the phase diagram of the theory at non-zero temperature. The chiral condensate develops an abrupt change of behavior as the gaugino mass becomes smaller, compatible with a second order phase transition. In comparison to our previous investigations without gradient flow in , the jump in the order parameter’s expectation value is clearer and thus more significant, which reflects the advantages of the method. Further, the study of the Binder cumulant of the Polyakov loop has been important to locate precisely the deconfinement phase transition. Through the identification of both phase transitions we found that the critical temperatures are in fact coincident. This leads us to the remarkable conclusion that chiral symmetry restoration and deconfinement are not independent and uncorrelated non-perturbative phenomena in supersymmetric Yang-Mills theory, but they may obey a common underlying dynamics. As pointed out in the end of the present work, this observation seems to be in agreement with previous semiclassical predictions found through certain brane configuration in IIA superstring/M-theory.
We are currently working towards new applications of the gradient flow. One immediate further step is to investigate thermal SYM for SU(3) gauge group. Another very interesting direction is to take profit of the method to study the vacuum’s structure of the theory at zero temperature. In addition we plan to investigate possible applications for the renormalization of the supercurrent in theories with extended supersymmetry .
The configurations and the parameters of this study have been created in the long term effort of the DESY-Münster collaboration to study the properties of SYM. We thank in particular Gernot Münster and Istvan Montvay for helpful comments and discussions. The authors gratefully acknowledge the Gauss Centre for Supercomputing e. V. (www.gauss-centre.eu) for funding this project by providing computing time on the GCS Supercomputers JUQUEEN, JURECA, and JUWELS at Jülich Supercomputing Centre (JSC) and SuperMUC at Leibniz Supercomputing Centre (LRZ). G. Bergner and C. López acknowledge support from the Deutsche Forschungsgemeinschaft (DFG) Grant No. BE 5942/2-1.
-  N. Seiberg and E. Witten, Electric - magnetic duality, monopole condensation, and confinement in N=2 supersymmetric Yang-Mills theory, Nucl. Phys. B426 (1994) 19-52, [arXiv: hep-th/9407087 ].
-  N. Seiberg, Electric-Magnetic Duality in Supersymmetric Non-Abelian Gauge Theories, Nucl. Phys. B435(1995) 129-146, [arXiv: hep-th/9411149 ].
-  K. Intriligator, N. Seiberg, Lectures on supersymmetric gauge theories and electric-magnetic duality, Nucl. Phys. Proc. Suppl. 45BC:1-28 (1996), [arXiv: hep-th/9509066 ].
-  N. M. Davies, T. J. Hollowood, V. V. Khoze, M. P. Mattis: Gluino Condensate and Magnetic Monopoles in Supersymmetric Gluodynamics, Nucl. Phys. B 559 (1999) 123, [arXiv: hep-th/9905015 ].
-  D. Amati, K. Konishi, Y. Meurice, G. C. Rossi and G. Veneziano: Nonperturbative Aspects in Supersymmetric Gauge Theories, Phys. Rept. 162 (1988) 169.
-  M. A. Shifman and A. I. Vainshtein, On Gluino Condensation in Supersymmetric Gauge Theories. SU(N) and O(N) Groups, Nucl. Phys. B 296 (1988) 445.
-  A. Y. Morozov, M. A. Olshanetsky and M. A. Shifman, Gluino condensate in supersymmetric gluodynamics, Nucl. Phys. B bf 304 (1988) 291.
-  M. Lüscher, Properties and uses of the Wilson flow in lattice QCD, JHEP 1008 (2010) 071 Erratum: [JHEP 1403 (2014) 092] [arXiv: 1006.4518 [hep-lat]].
-  M. Lüscher and P. Weisz, Perturbative analysis of the gradient flow in non-abelian gauge theories, JHEP 1102 (2011) 051 [arXiv: 1101.0963 [hep-th]].
-  M. Lüscher, Chiral symmetry and the Yang–Mills gradient flow, JHEP 1304 (2013) 123 [arXiv: 1302.5246 [hep-lat]].
-  Y. Taniguchi et al. [WHOT-QCD Collaboration], Energy-momentum tensor correlation function in full QCD at finite temperature, [arXiv: 1711.02262 [hep-lat]].
-  Y. Taniguchi, S. Ejiri, R. Iwami, K. Kanaya, M. Kitazawa, H. Suzuki, T. Umeda and N. Wakabayashi, Exploring = 2+1 QCD thermodynamics from the gradient flow, Phys. Rev. D 96 (2017) no.1, 014509
-  A. Carosso, A. Hasenfratz and E. T. Neil, Nonperturbative Renormalization of Operators in Near-Conformal Systems Using Gradient Flows, Phys. Rev. Lett. 121 (2018) no.20, 201601 [arXiv: 1806.01385 [hep-lat]].
-  K. Hieda, A. Kasai, H. Makino and H. Suzuki, 4D SYM supercurrent in terms of the gradient flow, PTEP 2017 (2017) no.6, 063B03 [arXiv: 1703.04802 [hep-lat]].
-  A. Kasai, O. Morikawa and H. Suzuki, Gradient flow representation of the four-dimensional super Yang–Mills supercurrent, PTEP 2018 (2018) no.11, 113B02 [arXiv: 1808.07300 [ hep-lat]].
-  K. Kikuchi and T. Onogi, Generalized Gradient Flow Equation and Its Application to Super Yang-Mills Theory, JHEP 1411 (2014) 094 [arXiv: 1408.2185 [hep-th]].
-  D. Kadoh and N. Ukita, Supersymmetric gradient flow in SYM, [arXiv: 1812.02351 [hep-th]].
-  S. R. Coleman and J. Mandula, All Possible Symmetries of the S Matrix, Phys. Rev. 159 (1967) 1251.
-  G. Veneziano and S. Yankielowicz, An Effective Lagrangian for the Pure N=1 Supersymmetric Yang-Mills Theory, Phys. Lett. 113B (1982) 231.
-  G. Münster, H. Stüwe, The mass of the adjoint pion in N=1 supersymmetric Yang-Mills theory, JHEP05(2014)034, [arXiv: 1402.6616 [hep-th]].
-  G. Bergner, P. Giudice, G. Münster, I. Montvay and S. Piemonte, The light bound states of supersymmetric SU(2) Yang-Mills theory, JHEP 1603 (2016) 080
-  E. Witten, Constraints on Supersymmetry Breaking, Nucl. Phys. B 202 (1982) 253.
-  T. R. Taylor, G. Veneziano and S. Yankielowicz, Supersymmetric QCD and Its Massless Limit: An Effective Lagrangian Analysis, Nucl. Phys. B 218 (1983) 493.
-  I. Affleck, M. Dine and N. Seiberg, Dynamical Supersymmetry Breaking in Supersymmetric QCD, Nucl. Phys. B 241 (1984) 493.
-  M. A. Shifman and A. I. Vainshtein, On Gluino Condensation in Supersymmetric Gauge Theories. SU(N) and O(N) Groups, Nucl. Phys. B 296 (1988) 445 [Sov. Phys. JETP 66 (1987) 1100].
-  A. Y. Morozov, M. A. Olshanetsky and M. A. Shifman, Gluino Condensate in Supersymmetric Gluodynamics, Sov. Phys. JETP 67 (1988) 222 [Zh. Eksp. Teor. Fiz. 94 (1988) 18].
-  D. Amati, K. Konishi, Y. Meurice, G. C. Rossi and G. Veneziano, Nonperturbative Aspects in Supersymmetric Gauge Theories, Phys. Rept. 162 (1988) 169.
-  J. Terning, Modern supersymmetry: Dynamics and duality, International series of monographs on physics. 132
-  M. Schwetz and M. Zabzine, Gaugino condensate and Veneziano-Yankielowicz effective Lagrangian, [arXiv: hep-th/9710125 ]
-  G. R. Dvali and M. A. Shifman, Domain walls in strongly coupled theories, Phys. Lett. B 396 (1997) 64 Erratum: [Phys. Lett. B 407 (1997) 452]
-  O. Aharony, J. Sonnenschein and S. Yankielowicz, A Holographic model of deconfinement and chiral symmetry restoration Annals Phys. 322 (2007) 1420 [arXiv: hep-th/0604161 ].
-  G. Bergner, P. Giudice, G. Münster, S. Piemonte and D. Sandbrink, Phase structure of the supersymmetric Yang-Mills theory at finite temperature, JHEP 1411 (2014) 049
-  M. Lüscher, Trivializing maps, the Wilson flow and the HMC algorithm, Commun. Math. Phys. 293 (2010) 899 [arXiv: 0907.5491 [hep-lat]].
-  R. Kirchner et al. [DESY-Munster Collaboration], Evidence for discrete chiral symmetry breaking in N=1 supersymmetric Yang-Mills theory, Phys. Lett. B 446 (1999) 209 [arXiv: hep-lat/9810062 ].
-  S. Ali, G. Bergner, H. Gerber, P. Giudice, I. Montvay, G. Münster, S. Piemonte and P. Scior, The light bound states of supersymmetric SU(3) Yang-Mills theory on the lattice, JHEP 1803 (2018) 113 [arXiv: 1801.08062 [hep-lat]].
-  G. T. Fleming, J. B. Kogut and P. M. Vranas, SuperYang-Mills on the lattice with domain wall fermions, Phys. Rev. D 64 (2001) 034510 doi:10.1103/PhysRevD.64.034510 [arXiv: hep-lat/0008009 ].
-  J. Giedt, R. Brower, S. Catterall, G. T. Fleming and P. Vranas, Lattice super-Yang-Mills using domain wall fermions in the chiral limit, Phys. Rev. D 79 (2009) 025015 [arXiv: 0810.5746 [hep-lat]].
-  M. G. Endres, Dynamical simulation of N=1 supersymmetric Yang-Mills theory with domain wall fermions, Phys. Rev. D 79 (2009) 094503 [arXiv: 0902.4267 [hep-lat]].
-  J. Giedt, Progress in four-dimensional lattice supersymmetry, Int. J. Mod. Phys. A 24 (2009) 4045 [arXiv: 0903.2443 [hep-lat]].
-  S. W. Kim et al. [JLQCD Collaboration], Lattice study of 4d N=1 super Yang-Mills theory with dynamical overlap gluino, PoS LATTICE 2011 (2011) 069 [arXiv: 1111.2180 [hep-lat]].
-  A. M. Ferrenberg, J. Xu and D. P. Landau, Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model, Phys. Rev. E 97 (2018) no.4, 043301 [arXiv: 1806.03558 [physics.comp-ph]].
-  F. Karsch and M. Lutgemeier, Deconfinement and chiral symmetry restoration in an SU(3) gauge theory with adjoint fermions, Nucl. Phys. B 550 (1999) 449
-  G. Cella, G. Curci, R. Tripiccione and A. Vicere, Scaling, asymptotic scaling and Symanzik improvement. Deconfinement temperature in SU(2) pure gauge theory, Phys. Rev. D 49 (1994) 511
-  P. Giudice and S. Piemonte, Improved thermodynamics of SU(2) gauge theory, Eur. Phys. J. C (2017) 77: 821, [arXiv: 1708.01216 [hep-lat]]
-  G. Lacroix, C. Semay and F. Buisseret, The deconfined phase of SUSY Yang-Mills: bound states and the equation of state, [arXiv: 1408.4979 [hep-th]].
-  E. Witten, Branes and the dynamics of QCD, Nucl. Phys. B 507 (1997) 658 [arXiv: hep-th/9706109 ].
-  A. Campos, K. Holland and U. J. Wiese, Complete wetting in supersymmetric QCD or why QCD strings can end on domain walls, Phys. Rev. Lett. 81 (1998) 2420 [arXiv: hep-th/9805086 ].