5 Fermions and Photons

The scalar particles we have studied so far had spin 0. What about higher spins? One can show a fundamental particle must be in the fundamental representation of the Lorentz group SO+⁢(1,3)\mathrm{SO}^{+}(1,3) which limits possible spins to

  • •
    ​

    j=0j=0: Higgs boson but also e.g. pion, Helium-4, Carbon-12;

  • •
    ​

    j=1/2j=1/2: quarks and leptons but also e.g. proton, neutron;

  • •
    ​

    j=1j=1: vector bosons like photons, gluons, ZZ and WW, but also deuteron, Nitrogen-14;

  • •
    ​

    j=3/2j=3/2: called a Rarita-Schwinger particle, no fundamental example has been discovered but composite particles like Lithium-7 or the Δ++\Delta^{++} exist;

  • •
    ​

    j=2j=2: any massless j=2j=2 particle can be shown to be a graviton of which we unfortunately do not have a consistent theory.

5.1 Fermions

Since fermions are particles of matter, let us consider them first. The main problem with the KG equation are its negative energy solutions. If viewed not in the context of a QFT but as a non-relativistic quantum theory, this would mean that no ground state could exist since it could always have less energy. These solutions appear because the KG is quadratic in ∂t\partial_{t}. This in turn was a consequence of the on-shell relation p2=m2p^{2}=m^{2}. Dirac’s approach was then to write the KG operator as a product

−(∂2+m2)=(i⁢γν⁢∂ν+m)⁢(i⁢γμ⁢∂μ−m).\displaystyle-(\partial^{2}+m^{2})=(i\gamma^{\nu}\partial_{\nu}+m)(i\gamma^{% \mu}\partial_{\mu}-m)\,. (263)

If we therefore choose our differential equation to be

(i⁢γμ⁢∂μ−m)⁢ψ=0,\displaystyle(i\gamma^{\mu}\partial_{\mu}-m)\psi=0\,, (264)

it is by construction linear in ∂t\partial_{t}, manifestly Lorentz invariant and fulfils the on-shell condition. Unfortunately, the γμ\gamma^{\mu} cannot be a mere number since we require

γμ⁢γν⁢∂μ∂ν=ημ⁢ν⁢∂μ∂ν.\displaystyle\gamma^{\mu}\gamma^{\nu}\partial_{\mu}\partial_{\nu}=\eta^{\mu\nu% }\partial_{\mu}\partial_{\nu}\,. (265)

Note that this does not mean that γμ⁢γν=ημ⁢ν\gamma^{\mu}\gamma^{\nu}=\eta^{\mu\nu} since the tensor ∂μ∂ν\partial_{\mu}\partial_{\nu} is symmetric under μ\mu-ν\nu exchange. Instead, we can write with the anti-commutator {A,B}=A⁢B+B⁢A\{A,B\}=AB+BA

γμ⁢γν⁢∂μ∂ν=12⁢(γμ⁢γν+γν⁢γμ)⁢∂μ∂ν=12⁢{γμ,γν}⁢∂μ∂ν.\displaystyle\gamma^{\mu}\gamma^{\nu}\partial_{\mu}\partial_{\nu}=\frac{1}{2}% \Big{(}\gamma^{\mu}\gamma^{\nu}+\gamma^{\nu}\gamma^{\mu}\Big{)}\partial_{\mu}% \partial_{\nu}=\frac{1}{2}\{\gamma^{\mu},\gamma^{\nu}\}\partial_{\mu}\partial_% {\nu}\,. (266)

This is now the defining property of the γ\gamma matrices

{γμ,γν}=2⁢ημ⁢ν.\displaystyle\{\gamma^{\mu},\gamma^{\nu}\}=2\eta^{\mu\nu}\,. (267a)
To ensure that the Hamiltonian is self-adjoint, we also require the following normalisation
(γ0)2=1and(γk)2=−1withk=1,2,3.\displaystyle\big{(}\gamma^{0}\big{)}^{2}=1\qquad\text{and}\qquad\big{(}\gamma% ^{k}\big{)}^{2}=-1\quad\text{with}\quad k=1,2,3\,. (267b)

These properties are actually sufficient for anything we may want to use γ\gamma matrices for, even without writing them down as explicit 4×44\times 4 objects.

Note that γμ\gamma^{\mu} is a Lorentz vector, i.e. a list of four 4×44\times 4 matrices. This can be made a bit clearer when using indices for this spinor space. For example with the identity matrix in spinor space II

(⁢264⁢)\displaystyle\eqref{eq:dirac} =∑b=14(i⁢γa⁢bμ−m⁢Ia⁢b)⁢ψb=0,\displaystyle=\sum_{b=1}^{4}(i\gamma^{\mu}_{ab}-m\,I_{ab})\psi_{b}=0\,, (268)
(⁢271a⁢)\displaystyle\eqref{eq:gammaident:2} =∑μ=03∑b=14γa⁢bμ⁢γμ,b⁢c=4⁢Ia⁢c,\displaystyle=\sum_{\mu=0}^{3}\sum_{b=1}^{4}\gamma^{\mu}_{ab}\gamma_{\mu,bc}=4% \,I_{ac}\,, (269)
(⁢271b⁢)\displaystyle\eqref{eq:gammaident:3} =∑μ=03∑b,c=14γa⁢bμ⁢γb⁢cν⁢γμ,c⁢d=−2⁢γa⁢dν.\displaystyle=\sum_{\mu=0}^{3}\sum_{b,c=1}^{4}\gamma^{\mu}_{ab}\gamma^{\nu}_{% bc}\gamma_{\mu,cd}=-2\gamma^{\nu}_{ad}\,. (270)

We will usually not write spinor indices and reserve the ⋅\cdot operator for a Lorentz or three-vector product.

Suggested Exercise

Show the following identities

γμ⁢γμ\displaystyle\gamma^{\mu}\gamma_{\mu} =4,\displaystyle=4\,, (271a)
γμ⁢γν⁢γμ\displaystyle\gamma^{\mu}\gamma^{\nu}\gamma_{\mu} =−2⁢γν,\displaystyle=-2\gamma^{\nu}\,, (271b)
γμ⁢γν⁢γρ⁢γμ\displaystyle\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\gamma_{\mu} =4⁢ην⁢ρ,\displaystyle=4\eta^{\nu\rho}\,, (271c)
γμ⁢γν⁢γρ⁢γσ⁢γμ\displaystyle\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\gamma^{\sigma}\gamma_{\mu} =−2⁢γσ⁢γρ⁢γν.\displaystyle=-2\gamma^{\sigma}\gamma^{\rho}\gamma^{\nu}\,. (271d)

For example

γμ⁢γμ=ημ⁢ν⁢γμ⁢γν=∗12⁢ημ⁢ν⁢(γμ⁢γν+γν⁢γμ)=ημ⁢ν⁢ημ⁢ν=4,\displaystyle\gamma^{\mu}\gamma_{\mu}=\eta_{\mu\nu}\gamma^{\mu}\gamma^{\nu}% \stackrel{{\scriptstyle*}}{{=}}\frac{1}{2}\eta_{\mu\nu}(\gamma^{\mu}\gamma^{% \nu}+\gamma^{\nu}\gamma^{\mu})=\eta_{\mu\nu}\eta^{\mu\nu}=4\,, (272)

where we have used at ∗* that the η\eta tensor is symmetric. For the next relation, we write using the anti-commutator and the previous result

γμ⁢γν⁢γμ=(2⁢ημ⁢ν−γν⁢γμ)⁢γμ=2⁢γν−4⁢γν=−2⁢γν.\displaystyle\gamma^{\mu}\gamma^{\nu}\gamma_{\mu}=(2\eta^{\mu\nu}-\gamma^{\nu}% \gamma^{\mu})\gamma_{\mu}=2\gamma^{\nu}-4\gamma^{\nu}=-2\gamma^{\nu}\,. (273)

The others follow exactly the same way.

Suggested Exercise

Using the fact that tr⁢(a⁢A+b⁢B)=a⁢tr⁢(A)+b⁢tr⁢(B){\rm tr}(aA+bB)=a{\rm tr}(A)+b{\rm tr}(B) and that tr⁢(A⋅B⁢⋯⁢C⋅D)=tr⁢(D⋅A⋅B⁢⋯⁢C){\rm tr}(A\cdot B\cdots C\cdot D)={\rm tr}(D\cdot A\cdot B\cdots C), show that

