On general ultraviolet properties of a class of confining propagators

We study the ultraviolet properties of theories whose fundamental fields display a confining, Gribov-type, propagator. These are propagators that exhibit complex poles and violate positivity, thus precluding a physical propagating particle interpretation. We show that the properties of this type of confining propagators do not change the ultraviolet behavior of the theory, in the sense that no new ultraviolet primitive divergences are generated, thus securing the renormalizability of these confining theories. We illustrate these properties by studying a variety of models, including bosonic and fermionic confined degrees of freedom. The more intricate case of Super Yang-Mills with ${\cal N} =1$ supersymmetries in the Wess-Zumino gauge is taken as example in order to prove these statements to all orders by means of the algebraic renormalization set up.


Introduction
The quantization of non-abelian gauge theories is still an open and rich subject for quantum field theorists. A general framework which takes into account the non-perturbative phenomena of gluon and quark confinement and of chiral symmetry breaking is still lacking. Needless to say, these issues represent a major challenge for our current understanding of non-abelian gauge theories in the non-perturbative infrared regime.
A successful approach to investigate these topics within the context of the Euclidean quantum field theory is provided by the Gribov framework 1 . In his seminal work [4], Gribov pointed out that the Faddeev-Popov gauge-fixing procedure is plagued by the existence of Gribov copies. In a path integral formulation this manifests in the fact that the Faddeev-Popov operator develops zero modes. For example, in the Landau gauge, ∂ µ A a µ = 0, the Faddeev-Popov operator is given by M ab = −(∂ 2 δ ab − gf abc A c µ ∂ µ ). For sufficiently enough strong coupling constant, zero modes of this operator start to appear, rendering ill-defined the Faddeev-Propov procedure. To deal with this issue, Gribov proposed to restrict the domain of integration in the path integral to a region in field space where the eigenvalues of the Faddeev-Popov operator are strictly positive. This region is known as the Gribov region Ω, being defined as It is useful to observe here that the inverse of the Faddeev-Popov operator, (M −1 ) ab , yields precisely the propagator G(k, A) of the Faddeev-Popov ghosts in the presence of an external gauge field A a µ [1,2,3,4], namely This property can be employed to implement the restriction to the Gribov region Ω by requiring the absence of poles in the ghost propagator for any non-zero value of the ghost external momentum k. This requirement is known as the Gribov no-pole condition. More precisely, following [4], one can always represent the exact ghost propagator in the presence of an external gauge field as where σ(k, A) is the ghost form factor. Using the general properties of the diagrammatic expansion of quantum field theory, we can write where "conn" stands for the connected set of diagrams and 1P I denotes the 1-particle irreducible ones. It can be shown that σ(k, A) 1P I is a decreasing function of k [1,2,3,4]. Therefore, the condition that the ghost propagator has no poles for any non-zero value of the ghost external momentum can be expressed as a condition for the maximum value of the ghost form factor, i.e.
Equation (5) expresses the no-pole condition. In [5], it has been shown that an exact closed expression for σ(0, A) can be obtained, being proportional to the horizon function H(A) of Zwanziger's [1] formalism: .
The no-pole condition (5) is thus equivalent to which is called the horizon condition. The relevance of the horizon function H(A) relies on the fact that the restriction of the domain of integration in the functional integral to the Gribov region Ω can be effectively implemented by adding to the Yang-Mills action the quantity H(A). More precisely, it turns out that, in the thermodynamic limit [1,6,7,8], the partition function of the theory with the cut-off at the Gribov region Ω is given by where the massive parameter γ is known as the Gribov parameter [1,2,3,4]. It is not a free parameter, being determined in a self-consistent way by the horizon condition (7), which can be rewritten as a stationary condition for the vacuum energy E, i.e. where Although expression (8) is non-local, it can be fully localized by introducing a suitable set of auxiliary fields. The Faddeev-Popov measure is localized as usual by means of the Faddeev-Popov ghosts (c a ,c a ) and of the Nakanishi-Lautrup field b a . Moreover, the horizon function H(A) can also be put in a local form [1,6,7,8] by using the auxiliary fields (ω ab µ , ω ab µ ,φ ab µ , ϕ ab µ ), where (φ ab µ , ϕ ab µ ) are a pair of bosonic fields, while (ω ab µ , ω ab µ ) are anti-commuting. The resulting local action is called the Gribov-Zwanziger action S GZ [1,6,7,8], i.e. where with S F P being the Faddeev-Popov action in the Landau gauge and S 0 , S γ given, respectively, by The GZ action (12) displays remarkable properties; it is renormalizable to all orders and has no extra free parameters with respect to the original Faddeev-Popov action (13); the parameter γ is completely determined by the gap equation (9) and does not renormalize independently. This means that the UV properties of the theory are not changed by restricting it to the Gribov region Ω. Only two renormalization factors are in fact needed to renormalize the action (12) [1,6,7,8].
As discussed in [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23], the GZ action (12) breaks the standard BRST symmetry of the Faddeev-Popov action in a soft way, i.e. where and s denotes the standard nilpotent BRST operator, defined by We notice that the breaking term ∆ is of dimension two in the quantum fields. As such, it is a soft breaking, which can be kept under control in the renormalization process [24].
Recently, it has been shown [25] that a non-perturbative nilpotent extension of the standard BRST operator can be constructed in such a way that it is an exact symmetry of the GZ action. In particular, the existence of the soft breaking ∆ turns out to be a consequence of the Ward identity stemming from the non-perturbative BRST exact symmetry [25]. Nevertheless, as far as the UV renormalization of the GZ action is concerned, the standard softly broken BRST symmetry turns out to be very helpful. In fact, using the tools of the algebraic renormalization [24], the softly broken identity (16) can be converted into useful Slavnov-Taylor identities which imply the all order UV renormalizability of expression (12), see for example [10,20,21,22,23].
We observe that, in the local formulation, the horizon condition (9) takes the form which expresses the condensation of the local dimension two operator gf abc A a µ (ϕ bc µ +φ bc µ ). As shown in [9,10,11], this is not the only dimension two condensate present in the theory. The condensation of other dimension two operators, A a µ A a µ and φ ab µ ϕ ab µ −ω ab µ ω ab µ , turns out to be energetically favoured [9,10,11]. The effective action which takes into account the formation of these condensates is known as the Refined-Gribov-Zwanziger (RGZ) action [9,10,11], being given by where the massive parameters (m 2 , µ 2 ) are not independent and have a dynamical origin, being related to the existence of the dimension two condensates A a µ A a µ and φ ab µ ϕ ab µ −ω ab µ ω ab µ . As the GZ action, also the RGZ action can be proven to be tenormalizable to all orders [9,10,11], while displaying the existence of a non-perturbative exact BRST symmetry [25].
The tree level gluon propagator obtained from the RGZ action (20) reads It is worth mentioning that the infrared behaviour of the RGZ gluon propagator (21) and of the corresponding ghost two-point function turns out to be in remarkable agreement with the recent numerical lattice simulations obtained on huge lattices [26,27,28]. Therefore, a numerical estimate of the non-perturbative parameters (m, µ, γ) can be obtained by fitting the lattice data by means of expression (21), see [28]. This leads to the presence of complex poles in the gluon propagator (21), as well as to a violation of reflection positivity, precluding thus a physical particle interpretation. As a consequence, gluons cannot belong to the physical spectrum of the theory. We see thus that the restriction to the Gribov region Ω captures non-trivial aspects of the gluon confinement.
Till now, the RGZ action has allowed for a variety of successful applications like: estimate of the masses of the first glueball states [29,30], yielding results which display the right mass hierarchy as observed in the available numerical simulations and whose accuracy is comparable to other non-perturbative approaches to the glueball spectrum (cf. e.g. [31] for a review), inclusion of quarks and estimate of the masses of meson states [32], study of the Casimir energy [33], finite temperature effects [34,35,36,37], study of the confinement/deconfinement transition in presence of Higgs fields [38,39], analysis of the relevance of the Gribov issue in supersymmetric theories [40,41].
The feature that we want to explore in the present work is the fact that both the GZ and the RGZ tree-level propagators hold the key for the good UV behavior of the theory. More precisely we note that the propagator D(k 2 ) (22) can be written as where The first term in (23) represents the usual propagator of a massive vector boson. The second term is the contribution coming from the restriction to the Gribov region. Notice the negative sign that points to an unphysical contribution that violates positivity requirements. The important feature we want to emphasize is the subleading contribution of the second term in the U V : it presents a ∼ 1/k 4 suppression with respect to the standard first term, which will always produce a UV convergent loop contribution in dimension 4. The renormalization of the RGZ and GZ (corresponding to µ = m = 0) follows from this important property and, as already mentioned, it is well known that γ does not renormalize independently and thus cannot be considered as an independent dynamically generated scale.
One is thus led to conjecture that this is a general property of theories displaying such confining propagators, with γ standing for a general mass scale associated with confinement of the fundamental fields; γ must be understood as a scale determined by other dynamically generated scales of the theory. More precisely, the second term in (23) cannot generate any new U V divergences in the theory and therefore cannot change the renormalization properties of the theory, which must be the same as with γ = 0. In a diagrammatic approach, only positive powers of propagators appear, so that it is clear that the highly-suppressed Gribov contribution (cf. (23), e.g.) will not influence the deep UV behavior of the theory. Furthermore, it follows that if the theory with γ = 0 does not generate a mass scale, then, since there can be no divergences proportional to γ, no mass scale will be generated in the γ = 0 theory. This in turn means that it is not possible to assign a dynamical meaning to the parameter γ in this case, i.e., the only possible solution is to have γ = 0 in these cases. An example of a theory displaying this feature is N = 4 supersymmetric Yang-Mills which, due to its conformal character, has vanishing β-function to all orders. As a consequence of the absence of a scale, it turns out that γ = 0, meaning that no mass scale associated to the Gribov copies is generated [41].
In the following sections we will study a variety of examples that support these claims. In section 2 we discuss the case of an interacting scalar field theory displaying a confining propagator. In section 3 we consider the inclusion of confined fermions interacting with the confined scalars through a Yukawa term. In section 4 we discuss the case of Super Yang-Mills with N = 1 supersymmetries and show to all orders via the algebraic renormalization approach that the adoption of Gribov-type propagators does not produce any new UV divergences, with the renormalization of the IR parameters being completely defined by the UV renormalization of the parameters of the original theory. Section 5 collects our summary and conclusions.

