Gribov ambiguities at the Landau -- maximal Abelian interpolating gauge

In a previous work, we presented a new method to account for the Gribov ambiguities in non-Abelian gauge theories. The method consists on the introduction of an extra constraint which directly eliminates the infinitesimal Gribov copies without the usual geometric approach. Such strategy allows to treat gauges with non-hermitian Faddeev-Popov operator. In this work, we apply this method to a gauge which interpolates among the Landau and maximal Abelian gauges. The result is a local and power counting renormalizable action, free of infinitesimal Gribov copies. Moreover, the interpolating tree-level gluon propagator is derived.


Introduction
One of the most important and challenging open problems in theoretical Physics is the full comprehension of the non-perturbative features of Yang-Mills theories. Responsible for describing the successful Standard Model at high energies, Yang-Mills theories still lack a complete consistent quantization. As pointed out by V. N. Gribov [1] at the Landau gauge, a residual gauge symmetry survives the Faddeev-Popov gauge fixing procedure [2]. It is a known fact that, to quantize a gauge theory, it is necessary to consistently eliminate the gauge freedom of the Yang-Mills action, see also [3]. The residual gauge symmetry is characterized by the presence of redundant configurations (called Gribov copies) which still contribute to the path integral. A very important remark is that it is not a particular defect of Landau gauge, but of all covariant gauges, as formally shown by I. M. Singer [4]. Since these configurations represent a redundancy in the theory, their elimination is an unavoidable requirement.
Still in [1], Gribov showed that copies which are related by infinitesimal gauge transformations are associated with the zero-modes of the Faddeev-Popov operator (or, equivalently, to the poles of the ghost propagator) for the Landau gauge. In fact, this is true at least for all gauges that depend exclusively on the gauge field, see [5]. Moreover, Gribov proposed the elimination of these infinitesimal copies by restricting the path integral to a region which is free of infinitesimal copies. This region is known as the first Gribov region or, simply, Gribov region. Essentially, it is defined as the region where the Faddeev-Popov operator is positive-definite, a property that ensures that no infinitesimal copies are present. A very important feature is that all gauge orbits actually cross the Gribov region [6]. Then, since all physical configurations have at least one representative inside the Gribov region, the restriction is a consistent improvement of the Faddeev-Popov trick. The restriction of the path integral to the Gribov region implies on a dramatic modification of the gluon and ghost propagators. In one hand, the gluon propagator is suppressed at the infrared regime and acquires imaginary poles, on the other hand, the ghost propagator is enhanced by an infrared behaviour of the type 1/k 4 . These properties show, in an explicit way, that the elimination of the infinitesimal copies is of great importance for a consistent quantization, deeply modifying the theory and providing evidences of confinement.
The solution proposed by Gribov works nicely when the Faddeev-Popov operator is hermitian. The reason is that hermiticity ensures that the spectrum of the Faddeev-Popov operator is real and, therefore, it is possible to establish an order relation between the eigenvalues of such operator. Hence, it is possible to define a region where the Faddeev-Popov operator is positivedefinite and the restriction of the path integral to this region ensures the absence of infinitesimal copies. However, if we desire to work with non-hermitean Faddeev-Popov operators, to perform such restriction is not a clear procedure because the order relation and, therefore, the definition of a region, do not make sense anymore. In this sense, hermiticity of the Faddeev-Popov operator plays a fundamental role for the elimination of copiesà la Gribov. extra constraint. Therefore, all infinitesimal Gribov copies are eliminated at the classical level. In a certain sense, this elimination is direct and does not require the construction of a geometric region to restrict the path integral. The only requirement is to avoid all zero-modes, which characterize the infinitesimal copies. Since the identification of Gribov copies and the zeromodes of the Faddeev-Popov operator is independent of the hermiticity of this operator, this method should also be employable to treat gauges with non-hermitian Faddeev-Popov operators. Thus, this method brings a new perspective on the elimination of copies. It is important to recall that, the method developed in [5] requires exclusively A-dependent gauge conditions, although any particular choice was imposed. Therefore, in principle, there is a large class of gauges for which the method would be applicable. In particular, still in [5], consistency tests were made by applying the method to the Landau and maximal Abelian gauges. It is worth mention that the method here described is not the only alternative to the Gribov and Zwanziger techniques. There are other alternative techniques to deal with the Gribov ambiguities, see for instance [20,21].
In the present work, we apply this new method [5] to a gauge with non-hermitian Faddeev-Popov operator: the Landau -maximal Abelian interpolating gauge (LMAIG) [22,23,24]. This gauge has, at least, three advantages to motivate the present investigation. The first one, as already mentioned, is that the traditional approaches are not able, in principle, to deal with this gauge or any other gauge with non-hermitian Faddeev-Popov operators. Second, it is a gauge that link the two gauges where the Gribov problem can be handled. Thus, it is possible to verify the consistency of the results by interpolating among both limits of the LMAIG. Third, this gauge can be defined through a minimizing function given by where η is the interpolating parameter and the indices refer to the non-Abelian and Abelian sectors of the SU (N ) group (see Sec. 3 for the conventions). The gauge conditions of the LMAIG can be obtained by the minimization of the operator (1) with respect to gauge transformations. This means that the LMAIG could be, in principle, implemented on the lattice. This is a very welcome feature because it can work as a test for the application of the method 1 .
This paper is organized as follows: in Sect. 2, a brief review of the method developed in [5] is given. In Sect. 3, we provide a review of the decomposition into diagonal and off-diagonal components of algebra-valued quantities, present the maximal Abelian gauge and make the explicit decomposition of the Landau gauge. After this, we introduce the LMAIG and discuss its important features for the analysis of Gribov copies. Then, in Sect. 4, we apply the method to the LMAIG, and construct an action free of infinitesimal copies. This is done in Sect. 4. In Sect. 5 we calculate the diagonal and off-diagonal gluon propagators and show how it is possible to deform it into Landau and maximal Abelian gauges propagators. In Sect. 6 we make some comments about the gap equation in this method. Finally, in Sect. 7, we provide our conclusions. Many algebraic details were left to appendices to avoid big interruptions along the text.