tr⁢(γμ)\displaystyle{\rm tr}(\gamma^{\mu}) =0,\displaystyle=0\,, (274a)
tr⁢(γμ1⁢⋯⁢γμk⏟odd)\displaystyle{\rm tr}(\underbrace{\gamma^{\mu_{1}}\cdots\gamma^{\mu_{k}}}_{% \text{odd}}) =0,\displaystyle=0\,, (274b)
tr⁢(γμ⁢γν)\displaystyle{\rm tr}(\gamma^{\mu}\gamma^{\nu}) =4⁢ημ⁢ν,\displaystyle=4\eta^{\mu\nu}\,, (274c)
tr⁢(γμ⁢γν⁢γρ⁢γσ)\displaystyle{\rm tr}\Big{(}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\gamma^{% \sigma}\Big{)} =4⁢(ημ⁢ν⁢ηρ⁢σ−ημ⁢ρ⁢ην⁢σ+ημ⁢σ⁢ην⁢ρ).\displaystyle=4\Big{(}\eta^{\mu\nu}\eta^{\rho\sigma}-\eta^{\mu\rho}\eta^{\nu% \sigma}+\eta^{\mu\sigma}\eta^{\nu\rho}\Big{)}\,. (274d)

We begin with

tr⁢(γμ)=1−2⁢tr⁢(γν⁢γμ⁢γν)=1−2⁢tr⁢(γν⁢γν⁢γμ)=4−2⁢tr⁢(γμ)\displaystyle{\rm tr}(\gamma^{\mu})=\frac{1}{-2}{\rm tr}(\gamma^{\nu}\gamma^{% \mu}\gamma_{\nu})=\frac{1}{-2}{\rm tr}(\gamma_{\nu}\gamma^{\nu}\gamma^{\mu})=% \frac{4}{-2}{\rm tr}(\gamma^{\mu}) (275)

which can only be satisfied if the trace is zero. Similarly, we can show γ5\gamma^{5} (cf. next section)

tr⁢(γμ1⁢γμ2⁢⋯⁢γμn)\displaystyle{\rm tr}(\gamma^{\mu_{1}}\gamma^{\mu_{2}}\cdots\gamma^{\mu_{n}}) =tr⁢(γμ1⁢γμ2⁢⋯⁢γμn⁢γ5⁢γ5)=tr⁢(γ5⁢γμ1⁢γμ2⁢⋯⁢γμn⁢γ5)=−tr⁢(γμ1⁢γ5⁢γμ2⁢⋯⁢γμn⁢γ5)\displaystyle={\rm tr}(\gamma^{\mu_{1}}\gamma^{\mu_{2}}\cdots\gamma^{\mu_{n}}% \gamma^{5}\gamma^{5})={\rm tr}(\gamma^{5}\gamma^{\mu_{1}}\gamma^{\mu_{2}}% \cdots\gamma^{\mu_{n}}\gamma^{5})=-{\rm tr}(\gamma^{\mu_{1}}\gamma^{5}\gamma^{% \mu_{2}}\cdots\gamma^{\mu_{n}}\gamma^{5})
=+tr⁢(γμ1⁢γμ2⁢γ5⁢⋯⁢γμn⁢γ5)=(−1)n⁢tr⁢(γμ1⁢γμ2⁢⋯⁢γμn⁢γ5⁢γ5).\displaystyle=+{\rm tr}(\gamma^{\mu_{1}}\gamma^{\mu_{2}}\gamma^{5}\cdots\gamma% ^{\mu_{n}}\gamma^{5})=(-1)^{n}{\rm tr}(\gamma^{\mu_{1}}\gamma^{\mu_{2}}\cdots% \gamma^{\mu_{n}}\gamma^{5}\gamma^{5})\,. (276)

If nn is odd, this means tr⁢(⋯)=−tr⁢(⋯){\rm tr}(\cdots)=-{\rm tr}(\cdots) which is only satisfied if the trace vanishes. Next,

tr⁢(γμ⁢γν)=12⁢(tr⁢(γμ⁢γν)+tr⁢(γν⁢γμ))=12⁢tr⁢({γμ,γν})=ημ⁢ν⁢tr⁢(1),\displaystyle{\rm tr}(\gamma^{\mu}\gamma^{\nu})=\frac{1}{2}\Big{(}{\rm tr}(% \gamma^{\mu}\gamma^{\nu})+{\rm tr}(\gamma^{\nu}\gamma^{\mu})\Big{)}=\frac{1}{2% }{\rm tr}\Big{(}\{\gamma^{\mu},\gamma^{\nu}\}\Big{)}=\eta^{\mu\nu}{\rm tr}(1)\,, (277)

with tr⁢(1)=4{\rm tr}(1)=4.

Suggested Exercise

Finally, show that

(γμ)†=γ0⁢γμ⁢γ0.\displaystyle(\gamma^{\mu})^{\dagger}=\gamma^{0}\gamma^{\mu}\gamma^{0}\,. (278)

We can write the Dirac equation (264) using a Hamiltonian

i⁢∂ψ∂t=(−i⁢γ0⁢γ→⋅∇→+γ0⁢m)⏟H⁢ψ.\displaystyle i\frac{\partial\psi}{\partial t}=\underbrace{\Big{(}-i\gamma^{0}% \vec{\gamma}\cdot\vec{\nabla}+\gamma^{0}m\Big{)}}_{H}\psi\,. (279)

Since we want H=H†H=H^{\dagger}, we need (γ0⁢γ→)(\gamma^{0}\vec{\gamma}) and (γ0)(\gamma^{0}) the be self-adjoint as well, justifying (267b).

5.1.1 Pauli’s fundamental theorem and basis of γ\gamma matrices

The γ\gamma matrices are fully defined through (267), i.e. we may use any set of matrices that fulfil these requirements. This means that if we have a different set of matrices (γ′)μ(\gamma^{\prime})^{\mu} that also fulfil (267), they must be related to γμ\gamma^{\mu} through a constant invertible matrix SS

(γ′)μ=S−1⁢γμ⁢S.\displaystyle(\gamma^{\prime})^{\mu}=S^{-1}\gamma^{\mu}S\,. (280)

This is called Pauli’s fundamental theorem. Its proof is not that important but it makes use of an important fact: we can write any product of γ\gamma matrices using a basis of 16=4×416=4\times 4 elements.

To pick these, it is customary to define a fifth γ\gamma matrix

γ5=i⁢γ0⁢γ1⁢γ2⁢γ3=−i4!⁢εμ⁢ν⁢ρ⁢σ⁢γμ⁢γν⁢γρ⁢γσ,\displaystyle\gamma^{5}=i\gamma^{0}\gamma^{1}\gamma^{2}\gamma^{3}=-\frac{i}{4!% }\varepsilon^{\mu\nu\rho\sigma}\gamma_{\mu}\gamma_{\nu}\gamma_{\rho}\gamma_{% \sigma}\,, (281)

with the totally anti-symmetric tensor ϵμ⁢ν⁢ρ⁢σ\epsilon^{\mu\nu\rho\sigma}, defined to be ε0123=1\varepsilon^{0123}=1.

Suggested Exercise

Show using the anti-commutation relations and the definition of γ5\gamma^{5}

(γ5)†\displaystyle(\gamma^{5})^{\dagger} =γ5,\displaystyle=\gamma^{5}\,, (282a)
(γ5)2\displaystyle(\gamma^{5})^{2} =1,\displaystyle=1\,, (282b)
{γ5,γμ}\displaystyle\{\gamma^{5},\gamma^{\mu}\} =0.\displaystyle=0\,. (282c)

For example, for the anti-commutator

{γ5,γλ}=−i4!⁢εμ⁢ν⁢ρ⁢σ⁢{γμ⁢γν⁢γρ⁢γσ,γλ},\displaystyle\{\gamma^{5},\gamma^{\lambda}\}=-\frac{i}{4!}\varepsilon_{\mu\nu% \rho\sigma}\{\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\gamma^{\sigma},\gamma^{% \lambda}\}\,, (283)

we write

