A non-perturbative study of matter field propagators in Euclidean Yang-Mills theory in linear covariant, Curci-Ferrari and maximal Abelian gauges

TL;DR

This paper investigates matter field propagators in Euclidean Yang-Mills theory using the Refined Gribov-Zwanziger approach, incorporating Gribov copies effects and providing explicit propagator calculations across various gauges with qualitative lattice agreement.

Contribution

It extends the Refined Gribov-Zwanziger framework to matter fields, deriving explicit propagators in multiple gauges considering Gribov effects.

Findings

Qualitative agreement with available lattice data

Explicit propagator expressions for scalar and quark fields

Incorporation of Gribov copies in matter sector

Abstract

In this work, we study the propagators of matter fields within the framework of the Refined Gribov-Zwanziger theory, which takes into account the effects of the Gribov copies in the gauge-fixing quantization procedure of Yang-Mills theory. In full analogy with the pure gluon sector of the Refined Gribov-Zwanziger action, a non-local long-range term in the inverse of the Faddeev-Popov operator is added in the matter sector. Making use of the recent BRST invariant formulation of the Gribov-Zwanziger framework achieved in [Capri et al 2016], the propagators of scalar and quark fields in the adjoint and fundamental representations of the gauge group are worked out explicitly in the linear covariant, Curci-Ferrari and maximal Abelian gauges. Whenever lattice data are available, our results exhibit good qualitative agreement.

Equations

$$S_{\mathrm{FP}}=\int {\text{d}}^{d}x\left(\frac{1}{4}{F}_{\mu \nu }^a{F}_{\mu \nu }^a+b^a{\partial }_{\mu }{A}_{\mu }^a+{\bar{c}}^a{\partial }_{\mu }{D}_{\mu }^{ab}(A)c^b\right),$$

$$F_{\mu \nu }^a$$

$$D_{\mu }^{ab}$$

$$A_{\mu }^{\prime a}=A_{\mu }^a-D_{\mu }^{ab}{\xi }^b,$$

$${\partial }_{\mu }{A}_{\mu }^{\prime a}={\partial }_{\mu }{A}_{\mu }^a-{\partial }_{\mu }{D}_{\mu }^{ab}{\xi }^b=-{\partial }_{\mu }{D}_{\mu }^{ab}{\xi }^b\ne 0,$$

$$\Omega=\left\{A^{a}_{\mu}\,,\,\partial_{\mu}A^{a}_{\mu}=0\,\Big{|}\,\mathcal{M}^{ab}(A)>0\right\}\,,$$

$$\mathcal{Z}={\int }_{\Omega }\left[\mathcal{D}A\right]\delta ({\partial }_{\mu }{A}_{\mu }^a)|\det ({\mathcal{M}}^{ab})|{\mathrm{e}}^{-S_{\mathrm{YM}}}=\int \left[\mathcal{D}A\right]\left[\mathcal{D}c\right]\left[\mathcal{D}\bar{c}\right]\left[\mathcal{D}b\right]{\mathrm{e}}^{-S_{\mathrm{FP}}-{\gamma }^{4}H(A)+dV{\gamma }^4(N^2-1)},$$

$$H(A)=g^2\int {\text{d}}^{d}x{\text{d}}^{d}y{f}^{abc}{A}_{\mu }^b(x){\left[{\mathcal{M}}^{-1}(A)\right]}^{ad}(x,y)f^{dec}{A}_{\mu }^e(y).$$

$$⟨H(A)⟩=dV(N^2-1),$$

$$S_{\mathrm{GZ}}=S_{\mathrm{FP}}-\int {\text{d}}^{d}x\left({\bar{\phi }}_{\mu }^{ac}{\mathcal{M}}^{ab}{\phi }_{\mu }^{bc}-{\bar{\omega }}_{\mu }^{ac}{\mathcal{M}}^{ab}{\omega }_{\mu }^{bc}+g{f}^{adl}{\bar{\omega }}_{\mu }^{ac}{\partial }_{\nu }\left({\phi }_{\mu }^{lc}{D}_{\nu }^{de}{c}^e\right)\right)+{\gamma }^2\int {\text{d}}^{d}xg{f}^{abc}{A}_{\mu }^a(\phi +\bar{\phi }{)}_{\mu }^{bc},$$

$$\mathcal{Z}=\int \left[\mathcal{D}{\mu }_{\mathrm{GZ}}\right]{\mathrm{e}}^{-S_{\mathrm{GZ}}+dV{\gamma }^4(N^2-1)},$$

$$\left[\mathcal{D}{\mu }_{\mathrm{GZ}}\right]=\left[\mathcal{D}A\right]\left[\mathcal{D}c\right]\left[\mathcal{D}\bar{c}\right]\left[\mathcal{D}b\right]\left[\mathcal{D}\phi \right]\left[\mathcal{D}\bar{\phi }\right]\left[\mathcal{D}\omega \right]\left[\mathcal{D}\bar{\omega }\right].$$

$$s{A}_{\mu }^a$$

$$s{\bar{c}}^a$$

$$s{\phi }_{\mu }^{ab}$$

$$s{\bar{\omega }}_{\mu }^{ab}$$

$$s{S}_{\mathrm{GZ}}={\gamma }^2\int {\text{d}}^{d}xg{f}^{abc}\left(-D_{\mu }^{ad}{c}^d(\phi +\bar{\phi }{)}_{\mu }^{bc}+A_{\mu }^a{\omega }_{\mu }^{bc}\right).$$

$$⟨A_{\mu }^a{A}_{\mu }^a⟩\propto {\gamma }^2,⟨{\bar{\phi }}_{\mu }^{ab}{\phi }_{\mu }^{ab}-{\bar{\omega }}_{\mu }^{ab}{\omega }_{\mu }^{ab}⟩\propto {\gamma }^2.$$

$$S_{\mathrm{RGZ}}=S_{\mathrm{GZ}}+\frac{m^2}{2}\int {\text{d}}^{d}x{A}_{\mu }^a{A}_{\mu }^a-M^2\int {\text{d}}^{d}x\left({\bar{\phi }}_{\mu }^{ab}{\phi }_{\mu }^{ab}-{\bar{\omega }}_{\mu }^{ab}{\omega }_{\mu }^{ab}\right),$$

$$⟨A_{\mu }^a(k)A_{\nu }^b(-k){⟩}_{d=3,4}={\delta }^{ab}\frac{k^2+M^2}{(k^2+m^2)(k^2+M^2)+2g^{2}N{\gamma }^4}\left({\delta }_{\mu \nu }-\frac{k_{\mu }{k}_{\nu }}{k^2}\right),$$

$$⟨A_{\mu }^a(k)A_{\nu }^b(-k){⟩}_{d=2}={\delta }^{ab}\frac{k^2}{k^4+2g^{2}N{\gamma }^4}\left({\delta }_{\mu \nu }-\frac{k_{\mu }{k}_{\nu }}{k^2}\right).$$

$${\partial }_{\mu }{A}_{\mu }^a=\alpha b^a,$$

$${\mathcal{M}}_{\mathrm{LCG}}^{ab}(A){\xi }^b=-{\delta }^{ab}{\partial }^2{\xi }^b+g{f}^{abc}{A}_{\mu }^c{\partial }_{\mu }{\xi }^b+g{f}^{abc}({\partial }_{\mu }{A}_{\mu }^c){\xi }^b=0,$$

$$A_{\mu }^h=\left({\delta }_{\mu \nu }-\frac{{\partial }_{\mu }{\partial }_{\nu }}{{\partial }^2}\right)\left(A_{\nu }-ig\left[\frac{1}{{\partial }^2}\partial A,A_{\nu }\right]+\frac{ig}{2}\left[\frac{1}{{\partial }^2}\partial A,{\partial }_{\nu }\frac{1}{{\partial }^2}\partial A\right]+\mathcal{O}(A^3)\right),$$

$${\tilde{S}}_{\mathrm{GZ}}^{\mathrm{LCG}}=S_{\mathrm{YM}}+\int {\text{d}}^{d}x\left(b^{h,a}{\partial }_{\mu }{A}_{\mu }^a-\frac{\alpha }{2}{b}^{h,a}{b}^{h,a}+{\bar{c}}^a{\partial }_{\mu }{D}_{\mu }^{ab}(A)c^b\right)+{\gamma }^{4}H(A^h),$$

$$H(A^h)=g^2\int {\text{d}}^{d}x{\text{d}}^{d}y{f}^{abc}{A}_{\mu }^{h,b}(x){\left[{\mathcal{M}}^{-1}(A^h)\right]}^{ad}(x,y)f^{dec}{A}_{\mu }^{h,e}(y),$$

$${\mathcal{M}}^{ab}(A^h)=-{\delta }^{ab}{\partial }^2+g{f}^{abc}{A}_{\mu }^{h,c}{\partial }_{\mu },\mathrm{with}{\partial }_{\mu }{A}_{\mu }^{h,a}=0.$$

