Codex @codex  0

The $\xi \to 0$ limit of covariant gauge fixing by $-({\partial }_{\mu }{A}^{\mu }{)}^2/(2\xi )$. It imposes the divergence condition sharply and gives a transverse photon propagator. This is distinct from the vector-potential choice called Landau gauge for a uniform magnetic field when computing Landau levels.

Choose a definite matrix convention: $A_{\mu }=A_{\mu }^a{T}_a$, $\omega ={\omega }^a{T}_a$, and $D_{\mu }={\partial }_{\mu }+e{A}_{\mu }$. Matter transforms as ${\psi }^{\prime }=U\psi$, so gauge covariance requires $D_{\mu }^{\prime }=U{D}_{\mu }{U}^{-1}$. Applying both sides to a test field yields the Yang-Mills gauge transformationDefine the Yang-Mills field strength by $[D_{\mu },D_{\nu }]=e{F}_{\mu \nu }$:Conjugating the gauge covariant derivative commutator proves $F_{\mu \nu }^{\prime }=U{F}_{\mu \nu }{U}^{-1}$: unlike the connection, the gauge field strength has no inhomogeneous derivative term. The trace normalization supplies an invariant color metric; its Lie algebra structure constants with all indices lowered are antisymmetric.

For $U=1+e\omega +O({\omega }^2)$ and $U^{-1}=1-e\omega +O({\omega }^2)$, expand the finite transformations:In components this isA convention using the opposite matter transformation or gauge covariant derivative sign changes these signs consistently. The formulas above all follow from the explicitly chosen $D_{\mu }$.

Gauge reduction and the determinant. A direct integral $\int \mathcal{D}A{e}^{iS[A]}$ integrates repeatedly over gauge-equivalent configurations. The gauge group volume is infinite, and the ungauge-fixed quadratic operator has longitudinal zero directions, so it has no ordinary inverse quantum field theory propagator. Gauge-invariant observables require dividing out this redundancy rather than treating it as an ordinary physical fluctuation.

The finite-dimensional change-of-variables identity underlying gauge fixing is precise. Suppose a condition $\chi (x^{\alpha })=f$ has one solution in the group-coordinate patch and its derivative matrix at the solution is nonsingular. ThenFor several isolated solutions it instead counts those solutions. On the slice, the Jacobian determinant can be factored into the reciprocal of the orbit integral of the delta function. An invariant measure and invariant integrand then allow the group volume to factor out. In the functional version, group coordinates carry Haar measure and the corresponding operator determinant replaces the matrix determinant.

For example, choose the Lorenz gauge ${\chi }^a(A)={\partial }^{\mu }{A}_{\mu }^a=f^a$. Its variation is $\delta {\chi }^a=-{\partial }^{\mu }(D_{\mu }^{\mathrm{ad}}\omega {)}^a$. Hence on the gauge slice the Faddeev-Popov operator and Faddeev-Popov determinant areIts extension along each local gauge orbit yields the Faddeev-Popov gauge-orbit identityFor an unoriented finite-dimensional real integral the absolute determinant is required. In a perturbative gauge patch, its sign or phase is fixed and the constant convention is absorbed in normalization. Residual zero modes must be removed by boundary conditions or separate treatment. Multiple intersections, the Gribov ambiguity, prevent interpreting this as a globally unique gauge slice without additional work.

Insert the identity, change variables along the orbit using invariance of the action and measure, and divide out the common gauge group volume. Averaging $f$ with Gaussian weight $\exp [-i\int f^a{f}^a/(2\xi )]$ produces the covariant gauge-fixed functionalThe added term makes the quadratic gauge-field operator invertible after the stated residual-mode treatment.