γμ⁢γν⁢γρ⁢γσ⁢γλ\displaystyle\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\gamma^{\sigma}\gamma^{\lambda} =γμ⁢γν⁢γρ⁢2⁢ηλ⁢σ−γμ⁢γν⁢γσ⁢2⁢ηλ⁢ρ+γμ⁢γρ⁢γσ⁢2⁢ην⁢λ−γν⁢γρ⁢γσ⁢2⁢ημ⁢λ+γλ⁢γμ⁢γν⁢γρ⁢γσ.\displaystyle=\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}2\eta^{\lambda\sigma}-% \gamma^{\mu}\gamma^{\nu}\gamma^{\sigma}2\eta^{\lambda\rho}+\gamma^{\mu}\gamma^% {\rho}\gamma^{\sigma}2\eta^{\nu\lambda}-\gamma^{\nu}\gamma^{\rho}\gamma^{% \sigma}2\eta^{\mu\lambda}+\gamma^{\lambda}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho% }\gamma^{\sigma}\,. (284)

Therefore

{γ5,γλ}=−2⁢i4!⁢εμ⁢ν⁢ρ⁢σ⁢(γμ⁢γν⁢γρ⁢ηλ⁢σ−γμ⁢γν⁢γσ⁢ηλ⁢ρ+γμ⁢γρ⁢γσ⁢ην⁢λ−γν⁢γρ⁢γσ⁢ημ⁢λ+γλ⁢γμ⁢γν⁢γρ⁢γσ)\displaystyle\{\gamma^{5},\gamma^{\lambda}\}=-\frac{2i}{4!}\varepsilon_{\mu\nu% \rho\sigma}\Big{(}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\eta^{\lambda\sigma}-% \gamma^{\mu}\gamma^{\nu}\gamma^{\sigma}\eta^{\lambda\rho}+\gamma^{\mu}\gamma^{% \rho}\gamma^{\sigma}\eta^{\nu\lambda}-\gamma^{\nu}\gamma^{\rho}\gamma^{\sigma}% \eta^{\mu\lambda}+\gamma^{\lambda}\gamma^{\mu}\gamma^{\nu}\gamma^{\rho}\gamma^% {\sigma}\Big{)} (285)

We can write any product of γ\gamma matrices using the following basis

Γ={1⏟1,γμ⏟4,i2⁢[γμ,γν]⏟6,γ5⏟1,γμ⁢γ5⏟4}.\displaystyle\Gamma=\Big{\{}\underbrace{1}_{1},\underbrace{\gamma^{\mu}}_{4},% \underbrace{\tfrac{i}{2}[\gamma^{\mu},\gamma^{\nu}]}_{6},\underbrace{\gamma^{5% }}_{1},\underbrace{\gamma^{\mu}\gamma^{5}}_{4}\Big{\}}\,. (286)

The numbers indicate that the number of basis elements of this form. We refer to these as scalar, vector, tensor, pseudo-scalar and pseudo-vector respectively.

If we simultaneously transform the spinor ψ\psi as ψ→S−1⁢ψ\psi\to S^{-1}\psi, the Dirac equation transforms to

(i⁢γ⋅∂−m)⁢ψ→(i⁢γ′⋅∂′−m)⁢ψ′=(i⁢(S−1⁢γ⁢S)⋅∂−m)⁢S−1⁢ψ=(i⁢γ⋅∂−m)\displaystyle(i\gamma\cdot\partial-m)\psi\to(i\gamma^{\prime}\cdot\partial^{% \prime}-m)\psi^{\prime}=(i(S^{-1}\gamma S)\cdot\partial-m)S^{-1}\psi=(i\gamma% \cdot\partial-m) (287)

This means physics will always be invariant under basis change.

Even if the exact form of the γ\gamma matrices does not matter, it is sometimes helpful to have one

γ0=(12×202×202×2−12×2),γi=(02×2σi−σi02×2),\displaystyle\gamma^{0}=\begin{pmatrix}\phantom{-}1_{2\times 2}&\phantom{-}0_{% 2\times 2}\\ \phantom{-}0_{2\times 2}&-1_{2\times 2}\end{pmatrix}\,,\qquad\gamma^{i}=\begin% {pmatrix}0_{2\times 2}&\sigma^{i}\\ -\sigma^{i}&0_{2\times 2}\end{pmatrix}\,, (288)

with the Pauli matrices σi\sigma^{i}

σ1=(0110),σ2=(0−ii0),σ3=(100−1).\displaystyle\sigma^{1}=\begin{pmatrix}0&1\\ 1&0\end{pmatrix}\,,\quad\sigma^{2}=\begin{pmatrix}0&-i\\ i&0\end{pmatrix}\,,\quad\sigma^{3}=\begin{pmatrix}1&0\\ 0&-1\end{pmatrix}\,. (289)

5.1.2 Transformations of the Dirac equation

If we want to view the spinor ψ\psi as physically meaningful, we need to understand how it transforms in different frames. Consider therefore two frames described by xx and x′x^{\prime} with (x′)μ=Λνμ⁢xν(x^{\prime})^{\mu}=\Lambda^{\mu}_{\phantom{\mu}\nu}x^{\nu}. If ψ⁢(x)\psi(x) is a solution in the xx frame and ψ′⁢(x′)\psi^{\prime}(x^{\prime}) is a solution in the x′x^{\prime} frame, we must have

(i⁢γμ⁢∂∂xμ−m)⁢ψ⁢(x)=(i⁢γμ⁢∂∂(x′)μ−m)⁢ψ′⁢(x′)=0.\displaystyle\bigg{(}i\gamma^{\mu}\frac{\partial}{\partial x^{\mu}}-m\bigg{)}% \psi(x)=\bigg{(}i\gamma^{\mu}\frac{\partial}{\partial(x^{\prime})^{\mu}}-m% \bigg{)}\psi^{\prime}(x^{\prime})=0\,. (290)

We can use (27) to transform the derivative

0=(i⁢γμ⁢∂∂xμ−m)⁢ψ⁢(x)=(i⁢γμ⁢Λμν⏟(γ′)ν⁢∂∂(x′)ν−m)⁢ψ⁢(Λ−1⁢x′).\displaystyle 0=\bigg{(}i\gamma^{\mu}\frac{\partial}{\partial x^{\mu}}-m\bigg{% )}\psi(x)=\bigg{(}i\underbrace{\gamma^{\mu}\Lambda^{\nu}_{\phantom{\nu}\mu}}_{% (\gamma^{\prime})^{\nu}}\frac{\partial}{\partial(x^{\prime})^{\nu}}-m\bigg{)}% \psi(\Lambda^{-1}x^{\prime})\,. (291)

The new γ\gamma matrices (γ′)μ=Λνμ⁢γν(\gamma^{\prime})^{\mu}=\Lambda^{\mu}_{\phantom{\mu}\nu}\gamma^{\nu} still need to fulfil (267a)

{(γ′)α,(γ′)β}={Λμα⁢γμ,Λνβ⁢γν}=Λμα⁢Λνβ⁢{γμ,γν}=Λμα⁢Λνβ⁢2⁢ημ⁢ν=2⁢ηα⁢β,\displaystyle\{(\gamma^{\prime})^{\alpha},(\gamma^{\prime})^{\beta}\}=\{% \Lambda^{\alpha}_{\phantom{\alpha}\mu}\gamma^{\mu},\Lambda^{\beta}_{\phantom{% \beta}\nu}\gamma^{\nu}\}=\Lambda^{\alpha}_{\phantom{\alpha}\mu}\Lambda^{\beta}% _{\phantom{\beta}\nu}\{\gamma^{\mu},\gamma^{\nu}\}=\Lambda^{\alpha}_{\phantom{% \alpha}\mu}\Lambda^{\beta}_{\phantom{\beta}\nu}2\eta^{\mu\nu}=2\eta^{\alpha% \beta}\,, (292)

since (23) guarantees ημ⁢ν⁢Λρμ⁢Λσν=ησ⁢ρ\eta_{\mu\nu}\Lambda^{\mu}_{\phantom{\mu}\rho}\Lambda^{\nu}_{\phantom{\nu}% \sigma}=\eta_{\sigma\rho}. Pauli’s fundamental theorem proves the existence of a spinor matrix S⁢(Λ)S(\Lambda) such that