$$S_{\mathrm{GZ}}^{\mathrm{LCG}}$$

$$s{A}_{\mu }^a$$

$$s{\bar{c}}^a$$

$$s{\phi }_{\mu }^{ab}$$

$$s{\bar{\omega }}_{\mu }^{ab}$$

$$s{A}_{\mu }^{h,a}$$

$$s{S}_{\mathrm{GZ}}^{\mathrm{LCG}}=0.$$

Full text

1. I Introduction

2. II Overview of the non-perturbative BRST invariant formulation of the Refined Gribov-Zwanziger framework

3. II.1 The Gribov problem in the Landau gauge

4. II.2 Going beyond the Landau gauge and the non-perturbative BRST symmetry: linear covariant gauges

5. II.3 Curci-Ferrari gauge

6. II.4 Maximal Abelian gauge (MAG)

7. III Non-perturbative coupling of scalar fields in the adjoint representation

8. III.1 Linear covariant and Curci-Ferrari gauges

9. III.2 Maximal Abelian gauge

10. IV Generalization of the non-perturbative matter coupling for quark fields

11. IV.1 Linear covariant and Curci-Ferrari gauges

12. IV.2 Maximal Abelian gauge

13. V Conclusions

14. A Conventions in Euclidean space

I Introduction

While quantum chromodynamics (QCD) is well understood at high energies, where perturbation theory is reliable due to asymptotic freedom, the low energy sector remains a challenging open problem in theoretical Physics. In the infrared region, perturbation theory breaks down and non-perturbative techniques are needed. A full control of the infrared regime of QCD would provide a fundamental understanding of the confinement of quarks and gluons, a goal not achieved till now.

Different approaches which take into account non-perturbative effects in QCD were devised in the last decades, see Alkofer:2000wg ; Binosi:2009qm ; Greensite:2011zz ; Brambilla:2014jmp ; Deur:2016tte . Up to now, the interplay of such approaches was able to produce non-trivial results. Though, a complete consistent picture of the mechanism behind colour confinement is still lacking.

One approach to deal with the confinement problem is the non-perturbative study of the correlation functions of the theory. Functional techniques based on the Dyson-Schwinger equations and on the functional renormalization group as well as numerical lattice simulations have been employed in the analysis of the infrared behavior of the correlation functions. In particular, the two-point gluon correlation function has been object of very intensive investigations. In fact, the infrared structure of two-point gluon correlation function, e.g. the gluon propagator, turns out to encode important features which are interpreted as signals of confinement. For instance, lattice numerical simulations as well as computations based on the Dyson-Schwinger equations show that the gluon propagator exhibits a violation of the reflection positivity. As such, it cannot be associated with a physical excitation of the spectrum of the theory. This property is interpreted as a manifestation of confinement, see vonSmekal:1997ohs ; vonSmekal:1997ern ; Cucchieri:2004mf ; Dudal:2013yva ; Cornwall:2013zra . One has to keep in mind that the gluon propagator is a gauge-dependent quantity. Nevertheless, it still contains important information about such (un)physical elementary fields, being the simplest correlation function one might compute. In the last decade, the gluon propagator has been studied in great detail in the Landau gauge, due to its special features, namely the transversality of the propagator itself and the important property of having a useful lattice formulation which has allowed for a numerical study of the gluon propagator on large lattices. More precisely, the most recent lattice simulations point towards an infrared suppressed gluon propagator which attains a finite value at zero momentum in four and three space-time dimensions, while it vanishes at zero momentum in two space-time dimensions. One says that in three and four dimensions the gluon propagator is of decoupling/massive type, while in two dimensions it is of scaling type, see Maas:2011se ; Cucchieri:2007rg ; Aguilar:2008xm ; Fischer:2008uz ; Maas:2008ri ; Cucchieri:2009zt ; Maas:2007uv .

Besides the aforementioned functional and numerical lattice approaches, an analytical framework which takes into account the existence of the Gribov copies Gribov:1977wm occurring in the Faddeev-Popov quantization of gauge theories has received increasing interest in the recent years. The so-called Refined Gribov-Zwanziger setup captures the effects of the spurious gauge copies as well as of additional non-perturbative effects related to the existence of dimension two-condensates, giving rise to an effective infrared action, the Refined Gribov-Zwanziger action, yielding a gluon propagator of the decoupling type which is in very good agreement with the most recent lattice data in both four and three space-time dimensions. In two dimensions, infrared singularities forbid the formation of the dimension two-condensates and the refinement does not take place. As a consequence, the gluon propagator turns out to be of the scaling type. In this paper, we focus on the Refined Gribov-Zwanziger formulation. An extensive review of the developments of this framework is presented in Sect. II.

Nonetheless, in QCD, in addition of the pure gluon sector, one has to face also the complex issue of quark confinement, to which different strategies have been devoted, see Greensite:2011zz ; Brambilla:2014jmp . As far as the Refined Gribov-Zwanziger framewrok is concerned, a possible mechanism to take into account matter confinement was proposed in Capri:2014bsa ; Dudal:2013vha ; Palhares:2016wqn ; Capri:2016aqq in the Landau gauge. More precisely, as it will be reviewed in Sect. II, within the Refined Gribov-Zwanziger approach in the Landau gauge, the non-perturbative effect of the Gribov copies is accounted for by restricting the domain of integration in the functional integral to a certain region $\Omega$, called the Gribov region, which is defined by demanding that the Faddeev-Popov operator $\mathcal{M}(A)$ is strictly positive, so that it is invertible within $\Omega$. Such a restriction enables us to eliminate a large set of copies. In practice, the restriction to $\Omega$ is achieved by adding to the starting Faddeev-Popov action an additional non-local term, known as the horizon function, which contains the inverse of the operator $\mathcal{M}(A)$. It is precisely the addition of this additional long-range term which is responsible for the infrared modifications of the gluon propagator, which turns out to be a confining propagator, exhibiting complex poles and lacking the Källén-Lehmann representation. Remarkably, the non-local horizon function can be cast in local form through the addition of a suitable set of auxiliary fields. The resulting action is multiplicatively renormalizable to all orders.

The proposal made in Capri:2014bsa ; Dudal:2013vha ; Palhares:2016wqn ; Capri:2016aqq consists in generalizing the introduction of the non-local horizon function to the matter sector, in complete analogy with the gluon sector. Moreover, as in the gluon sector, the non-local matter coupling term can be cast in local form, giving rise to a fully local and renormalizable action. In Capri:2014bsa ; Dudal:2013vha ; Palhares:2016wqn ; Capri:2016aqq , this prescription was implemented for scalar fields in the adjoint representation of the gauge group and for spinor fields in the fundamental representation. The whole procedure was carried out in the Landau gauge, for which the Refined Gribov-Zwanziger setup was well established, see also Capri:2015pfa and references therein for the maximal Abelian gauge.

Very recently, the Refined Gribov-Zwanziger framework has been extended to the class of the linear covariant and Curci-Ferrari gauges Capri:2015ixa ; Capri:2016aqq ; Capri:2015nzw ; Pereira:2016fpn ; Capri:2016gut , allowing, in particular, to establish the independence from the gauge parameter of the gauge invariant correlation functions as well of the poles of the transverse part of the gluon propagator. In the light of such developments, it seems natural to ask ourselves how matter fields should be coupled to the Refined Gribov-Zwanziger action in such gauges. This is precisely the aim of the present work.

The paper is organized is follows: Sect. II contains an overview of the Refined Gribov-Zwanziger action, covering the Landau, the linear covariant, the Curci-Ferrari as well as the maximal Abelian gauge. After that, in Sect. III, we describe the non-perturbative coupling of scalar fields in the adjoint representation of the gauge group within the Refined Gribov-Zwanziger action in the aforementioned gauges. Subsequently, in Sect. IV, we work out the non-perturbative coupling of quark fields and its consequences on the propagator. Finally, we collect our conclusions. To keep the paper self-contained as much as possible, we have added two appendices devoted to the details of our construction as well as to the conventions used throughout the paper.

II Overview of the non-perturbative BRST invariant formulation of the Refined Gribov-Zwanziger framework

In this section, we review the recently proposed BRST invariant formulation of the Refined Gribov-Zwanziger action in linear covariant Capri:2015ixa ; Capri:2015nzw , Curci-Ferrari Pereira:2016fpn and maximal Abelian gauges Capri:2015pfa . For the benefit of the reader, we start with a brief overview of the Gribov problem in the Landau gauge for which the Refined Gribov-Zwanziger action was originally constructed.

II.1 The Gribov problem in the Landau gauge

Let us consider Yang-Mills theory in $d$ Euclidean dimensions with $SU(N)$ gauge group quantized in the Landau gauge, namely $\partial_{\mu}A^{a}_{\mu}=0$. The Faddeev-Popov procedure results in the gauge-fixed action:

[TABLE]

with the field strength $F^{a}_{\mu\nu}$ and the covariant derivative $D^{ab}_{\mu}$ in the adjoint representation of the gauge group given by

