Infrared saturation and phases of gauge theories with BRST symmetry

We investigate the infrared limit of the quantum equation of motion of the gauge boson propagator in various gauges and models with a BRST symmetry. We find that the saturation of this equation at low momenta distinguishes between the Coulomb, Higgs and confining phase of the gauge theory. The Coulomb phase is characterized by a massless gauge boson. Physical states contribute to the saturation of the transverse equation of motion of the gauge boson at low momenta in the Higgs phase, while the saturation is entirely due to unphysical degrees of freedom in the confining phase. This corollary to the Kugo–Ojima confinement criterion in linear covariant gauges also is sufficient for confinement in general covariant gauges with BRST and anti-BRST symmetry, maximal Abelian gauges with an equivariant BRST symmetry, non-covariant Coulomb gauge and in the Gribov–Zwanziger theory.


Introduction
In their seminal work, [1], Kugo and Ojima develop the covariant operator formalism for gauge theories in linear covariant gauge. On the assumption of an unbroken BRST symmetry, they construct the physical Hilbert space of the theory and formulate a criterion for color confinement. The Hilbert space of a gauge theory is defined by the cohomology of s, the nilpotent generator of BRST transformations, The confining phase of a gauge theory according to [1] is characterized by an unbroken global color symmetry and J a μ is related to the canonical Noether current j a μ by the quantum equation of motion (QEoM) of the gluon, where S is the gauge-fixed action. According to Kugo and Ojima, color confinement is realized if neither of the two currents G a μ = −∂ ν F a νμ and N a μ = is(D μc ) a (4) create massless excitations. The corresponding charge operators are then both well defined. The global color charge Q a = G a +N a then also is well defined and vanishes on the physical Hilbert space. Following [1,2], we introduce the function u( p 2 ) by the correlator where L μν = p μ p ν / p 2 and T μν = δ μν − L μν are longitudinal and transverse projectors, and FT means Fourier transform. The BRST-exact charge N a is well defined only if u( p 2 ) Together with the absence of massless vector bosons, Eq. (7) is a sufficient condition for color confinement [2]. Here we express the confinement criterion of Kugo and Ojima in terms of the saturation of the gluonic QEoM, at vanishing momentum. With the classically conserved current of global color symmetry j a μ of Eq. (3), this equation in linear covariant gauge reads If G a μ does not create a massless vector boson, the first correlation function in Eq. (9) vanishes in the infrared limit p 2 → 0 and the current, must saturate Eq. (9) at vanishing momentum. The current j b μ is physically equivalent to the classical color current j b μ because they differ by a BRST-exact term only. Condition (7) for the correlator (6) implies that Thus: If the Kugo-Ojima criterion is fulfilled, the gluonic QEoM Eq. (9) in linear covariant gauge at vanishing momentum is saturated by BRST-exact states. Since the Physical Hilbert space does not include BRST-exact states, only unphysical states contribute to the saturation of the gluonic QEoM at vanishing momentum. We therefore have the following Proposition: In the confining phase of a gauge theory, unphysical states created by the color currentj a μ saturate the gluonic QEoM at vanishing momentum.
Since confinement allows only color singlet asymptotic states and any asymptotic state that contributes to the gluonic QEoM at vanishing momentum is in the adjoint color representation, the proposition clearly holds. An adjoint multiplet of physical asymptotic states on the other hand exists in Higgs and Coulomb phases. Thus one can discriminate between the Higgs and confinement phase in linear covariant gauge by whether or not physical states contribute to the matrix element of the (generalized) color currentj b μ at vanishing momentum.
It is another matter to prove that the gluonic QEoM indeed is saturated by unphysical states at vanishing momentum. Although we do not show this, we find that BRST-exact states in principle can saturate the gluonic QEoM at vanishing momentum in various gauges with a BRST or equivariant BRST symmetry. In particular, we verify in these gauges that (i) if the theory does not describe the Coulomb phase, the gluonic QEoM at vanishing momentum is saturated by the matrix element with a current that is physically equivalent to the conserved color current, and (ii) in non-Abelian gauge theories BRST-exact terms exist that can saturate the gluonic QEoM at low momentum.
Since physical states contribute to the gluonic QEoM at vanishing momentum in Higgs and Coulomb phases, we conclude that saturation of the gluonic QEoM at vanishing momentum by a BRST-exact term of the generalized current is a sufficient condition for confinement, provided that the BRST (or equivariant BRST) symmetry is unbroken.
In this article we identify the generalized current and the unphysical BRST-exact term in the gluonic QEoM in various gauges as well as in the Gribov-Zwanziger theory. The emergence of a similar pattern in all these models supports our proposition.
In a lattice theory with fundamental scalars, the Higgs and confining phases at finite lattice coupling β are analytically connected [3][4][5]. The situation is akin to a vapor-liquid transition below the critical point, and a (first-order) transition does seem to occur at sufficiently large β. It is therefore not clear whether the two phases remain analytically connected in the continuum limit. Since color charge is screened by a Higgs (or quark) field in the fundamental representation, the asymptotic behavior of the Wilson loop is expected to always be perimeter-like and cannot be used to distinguish between phases. As the liquid-vapor transition below the critical point demonstrates, the absence of an order parameter does not necessarily imply that the free energy is analytic everywhere. Gauge-dependent criteria for Higgs and confining phases give different critical curves below the critical βvalue, but are remarkably consistent above this critical point [5]. Analogous results have been obtained in a semi-classical continuum analysis of a non-Abelian Higgs model within the Gribov-horizon [6,7].
However, to unambiguously distinguish between a Higgs and a confining phase, we in this article consider only Yang-Mills theory without fundamental matter. This gauge theory is either in a Coulomb, a Higgs or a confining phase and the asymptotic behavior of the Wilson (or of the 't Hooft) loop [8] distinguishes between the last two. The investigation of transitions between these phases is beyond the scope of the present article.
We will consider Yang-Mills theory in linear covariant (LCG), generalized linear covariant (GLCG), maximal Abelian (MAG), Coulomb (CG) and minimal Landau (GZ) gauge. Some known results for these gauges are summarized below.
The best investigated Linear Covariant Gauge (LCG) is Landau gauge. Another corollary of the Kugo-Ojima confinement criterion [2] in this gauge relates the infrared limit of the ghost dressing function G( p 2 ) to u( p 2 ), lim p 2 →0 G( p 2 ) −1 = 1 + lim p 2 →0 u( p 2 ). When Eq. (7) is fulfilled, the ghost dressing function diverges in the infrared, and the ghost propagator at vanishing momentum is more singular than a massless pole. An exact infrared analysis of the whole tower of Dyson-Schwinger (DSE) and of Exact Renormalization Group Equations (ERGE) confirms the existence of solutions with this infrared behavior of the ghost propagator and a corresponding infrared suppressed gluon propagator; see for instance Refs. [9][10][11][12][13][14][15][16][17][18]. The solutions for which the Kugo-Ojima criterion is fulfilled have a power-like enhancement of the ghost propagator with related infrared exponents of other Green functions determined by an infinite tower of scaling relations [15,16]. According to numerical studies of Yang-Mills theory in 2d [19][20][21], the ghost and gluon propagators exhibit this scaling behavior, and, moreover, a strict analytic bound [22] implies that in the latter case the gluon propagator vanishes at zero momentum However, in 3d and 4d another one-parameter family of solutions to the DSEs and ERGEs also exists that is best parameterized by the value of G(0) −1 = 0. The ghost propagator of these solutions is only quantitatively enhanced and the gluon propagator is infrared finite [17,18,23,24]. Lattice gauge theory studies in three and four dimensions observe only this infrared finite behavior [25][26][27][28] of lattice propagators. But a whole one-parameter family of solutions can also be obtained on the lattice by tuning the lowest eigenvalue of the Faddeev-Popov operator [29] of the gauge fixing. The origin of this multitude of solutions is unresolved and it has been suggested [25] that the value of G(0) −1 could be considered an additional gauge parameter. Numerical solutions over the whole momentum range are only available for truncated DSEs and ERGEs.
Within the error of the employed approximation and/or truncation, both types of solutions lead to very similar phenomenology [30] and confine static quarks [31]. Unbroken BRST symmetry is essential for the Kugo-Ojima confinement criterion. Without recourse to a preserved nilpotent symmetry it is difficult to identify the unphysical sector of such truncated models. However, our proposition that confinement is a Higgs mechanism in the unphysical sector of the theory can be formulated in the absence of a nilpotent symmetry and may hold in all these scenarios.
In the minimal Coulomb gauge, the dressing function of the ghost propagator of Yang-Mills theory is more divergent than a simple massless pole, effectively leading to a confining color-Coulomb potential [32]. The infrared divergence of the instantaneous ghost propagator in this gauge has been verified by other calculations in the continuum [33][34][35] and on the lattice [36]. For a thorough discussion of the current status see [37].
The dual superconductor picture of the QCD vacuum [38][39][40] is the motivation for considering Yang-Mills in maximal Abelian gauge (MAG), [41][42][43][44][45]. This gauge discriminates between the Cartan subalgebra and the coset space of the gauge group and the partial gauge fixing breaks the local SU(N ) invariance down to the Cartan subgroup U (1) N −1 . The hypothesis of Abelian dominance [46] states that the Cartan gluons dominate long-range interactions. This has been observed in lattice simulations [47,48] and is also corroborated by an infrared analysis of the functional equations [49,50]. Furthermore, the Cartan gluons also dominate at large momenta [51,52]. A detailed understanding of the relation between the infrared behavior of Green's functions and confinement nevertheless is lacking in the MAG. It is of interest that a renormalization group analysis of interpolating gauges found that Abelian gauges form an invariant subspace that is not analytically connected to linear covariant gauges [53]. This explains why a literal interpretation of the Kugo-Ojima confinement criterion fails for this class of gauges [54]. Our proposition that unphysical states saturate the gluonic QEoM at vanishing momentum can nevertheless be realized. Contrary to Abelian gauge theories, the Abelian current of non-Abelian gauge theories includes an operator that only creates unphysical states that could saturate the gluonic QEoM.
The Gribov-Zwanziger framework [55][56][57][58][59] restricts the path integral to the first Gribov region of LCG and Coulomb gauge. Even though it drastically changes the gauge-fixed action and breaks BRST symmetry spontaneously, [60,61], it does not change the form of the DSEs in Landau gauge [32], and the infrared exponent of the scaling solution is that of ordinary Landau gauge [62,63]. We find that the color current of this model includes a BRST-exact contribution. The latter in fact saturates the gluonic QEoM at vanishing momentum already at tree level. However, it is difficult to verify that this BRST-exact operator only creates unphysical states since the BRST symmetry of this model is spontaneously broken.
Yang-Mills theory is expected to confine in all these gauges, and one hopes that the underlying mechanism can be characterized in a gauge-invariant fashion. Although the dynamics may be different, our proposition-concerning saturation of the gluonic QEoM at vanishing momentum by unphysical states-can be realized in all of them. Other similarities include that the dielectric function of the QCD vacuum appears to be related to the divergent dressing function of the ghost propagator in Coulomb gauge [64], and that in the Landau gauge the Kugo-Ojima criterion is related to the Gribov-Zwanziger scenario [65].
The present article is organized as follows: In Sect. 2, we examine the gluonic QEoM in Abelian gauge theory in linear covariant gauges, and we review the arguments that distinguish between Coulomb and Higgs phases in this specific case. In Sect. 2.1, we examine the Abelian Coulomb phase in greater detail, and in Sect. 2.2 we explicitly verify the implications of a spontaneously broken U (1)-symmetry in the t' Hooft gauge. In Sect. 3.1 the Kugo-Ojima confinement criterion for LCG is reviewed. We then proceed to generalize and adapt this criterion to other gauges: to generalized covariant gauges in Sect. 3.2, to covariant but nonlinear MAG in Sect. 3.3, and to the non-covariant Coulomb gauge in Sect. 3.4. In all these gauges a Kugo-Ojima-like confinement criterion is formulated. In Sect. 3.5 we examine signatures of confinement in the Gribov-Zwanziger (GZ) theory. Section 4 summarizes our results. Conventions and some technical details are deferred to three appendices.

Abelian gauge theories
In Abelian gauge theories one can, of course, dispense with a BRST construction of observables [66,67]. However, identifying the physical sector by a BRST symmetry is readily extended to non-Abelian gauge theories, and to noncanonical quantization.
To this end we consider Abelian gauge theories in general linear covariant gauges, Here ξ is a gauge parameter and b, c, andc are the Nakanishi-Lautrup (NL) and (anticommuting) ghost and antighost fields. The local function γ (. . . ) of canonical dimension 2 and vanishing ghost number is a polynomial of bosonic matter fields φ that does not depend on the gauge connection A μ or the NL field b. The matter part, L M , is invariant under U (1)-gauge transformations and may include covariantly coupled charged fermions and bosons. The variation s generates the nilpotent BRST symmetry [68,69] of L U (1) (12), Under the BRST variation s, charged matter fields vary by an infinitesimal U (1)-gauge transformation with the ghost field c(x) as variation. The longitudinal gauge field, ghost c, antighostc and NL field b form the elementary BRST quartet [1,70]. We note that one can define an anti-BRST variation in this Abelian setting bȳ with an obvious extension to matter fields. The generators of BRST and of anti-BRST transformations are nilpotent and anticommute, s 2 =s 2 = {s,s} = ss +ss = 0. The conserved BRST and anti-BRST charges corresponding to the transformations in Eqs. (13) and (14) in the Abelian case may be represented by the functional derivative operators, where the sums run over all matter fields φ M (x). The ghost number, also is conserved. The nilpotent BRST symmetry allows one to define the subset P of physical operators by the cohomology [68], Using a canonical construction, it was shown [1] that negative-norm states associated with asymptotic BRST quartets are unphysical. The elementary quartet thus is not observable. Transversely polarized photons, on the other hand, are physical. Instead of constructing the physical asymptotic Hilbert space directly, we prefer to define the space of physical operators, which is a slightly more flexible point of view that can be extended to spaces other than four-dimensional Minkowski spacetime. With a BRST invariant ground state, the two approaches are equivalent in Minkowski space. Since matter transforms covariantly under U (1)-gauge transformations, A μ → A μ + ∂ μ θ , and L M is an invariant, the conserved global U (1)-current, j U (1) μ , is obtained from the matter part of the action alone, The QEoM of the Abelian gauge boson propagator in linear covariant gauges reads It depends on the classically conserved current j μ U (1) (x) of Eq. (18) and holds for renormalized as well as for bare fields. The last, longitudinal, term on the r.h.s. in Eq. (19) arises from the linear covariant gauge fixing in Eq. (12) and is BRST-exact. It has no physically observable effects and may be included in the definition of the currentj In close analogy to the non-Abelian case discussed below one can use Eq. (14) to write, s∂ μc = ss A μ .
The first term on the r.h.s. of Eq. (19) is transverse due to the antisymmetry of the field strength tensor. Fourier transformation (denoted by . . . FT ) of Eq. (19) and longitudinal (L μν = p μ p ν / p 2 ) and transverse (T μν = δ μν − L μν ) projection give the identities, Using the equation of motion of the NL field, Eq. (20a) yields the Ward identity for the longitudinal photon propagator, where γ (x) = γ (φ(x), . . . ) is the local function of the fields in the BRST-exact term of Eq. (12). The correlator of the gauge field with the divergence of the field strength tensor is transverse and is described by a Lorentz invariant function f ( p 2 ), which for p 2 → 0 determines the phase of the model. f (0) = 0 implies a pole at p 2 = 0 due to a massless transverse vector boson in the correlator A model with f (0) = 0 thus has a massless photon and describes a Coulomb phase. The Abelian nature of the field strength is not essential for inferring a massless pole in Eq. (23). One merely exploits the Poincaré invariance and the antisymmetry of the curvature F μν .
In the Abelian case, f ( p 2 ) determines the transverse part of the vector boson propagator, and the photon is massless only if f (0) > 0. Insertion of the definition, Eq. (22), into Eq. (20b) shows that f ( p 2 ) also completely determines the transverse current matrix element, If the current saturates Eq. (20b) in the infrared, the vector boson is short ranged, and f (0) = 0. These relations hold for any Abelian gauge theory in linear covariant gauges. f (0) = 0 characterizes a model describing a Coulomb phase with a massless photon. We next examine Abelian gauge theories in the Coulomb and Higgs phase in more detail. In Sect. 3.3 we consider an Abelian gauge theory that confines color charge and investigate possible signatures of this phase.

The Coulomb phase
The Coulomb phase is characterized by an unbroken Abelian gauge symmetry and the failure of Eq. (26), i.e., failure of the current contribution to saturate Eq. (20b) in the infrared. Since the Abelian gauge symmetry is unbroken, the correlation function A σ (y) j ν U (1) (x) is transverse in any covariant gauge. Equation (20a) and the Ward identity Eq. (21) have the form The elementary quartet is free and massless:
It is interesting to note that in the canonical formalism f (0) = 0 implies that the electromagnetic charge operator, is not well defined. Up to terms proportional to the photon equation of motion, this charge is equivalent to Due to the antisymmetry of the field strength tensor, F νμ = −F μν , the current −∂ ν F νμ (x) and corresponding charge G are themselves conserved. Furthermore, the equal time commutator of G with any local physical operator (x) ∈ P vanishes, because causality requires operators with spatial separation to commute. One thus has is the canonical conjugate of the ghost operator at vanishing momentum and N = {Q BRST , C}. In deriving Eq. (33) one uses causality, the Jacobi identity and that N is BRST exact. All local physical operators (x) ∈ P thus are uncharged and physical operators creating charged particles like the electron necessarily are not local. (NB: To compare with non-Abelian gauge theories, note that Q BRST and C may be replaced byQ BRST andC = d 3 xπc(x) in Eq. (33).) One can construct non-local charged physical states in QED because infrared photon states are almost degenerate with the ground state [71][72][73][74]. The massless vector boson of the Coulomb phase thus prevents one from concluding from Eq. (33) that all physical states are uncharged.

The Abelian Higgs phase
A "spontaneously broken" Abelian gauge theory in the Higgs phase satisfies Eq. (19) differently. From the general discussion one expects that f (0) = 0, the vector boson is massive and the current saturates Eq. (20b) at low momenta, i.e., Eq. (26) holds. In the Higgs phase one also expects (unphysical) massless excitations. We explicitly verify this scenario in QED with charged scalar matter whose self-interactions are described by a Higgs potential with quartic coupling λ and a negative quadratic term proportional to −4λv 2 , In the Higgs phase with v > 0 we parameterize the fields by . The treelevel photon mass is m = gv and that of the Higgs field φ + is m 2 H = 8λv 2 . φ − is massless and couples to the longitudinal photon. The model is invariant under local gauge transfor- Replacing θ(x) by the anticommuting ghost field one arrives at the BRST variations A convenient gauge that eliminates the bilinear coupling of φ − to the longitudinal photon is given by the BRST-exact linear covariant gauge-fixing term [75], where in the second expression the NL field has been shifted: Note that the BRST-exact term L 't Hooft GF of Eq. (37) explicitly breaks not only local but also global U (1)-symmetry.
The QEoM of the photon propagator is given by Eq. (19) with the gauge-invariant and classically conserved matter current The current is BRST invariant, and its divergence is unphysical because the global gauge invariance of the model is broken by BRST exact terms only. To leading order in the loop expansion one has In fact, the divergence ∂ μ j μ U (1) is BRST exact up to equations of motion. In leading approximation where the tree-level QEoM of φ − and of the NL field b was used to obtain the intermediate expressions.
In the broken phase, the current contribution to Eq. (20a) does not vanish and in fact saturates it at low momenta. This is the signature of a "spontaneously broken" gauge theory. Since the divergence of the current is BRST exact up to equations of motion, it does not create physically observable Goldstone bosons and the gauge theory is in a Higgs phase. In 't Hooft gauges the mass of the unphysical scalar particle created by j μ U (1) depends on the gauge parameter ξ and vanishes for ξ = 0 only.
With γ (x) = ξ mφ − (x), the Ward identity of Eq. (21) to leading order asserts, Note that φ − does not contribute to the Ward identity at tree level. This cancelation of mixing terms is a feature of 't Hooft gauges. However, the current matrix element in Eq. (41) saturates Eq. (20a) for p 2 → 0 whereas it vanishes in the Coulomb phase in this limit. This is not a gauge artifact and for ξ = 0 does not depend on the gauge parameter.
Equation (41) gives the tree-level longitudinal propagator in the Higgs phase: which may be directly verified from the quadratic terms of the action Eq. (37). In the last expression of Eq. (42), the ξ -independent term at p 2 = 0 arises from the current matrix element. The second term is the ξ -dependent negative-norm It is one of the correlators of the elementary BRST quartet and for ξ m 2 = 0 vanishes as p 2 → 0, leaving the current to saturate the Eq. (20a) in the Higgs phase.
In the Higgs phase the tree-level approximation to the function f ( p 2 ) defined by Eq. (22) is Since f (0) = 0 the transverse vector boson is short ranged in this phase, The current of Eq. (38) thus saturates Eq. (20b) at low momenta, for any value of the gauge parameter ξ and Eq. (26) holds in the Higgs phase. Note that the physical Higgs particle and vector boson in this model are not charged. These examples illustrate (at tree level) the characteristics that distinguish the unbroken Coulomb and "spontaneously broken" Higgs phases of Abelian gauge theories. If the current contribution saturates the transverse and the longitudinal QEoM of the photon propagator at low momenta, the model describes a "spontaneously broken" Higgs phase. If the current contribution fails to saturate the transverse equation at low momenta, the Abelian gauge theory describes a Coulomb phase with a massless vector particle. The (conserved) transverse part of the Abelian current in our examples is BRST invariant and does not include BRST exact terms. At vanishing momentum it apparently creates only physical particles. This will change when we consider Abelian gauge theories that confine color charge in Sect. 3.3.
First however, let us revisit non-Abelian gauge theories in Linear Covariant Gauges (LCG) for which Kugo and Ojima originally formulated their confinement criterion.

The Kugo-Ojima confinement criterion in linear covariant gauges (LCG)
While the photon is massless and atoms are readily ionized, gluons are only of short range and no hadron has been colorionized. This confinement of color charge is one of the most prominent features of unbroken non-Abelian gauge theories. One expects the confinement of color and the absence of massless gluons to manifest themselves in solutions to the QEoM of the vector boson propagator. We here give a short review of Kugo and Ojima's analysis [1] of unbroken non-Abelian gauge theories in the linear covariant gauge (LCG). Yang-Mills theory in LCG is defined by the Lagrangian The nilpotent BRST transformation in the non-Abelian case is given by and is readily extended to covariantly coupled matter. The Lagrangian (45) is also invariant under an anti-BRST transformation generated bys, It is in addition invariant under global color rotations and preserves ghost number. The BRST and anti-BRST transformations are nilpotent and anticommute, s 2 =s 2 = {s,s} = 0. As in the Abelian case, one defines a BRST charge Q BRST and anti-BRST chargeQ BRST analogous to Eq. (15). The space P of physical operators then is the sector of vanishing ghost number of the BRST cohomology as in Eq. (17). Unphysical states of indefinite norm are associated with BRST quartets. To all orders in perturbation theory these do not contribute to the physical scattering matrix and cannot be created from physical states by physical operators [1,68]. The longitudinal gauge field, ghost c, antighostc and NL field b again form the elementary BRST quartet [1,70]. Contrary to the Abelian case, the transverse gluon of a non-Abelian gauge theory is part a non-perturbative BRST quartet [76][77][78] and not physical.
In the canonical formulation, the global color symmetry of this theory implies the conserved Noether currents, Up to the gluonic QEoM, is equivalent to j LCG a μ and the color charges G a and N a in defined by must vanish at p 2 = 0, implying the absence of a massless vector boson in the adjoint representation of the group. (ii) The function u LCG , defined by must assume the value, u( p 2 → 0) = −1, in the infrared limit.
As in the Abelian case, condition (52) with f LCG (0) = 0, also holds in a non-Abelian Higgs phase, in which the vector bosons is massive. The infrared behavior of u LCG ( p 2 ) thus distinguishes between the Higgs and confinement phase in LCG. In terms QEoM of the gauge boson propagator this distinction may be reformulated as: if the Kugo-Ojima criterion is fulfilled, that is, if the QEoM of the vector boson propagator, at vanishing momentum is saturated by unphysical states only. By contrast, for f LCG (0) = 0 and u LCG (0) = −1 the theory may describe a Higgs phase in which physical states contribute to the saturation of the transverse part of the current matrix element A a σ (y) j LCG b μ )(x) at vanishing momentum. Since one cannot be sure that physical states contribute when u LCG (0) = −1, the criterion of Kugo and Ojima of Eq. (54) is a sufficient criterion for confinement [2].
The longitudinal part of the correlation function in Eq. (53) is uniquely determined by the equation of motion of the ghost field. It saturates the longitudinal part of Eq. (56) for all momenta and implies that the current matrix element in LCG is transverse.
If the Kugo-Ojima criterion is fulfilled, the matrix element of the BRST-exact part of the (generalized) color current sat-urates the transverse gluonic QEoM of LCG at low momentum. Assuming the BRST symmetry is unbroken, Kugo and Ojima proved that the physical sector in this case is colorless [1,2,76]. In a Higgs phase, the physical, albeit massive, vector boson would contribute to the transverse part of the gluonic QEoM at vanishing momentum and it would not be saturated by unphysical states only.
Since color confinement is not compatible with physical states in the adjoint color representation, an alternative criterion for the confining phase, but one that includes the essential information, is the saturation of the gluonic QEoM by unphysical degrees of freedom in the infrared. We now investigate whether this proposal can be applied to a wider class of gauges.