(γ′)μ=Λνμ⁢γν=S⁢(Λ)−1⁢γμ⁢S⁢(Λ).\displaystyle(\gamma^{\prime})^{\mu}=\Lambda^{\mu}_{\phantom{\mu}\nu}\gamma^{% \nu}=S(\Lambda)^{-1}\gamma^{\mu}S(\Lambda)\,. (293)

Therefore,

0=(i⁢γμ⁢∂∂xμ−m)⁢ψ⁢(x)=S⁢(Λ)−1⁢(i⁢γμ⁢∂∂(x′)ν−m)⁢S⁢(Λ)⁢ψ⁢(Λ−1⁢x′).\displaystyle 0=\bigg{(}i\gamma^{\mu}\frac{\partial}{\partial x^{\mu}}-m\bigg{% )}\psi(x)=S(\Lambda)^{-1}\bigg{(}i\gamma^{\mu}\frac{\partial}{\partial(x^{% \prime})^{\nu}}-m\bigg{)}S(\Lambda)\psi(\Lambda^{-1}x^{\prime})\,. (294)

If we left-multiply with S⁢(Λ)S(\Lambda) and identify ψ′⁢(x)′=S⁢(Λ)⁢ψ⁢(Λ−1⁢x′)\psi^{\prime}(x)^{\prime}=S(\Lambda)\psi(\Lambda^{-1}x^{\prime}), we arrive at

0=(i⁢γμ⁢∂∂(x′)ν−m)⁢ψ′⁢(x′),\displaystyle 0=\bigg{(}i\gamma^{\mu}\frac{\partial}{\partial(x^{\prime})^{\nu% }}-m\bigg{)}\psi^{\prime}(x^{\prime})\,, (295)

the transformed Dirac equation (290)

To study S⁢(Λ)S(\Lambda) consider an infinitesimal transformation Λ\Lambda which should also be infinitesimal in SS, i.e. we have

Λνμ\displaystyle\Lambda^{\mu}_{\phantom{\mu}\nu} =δνμ+ωνμ,\displaystyle=\delta^{\mu}_{\phantom{\mu}\nu}+\omega^{\mu}_{\phantom{\mu}\nu}\,, (296)
S⁢(Λ)\displaystyle S(\Lambda) =1+i⁢ωμ⁢ν⁢Σμ⁢ν.\displaystyle=1+i\omega^{\mu\nu}\Sigma_{\mu\nu}\,. (297)

Since ωμ⁢ν+ων⁢μ=0\omega_{\mu\nu}+\omega_{\nu\mu}=0, we have the same anti-symmetry for Σ\Sigma. Let us now calculate what happens to γ\gamma

Λνμ⁢γν\displaystyle\Lambda^{\mu}_{\phantom{\mu}\nu}\gamma^{\nu} =!S⁢(Λ)−1⁢γμ⁢S⁢(Λ),\displaystyle\stackrel{{\scriptstyle!}}{{=}}S(\Lambda)^{-1}\gamma^{\mu}S(% \Lambda)\,, (298)
⇔γμ+ωνμ⁢γν+𝒪⁢(ω2)\displaystyle\Leftrightarrow\qquad\qquad\gamma^{\mu}+\omega^{\mu}_{\phantom{% \mu}\nu}\gamma^{\nu}+\mathcal{O}(\omega^{2}) =(1−i⁢ωλ⁢ρ⁢Σλ⁢ρ)⁢γμ⁢(1+i⁢ωλ⁢ρ⁢Σλ⁢ρ)+𝒪⁢(ω2),\displaystyle=\Big{(}1-i\omega^{\lambda\rho}\Sigma_{\lambda\rho}\Big{)}\gamma^% {\mu}\Big{(}1+i\omega^{\lambda\rho}\Sigma_{\lambda\rho}\Big{)}+\mathcal{O}(% \omega^{2})\,, (299)
⇔12⁢ωλ⁢ρ⁢(δλμ⁢γρ−δρμ⁢γλ)=ωνμ⁢γν\displaystyle\Leftrightarrow\qquad\qquad\frac{1}{2}\omega^{\lambda\rho}\Big{(}% \delta^{\mu}_{\phantom{\mu}\lambda}\gamma_{\rho}-\delta^{\mu}_{\phantom{\mu}% \rho}\gamma_{\lambda}\Big{)}=\omega^{\mu}_{\phantom{\mu}\nu}\gamma^{\nu} =i⁢ωλ⁢ρ⁢(γμ⁢Σλ⁢ρ−Σλ⁢ρ⁢γμ)=i⁢ωλ⁢ρ⁢[γμ,Σλ⁢ρ]\displaystyle=i\omega^{\lambda\rho}\Big{(}\gamma^{\mu}\Sigma_{\lambda\rho}-% \Sigma_{\lambda\rho}\gamma^{\mu}\Big{)}=i\omega^{\lambda\rho}[\gamma^{\mu},% \Sigma_{\lambda\rho}] (300)
⇔12⁢(δλμ⁢γρ−δρμ⁢γλ)\displaystyle\Leftrightarrow\qquad\qquad\frac{1}{2}\Big{(}\delta^{\mu}_{% \phantom{\mu}\lambda}\gamma_{\rho}-\delta^{\mu}_{\phantom{\mu}\rho}\gamma_{% \lambda}\Big{)} =i⁢[γμ,Σλ⁢ρ].\displaystyle=i[\gamma^{\mu},\Sigma_{\lambda\rho}]\,. (301)

This is satisfied by

Σλ⁢ρ=−i8⁢[γλ,γρ]\displaystyle\Sigma_{\lambda\rho}=-\frac{i}{8}[\gamma_{\lambda},\gamma_{\rho}] (302)

which can be integrated to

S⁢(Λ)=exp⁡(−i8⁢ωμ⁢ν⁢[γμ,γν])forΛ∈SO+⁢(1,3).\displaystyle S(\Lambda)=\exp\Big{(}-\frac{i}{8}\omega^{\mu\nu}[\gamma_{\mu},% \gamma_{\nu}]\Big{)}\qquad\text{for}\quad\Lambda\in\mathrm{SO}^{+}(1,3)\,. (303)
Suggested Exercise

Show that

[γμ,γλ⁢γρ]=2⁢(δλμ⁢γρ−δρμ⁢γλ),\displaystyle[\gamma^{\mu},\gamma_{\lambda}\gamma_{\rho}]=2(\delta^{\mu}_{% \phantom{\mu}\lambda}\gamma_{\rho}-\delta^{\mu}_{\phantom{\mu}\rho}\gamma_{% \lambda})\,, (304)

by using the anticommutator. Use this to show that (302) is a solution of (301).

The above discussion is only valid for proper transformations that have detΛ=+1\det\Lambda=+1. To also cover improper transformations, consider

P=(10000−10000−10000−1).\displaystyle P=\begin{pmatrix}1&0&0&0\\ 0&-1&0&0\\ 0&0&-1&0\\ 0&0&0&-1\end{pmatrix}\,. (305)

This transformation flips the spatial components, i.e. it is a parity transformation. We need to fulfil

S⁢(P)−1⁢γ0⁢S⁢(P)=γ0andS⁢(P)−1⁢γi⁢S⁢(P)=−γi.\displaystyle S(P)^{-1}\gamma^{0}S(P)=\gamma^{0}\qquad\text{and}\qquad S(P)^{-% 1}\gamma^{i}S(P)=-\gamma^{i}\,. (306)

There are two possible choices for S⁢(P)S(P)

S⁢(P)=S⁢(P)−1=±γ0\displaystyle S(P)=S(P)^{-1}=\pm\gamma^{0} (307)

correctly transforms γμ\gamma^{\mu}. When acting on the wavefunction, we can have two solutions as well

ψ⁢(t,x→)→ψ′⁢(t,x→)=±γ0⁢ψ⁢(t,−x→).\displaystyle\psi(t,\vec{x})\to\psi^{\prime}(t,\vec{x})=\pm\gamma^{0}\psi(t,-% \vec{x})\,. (308)

The sign is called intrinsic parity and only starts to matter once we consider system with changing numbers of particles.

5.1.3 Solutions of the Dirac equation

To eventually construct a field theory, we will need a basis of solutions to the wave equation (264). As always, we begin with a Fourier-transformation of ψ⁢(x)\psi(x)