Interacting scalar fields with confining propagators
Consider the theory of a real scalar field φ defined by the following action in D = 4 euclidean space where where m 2 is the mass of the scalar field in the deconfined (γ → 0) theory and λ is the quartic coupling. Here, γ is the confining parameter that shall play a similar role for the scalars as the Gribov mass does for the confined gluons. Our claim in this case is that the presence of the IR parameter γ does not affect the deep UV behavior of the theory at all.
The quadratic part of the total action furnishes the tree-level confining propagator for the scalar fields: where we have isolated the confining contribution to the scalar propagator, γ 4 ∆, with which is highly suppressed in the UV: ∆ ∼ 1/k 6 . The mass parameters M 2 ± are written in terms of γ and m being complex for large enough γ/m. The complexity of these IR mass parameters is closely related to positivity violation and the absence of a physical particle interpretation for these excitations, in line with confinement.
It is straightforward to see that there are no new U V divergences associated with the term S γ (28) by looking at the diagrams of primitive divergences of the theory.
In fact, the one-loop scalar selfenergy is The correction to the quartic coupling at one loop reads: As a representative example at two-loop order, we may look at the scalar selfenergy sunset diagram: In all examples above, the appearance of a general form for the contributions of the confining scale with increasingly UV convergent momentum integrals is clear. It is straightforward to realize then that this pattern will spread throughout all orders of the diagrammatic expansion, so that we are led to infer that contributions proportional to γ cannot give rise to new primitive divergences, besides the ones coming from S s (26) alone, i.e. the original theory.