A brief review of the method
The elimination method proposed in [5] is based on the introduction of an extra constraint that ruins the Gribov copies equation. In this section, we provide a brief review of the method in order to apply it to the interpolating LMAIG [23,24]. It is not our intent to be rigorous here. For any formal detail we refer to [5].
Let us consider Yang-Mills theory for a given semi-simple Lie group G. We choose a gauge condition ∆ A that depends exclusively on the gauge field, i.e. ∆ A = ∆ A (A), where the group indices vary as A, B, C, . . . ∈ {1, 2, . . . , dim G}. As pointed out by Gribov [1], the Faddeev-Popov gauge fixing procedure does not ruin completely the gauge symmetry. Thus, some redundant configurations, which are connected by gauge transformations, are still being considered in the path integral. The existence of these copies depends on the existence of solutions for the Gribov copies equation, and this is obtained by the requirement of gauge invariance of the gauge condition, i.e.
where U ∈ G and g is the coupling parameter. Besides the fact that we do not have much knowledge about the elimination of Gribov copies generated by large gauge transformations, we reasonably understand how to handle those copies generated by infinitesimal transformations. For this reason, we restrict ourselves to this case. The infinitesimal gauge transformation is then given by where ζ B is the infinitesimal gauge parameter. Thus, the copies equation (2) becomes where D AB µ ≡ δ AB ∂ µ − gf ABC A C µ is the covariant derivative and f ABC represents the structure constants. At first order in ζ, Eq. (4) can be written as where ∇ AB is the Faddeev-Popov operator, Summarizing, Eq. (5) is obtained by requiring infinitesimal gauge invariance of ∆(A).
The BRST transformation defined through the nilpotent operator s is given by where c A is the Faddeev-Popov ghost field, c A is the antighost field and b A is the auxiliary Nakanishi-Lautrup field. It is immediate to see that the first equation of (7) has the same form of (3). Of course, we have to understand that these are different transformations: The BRST, in particular, transforms a field with vanishing ghost number into a composite field with ghost number +1, while the gauge transformation does not change the ghost number. Nevertheless, it was proved in [5] that these transformations are homotopic. Thus, since they have the same formal structure, we can obtain the copies equation by requiring, not the gauge invariance of the gauge condition, but the BRST invariance 2 of the very same gauge condition. Hence, we can write the copies equation as The key point of the method resides at this stage: Since we want a theory free of copies, we have to ruin the copies equation. This might be seen as a new constraint for the theory. Thus, from Eq. (8), we can see that, to ruin this equation, we need to break the BRST invariance of the copies equation. In this sense, we want to write an equation such that where the term Ω A must prevent the theory to develop infinitesimal copies. Roughly speaking, this is the main idea behind the method. Now, in order to implement Eq. (9) in a gauge theory, we have to be careful to preserve all well established features of the perturbative regime. A very important requirement we have to make is that the BRST symmetry is restored at the perturbative regime. This means that the BRST breaking must be soft. Another requirement is that, since we do not want to affect the ghost sector, which is of great importance for the perturbative sector, we must introduce a set of trivial auxiliary fields to mimic Eq. (8). Finally, to ruin the copies equation and impose a consistent equation compatible with (9), the introduction of a soft BRST breaking term must be performed. As argued in [5], these goals are achieved by the introduction of two extra terms to the gauge fixed action, namely S triv and Ξ. These terms are responsible to implement a new constraint to the theory, satisfying the requirements mentioned before and reproducing Eq. (9). Thus, we impose the action where As stated before, the term S triv is introduced to mimic the copies equation. To do so, we introduce a BRST quartet in such a way that It is easy to see that the equation of motion for ϕ produced by S triv is exactly the copies equation. Moreover, the indices of the auxiliary fields are not arbitrary and describe the degeneracy of the copies equation. Since our point is precisely to ruin this equation, the term Ξ has the following general form where γ is a mass parameter introduced to fullfill the soft breaking requirement. With Ξ, we see that the equation of motion for ϕ is modified and represents a "ruined" copies equation which is the extra constraint that ensures the absence of infinitesimal copies.
The action given by Eq. (10) is then an action which satisfies the constraint given by Eq. (15). With this we eliminate all infinitesimal copies at the classical level. This result qualitatively coincides with the well establish refined Gribov-Zwanziger action, see [14]. It is worth mention that the form of the breaking term defined by eq. (14) has a kind of "freedom". To ruin the copies equation we must add a term which will be responsible for the breaking of BRST invariance of the gauge condition. This term must depend on ϕ for the obvious reason that, if it does not, the variation of the action (10) with respect to ϕ would not produce a "ruined" copies equation, as required. Moreover, the derivative of this term with respect to ϕ must depend exclusively on the gauge field A. The reason is that, if it depends on other fields, this term vanishes at their trivial vacua. Requiring the exclusive dependence on A, we ensure that the only copies that can be generated are related to A = 0. However, if they exist, they are necessarily different from zero and, therefore, the constraint will eliminate the copy A = 0 for the appropriate A. The conclusion is that the first term of Eq. (14), is sufficient for our requirements. In this sense, to ruin the copies equation in a minimal way, we could add only (16) to the original action and it will generate a theory free of infinitesimal copies. It this case, the extra terms can be included by the LCO technique, in the usual way [25,26,27,28,29]. On the other hand, once γ is at our disposal, the extra soft terms in (14) are allowed by power counting. What would decide if they are present or not are the Ward identities of the particular chosen gauge. In both cases, the effect is the obtention of the refined Gribov-Zwanziger action [14,15]. Furthermore, there is another possible freedom, for each term proportional to the mass parameter γ, we could replace it by different mass parameters. Essentially, this can also be obtained by the redefinition m i = ζ i γ, which means that the independent character of these coefficients are accounted by the parameters ζ i . Let us also comment on the term proportional to ζ 1 . As it is possible to see from eq.(14), a larger field combination is associated with the parameter ζ 1 . The reason is that we can introduce such combination as a BRST exact form, γ 2 s d 4 x ω AB µ ϕ AB µ , for instance. We could introduce this mass terms in a independent way, but following the idea of a minimal breaking of the BRST symmetry, a BRST exact term fits better for our plans.
Finally, it is important to understand that the inclusion of all extra terms proportional to γ implies on a deep modification of the so-called gap equation [5]. Until now, no results are known for this generalized gap equation and we are not able to decide if it is a better choice to follow it or not. However, we will keep these terms as in (14) because they are important to reproduce the refined Gribov-Zwanziger features and also because we will not deal with the gap equation in this work 3 . In fact, if the extra terms are not allowed for any reason, all we have to do is to set the corresponding ζ i to zero.

