5 Strong Interactions

In this section we will develop the theory of the strong interactions, QCD. The major difference between QED and QCD is that the gluons are self-interacting because they also carry colour charge (unlike the charge-neutral photon). Experimental evidence suggest that each quark comes in Nc=3N_{c}=3 different colours.

5.1 The QCD Lagrangian

In QCD we have two types of particles that carry colour

  • •
    ​

    Six types of quarks (up, charm, top with charge +2/3⁢e+2/3e and down, strange, bottom with charge −1/3⁢e-1/3e) with spin 1/21/2. For each flavour, there can be Nc=3N_{c}=3 of these.

  • •
    ​

    A massless, spin 11 gluon of which there are Nc2−1=8N_{c}^{2}-1=8 versions.

The QCD Lagrangian for a single quark with mass mm is

ℒQCD=−14⁢(Fa)μ⁢ν⁢(Fa)μ⁢ν+ψ¯i⁢(i⁢γ⋅Di⁢j−m⁢δi⁢j)⁢ψj.\displaystyle\mathcal{L}_{\rm QCD}=-\frac{1}{4}(F^{a})^{\mu\nu}(F^{a})_{\mu\nu% }+\bar{\psi}_{i}(i\gamma\cdot D_{ij}-m\delta_{ij})\psi_{j}\,. (151)

The indices a=1,…,Nc2−1a=1,...,N_{c}^{2}-1, and i,j=1,…,Nci,j=1,...,N_{c} are indices of this new colour charge and summed over. The covariant derivative and field-strength tensors are now

Di⁢jμ=∂μδi⁢j+i⁢gs⁢ti⁢ja⁢Aa⁢μ,(Fa)μ⁢ν=∂μAνa−∂νAμa+gs⁢fa⁢b⁢c⁢Aμb⁢Aνc.\displaystyle D_{ij}^{\mu}=\partial^{\mu}\delta_{ij}+ig_{s}t_{ij}^{a}A^{a\mu}% \,,\qquad(F^{a})_{\mu\nu}=\partial_{\mu}A^{a}_{\nu}-\partial_{\nu}A^{a}_{\mu}+% g_{s}f^{abc}A_{\mu}^{b}A_{\nu}^{c}\,. (152)

The tat^{a} are 3×33\times 3 matrices in colour space which fulfil a commutation relation

[ta,tb]=i⁢fa⁢b⁢c⁢tc.\displaystyle[t^{a},t^{b}]=if^{abc}t^{c}\,. (153)

This is similar to the angular moment [Ji,Jj]=i⁢εi⁢j⁢k⁢Jk[J_{i},J_{j}]=i\varepsilon_{ijk}J_{k}. In place of the anti-symmetric tensor we now have the structure constants fa⁢b⁢cf^{abc} that also appear in (Fa)μ⁢ν(F^{a})_{\mu\nu} which is also anti-symmetric under changes of indices but has a larger range of indices.

Just as the JiJ_{i} generate the rotation group SU⁢(2){\rm SU}(2), the tat^{a} generate the colour symmetry group SU⁢(3){\rm SU}(3). The tat^{a} are usually represented using Gell-Mann matrices λa=2⁢ta\lambda^{a}=2t^{a}

λ1=(010100000),λ2=(0−i0i00000),λ3=(1000−10000),λ4=(001000100),λ5=(00−i000i00),λ6=(000001010),λ7=(00000−i0i0),λ8=13⁢(10001000−2).\displaystyle\begin{split}\lambda^{1}=\begin{pmatrix}0&1&0\\ 1&0&0\\ 0&0&0\end{pmatrix}\,,\quad\lambda^{2}=\begin{pmatrix}0&-i&0\\ i&0&0\\ 0&0&0\end{pmatrix}\,,\quad\lambda^{3}=\begin{pmatrix}1&0&0\\ 0&-1&0\\ 0&0&0\end{pmatrix}\,,\\ \lambda^{4}=\begin{pmatrix}0&0&1\\ 0&0&0\\ 1&0&0\end{pmatrix}\,,\quad\lambda^{5}=\begin{pmatrix}0&0&-i\\ 0&0&0\\ i&0&0\end{pmatrix}\,,\quad\lambda^{6}=\begin{pmatrix}0&0&0\\ 0&0&1\\ 0&1&0\end{pmatrix}\,,\\ \lambda^{7}=\begin{pmatrix}0&0&0\\ 0&0&-i\\ 0&i&0\end{pmatrix}\,,\quad\quad\lambda^{8}=\frac{1}{\sqrt{3}}\begin{pmatrix}1&% 0&0\\ 0&1&0\\ 0&0&-2\end{pmatrix}\,.\end{split} (154)

