PRDPRVDAQPhysical Review DPhys. Rev. D2470-00102470-0029American Physical Society10.1103/PhysRevD.100.074501ARTICLESLattice field theories, lattice QCDStudy of center and chiral symmetry realization in thermal N=1 super Yang-Mills theory using the gradient flowSTUDY OF CENTER AND CHIRAL SYMMETRY REALIZATION …BERGNER GEORG, LÓPEZ CAMILO, AND PIEMONTE STEFANOhttps://orcid.org/0000-0002-7325-2220BergnerGeorg*https://orcid.org/0000-0001-5654-7804LópezCamilo†University of Jena, Institute for Theoretical Physics, Max-Wien-Platz 1, D-07743 Jena, GermanyPiemonteStefano‡Institute for Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany
georg.bergner@uni-jena.de
camilo.lopez@uni-jena.de
stefano.piemonte@ur.de
4October20191October2019100707450116July2019Published by the American Physical Society2019authorsPublished by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP3.
The realization of center and chiral symmetries in N=1 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 nonvanishing gaugino condensate, which breaks the nonanomalous Z2Nc chiral symmetry. In this work we find, for the first time, a nonvanishing condensate at zero temperature and zero gaugino mass with Wilson fermions. This is achieved by means of the renormalization properties of the Yang-Mills and fermion gradient flows, which are independent of the lattice regularization. At finite temperatures, we find that the phase transitions corresponding to deconfinement and chiral restoration occur at roughly the same critical temperature for SU(2) gauge group, implying the saturation of the bound Tchiralc≤Tdecc, recently obtained through anomaly matching. We furthermore discuss the agreement of our findings with conjectures from superstring theory.
Gauss Centre for Supercomputing e. V.Deutsche Forschungsgemeinschaft10.13039/501100001659BE 5942/2-1INTRODUCTION
Supersymmetric Yang-Mills theories (SYM) have been a useful laboratory to extend our understanding of nonperturbative phenomena. For instance, the study of supersymmetry combined with gauge symmetry has led to the discovery of duality as a feature of supersymmetric gauge theories. Besides the well-known gauge-gravity duality between the conformal N=4 SYM and string theory on the curved space AdS5×S5, electromagnetic duality has been conjectured for N=2 SYM and N=1 supersymmetric QCD [1–3]. The most interesting open challenge is the search for theories that are more similar to QCD and still allow for an analytical understanding of confinement and strong interactions due to supersymmetry. N=1 SYM is the supersymmetric extension of the pure gauge sector of QCD describing the strong interactions of gluons and gauginos. It shares many interesting aspects with QCD while retaining minimal supersymmetry. It is expected that this theory should exhibit confinement, chiral symmetry breaking and formation of a mass-gap [4–7]. N=1 SYM is therefore a perfect candidate for the analytical understanding of nonpertubative phenomena of non-Abelian gauge theories, like confinement and chiral symmetry breaking.
Numerical lattice simulations are the natural tool to explore nonperturbatively the phase diagram of strongly interacting gauge theories, such as N=1 SYM, and complement the analytical investigations. The aim of this work is to apply numerical methods in a search for interesting patterns in the phase transitions of N=1 SYM and to relate them to analytical arguments. In a previous study we have found an interesting indication for a coincidence of the chiral and deconfinement transition of this theory without an analytical explanation. In this work we find that this signal persists at finer lattice spacings and using refined numerical methods. Moreover we show that it is confirmed by analytical arguments from a string theory perspective. Another interesting open aspect of our previous study is the relation of the numerical data to the analytic prediction for the gluino condensate at zero temperature. The gradient flow methods that we apply in the current work allow the first step toward this relation: we are able to extrapolate the finite value of the gluino condensate in the chiral limit without additive renormalization. The remaining issue that is left open in our investigations is the multiplicative renormalization required to match the scheme of the analytical predictions.
The investigations of chiral symmetry breaking of N=1 SYM on the lattice requires additional considerations. We use the Wilson fermion discretization that breaks chiral symmetry explicitly. In addition, the lattice discretization also unavoidably breaks supersymmetry. We have seen in our previous studies that it is possible to recover both symmetries in the continuum limit if the gaugino mass is tuned to the chiral limit. However, the breaking still leads to a divergent additive renormalization of the chiral condensate.
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–10]. In particular, the flowed chiral condensate is free from the additive renormalization of Wilson fermions [10]. 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 [13], 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 work we present an extended study of the phase diagram of N=1 SYM with the gauge group SU(2) at zero and nonzero temperature using the gradient flow. This method is summarized in Sec. IV. The details of our lattice discretization as well as the expected phase transitions are explained in Secs. II and III. We measure the chiral condensate expectation value at positive flow time. We show strong evidences that chiral symmetry is spontaneously broken at zero temperature by a nonvanishing expectation value of the gaugino condensate in Sec. V. At high temperature chiral symmetry is restored and that a phase transition occurs in the massless limit. In Sec. VI, we provide evidence 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 indications of mixed phases where deconfinement occurs while chiral symmetry is broken in our numerical data. This section contains also an interpretation of our results from the point of view of superstring/M-theory.
N=1 SUPER YANG-MILLS THEORY ON THE LATTICE
N=1 super Yang-Mills theory is the minimal four-dimensional gauge theory consistent with supersymmetry. It is therefore the supersymmetric gauge theory most similar to Yang-Mills theory and QCD. The four-dimensional Euclidean on-shell action is S=∫d4x(14FaμνFμνa+12λ¯aDλa+mλ¯aλa+θ32π2εμνρσFaμνFaρσ),where F is the usual Yang-Mills (YM) field strength, λ the gaugino, a Majorana spinor in the adjoint representation, and D the covariant derivative. The mass term m breaks supersymmetry softly. The last term is a topological term, whose space-time integration yields the winding number Qtop of the gauge field, which we neglect in our investigations.
On the lattice, spinors are site variables and Yang-Mills fields Aμ are represented as links through the parallel transporter Uμ=eigaAμbτbF, being a the lattice spacing and τF the Lie generators in the fundamental representation. We employ the following SYM lattice action in our simulations [18]Slat=∑xRetr{βNc∑μ≠ν(53Pμν(x)-112Rμν(x))}+12∑x,yλ¯(y)Dw[Vμ](y,x)λ(y),where the Dirac-Wilson operator Dw(x,y)λ(y)=λ(x)-κ∑μ=14[(1-γμ)Vμ(x)λ(x+μ)+(1+γμ)Vμ†(x-μ)λ(x-μ)],has parallel transporters V in the adjoint representation, defined as Vμ(x)ab=2tr(Uμ(x)†τaFUμ(x)τbF). Here β is equal to 2Ncg2,Pμν and Rμν are the plaquette and the 2×1 rectangular Wilson loop, respectively. One-level stout smearing with parameter ρ=0.15 has been employed for the links in the fermion action. Finally, κ=12m+8 is the hopping parameter. Our simulations also include runs with a clover improved fermion action in order to test the lattice artefacts comparing two different discretizations.
A nonzero gaugino mass breaks supersymmetry softly and the bare parameter m 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 ma-π2∼m, as was shown in [19] within the framework of partially quenched perturbation theory. The lattice discretization breaks supersymmetry, but it is recovered in the massless and continuum limit, as confirmed by supersymmetric Ward identities [20].
THE PHASE DIAGRAM OF N=1 SUPER YANG-MILLS THEORY
N=1 SYM shares several features with QCD, as it is a confining theory without scalars and it includes spinor fields interacting with a gauge field. As explained in this section, these properties transform N=1 SYM to an important laboratory to gain a deeper understanding of nonperturbative phenomena of QCD such as chiral symmetry breaking and confinement. These phenomena are usually more accessible in supersymmetric theories due to the constraints that supersymmetry imposes.
Confinement of static fundamental charges
At zero temperature, the YM vacuum is expected a confining medium for external static color-electric charges. It can be probed through the Polyakov loop, which is the path ordered product of the links in the fundamental representation along a line which wraps in the compact direction PL=1NtV3∑x→Tr{∏t=1NtU4(x→,t)},where V3 denotes the three-dimensional lattice volume of the noncompact directions. The Polyakov loop can be related to the exponential of the free energy of an isolated static fundamental quark. Hence confinement is detected by a vanishing vacuum expectation value of PL, indicating that isolated fundamental quarks are states with infinite free energy. As the temperature is increased, the vacuum expectation value of the Polyakov loop becomes nonvanishing at some critical temperature, where quark deconfinement occurs.
In the case of QCD, 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-antiquark pair popping up from the vacuum. In contrast, the identification of the deconfinement phase transition is clear for N=1 SYM, as adjoint fermion fields are unable to provide string breaking to fundamental quark-antiquark pairs. The Polyakov loop represents a good order parameter for the deconfinement phase transition for N=1 SYM and a singularity of the partition function is expected at a critical nonzero temperature.
Chiral symmetry
N=1 supersymmetric Yang-Mills theory is invariant under 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
In this section λa is a Weyl spinor, as this notation is more natural for the massless theory.
λa→eiαλa.This symmetry is broken at the quantum level by instanton contributions, as the U(1) chiral rotation is equivalent to θ→θ-2Ncα,where Nc is the number of colors. The path integral is therefore invariant only for α=kπ/Nc, i.e., the chiral symmetry of the quantum theory is actually the discrete subgroup Z2Nc⊂U(1). As the gaugino condensate is not invariant with respect to the full Z2Nc group, i.e., ⟨λaλa⟩→e2iα⟨λaλa⟩, an interesting question is whether the quantum chiral symmetry is further spontaneously broken by a nonvanishing ⟨λaλa⟩. Indeed, analytical calculations show that the answer is affirmative [4–7,18,21–23]. Thanks to remarkable properties of supersymmetric theories like nonrenormalization and holomorphicity of the effective superpotential with respect to fields and couplings, it was proven [24,25] that if the effective theory at long distances is massive with color-singlets as degrees of freedom, the gaugino condensate takes on the form ⟨λaλa⟩=cΛ3.The global Z2Nc symmetry is therefore spontaneously broken down to its discrete subgroup Z2, i.e., to the sign flip λa→-λa. As a consequence, there are Nc degenerate vacua connected by ZNc chiral transformations and domain walls interpolating between different vacua are expected [26].
The exact value of the proportionality constant c in (3) has been computed based on strong coupling and weak coupling instanton analysis. The differences of these approaches have been addressed in compactified SYM, where the semiclassical analysis leads to a controlled approximation [27].
There have been many attempts to generalize the calculations of the condensate towards other gauge theories. The computed value has been compared to the results of one-flavour QCD based on the orientifold planar equivalence. This comparison is based on several assumptions that are discussed in Sec. V.
The phase diagram at nonzero temperatures
In thermal quantum field theory, temperature corresponds to the inverse radius of a compactified direction of the path integral with thermal boundary conditions, i.e., periodic for bosons and antiperiodic for fermions. The phases of N=1 SYM are related to the response of the mentioned properties of the zero temperature vacuum to the finite temperature. At some critical temperature the condensate is expected vanish and the spontaneously broken chiral symmetry is restored. In addition, center symmetry breaks spontaneously at the deconfinement temperature and the Polyakov loop should acquire a nonvanishing expectation value. An interesting question is whether both critical temperatures coincide, as it would be expected if confinement and chiral symmetry breaking share the same nonperturbative 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 [28]. In many 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, and no exact order parameter can be found which unambiguously signals the pseudocritical temperature.
Fortunately, both the Polyakov loop and the gaugino condensate are exact order parameters for the respective phase transitions in N=1 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 chiral condensate. The lattice realization with Wilson fermions on the lattice breaks chiral symmetry, which implies besides the multiplicative renormalization factor (Zλ¯λ) an additional additive renormalization (b0) of the gaugino condensate ⟨λ¯λ⟩R=Zλ¯λ(β)(⟨λ¯λ⟩B-b0(β)).
In previous investigations [29], the additive constant b0 was removed by subtracting the bare condensate at zero temperature, ⟨λ¯λ⟩S=⟨λ¯λ⟩BT=0-⟨λ¯λ⟩BTThe downside of this approach is the fact that the subtracted condensate is forced to vanish at T=0. Although this subtraction should preserve the behavior of the order parameter near the critical temperature, it is not possible to determine whether the renormalized condensate at zero temperature is nonzero. Hence, the picture of the realization of chiral symmetry in N=1 SYM is, following this method, incomplete. As it will be explained in the next section, the additive renormalization is not necessary when the gradient flow is used, allowing the computation of the condensate at zero temperature and, moreover, allowing a comparison to fermion discretizations which satisfy the Ginsparg-Wilson relation.
THE GRADIENT FLOWFlow equations
Motivated in the context of trivializing maps [30], Lüscher studied correlation functions of fields flowed through the equations [9]
The term Dμ∂νBν is a gauge parameter, which is included for mere technical reasons.
∂tBμ=DνGνμ+Dμ∂νBν,Bμ|t=0=Aμ,Gμν=∂μBν-∂νBμ+[Bμ,Bν]∂tχ=Δχ∂tχ¯=χ¯Δ←,χ|t=0=ψχ¯|t=0=ψ¯,Δ=DμDμ,where B (resp. A) and G are the Yang-Mills gauge field and field strength, respectively. χ and ψ are spinor fields and Dμ is the gauge-covariant derivative in the adjoint representation. The parameter t∈R describes a flow on the vector space of gauge fields.
The flow equations have a smoothening effect on the fields, which are Gaussian-like smeared over an effective radius rt=8t, as it can be easily seen by integrating equation (4) in the noninteracting limit [B,B]∼0, Bμ,1(t,x)=∫dDyKt(x-y)Aμ(y),Kt(z)=∫dDp(2π)Deipze-tp2=e-z2/4t(4πt)D/2,where D is the number of space-time dimensions and K the heat-kernel. The term e-tp2 regularizes the integral in momentum space when t>0, removing the UV divergences at large momenta. Some years ago it was shown that through the gradient flow the smearing property remains at all orders in perturbation theory [9]. Correlation functions of monomials of flowed gauge-fields are renormalized without extra counterterms, while spinor operators renormalize multiplicatively according to the field content [10]. An advantage of the gradient flow is its intrinsic regularization-scheme independence, and the renormalization properties of the flow hold also once discretized on the lattice. Therefore, currents and densities which are explicitly broken by the lattice discretization should be easily accessible at positive flow-time. In particular, the complications of explicit chiral symmetry breaking of the Wilson discretization of the fermion action are solved in the sense that the flowed chiral condensate renormalizes only multiplicatively.
Written explicitly, the flowed bare gaugino condensate is [10]⟨λ¯(x,t)λ(x,t)⟩=-∫dDvdDwK(t,x;0,v)S(v,w)K(t,x;0,w)†.The flowed condensate is equivalent to the action of the heat-kernel on the fermion propagator S. The condensate at a fixed flow time that is computed in this way is proportional to any other renormalized condensate. It defines a correct order parameter for the chiral symmetry breaking even for Wilson fermions. The multiplicative factor in comparison to the MS¯ scheme can be determined from a small flow-time expansion. Consequently, it is possible to study if ⟨λ¯λ⟩≠0 at T=0 even with Wilson fermions. The value of Zλ¯λ is irrelevant if the bare lattice gauge coupling is kept fixed.
The numerical integration of the flow is relatively straightforward, and the methods we employ to compute the flowed gluino condensate are described in Appendix A.
GAUGINO CONDENSATE AT ZERO TEMPERATURE
The gradient flow method allows for the computation of the gaugino condensate at zero temperature, without an additive renormalization. We consider four ensembles from our previous investigations of SU(2) SYM at β=1.75 on a 323×64 lattice with hopping-parameter values κ∈{0.1490,0.1492,0.1494,0.1495}.These four points are used to extrapolate to the chiral point, i.e., to vanishing renormalized gaugino mass. The scale in physical units is set by the gradient flow observables t0, which is defined as t2⟨ϵ(t)⟩|t=t0=0.3,ϵ(t)=14Gμνa(t)Gμνa(t),with ϵ the field energy density and Gμνa the flowed Yang-Mills field strength. The t0 values for each κ value are summarized in Table I. For the computation of the flowed condensate, a mass-independent scheme was used, i.e., a common flow time was chosen corresponding to the chiral extrapolated t0chiral=12.81(35). The flow equations were integrated up to t0 and the condensate obtained from Eq. (A1). In the following we use the condensate in units of t0, defined as: ⟨λ¯λ⟩≡t03/2⟨λ¯λ(t)⟩|t=t0.
I10.1103/PhysRevD.100.074501.t1
The gradient flow scale t0 and the values of the gluino condensate for each of the ensembles.
As shown in Fig. 1, we obtain a nonzero value of the gluino condensate in the chiral limit. For the numerical simulations with Wilson fermions, this is the first clear evidence 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 nonzero condensate has been possible from a double peak of the histogram that appears in rather unstable simulations on small volumes close to the chiral point [31,32].
110.1103/PhysRevD.100.074501.f1
Extrapolation of the chiral condensate to vanishing renormalized gaugino mass. The gaugino condensate scales almost linearly with the lattice spacing a, while some mild quadratic dependence is noticeable as the gaugino mass grows. Most importantly, the condensate does not vanish at the chiral point, which implies the spontaneously symmetry breaking of the Z2Nc symmetry. The extrapolated values are ⟨λ¯λ⟩chiral=0.0002304(19) (linear) and ⟨λ¯λ⟩chiral=0.0002235(36) (quadratic).
Our results are of the same significance as the ones obtained with domain-wall and overlap fermions in [33–37]. In these simulations the renormalization factor of the condensate has not been determined and therefore the results are, as in our case, proportional to the renormalized condensate in MS¯ scheme. The gradient flow method opens up the possibility to study the different phases at zero temperature without the computational more expensive Ginsparg-Wilson fermions. A direct comparison of the results with different lattice actions would have been possible from the small flow-time expansion combined with finite-volume methods.
The comparison of our numerical result with strong and weak coupling instanton calculations (3) requires two further important steps, as already outlined in the comparison of one-flavor QCD with numerical results [38]. The first step is the matching of scale Λ, which, according to [38], can be set to the value in the MS¯ scheme. We have already done this analysis and found a value of ΛMS¯ in SYM that is quite compatible to QCD [39]. The second step requires the determination of the renormalization factor Zλ¯λ of the condensate for the exact result. The analysis of [38] shows that the renormalization scheme dependence can be rephrased in terms of a second proportionality factor c˜ multiplied to c in (3). This factor is assumed to have a subleading correction in the large Nc expansion, but there is so far no clear determination at small Nc. Without an analytic prediction for the additional factor, it is not possible to match the semiclassical predictions to any numerical value.
DISCRETE CHIRAL SYMMETRY RESTORATION AND QUARK DECONFINEMENT
The second main purpose of this study is to investigate the realization of chiral and center symmetries in SU(2) SYM at finite temperature. Some first results in this direction were obtained in [29], and now we revise our early investigation using the gradient flow on four different new lattice ensembles at β=1.75 and hopping parameter κ∈{0.1480,0.1490,0.14925}. To further cross-check the validity of our results, we have also analyzed a set of ensembles generated at β=1.65 and κ=0.175 using the tree-level clover improved fermion action with unsmeared links. Nonzero 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 243×Nt for Nt∈{4,5,6,…,48}, with the upper bound parametrizing the zero-temperature limit. For this finite-temperature analysis, a mass-dependent scale setting was used. This means that the flow time of the gaugino condensate is different for each κ value, in contrast to the zero-temperature measurements. One reason for this choice is that the integration of the adjoint flow equation gets very expensive as the flow time grows. A second reason is that the possible range of t is limited at smaller Nt. A reliable value can not be obtained once the smearing radius becomes compatible with Nt. We have checked that around the mass dependent value of t0 the condensate shows only a very weak t dependence. Therefore we expect no relevant change of the results for a mass independent choice of t0.
The results are summarized in Figs. 2–4 and Table II. 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 pseudocritical temperature is less pronounced for a larger explicit symmetry breaking. At the smallest gaugino mass (κ=0.14925) 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 Fig. 4.
The disconnected part of the susceptibility is expected to represent the largest contribution to the phase transition’s peak [29].
It is remarkable that the phase transition appears at T=t0/Nt∼0.25 also with the improved lattice action, see Fig. 4.
210.1103/PhysRevD.100.074501.f2
Temperature dependence of the bare gaugino condensate is shown for κ=0.1480,0.1490 and 0.14925, with the smallest κ in the uppermost plot. On the undermost graphic, the point corresponding to the highest temperature appears to show a growth in the condensate. This is however a nonphysical oversmoothing artefact due to the fact that 8t0>Nt in that region.
310.1103/PhysRevD.100.074501.f3
Binder cumulant of the Polyakov loop for the same κ values of Fig. 2.
410.1103/PhysRevD.100.074501.f4
At the top: Chiral susceptibility for κ=0.14925. The gray band denotes the critical temperature of the Polyakov loop. At the bottom: gaugino condensate at β=1.65 and κ=0.175. The phase transition occurs at Ntc=7, which roughly agrees in dimensionless units with the critical temperature for β=1.75.
This temperature approximately coincides with the deconfinement phase transition, which was found in [29] to be second order and to have the critical behavior of the Z2 Ising model. The deconfinement transition occurs thus at the point where the Binder cumulant of the Polyakov loop B4(PL)=1-13⟨PL4⟩⟨PL2⟩2reaches the critical value B4c=0.46548(5)[40]. This point is shown in the plots of Fig. 3. In the fixed scale approach, the temperature can be changed only by discrete steps. In order to precisely estimate the critical point of the deconfinement phase transition, we perform a linear interpolation of the Binder cumulant at two temperatures where B4 approaches B4c both from the left and from the right. After finding the critical temperature for each κ value, the results were extrapolated to the chiral limit by means of a linear fit (see Fig. 5 and Table III).
510.1103/PhysRevD.100.074501.f5
Chiral extrapolation of the deconfinement phase transition. The deconfinement temperature extrapolated to the chiral limit is Tc|ma–π=0=0.2574(26).
Our results show that, up to numerical uncertainties, in N=1 SU(2) supersymmetric Yang-Mills theory, the deconfinement phase transition and restoration of the discrete Z2Nc chiral symmetry occur simultaneously. This fact is far from trivial. Indeed, some nonsupersymmetric QCD-like theories, like YM with several adjoint spinors, have been observed to exhibit phases with mixed deconfining and broken chiral symmetry phases [41]. In our previous investigations, we have not provided analytical confirmation for our finding. Here we provide arguments that predict the coincidence of the transitions in Sec. VI A.
The agreement of the two transitions is obtained up to the systematic uncertainties of our numerical determination which includes possible supersymmetry breaking lattice artefacts. Based on our previous investigations of particle spectrum and Ward identities [20], the role of these lattice artefacts is expected to be small. Furthermore our current results obtained with an improved lattice action and at a smaller lattice spacing are consistent with our earlier findings [29]. Consequently, we expect any systematic effect to be much smaller than the large difference of the two transition observed in [41].
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 [42] through the ratio Tc(SYM)Tc(YM)=0.2574(26)0.3082(2)=0.8352(90),where the two temperatures are compared in dimensionless units as t0/Ntc, and the scale t0 for pure gauge has been computed in [43]. This value is compatible with our previous investigations in [29] at a smaller β. The ratio furthermore roughly agrees with an analytical prediction found in [44], where it is claimed that Tc(SYM)Tc(YM)=23∼0.82. It is important to stress that the ratio of Eq. (8) depends on the observable used to set the scale of the two different theories, in our case t0, however a general agreement is expected among scales defined from purely gluonic operators.
Prediction for the relation of chiral and deconfinement transition from string theory and t’Hooft anomaly matching
The result that deconfinement and restoration of chiral symmetry occur at the same critical point is far from trivial and can only be confirmed by nonperturbative 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. We provide here some rough arguments based on t’Hooft anomaly matching and a more detailed picture derived from the string theory perspective.
The coincidence of the transitions is in agreement with the predictions of [45,46]. As shown in these references based on anomaly matching arguments, the unbroken center symmetry must imply broken Z2Nc chiral symmetry in N=1 supersymmetric Yang-Mills, meaning that Tchiralc≥Tdeconfinementc. The same constraint was also found from anomaly matching for SU(Nc) with adjoint and fundamental matter in [47]. Our results indicate that this inequality is saturated, as predicted remarkably for adjoint QCD and SYM on R3×S1 in [48–51].
In [52] Witten considered certain brane configuration consisting of two differently oriented NS five-branes and Nc D-four-branes stretching between them in weakly coupled IIA superstring theory.
The 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(Nc) gauge theory with N=1 supersymmetry. This theory has Nc 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 [53] by means of an effective field theory for the gaugino condensate and the Polyakov loop. The authors found that, for Nc=3,Z3χ restoration implies Z3c 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.
CONCLUSIONS AND OUTLOOK
The gradient flow has enabled us to explore the confining and the chiral properties of the N=1 SYM both at zero and finite temperature. We have been able to extrapolate the chiral condensate at zero temperature to the massless 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 is possible once a common renormalization scheme is chosen to fix the multiplicative renormalization constant. The matching to the semiclassical predictions of the condensate would be interesting from the theoretical point of view, however a deeper understanding of the exact calculations of the gluino condensate is required. Any prediction or numerical value of ⟨λ¯λ⟩ holds only in a given scheme, i.e., up to a multiplicative factor, and provides the same information as the condensate at a finite flow time. A comparison to the semiclassical calculations of the gluino condensate would need to find first the missing perturbative connection of the exact formulas of [6] to the MS¯ scheme.
We have also explored the phase diagram of the theory at nonzero temperature. The chiral condensate provides a signal compatible with a second order phase transition at smaller gaugino masses. In comparison to our previous investigations without gradient flow in [29], the signal of the order parameter’s expectation value is clearer and thus more significant, which reflects the advantages of the method. 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 nonperturbative phenomena in N=1 supersymmetric Yang-Mills theory, but they may obey a common underlying dynamics. As pointed out in the end of the present work, this observation is in agreement with 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 [15].
ACKNOWLEDGMENTS
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.
NUMERICAL INTEGRATION OF THE GRADIENT FLOW AND COMPUTATION OF THE GLUINO CONDENSATE
In this Appendix we describe the path for the numerical computation of the gradient flow [10,12]. First, one initial gauge configuration Vt=0 is picked and flowed by integrating the discrete version of the gauge part of Eq. (4)V˙t(x,μ)=-{∂x,μSw(Vt)}Vt(x,μ),U(x,μ)=Vt(x,μ)|t=0,up to some t by means of the following Runge-Kutta integrator with step-size εW0=Vt,W1=exp(1/4)Z0W0,W2=exp(8/9)Z1-(17/36)Z0W1,Vt+ε=exp(3/4)Z2-(8/9)Z1+(17/36)Z0W2,where Zi=-ε∂x,μSw(Wi) and Sw 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 ξs=t=η is generated on the lattice and then integrated by means of the Runge-Kutta integrator down to s=0λ3=ξs+εε,λ2=34Δ2λ3,λ1=λ3+89Δ1λ2,λ0=λ1+λ2+14Δ0(λ1-89λ2),withξsε=λ0.
Further, the vector W(v)=∑wD(v,w)-1ξ(t;0,w) is computed, e.g., through conjugate gradient. Subsequently it is contracted with ξ(t;0,v)† and the result averaged over the lattice sites, i.e., one calculates -1Nlat∑vξk(t;0,v)†W(v).
The discrete version of the gluino condensate ⟨λ¯(x,t)λ(x,t)⟩ in Eq. (7) can be straightforwardly written: -∑v,w⟨tr{K(t,x;0,v)(DW(v,w))-1K(t,x;0,w)†}⟩,where DW is the Wilson-Dirac operator and the trace runs over space-time, spinor and color indices. Following [10], the trace is estimated stochastically by inserting a complete set of random complex vectors η(x) with ⟨η(x)⟩=0 and ⟨η(x)η(y)†⟩=δx,y: ⟨λ¯λ(t)⟩=1NΓ∑x∈Γ⟨λ¯λ(x,t)⟩=-1NΓ∑v,w⟨ξ(t;0,v)†DW-1ξ(t;0,w)⟩,ξ(t;s,w)=∑xK(t,x;s,w)†η(x).
Here ⟨…⟩ denotes the average with respect to both the Monte-Carlo time and any internal group-symmetry representations. Finally, to compute the new vectors ξ, the so-called adjoint flow equation ∂sξ(t;s,w)=-Δ(Vs)ξ(t;s,w),ξ(t;t,w)=η(w),must be integrated from s=t to s=0, i.e., backwards in comparison with the flow equations presented above. Here, the gauge connection Vs in the covariant four-dimensional Laplacian Δ is flowed up to the same t as η.
FINITE TEMPERATURE DATAII10.1103/PhysRevD.100.074501.t2
Condensate, chiral susceptibility and Binder cumulant of the Polyakov loop at β=1.75,243×Nt.
κama-π0.14800.4119(39)0.14900.23780(97)0.149250.1896(17)1N. Seiberg and E. Witten, Electric—magnetic duality, monopole condensation, and confinement in N=2 supersymmetric Yang-Mills theory, B426, 19 (1994).NUPBBO0550-321310.1016/0550-3213(94)90124-42N. Seiberg, Electric-magnetic duality in supersymmetric non-Abelian gauge theories, B435, 129 (1995).NUPBBO0550-321310.1016/0550-3213(94)00023-83aK. Intriligator and N. Seiberg, Lectures on supersymmetric gauge theories and electric-magnetic duality, 45, 1 (1996); NPBSE70920-563210.1016/0920-5632(95)00626-53bK. Intriligator and N. Seiberg34, 237 (1997).SUSEE44N. M. Davies, T. J. Hollowood, V. V. Khoze, and M. P. Mattis, Gluino condensate and magnetic monopoles in supersymmetric gluodynamics, B559, 123 (1999).NUPBBO0550-321310.1016/S0550-3213(99)00434-45D. Amati, K. Konishi, Y. Meurice, G. C. Rossi, and G. Veneziano, Nonperturbative aspects in supersymmetric gauge theories, 162, 169 (1988).PRPLCM0370-157310.1016/0370-1573(88)90182-26M. A. Shifman and A. I. Vainshtein, On Gluino condensation in supersymmetric gauge theories. SU(N) and O(N) Groups, B296, 445 (1988).NUPBBO0550-321310.1016/0550-3213(88)90680-37A. Y. Morozov, M. A. Olshanetsky, and M. A. Shifman, Gluino condensate in supersymmetric gluodynamics, B304, 291 (1988).NUPBBO0550-321310.1016/0550-3213(88)90628-18aM. Lüscher, Properties and uses of the Wilson flow in lattice QCD, 08 (2010) 071; JHEPFG1029-847910.1007/JHEP08(2010)0718bM. LüscherErratum, 03 (2014) 92.JHEPFG1029-847910.1007/JHEP03(2014)0929M. Lüscher and P. Weisz, Perturbative analysis of the gradient flow in non-Abelian gauge theories, 02 (2011) 051.JHEPFG1029-847910.1007/JHEP02(2011)05110M. Lüscher, Chiral symmetry and the Yang–Mills gradient flow, 04 (2013) 123.JHEPFG1029-847910.1007/JHEP04(2013)12311Y. Taniguchi (WHOT-QCD Collaboration), Energy-momentum tensor correlation function in Nf=2+1 full QCD at finite temperature, arXiv:1711.02262.12Y. Taniguchi, S. Ejiri, R. Iwami, K. Kanaya, M. Kitazawa, H. Suzuki, T. Umeda, and N. Wakabayashi, Exploring Nf=2+1 QCD thermodynamics from the gradient flow, 96, 014509 (2017).PRVDAQ2470-001010.1103/PhysRevD.96.01450913A. Carosso, A. Hasenfratz, and E. T. Neil, Nonperturbative Renormalization of Operators in Near-Conformal Systems Using Gradient Flows, 121, 201601 (2018).PRLTAO0031-900710.1103/PhysRevLett.121.20160114K. Hieda, A. Kasai, H. Makino, and H. Suzuki, 4D N=1 SYM supercurrent in terms of the gradient flow, 2017, 063B03 (2017).PTEPCR2050-391110.1093/ptep/ptx07315A. Kasai, O. Morikawa, and H. Suzuki, Gradient flow representation of the four-dimensional N=2 super Yang–Mills supercurrent, 2018, 113B02 (2018).PTEPCR2050-391110.1093/ptep/pty11716K. Kikuchi and T. Onogi, Generalized gradient flow equation and its application to super Yang-Mills theory, 11 (2014) 094.JHEPFG1029-847910.1007/JHEP11(2014)09417D. Kadoh and N. Ukita, Supersymmetric gradient flow in N=1 SYM, arXiv:1812.02351.18G. Veneziano and S. Yankielowicz, An effective Lagrangian for the pure N=1 supersymmetric Yang-Mills theory, 113B, 231 (1982).PYLBAJ0370-269310.1016/0370-2693(82)90828-019G. Münster and H. Stüwe, The mass of the adjoint pion in N=1 supersymmetric Yang-Mills theory, 05 (2014) 034.JHEPFG1029-847910.1007/JHEP05(2014)03420G. Bergner, P. Giudice, G. Münster, I. Montvay, and S. Piemonte, The light bound states of supersymmetric SU(2) Yang-Mills theory, 03 (2016) 080.JHEPFG1029-847910.1007/JHEP03(2016)08021E. Witten, Constraints on supersymmetry breaking, B202, 253 (1982).NUPBBO0550-321310.1016/0550-3213(82)90071-222T. R. Taylor, G. Veneziano, and S. Yankielowicz, Supersymmetric QCD and its massless limit: An effective Lagrangian analysis, B218, 493 (1983).NUPBBO0550-321310.1016/0550-3213(83)90377-223I. Affleck, M. Dine, and N. Seiberg, Dynamical supersymmetry breaking in supersymmetric QCD, B241, 493 (1984).NUPBBO0550-321310.1016/0550-3213(84)90058-024J. Terning, (Oxford University Press, Oxford, 2006), Vol. 132.25M. Schwetz and M. Zabzine, Gaugino condensate and Veneziano-Yankielowicz effective Lagrangian, arXiv:hep-th/9710125.26aG. R. Dvali and M. A. Shifman, Domain walls in strongly coupled theories, 396, 64 (1997); PYLBAJ0370-269310.1016/S0370-2693(97)00131-726bG. R. Dvali and M. A. ShifmanErratum, 407, 452(E) (1997).PYLBAJ0370-269310.1016/S0370-2693(97)00808-327N. M. Davies, T. J. Hollowood, V. V. Khoze, and M. P. Mattis, Gluino condensate and magnetic monopoles in supersymmetric gluodynamics, B559, 123 (1999).NUPBBO0550-321310.1016/S0550-3213(99)00434-428O. Aharony, J. Sonnenschein, and S. Yankielowicz, A Holographic model of deconfinement and chiral symmetry restoration322, 1420 (2007).APNYA60003-491610.1016/j.aop.2006.11.00229G. Bergner, P. Giudice, G. Münster, S. Piemonte, and D. Sandbrink, Phase structure of the N=1 supersymmetric Yang-Mills theory at finite temperature, 11 (2014) 049.JHEPFG1029-847910.1007/JHEP11(2014)04930M. Lüscher, Trivializing maps, the Wilson flow and the HMC algorithm, 293, 899 (2010).CMPHAY0010-361610.1007/s00220-009-0953-731R. Kirchner, I. Montvay, J. Westphalen, S. Luckmann, and K. Spanderen (DESY-Munster Collaboration), Evidence for discrete chiral symmetry breaking in N=1 supersymmetric Yang-Mills theory, 446, 209 (1999).PYLBAJ0370-269310.1016/S0370-2693(98)01523-832S. Ali, G. Bergner, H. Gerber, P. Giudice, I. Montvay, G. Münster, S. Piemonte, and P. Scior, The light bound states of N=1 supersymmetric SU(3) Yang-Mills theory on the lattice, 03 (2018) 113.JHEPFG1029-847910.1007/JHEP03(2018)11333G. T. Fleming, J. B. Kogut, and P. M. Vranas, SuperYang-Mills on the lattice with domain wall fermions, 64, 034510 (2001).PRVDAQ0556-282110.1103/PhysRevD.64.03451034J. Giedt, R. Brower, S. Catterall, G. T. Fleming, and P. Vranas, Lattice super-Yang-Mills using domain wall fermions in the chiral limit, 79, 025015 (2009).PRVDAQ1550-799810.1103/PhysRevD.79.02501535M. G. Endres, Dynamical simulation of N=1 supersymmetric Yang-Mills theory with domain wall fermions, 79, 094503 (2009).PRVDAQ1550-799810.1103/PhysRevD.79.09450336J. Giedt, Progress in four-dimensional lattice supersymmetry, 24, 4045 (2009).IMPAEF0217-751X10.1142/S0217751X0904549237S. W. Kim (JLQCD Collaboration), Lattice study of 4d N=1 super Yang-Mills theory with dynamical overlap gluino, , LATTICE2011 (2011) 069 [arXiv:1111.2180].38A. Armoni, M. Shifman, and G. Veneziano, QCD quark condensate from SUSY and the orientifold large N expansion, 579, 384 (2004).PYLBAJ0370-269310.1016/j.physletb.2003.10.09439G. Bergner and S. Piemonte, Running coupling from gluon and ghost propagators in the Landau gauge: Yang-Mills theories with adjoint fermions, 97, 074510 (2018).PRVDAQ2470-001010.1103/PhysRevD.97.07451040A. M. Ferrenberg, J. Xu, and D. P. Landau, Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model, 97, 043301 (2018).PRESCM2470-004510.1103/PhysRevE.97.04330141F. Karsch and M. Lutgemeier, Deconfinement and chiral symmetry restoration in an SU(3) gauge theory with adjoint fermions, B550, 449 (1999).NUPBBO0550-321310.1016/S0550-3213(99)00129-742G. Cella, G. Curci, R. Tripiccione, and A. Vicere, Scaling, asymptotic scaling and Symanzik improvement. Deconfinement temperature in SU(2) pure gauge theory, 49, 511 (1994).PRVDAQ0556-282110.1103/PhysRevD.49.51143P. Giudice and S. Piemonte, Improved thermodynamics of SU(2) gauge theory, 77, 821 (2017).EPCFFB1434-604410.1140/epjc/s10052-017-5392-644G. Lacroix, C. Semay, and F. Buisseret, The deconfined phase of N=1 SUSY Yang-Mills: bound states and the equation of state, arXiv:1408.4979.45D. Gaiotto, A. Kapustin, Z. Komargodski, and N. Seiberg, Theta, time reversal, and temperature, 05 (2017) 091.JHEPFG1029-847910.1007/JHEP05(2017)09146Z. Komargodski, T. Sulejmanpasic, and M. Ünsal, Walls, anomalies, and deconfinement in quantum antiferromagnets, 97, 054418 (2018).PRBMDO2469-995010.1103/PhysRevB.97.05441847H. Shimizu and K. Yonekura, Anomaly constraints on deconfinement and chiral phase transition, 97, 105011 (2018).PRVDAQ2470-001010.1103/PhysRevD.97.10501148M. M. Anber, E. Poppitz, and M. Unsal, 2d affine XY-spin model/4d gauge theory duality and deconfinement, 04 (2012) 040.JHEPFG1029-847910.1007/JHEP04(2012)04049M. M. Anber, S. Collier, and E. Poppitz, The SU(3)/Z3 QCD(adj) deconfinement transition via the gauge theory/‘affine’ XY-model duality, 01 (2013) 126.JHEPFG1029-847910.1007/JHEP01(2013)12650M. M. Anber, S. Collier, E. Poppitz, S. Strimas-Mackey, and B. Teeple, Deconfinement in N=1 super Yang-Mills theory on R3×S1 via dual-Coulomb gas and “affine” XY-model, 11 (2013) 142.JHEPFG1029-847910.1007/JHEP11(2013)14251M. M. Anber and B. J. Kolligs, Entanglement entropy, dualities, and deconfinement in gauge theories, 08 (2018) 175.JHEPFG1029-847910.1007/JHEP08(2018)17552E. Witten, Branes, and the dynamics of QCD, B507, 658 (1997).NUPBBO0550-321310.1016/S0550-3213(97)00648-253A. Campos, K. Holland, and U. J. Wiese, Complete Wetting in Supersymmetric QCD or Why QCD Strings Can End on Domain Walls, 81, 2420 (1998).PRLTAO0031-900710.1103/PhysRevLett.81.2420