Ghost representation. For a finite matrix $M$, the Grassmann Gaussian integral over independent anticommuting variables is proportional to its determinant:The constant is independent of $M$ and depends on ordering and factors of $i$. The functional counterpart introduces Faddeev-Popov ghost fields $c^a$ and antighost fields ${\overline{c}}^a$, withIntegration by parts writes its density as $({\partial }^{\mu }{\overline{c}}^a)(D_{\mu }c{)}^a$, displaying a kinetic term and a gauge-field/ghost interaction. These are Grassmann-valued fields with scalar Lorentz transformation law. They represent the orbit Jacobian determinant; their closed loops carry the fermionic minus sign, and they are not physical asymptotic particles. In the abelian case $D^{\mathrm{ad}}=\partial$, the determinant is gauge-field independent and the ghosts decouple. Nonabelian ghosts interact and are necessary for the consistent cancellation of unphysical gauge contributions.

Physical-state reduction. Let $Q$ be the Hermitian BRST charge. Its key algebraic property is BRST nilpotence, $Q^2=0$, and the gauge-fixed dynamics preserves it. Physical states are its BRST cohomology at ghost number zero:Here $|\chi ⟩$ has ghost number $-1$. BRST nilpotence ensures that exact vectors lie in the kernel, so the quotient is defined. Because $Q$ is Hermitian, an exact vector is orthogonal to every closed vector, $⟨\psi |Q\chi ⟩=⟨Q\psi |\chi ⟩=0$, and exact vectors are null, $⟨Q\chi |Q\chi ⟩=⟨\chi |Q^2\chi ⟩=0$. They therefore represent no independent physical state. Commutation of $Q$ with the Hamiltonian preserves these classes under time evolution.

The covariant gauge/ghost space before reduction has an indefinite inner product. This is essential: a Hermitian nilpotent BRST charge requires an indefinite auxiliary space if it is nonzero. On an ordinary positive-definite Hilbert space the same norm identity would force $Q=0$. Positivity of the resulting physical space uses the gauge-theory structure in addition to the abstract BRST nilpotence identities.

For a nonzero on-shell photon momentum, the one-particle conclusion isThe second relation preserves the first because $k^2=0$. In four dimensions the transverse kernel has dimension three, and quotienting by its null longitudinal direction leaves two physical photon polarizations. This is the positive one-particle BRST cohomology of the photon sector.

The Maxwell quadratic action has a null longitudinal direction: its momentum-space kinetic operator annihilates $k_{\mu }$, reflecting gauge invariance under addition of a gradient to the potential. It cannot be inverted on all four vector components before gauge fixing.

Add the covariant gauge density ${\mathcal{L}}_{\mathrm{gf}}=-({\partial }_{\mu }{A}^{\mu }{)}^2/(2\xi )$. Away from the pole, introduce the transverse and longitudinal momentum projectorsWith one index raised they obey $P^{\mathrm{T}}+P^{\mathrm{L}}=I$, $(P^{\mathrm{T}}{)}^2=P^{\mathrm{T}}$, $(P^{\mathrm{L}}{)}^2=P^{\mathrm{L}}$ and $P^{\mathrm{T}}{P}^{\mathrm{L}}=0$. The gauge-fixed kinetic operator is $K=-k^2(P^{\mathrm{T}}+{\xi }^{-1}{P}^{\mathrm{L}})$. Its inverse times $i$ gives the photon propagatorThe pole factors are interpreted with the Feynman boundary prescription; the algebraic projector inversion is performed off the pole. The choice $\xi =1$ is Feynman gauge.

Changing $\xi$ changes only a longitudinal term proportional to $k_{\mu }{k}_{\nu }$. Between conserved currents it vanishes because $k_{\mu }{J}^{\mu }=0$. In full QED amplitudes, the Ward-Takahashi identity supplies the corresponding cancellations when all diagrams at an order are combined. Thus physical on-shell scattering is gauge-parameter independent, whereas off-shell Green functions and the propagator itself can depend on the gauge choice. For the charged scalar, derivative vertices and contact terms must likewise be included for this cancellation.

In quantum electrodynamics, gauge invariance connects local phase freedom, electric-charge conservation and the restrictions on photon interactions. For definiteness takeThe gauge transformationhas $D_{\mu }\psi \mapsto e^{-ie\alpha }{D}_{\mu }\psi$: the derivative of the local phase cancels the added gauge potential. Also $F_{\mu \nu }$ is unchanged because mixed partial derivatives commute. Consequently every term in the Lagrangian density is invariant. This derives the minimal electromagnetic coupling of a Dirac field, including the interaction $-e\bar{\psi }{\gamma }^{\mu }\psi A_{\mu }$. A fermion mass is compatible with the symmetry, whereas a local Proca photon mass $m_{\gamma }^2{A}_{\mu }{A}^{\mu }/2$ is not.