[TABLE]

The parameter $g$ stands for the gauge coupling111The coupling $g$ is dimensionless in $d=4$. and $f^{abc}$ are the real and totally antisymmetric structure constants of the gauge group. The fields $(\bar{c}^{a},c^{a})$ denote the Faddeev-Popov ghosts, while $b^{a}$ is the Lagrange multiplier implementing the Landau gauge condition. Nevertheless, as investigated by Gribov in Gribov:1977wm , the action (1) is plagued by the existence of gauge copies, i.e. equivalent gauge configurations which still obey the gauge condition. This fact can be observed very concretely by considering a gauge field configuration $A^{a}_{\mu}$ which satisfies the Landau gauge condition and another gauge field $A^{\prime a}_{\mu}$ which is connected to $A^{a}_{\mu}$ by an infinitesimal gauge transformation:

[TABLE]

with $\xi^{a}$ being the infinitesimal gauge parameter of the transformation. If the Landau gauge condition were ideal, i.e. if it were selecting only one representative $A^{a}_{\mu}$ per gauge orbit222A gauge orbit of a given configuration $A^{a}_{\mu}$ is the set of all gauge fields related to $A^{a}_{\mu}$ by a gauge transformation., then $A^{\prime a}_{\mu}$ would not obey anymore the Landau gauge condition, $\partial_{\mu}A^{\prime a}_{\mu}\neq 0$. Therefore, from eq.(3) we should have

[TABLE]

where the condition $\partial_{\mu}A^{a}_{\mu}=0$ was employed. Hence, eq.(4) shows that if the Faddeev-Popov operator $\mathcal{M}^{ab}(A)\equiv-\partial_{\mu}D^{ab}_{\mu}(A)$ develops zero-modes, then the Landau gauge is not ideal. Gribov proved in Gribov:1977wm that the operator $\mathcal{M}^{ab}(A)\equiv-\partial_{\mu}D^{ab}_{\mu}(A)$ does exhibit in fact zero-modes. As a consequence, a residual gauge symmetry remains even after the implementation of (1). The existence of such spurious configurations known as Gribov copies is the so-called Gribov problem. For a pedagogical review of this subject, we refer to Sobreiro:2005ec ; Vandersickel:2012tz ; Vandersickel:2011zc ; Pereira:2016inn . Let us emphasize that the previous argument is restricted to Gribov copies generated by infinitesimal gauge transformations. Finite gauge transformations were considered in vanBaal:1991zw .

The aforementioned discussion on the existence of the Gribov problem might induce the reader to think that this is a particular pathology of the Landau gauge or of a subclass of gauges, which can be circumvented by a suitable choice of a more appropriate gauge condition. Nevertheless, it was proved by Singer in Singer:1978dk that this is not the case. In fact, it turns out that the Gribov problem has to do with the non-trivial topological structure of Yang-Mills theories.

In order to deal with the Gribov copies in the path integral measure, Gribov proposed in Gribov:1977wm the restriction of the path integral domain to a smaller region $\Omega$ in field space, known as the Gribov region, defined as

[TABLE]

namely, $\Omega$ is the set of gauge fields which satisfy the Landau gauge condition and for which the Faddeev-Popov operator is strictly positive. The boundary $\partial\Omega$, where the first vanishing eigenvalue of the Faddeev-Popov operator shows up, is called the first Gribov horizon. We should mention that, although the region $\Omega$ is free from infinitesimal Gribov copies, it still contains additional copies vanBaal:1991zw related to finite gauge transformations. A smaller region, contained inside $\Omega$ and known as the Fundamental Modular Region, exists and turns out to be fully free from Gribov ambiguities. However, unlike the Gribov region $\Omega$, a practical way to implement the restriction of the domain of integration in the functional integral to the Fundamental Modular Region has not yet been achieved so far. Therefore, we stick to the Gribov region $\Omega$ which displays important properties: i) it is bounded in all directions in field space; ii) it is convex; iii) all gauge orbits cross it at least once. These properties were proven in a rigorous fashion in Dell'Antonio:1991xt and give a well defined support to original Gribov’s proposal of restricting the functional integral to $\Omega$. Such restriction is effectively implemented through the addition of an extra term into the action (1), as shown independently by Gribov and Zwanziger, Gribov:1977wm ; Zwanziger:1989mf ; Capri:2012wx , i.e.

[TABLE]

where $H(A)$ is known as the horizon function, being given by

[TABLE]

In expression (6), $V$ is the space-time volume and $\gamma$ is a parameter with the dimension of a mass, known as the Gribov parameter. It is not a free parameter, being determined in a self-consistent way through the so-called horizon condition, i.e.

[TABLE]

where the expectation value $\langle\ldots\rangle$ is taken with respect to the modified measure given by eq.(6).

As is apparent from eq.(7), the presence of the inverse of the Faddeev-Popov operator makes the horizon function a non-local quantity. Nevertheless, it can be cast in local form by the introduction of a suitable set of auxiliary fields, namely a pair of commuting $(\varphi^{ab}_{\mu},\bar{\varphi}^{ab}_{\mu})$ and of anticommuting $(\omega^{ab}_{\mu},\bar{\omega}^{ab}_{\mu})$ fields. Written in terms of these new fields, the Gribov-Zwanziger action $S_{\mathrm{GZ}}$ is expressed as

[TABLE]

and eq.(6) takes the form

[TABLE]

with

[TABLE]

Remarkably, the action $S_{\mathrm{GZ}}$ is local and renormalizable to all orders in perturbation theory Zwanziger:1989mf . Hence, the Gribov-Zwanziger action is an effective framework which implements the restriction of the domain of integration in the path integral to the Gribov region $\Omega$ in a renormalizable and local way.

The Gribov-Zwanziger action has many interesting and non-trivial properties. For our present purposes, we focus on a few of them. First, the gluon propagator computed out of (9) is suppressed in the deep infrared regime and attains a vanishing value at zero-momentum, a result which is at odds with the divergent perturbative behavior. This propagator is said to be of the scaling type. Also, it violates reflection positivity and, therefore, gluons cannot be interpreted as excitations of the physical spectrum, being thus confined. The ghost propagator, however, is enhanced in the strong coupling regime, diverging as $1/k^{4}$ for $k\approx 0$. Another property of the Gribov-Zwanziger action (9) is that it breaks the BRST symmetry, given by the following transformations,

[TABLE]

in an explicit way, namely

[TABLE]

Being proportional to $\gamma^{2}$, the BRST breaking is a soft breaking. It becomes relevant in the non-perturbative infrared region. Though, it does not affect the deep ultraviolet region, so that the perturbative results are recovered.

More recently, it has been realized that the localizing fields $(\varphi,\bar{\varphi},\omega,\bar{\omega})$ develop their own dynamics and non-trivial additional effects are generated. In particular, it has been shown that dimension-two condensates, $\langle A^{a}_{\mu}A^{a}_{\mu}\rangle$ and $\langle\bar{\varphi}^{ab}_{\mu}\varphi^{ab}_{\mu}-\bar{\omega}^{ab}_{\mu}\omega^{ab}_{\mu}\rangle$, are dynamically generated Dudal:2007cw ; Dudal:2008sp ; Dudal:2011gd ; Gracey:2010cg , i.e.

[TABLE]

Taking into account the existence of such dimension-two condensates from the beginning, gives rise to the so-called Refined Gribov-Zwanziger action, which is expressed as

[TABLE]

where, much alike the Gribov parameter $\gamma^{2}$, the massive parameters $(m^{2},M^{2})$ are not independent, being determined by suitable gap equations obtained through the evaluation of the effective potential for the condensates $\langle A^{a}_{\mu}A^{a}_{\mu}\rangle$ and $\langle\bar{\varphi}^{ab}_{\mu}\varphi^{ab}_{\mu}-\bar{\omega}^{ab}_{\mu}\omega^{ab}_{\mu}\rangle$, see Dudal:2011gd .

The addition of the dimension-two operators, $A^{a}_{\mu}A^{a}_{\mu}$ and $(\bar{\varphi}^{ab}_{\mu}\varphi^{ab}_{\mu}-\bar{\omega}^{ab}_{\mu}\omega^{ab}_{\mu})$, does not spoil the renormalizability of the refined action (15). Notably, taking into account these additional non-perturbative effects, changes the gluon and ghost propagators. For instance, the gluon propagator displays now a decoupling/massive behavior, exhibiting a finite non-vanishing value at zero-momentum, while being still suppressed in the deep infrared sector. The ghost propagator, however, is not enhanced anymore in the strong coupling and, for $k\approx 0$, it behaves as $1/k^{2}$. Such behavior of the gluon and ghost propagator is in very good agreement with the most recent lattice simulations in the Landau gauge, see Cucchieri:2007rg ; Cucchieri:2008fc ; Maas:2008ri ; Cucchieri:2011ig ; Oliveira:2012eh ; Duarte:2016iko .