The Landau and maximal Abelian gauges and their interpolation
From now on, we restrict ourselves to the SU (N ) gauge group. In [23,24], a gauge fixing which interpolates among Landau, Coulomb and maximal Abelian gauges was studied. In the present work, we will analyze the Gribov problem in this gauge. However, we will consider only the interpolation between Landau and the maximal Abelian gauges and avoid the Coulomb sector of the gauge. Since the MAG is characterized by imposing different gauge conditions to the diagonal and off-diagonal components of the Lie algebra-valued fields, we will decompose the Landau gauge in order to provide an explicit comparison with the reduction of the interpolating gauge to the Landau case 4 . To fix notation and conventions, we will briefly review this kind of decomposition, called Abelian decomposition [18]. Essentially, the SU (N ) group is dismembered into its Abelian and non-Abelian sectors where the Abelian sector is recognized as the Cartan subgroup. The gauge field decomposition is taken as where G A correspond to the (N 2 − 1) generators of the SU (N ) group; G a are the N (N − 1) off-diagonal generators of the gauge group; and G i represent the (N − 1) Cartan subgroup generators. The indices {a, b, c, . . . , h} run from 1 to N (N − 1) and the indices {i, j, k, . . .} run from 1 to (N − 1). As a consequence of this decomposition, we can write the decomposed BRST transformations (7) as where the covariant derivative is defined with respect to the Abelian sector and acts on non-Abelian quantities, We can now write the gauge fixed Yang-Mills action (11) as where F a µν and F i µν are the components of the field strength, which are, explicitly, and ∆ a (A) and ∆ i (A) are related to the gauge condition components. To complete the Abelian decomposition we can write the Jacobi identities as