ψ⁢(x)=∫d3⁢p→(2⁢π)3⁢2⁢Ep→⁢(u⁢(p)⁢e−i⁢p⋅x+v⁢(p)⁢ei⁢p⋅x).\displaystyle\psi(x)=\int\frac{{\rm d}^{3}\vec{p}}{(2\pi)^{3}\sqrt{2E_{\vec{p}% }}}\Big{(}u(p){\rm e}^{-ip\cdot x}+v(p){\rm e}^{ip\cdot x}\Big{)}\,. (309)

Here the uu and vv objects are vectors in spinor space. They fulfil the momentum-space Dirac equation

(γ⋅p−m)⁢u⁢(p)=(γ⋅p+m)⁢v⁢(p)=0.\displaystyle(\gamma\cdot p-m)u(p)=(\gamma\cdot p+m)v(p)=0\,. (310)

In the restframe of the particle where p=(E,0,0,0)p=(E,0,0,0), we have

(γ0−1)⁢u⁢(0)=(γ0+1)⁢v⁢(0)=0.\displaystyle(\gamma^{0}-1)u(0)=(\gamma^{0}+1)v(0)=0\,. (311)

We can find explicit answers for the spinors using the explicit representation of (288)

u1⁢(0)=𝒩⁢(1000),u2⁢(0)=𝒩⁢(0100),v1⁢(0)=𝒩⁢(0010),v2⁢(0)=𝒩⁢(0001).\displaystyle u_{1}(0)=\mathcal{N}\begin{pmatrix}1\\ 0\\ 0\\ 0\end{pmatrix}\,,\quad u_{2}(0)=\mathcal{N}\begin{pmatrix}0\\ 1\\ 0\\ 0\end{pmatrix}\,,\quad v_{1}(0)=\mathcal{N}\begin{pmatrix}0\\ 0\\ 1\\ 0\end{pmatrix}\,,\quad v_{2}(0)=\mathcal{N}\begin{pmatrix}0\\ 0\\ 0\\ 1\end{pmatrix}\,. (312)

Note that we have not one but two solution for each direction of pp. One can show that these correspond to the two spin directions. At this point the normalisation factor 𝒩\mathcal{N} is a free parameter. It turns out that 𝒩=2⁢m\mathcal{N}=\sqrt{2m} but it is possible to modify the following equations to accommodate another normalisation.

To turn these into solutions for u⁢(p)u(p), we could perform a Lorentz boost. Alternatively, we can note that

(γ⋅p−m)⁢(γ⋅p+m)=(p2−m2),\displaystyle(\gamma\cdot p-m)(\gamma\cdot p+m)=(p^{2}-m^{2})\,, (313)

to write

ur⁢(p)∝(γ⋅p+m)⁢ur⁢(0),vr⁢(p)∝(−γ⋅p+m)⁢vr⁢(0).\displaystyle u_{r}(p)\propto(\gamma\cdot p+m)u_{r}(0)\,,\qquad v_{r}(p)% \propto(-\gamma\cdot p+m)v_{r}(0)\,. (314)
Suggested Exercise

(313) relies on the fact that (γ⋅a)⁢(γ⋅a)=a2(\gamma\cdot a)(\gamma\cdot a)=a^{2}. Proof this.

A suitable normalisation would be

ur⁢(p)=γ⋅p+m2⁢m⁢(m+Ep→)⁢ur⁢(0),vr⁢(p)=−γ⋅p+m2⁢m⁢(m+Ep→)⁢vr⁢(0),\displaystyle u_{r}(p)=\frac{\gamma\cdot p+m}{\sqrt{2m(m+E_{\vec{p}})}}u_{r}(0% )\,,\qquad v_{r}(p)=\frac{-\gamma\cdot p+m}{\sqrt{2m(m+E_{\vec{p}})}}v_{r}(0)\,, (315)

since it leads to be the “boost operator” to be normalised, i.e. for p→→0\vec{p}\to 0, the fraction has components of one. This choice results in the

ur⁢(p)†⁢γ0⁢us⁢(p)=2⁢m⁢δr⁢s,vr⁢(p)†⁢γ0⁢vs⁢(p)=−2⁢m⁢δr⁢s,ur⁢(p)†⁢γ0⁢vs⁢(p)=vs⁢(p)†⁢γ0⁢ur⁢(p)=0.\displaystyle u_{r}(p)^{\dagger}\gamma^{0}u_{s}(p)=2m\delta_{rs}\,,\qquad v_{r% }(p)^{\dagger}\gamma^{0}v_{s}(p)=-2m\delta_{rs}\,,\qquad u_{r}(p)^{\dagger}% \gamma^{0}v_{s}(p)=v_{s}(p)^{\dagger}\gamma^{0}u_{r}(p)=0\,. (316)

Further, we can show that (completeness relation)

∑r=1,2ur⁢(p)⁢ur⁢(p)†⁢γ0=γ⋅p+mand∑r=1,2vr⁢(p)⁢vr⁢(p)†⁢γ0=γ⋅p−m.\displaystyle\sum_{r=1,2}u_{r}(p)u_{r}(p)^{\dagger}\gamma^{0}=\gamma\cdot p+m% \qquad\text{and}\qquad\sum_{r=1,2}v_{r}(p)v_{r}(p)^{\dagger}\gamma^{0}=\gamma% \cdot p-m\,. (317)

5.1.4 Quantisation of the free Dirac field

To define the Lagrangian of the free Dirac field, it is helpful to first define the adjoint spinor

ψ¯=ψ†⁢γ0.\displaystyle\bar{\psi}=\psi^{\dagger}\gamma^{0}\,. (318)

With this, the Lagrangian can be written as

ℒ=ψ¯⁢(i⁢γμ⁢∂μ−m)⁢ψ,\displaystyle\mathcal{L}=\bar{\psi}(i\gamma^{\mu}\partial_{\mu}-m)\psi\,, (319)

since it results in the correct Euler-Lagrange equation for ψ¯\bar{\psi}. To see this, we calculate

∂ℒ∂(∂μψ)=i⁢ψ¯⁢γμand∂ℒ∂ψ=−ψ¯⁢m,\displaystyle\frac{\partial\mathcal{L}}{\partial(\partial_{\mu}\psi)}=i\bar{% \psi}\gamma^{\mu}\quad\text{and}\quad\frac{\partial\mathcal{L}}{\partial\psi}=% -\bar{\psi}m\,, (320)

and write

0=−i⁢ψ¯⁢∂μ←⁢γμ−m⁢ψ¯.\displaystyle 0=-i\bar{\psi}\overleftarrow{\partial_{\mu}}\gamma^{\mu}-m\bar{% \psi}\,. (321)

Here we have used the notation ∂μ←\overleftarrow{\partial_{\mu}} to indicate that the derivative is acting to the left rather than the right. We can Hermitian-conjugate this to arrive at (264).

The assosciated Hamiltonian is

H=∫d3⁢x⁢ℋ=∫d3⁢x⁢ψ¯⁢(−i⁢γ→⋅∇→+m)⁢ψ.\displaystyle H=\int{\rm d}^{3}x\ \mathcal{H}=\int{\rm d}^{3}x\ \bar{\psi}(-i% \vec{\gamma}\cdot\vec{\nabla}+m)\psi\,. (322)

Since uu and vv are vectors, we will extend our Fourier-decomposition of the fermion field slightly to split the creation and annihilation operators from the spinor vectors, i.e. we write

ψ⁢(x)=∫d3⁢p→(2⁢π)3⁢2⁢Ep→⁢∑r=1,2(ar⁢(p→)⁢ur⁢(p)⁢e−i⁢p⋅x+br⁢(p→)⁢vr⁢(p)⁢e+i⁢p⋅x).\displaystyle\psi(x)=\int\frac{{\rm d}^{3}\vec{p}}{(2\pi)^{3}\sqrt{2E_{\vec{p}% }}}\sum_{r=1,2}\Big{(}a_{r}(\vec{p})u_{r}(p){\rm e}^{-ip\cdot x}+b_{r}(\vec{p}% )v_{r}(p){\rm e}^{+ip\cdot x}\Big{)}\,. (323)

We have further flipped the propagation direction of the vv spinors to ensure positive energy. Recall how we used CCRs to quantise the KG field