An interesting property of the refinement of the Gribov-Zwanziger action is that its occurrence depends on the space-time dimension $d$. In particular, for $d=3,4$, the formation of dimension-two condensates is dynamically favoured and the Gribov-Zwanziger action is naturally refined Dudal:2008sp ; Dudal:2008rm . Nevertheless, in $d=2$, infrared singularities prevent the introduction of such operators and the refinement does not take place. In particular, this implies that, for $d=3,4$, the gluon propagator is of decoupling type, while in $d=2$, it is of scaling type Dudal:2008xd . Remarkably, this phenomenon was observed by recent lattice numerical simulationsn Cucchieri:2009zt ; Dudal:2012hb . It is worth mentioning that, considering the Gribov-Zwanziger action as an effective action with an energy scale ultraviolet cutoff, it is possible to show that, at the strong coupling, the refinement is also favored in $d>4$, Guimaraes:2016okb .

For completeness, we display the form of the tree-level gluon propagator in333Due to the different values of $d$, one should keep in mind the different meanings of the space-time indices and mass dimensions. $d=3,4$

[TABLE]

and in $d=2$,

[TABLE]

II.2 Going beyond the Landau gauge and the non-perturbative BRST symmetry: linear covariant gauges

Although not peculiar to the Landau gauge, Gribov copies are very difficult to be handled when one chooses a different gauge condition. The main reason is that in the Landau gauge, the transversality of the gauge field ensures that the Faddeev-Popov operator $\mathcal{M}^{ab}$ is Hermitean. As such, this operator has a real spectrum which meaningfully allows for a definition of a Gribov region $\Omega$, where it is positive. Nevertheless, in general, the Faddeev-Popov operator is not Hermitean. The lack of such a property hinders a direct and clear definition of what would be the Gribov region in gauges different from the Landau gauge. There are two notable examples of gauge choices which also possess a Hermitean Faddeev-Popov operator: the maximal Abelian and Coulomb gauges. For such gauges, an explicit construction of the Gribov-Zwanziger action and its refinement was performed, see Dudal:2006ib ; Capri:2006cz ; Capri:2008ak ; Capri:2008vk ; Capri:2010an ; Gongyo:2013rua ; Capri:2015pfa ; Grotowski:1999ay ; Zwanziger:2006sc ; Burgio:2008jr ; Reinhardt:2008pr ; Golterman:2012dx ; Guimaraes:2015bra ; Burgio:2016nad . Despite of this fact, these gauges have their own peculiarities and the development of the Refined Gribov-Zwanziger scenario for them is not at the same level as in the Landau gauge.

A natural extension of the Landau gauge, which preserves Lorentz and color covariance, is given by the so-called linear covariant gauges, whose corresponding gauge condition is written as

[TABLE]

with $\alpha$ a non-negative gauge parameter and $b^{a}$ being, at this level, a given function. Clearly, if one sets $\alpha=0$, the Landau gauge is recovered. Infinitesimal Gribov copies in these gauges are characterized by the zero-modes equation,

[TABLE]

where, in contrast to Landau gauge, $\partial_{\mu}A^{a}_{\mu}\neq 0$ in general. It is precisely the fact that the gauge field is not purely transverse in these gauges that spoils the Hermiticity of $\mathcal{M}^{ab}_{\mathrm{LCG}}$.

The lack of Hermiticity makes the definition of the analogue of the Gribov region in linear covariant gauges very difficult. The first strategy to circumvent this technical difficulty was to take $\alpha$ as an infinitesimal parameter Sobreiro:2005vn , namely, the linear covariant gauge is taken as a small perturbation of the Landau gauge. As a consequence, it was proven in Sobreiro:2005vn that, in this situation, one can restrict the transverse component of the gauge field $A^{T,a}_{\mu}=(\delta_{\mu\nu}-\partial_{\mu}\partial_{\nu}/\partial^{2})A^{a}_{\nu}$ to the Gribov region $\Omega$ in the domain of integration in the path integral and all infinitesimal Gribov copies are removed. In Capri:2015pja , it was pointed out that the same strategy works for finite values of $\alpha$ with the exception of pathological infinitesimal Gribov copies, corresponding to zero modes which are not smooth functions of the gauge parameter $\alpha$. Hence, modulo a certain subclass of pathological copies, the restriction of the domain of integration in the path integral to the region where the transverse component belongs to $\Omega$ removes the infinitesimal Gribov copies. The resulting action was expressed in local form Capri:2015pja and its renormalizability proof to all orders in perturbation theory was achieved in Capri:2016aif . We also refer to Moshin:2015gsa .

Nevertheless, the presence of the gauge parameter $\alpha$ allows for an explicit check of the gauge independence of correlation functions of gauge invariant operators. In standard perturbation theory, this is controlled by the BRST symmetry. However, the soft breaking of the BRST symmetry in the (refined) Gribov-Zwanziger setup gives rise to non-trivial complications for such a task. Nevertheless, recently, in Capri:2015ixa , a reformulation of the Gribov-Zwanziger action in the Landau gauge in terms of a transverse and gauge invariant field444We refer to Appendix A of Capri:2015ixa for the construction of the gauge invariant field $A^{h,a}_{\mu}$., see Zwanziger:1990tn ; Lavelle:1995ty ; Lavelle:2011yc , $A^{h,a}_{\mu}$, with $\partial_{\mu}A^{h,a}_{\mu}=0$, and its generalization to linear covariant gauges was proposed. In this new formulation, the Gribov-Zwanziger enjoys an exact nilpotent BRST symmetry, which is a direct consequence of the gauge invariance of $A^{h,a}_{\mu}$ and which enables us to establish the independence from the parameter $\alpha$ of the gauge invariant correlation functions, and this even in the presence of the Gribov horizon.

As shown in Appendix A of Capri:2015ixa , the gauge invariant field555We write it in the matrix notation $A^{h}_{\mu}=A^{h,a}_{\mu}T^{a}$, with $T^{a}$ the generators of $SU(N)$. $A^{h}_{\mu}$ is expressed as an infinite series in powers of $A_{\mu}$, namely

[TABLE]

which, albeit transverse and gauge invariant, is a non-local expression. Upon a suitable redefinition of the field $b^{a}$, $b^{a}\rightarrow b^{h,a}$ Capri:2015ixa , with the introduction of the gauge invariant field $A^{h,a}_{\mu}$, the resulting Gribov-Zwanziger action in linear covariant gauges is written as Capri:2015ixa

[TABLE]

with

[TABLE]

and

[TABLE]

Before proceeding, one should note that the horizon function $H(A^{h})$ has now two sources of non-localities: the first one is related to the inverse of the operator $\mathcal{M}(A^{h})$, which is similar to the non-locality of the horizon function in the Landau gauge, see eq.(7). The second source of non-locality is associated with the field $A^{h}_{\mu}$ itself, see eq.(20). In order to localize the first type of non-locality present in (22), one proceeds as in the Landau gauge and introduces the set of auxiliary fields $(\bar{\varphi},\varphi,\bar{\omega},\omega)^{ab}_{\mu}$, which gives rise to the following action666We omit the vacuum term $-dV\gamma^{4}(N^{2}-1)$.

[TABLE]

Some properties of (24) are listed: i) The action $S^{\mathrm{LCG}}_{\mathrm{GZ}}$ is non-local due to the presence of the field $A^{h}_{\mu}$; ii) In the limit $\alpha\rightarrow 0$, i.e. $\partial_{\mu}A^{a}_{\mu}=0$, one has that $A^{h}_{\mu}\rightarrow A^{T}_{\mu}$ and the action (25) is equivalent to (9), the Gribov-Zwanziger action in the Landau gauge; iii) The action (24) enjoys an exact nilpotent BRST symmetry defined by the following transformations,

[TABLE]

with

[TABLE]

Up to now, we have presented a BRST invariant non-local action, eq.(24), which restricts the domain of integration in the path integral to a region free from a large set of Gribov copies. Moreover, as reported in Capri:2016aqq , this action can be fully localized by means of the introduction of additional auxiliary fields. In particular, the localization procedure worked out in Capri:2016aqq relies on the introduction of an auxiliary Stueckelberg-type field $\xi^{a}$, namely

[TABLE]

The field $A^{h}_{\mu}=A^{h,a}_{\mu}T^{a}$ is expressed in terms of the local field $\xi^{a}$ as

[TABLE]

An important feature of $A^{h}$, as defined by (28), is that it is gauge invariant, that is

[TABLE]

as can be explicitly seen through a gauge transformation parametrized by the $SU(N)$ matrix $V$

[TABLE]

Although non-polynomial, the field $A^{h}_{\mu}$ (28) is now a local field and can be expanded in terms of $\xi^{a}$, yielding

[TABLE]

Also, we must impose that the local field $A^{h}_{\mu}$, eq.(28), is transverse, namely, $\partial_{\mu}A^{h}_{\mu}=0$. Solving the transversality condition for the local field $\xi^{a}$ field, we obtain back the non-local expression for $A^{h}_{\mu}$ of eq.(20), see Appendix A of Capri:2015ixa . Therefore, besides introducing the field $\xi^{a}$, we should enforce the transversality of $A^{h}_{\mu}$ by means of a Lagrange multiplier $\tau^{a}$, a task which can be accomplished by introducing in the action the term