The maximal Abelian gauge
The maximal Abelian gauge imposes different gauge conditions to the diagonal and off-diagonal sectors of the gauge fields, namely and the corresponding gauge fixing action is given by where the operator ∇ ab is the Faddeev-Popov operator, The gauge conditions (23) can be obtained from the equations of b a and b i . If we think of the conditions for the existence of Gribov copies, we can derive the copies equation requiring the gauge/BRST invariance of the gauge condition [5]. Hence, since we have two different gauge conditions in MAG, it is natural to expect two copies equations. In fact, if we calculate directly the copies equation from the gauge conditions (23), we obtain the following equations where ζ a and ζ i are the off-diagonal and diagonal components of the infinitesimal gauge parameter, respectively. From (26) we see that the first equation just involves the off-diagonal component of the gauge parameter while the second involves both. Simple manipulations of the second equation provide which shows that, once one has solved the first equation of (26), the second does not contribute with any extra information. This redundancy is the reason why only the first equation of (26) is considered as the copies equation for the MAG. A final comment is that the Faddeev-Popov operator is hermitian in this case, see [16,17] and references therein for more details.

The decomposed Landau gauge
The Landau gauge condition does not distinct the diagonal and off-diagonal sectors of the gauge connection. However, since we will work with decomposed fields, we also write the Landau gauge fixing relevant expressions in the Abelian decomposition. The result is the decomposed Landau gauge fixing action, given by It is immediate to obtain the Faddeev-Popov operator from (29). In components, it is given by Unlike the case of MAG, all components above in (30) contribute to the copies equations. If we again consider the diagonal and off-diagonal components of the infinitesimal gauge parameter, ζ i and ζ a , respectively, we can write the following copies equations In this case, both equations encompass all components of the infinitesimal gauge parameter. Hence, we cannot put away any of them and all components of the Faddeev-Popov operator are essential to the analysis.
Another important remark is that the full Faddeev-Popov operator for the Landau gauge is also hermitian, see for instance [1,3]. If we adopt a matrix viewpoint, an hermitian operator is such that the elements of its diagonal are hermitian operators and all elements above the diagonal are hermitian conjugate of the elements below it. If we analyze the mixed components of (30) we can see that (∇ ai ) T * = ∇ ia .