Confined fermions and scalars with Yukawa interaction
The same reasoning can be applied when Dirac fermions are added to the theory, with an Yukawa coupling and a fermionic Gribov-type term rendering the fermionic excitations also confined.
We consider here the theory in the absence of scalar condensates. In this case, the full action where where M is the mass of the original fermion field (i.e. for Γ → 0) and g is the Yukawa coupling.
In the fermionic sector the IR mass scale analogous to the Gribov parameter is Γ.
Analogously to the purely scalar case, it is easily seen that there are no U V divergences associated to the whole term S γ,Γ (39). The scalar excitations display the confining propagator of the last section, (29), while for the confining fermion propagator, we have Again, the isolated confining contribution to the propagator is highly suppressed in the UV with respect to the standard massive Dirac term (∼ 1/k): and we anticipate that the primitive divergences of the theory with confined propagators will be exactly the ones coming from S s (26) and S f (37) alone, since any contribution proportional to γ or Γ will be strongly suppressed in the UV.
At one loop order, besides the diagrams already analyzed in the previous section, new diagrams contributing to primitive divergences appear, due to the presence of fermion lines (dashed ones): It should be noticed that the Yukawa coupling breaks the discrete symmetry φ → −φ originally present in the scalar sector, generating at the quantum level a cubic scalar interaction. This means that the renormalizable version of this theory requires a counterterm for the cubic scalar interaction, even if the physical value of this coupling is set to zero. In the case of a pseudoscalar Yukawa coupling (i.e. gφψψ → gφψγ 5 ψ), parity symmetry guarantees that the cubic terms vanish identically. We emphasize, however, that our statement concerning the UV properties of Gribovtype confining propagators remains valid in any case, as will be made explicit below via the whole set of primitive divergences at one loop order.
In order to investigate the influence of the confining propagators in the UV regime, we may isolate the free fermion and scalar propagators from the confining contributions, namely Σ(k) (a) the one-loop fermion self energy: (42) (b) the fermion loop contributing to the scalar self energy: (c) the triangular diagram contributing to the scalar cubic interaction: (d) the fermion loop correction to the φ 4 vertex: (e) the modification of the Yukawa coupling: As already occurred for the confining scalar theory in the previous section, the highly suppressed UV behavior of the confining pieces Σ(k) U V ∼ 1/k 4 and ∆(k 2 ) U V ∼ 1/k 6 enforces the convergence of all terms proportional to the new massive parameters introduced (γ and Γ). The divergent integrals in all diagrams above are exactly the ones coming from the original action, i.e. S s + S f + S int (cf. (35)). In the theory including the confining quadratic nonlocal terms in S γ,Γ , the absence of new primitive divergences then guarantees that the parameters γ and Γ can be consistently related to dynamically generated scales and do not affect the UV regime of the theory.
Realizing that any diagrammatic expression at higher loops will involve higher powers of the propagators, it becomes straightforward to envision the generalization of our claim in the full diagrammatic expansion of this general Yukawa theory. Therefore, given the renormalizability of the original theory, one concludes that the resulting action with confining, Gribov-type propagators is renormalizable and the IR confining parameters in both fermionic and bosonic sectors do not display an independent renormalization, being thus consistent with dynamically generated mass scales.