[a^⁢(p→),a^†⁢(q→)]\displaystyle[\hat{a}(\vec{p}),\hat{a}^{\dagger}(\vec{q})] =(2⁢π)3⁢δ(3)⁢(p→−q→),\displaystyle=(2\pi)^{3}\delta^{(3)}(\vec{p}-\vec{q})\,, (83a)
[a^⁢(p→),a^⁢(q→)]\displaystyle[\hat{a}(\vec{p}),\hat{a}(\vec{q})] =[a^†⁢(p→),a^†⁢(q→)]=0.\displaystyle=[\hat{a}^{\dagger}(\vec{p}),\hat{a}^{\dagger}(\vec{q})]=0\,. (83b)

Here we would have

[ar⁢(p→),as†⁢(q→)]=[br⁢(p→),bs†⁢(q→)]=(2⁢π)3⁢δ(3)⁢(p→−q→)⁢δr⁢sand all other commutators zero.\displaystyle[a_{r}(\vec{p}),a_{s}^{\dagger}(\vec{q})]=[b_{r}(\vec{p}),b_{s}^{% \dagger}(\vec{q})]=(2\pi)^{3}\delta^{(3)}(\vec{p}-\vec{q})\delta_{rs}\qquad% \text{and all other commutators zero}\,. (324)

However, this will lead to a contradiction. We can show that the Hamiltonian of this theory is

H=∫d3⁢p(2⁢π)3⁢Ep→⁢∑r(ar†⁢(p→)⁢ar⁢(p→)−br†⁢(p→)⁢br⁢(p→)).\displaystyle H=\int\frac{{\rm d}^{3}p}{(2\pi)^{3}}E_{\vec{p}}\sum_{r}\Big{(}a% _{r}^{\dagger}(\vec{p})a_{r}(\vec{p})-b_{r}^{\dagger}(\vec{p})b_{r}(\vec{p})% \Big{)}\,. (325)

Since we can view ar†⁢ara_{r}^{\dagger}a_{r} and br†⁢brb_{r}^{\dagger}b_{r} as particle numbers for particles of type aa and bb, this would mean that creating more bb-type particles decreases the energy of the system.

If we instead swap bb and b†b^{\dagger}

ψ⁢(x)\displaystyle\psi(x) =∫d3⁢p→(2⁢π)3⁢2⁢Ep→⁢∑r=1,2(ar⁢(p→)⁢ur⁢(p)⁢e−i⁢p⋅x+br†⁢(p→)⁢vr⁢(p)⁢e+i⁢p⋅x),\displaystyle=\int\frac{{\rm d}^{3}\vec{p}}{(2\pi)^{3}\sqrt{2E_{\vec{p}}}}\sum% _{r=1,2}\Big{(}a_{r}(\vec{p})u_{r}(p){\rm e}^{-ip\cdot x}+b^{\dagger}_{r}(\vec% {p})v_{r}(p){\rm e}^{+ip\cdot x}\Big{)}\,, (326)
ψ¯⁢(x)\displaystyle\bar{\psi}(x) =∫d3⁢p→(2⁢π)3⁢2⁢Ep→⁢∑r=1,2(ar†⁢(p→)⁢u¯r⁢(p)⁢e+i⁢p⋅x+br⁢(p→)⁢v¯r⁢(p)⁢e−i⁢p⋅x),\displaystyle=\int\frac{{\rm d}^{3}\vec{p}}{(2\pi)^{3}\sqrt{2E_{\vec{p}}}}\sum% _{r=1,2}\Big{(}a_{r}^{\dagger}(\vec{p})\bar{u}_{r}(p){\rm e}^{+ip\cdot x}+b_{r% }(\vec{p})\bar{v}_{r}(p){\rm e}^{-ip\cdot x}\Big{)}\,, (327)

and chose anitcommutation relations

{ar⁢(p→),as†⁢(q→)}={br⁢(p→),bs†⁢(q→)}=(2⁢π)3⁢δ(3)⁢(p→−q→)⁢δr⁢sand all other anticommutators zero,\displaystyle\{a_{r}(\vec{p}),a_{s}^{\dagger}(\vec{q})\}=\{b_{r}(\vec{p}),b_{s% }^{\dagger}(\vec{q})\}=(2\pi)^{3}\delta^{(3)}(\vec{p}-\vec{q})\delta_{rs}% \qquad\text{and all other anticommutators zero}\,, (328)

we would find

H=∫d3⁢p(2⁢π)3⁢Ep→⁢∑r(ar†⁢(p→)⁢ar⁢(q→)+br†⁢(p→)⁢br⁢(p→)).\displaystyle H=\int\frac{{\rm d}^{3}p}{(2\pi)^{3}}E_{\vec{p}}\sum_{r}\Big{(}a% _{r}^{\dagger}(\vec{p})a_{r}(\vec{q})+b_{r}^{\dagger}(\vec{p})b_{r}(\vec{p})% \Big{)}\,. (329)

which means that both a†a^{\dagger} and b†b^{\dagger} create a particle of mass m2=p2m^{2}=p^{2}.

Suggested Exercise

Show (325) and (329).

Suggested Exercise

You may find the following quantum mechanics problem instructive. We normally consider the bosonic harmonic oscillator defined as

HB=ω2⁢(a†⁢a+a⁢a†)with[a,a†]=1,[a,a]=[a†,a†]=0.\displaystyle H_{B}=\frac{\omega}{2}(a^{\dagger}a+aa^{\dagger})\qquad\text{% with}\qquad[a,a^{\dagger}]=1\,,\quad[a,a]=[a^{\dagger},a^{\dagger}]=0\,. (330)

Now define the fermionic oscillator with

HF=ω2⁢(b†⁢b−b⁢b†)with{b,b†}=1,{b,b}={b†,b†}=0.\displaystyle H_{F}=\frac{\omega}{2}(b^{\dagger}b-bb^{\dagger})\qquad\text{% with}\qquad\{b,b^{\dagger}\}=1\,,\quad\{b,b\}=\{b^{\dagger},b^{\dagger}\}=0\,. (331)

Write HFH_{F} and HBH_{B} in terms of number operators NF=b†⁢bN_{F}=b^{\dagger}b and NB=a†⁢aN_{B}=a^{\dagger}a. What are the allowed eigenvalues of NFN_{F} and NBN_{B}?

You can also define a combined system H=HB+HFH=H_{B}+H_{F} with |n⟩=|nB⟩⊗|nF⟩≡|nB,nF⟩|{n}\rangle=|{n_{B}}\rangle\otimes|{n_{F}}\rangle\equiv|{n_{B},n_{F}}\rangle. This system treats bosons and fermions the same and is therefore supersymmetric. Show that the supercharge operator Q=a⁢b†Q=ab^{\dagger} fulfils

{Q,Q}={Q†,Q†}=0,ω⁢{Q,Q†}=H,[Q,H]=[Q†,H]=0.\displaystyle\{Q,Q\}=\{Q^{\dagger},Q^{\dagger}\}=0\,,\quad\omega\{Q,Q^{\dagger% }\}=H\,,\quad[Q,H]=[Q^{\dagger},H]=0\,. (332)

Therefore, QQ is a conserved quantity. Finally, apply Q⁢|nB,nF⟩Q|{n_{B},n_{F}}\rangle and explain what the operator does to a fermion or boson.

We can now define a few states. As before, we define |0⟩|{0}\rangle as the state destroyed by aa and bb

ar⁢(p→)⁢|0⟩=br⁢(p→)⁢|0⟩=0.\displaystyle a_{r}(\vec{p})|{0}\rangle=b_{r}(\vec{p})|{0}\rangle=0\,. (333)

We can also define two different one-particles states

|p→,s,+⟩=2⁢Ep→⁢as†⁢(p→)⁢|0⟩and|p→,s,−⟩=2⁢Ep→⁢bs†⁢(p→)⁢|0⟩.\displaystyle|{\vec{p},s,+}\rangle=\sqrt{2E_{\vec{p}}}a_{s}^{\dagger}(\vec{p})% |{0}\rangle\qquad\text{and}\qquad|{\vec{p},s,-}\rangle=\sqrt{2E_{\vec{p}}}b_{s% }^{\dagger}(\vec{p})|{0}\rangle\,. (334)

These states are properly normalised such that

⟨p→,s,±|q→,r,±⟩=2Ep→(2π)3δ(3)(p→−q→)δr⁢s.\displaystyle\langle{\vec{p},s,\pm}|{\vec{q},r,\pm}\rangle=2E_{\vec{p}}\,(2\pi% )^{3}\delta^{(3)}(\vec{p}-\vec{q})\delta_{rs}\,. (335)