Interpolating gauge
In order to provide a gauge fixing which interpolates among Landau and maximal Abelian gauges [23], we introduce a real interpolating parameter η and write the following gauge conditions Thus, it is clear that the gauge condition for the diagonal component of the gauge field is identical for the Landau and maximal Abelian gauges cases. Moreover, for the first equation of (32), the case η = 1 gives the Landau gauge condition while for η = 0, the MAG condition is achieved. Consequently, we can write the gauge fixing term as It is a simple exercise to verify that for η = 0, Eq. (33) reduces to Eq. (24), and for η = 1, it reduces to Eq. (29).
4 Eliminating the Gribov copies

The Faddeev-Popov operator and Gribov ambiguities
The gauge conditions presented in the last Section provide a way to interpolate among the Landau and maximal Abelian gauges. To analyze the Gribov problem in this gauge, it is fundamental to study the Faddeev-Popov operator in order to establish the copies equation and its main properties. The way we do this is completely analogous to the MAG (as showed in Sect. 3.1). Thus, by requiring gauge/BRST invariance of (32), the following operators are obtained These operators act on a gauge parameter pair (ζ a , ζ i ), exactly as in (31). First of all, as a consistency check, we must verify if, for the suitable choices of the parameter η, this operator reduces to the previous operators for Landau and maximal Abelian gauges. Hence, starting with η = 0, the first equation of (34) reduces immediately to (25), the second turns out to be the null operator and the last remain unaffected. Actually, the two last equations simply define the redundant condition (27) for the MAG. We conclude then that the choice η = 0 returns the Faddeev-Popov operator of the MAG, as expected. On the other hand, if we choose η = 1, we can see that, after some simple manipulations, the first equation of (34) reduces to ∇ ab = −∂ µ D ab µ − gf acb A c µ ∂ µ , which is exactly the purely off-diagonal components of the Faddeev-Popov operator for the Landau gauge. The second equation of (34) reduces to the second equation of (30). The third equation of (34) does not involve the interpolating parameter η, but once we choose η = 1, the first gauge condition of (32) turns out to be ∂ µ A a µ = 0, which means that we can substitute this in the third equation of (34) and obtain the same result for the Landau gauge (30). Finally, the fourth equation is unchanged. Summarizing, these are the components of the Faddeev-Popov operator for the Landau gauge. This concludes our consistency checks for now.
Let us make a quick remark about the Faddeev-Popov operator. It is widely known that, in standard techniques employed to deal with the Gribov problem [1,3], the hermiticity of the Faddeev-Popov operator is essential. Since we can associate Gribov copies with zero-modes of the Faddeev-Popov operator, the knowledge about its spectrum is very welcome. A hermitian operator has only real eigenvalues, allowing to establish an order relation between them. In the Landau gauge, for instance, is through this analysis that it is possible to construct a region where the Faddeev-Popov operator is positive-definite. For this reason, it is possible to eliminate all infinitesimal Gribov copies from the path integral by the restriction of the integration to this domain. This technique makes the analysis of the Gribov problem a geometrical problem. On the other hand, this method cannot be employed to non-hermitian Faddeev-Popov operators. It is then not clear how to generate a region that will restrict the path integral. Nevertheless, the method developed in [5] and briefly reviewed in Sect. 2 does not require the definition of a region to perform the functional integration. Thus, we can apply it for gauges with non-hermitian Faddeev-Popov operators. In fact, as discussed in [5], in the case of hermitian Faddeev-Popov operators, the new method is equivalent to restrict the path integral to a region defined by the zero-modes of the corresponding FaddeevPopov operator.
Getting back to the LMAIG Faddeev-Popov operator (34), its first decomposed operator can be rewritten as is the Faddeev-Popov operator for the MAG. As mentioned in Sect. 3.1, the operator defined by Eq. (36) is hermitian. It is possible to show that the operator ∇ ab −∇ ab is not hermitian, the details of the proof can be found at Ap. A. The purely diagonal component of the Faddeev-Popov operator is trivially hermitian. Now, only the mixed components are left and we will follow an analogous idea presented in Sect. 3.2. The whole idea is based on the fact that we can write the Faddeev-Popov operator in a matrix form and it has two blocks formed by the purely off-diagonal and diagonal components. Terms outside these blocks are the mixed ones and to analyze their hermiticity we must take their hermitian conjugate. Thus, transposing and taking the complex conjugate of the matrix it is clear that (∇ ai ) † = ∇ ia . In fact, since one involves the parameter η and the other does not, it is impossible to establish a hermiticity relation between these components. Thus, pictorially, we have that the full Faddeev-Popov operator, is not hermitian. Moreover, unlikely the MAG, it is not possible to eliminate some components of this operator to analyze the Gribov copies. See (31). This, and the fact that the LMAIG is an exclusively A-dependent gauge, are the leads allowing the direct elimination of the zero-modes within the method developed in [5].

