# Non-perturbative analysis of the Gribov-Zwanziger action

## Abstract

In the non-perturbative regime the usual gauge fixing is not sufficient due to the Gribov problem. To deal with it one can restrict the integration in the path integral to the first Gribov region by using the Gribov-Zwanziger action. In its local form it features additional auxiliary fields which mix with the gluon at the two-point level. We present an explicit infrared analysis of this action. We show that from the two possible scaling solutions obtained previously only one remains: It coincides exactly with the results from the Faddeev-Popov action, i.e., the ghost propagator is infrared enhanced and the gluon propagator infrared suppressed and the corresponding power law behavior is described by only one parameter . This corroborates the argument by Zwanziger that for functional equations it suffices to take into account the appropriate boundary conditions and no explicit restriction in the path integral measure is required.

###### :

11.10.-z,03.70.+k,11.15.Tk8x11double

address=Physikalisch-Theoretisches Institut, Friedrich-Schiller-Universität Jena, Max-Wien-Platz 1, 07743 Jena, Germany

address=Institut für Physik, Universität Graz, Universitätsplatz 5, 8010 Graz, Austria

address=Departamento de Física Teórica, Instituto de Física, UERJ - Universidade do Estado do Rio de Janeiro, Rua São Francisco Xavier 524, 20550-013 Maracanã, Rio de Janeiro, Brazil

Understanding confinement, i.e., the absence of quarks and gluons from the physical spectrum, has been a long-standing challenge for the theory of strong interations. In the course of time several seemingly different confinement scenarios have emerged. By now it has become clear that different explanations are not mutually exclusive but are only different manifestations. Especially in the last fifteen years functional methods have provided new insights into some of these scenarios, as the basic quantities in these approaches, the Green functions, give direct access to some of them, see, e.g., Alkofer and von Smekal (2001); Braun et al. (2010). Because of its accessibility mostly the Landau gauge has been investigated in this context, see, for example, von Smekal et al. (1998); Zwanziger (2002); Lerche and von Smekal (2002); Pawlowski et al. (2004); Fischer et al. (2009).

A particular gluon confinement picture is due to Gribov and Zwanziger Gribov (1978); Zwanziger (1989), where the gluon is removed from the physical spectrum due to a strong suppression of its propagator at low momenta. The argument is based on the insufficient gauge fixing implemented by the Faddeev-Popov procedure Gribov (1978). While this is not harmful for perturbation theory it is important for the non-perturbative regime. The so-called Gribov-Zwanziger action incorporates a refined restriction in field configuration space in an attempt to avoid this problem Gribov (1978); Zwanziger (1989): Instead of integrating over all configurations fulfilling the Landau gauge condition one takes into account only configurations from the first Gribov region defined as . As a result already the tree level gluon propagator is confined because it vanishes at zero momentum: In this case positivity, a requirement for any physical particle, is violated maximally Zwanziger (1991). The ghost propagator on the other hand is found to diverge more strongly than the tree-level part, namely like Gribov (1978); Zwanziger (1993). This result is only obtained when enforcing the horizon condition, which arises as an additional condition in the gauge fixing procedure. Also in our approach it plays a pivotal role.

Although the insufficiency of the standard gauge fixing has been long known, Dyson-Schwinger equations (DSEs) relied up to now on the Faddeev-Popov action. To understand why this is working we refer to a conjecture of Zwanziger Zwanziger (2002, 2004): Deriving DSEs with an integration in field configuration space formally restricted to the Gribov region does not change their form. The reason is that cutting the integration at the Gribov horizon does not introduce any boundary terms as the Faddeev-Popov determinant, which appears in the integrand of the path integral, vanishes on the horizon. Hence one can use the Faddeev-Popov action but has to take into account a proper boundary condition.

The original Gribov-Zwanziger action contains a non-local term Zwanziger (1989), which can be localized with auxiliary fields , , and Zwanziger (1989). The statistics obeyed by these fields is important: The first two are anti-commuting and the other two commuting. The additional adjoint color index accounts for the correct number of degrees of freedom. For our purposes the local action is rewritten as

(1) |

where the field is the imaginary part of and :

(2) |

