Skip to content

Dirac Equation and Spinor Solutions

The previous page constructed the spinor representations of the Lorentz group and introduced gamma matrices as the algebraic bridge between vectors and spinors. We now turn that representation theory into dynamics. The result is the Dirac equation: a first-order relativistic equation whose solutions describe spin-one-half modes.

The main lesson is that the Dirac equation is not a replacement for the mass shell p2=m2p^2=m^2. It is a refined, spinorial way of imposing it. The Clifford algebra makes the operator γpm\gamma\cdot p-m a Lorentz-covariant square root of p2m2p^2-m^2, while the spinor components encode the two spin states of a massive fermion and the two chiralities that become independent in the massless limit.

For calculations, the page has a simple payoff: learn to recognize which equations define the us(p)u_s(p) spinors, which equations define the vs(p)v_s(p) spinors, and why the spin sums susuˉs\sum_su_s\bar u_s and svsvˉs\sum_sv_s\bar v_s are more useful than explicit spinor columns.

The square root of the Klein–Gordon operator

Section titled “The square root of the Klein–Gordon operator”

A free scalar mode satisfies

(p2m2)ϕ(p)=0.(p^2-m^2)\phi(p)=0.

In position space, with pμiμp_\mu\to i\partial_\mu, this is the Klein–Gordon equation

(2+m2)ϕ(x)=0,2=μμ.(\partial^2+m^2)\phi(x)=0, \qquad \partial^2=\partial_\mu\partial^\mu.

For spinors, we can impose a stronger first-order equation. Define

γp=γμpμ.\gamma\cdot p=\gamma^\mu p_\mu.

Using the Clifford algebra,

(γp)2=γμγνpμpν=12{γμ,γν}pμpν=p21,(\gamma\cdot p)^2 =\gamma^\mu\gamma^\nu p_\mu p_\nu ={1\over2}\{\gamma^\mu,\gamma^\nu\}p_\mu p_\nu =p^2\mathbf 1,

because pμpνp_\mu p_\nu is symmetric and therefore kills the commutator part of γμγν\gamma^\mu\gamma^\nu. Hence

(γpm)(γp+m)=(p2m2)1.(\gamma\cdot p-m)(\gamma\cdot p+m)=(p^2-m^2)\mathbf 1.

The free Dirac equation in momentum space is

(γpm)u(p)=0.(\gamma\cdot p-m)u(p)=0.

In position space it is

(iγμμm)ψ(x)=0.(i\gamma^\mu\partial_\mu-m)\psi(x)=0.

Acting once more with iγμμ+mi\gamma^\mu\partial_\mu+m gives

(iγμμ+m)(iγννm)ψ=(2+m2)ψ=0.(i\gamma^\mu\partial_\mu+m)(i\gamma^\nu\partial_\nu-m)\psi =-(\partial^2+m^2)\psi=0.

The overall minus sign is harmless because the expression is set to zero. The important content is that the Dirac equation implies the Klein–Gordon equation with the course-wide mostly-minus convention (+m2)ψ=0(\Box+m^2)\psi=0.

Thus every component of a Dirac spinor obeys the Klein–Gordon equation. The converse is not true: a collection of four Klein–Gordon solutions is not automatically a Dirac spinor solution. The Dirac equation keeps only the spinor configurations whose components are related by the first-order constraint.

The Dirac operator as a square root of the Klein–Gordon operator

The Clifford algebra factorizes the relativistic mass-shell polynomial. The first-order equation (γpm)u=0(\gamma\cdot p-m)u=0 implies the Klein–Gordon condition (p2m2)u=0(p^2-m^2)u=0, while preserving spinor transformation properties.

The size of the spinor is fixed by the Clifford algebra. In dd spacetime dimensions, an irreducible complex spinor has dimension 2d/22^{\lfloor d/2\rfloor}. In even dimension the complex Clifford algebra has one irreducible matrix representation up to equivalence; in odd dimension there are two choices distinguished by the sign of the highest-grade gamma product. In four dimensions this gives four complex components. On shell, the Dirac equation reduces the independent positive-energy solutions to two spin states, as required for a massive spin-one-half particle.