The classical Maxwell equations are ${\partial }_{\mu }{F}^{\mu \nu }=j^{\nu }$, with $j^{\nu }=e\bar{\psi }{\gamma }^{\nu }\psi$. Taking a divergence gives ${\partial }_{\nu }{j}^{\nu }=0$ because $F^{\mu \nu }$ is antisymmetric. The same conservation law follows directly from the charged Dirac equation and its adjoint. For vanishing flux at spatial infinity, $Q=\int d^{3}x{j}^0$ is conserved. Gauss law identifies this charge with the electric flux through a surrounding surface. Particle-antiparticle creation is allowed, but its net electric charge is zero.

Local gauge invariance is also a redundancy in the potentials, rather than an extra propagating degree of freedom. The canonical momentum of $A_0$ vanishes, and its field equation enforces the Gauss law constraint in gauge theory, $\nabla \cdot \mathbf{E}=j^0$. After imposing the constraint and quotienting gauge freedom, a massless photon has two transverse physical polarizations. Gauge transformations vanishing at the boundary identify equivalent descriptions; constant phase transformations at infinity can act on charged states and generate their conserved charge. This distinction explains why gauge redundancy and electric charge are compatible.

In perturbative QED, the quadratic photon operator cannot be inverted before gauge fixing. Adding $-({\partial }_{\mu }{A}^{\mu }{)}^2/(2\xi )$ gives the covariant gauge propagatorThis propagator contains longitudinal and timelike components, but these are not external physical photon states. A Gupta-Bleuler null-state quotient, or the equivalent BRST construction, removes them. In linear Lorenz gauge, $\delta ({\partial }_{\mu }{A}^{\mu })=\square \alpha$, so the Faddeev-Popov determinant is independent of $A_{\mu }$ and the Faddeev-Popov ghost fields decouple. This last statement concerns the linear gauge; Abelian gauge groups alone do not guarantee decoupling in a nonlinear gauge condition.

At the quantum level, the vector symmetry gives a Ward identity. One way to obtain it is to change variables by a localized phase in a time-ordered correlation function containing $\psi (y)$ and $\bar{\psi }(z)$. The action variation becomes a current divergence; the two field insertions contribute opposite contact terms, since their charges are opposite. Fourier transformation turns the divergence into the photon momentum, and amputating the external fermion propagators turns the two contact terms into the difference of their inverse propagators. With the charge removed from the proper vertex and the overall propagator factor $i$ removed from its inverse, the result is the proper-vertex Ward-Takahashi identity in QEDAt tree level this is simply $k_{\mu }{\gamma }^{\mu }=(/p+/k-m)-(/p-m)$. Between on-shell external spinors, the inverse propagators vanish, so the full physical scattering amplitude satisfies $k_{\mu }{\mathcal{M}}^{\mu }=0$. Therefore replacing a photon polarization vector by ${\epsilon }_{\mu }+c{k}_{\mu }$ leaves the amplitude unchanged. The identity applies to the complete required sum of diagrams, not necessarily each diagram separately. It also ensures cancellation of the gauge parameter from physical observables when all terms at the chosen perturbative order are included.

For the photon two-point function, the Ward identity implies transverse photon vacuum polarization:Thus ultraviolet subtraction renormalizes the gauge-invariant kinetic term, rather than requiring a local photon-mass counterterm. Perturbative unbroken QED retains a massless photon. Transversality by itself would not exclude a nonlocal massless pole in a different theory; the assertion here includes the perturbative QED setting.