[TABLE]

We are now ready to write down the local and non-perturbative BRST invariant Gribov-Zwanziger action in the linear covariant gauges, i.e.

[TABLE]

The local action turns out to be renormalizable to all orders in perturbation theory Capri:2017 , while implementing the restriction of the domain of integration in the path integral to a region free from a large set of Gribov copies in the linear covariant gauges in a BRST-invariant way. Such a feature allows for a well-defined Slavnov-Taylor identity, through which the gauge parameter independence of gauge-invariant correlation functions can be established. An extensive analysis of these properties was carried out in Capri:2016aqq and Capri:2016gut .

As discussed in Subsect. II.1, the action (33) needs to be further refined, due to the dynamical formation of dimension two condensates. This fact was exploited in Capri:2015nzw where it was verified that, as in the Landau gauge, the refinement of the Gribov-Zwanziger action occurs in $d=3,4$, while in $d=2$ it is forbidden due to the presence of infrared singularities which prevent the formation of the dimension two condensates. Hence, in $d=3,4$, the action (33) is replaced by its refined version

[TABLE]

The tree-level gluon propagator computed out of (34) is given by

[TABLE]

being in very good agreement with the most recent lattice data Cucchieri:2009kk ; Cucchieri:2011pp ; Bicudo:2015rma . Although the transverse part of the propagator might acquire loop corrections, the longitudinal sector is exact to all orders, a consequence of the BRST symmetry. It is worth mentioning that this propagator has a decoupling/massive behavior in $d=3,4$, while in $d=2$ it is of scaling type due to the absence of refinement, see Capri:2015nzw . For completeness, the local Refined Gribov-Zwanziger action, eq.(34), is invariant under the nilpotent BRST transformations

[TABLE]

from which the BRST transformation of the field $\xi^{a}$, eq.(27), can be evaluated iteratively, giving

[TABLE]

It is instructive to check here explicitly the BRST invariance of $A^{h}$. For this, it is better to employ a matrix notation for the fields, namely

[TABLE]

with $A_{\mu}=A^{a}_{\mu}T^{a}$, $c=c^{a}T^{a}$, $\xi=\xi^{a}T^{a}$. From expression (28) we get

[TABLE]

Finally, we have

[TABLE]

It is important to emphasize that, in the action (34), the massive parameters $(\gamma,m,M)$ are coupled to BRST invariant expressions which are easily verified to be not BRST exact, i.e. cannot be expressed as pure $s$-variations. This fact ensures that these parameters are not akin to gauge parameters, having a physical meaning. As such, they will be present in the gauge-invariant correlation functions. Also, they are not free, being determined by their own gap equations as discussed in Dudal:2008sp ; Dudal:2011gd .

II.3 Curci-Ferrari gauge

In Pereira:2016fpn , it was argued that the Gribov problem in the Curci-Ferrari gauge is intimately related to the existence of copies in the linear covariant gauges. By a suitable shift of the $b$-field, it was shown that the copies equation is the same in both gauges. As such, the issue of the Gribov copies can be handled in the same way and the implementation of the restriction of the domain of integration in the path integral is obtained by the introduction of the same horizon function (22). As discussed in Pereira:2016fpn , the Gribov-Zwanziger action in the Curci-Ferrari gauge is

[TABLE]

As in the case of the linear covariant gauges, this theory suffers from non-perturbative instabilities which give rise to the dynamical formation of condensates in $d=3,4$. Therefore, expression (41) is refined by the inclusion of the same operators as in eq.(34) i.e.

[TABLE]

This action is invariant under the BRST transformations of eq.(36). The resulting tree-level gluon propagator coincides with that given in expression (35). Nevertheless, since the Curci-Ferrari gauge is non-linear777The non-linearity of the Curci-Ferrari gauge can be appreciated through the fact that, upon elimination of the Lagrange multiplier field $b^{a}$, a quartic ghost interaction term shows up., it does not enjoy the same set of Ward identities as the linear covariant gauges. A particular consequence of this fact is that, unlike the case of the linear covariant gauge, the longitudinal part of the propagator is now affected by quantum corrections.

II.4 Maximal Abelian gauge (MAG)

In order to construct the BRST-invariant (Refined) Gribov-Zwanziger action in the MAG, let us first set our conventions for this gauge. To avoid unnecessary complications, we restrict ourselves to the case of the gauge group $SU(2)$. In this case, the gauge field $A_{\mu}=A^{a}_{\mu}T^{a}$ can be decomposed into diagonal and off-diagonal components, as

[TABLE]

with $\alpha=\left\{1,2\right\}$ denoting the indices corresponding to the off-diagonal components. The diagonal generator $T^{3}\equiv T$ belongs to the Cartan subalgebra of $SU(2)$. Therefore, the following commutation relations hold:

[TABLE]

with $\epsilon^{\alpha\beta}=\epsilon^{\alpha\beta 3}$ being the totally antisymmetric symbol. The explicit decomposition of the field strength $F^{a}_{\mu\nu}$ yields

[TABLE]

with $\mathcal{D}^{\alpha\beta}_{\mu}$ being the covariant derivative defined with respect to the Abelian component $A_{\mu}=A_{\mu}^{3}$, namely

[TABLE]

By means of eq.(45), we can express the Yang-Mills action as

[TABLE]

which is left invariant under the following infinitesimal gauge transformations,

[TABLE]

The MAG is defined by the gauge conditions,

[TABLE]

giving rise to the following Faddeev-Popov operator $\mathcal{M}^{\alpha\beta}(A)$:

[TABLE]

The gauge fixed Yang-Mills action in the MAG is written as

[TABLE]

As discussed in Capri:2015pfa , an analogous of the Gribov region $\Omega$ of the Landau gauge can be introduced in the MAG. More precisely, the Gribov region $\Omega_{\mathrm{MAG}}$ for the MAG is defined by

[TABLE]

As in the case of the Landau gauge, the restriction of the domain of integration in the path integral to the region $\Omega_{\mathrm{MAG}}$ can be achieved in a BRST-invariant way by the introduction of the following horizon function

[TABLE]

where $\mathcal{M}^{\alpha\beta}(A^{h})$ means

[TABLE]

and $\mathcal{D}^{\alpha\beta}_{\mu}(A^{h})=\delta^{\alpha\beta}\partial_{\mu}-g\epsilon^{\alpha\beta}A^{h,3}_{\mu}$. The Gribov-Zwanziger action in the MAG is thus given by

[TABLE]

As before, expression (55) has two sort of non-localities encoded in the horizon function $H_{\mathrm{MAG}}(A^{h})$. In complete analogy with the procedure described in Subsect. II.2, it is possible to cast the action (55) in a local fashion. The resulting local action is expressed by

[TABLE]

This action is invariant under the following BRST tranformations,

[TABLE]

with

[TABLE]

As in the case of the gauges discussed before, the Gribov-Zwanziger action in the MAG also suffers from non-perturbative instabilities and dimension two condensates are dynamically generated in $d=3,4$, while, in $d=2$, their formation is invalidated by infrared singularities. Therefore, as in the case of the previous gauges, the Gribov-Zwanziger action in the MAG does not refine in $d=2$. We refer to Capri:2015pfa for a detailed discussion of this feature. The Refined Gribov-Zwanziger action in $d=3,4$ in the MAG is written as

[TABLE]

where the mass parameters $(m^{2}_{\mathrm{diag}},m^{2}_{\mathrm{off}},M^{2})$ reflect the existence of the dimension-two condensates $\langle A^{h,3}_{\mu}A^{h,3}_{\mu}\rangle$, $\langle A^{h,\alpha}_{\mu}A^{h,\alpha}_{\mu}\rangle$, $\langle\bar{\varphi}^{\alpha\beta}_{\mu}\varphi^{\alpha\beta}_{\mu}-\bar{\omega}^{\alpha\beta}_{\mu}\omega^{\alpha\beta}_{\mu}\rangle$.

The diagonal gluon propagator is given by

[TABLE]

while the off-diagonal gluon propagator is

[TABLE]

From expression (61), we see that the off-diagonal gluon propagator displays a Yukawa type behavior. Lattice simulations give support to this result, see Bornyakov:2003ee ; Mendes:2008ux ; Gongyo:2012jb ; Gongyo:2013sha . Moreover, this behavior is in agreement with the Abelian dominance scenario Ezawa:1982bf , where off-diagonal gluons should acquire a dynamical mass, reponsible for their decoupling at low energy. On the other hand, the diagonal gluon propagator (60) is of the refined Gribov type. As such, it is infrared suppressed and attains a non-vanishing value for $k=0$, in agreement with the lattice studies Bornyakov:2003ee ; Mendes:2008ux ; Gongyo:2012jb ; Gongyo:2013sha . The diagonal gluon propagator also displays reflection positivity violation, a feature which is interpreted as a signal of confinement. Again, this result is in agreement with the Abelian dominance scenario.