Saturation and confinement in generalized linear covariant gauges (GLCG)
Let us therefore next examine saturation of the gluonic QEoM in Generalized Linear Covariant gauges (GLCG) given by the Lagrangian [79,80] with α +ᾱ = 1. The Lagrangian of Eq. (57) interpolates between LCG (α = 0), its Faddeev-Popov conjugate (α = 1), and the ghost-antighost symmetric gauge at α =ᾱ = 1 2 . The generalized Kugo-Ojima confinement scenario for this Lagrangian is discussed in [81]. For any value of the gauge parameters α and ξ , L GLCG is globally color symmetric and invariant under the nilpotent BRST and anti-BRST transformations, The ghost-antighost symmetric gauge α =ᾱ = 1 2 possesses an additional continuous global SL(2, R) symmetry generated by ∓ and the ghost number , One verifies that the BRST and anti-BRST charges Q BRST andQ BRST , which generate the transformations of Eqs. (58) and (59), anticommute [79], and that the graded algebra of {Q BRST ,Q BRST , + , − , } closes in ghost-antighost symmetric gauges. The space of physical operators P is defined as in Eq. (17). Due to the invariance of Eq. (57) under global color transformations, the Noether currents, are conserved. These currents depend explicitly on the gauge parameter α. They are part of the gluonic QEoM, which here takes the form, The last term in Eq. (62) again involves only unphysical excitations and is of the same form as in the LCG studied above, The equation of motion of the ghost by itself does not suffice to determine the longitudinal part of Eq. (63). Instead one has Since A a μ = 0 and the BRST transformation is nilpotent, Eq. (64) determines the longitudinal part of the correlation function in Eq. (63). As in LCG, unphysical degrees of freedom saturate the longitudinal part of the gluonic QEoM in Eq. (62), and the current matrix element is transverse.
The form factor, f GLCG ( p 2 ), is defined as in Eq. (52), and the same discussion as in Sect. 3.1 applies. The theory describes a Coulomb phase with a massless vector boson only if f GLCG (0) = 0. The transverse part of the correlator in Eq. (63) defines a form factor u GLCG ( p 2 ) whose value in the infrared can be used as a criterion in these gauges. The transverse gluonic QEoM is saturated by unphysical degrees of freedom in the infrared and the phase is confining if u GLCG (0) = −1. One thus formally has the same confinement criterion as in LCG [81]. However, although the unphysical correlation functions in Eqs. (63) and (53) formally look similar, unphysical correlations differ, and in general u GLCG ( p 2 ) = u LCG ( p 2 ) if α = 0. The confinement criterion stated above asserts that these functions in the confining phase may coincide at p 2 = 0 in any gauge parametrized by (α, ξ ). We next turn to covariant gauges that explicitly break global color invariance to the Cartan subgroup.