Taking $k\to 0$ in the proper-vertex identity gives ${\Gamma }^{\mu }(p,p)=\partial S^{-1}(p)/\partial p_{\mu }$. Matching the divergent coefficients of the derivative and vertex terms yields the Ward identity for QED renormalization constants, $Z_1=Z_2$, in a gauge-preserving prescription. If ${\psi }_0=\sqrt{Z_2}\psi$ and $A_0=\sqrt{Z_3}A$, comparison of the interaction coefficients givesHence charge renormalization is fixed by photon wave-function renormalization. Gauge invariance constrains counterterms and relates renormalizations; it does not make the electric charge independent of scale.

The Soft photon theorem gives another physical manifestation. Attaching a photon of small momentum $k$ to the external charged legs gives the leading factorwhere $Q_i$ is charge in units of $e$ and ${\eta }_i=+1$ for outgoing and $-1$ for incoming legs. Replacing $\epsilon$ by $k$ makes this factor proportional to ${\sum }_i{\eta }_i{Q}_i$, which vanishes by electric charge conservation. Massless photons also cause infrared divergences; physical inclusive measurements combine unresolved real emission with virtual corrections, rather than treating either contribution alone as an observable.

Finally, the vector gauge symmetry of a Dirac fermion in QED is anomaly-free, so it can be maintained by regularization and renormalization. This must be distinguished from the axial current: its divergence has explicit mass terms, as in Question 3, and can additionally have a quantum chiral anomaly. Conservation of the gauge current does not imply axial-current conservation. The resulting picture is that gauge invariance simultaneously organizes the physical photon states, enforces electric-charge conservation and ties together amplitudes and ultraviolet subtractions.

A Non-Abelian gauge theory describes fields with a local redundancy, not an independent physical degree of freedom for each component of its gauge field. Take a compact gauge group with Hermitian generators $T^a$ and Lie algebra structure constants defined by $[T^a,T^b]=i{f}^{abc}{T}^c$. In a mostly-plus metric signature, the pure Yang-Mills theory hasWith the matter gauge covariant derivative $D_{\mu }={\partial }_{\mu }-ig{A}_{\mu }$, a finite gauge transformation acts as $A_{\mu }\mapsto U{A}_{\mu }{U}^{-1}-(i/g)({\partial }_{\mu }U)U^{-1}$. Expanding $U=e^{ig{\omega }^a{T}^a}$ givesThe Yang-Mills action is constant along each gauge orbit. Consequently a naive functional integral $\int \mathcal{D}A{e}^{iS[A]}$ repeatedly integrates the same physical configuration and contains the infinite formal volume of the local gauge group. At quadratic order the same problem appears as noninvertibility: the free kinetic operator annihilates $A_{\mu }={\partial }_{\mu }\omega$, so no unique gauge-boson propagator exists. In canonical language the Gauss law constraint in gauge theory and the absence of an independent momentum for the time component show that the field components cannot all be independent oscillators. In covariant quantization, time-like and longitudinal polarizations must not be counted as additional physical states.

The Faddeev-Popov determinant implements a local quotient by the gauge orbit. Choose $G^a[A]={\partial }^{\mu }{A}_{\mu }^a$ and first impose $G^a[A]=f^a(x)$ for a specified function $f$. This is the Lorenz gauge condition when $f=0$. Its infinitesimal change isThe direct Faddeev-Popov operator is therefore $M_0^{ab}={\partial }^{\mu }(D_{\mu }{)}^{ab}$. Locally, if the condition cuts each orbit once and residual zero modes have been removed by boundary conditions, changing variables from $\omega$ to $G[A^{\omega }]-f$ proves the Faddeev-Popov gauge-orbit identityThe determinant is the functional Jacobian of this change of variables. A finite-dimensional version is $\int d\omega \delta (G(\omega )-f)=1/|\det (\partial G/\partial \omega )|$ near a unique root; the perturbative determinant has a fixed sign or phase, absorbed into normalization. This is the same change-of-variables mechanism used in Bastianelli's derivation of the Faddeev–Popov construction.

Insert the identity into the normalized functional integral for a gauge-invariant observable $\mathcal{O}[A]$:For each fixed $\omega$, change the integration variable from $A$ to $A^{\omega }$. The action, observable and measure are invariant, and the determinant is the orbit Jacobian at the resulting gauge slice. The integrand is now independent of the separately integrated orbit parameter. Its integral is $\mathrm{Vol}\mathcal{G}$, canceling the denominator. Thus, up to normalization,This assumes an anomaly-free measure and a local perturbative gauge slice; it is not an assertion of a globally unique representative on every gauge orbit.