The anti-commuting fields and comprise the Faddeev-Popov ghosts, the real part of and and the fields and . The index runs from . Treating all these fields with only one pair of fields is possible, because they all interact only via the Faddeev-Popov operator . The massive parameter is the Gribov parameter. Its value is fixed by the horizon condition. Although the part proportional to breaks BRST invariance, see, e.g., Zwanziger (1993); Capri et al. (2010), this action is renormalizable. See Dudal et al. (2010) and references therein for more details on the renormalization and the horizon condition.

DSEs for this action can be derived with the Mathematica package DoDSE Alkofer et al. (2009) which is of great help here, as due to the existence of mixed propagators even the two-point DSEs become quite complicated: The , , , and two-point function DSEs have , , , and terms. The equations can be found in Huber et al. (2010a).

For the analysis of the equations we write the full two-point function as

(3) |

For the two-point functions and we make ansätze which take only the tensors from the bare action into account:

(4) | ||||

(5) |

The functions , , are non-perturbative dressing functions for the gluon, field and mixed two-point functions. In principle more tensors are possible. In order to obtain the corresponding propagators we have to invert the two-point matrix , :

(6) |

(7) | |||

(8) |

The propagator of the field is unaffected by the mixing:

(9) |

For determining the behavior of Green functions at low momenta we study the DSEs in this limit. Then the integrand is dominated by small momenta and all dressing functions can be replaced by their infrared (IR) expressions. For the dressing functions we make the ansätze and similarly for the vertices. The complete qualitative behavior is now described by the IR exponents . For the propagators we cannot say which part will be dominant in the IR, so we have to distinguish four scenarios. For every case one can then perform an analysis along the lines of ref. Huber et al. (2010b), where a general method for the calculation of possible IR scaling solutions is described. It was shown in ref. Huber et al. (2010a) that two scenarios are inconsistent, as the IR analyses yield conflicting conditions for the IR exponents. Hence we concluded that only two possible solutions remain, which are, however, both in qualitative agreement with the results from the Faddeev-Popov action and the IR leading diagrams are given by those which have a bare gluon- or gluon- vertex. Here we will present evidence that one of those solutions can also be discarded.

This case was denoted III in Huber et al. (2010a) and is given by , i.e., both parts of the determinant appearing in the propagators, , scale equally. The scaling analysis resulted in , and , where is a constant to be determined. For this one has to solve a system of equations consisting of the four two-point DSEs truncated to their IR leading diagrams. The unknowns are the combinations of the dressing function constants , , and . As the determinant remains unchanged here, the variables appear in a more involved way than for the second case treated below, where the propagators have a simple dependence on the dressing functions. Consequently a complete analytic solution is not possible. However, one can solve for two variables analytically and try to find a solution numerically for the remaining equations. The required analytic expressions can be obtained with the current version of DoDSE Alkofer et al. (2009). is then expressed analytically from one of the equations as is or . The results are plugged into the remaining two equations. Finally, they are plotted for the allowed values of the variables: is known to be positive from the scaling analysis and it has to be smaller than in order for the Fourier transformations of the propagators to exists. The values of and must be positive as the dressings and are necessarily positive. It turns out that no solution exists. Hence this case is ruled out for the chosen ansatz.

Therefore the only solution in our truncation is given by , and and the dominant part in the determinant is . However, for the propagator the second term is suppressed compared to the first one, so that its final form in the IR is

(10) |

As a consequence the dependence of the IR leading diagrams in the two-point DSEs of the and field on the mixed propagator vanishes. Even more, the equations turn out to be completely the same in the IR as those obtained from the Faddeev-Popov action, because the contributions from the field and the fields cancel each other such that only contributions equivalent to those from the Faddeev-Popov ghosts remain. Hence also the result for , which can be calculated analytically, is the same: Zwanziger (2002); Lerche and von Smekal (2002). Furthermore, the IR exponents of all vertices are determined as Alkofer et al. (2005); Huber et al. (2008). With this information the IR exponent of the mixed two-point function can be calculated. It turns out that there are several possible solutions for , the smallest one being .