III Non-perturbative coupling of scalar fields in the adjoint representation

In this section, we generalize the construction of Capri:2014bsa to linear covariant, Curci-Ferrari and maximal Abelian gauges. To begin with, we consider scalar fields in the adjoint representation of888In the case of the MAG, we restrict ourselves to $SU(2)$ for simplicity. $SU(N)$. The idea proposed in Capri:2014bsa consists in the introduction of a term akin to the horizon function for the matter sector, which provides a non-perturbative coupling between matter fields and the gauge sector. Although for the gluon sector the horizon function has a clear geometrical meaning, implementing the restriction of the domain of integration in the path integral to the Gribov region, the introduction of an analogous term in the matter sector does not yet exhibit the same well defined geometric support. Nevertheless, recently, it was observed that such a non-perturbative coupling between matter and gauge fields could be motivated through the dimensional reduction of higher-dimensional Yang-Mills theory, see Guimaraes:2016okb . More precisely, upon reduction of a five dimensional Yang-Mills to the four dimensional theory Guimaraes:2016okb , a non-perturbative coupling between the scalar field corresponding to the fifth component of the gauge connection and the four dimensional gauge field shows up, being precisely of the type introduced in Capri:2014bsa . As we shall see, this prescription gives rise to non-perturbative matter fields propagators which turn out to be in good agreement with lattice data, whenever available.

III.1 Linear covariant and Curci-Ferrari gauges

Let us consider the standard action of scalar fields in the adjoint representation of $SU(N)$, minimally coupled with the gauge sector, i.e.

[TABLE]

Of course, expression (62) is left invariant by BRST transformations (25), with the scalar field $\phi^{a}$ transforming as

[TABLE]

where $\{T^{a}\}$ stand for the generators of $SU(N)$ in the adjoint representation.

Making use of the Stueckelberg field $\xi$, eq.(27), a BRST invariant scalar field is constructed as follows Capri:2016aqq :

[TABLE]

To first order, we get

[TABLE]

It is easy to verify that $\phi^{h}$ is left invariant by the BRST transformations, i.e.

[TABLE]

The prescription introduced in Capri:2014bsa amounts to introduce the following non-local BRST invariant term to the scalar action (62),

[TABLE]

where $\left(\mathcal{M}(A^{h})\right)^{ad}$ stands for the Faddeev-Popov operator of eq.(23). It is almost immediate to realise that expression (67) shares great similarity with the horizon function of the gluon sector, eq.(22). In fact, as already mentioned, expression (67) can be obtained through the dimensional reduction of higher-dimensional Yang-Mills theory Guimaraes:2016okb .

The action of the scalar field with the addition of the non-perturbative coupling (67) is given by

[TABLE]

where the massive parameter $\sigma$ plays the same role of the Gribov parameter $\gamma$. Again, due to the presence of the operator $\mathcal{M}^{-1}$ in expression (67), the action (68) is non-local. Moreover, it turns out to be possible to cast the action $S^{\phi}$ in local form following the same procedure adopted in the previous sections for the localization of the Gribov-Zwanziger action. To that purpose, we introduce a set of auxiliary fields $(\bar{\eta},\eta,\bar{\theta},\theta)^{ab}$ akin to Zwanziger’s localizing fields in such a way that

[TABLE]

The fields $(\bar{\eta},\eta)$ are commuting while $(\bar{\theta},\theta)$ are anti-commuting. Integrating out these fields in the functional integration gives back the non-local expression (67).

Therefore, the local scalar field action non-perturbatively coupled to the gauge sector is expressed by

[TABLE]

One should keep in mind that in expression (70), both $A^{h}_{\mu}$ and $\phi^{h}$ are expressed in terms of the Stueckelberg field $\xi^{a}$ and are thus local fields, albeit non-polynomial.

As it happens in the gauge sector of the Gribov-Zwanziger action, the non-local mass term (67) entails non-perturbative instabilities which give rise to the dimension two condensates, $\langle\phi^{h,a}\phi^{h,a}\rangle$ and $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$, akin to those of the Refined Gribov-Zwanziger action. It is worth to proceed by evaluating those condensates at first order, a task which can be accomplished by introducing the operators

[TABLE]

in expression (70), where $(J,\tilde{J})$ are constant sources. Thus, we define the action $\Sigma(J,\tilde{J})$ by

[TABLE]

To first order, the condensates $\langle\phi^{h,a}\phi^{h,a}\rangle$ and $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$ can be obtained by taking the derivatives of the one-loop vacuum energy $\mathcal{E}^{(1)}$ with respect to the sources $(J,\tilde{J})$, and setting them to zero, where

[TABLE]

$\Sigma^{(2)}(J,\tilde{J})$ denotes the quadratic part of (72), while the path integral measure is expressed as

[TABLE]

At one-loop order, the vacuum energy is easily evaluated, being given by

[TABLE]

where dimensional regularization has been employed. Therefore, at first order, for the condensates we get

[TABLE]

One sees that the one-loop result already shows non-vanishing expressions for the condensates $\langle\phi^{h,a}\phi^{h,a}\rangle$ and $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$. Remarkably, the contributions coming from the introduction of the non-perturbative mass term (67) to the standard scalar field action are ultraviolet convergent. Interesting to note, very much alike the refinement of the Gribov-Zwanziger action, infrared singularities show up in the integrals (76), preventing the formation of such condensates in $d=2$. As in the case of the gluon sector, in $d=3,4$, the effects of the existence of the condensates $\langle\phi^{h,a}\phi^{h,a}\rangle$ and $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$ can be taken into account by refining the matter action as:

[TABLE]

where the parameters $(\tilde{m}^{2}_{\phi},\rho^{2})$ have dynamical origin and can be obtained through the evaluation of the effective potential for $\langle\phi^{h,a}\phi^{h,a}\rangle$ and $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$. Furthermore, in the case of $d=2$, due to the absence of condensates, the action remains the original one given by eq.(70).

After these considerations, we can compute the tree-level scalar field propagator for different values of $d$. In $d=2$, we have

[TABLE]

while in $d=3,4$,

[TABLE]

In analogy with the case of gluon propagator, the scalar field propagator attains a finite value at zero momentum in $d=3,4$ while in $d=2$ it vanishes at $k=0$. In both cases the scalar propagator is infrared suppressed. Also, at the tree-level, there is no $\alpha$-dependence as it is apparent from eqs.(78) and (79). Hence, the Landau limit $\alpha=0$ is trivial and agrees with the results reported in Dudal:2008sp . Also, the propagators (78) and (79) violate reflection positivity, a feature which is interpreted as a signal of confinement. We see thus that the introduction of the non-perturbative matter coupling (67) has the effect of confinining the scalar matter fields.

It is worth here to add some further remarks on the specific case of $d=2$, eq.(78). One should keep in mind that expression (78) is a consequence of the first-order absence of the condensate $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$, as it follows from eq.(76). Though, we underline that this is only a first order analysis. As such, expression (78) would retain its validity at this order. Willing to make an all order statement, a higher loop analysis of the condensate $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$ would be required, a matter which is well beyond the aim of the present paper. Although the available lattice simulations Maas:2011yx point towards a similar behavior for the scalar field propagator in the infrared for different values of $d=4,3,2$ in the Landau gauge, in order to make a comparison with the lattice data in $d=2$ a detailed analysis of the higher order condensate $\langle\bar{\eta}^{ab}\eta^{ab}-\bar{\theta}^{ab}\theta^{ab}\rangle$ would definitively be needed.

Finally, as dicussed in Subsect. II.3, the Gribov problem in the Curci-Ferrari and linear covariant gauges can be treated by means of a formal equivalence. As such, the non-perturbative matter term, eq.(67), in both gauges is the same. Therefore, the scalar field action non-perturbatively coupled to the gauge sector is given by (70). Clearly, at first order, all the computations presented in this section remain valid for the Curci-Ferrari gauges, namely, the calculation of the vacuum energy and of the scalar field propagator. Of course, taking into account higher loops contributions, the non-linear character of the Curci-Ferrari gauges will show up giving results which will differ from those of the linear covariant gauges. Since this is beyond the scope of the present work, we limit ourselves to the first order computations already presented in the case of linear covariant gauges, which retain their validity also in the Curci-Ferrari gauges.

III.2 Maximal Abelian gauge

In the case of the MAG, although the prescription is the same, care is due to the decomposition of color indices into diagonal and off-diagonal ones. Firstly, we express the minimally coupled scalar field action in a color decomposed fashion, namely

[TABLE]

with $\phi\equiv\phi^{3}$.

As in the case of linear covariant and Curci-Ferrari gauges, the non-perturbative matter coupling is obtained through the addition, in the scalar field action, of a non-local term which shares great similarity with the corresponding horizon function of the gluon sector in the MAG, eq.(53), i.e.