Saturation and confinement in maximal Abelian gauges (MAG)
An equivariant BRST construction allows one to partially localize a non-Abelian gauge theory to an equivalent Abelian model with the same gauge-invariant correlators. This partial localization is possible on the lattice [41,82,83] as well as in the continuum [44,84,85] and defines the continuum theory as the critical limit of a lattice model with the same global symmetries. The equivariant construction may be viewed as a partial gauge fixing that leaves the Abelian Cartan subgroup free, hence Maximal Abelian Gauge. The resulting Abelian gauge theory has the same gauge-invariant physical correlation functions as the non-Abelian model and the residual Abelian gauge symmetry of MAG can be dealt with as in any Abelian gauge theory. An SU(N ) gauge theory in MAG presents itself as an Abelian (U (1)) N −1 gauge theory that is asymptotically free but retains typical Abelian Ward identities. In the following we investigate how a confining phase may manifest itself in this Abelian gauge theory. In MAG one discriminates between the Cartan subgroup and the coset space. For an SU(N ) gauge theory the mutually commuting (Hermitian) generators of the Cartan subgroup will be denoted by {T i ; [T i , T j ] = 0 ∀i, j = 1, . . . , N − 1}, whereas the remaining N (N − 1) coset generators carry Latin indices from the beginning of the alphabet {T a ; a = 1, . . . , N (N − 1)}. The non-Abelian SU(N ) connection is decomposed as A μ = A i μ T i + B a μ T a and the field strength tensor may similarly be decomposed into Cartan and coset components, Here and in the following the covariant derivative with respect to the Cartan gluons in the adjoint representation is denoted by D ab μ ; see Appendix A. The SU(N ) Yang-Mills Lagrangian in the components of Eq. (65) reads Although a more general discussion is possible, we for simplicity consider the gauge group SU(2) in the following. It illustrates our main points and connects to our considerations in Sect. 2. The Cartan subalgebra in this case is one dimensional and we suppress its index. The coset space is two dimensional with components a = 1, 2.
The distinction between Abelian and coset degrees of freedom is accomplished by the "gauge-fixing" part of the MAG Lagrangian, L MAG GF is obtained using the equivariant BRST and anti-BRST transformations [41,44], Analogous to Eq. (15) one can define the equivariant BRST and anti-BRST charges, Q ε B RST andQ ε B RST . The Lagrangian (67) and the algebra of Q ε B RST andQ ε B RST are ghost-antighost symmetric. In addition to the ghost number , one thus has the complete SL(2, R) algebra of charges typical of ghost-antighost symmetric gauges. However, in MAG the ∓ of Eq. (60) involve only two coset (anti-)ghosts rather than the three of an SU(2) gauge theory in ghost-antighost symmetric GLCG.
The equivariant BRST transformations (68) with gauge-fixing parameter ξ and (uncharged) auxiliary field b. b is taken to be invariant under s ε ands ε , s ε b = s ε b = 0. We thus consider the SU(2)-Yang-Mills theory in MAG specified by the action S MAG = d 4 xL MAG (x) with, The where generates a local U (1) gauge transformation at x. Defining the set W of U (1)-invariant operators, we show in Appendix B that On the set W of U (1) invariant operators of the equivalent Abelian gauge theory one thus recovers s ε ands ε as nilpotent BRST symmetries. One then can define the set of physical operators P ⊆ W of the underlying non-Abelian SU (2) gauge theory by the equivariant cohomology, The conserved U (1) current of the Cartan subalgebra of SU (2) in MAG is and can be rewritten in the form, In fact, each term in Eq. (78b) separately is an element of W, but i [s ε ,s ε ] A μ does not create physical states. The QEoM of the Cartan gluon propagator depends on the conserved Abelian Noether current of Eq. (78b) as in the Abelian gauge theory studied in Sect. 2, As for an unbroken Abelian gauge theory, the last term of Eq. (79) saturates the longitudinal part of the gluonic QEoM because of the Abelian Ward identity of Eq. (73a), The first correlator in Eq. (79) is transverse due to the antisymmetry of f μν and the current matrix element thus is transverse as well, The functions f MAG , u MAG , h MAG , and MAG are defined from the correlators and the transverse part of Eq. (79) yields the constraint As in LCG and GLCG, the transverse gluonic QEoM is saturated by unphysical degrees of freedom and the Cartan color charge of physical states vanishes when the Kugo-Ojima-like criterion holds since it implies that h MAG (0) = 0. The conditions (84) guarantee saturation of the gluonic QEoM in the infrared by unphysical degrees of freedom in MAG and implies that physical states are colorless. Equation (81) together with Eq. (84) imply that this scenario can only be realized in MAG if unphysical degrees of freedom created by the conserved Abelian current j MAG μ saturate the QEoM of the Abelian propagator at low momenta. From the point of view of the Abelian gauge theory, saturation of the transverse equation at low momenta in confinement and Higgs phases thus is similar. In MAG the only difference is that whereas some physical degrees of freedom contribute to the matrix element of the current at vanishing momentum in the Higgs phase, only unphysical states contribute in the confinement phase. The current saturates the QEoM of the Abelian propagator at low energies in the Higgs phase described by Eq. (37) as well as in the confinement phase of the SU(2) gauge theory in MAG defined by Eq. (69). This supports the idea that the phases and the two Abelian models describing them are dual [8].
In this context it is interesting to consider the condition f MAG (0) = 0 more closely. The analogous condition implies a massive physical vector boson in the Abelian theory described by Eq. (34). The Cartan gluon on the other hand is not a physical asymptotic state in the confinement phase and it has been conjectured [49] that the Abelian propagator in this case may even be enhanced in the infrared. That this scenario can be reconciled with the criterion of Eq. (84) rests on the definition Eq. (65a) of the Abelian field strength tensor f μν . Equation (82a) implies Although a massive Abelian vector boson allows one to fulfill f MAG (0) = 0, the last correlator of Eq. (85) prohibits one from asserting that the diagonal gluon propagator has to be suppressed at low momenta. Introducing the function α MAG ( p 2 ), Equation (85) states that If the Cartan gluon correlator is infrared enhanced, Eq. (87) determines only the infrared singular behavior of α MAG ( p 2 ) when f MAG (0) = 0. To gain some more information about these functions, we define a U (1)-invariant transverse field strength, It is not an invariant of the equivariant BRST (or anti-BRST) and, in contrast to the U (1) gauge theory considered in Sect. 2, is not a physical operator of the SU(2) gauge theory. The function u MAG ( p 2 ) defined in Eq. (82b) also describes the correlation functions, where Eq. (76)   also need not be related to MAG ( p 2 ) defined by Eq. (82b). 1 We next study signatures of the confining phase in a gauge that is not covariant.