The orbital Lorentz generators acting on functions of spacetime are

Lμν=i(xμνxνμ).L_{\mu\nu}=i(x_\mu\partial_\nu-x_\nu\partial_\mu).

A spinor also needs intrinsic Lorentz generators acting on its components. In the global convention

σμν=i2[γμ,γν],\sigma^{\mu\nu}={i\over2}[\gamma^\mu,\gamma^\nu],

while the actual generator matrix in the spinor representation is

Σμν=12σμν=i4[γμ,γν].\Sigma^{\mu\nu}={1\over2}\sigma^{\mu\nu} ={i\over4}[\gamma^\mu,\gamma^\nu].

From the Clifford algebra one obtains

[Σμν,γρ]=i(ηνργμημργν).[\Sigma^{\mu\nu},\gamma^\rho] =i(\eta^{\nu\rho}\gamma^\mu-\eta^{\mu\rho}\gamma^\nu).

This is the algebraic statement that the gamma matrices carry a Lorentz vector index while acting on spinor indices. For a finite Lorentz transformation,

S(Λ)=exp(i2ωμνΣμν),S(\Lambda)=\exp\left(-{i\over2}\omega_{\mu\nu}\Sigma^{\mu\nu}\right),

and

S(Λ)1γρS(Λ)=Λρσγσ.S(\Lambda)^{-1}\gamma^\rho S(\Lambda)=\Lambda^\rho{}_{\sigma}\gamma^\sigma.

If x=Λxx'=\Lambda x and the spinor transforms as

ψ(x)=S(Λ)ψ(x),\psi'(x')=S(\Lambda)\psi(x),

then