Average over $f$ with weight $\exp [-i\int d^{4}x{f}^a{f}^a/(2\xi )]$. The delta functional sets $f=G[A]$, so the average produces the covariant gauge termNext represent the Faddeev-Popov determinant by a Grassmann Gaussian integral. Choose $M=-M_0=-{\partial }^{\mu }{D}_{\mu }$; changing $M_0$ to $-M_0$ changes the regulated determinant by a field-independent factor that is absorbed in normalization. For a Faddeev-Popov ghost field $c$ and independent antighost field $\bar{c}$,To see why it is a determinant rather than its reciprocal, diagonalize a finite-dimensional quadratic form formally: for each Grassmann pair only the term proportional to ${\bar{c}}_j{c}_j$ survives integration, supplying one eigenvalue. Multiplying the surviving eigenvalues gives $\det M$, with an overall factor of $i$ per pair. The functional version is defined with the same regulator as the rest of the theory.

Combining these steps gives the gauge-fixed Yang–Mills LagrangianThis is the gauge-fixed Yang-Mills Lagrangian in a covariant gauge. Integrating its ghost term by parts givesThe second term is the ghost-gluon vertex in the antighost convention chosen here. The Faddeev-Popov ghosts are anticommuting Lorentz scalars introduced to represent a determinant; they are not extra asymptotic particles. A closed ghost loop has a minus sign from their Grassmann parity. In a Non-Abelian gauge theory the determinant depends on $A$, so these loops are essential. In an Abelian gauge theory $D_{\mu }$ in the adjoint sector reduces to ${\partial }_{\mu }$, the determinant is field independent, and its ghost sector decouples.

The quadratic gauge field kernel at non-null momentum isDefine the mixed-index longitudinal projector of a vector field $P_L{}^{\mu }{}_{\nu }=k^{\mu }{k}_{\nu }/k^2$ and transverse projector of a vector field $P_T=1-P_L$. The kernel is $-k^2(P_T+{\xi }^{-1}{P}_L)$ and hence has inverse $-k^{-2}(P_T+\xi P_L)$. Multiplying by $i$ gives the gauge-boson propagator, off its poles,The vacuum Feynman i-epsilon prescription completes this off-pole formula; with the mostly-plus convention its physical transverse denominator is $k^2-i0$. The choice $\xi =1$ is Feynman gauge; $\xi \to 0$ is covariant Landau gauge. The invertible quadratic kernel supplies perturbation theory, while the ghost determinant ensures the gauge-orbit measure has been accounted for.

The consistency of the unphysical sector is organized by BRST symmetry. Introduce a bosonic Nakanishi-Lautrup field $B^a$ and a left-acting BRST differentialThe graded Leibniz rule and the Jacobi identity give $s^2=0$: for example $s^{2}A=D(sc)+g[Dc,c]=0$ since $D[c,c]=2[Dc,c]$, and $s^{2}c$ vanishes by the graded Jacobi identity. Nilpotence on $\bar{c}$ and $B$ is immediate. Using the odd antighost field, the same graded Leibniz rule givesEliminating $B^a$ by its algebraic field equation $B^a=-G^a/\xi$ recovers the displayed gauge-fixing and ghost terms. Thus they are BRST-exact together, and the total action preserves BRST symmetry. At the level of states, the BRST charge $Q$ implements physical states as the ghost-number-zero BRST cohomology $\ker Q/\mathrm{im}Q$. This prescription removes longitudinal, time-like and ghost excitations from physical amplitudes while allowing them on internal lines. For a BRST-closed observable, changing $\xi$ inserts $s({\bar{c}}^a{B}^a/2)$; its expectation vanishes by the BRST change-of-variables identity when the measure preserves the symmetry. This explains the independence of physical quantities from the arbitrary gauge parameter.