In practice, we are not interested in calculating one particular colour component and instead work with sums over all colours which ultimately leads to traces over the tat^{a}-matrices. We will see explicit examples of this in the sections that follow and here just collect some useful identities

tr⁢(ta)\displaystyle{\rm tr}\big{(}t^{a}\big{)} =0,\displaystyle=0\,, (155a)
tr⁢(ta⁢tb)\displaystyle{\rm tr}\big{(}t^{a}t^{b}\big{)} =12⁢δa⁢b,\displaystyle=\frac{1}{2}\delta^{ab}\,, (155b)
∑ati⁢ja⁢tj⁢ka\displaystyle\sum_{a}t_{ij}^{a}t_{jk}^{a} =CF⁢δi⁢k,\displaystyle=C_{F}\delta_{ik}\,, (155c)
∑a⁢bfa⁢b⁢c⁢fa⁢b⁢d\displaystyle\sum_{ab}f^{abc}f^{abd} =CA⁢δc⁢d,\displaystyle=C_{A}\delta^{cd}\,, (155d)

where CFC_{F} and CAC_{A} are the Casmir operators of the fundamental representation, i.e. the quarks, and the adjoint representation, i.e. the gluons, respectively. They evaluate to

CF=Nc2−12⁢NcandCA=Nc.\displaystyle C_{F}=\frac{N_{c}^{2}-1}{2N_{c}}\qquad\text{and}\qquad C_{A}=N_{% c}\,. (156)

Note that, the label of the index matters. We have used ii, jj, kk for the fundamental representation and aa, bb, cc, dd for the adjoint. This is particularly important when calculating the trace of the identity matrix δ\delta

δa⁢a=Nc2−1=8,\displaystyle\delta^{aa}=N_{c}^{2}-1=8\,, (157)
δi⁢i=Nc=3.\displaystyle\delta_{ii}=N_{c}=3\,. (158)

As in QED the QCD Lagrangian follows from a gauge symmetry and requires a gauge fixing term which we will discuss later.

5.2 Feynman rules

The Feynman rules of QCD are similar to those of QED (73) with a few simple modifications.