N = 1 Super Yang-Mills in Wess-Zumino gauge within the Gribov-Zwanziger approach
Let us now investigate a more intricate theory with confining propagators, including gauge interactions as well as Majorana fermions. We consider here Yang-Mills theory in D = 4 spacetime dimensions with N = 1 supersymmetry in the presence of the Gribov horizon. We shall use this (most complicated) example to prove, to all-orders in the loop expansion, our claim concerning the good UV behavior of Gribov-type propagators. The IR parameters introduced will be shown to have renormalization parameters that are completely determined by the renormalization of the original theory.
This theory has already been put forward and investigated in Ref. [40]. There, the extension of the Gribov-Zwanziger framework to N = 1 Super-Yang-Mills (SYM) theories quantized in the Wess-Zumino gauge by imposing the Landau gauge condition was presented. The resulting effective action is Here, the operator Q stands for the generalized BRST operator which encodes both gauge and supersymmetry transformations 2 ; S γ is the horizon term in its local form, eq.(15), namely and the term S G is given by which also has a new massive constant M . This quantum action takes into account the existence of Gribov copies in the path-integral quantization of the theory. It encodes the restriction to the first Gribov horizon while keeping full compatibility with non-perturbative supersymmetric features, such as the exactly vanishing vacuum energy.
Even though this non-perturbative framework has been constructed through the introduction of two massive parameters γ, M which are not present in the classical action, those new parameters are determined in a dynamical, self-consistent way via two non-perturbative conditions: (i) the Gribov gap equation, that fixes γ by imposing the positivity of the Faddeev-Popov operator and eliminating a large set of Gribov copies from the functional integral, and (ii) the vanishing of the vacuum energy, which determines the parameter M that plays the role of a supersymmetric counterpart of the Gribov parameter γ, guaranteeing a consistent non-perturbative fermion sector. Interestingly, the appearance of the dynamical fermionic scale M has been shown to be directly related to the formation of a gluino condensate, a well-known non-perturbative property of N = 1 SYM theories. For further details, the reader is referred to Ref. [40]. A brief summary of the notation adopted may also be found in the Appendix B.
The propagators of the theory (47) can be straightforwardly shown to be of the Gribov type. The gauge field propagator is: which, apart from the more complicated tensorial structure, is equivalent to the Gribov scalar propagator studied above in section 2. The gauge field propagator in this Gribov-extended N = 1 SYM theory displays thus a confining contribution that is suppressed by an extra 1/p 4 factor in the UV as compared to the free term.
For gluino fields we have: where C αβ is the charge conjugation matrix and The presence of three two-point correlation functions involving gluino fields is a result of the lack of charge conservation for Majorana fermions. One verifies however that all of them have the form of Gribov propagators with M playing an analogous role as the Gribov parameter in the gluino sector. In particular, one can easily check that the same structure observed for the Gribov fermion propagator in the previous section (cf. Eq. (40)) is found here: where the isolated confining contribution Σ λ to the gluino propagator is again highly suppressed in the UV with respect to the leading term (∼ 1/k): The same reasoning applied in the scalar and Yukawa theories above may be followed here in order to prove that the UV regime of the theory remains the same even after the inclusion of nonlocal confining terms in the propagators. One may compute the one-loop primitive divergences and show that the confining parameters γ, M will not affect the UV divergent pieces, due to the high suppression observed in the Gribov-type propagators. We shall, however, use this most complicated theory analyzed in the current section to present an all-order algebraic proof of renormalizability and of the fact that the confining parameters γ, M do not display independent renormalization.
The non-local action (47) is, however, not helpful in the algebraic renormalization procedure. Fortunately we are able to write its local form with the insertion of auxiliary fields.
The whole action which describes our model can then be written in its local form as, where the set of auxiliary fields (θ aα , θ aα ,ζ aα , ζ aα ) has an analogous role of the set of auxiliary localizing fields (φ a µ , ϕ a µ ,ω a µ , ω a µ ) introduced by Zwanziger in GZ model, i.e. it allows to describe the non-local gluino term (49) in a local fashion.
Applying the algebraic renormalization procedure to the local action (57) above we are able to prove that: (i) the Gribov-extended SYM theory is renormalizable; and (ii) the massive parameters γ, M introduced in the infrared action do not renormalize independently, meaning that they are consistent with dynamically generated mass scales, produced by nonperturbative interactions in the original theory. All details of the proof may be found in the Appendix A.
The final results for the renormalization factors related to the confining parameters M, γ may be read from the renormalization of external sources conveniently introduced in the algebraic procedure (cf. Appendix A). The renormalization of the sources M andM give us the renormalization factor of the Gribov parameter γ 2 , while the renormalization of V andV give us the renormalization of M 3/2 , when every source assumes its physical value stated at (64). We have: which clearly prove that the renormalization of the infrared parameters M, γ is fixed by the renormalization factor of the original SY M theory: the renormalization of the gauge coupling, Z g , the wave function renormalization of the gauge field, Z A , and and the wave function renormalization of the gluing field, Z λ .
Therefore we conclude that this action is indeed a suitable nonperturbative infrared action for N = 1 SYM theories, reducing consistently to the ultraviolet original action. Moreover, even in this very intricate non-Abelian gauge theory with matter fields, the good UV behavior in the presence of confining propagators of the Gribov type shows up at all orders.