5.1.5 Charge of the Dirac field and bilinear forms

The Dirac Lagrangian (319) has a symmetry ψ→ei⁢α⁢ψ\psi\to{\rm e}^{i\alpha}\psi which means that there must be conserved current.

This current is
jVμ=ψ¯⁢(x)⁢γμ⁢ψ⁢(x).\displaystyle j_{V}^{\mu}=\bar{\psi}(x)\gamma^{\mu}\psi(x)\,. (336a)
It is customary to also define
jS\displaystyle j_{S} =ψ¯⁢(x)⁢ψ⁢(x),\displaystyle=\bar{\psi}(x)\psi(x)\,, (336b)
j5\displaystyle j_{5} =ψ¯⁢(x)⁢γ5⁢ψ⁢(x),\displaystyle=\bar{\psi}(x)\gamma^{5}\psi(x)\,, (336c)
j5⁢Vμ\displaystyle j_{5V}^{\mu} =ψ¯⁢(x)⁢γ5⁢γμ⁢ψ⁢(x).\displaystyle=\bar{\psi}(x)\gamma^{5}\gamma^{\mu}\psi(x)\,. (336d)

It is easy to see that jVμj_{V}^{\mu} and j5⁢Vμj_{5V}^{\mu} are conserved

∂μjVμ=(∂μψ¯)⁢γμ⁢ψ+ψ¯⁢γμ⁢(∂μψ)=(i⁢m⁢ψ¯)⁢ψ+ψ¯⁢(−i⁢m⁢ψ)=0,\displaystyle\partial_{\mu}j_{V}^{\mu}=\big{(}\partial_{\mu}\bar{\psi}\big{)}% \gamma^{\mu}\psi+\bar{\psi}\gamma^{\mu}\big{(}\partial_{\mu}\psi\big{)}=(im% \bar{\psi})\psi+\bar{\psi}(-im\psi)=0\,, (337)
Suggested Exercise

Show that j5⁢Vj_{5V} is conserved as well as long as m=0m=0.

Let us also see how jj transforms under Lorentz transformation. For example,

jS→ψ¯′⁢(x′)⁢ψ′⁢(x′)=ψ¯⁢(x)⁢S⁢(Λ)−1⁢S⁢(Λ)⁢ψ⁢(x)=ψ¯⁢(x)⁢ψ⁢(x)=jS⁢(x).\displaystyle j_{S}\to\bar{\psi}^{\prime}(x^{\prime})\psi^{\prime}(x^{\prime})% =\bar{\psi}(x)S(\Lambda)^{-1}S(\Lambda)\psi(x)=\bar{\psi}(x)\psi(x)=j_{S}(x)\,. (338)
Suggested Exercise

Use the proof to show that

jVμ⁢(x)→(jV′)μ⁢(x′)=Λνμ⁢jVν⁢(x),\displaystyle j_{V}^{\mu}(x)\to(j_{V}^{\prime})^{\mu}(x^{\prime})=\Lambda^{\mu% }_{\phantom{\mu}\nu}j_{V}^{\nu}(x)\,, (339)

and similarly for j5j_{5} and j5⁢Vμj_{5V}^{\mu}.

The effect of parity is slightly more interesting. For example,

jS\displaystyle j_{S} →ψ¯′⁢(x′)⁢ψ′⁢(x′)=ψ¯⁢(x)⁢S⁢(P)−1⁢S⁢(P)⁢ψ⁢(x)=ψ¯⁢(x)⁢ψ⁢(x)=jS,\displaystyle\to\bar{\psi}^{\prime}(x^{\prime})\psi^{\prime}(x^{\prime})=\bar{% \psi}(x)S(P)^{-1}S(P)\psi(x)=\bar{\psi}(x)\psi(x)=j_{S}\,, (340)
j5\displaystyle j_{5} →ψ¯′⁢(x′)⁢(γ′)5⁢ψ′⁢(x′)=ψ¯⁢(x)⁢S⁢(P)−1⁢γ5⁢S⁢(P)⁢ψ⁢(x)=−ψ¯⁢(x)⁢γ5⁢ψ⁢(x)=−j5,\displaystyle\to\bar{\psi}^{\prime}(x^{\prime})(\gamma^{\prime})^{5}\psi^{% \prime}(x^{\prime})=\bar{\psi}(x)S(P)^{-1}\gamma^{5}S(P)\psi(x)=-\bar{\psi}(x)% \gamma^{5}\psi(x)=-j_{5}\,, (341)

and similarly for the vector currents. We have used that

S⁢(P)−1⁢γ5⁢S⁢(P)=γ0⁢γ5⁢γ0=−γ0⁢γ0⁢γ5=−γ5.\displaystyle S(P)^{-1}\gamma^{5}S(P)=\gamma^{0}\gamma^{5}\gamma^{0}=-\gamma^{% 0}\gamma^{0}\gamma^{5}=-\gamma^{5}\,. (342)

This is the original of the labels we have used for the different basis elements in (286): since jSj_{S} (jVj_{V}) transforms like a Lorentz scalar (Lorentz vector) we call it a scalar (vector) current. The “pseudo” prefix indicates that the current picks up a sign under parity conservation, the same way that e.g. the angular momentum L→=x→×p→\vec{L}=\vec{x}\times\vec{p} does.

Classically, the vector current jVj_{V} corresponds to the electromagnetic current with the charge density as ρ=jV0\rho=j_{V}^{0}. Let us calculate this current for our QFT

Q\displaystyle Q =∫d3⁢x⁢ρ⁢(x)=∫d3⁢x⁢ψ†⁢(x)⁢ψ⁢(x)\displaystyle=\int{\rm d}^{3}x\ \rho(x)=\int{\rm d}^{3}x\ \psi^{\dagger}(x)% \psi(x) (343)
=∫d3⁢x⁢d3⁢q→⁢d3⁢p→(2⁢π)6⁢2⁢Ep→⁢2⁢Eq→⁢∑r,s(ar†⁢(p→)⁢ur⁢(p)†⁢e+i⁢p⋅x+br⁢(p→)⁢vr⁢(p)†⁢e−i⁢p⋅x)⁢(ar⁢(p→)⁢ur⁢(p)⁢e−i⁢p⋅x+br†⁢(p→)⁢vr⁢(p)⁢e+i⁢p⋅x)\displaystyle=\int\frac{{\rm d}^{3}x\ {\rm d}^{3}\vec{q}\ {\rm d}^{3}\vec{p}}{% (2\pi)^{6}\sqrt{2E_{\vec{p}}2E_{\vec{q}}}}\sum_{r,s}\Big{(}a_{r}^{\dagger}(% \vec{p})u_{r}(p)^{\dagger}{\rm e}^{+ip\cdot x}+b_{r}(\vec{p})v_{r}(p)^{\dagger% }{\rm e}^{-ip\cdot x}\Big{)}\Big{(}a_{r}(\vec{p})u_{r}(p){\rm e}^{-ip\cdot x}+% b_{r}^{\dagger}(\vec{p})v_{r}(p){\rm e}^{+ip\cdot x}\Big{)} (344)
=∫d3⁢p→(2⁢π)3⁢∑r(ar†⁢(p→)⁢ar⁢(p→)+br⁢(p→)⁢br†⁢(p→))=∫d3⁢p→(2⁢π)3⁢∑r(ar†⁢(p→)⁢ar⁢(p→)−br†⁢(p→)⁢br⁢(p→)),\displaystyle=\int\frac{{\rm d}^{3}\vec{p}}{(2\pi)^{3}}\sum_{r}\Big{(}a_{r}^{% \dagger}(\vec{p})a_{r}(\vec{p})+b_{r}(\vec{p})b_{r}^{\dagger}(\vec{p})\Big{)}=% \int\frac{{\rm d}^{3}\vec{p}}{(2\pi)^{3}}\sum_{r}\Big{(}a_{r}^{\dagger}(\vec{p% })a_{r}(\vec{p})-b_{r}^{\dagger}(\vec{p})b_{r}(\vec{p})\Big{)}\,, (345)

where we have dropped yet another ∫d3⁢p⁢1\int{\rm d}^{3}p1 constant in the last step. This proves that the particles created by a†a^{\dagger} have charge Q=+1Q=+1 and the those created by b†b^{\dagger} have Q=−1Q=-1 Therefore, we call the former particles and the latter antiparticles.