The propagators get a δ\delta function in the correct representation (dropping spinor indices)
For each internal fermion
A quark of momentum p flowing from j to k
=δk⁢j⁢iγ⋅p−m+i⁢ϵ,\displaystyle=\delta_{kj}\frac{i}{\gamma\cdot p-m+i\epsilon}\,, (159a)
For each internal gluon
A gluon of momentum p flowing from μ,a to ν,b
=δa⁢b⁢−i⁢gμ⁢νp2+i⁢ϵ.\displaystyle=\delta^{ab}\frac{-ig^{\mu\nu}}{p^{2}+i\epsilon}\,. (159b)
For the quark-quark-gluon vertex, we also just add a tt matrix
For each q⁢q¯⁢gq\bar{q}g vertex
A vertex connecting two quarks (indices j and k) with a gluon (index μand a)
=−i⁢gs⁢γμ⁢tk⁢ja.\displaystyle=-ig_{s}\gamma^{\mu}\,t_{kj}^{a}\,. (159c)
However, we now get new vertices from the fa⁢b⁢c⁢Aμb⁢Aνcf^{abc}A_{\mu}^{b}A_{\nu}^{c} term in (Fa)μ⁢ν(F^{a})_{\mu\nu}.
For each g⁢g⁢gggg vertex
A vertex connecting three gluons with momenta p^μ, qν, and r^ρ. The indices are a, b, and c.
=−gsfa⁢b⁢c(gμ⁢ν(p−q)ρ\displaystyle=-g_{s}f^{abc}\Big{(}g^{\mu\nu}(p-q)^{\rho}
+gν⁢ρ⁢(q−r)μ\displaystyle\qquad+g^{\nu\rho}(q-r)^{\mu}
+gρ⁢μ(r−p)ν),\displaystyle\qquad+g^{\rho\mu}(r-p)^{\nu}\Big{)}\,, (159d)
For each g⁢g⁢g⁢ggggg vertex
A vertex connecting four gluons with indices are a and μ, b and ν, c and ρ, and d and σ.
=−igs2(fe⁢a⁢bfe⁢c⁢d(gμ⁢ρgν⁢σ−gμ⁢σgν⁢ρ)\displaystyle=-ig_{s}^{2}\Big{(}f^{eab}f^{ecd}(g^{\mu\rho}g^{\nu\sigma}-g^{\mu% \sigma}g^{\nu\rho})
+fe⁢a⁢c⁢fe⁢b⁢d⁢(gμ⁢ν⁢gρ⁢σ−gμ⁢σ⁢gν⁢ρ)\displaystyle\qquad+f^{eac}f^{ebd}(g^{\mu\nu}g^{\rho\sigma}-g^{\mu\sigma}g^{% \nu\rho})
+fe⁢a⁢dfe⁢b⁢c(gμ⁢νgρ⁢σ−gμ⁢ρgν⁢σ))\displaystyle\qquad+f^{ead}f^{ebc}(g^{\mu\nu}g^{\rho\sigma}-g^{\mu\rho}g^{\nu% \sigma})\Big{)} (159e)
The rules for the external particles remain unchanged since we treat these indices as ‘open’
For each external gluon
An external gluon leaving a blob with momentum p and index μ
=ϵμ∗⁢(p)(final),\displaystyle=\epsilon^{*}_{\mu}(p)\qquad\quad\ \text{(final)}\,, (159f)
For each external gluon
An external gluon entering a blob with momentum p and index μ
=ϵμ⁢(p)(initial),\displaystyle=\epsilon_{\mu}(p)\qquad\quad\ \text{(initial)}\,, (159g)
For each external quark
An external quark leaving a blob with momentum p. The momentum flow and spinor flow are aligned
=u¯s⁢(p)(final),\displaystyle=\bar{u}_{s}(p)\qquad\quad\ \text{(final)}\,, (159h)
For each external quark
An external quark entering a blob with momentum p. The momentum flow and spinor flow are aligned
=us⁢(p)(initial),\displaystyle=u_{s}(p)\qquad\quad\ \text{(initial)}\,, (159i)
For each external antiquark
An external quark leaving a blob with momentum p. The momentum flow and spinor flow are not aligned
=vs⁢(p)(final),\displaystyle=v_{s}(p)\qquad\quad\ \text{(final)}\,, (159j)
For each external antiquark
An external quark entering a blob with momentum p. The momentum flow and spinor flow are not aligned
=v¯s⁢(p)(initial).\displaystyle=\bar{v}_{s}(p)\qquad\quad\ \text{(initial)}\,. (159k)

A few comments are now in order:

  • •
    ​

    The indices aa, bb, cc, dd, and ee are always in the adjoint representation (upper indices in tt, ff, and δ\delta). This means that they refer to the gluons. The lower indices jj and kk are in the fundamental representation (lower indices in tt and δ\delta) and refer to quarks.

  • •
    ​

    The Dirac matrix γμ\gamma^{\mu} still appears. Even though we have dropped the explicit spinor indices, these are still present. However, γμ\gamma^{\mu} and tat^{a} act on different spaces and therefore commute.

  • •
    ​

    This means that now have four kinds of indices: spinor indices, Lorentz indices, and two types of colour indices. Make sure to not mix them up!

  • •
    ​

    The new vertices arise when multiplying out the fa⁢b⁢cf^{abc} term in the field-strength tensor. In principle, this also exists in QED but there structure constant fa⁢b⁢c=0f^{abc}=0 so we do not usually write it.

5.3 Gauge invariance

Just as in QED, the Lagrangian is invariant under a gauge transformation

ψi⁢(x)→(δi⁢j−i⁢gs⁢θa⁢(x)⁢ti⁢ja)⁢ψj⁢(x),Aμa⁢(x)→Aμa⁢(x)+Dμa⁢b⁢θb⁢(x),\displaystyle\begin{split}\psi_{i}(x)&\to\big{(}\delta_{ij}-ig_{s}\theta^{a}(x% )t_{ij}^{a}\big{)}\psi_{j}(x)\,,\\ A_{\mu}^{a}(x)&\to A_{\mu}^{a}(x)+D_{\mu}^{ab}\theta^{b}(x)\,,\end{split} (160)