We want to stress that the solution obtained here solves the complete tower of DSEs in the IR and is the only possible solution. An important point is how the horizon condition enters. Without it we only could get the trivial result of canonical scaling of all Green function because of the bare two-point functions appearing on the right-hand sides of the two-point DSEs. However, solving the DSEs requires a renormalization condition which can be chosen as the value of the ghost dressing function at zero momentum. It was shown in ref. Fischer et al. (2009) that this value decides between the scaling solution (when it is chosen to be 0) and a family of decoupling solutions. Here the horizon condition explicitly forces us to use the former renormalization Zwanziger (2002); Gribov (1978); Zwanziger (1993). This is how this additional condition enters and why we get the scaling solution and no decoupling type of solution like, for example, in Fischer et al. (2009); Aguilar et al. (2008); Alkofer et al. (2010). In principle the analysis also allows for other possibilities, e.g., for an IR divergent gluon propagator, but they are not realized. Discussions about the different solutions can be found, e.g., in Fischer et al. (2009); Alkofer et al. (2010); Sternbeck and von Smekal (2010); Maas (2010).

Our results are in qualitative agreement with those obtained in refs. Gracey (2010); Zwanziger (2010), where some propagators of the auxiliary fields are found to diverge like . This corresponds to the value , which was already obtained by Gribov Gribov (1978) in a perturbative analysis. The agreement we find in the IR with the scaling solution of the Faddeev-Popov action is another indication of the solidity of this confinement scenario, since it corroborates the conjecture of Zwanziger Zwanziger (2002, 2004) about cutting the integration at the Gribov horizon.

Acknowledgments: Discussions with Daniel Zwanziger are gratefully acknowledged. M. Q. H. was supported by DFG projects Gi 328/1-4 and GRK 1523 and the FWF project W1203.

### Footnotes

- Speaker at the Quark Confinement and the Hadron Spectrum IX Conference, August 30 - September 3, 2010, Madrid

### References

- R. Alkofer, and L. von Smekal, Phys. Rept. 353, 281 (2001), hep-ph/0007355.
- J. Braun, H. Gies, and J. M. Pawlowski, Phys. Lett. B684, 262–267 (2010), 0708.2413.
- L. von Smekal, A. Hauck, and R. Alkofer, Ann. Phys. 267, 1 (1998), hep-ph/9707327.
- D. Zwanziger, Phys. Rev. D65, 094039 (2002), hep-th/0109224.
- C. Lerche, and L. von Smekal, Phys. Rev. D65, 125006 (2002), hep-ph/0202194.
- J. M. Pawlowski, et al., Phys. Rev. Lett. 93, 152002 (2004), hep-th/0312324.
- C. S. Fischer, A. Maas, and J. M. Pawlowski, Annals Phys. 324, 2408–2437 (2009), 0810.1987.
- V. N. Gribov, Nucl. Phys. B139, 1 (1978).
- D. Zwanziger, Nucl. Phys. B323, 513–544 (1989).
- D. Zwanziger, Nucl. Phys. B364, 127–161 (1991).
- D. Zwanziger, Nucl. Phys. B399, 477–513 (1993).
- D. Zwanziger, Phys. Rev. D69, 016002 (2004), hep-ph/0303028.
- M. A. L. Capri, et al. (2010), to be published in PRD, 1009.4135.
- D. Dudal, S. P. Sorella, and N. Vandersickel, Eur. Phys. J. C68, 283–298 (2010), 1001.3103.
- R. Alkofer, M. Q. Huber, and K. Schwenzer, Comput. Phys. Commun. 180, 965–976 (2009), 0808.2939.
- M. Q. Huber, R. Alkofer, and S. P. Sorella, Phys. Rev. D81, 065003 (2010a), 0910.5604.
- M. Q. Huber, K. Schwenzer, and R. Alkofer, Eur. Phys. J. C68, 581–600 (2010b), 0904.1873.
- R. Alkofer, C. S. Fischer, and F. J. Llanes-Estrada, Phys. Lett. B611, 279–288 (2005), hep-th/0412330.
- M. Q. Huber, R. Alkofer, C. S. Fischer, and K. Schwenzer, Phys. Lett. B659, 434–440 (2008), 0705.3809.
- A. C. Aguilar, D. Binosi, and J. Papavassiliou, Phys. Rev. D78, 025010 (2008), 0802.1870.
- R. Alkofer, M. Q. Huber, and K. Schwenzer, Phys. Rev. D81, 105010 (2010), 0801.2762.
- A. Sternbeck, and L. von Smekal, Eur. Phys. J. C68, 487–503 (2010), 0811.4300.
- A. Maas, Phys. Lett. B689, 107–111 (2010), 0907.5185.
- J. A. Gracey, JHEP 02, 009 (2010), 0909.3411.
- D. Zwanziger, Phys. Rev. D81, 125027 (2010), 1003.1080.