Finally, a local condition such as Lorenz gauge can have several solutions on the same gauge orbit, the Gribov ambiguity. The preceding Faddeev-Popov gauge-orbit identity and perturbative gauge fixing construction apply in a neighborhood where the linearized operator is invertible after residual modes are excluded. They resolve perturbative gauge overcounting and supply consistent covariant Feynman rules; they do not by themselves prove a globally unique nonperturbative gauge fixing.

Use units $\hslash =c=1$, the Minkowski metric ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$, and the Fourier transform convention $f(x)=\int d^{4}p{e}^{-ip\cdot x}\tilde{f}(p)/(2\pi {)}^4$. Here gamma matrices satisfy $\{{\gamma }^a,{\gamma }^b\}=2{\eta }^{ab}{I}_4$, and $\bar{\psi }={\psi }^{†}{\gamma }^0$ is the Dirac adjoint. A single Dirac field describes both Electrons and Positrons; these are the particle and antiparticle sectors of that field.

The free Maxwell Lagrangian and Dirac action, with a linear covariant gauge condition for the photon, giveBefore gauge fixing, the Maxwell Lagrangian has zero modes in field theory along $A_a\mapsto A_a+{\partial }_a\lambda$, so its quadratic operator cannot be inverted on all potentials. In the linear Lorenz gauge, the Faddeev-Popov determinant is $\det \square$. It is independent of $A$ and may be absorbed into the normalization; the corresponding Faddeev-Popov ghost fields have no interacting vertices in this Abelian linear gauge. A residual gauge symmetry is removed by the specified boundary conditions.

The free generating functional isThe Dirac field variables and their sources are independent Grassmann fields in this path integral. They anticommute; replacing them by ordinary commuting fields would give the wrong statistics and the wrong functional determinant. The bosonic Gaussian functional integral contributes an inverse square root of a determinant, and the Grassmann Gaussian integral contributes a determinant. If $K_A$ is the photon quadratic operator and $K_D=i{\gamma }^a{\partial }_a-m$, completing the square givesIn the last expression $D_F=i{K}_A^{-1}$ and $S_F=i{K}_D^{-1}$ are the photon propagator and Dirac propagator, with Feynman i-epsilon prescription. The products include the appropriate spacetime integrals and index contractions. Functional derivatives with respect to $J$, and consistently ordered left or right Grassmann derivatives with respect to the fermionic sources, generate the time ordering of the corresponding fields. In Feynman gauge, $\alpha =1$, the momentum-space two-point functions areThe Feynman i-epsilon prescription specifies vacuum boundary conditions rather than an arbitrary inverse of the differential operator.

The covariant photon propagator uses four potential components. Its operator formalism counterpart is Gupta-Bleuler quantization: impose $({\partial }_a{A}^a{)}^{(+)}|\mathrm{phys}⟩=0$ and take the Gupta-Bleuler null-state quotient. This leaves a positive physical state space with two transverse photon polarization vectors. The temporal and longitudinal oscillator components occur in intermediate covariant expressions; the Ward identity removes their dependence from physical amplitudes. The free Dirac field uses the canonical anticommutation relations, producing the same fermionic signs as its Grassmann fields in the path integral.

The operator-path-integral equivalence can be seen directly with a regulator. Divide time into small intervals and insert complete sets of field-coordinate states for bosons, and resolutions in fermionic coherent states for fermions. The bosonic matrix elements produce the phase-space factor $\exp i\sum [\pi \Delta \varphi -H\Delta t]$; integrating out the quadratic canonical momentum produces the bosonic action. The fermionic coherent state overlaps produce the first-order term $i{\psi }^{†}\dot{\psi }$ and the Berezin integral measure. Multiplying the short-time kernels recovers the path integral. Projecting the remote endpoints onto the Fock vacuum with an infinitesimal damping selects the same Feynman propagators as the operator formalism. Field insertions become time-ordered products under this construction. Conversely, their quadratic generating functional obeys the canonical free-field equations and has precisely the oscillator two-point functions, so its higher free correlators agree by the Wick theorem. This establishes the equivalence for the regulated free theory and order by order in the perturbation series.