Here we have introduced another form of the covariant derivative. Unlike in (152) where we had it in the fundamental representation, here it is in the adjoint representation

Dμa⁢b=∂μδa⁢b+i⁢gs⁢Aμc⁢(Tc)a⁢b.\displaystyle D_{\mu}^{ab}=\partial_{\mu}\delta^{ab}+ig_{s}A_{\mu}^{c}\ (T^{c}% )^{ab}\,. (161)

The TcT^{c} is again a generator of the SU⁢(3){\rm SU}(3) colour algebra but for a different representation. In fact we have already encountered the matrices for this

(Tc)a⁢b=i⁢fa⁢c⁢b=−i⁢fa⁢b⁢c.\displaystyle(T^{c})^{ab}=if^{acb}=-if^{abc}\,. (162)

These need to fulfil the commutation relations of the algebra (153)

[Ta,Tb]=i⁢fa⁢b⁢c⁢Tc,\displaystyle[T^{a},T^{b}]=if^{abc}T^{c}\,, (163)

which is equivalent to the Jacobi identity

fa⁢b⁢c⁢fc⁢d⁢e+fb⁢d⁢c⁢fc⁢a⁢e+fd⁢a⁢c⁢fc⁢b⁢e=0.\displaystyle f^{abc}f^{cde}+f^{bdc}f^{cae}+f^{dac}f^{cbe}=0\,. (164)

As we have seen in (160) the gluon transforms in the adjoint representation while the quark transforms in the fundamental representation. Notice that the gauge transformation for AμaA_{\mu}^{a} involves the strong coupling gsg_{s}

Aμa→Aμa+∂μθa+gs⁢fa⁢b⁢c⁢θb⁢Aμc=Aμa+∂μaθa+𝒪⁢(gs),\displaystyle A_{\mu}^{a}\to A_{\mu}^{a}+\partial_{\mu}\theta^{a}+g_{s}f^{abc}% \theta^{b}A^{c}_{\mu}=A_{\mu}^{a}+\partial_{\mu}^{a}\theta^{a}+\mathcal{O}(g_{% s})\,, (165)

and only at lowest order in gsg_{s} does it reduce to the analogous transformation for QED (63). The gauge fixing works exactly like in QED, cf. (68)

ℒgf=−12⁢ξ⁢(∂μAaμ)2,\displaystyle\mathcal{L}_{\rm gf}=-\frac{1}{2\xi}(\partial_{\mu}A_{a}^{\mu})^{% 2}\,, (166)

5.3.1 Compton scattering in QCD

Let us calculate the equivalent of Compton scattering in QCD, namely

q⁢(p1)⁢q¯⁢(p2)→g⁢(p3)⁢g⁢(p4).\displaystyle q(p_{1})\bar{q}(p_{2})\to g(p_{3})g(p_{4})\,. (167)

Up to the crossing, this is equivalent to what we have calculated in Section 4.4 except that we have now an additional diagram

ℳ\displaystyle\mathcal{M} =

An quark comes in with momentum p1 and emits a gluon with momentum p3 going out. The quark line continues downwards where it joins with a anti-quark coming in (momentum p2) that has emitted a gluon with momentum p4.
+

An quark comes in with momentum p1 and emits a gluon with momentum p4 going out. The quark line continues downwards where it joins with a anti-quark coming in (momentum p2) that has emitted a gluon with momentum p3. The outgoing gluons are visibly swapped.
+

An quark and an anti-quark come in with momenta p1 and p2 and annihilate into a gluon with momenta p1+p2. The gluon continues and turns into two gluons with momenta p3 and p4.
.
\displaystyle=\begin{gathered}\includegraphics{bmlimages/notes-27.svg}% \bml@image@depth{27}\bmlDescription{ An quark comes in with momentum p1 and emits a gluon with momentum p3 going % out. The quark line continues downwards where it joins with a anti-quark % coming in (momentum p2) that has emitted a gluon with momentum p4. }\end{gathered}\quad+\quad\begin{gathered}\includegraphics{bmlimages/notes-28.% svg}\bml@image@depth{28}\bmlDescription{ An quark comes in with momentum p1 and emits a gluon with momentum p4 going % out. The quark line continues downwards where it joins with a anti-quark % coming in (momentum p2) that has emitted a gluon with momentum p3. The % outgoing gluons are visibly swapped. }\end{gathered}\quad+\quad\begin{gathered}\includegraphics{bmlimages/notes-29.% svg}\bml@image@depth{29}\bmlDescription{ An quark and an anti-quark come in with momenta p1 and p2 and annihilate into % a gluon with momenta p1+p2. The gluon continues and turns into two gluons with% momenta p3 and p4. }\end{gathered}\,.
(171)