Saturation and confinement in non-Abelian coulomb gauge (CG)
Coulomb gauge breaks manifest Lorentz covariance by treating timelike and spacelike gluons differently. We here study to what extent the Kugo-Ojima criterion depends on a manifestly Lorentz-invariant gauge condition. Coulomb gauge is described by the Lagrangian, where Latin indices denote spatial components of a Lorentzvector, i, j, . . . = 1, 2, 3. The BRST transformations are the same as in Eq. (46), 1 Note that the definition of the function u MAG ( p 2 ) in MAG apparently differs by a factor of (−2) from that of GLCG given by Eq. (63), since (D μ c) 3 = ∂ μ c 3 + B μ ×c 3 formally differs from B μ ×c in Eq. (90) by a longitudinal contribution only. However, Eqs. (63) and (90) are gauge dependent correlation functions that are not required to coincide in two different gauges. In addition, as LCG and MAG are not analytically connected, [53], no quantitative relation should be expected between the functions u LCG and u MAG .
The BRST charge in Coulomb gauge can be written in terms of Gauss's law [86] The anti-BRST transformations and the corresponding charge may be defined analogously and the set of physical operators again is given by the BRST cohomology of Eq. (17). Coulomb gauge manifestly preserves global color symmetry and the color currents are conserved. The absence of manifest Lorentz invariance in Coulomb gauge implies two distinct gluonic QEoMs. The QEoM of the time component reads Since all physical states satisfy Gauss's Law in Coulomb gauge, this equation of motion is saturated by unphysical states only, whether the model confines or not. To see that all states that contribute to Eq. (95) are unphysical note that physical states | phys are created by physical operators defined in Eq. (17). They have vanishing ghost number and are annihilated by the "Gauss-BRST" charge of Eq. (93), The ghost field c does not annihilate | phys , since its only effect is to create a ghost. Equation (96) thus has to be ensured by gluonic contributions only, and one gets back Gauss's law as the subsidiary condition, Any non-vanishing contribution to Eq. (95) thus must be due to unphysical |ψ ∈ {| phys }. In Appendix C we give an explicit calculation of the r.h.s. of Eq. (95), and we relate the propagator of the temporal gluon to the Faddeev-Popov operator to show that it is saturated by instantaneous contributions only. The discussion of the spatial components of the gluonic QEoM is very similar to that in LCG. The equation of motion for the spatial part of the gluon propagator is given by The first matrix element necessarily is spatially transverse in Coulomb gauge with ∂ i A i = 0. The Faddeev-Popov operator of Coulomb gauge is instantaneous, and the contribution of the last term in Eq. (98) therefore is instantaneous, Its Fourier transform depends on spatial momenta only. The QEoM of the ghost gives the longitudinal part of the correlation function where t i j and l i j are the longitudinal and transverse spatial projectors. The confinement criterion in Coulomb gauge reads lim p 2 →0 u C (p 2 ) = −1 and lim where the function f C ( p 0 , p) is defined by The conditions of Eq. (103) ensure that the spatial gluonic QEoM is saturated by unphysical degrees of freedom in Coulomb gauge. It is interesting that the correlation function in Eq. (101) is related to the horizon function of minimal Coulomb gauge in a finite quantization volume V [56][57][58][59], with In minimal Coulomb gauge, the configuration space is constrained to the first Gribov region by imposing the constraint which in fact is equivalent to the condition u C (0) = −1 for color confinement of Eq. (106). A similar relation between the Kugo-Ojima criterion and the horizon condition also holds in minimal Landau gauge [58,65], to which we now turn.