Conclusion
In this paper we have studied the UV behaviour of quantum field theory models in which the two-point correlation functions of the elementary fields are described by confining propagators of the Gribov type.
Relying on the decompositions (29), (40), we have been able to show that the UV divergent behaviour of the Feynman diagrams is not affected by the infrared parameters, e.g. (γ, Γ), encoded in the aforementioned confining propagators.
From this property, it follows that no new UV divergences in the infrared parameters can arise. Otherwise said, the only UV divergences affecting the 1PI Green's functions of the theory are those present when the infrared parameters are set to zero. As a consequence, the infrared parameters do not renormalize independently, as explicitly shown in the case of N = 1 supersymmetric Yang-Mills theory.
In particular, in the case of a generic non-Abelian theory, the implementation of the restriction of the domain of integration to the Gribov region Ω has no consequences on the UV renormalisation properties of the theory.
In conclusion, the main result of the present work can be stated as follows: given a multiplicatively renormalizable Faddeev-Popov action, adding a Gribov horizon term in both gluon and matter sectors will not affect the ultraviolet properties of the theory at all. The resulting action remains multiplicatively renormalizable, with the same counterterms as the original theory. Moreover, the IR parameters originated by the horizon terms do not renormalize independently, being thus consistent with dynamically generated mass scales.