The first two diagrams evaluate just like in (141) but with colour matrices

ℳ|Abelian=−i⁢gs2⁢ϵ∗μ⁢(p3)⁢ϵ∗ν⁢(p4)⁢v¯⁢(p2)⁢(γν⁢tb⁢γ⋅(p1−p3)(p1−p3)2⁢γμ⁢ta⁢γμ⁢ta⁢γ⋅(p1−p4)(p1−p4)2⁢γν⁢tb)⁢u⁢(p1).\displaystyle\mathcal{M}\Big{|}_{\rm Abelian}=-ig_{s}^{2}\epsilon^{*\mu}(p_{3}% )\,\epsilon^{*\nu}(p_{4})\,\,\bar{v}(p_{2})\Bigg{(}\gamma_{\nu}t^{b}\frac{% \gamma\cdot(p_{1}-p_{3})}{(p_{1}-p_{3})^{2}}\gamma_{\mu}t^{a}\gamma_{\mu}t^{a}% \frac{\gamma\cdot(p_{1}-p_{4})}{(p_{1}-p_{4})^{2}}\gamma_{\nu}t^{b}\Bigg{)}u(p% _{1})\,. (172)

If we just this Abelian part of ℳ\mathcal{M} and replaced ϵ⁢(p)\epsilon(p) with pp to verify the Ward identity (97), i.e. testing whether the replacement ϵμ→ϵμ+λ⁢p\epsilon_{\mu}\to\epsilon_{\mu}+\lambda p is valid, we would find

p3μ⁢ℳμ⁢ν|Abelian=−i⁢gs2⁢[ta,tb]⁢v¯⁢(p2)⁢γν⁢u⁢(p1)≠0.\displaystyle p_{3}^{\mu}\ \mathcal{M}_{\mu\nu}\Big{|}_{\rm Abelian}=-ig_{s}^{% 2}\Big{[}t^{a},t^{b}\Big{]}\bar{v}(p_{2})\gamma_{\nu}u(p_{1})\neq 0\,. (173)

The non-zero commutator makes these diagrams alone not gauge-invariant. Adding the final diagram gives a contribution which exactly cancels this (try this!) but yields another term proportional to p4νp_{4}^{\nu}. This vanishes when we remember the whole expression is contracted with ϵ∗ν⁢(p4)\epsilon^{*\nu}(p_{4}), and so gauge invariance is only obeyed once we project onto physical polarisations. This was not necessary in QED.

This problem is due to the sum over polarisations (102) where we did the replacement

∑α=12ϵμα⁢ϵν∗α→−gμ⁢ν.\displaystyle\sum_{\alpha=1}^{2}\epsilon_{\mu}^{\alpha}\epsilon_{\nu}^{*\alpha% }\to-g_{\mu\nu}\,. (102)

Although the right-hand summed all polarisations and not only the physical transverse ones, in actual calculations the unphysical longitudinal gluon polarisations automatically cancelled. This is no longer the case in QCD, where one has to sum strictly over physical polarisations. However, this can make calculations more cumbersome, so it might still be useful to sum over all polarisations, and to cancel in some way the unphysical degrees of freedom. How this cancellation is performed depends on the gauge. In covariant gauges, like the Feynman gauge, this is done by introducing extra fields, called the ghost fields. The alternative is to use the so-called physical gauges, that ensure that that only physical degrees of freedom propagate on shell.

5.3.2 Unitarity

Before we can understand this, we need to return to e⁢e→γ⁢γee\to\gamma\gamma which we have calculated in Section 4.4. It is possible to relate the amplitude for e⁢e→γ⁢γee\to\gamma\gamma to the process e⁢e→e⁢eee\to ee using unitarity. One can show that the imaginary part of forward scattering of e⁢e→e⁢eee\to ee is related to the amplitude square of e⁢e→γ⁢γee\to\gamma\gamma, i.e.