(iγμμm)ψ(x)=S(Λ)(iγμμm)ψ(x).(i\gamma^\mu\partial'_\mu-m)\psi'(x') =S(\Lambda)(i\gamma^\mu\partial_\mu-m)\psi(x).

So solutions are carried into solutions. Lorentz covariance is not automatic merely because the wavefunction has four components; it follows from the compatibility between S(Λ)S(\Lambda) and the gamma matrices.

The full Lorentz generator acting on a spinor field is

Jμν=Lμν+Σμν.J_{\mu\nu}=L_{\mu\nu}+\Sigma_{\mu\nu}.

The orbital part moves the argument of the field. The spin part rotates and boosts the spinor indices.

Write a Dirac spinor in the chiral basis as

ψ=(ψLψR),\psi= \begin{pmatrix} \psi_L\\ \psi_R \end{pmatrix},

where ψL=PLψ\psi_L=P_L\psi and ψR=PRψ\psi_R=P_R\psi. On the previous page, AA was fixed by the Hermitian-vector map X=AXAX'=AXA^\dagger. With that choice, ψR\psi_R transforms with AA and ψL\psi_L with (A)1(A^\dagger)^{-1}. The important invariant statement is that the mass term couples the two inequivalent Weyl representations. The momentum-space Dirac equation becomes

(mpσpσˉm)(ψLψR)=0,\begin{pmatrix} -m&p\cdot\sigma\\ p\cdot\bar\sigma&-m \end{pmatrix} \begin{pmatrix} \psi_L\\ \psi_R \end{pmatrix}=0,

or

(pσ)ψR=mψL,(pσˉ)ψL=mψR.(p\cdot\sigma)\psi_R=m\psi_L, \qquad (p\cdot\bar\sigma)\psi_L=m\psi_R.

Equivalently,

(E1pσ)ψR=mψL,(E1+pσ)ψL=mψR.(E\mathbf 1-\mathbf p\cdot\boldsymbol\sigma)\psi_R=m\psi_L, \qquad (E\mathbf 1+\mathbf p\cdot\boldsymbol\sigma)\psi_L=m\psi_R.

The mass term couples the two chiralities. If m=0m=0, the equations decouple:

(pσˉ)ψL=0,(pσ)ψR=0.(p\cdot\bar\sigma)\psi_L=0, \qquad (p\cdot\sigma)\psi_R=0.

This is the Weyl limit. A massless left-handed spinor and a massless right-handed spinor are independent Lorentz-covariant fields; a massive Dirac fermion needs both.

The Dirac equation couples left- and right-handed Weyl spinors through the mass

The momentum matrices pσp\cdot\sigma and pσˉp\cdot\bar\sigma map between the two Weyl representations. The mass term identifies the two chiralities inside one massive Dirac spinor.

At rest, p=(m,0,0,0)p=(m,0,0,0), the positive-frequency equation gives

ψL=ψR.\psi_L=\psi_R.

Thus the rest-frame positive-frequency spinors may be written as

us(0)=m(χsχs),s=1,2,u_s(0)=\sqrt m \begin{pmatrix} \chi_s\\ \chi_s \end{pmatrix}, \qquad s=1,2,

where χs\chi_s are ordinary two-component Pauli spinors. The spin label ss is a rotation spin label, not a chirality label.

For positive-frequency solutions, take

ψ(x)=us(p)eipx,p0=Ep=p2+m2.\psi(x)=u_s(p)e^{-ip\cdot x}, \qquad p^0=E_{\mathbf p}=\sqrt{\mathbf p^2+m^2}.

Then

(γpm)us(p)=0.(\gamma\cdot p-m)u_s(p)=0.

For negative-frequency solutions, write

ψ(x)=vs(p)e+ipx,p0=Ep>0,\psi(x)=v_s(p)e^{+ip\cdot x}, \qquad p^0=E_{\mathbf p}>0,

which gives

(γp+m)vs(p)=0.(\gamma\cdot p+m)v_s(p)=0.

Here p0p^0 in the label of vs(p)v_s(p) is still positive. The sign change comes from the phase: iμe+ipx=pμe+ipxi\partial_\mu e^{+ip\cdot x}=-p_\mu e^{+ip\cdot x}. This convention is what later lets the vv branch appear as antiparticle creation rather than as a negative-energy particle.

A compact chiral-basis solution uses the unique positive-semidefinite Hermitian square roots of pσp\cdot\sigma and pσˉp\cdot\bar\sigma. They exist because, for p0=Ep>0p^0=E_{\mathbf p}>0 on shell, the eigenvalues are Epp0E_{\mathbf p}\mp|\mathbf p|\geq0:

us(p)=(pσχspσˉχs),u_s(p)= \begin{pmatrix} \sqrt{p\cdot\sigma}\,\chi_s\\ \sqrt{p\cdot\bar\sigma}\,\chi_s \end{pmatrix},

and

vs(p)=(pσηspσˉηs).v_s(p)= \begin{pmatrix} \sqrt{p\cdot\sigma}\,\eta_s\\ -\sqrt{p\cdot\bar\sigma}\,\eta_s \end{pmatrix}.

Here χs\chi_s and ηs\eta_s are orthonormal two-component spinor bases. On shell,

(pσ)(pσˉ)=m21,(p\cdot\sigma)(p\cdot\bar\sigma)=m^2\mathbf 1,

so

(pσ)pσˉ=mpσ,(pσˉ)pσ=mpσˉ.(p\cdot\sigma)\sqrt{p\cdot\bar\sigma}=m\sqrt{p\cdot\sigma}, \qquad (p\cdot\bar\sigma)\sqrt{p\cdot\sigma}=m\sqrt{p\cdot\bar\sigma}.

This verifies the two coupled Weyl equations directly.

Positive- and negative-frequency spinor branches on the mass shell

Writing a plane wave as eiqxe^{-iq\cdot x}, the uu mode has Fourier momentum q=pq=p while vs(p)e+ipxv_s(p)e^{+ip\cdot x} has q=pq=-p. The label pp in both spinors nevertheless has p0=Ep>0p^0=E_{\mathbf p}>0; after quantization the lower-frequency branch becomes the antiparticle creation term.

For explicit calculations it is often useful to choose Pauli helicity eigenspinors

σp^χλ=λχλ,λ=±1.\boldsymbol\sigma\cdot\widehat{\mathbf p}\,\chi_\lambda=\lambda\chi_\lambda, \qquad \lambda=\pm1.

The physical helicity carried by the two-component spinor is λ/2\lambda/2. With this notation,

(pσ)χλ=(Eλp)χλ,(pσˉ)χλ=(E+λp)χλ.(p\cdot\sigma)\chi_\lambda=(E-\lambda|\mathbf p|)\chi_\lambda, \qquad (p\cdot\bar\sigma)\chi_\lambda=(E+\lambda|\mathbf p|)\chi_\lambda.

The positive-frequency spinor becomes

uλ(p)=(EλpχλE+λpχλ).u_\lambda(p)= \begin{pmatrix} \sqrt{E-\lambda|\mathbf p|}\,\chi_\lambda\\ \sqrt{E+\lambda|\mathbf p|}\,\chi_\lambda \end{pmatrix}.

The square roots are exactly what makes the two first-order equations true, because

(E+λp)(Eλp)=E2p2=m2.(E+\lambda|\mathbf p|)(E-\lambda|\mathbf p|)=E^2-\mathbf p^2=m^2.

The standard spinor normalization used in the next page is

uˉr(p)us(p)=2mδrs,vˉr(p)vs(p)=2mδrs,\bar u_r(p)u_s(p)=2m\delta_{rs}, \qquad \bar v_r(p)v_s(p)=-2m\delta_{rs},

and

ur(p)us(p)=2Epδrs,vr(p)vs(p)=2Epδrs.u_r(p)^\dagger u_s(p)=2E_{\mathbf p}\delta_{rs}, \qquad v_r(p)^\dagger v_s(p)=2E_{\mathbf p}\delta_{rs}.

The first line is Lorentz invariant. The second line is the equal-time Hilbert-space normalization used in canonical quantization.

The spin sums are

s=12us(p)uˉs(p)=γp+m,s=12vs(p)vˉs(p)=γpm.\sum_{s=1}^2 u_s(p)\bar u_s(p)=\gamma\cdot p+m, \qquad \sum_{s=1}^2 v_s(p)\bar v_s(p)=\gamma\cdot p-m.

These identities are usually more important than the explicit spinors. They convert sums over external spin states into traces of gamma matrices, and they are independent of the particular spin basis used for the two states. In practice, one often chooses helicity spinors for external particles but immediately replaces spin sums by the covariant numerators above. Whenever a calculation seems to differ by a factor of 2m2m or 2Ep2E_{\mathbf p}, check the spinor normalization before changing the physics.

For m>0m>0, the positive-frequency numerator defines a projector:

Λ+(p)=γp+m2m,Λ+2=Λ+.\Lambda_+(p)={\gamma\cdot p+m\over2m}, \qquad \Lambda_+^2=\Lambda_+.

Indeed,

(γp+m)2=p2+2mγp+m2=2m(γp+m),(\gamma\cdot p+m)^2=p^2+2m\gamma\cdot p+m^2=2m(\gamma\cdot p+m),

where p2=m2p^2=m^2 was used. Since trΛ+=2\operatorname{tr}\Lambda_+=2, this projector has rank two, matching the two spin states of a massive positive-energy fermion. The corresponding projector onto solutions of (γp+m)v=0(\gamma\cdot p+m)v=0 at the same positive-energy momentum label is

Λ(p)=mγp2m.\Lambda_-(p)={m-\gamma\cdot p\over2m}.

It obeys Λvs=vs\Lambda_-v_s=v_s and Λ2=Λ\Lambda_-^2=\Lambda_-. Because the covariant normalization of the vv spinors is negative,

svs(p)vˉs(p)=γpm=2mΛ(p).\sum_s v_s(p)\bar v_s(p) =\gamma\cdot p-m =-2m\Lambda_-(p).

Neither Λ+\Lambda_+ nor Λ\Lambda_- should be used at m=0m=0, where division by 2m2m is singular; the chiral and helicity projectors are the appropriate massless objects.

Spin sums as covariant projectors onto Dirac solution spaces

Spin sums turn basis-dependent spinor products into covariant matrices. With the stated normalization, susuˉs=2mΛ+\sum_su_s\bar u_s=2m\Lambda_+ while svsvˉs=2mΛ\sum_sv_s\bar v_s=-2m\Lambda_-; the minus sign in the latter is tied to vˉsvs=2m\bar v_sv_s=-2m.

The same numerator appears in the free Dirac propagator,

SF(p)=i(γp+m)p2m2+iϵ,S_F(p)={i(\gamma\cdot p+m)\over p^2-m^2+i\epsilon},

which will be derived after quantizing the Dirac field.

Chirality and helicity in the massless limit

Section titled “Chirality and helicity in the massless limit”

When m=0m=0, the Dirac equation separates into two Weyl equations:

iσˉμμψL=0,iσμμψR=0.i\bar\sigma^\mu\partial_\mu\psi_L=0, \qquad i\sigma^\mu\partial_\mu\psi_R=0.

For a positive-energy mode with p0=pp^0=|\mathbf p|,

pσˉ=p(1+σp^),pσ=p(1σp^).p\cdot\bar\sigma=|\mathbf p|(1+\boldsymbol\sigma\cdot\widehat{\mathbf p}), \qquad p\cdot\sigma=|\mathbf p|(1-\boldsymbol\sigma\cdot\widehat{\mathbf p}).

Therefore

σp^ψL=ψL,σp^ψR=+ψR.\boldsymbol\sigma\cdot\widehat{\mathbf p}\,\psi_L=-\psi_L, \qquad \boldsymbol\sigma\cdot\widehat{\mathbf p}\,\psi_R=+\psi_R.

Since the spin operator is S=σ/2\mathbf S=\boldsymbol\sigma/2, a positive-energy left-handed Weyl mode has helicity 1/2-1/2, while a positive-energy right-handed Weyl mode has helicity +1/2+1/2.

Massless chirality and helicity for positive-energy spinors

For m=0m=0, chirality and helicity lock together for positive-energy particle modes. Left-handed modes have spin antiparallel to momentum; right-handed modes have spin parallel to momentum.

For massive fermions the distinction matters. Chirality is a Lorentz-representation label, implemented by

PL=1γ52,PR=1+γ52.P_L={1-\gamma^5\over2}, \qquad P_R={1+\gamma^5\over2}.

Helicity is the projection of spin along momentum. A massive particle can be viewed from a frame that overtakes it and reverses its momentum, so helicity is not Lorentz invariant for massive particles. A massless particle cannot be overtaken, and helicity becomes invariant. This is why weak-interaction phrases such as “left-handed fermion” must be read carefully: for a massive fermion, a chirality projection and a helicity measurement are related but not identical.

The ordinary Hermitian inner product ψψ\psi^\dagger\psi is invariant under spatial rotations, but not under boosts, because finite-dimensional spinor boost matrices are not unitary. The gamma matrices instead obey

(γ0)=γ0,(γi)=γi,(\gamma^0)^\dagger=\gamma^0, \qquad (\gamma^i)^\dagger=-\gamma^i,

or, equivalently,

(γμ)=γ0γμγ0.(\gamma^\mu)^\dagger=\gamma^0\gamma^\mu\gamma^0.

This identity implies

S(Λ)γ0S(Λ)=γ0.S(\Lambda)^\dagger\gamma^0S(\Lambda)=\gamma^0.

The invariant matrix is therefore γ0\gamma^0, not the identity. This motivates the Lorentz-covariant adjoint

ψˉ=ψγ0.\bar\psi=\psi^\dagger\gamma^0.

Indeed, if ψ=Sψ\psi'=S\psi, then

ψˉ=ψSγ0=ψˉS1.\bar\psi' =\psi^\dagger S^\dagger\gamma^0 =\bar\psi S^{-1}.

It follows immediately that ψˉψ\bar\psi\psi is a Lorentz scalar and ψˉγμψ\bar\psi\gamma^\mu\psi is a Lorentz vector.

The free Dirac Lagrangian is

L=ψˉ(iγμμm)ψ.\mathcal L=\bar\psi(i\gamma^\mu\partial_\mu-m)\psi.

Varying ψˉ\bar\psi gives the Dirac equation,

(iγμμm)ψ=0.(i\gamma^\mu\partial_\mu-m)\psi=0.

Varying ψ\psi and integrating by parts gives the adjoint equation

i(μψˉ)γμ+mψˉ=0.i(\partial_\mu\bar\psi)\gamma^\mu+m\bar\psi=0.

The Lagrangian is invariant under the global phase symmetry

ψeiαψ,ψˉψˉeiα.\psi\mapsto e^{i\alpha}\psi, \qquad \bar\psi\mapsto \bar\psi e^{-i\alpha}.

The associated current is

jμ=ψˉγμψ.j^\mu=\bar\psi\gamma^\mu\psi.

Using the Dirac equation and its adjoint,

μjμ=(μψˉ)γμψ+ψˉγμμψ=imψˉψimψˉψ=0.\partial_\mu j^\mu =(\partial_\mu\bar\psi)\gamma^\mu\psi+\bar\psi\gamma^\mu\partial_\mu\psi =im\bar\psi\psi-im\bar\psi\psi=0.

The time component is

j0=ψˉγ0ψ=ψψ.j^0=\bar\psi\gamma^0\psi=\psi^\dagger\psi.

As a classical wave equation this gives a positive density. In quantum field theory the deeper interpretation is that jμj^\mu is a charge current, while the negative-frequency modes become antiparticle creation modes.

The same gamma algebra organizes all local Dirac bilinears. A basis for the 4×44\times4 complex matrices contains

1+4+6+4+1=161+4+6+4+1=16

independent elements:

Gamma structureBilinearProper-Lorentz type
1\mathbf 1ψˉψ\bar\psi\psiscalar
γμ\gamma^\muψˉγμψ\bar\psi\gamma^\mu\psivector
σμν\sigma^{\mu\nu}ψˉσμνψ\bar\psi\sigma^{\mu\nu}\psiantisymmetric tensor
γμγ5\gamma^\mu\gamma^5ψˉγμγ5ψ\bar\psi\gamma^\mu\gamma^5\psiaxial vector
γ5\gamma^5ψˉγ5ψ\bar\psi\gamma^5\psipseudoscalar

The adjectives “axial” and “pseudo” distinguish behavior under parity; under the proper orthochronous Lorentz group they transform like a vector and scalar, respectively. This 1+4+6+4+11+4+6+4+1 decomposition is also a useful completeness check for gamma-matrix manipulations.

The Dirac equation

(iγμμm)ψ=0(i\gamma^\mu\partial_\mu-m)\psi=0

combines Lorentz covariance, the Clifford algebra, and spin-one-half representation theory. Squaring the operator gives the Klein–Gordon equation, so Dirac spinors live on the same mass shell as scalar particles, but with spinor polarizations attached to each momentum.

A massive Dirac spinor contains left- and right-handed Weyl components coupled by the mass. In the massless limit the components decouple, and chirality becomes tied to helicity for positive-energy particle modes.

The positive-frequency spinors us(p)u_s(p) obey (γpm)us=0(\gamma\cdot p-m)u_s=0, while the negative-frequency spinors vs(p)v_s(p) obey (γp+m)vs=0(\gamma\cdot p+m)v_s=0. Their spin sums,

sus(p)uˉs(p)=γp+m,svs(p)vˉs(p)=γpm,\sum_s u_s(p)\bar u_s(p)=\gamma\cdot p+m, \qquad \sum_s v_s(p)\bar v_s(p)=\gamma\cdot p-m,

are the practical heart of spinor calculations and foreshadow the numerator of the fermion propagator. In both sums p0p^0 is the positive on-shell energy.

Chirality and helicity differ at nonzero mass. They coincide for positive-energy massless particle modes, but a boost can reverse the momentum of a massive particle without changing its chirality.

Spinor boost matrices need not be unitary. The finite-dimensional matrix S(Λ)S(\Lambda) is not generally unitary for boosts. This does not contradict unitarity of the Lorentz transformation on the Hilbert space.

The vs(p)v_s(p) spinors are not negative-probability states. Their label still has p0=+Epp^0=+E_{\mathbf p}. The negative Fourier frequency is carried by e+ipxe^{+ip\cdot x}, and after quantization this mode multiplies an antiparticle creation operator.

The momentum label in the vv spin sum has positive energy. In svs(p)vˉs(p)=γpm\sum_s v_s(p)\bar v_s(p)=\gamma\cdot p-m, do not silently replace p0=+Epp^0=+E_{\mathbf p} by Ep-E_{\mathbf p}.

Spinor normalizations fix the completeness relations. The displayed spin sums assume uˉsus=2mδss\bar u_su_{s'}=2m\delta_{ss'} and vˉsvs=2mδss\bar v_sv_{s'}=-2m\delta_{ss'}. Changing normalization changes the prefactors.

The Dirac adjoint is essential. The scalar ψˉψ\bar\psi\psi and vector ψˉγμψ\bar\psi\gamma^\mu\psi are Lorentz covariant; ψψ\psi^\dagger\psi alone is neither a Lorentz scalar nor boost invariant.

Explicit columns are basis dependent. Spin sums and bilinears are usually safer covariant objects. Projectors divided by 2m2m are only defined for m>0m>0.

Show directly that the Dirac equation implies the Klein–Gordon equation for each spinor component.

Solution

Start from

(iγμμm)ψ=0.(i\gamma^\mu\partial_\mu-m)\psi=0.

Act on the left with iγνν+mi\gamma^\nu\partial_\nu+m:

0=(iγνν+m)(iγμμm)ψ.0=(i\gamma^\nu\partial_\nu+m)(i\gamma^\mu\partial_\mu-m)\psi.

The cross terms cancel, leaving

0=γνγμνμψm2ψ.0=-\gamma^\nu\gamma^\mu\partial_\nu\partial_\mu\psi-m^2\psi.

Because νμ\partial_\nu\partial_\mu is symmetric, the commutator part of γνγμ\gamma^\nu\gamma^\mu drops out:

γνγμνμ=12{γν,γμ}νμ=2.\gamma^\nu\gamma^\mu\partial_\nu\partial_\mu ={1\over2}\{\gamma^\nu,\gamma^\mu\}\partial_\nu\partial_\mu =\partial^2.

Therefore

(2+m2)ψ=0.(\partial^2+m^2)\psi=0.

Exercise 2: chiral equations imply the mass shell

Section titled “Exercise 2: chiral equations imply the mass shell”

Using the chiral-basis equations

(E1pσ)ψR=mψL,(E1+pσ)ψL=mψR,(E\mathbf 1-\mathbf p\cdot\boldsymbol\sigma)\psi_R=m\psi_L, \qquad (E\mathbf 1+\mathbf p\cdot\boldsymbol\sigma)\psi_L=m\psi_R,

verify the helicity-basis spinor

uλ(p)=(EλpχλE+λpχλ),σp^χλ=λχλ,λ=±1.u_\lambda(p)= \begin{pmatrix} \sqrt{E-\lambda|\mathbf p|}\,\chi_\lambda\\ \sqrt{E+\lambda|\mathbf p|}\,\chi_\lambda \end{pmatrix}, \qquad \boldsymbol\sigma\cdot\widehat{\mathbf p}\,\chi_\lambda=\lambda\chi_\lambda, \qquad \lambda=\pm1.
Solution

Let

ψL=Aλχλ,ψR=Bλχλ,\psi_L=A_\lambda\chi_\lambda, \qquad \psi_R=B_\lambda\chi_\lambda,

with

Aλ=Eλp,Bλ=E+λp.A_\lambda=\sqrt{E-\lambda|\mathbf p|}, \qquad B_\lambda=\sqrt{E+\lambda|\mathbf p|}.

The first Weyl equation gives

(E1pσ)ψR=(Eλp)Bλχλ.(E\mathbf 1-\mathbf p\cdot\boldsymbol\sigma)\psi_R =(E-\lambda|\mathbf p|)B_\lambda\chi_\lambda.

Using E2p2=m2E^2-\mathbf p^2=m^2,

(Eλp)Bλ=(Eλp)E+λp=mEλp=mAλ.(E-\lambda|\mathbf p|)B_\lambda =(E-\lambda|\mathbf p|)\sqrt{E+\lambda|\mathbf p|} =m\sqrt{E-\lambda|\mathbf p|}=mA_\lambda.

So the first equation holds. Similarly,

(E1+pσ)ψL=(E+λp)Aλχλ=mBλχλ,(E\mathbf 1+\mathbf p\cdot\boldsymbol\sigma)\psi_L =(E+\lambda|\mathbf p|)A_\lambda\chi_\lambda=mB_\lambda\chi_\lambda,

so the second equation also holds.

Exercise 3: the positive-energy spin-sum projector

Section titled “Exercise 3: the positive-energy spin-sum projector”

Show that

Λ+(p)=γp+m2m\Lambda_+(p)={\gamma\cdot p+m\over2m}

is a projector on shell and acts as the identity on any us(p)u_s(p).

Solution

Since (γpm)us(p)=0(\gamma\cdot p-m)u_s(p)=0,

γpus(p)=mus(p).\gamma\cdot p\,u_s(p)=m u_s(p).

Thus

Λ+(p)us(p)=γp+m2mus(p)=us(p).\Lambda_+(p)u_s(p) ={\gamma\cdot p+m\over2m}u_s(p)=u_s(p).

On shell, p2=m2p^2=m^2, so

(γp+m)2=(γp)2+2mγp+m2=p2+2mγp+m2=2m(γp+m).(\gamma\cdot p+m)^2 =(\gamma\cdot p)^2+2m\gamma\cdot p+m^2 =p^2+2m\gamma\cdot p+m^2 =2m(\gamma\cdot p+m).

Therefore

Λ+2=γp+m2m=Λ+.\Lambda_+^2={\gamma\cdot p+m\over2m}=\Lambda_+.

Exercise 4: chirality versus helicity at zero mass

Section titled “Exercise 4: chirality versus helicity at zero mass”

For a massless positive-energy spinor, show that ψL\psi_L has helicity 1/2-1/2 and ψR\psi_R has helicity +1/2+1/2.

Solution

For m=0m=0,

(pσˉ)ψL=0,(pσ)ψR=0.(p\cdot\bar\sigma)\psi_L=0, \qquad (p\cdot\sigma)\psi_R=0.

With p0=pp^0=|\mathbf p|,

pσˉ=p(1+σp^),pσ=p(1σp^).p\cdot\bar\sigma=|\mathbf p|(1+\boldsymbol\sigma\cdot\widehat{\mathbf p}), \qquad p\cdot\sigma=|\mathbf p|(1-\boldsymbol\sigma\cdot\widehat{\mathbf p}).

Hence

σp^ψL=ψL,σp^ψR=+ψR.\boldsymbol\sigma\cdot\widehat{\mathbf p}\,\psi_L=-\psi_L, \qquad \boldsymbol\sigma\cdot\widehat{\mathbf p}\,\psi_R=+\psi_R.

The spin operator is S=σ/2\mathbf S=\boldsymbol\sigma/2, so the helicities are 1/2-1/2 and +1/2+1/2.

Use the Dirac equation and its adjoint to prove current conservation for jμ=ψˉγμψj^\mu=\bar\psi\gamma^\mu\psi.

Solution

The Dirac equation gives

γμμψ=imψ.\gamma^\mu\partial_\mu\psi=-im\psi.

The adjoint equation gives

(μψˉ)γμ=imψˉ.(\partial_\mu\bar\psi)\gamma^\mu=im\bar\psi.

Therefore

μjμ=(μψˉ)γμψ+ψˉγμμψ=imψˉψimψˉψ=0.\partial_\mu j^\mu =(\partial_\mu\bar\psi)\gamma^\mu\psi +\bar\psi\gamma^\mu\partial_\mu\psi =im\bar\psi\psi-im\bar\psi\psi=0.
  • Coleman, Sidney. Lectures on Quantum Field Theory. World Scientific, 2019, chapters 19–21.
  • Srednicki, Mark. Quantum Field Theory. Cambridge University Press, 2007, sections 34–38.
  • Weinberg, Steven. The Quantum Theory of Fields, Volume I: Foundations. Cambridge University Press, 1995, sections 5.4–5.6 and 14.1.
  • Zee, A. Quantum Field Theory in a Nutshell. 2nd ed., Princeton University Press, 2010, chapters II.1–II.2.