[TABLE]

where the Faddeev-Popov operator $\mathcal{M}(A^{h})$ is now given by eq.(54). The scalar field action supplemented with the non-perturbative coupling (81) becomes

[TABLE]

The parameter $\sigma$ has mass dimension and is the analogue of the Gribov parameter $\gamma$ in the matter sector. As before, the non-local action (82) can be cast in local form by means of the introduction of auxiliary fields and of a Stueckelberg field, also used to localize $A^{h}_{\mu}$. In local form, the action (82) is written as

[TABLE]

with $\phi^{h}=h^{\dagger}\phi h$.

As pointed out in Subsect. III.1, the auxiliary localizing fields $(\bar{\eta},\eta,\bar{\theta},\theta)^{\alpha\beta}$ develop their own dynamics and give rise to the dynamical formation of condensates. This is in very much analogy with the refinement of the Gribov-Zwanziger action. In order to explicitly check the existence of such condensates to first order, we proceed as before and introduce the following operators to (83),

[TABLE]

where $J$ and $\tilde{J}$ are constant sources. This gives rise to

[TABLE]

Our aim is to compute the following condensates at one-loop order:

[TABLE]

This is achieved by taking the derivatives with respect to $J$ and $\tilde{J}$ of the vacuum energy $\mathcal{E}$, defined by

[TABLE]

and setting the sources to zero at the end, namely

[TABLE]

The measure of the path integral (87) is written as

[TABLE]

At one-loop order, we should take the quadratic part of999For the moment, we can ignore the contribution from $S^{\mathrm{loc}}_{\mathrm{RGZ}}$, which is $(J,\tilde{J})$-independent. $\Sigma(J,\tilde{J})$,

[TABLE]

Integrating the auxiliary fields $(\tau^{\alpha},\tau)$ which enforce the transversality condition of $A^{h,a}=(A^{h,\alpha}_{\mu},A^{h}_{\mu})$, we see that the gauge-invariant scalar field $\phi^{h,a}=(\phi^{h,\alpha},\phi^{h,3})$ can be expressed as in eq.(65) with $\xi^{a}=\frac{\partial A^{a}}{\partial^{2}}$.

Since we want to maintain the action $\Sigma^{(2)}$ to the quadratic order in the fields, we see that $\phi^{h,a}\approx\phi^{a}$. Hence, the $(J,\tilde{J})$-dependent part of the one-loop order vacuum energy $\mathcal{E}$ is given by

[TABLE]

This implies,

[TABLE]

and

[TABLE]

Eq.(93) shows that, at the one-loop level, the condensate of auxiliary fields, $\langle\bar{\eta}^{\alpha\beta}\eta^{\alpha\beta}-\bar{\theta}^{\alpha\beta}\theta^{\alpha\beta}\rangle_{\mathrm{1-loop}}$, is ultraviolet convergent. For $d=3,4$, such a condensate is perfectly well-defined in the infrared region and can be safely introduced. In $d=2$, an infrared singularity at $k=0$ turns out to appear. This is in agreement with the refining condensates in the Gribov-Zwanziger setup. From eq.(92), we see that the condensate $\langle\phi^{h,3}\phi^{h,3}\rangle$ has two contributions: one proportional to $\tilde{m}^{2}_{\mathrm{diag}}$ which exists irrespective of the presence of $\sigma$ and the other one proportional to $\sigma^{4}$. The former contains an ultraviolet divergence which can be taken into account by the standard renormalization techniques while the latter is ultraviolet convergent and free from infrared divergences in $d=3,4$. In $d=2$, an infrared singularity appears preventing the introduction of this condensate. We must emphasize that this condensate does not affect the qualitative behavior of the initial theory.

Therefore, in $d=2$, the scalar field action non-perturbatively coupled with the gauge sector is given by (83), while in $d=3,4$ the condensates $\langle\bar{\eta}^{\alpha\beta}\eta^{\alpha\beta}-\bar{\theta}^{\alpha\beta}\theta^{\alpha\beta}\rangle$ and $\langle\phi^{h,3}\phi^{h,3}\rangle$ have to be taken into account, giving rise to the following refined action

[TABLE]

From the actions (94) and (83) we can compute the tree-level Abelian component of the scalar field propagator. The expressions in $d=2$ and $d=3,4$ are, respectively,

[TABLE]

and

[TABLE]

From eq.(95) and (96), we see that the propagator of the Abelian component of the scalar field displays the same features observed for the tree-level propagator of the scalar field in the linear covariant and Curci-Ferrari gauges, eqs.(78),(78).

IV Generalization of the non-perturbative matter coupling for quark fields

In the previous section, we have presented a prescription for the non-perturbative coupling of scalar fields in the adjoint representation of the gauge group with the gauge sector. Such a coupling arises from the introduction of a non-local term which shares great similarity with the corresponding horizon term introduced in the gluon sector to implement the restriction of the domain of integration to the Gribov region. Interestingly, this term naturally appears through the dimensional reduction of higher-dimensional Yang-Mills theory Guimaraes:2016okb .

In the present section we follow the same reasoning for the case of fermionic matter fields in the fundamental representation of the gauge group. This case is particularly important since it allows us to obtain an analytic non-perturbative expression of the quark field propagator. As before, we divide the analysis in two subsections for linear covariant/Curci-Ferrari gauges and for the maximal Abelian gauge. We have collected our conventions regarding spinors and related issues in Appendix A.

IV.1 Linear covariant and Curci-Ferrari gauges

Let us begin by considering the Dirac action in Euclidean space minimally coupled with the gauge sector,

[TABLE]

where capital latin indices $\left\{I,J,\ldots\right\}$ stand for the fundamental representation of $SU(N)$. The covariant derivative $D^{IJ}_{\mu}$ is defined by

[TABLE]

with $T^{a}$ the generators of $SU(N)$ in the fundamental representation. In strict analogy to what has been proposed in Sect. III, the non-perturbative fermion matter coupling is introduced by adding to the Dirac action the non-local term

[TABLE]

with $\left(\mathcal{M}(A^{h})\right)^{ad}$ given by eq.(23) and where the gauge-invariant spinor $\psi^{h}$ is defined as

[TABLE]

Employing the Stueckelberg field $\xi^{a}$, the all order BRST invariant spinor field $\psi^{h}$ is obtained as

[TABLE]

From

[TABLE]

it immediately follows that $\psi^{h}$ is BRST invariant, namely

[TABLE]

Solving the transversality condition $\partial_{\mu}A^{h,a}_{\mu}=0$ for the Stueckelberg field $\xi^{a}$ and plugging it in eq.(101), see Appendix A of Capri:2015ixa , we reobtain expression (100). Hence, following the prescription discussed in Capri:2014bsa ; Capri:2016aqq , the fermionic action non-perturbatively coupled to the gauge sector is given by

[TABLE]

where $M$ is the analogue of the Gribov parameter $\gamma$ for the fermionic sector.

The term $\mathcal{H}(\psi^{h})$ is non-local due to the inverse of $\mathcal{M}(A^{h})$, eq.(23). Nevertheless, the action (104) can be localized in complete analogy with the localization of the Gribov-Zwanziger action by means of the introduction of commuting spinor fields $(\bar{\theta},\theta)^{aI}$ as well as of anti-commuting ones $(\bar{\lambda},\lambda)^{aI}$. The local form of expression (99) is given by

[TABLE]

which, upon integration over the auxiliary fields $(\bar{\theta},\theta)^{aI}$ and $(\bar{\lambda},\lambda)^{aI}$, gives back the non-local quantity of eq.(99)

Therefore, the local action with the non-perturbative coupling between fermionic matter and the gauge sector is expressed as

[TABLE]

As extensively discussed in the present work, the presence of the parameter $M$, akin to the Gribov parameter $\gamma$, and of the quadratic coupling between the auxiliary localizing fields and the corresponding matter field give rise to a dynamical and non-perturbative instability, resulting in the formation of condensates. Again, we present the one-loop computation which hints the existence of such condensates. To do so, we introduce the following operators

[TABLE]

into the action (106), yielding

[TABLE]

We aim at computing the following condensates:

[TABLE]

which can be obtained by taking taking the derivatives with respect to $(J,\tilde{J})$ of the vacuum energy101010We restrict ourselves to the contributions relevant for our purposes. $\mathcal{E}(J,\tilde{J})$ at one-loop order,

[TABLE]

with $\Sigma^{(2)}(J,\tilde{J})$ the quadratic part of $\Sigma(J,\tilde{J})$, namely

[TABLE]

Explicitly, $\Sigma^{(2)}(J,\tilde{J})$ is written as

[TABLE]

Performing the path integral over the auxiliary localizing fields yields the following expression

[TABLE]

Making use of the relation

[TABLE]

and performing the path integral over $(\bar{\psi},\psi)$, one obtains

[TABLE]

After simple manipulations and employing the identity

[TABLE]

one ends up with

[TABLE]