Trivial set of auxiliary fields
According to the method, it is possible to eliminate Gribov copies directly, by imposing a new constraint to the theory. This constraint is, essentially, the requirement that the copies equation is not obeyed. This is done by the introduction of a set of auxiliary fields, forming a BRST quartet, through a trivial term and a soft BRST breaking term. From now on we will deal exclusively with the interpolating gauge, so when we refer to the Faddeev-Popov operator, we are talking about the operator defined by Eq. (34).
The trivial term is given by 5 where capital latin indices refer to the complete Lie algebra. The full decomposition of (38) in off-diagonal and diagonal components results in where terms involving the functional derivative with respect to A of ∇ ij and terms involving the functional derivative with respect to A j of ∇ ia and ∇ ai are not present because they vanish.
The explicit form of (39) is A consistency check must be performed to verify if this trivial term interpolates among Landau and maximal Abelian gauges trivial terms [5]. To avoid many tedious algebraic steps along the text, we leave this proof for the Ap. B.

Breaking Term
We have now to introduce a soft BRST breaking term at the original action S 0 , given by The reason is the following: Since we can obtain the Gribov copies equation requiring the BRST invariance of the gauge condition, to ruin this equation we must break the BRST invariance. This breaking, however, is not arbitrary. Therefore, when we look for perturbative effects of the theory, the BRST invariance must be restored, and for this reason, we call this a soft breaking [30,31]. A soft breaking can be obtained by the introduction of a mass parameter γ which makes the dimension of this term lower than the spacetime dimension. See [5] for more details about this construction. The general form of the soft breaking term is where the parameters ξ, θ, χ and ζ i must be η-dependent in order to permit the interpolation of the breaking term of Landau and maximal Abelian gauges. Such interpolation is presented at Ap. C. Comparing eq.(42) with eq.(67) and eq.(68), we obtain with ζ 1 and ζ 2 being independent from η. We remark that, depending on the Ward identities, some of these parameters might be zero. Hence, the breaking term Ξ with the appropriated fixed parameters is We remark that this term not only breaks the BRST in a soft manner, but also ensures that the copies equation is ruined. As discussed in Sect. 2, this breaking term has some sort of "freedom". In order to write an action free of infinitesimal copies, we could just introduce the following term The point here is that the terms given by Ξ −Ξ permit a construction very close to the refined Gribov-Zwanziger action [14]. In fact, to make contact between Eq.(44) and the LCO formalism, we should write the mass terms as independent masses m i = ζ i γ 2 and deal with local composite operators and their condensation. An immediate difference between both methods relies on the gap equation, see Sect. 6.