Construction of a complete invariant action
In order to prove the renormalizability of the action (57) introduced in Sect.4, we follow the procedure already employed in [7,10,21,22,23,40] and embed the action (57) into a more general one displaying a huge set of symmetries and Ward identities. In the present case, it turns out that the action (57) can be recovered as a particular case of the following expression Let us proceed by specifying the various terms appearing in the action (59). The first term, S N =1 SY M , is the N = 1 Euclidean Super Yang-Mills action with Majorana fermions in the Wess-Zummino gauge, without matter fields, namely, with D ab µ = δ ab ∂ µ − gf abc A c µ denoting the covariant derivative in the adjoint representation of the SU (N ); λ α being a four component Majorana spinor; and D a standing for a dimension two auxiliary fields needed in order to close the algebra of N = 1 supersymmetry. Also,λ = λ T C with C being the charge conjugation matrix, which is defined together with the Euclidean gamma matrices γ µ in Appendix B.
The second term in eq.(59) is the gauge fixing term in the Landau gauge, given by In this term, b a is the Lagrange multiplier enforcing the Landau gauge condition, ∂ µ A a µ = 0, while (c a ,č a ) are the Fadeev-Popov ghost fields and¯ = T C is a constant ghost. As we shall see later, the constant ghost¯ is needed to encode the supersymmetric transformations into a unique generalized BRST operator.
The third term of eq.(59) corresponds to the local and invariant Gribov-Zwanziger term 3 being given by In the expression above, (M ab µν , M ab µν ,Ñ ab µν , N ab µν ) are external sources which will be set equal to their physical values after the renormalization procedure, i.e. after removing the UV divergencies.
Finally, the fourth and last term of the action (59) is the local and invariant confining term for the gluino sector. This term is the analogous of the GZ term (62). It can be seen as the supersymmetric counterpart of expression (62). It reads Analogously to the term (62), this term depends on the external sources (V ab αβ , V ab αβ ,Û ab αβ , U ab αβ ).
As already mentioned, the original action (57) can be re-obtained from expression (59) when the external sources attain the following physical values: We see thus that the action (57) is a particular case of the more general expression (59). Therefore, we will turn our attention to the action (59), keeping in mind that we can always go back to the action (57) by taking the limit (64).
The advantage of working with the most general action (59) is that it is left invariant by both SUSY and BRST transformations, which can be embedded into a unique generalized BRST operator, see [40,42], where s is the usual BRST operator and δ α are the SUSY generators, with α being the constant ghost. More precisely, it turns out that where the action of the operator Q on the fields and external sources is defined as The operator Q has the pleasant property that its square gives the generator of the translations [40,42], i.e. where Q 2 is defined as the translation operator Notice that the transformations QA a µ , Qλ aβ , QD a and Qc a are non-linear in the quantum fields, meaning that they have to be treated as composite operators. Therefore, following the algebraic renormalization procedure [24], we introduce them into the starting action coupled to suitable external sources (K a µ , L a , T a , Y aα ) : with the following transformations Finally, for the complete starting invariant action suitable to study the symmetry content and renormalizability of the theory, we have where, following [7,10,21,22,23,40], we have introduced the composite index notation: and θ a α ,θ a α , ζ a α ,ζ a α → θ I ,θ I , ζ I ,ζ I , with (a, µ) ≡ i, j, k, l, · · · ∈ {1, . . . , 4(N 2 − 1)} , (a, α) ≡ I, J, K, L, · · · ∈ {1, . . . , 4(N 2 − 1)} .
A summary of all indices used here can be found in Appendix B, where we also display the quantum numbers, i.e. mass dimensions and charges of all fields and sources of the model, see Tables 1 and 2. According to those Tables, the fields and sources denoted with a (˜) have charge q f equal to (−1), while the fields and sources denoted with a (ˆ) have charge q f equal to (−1). That is why in Sect.4, a slightly different notation with respect to that used in the Introduction has been adopted.