2⁢ℑ⁡(

A fermion coming in, emitting a photon, going down and emitting another photon before meeting up with the other incoming anti-fermion. The two photons travel horizontally and then meet up with another fermion line that is outgoing. The diagram is vaguely box-shaped.
+

A fermion coming in, emitting a photon, going down and emitting another photon before meeting up with the other incoming anti-fermion. The two photons cross over without interacting and then meet up with another fermion line that is outgoing. The diagram is shaped liked a twisted box.
+diags w/oimaginary
)
=∫dΦ⁢∑|

An quark comes in and emits a photon going out. The quark line continues downwards where it joins with a anti-quark coming in that has emitted a photon.
+

An quark comes in and emits a photon going out. The quark line continues downwards where it joins with a anti-quark coming in that has emitted a photon. The outgoing photons are visibly swapped.
|
2
.
\displaystyle 2\Im\Bigg{(}\ \begin{gathered}\includegraphics{bmlimages/notes-3% 0.svg}\bml@image@depth{30}\bmlDescription{ A fermion coming in, emitting a photon, going down and emitting another photon% before meeting up with the other incoming anti-fermion. The two photons % travel horizontally and then meet up with another fermion line that is % outgoing. The diagram is vaguely box-shaped. }\end{gathered}\ +\ \begin{gathered}\includegraphics{bmlimages/notes-31.svg}% \bml@image@depth{31}\bmlDescription{ A fermion coming in, emitting a photon, going down and emitting another photon% before meeting up with the other incoming anti-fermion. The two photons cross% over without interacting and then meet up with another fermion line that is % outgoing. The diagram is shaped liked a twisted box. }\end{gathered}\ +\genfrac{}{}{0.0pt}{0}{\text{diags w/o}}{\text{imaginary}}% \Bigg{)}=\int{\rm d}\Phi\sum\Bigg{|}\ \begin{gathered}\includegraphics{% bmlimages/notes-32.svg}\bml@image@depth{32}\bmlDescription{ An quark comes in and emits a photon going out. The quark line continues % downwards where it joins with a anti-quark coming in that has emitted a photon% . }\end{gathered}+\begin{gathered}\includegraphics{bmlimages/notes-33.svg}% \bml@image@depth{33}\bmlDescription{ An quark comes in and emits a photon going out. The quark line continues % downwards where it joins with a anti-quark coming in that has emitted a photon% . The outgoing photons are visibly swapped. }\end{gathered}\ \Bigg{|}^{2}\,.
(178)

This result is known as the optical theorem and it can be extremely powerful. For our purposes, it is enough to note that the sum on the right-hand side goes onlt over physical polarisations.

The way this is done in practice is by cutting propagators, i.e. replacing i/(p2−m2+i⁢ϵ)i/(p^{2}-m^{2}+i\epsilon) with (2⁢π)⁢Θ⁢(p0)⁢δ⁢(p2−m2)(2\pi)\Theta(p^{0})\delta(p^{2}-m^{2}). Everything on the right-hand side of the cut then gets complex-conjugated. Cuts that violate momentum conservation do not contribute. By drawing the cut as a dashed line, we can write the optical theorem as

2⁢ℑ⁡(

A fermion and an anti-fermion enter a blob. Two photons leave the blob before entering another blob. A fermion and an anti-fermion leave the blob.
)
=

A fermion and an anti-fermion enter a blob. Two photons leave the blob before entering another blob. A fermion and an anti-fermion leave the blob.
There is a dashed cut through the photons
,
\displaystyle 2\Im\Bigg{(}\ \begin{gathered}\includegraphics{bmlimages/notes-3% 4.svg}\bml@image@depth{34}\bmlDescription{ A fermion and an anti-fermion enter a blob. Two photons leave the blob before % entering another blob. A fermion and an anti-fermion leave the blob. }\end{gathered}\ \Bigg{)}={\ \begin{gathered}\includegraphics{bmlimages/notes-% 35.svg}\bml@image@depth{35}\bmlDescription{ A fermion and an anti-fermion enter a blob. Two photons leave the blob before % entering another blob. A fermion and an anti-fermion leave the blob. There is a dashed cut through the photons }\end{gathered}\ }\,,
(181)