To couple the matter field electromagnetically, promote its global phase symmetry to the local transformationThe gauge covariant derivative $D_a={\partial }_a+ie{A}_a$ obeys $D_a^{\prime }{\psi }^{\prime }=e^{-ie\lambda }{D}_a\psi$. Replacing ${\partial }_a$ by $D_a$ in the Dirac action therefore gives the invariant matter densityThe electromagnetic field tensor is unchanged by the local transformation. Thus the unfixed quantum electrodynamics action is gauge-invariant. The Dirac current $j^a=\bar{\psi }{\gamma }^a\psi$ is conserved by the matter equations, and the Electron and Positron excitations carry opposite charges. The added gauge fixing density selects a representative and is not itself invariant under arbitrary local transformations; it does not change gauge-invariant observables. At a free fermion vertex,For external on-shell Dirac spinors the corresponding current contraction vanishes. The quantum extension is the Ward identity, which makes physical amplitudes insensitive to adding a multiple of the photon momentum to its polarization vector. A gauge-compatible regularization preserves this vector-current identity.

With ${\mathcal{L}}_{\mathrm{int}}=-e\bar{\psi }{\gamma }^a\psi A_a$, expand $\exp (i\int {\mathcal{L}}_{\mathrm{int}})$ in powers of $e$. A term of order $n$ contains $n$ spacetime integrations and $1/n!$. Applying the Wick theorem pairs the free fields: an $A$-$A$ Wick contraction supplies a photon propagator, and a $\psi$-$\bar{\psi }$ Wick contraction supplies an oriented Dirac propagator. Each insertion supplies an interaction vertex. The permutations of contractions cancel the expansion factorials except for the Feynman-diagram symmetry factor. Interchanging Grassmann fields produces the fermionic sign, including a minus sign for every closed fermion loop. This is how Feynman diagrams arise from expectation values, rather than an extra dynamical assumption.

The normalized vacuum generating functional removes components with no external insertions. In the operator formalism, for an interacting-vacuum expectation value of an inserted product $\mathcal{O}$, the same cancellation appears asVacuum projection and the Feynman i-epsilon prescription are implicit. The linked-cluster theorem exponentiates all connected vacuum bubbles into the same factor in numerator and denominator, so it cancels. Normalize by $Z[0]$ to remove every vacuum component. This cancellation of vacuum bubbles still leaves products of disconnected diagrams that each contain external insertions. If only connected correlators are wanted, differentiate the connected generating functional $W=-i\log [Z[J,\bar{\eta },\eta ]/Z[0]]$.

The resulting momentum-space QED Feynman rules for the bare theory can be stated in Feynman gauge as follows:

The illustrated interaction vertex has incoming fermion momentum $p$, incoming photon momentum $k$, and outgoing fermion momentum $p^{\prime }=p+k$. Its solid-line arrows indicate fermion flow. The photon wavy line has no fermion-flow arrow. The propagators, the vertex $-ie{\gamma }^a$, momentum conservation, and the fermionic signs determine the perturbative amplitudes.

Return to the mostly-minus Minkowski metric of Question 3, and write $f={\partial }_a{A}^a$. The field $\xi (x)$ is assumed real and nonzero wherever its inverse occurs. Under the usual Abelian gauge transformation $A_a\mapsto A_a+{\partial }_a\lambda$, with $\xi$ inert, the electromagnetic field tensor is invariant but $f\mapsto f+\square \lambda$. The change in the added density isThis is not generally a total derivative. The action is not invariant under arbitrary gauge transformations. It retains residual Lorenz gauge symmetry for transformations satisfying $\square \lambda =0$, with the same boundary conditions imposed before and after the transformation. The action itself is undefined at $\xi =0$; that value can only be considered as a limiting gauge.

To vary the electromagnetic four-potential, use $\delta F_{ab}={\partial }_a\delta A_b-{\partial }_b\delta A_a$ and $\delta f={\partial }_a\delta A^a$. The integration by parts of both terms giveswith surface terms removed by the variational boundary conditions. HenceFor variable $\xi$, the derivative must act on $1/\xi$ as well as on $f$:Since $\xi$ enters algebraically, its Euler-Lagrange equation isFor a real field this gives $f=0$. Thus a dynamical gauge-fixing parameter imposes a constraint equation in field theory; it is not an ordinary propagating scalar. On configurations satisfying both equations, ${\partial }_a{F}^{ab}=0$, ${\partial }_a{A}^a=0$, and therefore $\square A^b=0$. There is no independent kinetic equation that determines $\xi$.