Saturation and confinement in the covariant Gribov-Zwanziger(GZ) theory
To avoid summation over gauge equivalent configurations, the GZ approach seeks to dynamically restrict the path integral to the first Gribov region [55][56][57][58][59]. The restriction leads to a horizon condition similar to Eq. (107) and can be implemented in a local renormalizable field theory with additional auxiliary ghosts. It was shown [59,60] that the GZ Lagrangian differs from L YM by BRST-exact terms only. An infrared analysis of the GZ action reveals that its scaling solution coincides exactly with the solution calculated from the Faddeev-Popov action for Landau gauge [62,63]: the ghost propagator is infrared enhanced and the gluon propagator infrared suppressed, the respective infrared exponents are identical. This corroborates the argument that for functional equations it suffices to take into account the appropriate boundary conditions, and no explicit restriction in the pathintegral measure is required. However, the horizon condition implies that the BRST symmetry of the GZ action is spontaneously broken. In this last section we want to investigate how the gluonic QEoM in minimal Landau gauge is saturated in the infrared, even though the spontaneously broken BRST symmetry prohibits a definition of physical operators as in the foregoing sections. The auxiliary ghost-fields, φ a μb ,φ a μb , ω a μb , andω a μb are vector fields with two color indices that transform under the adjoint representation of the global color group (in SU(3) they are reducible 8 ⊗ 8 = 1 ⊕ 8 ⊕ 27 ⊕ 8 ⊕ 10 ⊕ 10 color tensors), for any a μb ∈ {φ a μb ,φ a μb , ω a μb ,ω a μb }. While φ a μb andφ a μb are bosonic, ω a μb andω a μb are fermionic. The auxiliary ghosts form a BRST quartet, sω a μb =φ a μb , sφ a μb = 0.
Including this auxiliary quartet in the BRST transformations of Eq. (46), the BRST-exact part of the GZ Lagrangian is The restriction of configuration space to the first Gribov region can be interpreted as a spontaneous breakdown of this BRST symmetry. As in any instance of a spontaneously broken symmetry it is advantageous to express the Lagrangian in terms of fluctuations about the symmetry breaking ground state. In the GZ framework this amounts to a shift of the fields by where tr a μ = g f abc b μc . 2 This change of variables in Eq. (111) gives the GZ Lagrangian of minimal Landau gauge [59,60], Although the shift (112) The Gribov parameter γ is found by demanding that the model is quantized about an extremum of the quantum effective action , δ δγ = 0.
The inhomogeneous term of Eq. (114d) causes the BRST symmetry of the local Lagrangian to be spontaneously broken for any extremum of the quantum action with non-vanishing γ . It is perhaps worth noting that for γ = 0 the Poincaré generators do not commute with the BRST charge even though Poincaré invariance is not spontaneously broken.
To proceed with our investigation of confinement criteria in various gauges, we consider the gluonic QEoM implied by the Lagrangian where we suppressed all indices that are summed. Contractions with structure constants in the "covariant" and "contravariant" color indices are denoted by and Equation (116) includes the global color current j LCG a μ of LCG given in Eq. (48). However, j LCG a μ is not the conserved color current of the GZ action since the auxiliary fields transform according to Eq. (108). The corresponding conserved color current is Using Eq. (119) the QEoM of Eq. (116) may be rewritten with The gluonic QEoM of the GZ action therefore has the now already familiar form Color transport is short ranged and the current matrix element does not contribute in the infrared if the functions f GZ ( p 2 ) and u GZ ( p 2 ) defined by satisfy the criteria However, in this case of a spontaneously broken BRST symmetry it is not entirely clear that holds, which would imply that only (unphysical) quartet states contribute to the matrix element of Eq. (123b). Due to the equations of motion of the antighostc and of the NL field b a the longitudinal part of Eq. (125) is satisfied. Although a proof is lacking, it therefore is at least plausible that the transverse part of Eq. (125) also holds. The GZ action incorporates non-perturbative features and in fact satisfies the criteria (124) already at tree level. Expanding the gluonic QEoM (120) to tree level yields We again have that the longitudinal part of the gluon propagator is saturated by the NL field as in the foregoing investigations. The transverse part of Eq. (126) is satisfied by the tree-level propagators, given for example in [59,87] (with and, The GZ gluon propagator vanishes in the infrared, and f GZ (0) = 0. In addition the last term in Eq. (126), derived entirely from the BRST-exact term in Eq. (120), saturates the transverse part of the gluonic QEoM at tree level for vanishing momenta, Quite strikingly, both criteria of Eq. (124) for a confining phase are thus satisfied by the GZ Lagrangian already at tree level. Perturbative calculations to one-and two-loop order in 3 [88] and 4 [87,89,90] Euclidean dimensions as well as a non-perturbative infrared analysis [62,63] of the GZ action show that in the infrared the gluon propagator remains suppressed and the ghost propagator diverges more strongly than a massless pole [58,62,63]. This infrared behavior agrees with the original Kugo-Ojima scenario [2].

Conclusion
In summary, we have formulated as confinement criterion that the gluonic QEoM be saturated by unphysical states in the infrared. In the Higgs and Coulomb phases this is not the case. These conditions thus are sufficient for distinguishing a color-confining phase from a Higgs and a Coulomb phase in linear covariant (LCG), generalized linear covariant(GLCG), maximal Abelian (MAG), and Coulomb (CG) gauges. Although the details depend somewhat on the chosen gauge, a universal qualitative criterion emerges in theories with an unbroken BRST or equivariant BRST symmetry that distinguishes between physical and unphysical states.
In the considered gauges the QEoM of the gauge boson propagator is of the form where the local currentj a μ (x) differs from the canonical Noether current j b μ (x) of the model by a BRST-exact contribution only, In models with unbroken BRST symmetry,j a μ thus is physically equivalent to the conserved Noether current j a μ . The criteria distinguish the phases depending upon which term on the right-hand side of Eq. (130) saturates the unity on the left in the infrared.
For the models we considered, the generalized color currentj a μ (x) is given bỹ In all gauges with unbroken BRST or equivariant BRST symmetryj μ (x) is physically equivalent to the conserved current j μ (x) because the additional terms either are BRST-exact or vanish on the physical Hilbert space due to subsidiary conditions.
In the Coulomb phase, the first matrix element in Eq. FT , does not vanish in the infrared. f (0) = 0 implies the existence of a massless vector boson.
In the Higgs and color-confining phases on the other hand, f (0) = 0 and the current matrix element saturates Eq. (130) in the infrared limit p 2 → 0.
Since no gauge-invariant order parameter discriminates between the Higgs and color-confining phases [3], the question arises whether one can distinguish between them at all. All physical states are colorless in both phases [8,91]. However, the states contributing to A a σ (y)j b μ (x) FT at low momentum differ in these two phases. In the confining phase the infrared limit of the transverse part of this current matrix element is entirely saturated by unphysical states. One in particular can be sure that the phase is confining if a BRSTexact part ofj a μ (x) saturates the current matrix element. In the Higgs phase physical states contribute to this matrix element.
We found that such criteria distinguishing the confining from the Higgs and Coulomb phases also exist for Abelian gauge theories. In the Abelian Coulomb phase, discussed in Sect. 2, the current does not saturate the QEoM of the photon propagator at low momenta and the photon is massless. At tree level in the Abelian Higgs model of Sect. 2.2, f (0) = 0, and the massive physical vector boson saturates the transverse part of the current matrix element. The only BRST-exact contribution to the generalized current in this case is longitudinal. The SU(2) gauge theory in MAG, discussed in Sect. 3.3, can be viewed as an Abelian U (1) gauge theory in a confining phase. If unphysical states created by the [s ε ,s ε ]A μ -part of the Noether current saturate the current matrix element at vanishing momentum, the theory describes a confining phase. Only the mutually commuting color charges of the Cartan subgoup are conserved in non-Abelian gauge theories in MAG and an unphysical part of the corresponding Abelian Noether current can saturate the gluonic QEoM in the infrared.
The saturation in GLCG, considered in Sect. 3.2, resembles that in LCG originally discussed by Kugo and Ojima. In these gauges the model confines color if the BRST-exact term, −is α (D μc ) a , ofj a μ saturates the gluonic QEoM in the infrared.
In the non-Abelian Coulomb gauge studied in Sect. 3.4, only unphysical states that do not satisfy Gauss's Law contribute to the temporal part of the gluonic QEoM in all phases (and at all momenta). The temporal part of the gluonic QEoM thus cannot discriminate between phases. However, the theory is again confining if the spatial part of the gluonic QEoM in Coulomb gauge is saturated by the BRST-exact s(D kc ) a term ofj a k . The corresponding confinement criterion of the Coulomb gauge was in addition found to be identical to the horizon condition of minimal Coulomb gauge.
An equivalence between the Kugo-Ojima confinement criterion and the horizon condition of minimal Landau gauge has also been established in [58,65]. The auxiliary fields also contribute to the conserved Noether currents j GZ a μ of the GZ action (see Eq. (119)), but the gluonic QEoM retains the form of Eq. (130) withj a μ given by Eq. (132f). In this model the BRST-exact part sχ a μ ofj a μ saturates the gluonic QEoM at p 2 = 0 already at tree level. There is no massless vector boson, and the gluon propagator at low momentum is suppressed. The GZ-theory in this sense satisfies all the confinement criteria for gauge theories with BRST symmetry. However, at present it is not known how to define a physical Hilbert space in this model with a spontaneously broken BRST symmetry [59][60][61]92] and one has to prove that the BRST-exact contribution of the generalized current does not create physical states.
The non-Abelian field strength is defined by the relation We use an equivalent sign ≡ between expressions that differ by terms that vanish when the classical equations of motion are satisfied and ≈ if expressions coincide to leading order, usually tree level.