where the blob indicates all possible ways to connect the photon, i.e. the two diagrams in tt and uu channel. In Feynman gauge we replace

−i⁢gμ⁢νk2+i⁢ϵ→−gμ⁢ν⁢(2⁢π)⁢Θ⁢(k0)⁢δ⁢(k2).\displaystyle-\frac{ig_{\mu\nu}}{k^{2}+i\epsilon}\to-g_{\mu\nu}\ (2\pi)\Theta(% k^{0})\delta(k^{2})\,. (182)

Let us call the diagram on the left side of the cut ℳμ⁢ν\mathcal{M}^{\mu\nu}. Then the Ward identity kμ⁢ℳμ⁢ν=0k_{\mu}\mathcal{M}^{\mu\nu}=0, means we can obtain the contribution for a photon kk to the imaginary part as

∫d4⁢k(2⁢π)4⁢−i⁢gμ⁢νk2+i⁢ϵ⁢ℳμ⁢ν→∫d4⁢k(2⁢π)4⁢(−gμ⁢ν)⁢(2⁢π)⁢Θ⁢(k0)⁢δ⁢(k2)⁢ℳμ⁢ν=∫d3⁢k(2⁢π)3⁢2⁢k0⁢ℳμ⁢ν⁢∑α=12ϵμα⁢(k)⁢ϵν∗α⁢(k).\displaystyle\int\frac{{\rm d}^{4}k}{(2\pi)^{4}}\frac{-ig_{\mu\nu}}{k^{2}+i% \epsilon}\mathcal{M}^{\mu\nu}\to\int\frac{{\rm d}^{4}k}{(2\pi)^{4}}(-g_{\mu\nu% })\ (2\pi)\Theta(k^{0})\delta(k^{2})\mathcal{M}^{\mu\nu}=\int\frac{{\rm d}^{3}% k}{(2\pi)^{3}2k^{0}}\mathcal{M}^{\mu\nu}\sum_{\alpha=1}^{2}\epsilon_{\mu}^{% \alpha}(k)\epsilon_{\nu}^{*\alpha}(k)\,. (183)

This verifies the explicit relation (102) and means that, in QED, we obtain (102) from unitarity.

In QCD, the fact that the Ward identity is partially broken in (172) implies that the amplitude for the process q⁢q¯→g⁢gq\bar{q}\to gg is not given by the imaginary part of the forward amplitude for q⁢q¯→q⁢q¯q\bar{q}\to q\bar{q} if we only consider gluons as internal particles. In fact, the cut forward amplitude contains the contribution of non-physical longitudinal polarisations, which do not contribute to the amplitude squared for q⁢q¯→g⁢gq\bar{q}\to gg. This would violate unitarity, so there has to be additional fields that are responsible for the cancellation of the contribution of non-physical polarisations in the imaginary part of the forward amplitude. These new fields are called ghosts.

5.3.3 Ghosts

Ghosts are scalar fields but follow Pauli’s exclusion principle and pick up a sign under swapping, just like fermions. Unlike fermions, however, they transform like gluons in the adjoint representation of SU⁢(3){\rm SU}(3). The Feynman rules for the ghosts are

For each internal ghost
A ghost of momentum p flowing from a to b
=δa⁢b⁢ip2+i⁢ϵ,\displaystyle=\delta^{ab}\frac{i}{p^{2}+i\epsilon}\,, (184a)
For each c⁢c¯⁢gc\bar{c}g vertex
A vertex connecting two ghosts (indices a and b) with a gluon (index μ)
=gs⁢fa⁢b⁢c⁢qμ.\displaystyle=g_{s}f^{abc}q^{\mu}\,. (184b)

Note that ghost only appear in internal lines, never as final states.

Consider now again the relation (181) but in QCD. We now know that there is an additional diagram with a ghost


A fermion and an anti-fermion enter a blob. Two gluons leave the blob before entering another blob. A fermion and an anti-fermion leave the blob.
There is a dashed cut through the gluons
+