From (117) it is immediate to extract the vacuum energy $\mathcal{E}^{(1)}$, which is written as

[TABLE]

Finally, we are ready to compute the expectation values (111), by differentiating (118) with respect to the sources $(J,\tilde{J})$ as in (111). One obtains,

[TABLE]

and

[TABLE]

where the prescriptions of the dimensional regularization were employed. In $d=4$, we see from eq.(119) that the contribution which is directly proportional to the parameter $M$ is perfectly ultraviolet convergent. This is in agreement with the fact that the introduction of the non-local term of the type of eq.(99) does not introduce any new ultraviolet divergence Capri:2015mna . From eq.(120), we easily see that the one-loop contribution to the condensate $\langle\bar{\theta}^{aI}\theta^{aI}-\bar{\lambda}^{aI}\lambda^{aI}\rangle$ is non-vanishing and ultraviolet convergent. These results show explicitly, already at one-loop order, that the introduction of the non-perturbative matter coupling (120) contributes definitively to the formation of such condensates.

As usual, the dynamical formation of those condensates can be taken into account from the beginning by refining the matter sector in the following way

[TABLE]

Finally, one can compute the quark field propagator at tree level from the refined action (121). The result is

[TABLE]

with $M_{\psi}=m_{\psi}+\tilde{m}_{\psi}$.

The propagator (122) is the same as the one computed in the Landau gauge Capri:2014bsa , i.e. $\alpha=0$. Of course, higher orders correction will, eventually, introduce some $\alpha$-dependence in (122). In the particular case of $\alpha=0$, the propagator (122) fits well recent lattice data, see Capri:2014bsa and references therein. To the best of our knowledge, there are no available numerical simulations of the quark propagator in linear covariant gauges. Hence, our result could be a motivation for such an endeavour in the near future.

As described in the case of scalar fields, the generalization of the present construction to the case of the Curci-Ferrari gauges is straightforward. In particular, the results obtained here also hold in the Curci-Ferrari gauge, which differs from the linear covariant gauges by non-linear terms which do not contribute to the order we are dealing with. In particular, the quark propagator at the tree-level remains the same as in eq.(122).

IV.2 Maximal Abelian gauge

In this subsection, we proceed with the analysis of the non-perturbative coupling of quark matter fields in the maximal Abelian gauge case. In full analogy with the case of the scalar matter field, eq.(81), for the non-perturbative BRST invariant coupling in the quark sector we write

[TABLE]

where the gauge-invariant field $\psi^{h}$ is defined by eq.(101), while the Faddeev-Popov operator in the maximal Abelian gauge, $\mathcal{M}^{\alpha\beta}$, is given by eq.(50). The non-perturbative coupling of quark matter fields with the gauge sector in the maximal Abelian gauge is thus given by

[TABLE]

where, as before, the parameter $M$ plays an analogue role of the Gribov parameter $\gamma$ in the matter sector. As exhaustively discussed in the previous sections, the non-local quark matter term (123) can be localized by means of auxiliary fields. The gauge-invariant field $\psi^{h}$ can be written in local form in the same manner described in Subsect. IV.1. On the other hand, a pair of commuting $(\bar{\theta},\theta)^{\alpha I}$ and anticommuting $(\bar{\lambda},\lambda)^{\alpha I}$ fields are introduced in order to localize $\mathcal{H}_{\mathrm{MAG}}(\psi^{h})$, namely,

[TABLE]

Therefore, the action of quark matter fields coupled with the gauge sector in a non-perturbative way is expressed, in local form, as

[TABLE]

At this stage, it is not unexpected to predict that, again, the action (126) suffers from dynamical non-perturbative instabilities, giving rise to the formation of condensates. The procedure to explicit check the existence of such codensates goes exactly along the same lines of the previous case, namely: constant sources $J$ and $\tilde{J}$ are coupled to the composite operators $\bar{\psi}^{h,I}\psi^{h,I}$ and $(\bar{\theta}^{\alpha I}\theta^{\alpha I}-\bar{\lambda}^{\alpha I}\lambda^{\alpha I})$, i.e.

[TABLE]

which are introduced in the action (126), giving rise to

[TABLE]

The condensates are obtained by taking derivatives of the vacuum energy $\mathcal{E}$ corresponding to the action (128) with respect to the sources $J$ and $\tilde{J}$, and setting them to zero, i.e.

[TABLE]

with

[TABLE]

At one-loop order, using the same techniques presented in Sect. III and Subsect. IV.1, one obtains

[TABLE]

Plugging eq.(131) into eq.(129), one immediately gets

[TABLE]

and

[TABLE]

Once again, one notices that the contributions proportional to $M$ are ultraviolet finite. As such, we find already at one-loop order that such condensates are non-vanishing, due to the introduction of the non-perturbative coupling (123) in the quark matter sector. We should emphasize that, unlike the case of the linear covariant and Curci-Ferrari gauges, the condensate of the auxiliary fields, $\langle\bar{\theta}^{\alpha I}(x)\theta^{\alpha I}(x)-\bar{\lambda}^{\alpha I}(x)\lambda^{\alpha I}(x)\rangle$, is purely diagonal. This is a direct consequence of the decomposition into diagonal and off-diagonal indices of the maximal Abelian gauge. Finally, as before, the dynamical generation of the condensates (129) can be taken into account by the refinement of the quark action, i.e.

[TABLE]

Out of the action (134), one can compute the tree-level quark propagator, which is given by

[TABLE]

Quite importantly, the tree-level quark propagator (135) is in qualitative agreement with the very recent lattice results reported in Schrock:2015pna . Such an agreement works as a highly non-trivial check of the non-perturbative matter coupling proposed here.

V Conclusions

In this work, we have extended the non-perturbative gauge-matter coupling proposed in Capri:2016aqq ; Capri:2014bsa to linear covariant, Curci-Ferrari and maximal Abelian gauges. In particular, we have investigated the coupling of scalar fields in the adjoint representation of the gauge group as well as of quark fields in the fundamental representation.

The non-perturbative nature of the proposal relies on the introduction of an additional term in which the matter fields are coupled to the inverse of the operator $\mathcal{M}(A^{h})$, whose existence is ensured by the restriction of the domain of integration in the functional integral to the Gribov region. As discussed in details throughout the paper, this additional term in the matter fields shares great similarity with the horizon function introduced in the pure gauge sector in order to implement the restriction to the Gribov region. Albeit non-local, the resulting action can be cast in local form by the introduction of auxiliary fields which, as in the case of the localizing Zwanziger fields of the pure gauge sector, develop their own dynamics giving rise to the formation of condensates, as explicitly checked through one-loop computations. Moreover, the condensates arising in the matter sector can be taken into account through an effective action which looks much alike the Refined Gribov-Zwanziger action which accounts for the existence of similar condensates in the gluon sector. Out of this action, the tree-level propagators for matter fields were analysed, giving rise to reflection positivity violating propagators. As in the case of the gluon propagator, the positivity violation is taken as a signal that colored matters fields are confined too.

We emphasize that the final effective action which encodes the non-perturbative effects of the matter sector is invariant under BRST transformations. This was achieved by the introduction of the suitable gauge-invariant fields $A^{h}$, $\phi^{h}$ and $\psi^{h}$, see Capri:2015ixa ; Capri:2016aqq ; Capri:2015nzw ; Pereira:2016fpn ; Capri:2016gut , which, albeit local, are non-polynomial in the auxiliary Stueckelberg type field $\xi^{a}$. Nevertheless, such variables as well as the proposed non-perturbative matter coupling give rise to a local ation which can be proven to be renormalizable to all orders, see Fiorentini:2016rwx ; Capri:2017 .

The present work can give rise to several future investigations among which we quote: i) as done in the pure gauge sector Capri:2016gut , we are now ready for a detailed analysis of the Nielsen identities, in the case of linear covariant gauges, to investigate the independence of the poles of the matter field propagator from the gauge parameter $\alpha$; ii) use the gauge-invariance of $A^{h}$, $\phi^{h}$ and $\psi^{h}$ to explore the Landau-Khalatnikov-Fradkin tranformations, as briefly discussed in Capri:2016gut for the gluon sector, and analyse how gauge-matter correlators depend on the gauge parameter $\alpha$, while checking out how the results compare with those obtained through the aforementioned Nielsen identitites; iii) study of how the presence of the Higgs mechanism can drive the transition between the confining and de-confining regimes in a BRST invariant fashion, iv) investigate how the present proposal generalizes to supersymmetric gauge theories, v) stimulate different groups from other approaches such as lattice simulations and Dyson-Schwinger equations to study two-point functions of matter fields away from Landau gauge. As in the gluon sector, the interplay between different approaches in the study of non-perturbative correlation functions will certainly be very successful.

Acknowledgements

The Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-Brazil) and The Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) are acknowledged.

Appendix A Conventions in Euclidean space

The gamma matrices $\gamma_{\mu}$ obey the Clifford algebra

[TABLE]

with

[TABLE]

and

[TABLE]