Appendix B: Proof of restored BRST symmetries in MAG
Here we prove that the expectation values of equivariant BRST and anti-BRST variations (given in Eq. (68)) of U (1)invariant operators vanish for an SU(2) gauge theory in MAG, cf. Eq. (76), where Eq. (75) defines the space W of U (1)-invariant operators. With mild restrictions on the topology of spacetime, i.e. the Laplace-operator has to have an inverse, Eq. (73a) for any O ∈ W implies that The variations of the U (1)-gauge field A μ satisfy δ x s ε A μ = δ x B μ ×c = 0 (140) and δ xsε A μ = δ x B μ ×c = 0.

Appendix C: The gluonic QEoM of the A a 0 field
In this appendix we integrate out the A 0 field in Eq. (95) and show that the equation is saturated by instantaneous contributions only. With the action S C given by the Lagrangian Eq. (91), the QEoM for the A 0 -field is We decompose the right-hand side into two terms, so the gluonic QEoM reads δ(x − y) = I 1 (x − y) + I 2 (x − y), where I 1 (x − y)δ ab := − A b 0 (y) (D 2 i A 0 ) a (x) Here we used (D iȦi ) a = A i ×Ȧ i a , which follows from the transverse Coulomb gauge condition ∂ i A a i = 0. To improve readability in the following we suppress color indices and add them only where necessary. We wish to express these expectation values in terms of an integral over the canonical variables and make use of an identity proven in [93]: where ρ is a source for A 0 . To obtain this formula, one introduces the color-electric field E i by an auxiliary integration, after which one integrates out A 0 and the longitudinal part of E i . This takes one from the Faddeev-Popov formula for integrating over d 4 A = d 3 A i d A 0 to an integration over the canonical variables of the Coulomb gauge, A tr b i and E tr b i , which are the three-dimensionally transverse vector potential and chromoelectric field. In the last formula, the Hamiltonian density is given by where and etc., where (Kρ coul )(x) ≡ d 4 y K(x, y)ρ coul (y), and the color-Coulomb kernel is given by The identity Eq. (150), when applied to Eq. (145), gives We next separate the instantaneous and non-instantaneous parts of these expressions. The kernel K (x, y) = K (x, y)δ (x 0 − y 0 ) is instantaneous, so the first term of I 1 is purely instantaneous. The second term of I 1 involves the canonical fields E tr and A tr at time x 0 and the canonical fields at time y 0 . These are the dynamical degrees of freedom so their correlators are non-instantaneous.
Keeping only the instantaneous part in Eq. (152) one gets To separate the instantaneous part in I 2 , we shall expressȦ tr i in terms of the canonical fields A tr i and E tr i . For this purpose we use the fact that the integral of a derivative vanishes,