An quark and an anti-quark come in with and annihilate into a gluon. The gluon continues and turns into a ghost bubble before re-combining back into a gluon. The gluon continues further and turns into two gluons.
\displaystyle{\ \begin{gathered}\includegraphics{bmlimages/notes-38.svg}% \bml@image@depth{38}\bmlDescription{ A fermion and an anti-fermion enter a blob. Two gluons leave the blob before % entering another blob. A fermion and an anti-fermion leave the blob. There is a dashed cut through the gluons }\end{gathered}\ }+\begin{gathered}\includegraphics{bmlimages/notes-39.svg}% \bml@image@depth{39}\bmlDescription{ An quark and an anti-quark come in with and annihilate into a gluon. The gluon% continues and turns into a ghost bubble before re-combining back into a gluon% . The gluon continues further and turns into two gluons. }\end{gathered}
=∫dΦ⁢∑|

An quark and an anti-quark come in and enter a blob. Out of the blob come two gluon.
|
2
\displaystyle=\int{\rm d}\Phi\sum\Bigg{|}\ \begin{gathered}\includegraphics{% bmlimages/notes-40.svg}\bml@image@depth{40}\bmlDescription{ An quark and an anti-quark come in and enter a blob. Out of the blob come two % gluon. }\end{gathered}\ \Bigg{|}^{2}
(188)

The ghost-antighost loop contributes to the imaginary part of the forward amplitude with a factor (−1)(-1), just like a normal fermion loop, so as to cancel the contribution of the unphysical longitudinal gluon polarisations when summing over all diagrams. The resulting imaginary part equals the amplitude squared for the process q⁢q¯→g⁢gq\bar{q}\to gg, integrated over the gluon phase space and summed over physical gluon polarisations, as required by unitarity of QCD.

5.3.4 Physical gauges

Alternatively, we can impose a so-called physical gauge condition on the gluon fields to eliminate unphysical polarisations from the start. This eliminates the need for ghosts, which do not interact with gluons anymore, but complicates the gluon propagator. In place of the Lorent gauge condition ∂μAμa=0\partial^{\mu}A_{\mu}^{a}=0, we impose

Aμa⁢nμ=0\displaystyle A_{\mu}^{a}n^{\mu}=0 (189)

for some arbitrary reference vector nμn^{\mu}. This is also sometimes called axial gauge. The gauge fixing Lagrangian is

ℒgf=−limξ→012⁢ξ⁢(Aμa⁢nμ)2.\displaystyle\mathcal{L}_{\rm gf}=-\lim_{\xi\to 0}\frac{1}{2}\xi\Big{(}A_{\mu}% ^{a}n^{\mu}\Big{)}^{2}\,. (190)

The gluon propagator changes to

For each internal gluon
A gluon of momentum p flowing from μ,a to ν,b
=δa⁢b⁢ip2+i⁢ϵ⁢(−gμ⁢ν+pμ⁢nν+pν⁢nμp⋅n+n2⁢pμ⁢pν(p⋅n)2).\displaystyle=\delta^{ab}\frac{i}{p^{2}+i\epsilon}\Bigg{(}-g^{\mu\nu}+\frac{p^% {\mu}n^{\nu}+p^{\nu}n^{\mu}}{p\cdot n}+n^{2}\frac{p^{\mu}p^{\nu}}{(p\cdot n)^{% 2}}\Bigg{)}\,. (191)

And the polarisation sum becomes

∑α=12ϵμα⁢ϵν∗α→−gμ⁢ν+kμ⁢nν+kν⁢nμk⋅n+n2⁢kμ⁢kν(k⋅n)2.\displaystyle\sum_{\alpha=1}^{2}\epsilon_{\mu}^{\alpha}\epsilon_{\nu}^{*\alpha% }\to-g_{\mu\nu}+\frac{k_{\mu}n_{\nu}+k_{\nu}n_{\mu}}{k\cdot n}+n^{2}\frac{k_{% \mu}k_{\nu}}{(k\cdot n)^{2}}\,. (192)

Physical quantities cannot depend on the choice of nn as long as n≠kn\neq k and one can check this explicitly.

A relevant example of a physical gauge is the light-cone gauge, in which n2=0n^{2}=0. We would for example choose k=(ω,k→)k=(\omega,\vec{k}) and n=(1,−k→/ω)n=(1,-\vec{k}/\omega) so that the n2n^{2} term drops out and we have an exact equality in (192) (rather than the “→\to” we have used previously).