A.2 Symmetry content of the model
It turns out that the the complete action (74) displays a huge set of Ward identities, which we enlist below:

Ward identities
• The Slavnov-Taylor identity: • The gauge-fixing condition and anti-ghost equation: • The equations of motion of the auxiliary fields: • The identities in the external BRST sources: • The U (f = 4(N 2 − 1)) invariance: The trace of this symmetry defines a q f charge and the composite index (a, µ) ≡ i, j, k, l, . . . .
• The ghost equation: • The equation of the source T a : Let us also remark that some of the identities enlisted above are linearly broken, i.e. they display a breaking term which is linear in the quantum fields. Such a breaking is a classical breaking, not affected by the renormalization process [24]. order in the parameter η,. This gives the following constraints for the counterterm Σ count : In particular, thanks to the property B Σ B Σ = ∇, the general solution of the constraint B Σ Σ count = 0, i.e. the first eq. of (101), can be written as wherea 0 is a free parameter and ∆ (−1) is an integrated polynomial in the fields and sources of dimension 3, ghost number −1, and q f = q f = 0. Taking into account the remaining constraints and the discrete symmetries (97) and (100), it follows that the most general expression for ∆ (−1) turns out to be with a 1 , a 2 free parameters. We see therefore that Σ count depends on three arbitrary coefficients, i.e. (a 0 , a 1 , a 2 ).

A.4 Renormalization factors
In order to complete the analysis of the algebraic renormalization, we still need to show that the counterterm Σ count can be reabsorbed into the starting action Σ through a redefinition of the fields {φ}, sources {S} and parameters {p} = {g, }, namely, where (φ 0 , S 0 , p 0 ) are the so-called bare quantities, defined through the renormalization factors as where with the {z} being linear combinations of the coefficients (a 0 , a 1 , a 2 ). Moreover, in the present case, a little care has to be taken due to the potential mixing of quantities which have the same quantum numbers, see also [40,42]. In fact, as it can be checked from Table 1 and Table 2, one notices that the field λ aα and the combination γ 5 T a have the same dimension and quantum numbers as well as the field D a and the combination (Y a −Û Ia Cθ I − U Iaθ I )γ 5 . As a consequence, these quantities can mix at quantum level, a well known property of renormalization theory. This feature can be properly taken into account by writing the renormalization of the fields λ and D as λ aα and while the remaining fields, sources and parameters still obey (106).
By direct inspection of eq.(105), we find All other remaining renormalization factors can be expressed in terms of the tree independent quantities (Z A , Z λ , Z g ), namely In particular, from eq.(64), it follows that the renormalization of the sources M andM gives us the renormalization factor of the Gribov parameter γ 2 , while the renormalization of V andV yields the renormalization of M 3/2 , i.e.
The other sources renormalize as Finally, the renormalization factor of the supersymmetric ghost parameter is while we also have This concludes the proof of the algebraic renormalizability of the model.

B Notations and conventions in Euclidean space-time
Units: = c = 1.
Wick rotations: Gamma matrices: Pauli matrices: The gamma matrices obey the following relations: We also define the matrix γ 5 as with the following properties: The charge conjugation matrix is with The matrices σ µν are defined as (σ µν ) β α ≡ and have the property σ † µν = −σ µν .
Table of quantum numbers: We display below the quantum numbers of the fields and sources appearing in the action (74). Notice that by "nature" we mean "C" for commuting (or bosonic) and "A" for anti-commuting (or fermionic).