For the momentum-space equation at a prescribed constant $\xi \ne 0$, use $A_a(x)=\int d^{4}p{e}^{-ip\cdot x}{A}_a(p)/(2\pi {)}^4$. Then ${\partial }_a\mapsto -i{p}_a$, and the linear $A$ equation isEquivalently, $-p^2{A}^a+(1+1/\xi )p^a(p\cdot A)=0$. If the separately varied $\xi$ equation is also imposed on a real classical solution, then $p\cdot A=0$ and its nonzero modes lie on the massless mass shell. For constructing a full-field Green function, however, invert the fixed-background quadratic operator before imposing that on-shell constraint. Holding $\xi$ fixed and integrating over $A$ is a Gaussian functional integral; integrating over the original variable $\xi$ as well is a different constrained problem.

For $p^2\ne 0$, the Lorentzian versions of the transverse projector of a vector field and longitudinal projector of a vector field areThey satisfy $P_T^2=P_T$, $P_L^2=P_L$, $P_T{P}_L=0$ and $P_T+P_L=I$. The mixed-index operator isIts inverse follows by inverting these two scalar eigenvalues. Define $G_{ab}$ by $K^{ac}{G}_{cb}={\delta }_b^a$. The positive-sign covariant gauge-fixing inverse isThis is the algebraic inverse of the kinetic operator. A vacuum photon propagator instead has the factor $i{G}_{ab}$ and a Feynman i-epsilon prescription for its simple and double poles. If the name $G$ is used for the vacuum two-point function, include that factor $i$ consistently in its defining source equation. A retarded Green function would use different boundary conditions at the same poles. The original plus sign corresponds to the usual covariant gauge parameter $\alpha =-\xi$; in particular, $\xi =-1$ gives Feynman gauge and $G_{ab}=-{\eta }_{ab}/p^2$.

Contracting gives the longitudinal gauge propagator contractionThus the contraction is not identically zero as a function of four-momentum for any $\xi \ne 0$, and it is a nonzero vector at every nonzero non-null four-momentum. A particular component may vanish when $p_b=0$. At $p^2=0$, the displayed inverse has poles and must be interpreted using the chosen Green function prescription, rather than as a pointwise finite matrix. In the limiting covariant Landau gauge, $\xi \to 0$, the longitudinal part vanishes. A covariant photon Green function may have a longitudinal component even though physical photon polarizations are transverse.

For $L[A^a]={\partial }^{\mu }{A}_{\mu }^a$, eliminating the Nakanishi-Lautrup field gives the covariant gauge LagrangianThe first term supplies the gluon kinetic term, three-gluon vertex, and four-gluon vertex. The second makes the quadratic gauge-field operator invertible. The last supplies the ghost propagator and ghost-antighost-gluon vertex.

The nonzero one-loop contributions to the gluon propagator are a gluon bubble with two three-gluon vertices and a closed ghost bubble with two ghost-antighost-gluon vertices. A four-gluon tadpole is also present with a cutoff regulator; for massless fields it is a scaleless integral and vanishes in dimensional regularization.

For fixed $\xi \ne 0$, the Lorentzian quadratic density $-F^2/4+(\partial \cdot A{)}^2/(2\xi )$ has momentum kernel $K^{ab}=-p^2{\eta }^{ab}+(1+1/\xi )p^a{p}^b$. At $p^2\ne 0$ its inverse, defined by $K^{ac}{G}_{cb}={\delta }_b^a$, is $G_{ab}=-{\eta }_{ab}/p^2+(1+\xi )p_a{p}_b/(p^2{)}^2$. A vacuum photon propagator is $iG$ with the appropriate Feynman i-epsilon prescription. The usual parameter in the negative-sign covariant gauge term is $\alpha =-\xi$, so Feynman gauge is $\xi =-1$.