Suggested Exercise

Show the above.

5.1.6 Dirac propagator

Finally, we should calculate the propagator of a fermion. To do this, we simply write

⟨0|ψa⁢(x)⁢ψ¯b⁢(y)|0⟩\displaystyle\langle{0}|\psi_{a}(x)\bar{\psi}_{b}(y)|{0}\rangle =∫d3⁢p(2⁢π)3⁢2⁢Ep→⁢ei⁢p⋅(x−y)⁢∑rur⁢(p)a⁢u¯r⁢(p)b⏟(γ⋅p+m)a⁢b,\displaystyle=\int\frac{{\rm d}^{3}p}{(2\pi)^{3}2E_{\vec{p}}}{\rm e}^{ip\cdot(% x-y)}\underbrace{\sum_{r}u_{r}(p)_{a}\bar{u}_{r}(p)_{b}}_{(\gamma\cdot p+m)_{% ab}}\,, (346)
⟨0|ψ¯b⁢(y)⁢ψa⁢(x)|0⟩\displaystyle\langle{0}|\bar{\psi}_{b}(y)\psi_{a}(x)|{0}\rangle =∫d3⁢p(2⁢π)3⁢2⁢Ep→⁢ei⁢p⋅(y−y)⁢∑rvr⁢(p)a⁢v¯r⁢(p)b⏟(γ⋅p−m)a⁢b.\displaystyle=\int\frac{{\rm d}^{3}p}{(2\pi)^{3}2E_{\vec{p}}}{\rm e}^{ip\cdot(% y-y)}\underbrace{\sum_{r}v_{r}(p)_{a}\bar{v}_{r}(p)_{b}}_{(\gamma\cdot p-m)_{% ab}}\,. (347)

Keep in mind that ψ\psi and ψ¯\bar{\psi} are operator-valued vector fields and therefore have spinor indices aa and bb, not to be confused with the operators aa and bb. By writing pμ→i⁢∂μp_{\mu}\to i\partial_{\mu}, we can pull the Dirac structure out and are left with D⁢(x−y)D(x-y) of (116)

⟨0|ψa⁢(x)⁢ψ¯b⁢(y)|0⟩\displaystyle\langle{0}|\psi_{a}(x)\bar{\psi}_{b}(y)|{0}\rangle =+(i⁢γ⋅∂x+m)a⁢b⁢D⁢(x−y),\displaystyle=+(i\gamma\cdot\partial_{x}+m)_{ab}D(x-y)\,, (348)
⟨0|ψ¯b⁢(y)⁢ψa⁢(x)|0⟩\displaystyle\langle{0}|\bar{\psi}_{b}(y)\psi_{a}(x)|{0}\rangle =−(i⁢γ⋅∂x−m)a⁢b⁢D⁢(y−x).\displaystyle=-(i\gamma\cdot\partial_{x}-m)_{ab}D(y-x)\,. (349)

Similarly, the delayed Green’s function can be written as

SRa⁢b⁢(x−y)=θ⁢(x0−y0)⁢⟨0|{ψa⁢(x),ψ¯b⁢(y)}|0⟩=(i⁢γ⋅∂x+m)a⁢b⁢DR⁢(x−y).\displaystyle S_{R}^{ab}(x-y)=\theta(x^{0}-y^{0})\langle{0}|\{\psi_{a}(x),\bar% {\psi}_{b}(y)\}|{0}\rangle=(i\gamma\cdot\partial_{x}+m)_{ab}D_{R}(x-y)\,. (350)
Suggested Exercise

Verify that SRS_{R} is indeed a Green’s function of the Dirac operator i⁢∂−mi\partial-m.

For the Fourier-transformed Green’s function, we find

S~R⁢(x−y)=i⁢(γ⋅p+m)p2−m2=iγ⋅p−m.\displaystyle\tilde{S}_{R}(x-y)=\frac{i(\gamma\cdot p+m)}{p^{2}-m^{2}}=\frac{i% }{\gamma\cdot p-m}\,. (351)

We can use the same constructions to define the Feynman propagator

SF⁢(x−y)=∫d4⁢p(2⁢π)4⁢iγ⋅p−m+i⁢ϵ⁢e−i⁢p⋅(x−y)=⟨0|T⁢ψ⁢(x)⁢ψ¯⁢(y)|0⟩.\displaystyle S_{F}(x-y)=\int\frac{{\rm d}^{4}p}{(2\pi)^{4}}\frac{i}{\gamma% \cdot p-m+i\epsilon}{\rm e}^{-ip\cdot(x-y)}=\langle{0}|T\psi(x)\bar{\psi}(y)|{% 0}\rangle\,. (352)

5.2 Vector fields

The last particle we will consider is the photon, i.e. the particle of the electromagnetic field. The classical Lagrangian is just

ℒ=−14⁢Fμ⁢ν⁢Fμ⁢νwithFμ⁢ν=∂μAν−∂νAμ.\displaystyle\mathcal{L}=-\frac{1}{4}F^{\mu\nu}F_{\mu\nu}\quad\text{with}\quad F% ^{\mu\nu}=\partial^{\mu}A^{\nu}-\partial^{\nu}A^{\mu}\,. (353)

Fμ⁢νF^{\mu\nu} is called the field-strength tensor and AμA^{\mu} the vector potential. We will not discuss how to quantise this field because AA is a gauge field which makes canonical quantisation much more complicated. This is because the conjugate momentum to A0A^{0} is zero as ℒ\mathcal{L} does not contain A˙0\dot{A}^{0} (the term ∂0A0\partial^{0}A^{0} would be in F00=∂0A0−∂0A0=0F^{00}=\partial^{0}A^{0}-\partial^{0}A^{0}=0). The way to circumvent this problem is add a gauge-fixing term to ℒ\mathcal{L} that forces a specific gauge, e.g. Lorentz gauge, i.e. ∂μAμ=0\partial_{\mu}A^{\mu}=0, in which we can write down CCRs

[Aμ⁢(x),A˙ν⁢(y)]=−i⁢ημ⁢ν⁢δ⁢(x−y).\displaystyle[A^{\mu}(x),\dot{A}^{\nu}(y)]=-i\eta^{\mu\nu}\delta(x-y)\,. (354)

This is very similar the KG field and we can almost proceed along the same lines88 8 There is one more problem related to the fact that states created by a0⁣†a^{0{\dagger}} have negative norm.. The photon field can be written as

Aμ⁢(x)=∫d3⁢p(2⁢π)3⁢2⁢Ep→⁢∑λϵλμ⁢(p)⁢[ap,λ⁢ei⁢p⋅x+ap,λ†⁢e−i⁢p⋅x].\displaystyle A^{\mu}(x)=\int\frac{{\rm d}^{3}p}{(2\pi)^{3}\sqrt{2E_{\vec{p}}}% }\sum_{\lambda}\epsilon^{\mu}_{\lambda}(p)\Big{[}a_{p,\lambda}{\rm e}^{ip\cdot x% }+a_{p,\lambda}^{\dagger}{\rm e}^{-ip\cdot x}\Big{]}\,. (355)

Here we use λ\lambda to sum over polarisations and ϵ\epsilon to denote the polarisation vector itself. Naively we would expect two polarisations. However, the gauge fixing leads to two unphysical polarisations that we also need to sum over. Like the uu and vv spinors, the polarisation vector ϵ\epsilon has a complenetess relation

∑λϵλμ⁢(p)⁢ϵλ∗ν⁢(p)=−ημ⁢ν.\displaystyle\sum_{\lambda}\epsilon_{\lambda}^{\mu}(p)\epsilon_{\lambda}^{*\nu% }(p)=-\eta^{\mu\nu}\,. (356)

The Feynman propagator of this theory is

DFμ⁢ν⁢(x−y)=∫d4⁢p(2⁢π)4⁢ημ⁢νp2+i⁢ϵ⁢e−i⁢p⋅(x−y).\displaystyle D_{F}^{\mu\nu}(x-y)=\int\frac{{\rm d}^{4}p}{(2\pi)^{4}}\frac{% \eta^{\mu\nu}}{p^{2}+i\epsilon}{\rm e}^{-ip\cdot(x-y)}\,. (357)