Gluon Propagator
It is remarkable that the elimination of infinitesimal Gribov copies, as an apparent technicality, brings rich effects to the physical properties of non-abelian gauge theories. Of course, they play a fundamental role for a consistent quantization, but their elimination provides a profound change in the gluon and ghost propagators, specially at the infrared regime. This is a well-known feature for Landau and maximal Abelian gauges, see for instance [1,16]. In fact, the inclusion of dimension 2 condensates makes the analytic result of the propagators to stay in harmony with lattice results [12,13,34,35]. In this section, we compute the off-diagonal and diagonal gluon propagators for the interpolating gauge. As mentioned before, this provides a good way to test the free of copies action presented here, since this gauge can be implemented in the lattice. An interesting feature to study is the deformation of the propagators of the Landau gauge into the propagators of MAG, a property that could be investigated in the lattice.
The full action S is given by For the gluon propagator, only the quadratic action S q (A) is required 6 , where α and β are gauge parameters. The actual interpolating gauge is obtained at the limit α = β = 0. Taking the Fourier transform of Eq. (47), we obtain the following expression, It is not a difficult task to obtain the diagonal and off-diagonal gluon propagators from (48). One has only to invert the corresponding wave operators in the usual way. The expressions for the diagonal and off-diagonal gluon propagators are, respectively, where the limits α = β = 0 were already taken. If we rename each mass term 7 which appears with a ζ i parameter in terms of independent masses, we can write these propagators as Since we are dealing with the interpolating gauge, we have to check the deformation of the propagators (49) and (50) among Landau and maximal Abelian gauges for η = 1 and η = 0, respectively. It is obvious that the diagonal gluon propagator is independent of the parameter η, which means that its form does not change for intermediary choices of η. It is known that this propagator is identical in Landau and maximal Abelian gauges [14,17], and the reason why it has the same form for a interpolation among them is that the gauge fixing condition for the diagonal components remains invariant as η varies. So, we have to analyze only the off-diagonal sector in order to provide an explicit test of its reduction to the known cases. For η = 0, the propagator (50) reduces to which is exactly the off-diagonal gluon propagator for the maximal Abelian gauge [17]. Choosing η = 1, we have which coincides with (49), as expected. Of course, the new features that all these calculations bring reside on different values for η than 0 and 1. Hence, since η ∈ [0, 1], we can choose η in such a way that an explicit continuous deformation of the off-diagonal propagator can be seen. Another remark is that for an arbitrary value of η different from 0 and 1, we have the propagators for a gauge with non-hermitian Faddeev-Popov operator. This is a very important feature because, in usual approaches [1,7,14], the construction of these propagators would be very difficult (if not impossible). We also remark that expressions (51) and (52) are much more related to the usual results of the refined Gribov-Zwanziger approach where the mass parameter are independent of the Gribov parameter (at least in a tree-level analysis).
The deformation of the off-diagonal propagator is displayed in Fig. 1 (Both pictures are exactly the same, except for the viewpoint angle). From this picture, it is possible to see the 7 In order to match with the usual conventions, we must rename the terms involving ζ1 with a minus, i.e. −γ 2 ζ1 = m 2 1 . This is very important, because the opposite choice will lead to negative m 2 1 values. On the other hand, the terms γ 2 ζ2 are correctly identified with m 2 2 . explicit continuous deformation from the MAG off-diagonal gluon propagator to the Landau offdiagonal gluon propagator. A very well-known feature is that, for k = 0, the gluon propagator does not vanish, at both gauges. It is interesting, however, that this value has some kind of oscillatory behavior in the intermediary choices of η and this is explicit in the right plot of Fig. 1. There is an absolute maximum at η = 1/3 while η = 1 is the absolute minimum (obviously, the domain we are considering is η ∈ [0, 1]), which corresponds to the Landau gauge at zero momentum. As a curious fact, in the case of N = 2, there is a term in the denominator of Eq.(52) which vanishes. The consequence is that the the curve on the k = 0 plane loses the oscillatory behavior and becomes a constant with respect to η, as shown in Fig. 2

A few words about the gap equation
As discussed in [5], the method applied in this work can provide a generalized gap equation. The reason why it "can" lies on the fact that we have some freedom in the choice of the breaking term, as discussed in Sect. 2. First, if we do not choose the mass terms which give rise to the refined Gribov-Zwanziger action, no modification to the gap equation emerges. We must remind that, if we opt for a minimal breaking of the BRST symmetry, these terms are excluded (at least the A 2 and cc which are not BRST invariant) and can be considered through the LCO formalism [14,15]. On the other hand, we can include these terms for two reasons: (i) they are permitted by power counting and dimensional analysis and (ii) we recover the refined Gribov-Zwanziger propagators independently from the gap equation form. We could include them with an explicit dependence on γ and, since the gap equation is obtained by minimizing the quantum action with respect to γ 2 , all these terms will contribute to it, see [5]. This is the difference to the usual gap equation which does not contain these terms. Alternatively, these terms could also be included with independent mass parameters (i.e., with no explicit dependence on γ) and the gap equation would not be affected. The main importance of these possibilities is that they can be very welcome, since the usual gap equation (at the Landau gauge, for instance) throws the theory right at the horizon, which is precisely the place where infinitesimal copies start to appear. The decision of how good is an alternative gap equation will rise with its physical effects and consistency checks. For obvious reasons, this analysis is left for future investigation.

Conclusions
In this work we applied the method developed in [5] to eliminate infinitesimal Gribov copies from the interpolating Landau -maximal Abelian gauge. We obtained an action free of copies, given by Eq. (46), which has the same structure of the refined Gribov-Zwanziger-type actions [14,15]. After that, with a suitable choice of the interpolating parameter, we extracted the diagonal and off-diagonal gluon propagators and showed that the results reduce to the well-known propagators for the Landau and maximal Abelian gauge fixings.
Although the elimination of infinitesimal Gribov copies for the interpolating gauge is important and interesting by its own merit, it brings some new features for the general problem of dealing with Gribov copies. As discussed throughout Sect. 4.1, this gauge has a non-hermitian Faddeev-Popov operator, which means that no order relation can be established between its eigenvalues. Consequently, the possibility of constructing a well defined free of copies region in functional space is not evident. In this sense, the elimination performed here through the method developed in [5] opens a new door to the understanding of non-hermitian Faddeev-Popov operators. Moreover, the propagators computed in Sect. 5 are already in a form to be compared with a possible lattice simulation of the LMAIG.
There are many issues that should be addressed now. All of them deserves investigation. However, each of them is quite extend and intricate, and therefore are beyond the goals of this work. Nevertheless, they are left for future investigation. To cite a few interesting topics to be investigated, we can start with the renormalizability problem of action (46). As in the case of the MAG, many complications and extended expressions, due to extra quartic ghost interacting terms, are expected. Another important issue is the Abelian and non-Abelian ghost propagators, another task that should demand a laborious amount of computations (a smart attack would be to start with the SU (2) case). A third problem to be studied is the comprehension of what could be, if any, an analogous Gribov region for this gauge and the interpolation between the known regions of Landau and maximal Abelian gauges. Finally, as discussed in Sect. 6, the effect of a possible alternative gap equation has to be taken under consideration. This last question opens the possibility of introducing the refined mass parameters directly on the gap equation as a function of the Gribov parameter instead of a independent local composite operators condensation. Obviously, we would start this study at the Landau gauge, which is the gauge where the Gribov problem and its effects is most understood.

A (non-)Hermiticity of ∇ ab
Let's define, using Eq. (35), an operator ∇ given by Since∇ is the MAG Faddeev-Popov operator, it is hermitian. Thus, if we prove that ∇ is hermitian, it is sufficient to say that ∇ is also hermitian. To study this possibility, let us consider the following expression 9 Now, we will consider the three terms of the rhs of Eq.(56) separately. The first term is which proves its hermiticity. The second term is given by from where we can see that it is also hermitian. Finally, the third term is written as Clearly, the third term is not hermitian. Obviously, the reason is the presence of the piece in Eq. (59). Since the sum of hermitian operators is a hermitian operator, we can see from expression (56) that ∇ would be hermitian only if (60) vanishes. The conclusion is that the purely off-diagonal components of the Faddeev-Popov operator of the Landau -MAG interpolating gauge is not hermitian, except for the MAG limit. In the case of the Landau gauge, this operator combines itself with the other sectors in order to provide another hermitian Faddeev-Popov operator.

B Trivial terms -Maximal Abelian Gauge
As discussed in Sect. 4.1, the Faddeev-Popov operator of the interpolating gauge has four different sectors: One purely off-diagonal, one purely diagonal and two mixed operators. The Landau gauge has this same feature and the comparison between the trivial terms of these gauges can be done term by term. Since we have seen that the operator given by Eq. (34) reduces to the Faddeev-Popov operator at the Landau gauge for η = 1, the expressions of the trivial term of the interpolating gauge, for η = 1, must also coincide with the trivial term of the Landau gauge (which must be decomposed for a comparison). This is easy to verify. On the other hand, the usual MAG [16,17] has only off-diagonal components for the Faddeev-Popov operator. In this sense, we have to be careful with the trivial term because, when we choose η = 0 for the interpolating parameter, the mixed and purely diagonal components of the Faddeev-Popov operator will not vanish and this will provide a different trivial term with respect to the usual MAG case. Here, we have to remember that all mixed components of the Faddeev-Popov operator of the usual MAG are eliminated from the very beginning of the analysis of the Gribov problem (due to the redundant condition), see Sect. 3.1. Hence, the trivial term of the MAG must have the form S MAG triv = − d 4 x ϕ ac µ ∇ ab ϕ bc µ − ω ac µ ∇ ab ω bc µ − ω ac µ (sA d ν )
(68) Both expressions are consistent with those obtained in [5].