Spinor normalisation & completeness

Once the Dirac equation is solved, you hold four plane-wave spinors in your hand — two describing a spin-½ particle, two describing its antiparticle — and you must answer two questions before a single cross section can be computed. First, how large is each spinor? The equation fixes their direction in the four-dimensional spinor space but leaves an overall constant NN free, exactly as an eigenvalue equation fixes an eigenvector only up to scale. Second, how do the four together fill that space? A matrix element is never evaluated on one spin state in isolation; real experiments average over initial spins and sum over final spins, and that sum is only tractable if the spinors close into something simple. Both questions have the same answer, and it is dictated not by taste but by relativity. The normalisation is fixed to uˉu=2m\bar u u = 2m and u†u=2Eu^\dagger u = 2E — the spin-½ echo of the scalar field's "2E2E particles per unit volume" — so that spinor states dovetail with Lorentz-invariant phase space and the box volume cancels from every rate. The completeness relations then read ∑su(s)uˉ(s)=p̸+m\sum_s u^{(s)}\bar u^{(s)} = \slashed p + m and ∑sv(s)vˉ(s)=p̸−m\sum_s v^{(s)}\bar v^{(s)} = \slashed p - m: the spin sum collapses to a matrix built from the momentum alone. This single fact is what turns the forbidding-looking ∣M∣2‾\overline{|\mathcal M|^2} into a trace of gamma matrices — Casimir's trick — and is executed inside every QED and QCD calculation ever done. This article fixes NN, proves the two completeness relations, builds the energy projection operators Λ±\Lambda_\pm they generate, and shows why the whole machinery of spin-averaged scattering rests on them.

Key takeaways
  • The free-particle Dirac equation becomes algebraic, (p̸−m)u=0(\slashed p - m)u = 0 and (p̸+m)v=0(\slashed p + m)v = 0; it fixes each spinor's direction but leaves an overall constant NN free.
  • The relativistic convention sets N=E+mN = \sqrt{E+m}, giving uˉ(r)u(s)=2m δrs\bar u^{(r)} u^{(s)} = 2m\,\delta^{rs} (Lorentz scalar) and u(r)†u(s)=2E δrsu^{(r)\dagger} u^{(s)} = 2E\,\delta^{rs} (a density, the spin-½ version of "2E2E per unit volume").
  • Particle and antiparticle spinors are orthogonal, uˉ(r)v(s)=vˉ(r)u(s)=0\bar u^{(r)} v^{(s)} = \bar v^{(r)} u^{(s)} = 0, and the antiparticle scalar flips sign, vˉ(r)v(s)=−2m δrs\bar v^{(r)} v^{(s)} = -2m\,\delta^{rs}.
  • The completeness relations ∑su(s)uˉ(s)=p̸+m\sum_s u^{(s)}\bar u^{(s)} = \slashed p + m and ∑sv(s)vˉ(s)=p̸−m\sum_s v^{(s)}\bar v^{(s)} = \slashed p - m close the four spinors and generate the energy projection operators Λ±=(±p̸+m)/2m\Lambda_\pm = (\pm\slashed p + m)/2m.
  • These sums are what convert a spin-averaged ∣M∣2‾\overline{|\mathcal M|^2} into a trace of gamma matrices — Casimir's trick — the backbone of every unpolarised cross-section calculation.
  • Key idea: the Dirac equation fixes each spinor only up to a constant; setting it by uˉu=2m\bar u u = 2m (so u†u=2Eu^\dagger u = 2E) matches relativistic phase space, and the spin sums ∑suuˉ=p̸+m\sum_s u\bar u = \slashed p + m, ∑svvˉ=p̸−m\sum_s v\bar v = \slashed p - m let every spin-averaged matrix element be evaluated as a gamma-matrix trace.
The core insight

The Dirac equation fixes each spinor only up to an overall constant; choosing it so that uˉu=2m\bar u u = 2m and u†u=2Eu^\dagger u = 2E matches relativistic phase space, and the completeness relations ∑su(s)uˉ(s)=p̸+m\sum_s u^{(s)}\bar u^{(s)} = \slashed p + m, ∑sv(s)vˉ(s)=p̸−m\sum_s v^{(s)}\bar v^{(s)} = \slashed p - m turn every spin-averaged matrix element into a trace of gamma matrices.

Where this lives today

The completeness relations are not a historical curiosity — they run inside every calculation of a scattering rate. The moment an experiment at the LHC or a lepton collider does not measure the spin of its particles (which is nearly always), the theorist must average over initial spins and sum over final ones, and Casimir's trick — ∣M∣2‾=14∑spins∣M∣2=14Tr[⋯ ]\overline{|\mathcal M|^2} = \tfrac{1}{4}\sum_{\text{spins}}|\mathcal M|^2 = \tfrac14\mathrm{Tr}[\cdots] built from ∑uuˉ=p̸+m\sum u\bar u = \slashed p + m — is the only practical way to do it. Every textbook QED cross section (Møller, Bhabha, Compton, pair annihilation) and every parton-level QCD rate is a gamma-matrix trace generated this way, and symbolic engines such as FORM, FeynCalc, and the automated pipelines inside MadGraph, Sherpa, and WHIZARD implement exactly these spin sums and the Clifford-algebra trace identities to churn out ∣M∣2|\mathcal M|^2 for processes with dozens of external legs. When physics is spin-dependent, the same projection-operator technology carries over: polarised-beam programmes at a future ILC or FCC-ee insert the spin projector 12(1+γ5s̸)\tfrac12(1 + \gamma^5\slashed s) alongside Λ±\Lambda_\pm to compute helicity-dependent asymmetries, and lattice and effective-theory calculations use the massless limit p̸\slashed p and the chirality projectors 12(1±γ5)\tfrac12(1\pm\gamma^5) that Λ±\Lambda_\pm degenerate into. The tidy fact that four normalised spinors close into p̸±m\slashed p \pm m is, quietly, one of the most-executed pieces of mathematics in all of particle physics.

The physical picture

Picture the four solutions as four arrows spanning the four-dimensional space the Dirac matrices act on. In the particle's rest frame the picture is at its simplest: a positive-energy spinor is pure "upper block", u=2m (χ0)u = \sqrt{2m}\,\binom{\chi}{0}, and its two-component top χ\chi is nothing but the ordinary Pauli spin-up or spin-down state. Boost the particle and the lower block fills in, tied to the upper one by σ⋅p/(E+m)\boldsymbol\sigma\cdot\mathbf p/(E+m); the arrow tilts into the "small components", and its length as measured by u†uu^\dagger u grows from 2m2m to 2E2E. But one measure of the arrow does not change — the Lorentz-scalar combination uˉu=u†γ0u\bar u u = u^\dagger\gamma^0 u stays pinned at 2m2m in every frame. The first figure shows this split personality: a density that grows with energy, and an invariant that does not.

Figure
A rest-frame positive-energy spinor is pure upper block, u = √(2m)(χ, 0), with u†u = 2m. Under a boost to energy E the lower block fills in through σ·p ⁄ (E+m) and the density grows to u†u = 2E, yet the Lorentz scalar ū u stays fixed at 2m in every frame — the invariant that fixes the normalisation constant N = √(E+m).
The density u†u grows from 2m to 2E under a boost, but the scalar ū u stays 2m rest frame (E = m) χ↑ χ↓ 0 0 u†u = 2m ū u = 2m boost γ moving (E > m) χ↑ χ↓ σ·p E+m χ u†u = 2E ū u = 2m (unchanged)

The four spinors do more than sit there — together they close. Sum the two particle spinors over spin and you do not get the identity; you get a specific matrix, p̸+m\slashed p + m, that projects onto the positive-energy subspace and nothing else. The antiparticle pair sums to p̸−m\slashed p - m. These two completeness relations are the entire practical payoff of the chapter: they let an experimentalist who never measures a single spin still compute the cross section, because the messy sum over unobserved spin states collapses into a trace. The second figure makes the normalisation quantitative for a real electron — the invariant uˉu=2m\bar u u = 2m is a flat line, while the density u†u=2Eu^\dagger u = 2E climbs without bound.

Figure
The two spinor invariants for an electron (m = 0.511 MeV) versus momentum. The Lorentz scalar ū u = 2m = 1.022 MeV is flat, independent of frame; the density u†u = 2E = 2√(p² + m²) climbs without bound. At p = 0.885 MeV the electron has E = 2m, so u†u = 2.044 MeV while ū u is still 1.022 MeV.
The punchline: ū u = 2m is invariant, u†u = 2E grows with energy 0 2 4 6 8 10 0 1 2 3 4 5 electron momentum p (MeV) invariant value (MeV) ū u = 2m = 1.022 MeV u†u = 2E E = 2m: u†u = 2.044 MeV p = 0.885 MeV, γ = 2

The core equations

Substituting a plane wave ψ=u(p) e−ip⋅x\psi = u(p)\,e^{-ip\cdot x} into (i∂̸−m)ψ=0(i\slashed\partial - m)\psi = 0 sends every derivative to −ipμ-ip_\mu, converting the differential equation into a purely algebraic one for the constant spinor. The antiparticle wave ψ=v(p) e+ip⋅x\psi = v(p)\,e^{+ip\cdot x} carries the opposite sign in the exponent, flipping the mass term.

(p̸−m) u(p)=0,(p̸+m) v(p)=0,p̸≡γμpμ(1)(\slashed p - m)\,u(p) = 0,\qquad (\slashed p + m)\,v(p) = 0,\qquad \slashed p \equiv \gamma^\mu p_\mu \qquad (1)

where: u(p)u(p) — particle (positive-energy) spinor; v(p)v(p) — antiparticle spinor; p̸\slashed p — the momentum contracted with the gamma matrices; pμ=(E,p)p^\mu = (E,\mathbf p) on shell with E=+p2+m2E = +\sqrt{\mathbf p^2 + m^2}. Reading: this is an eigenvalue-type equation — it fixes the direction of uu in spinor space but not its length, leaving the constant to be set below. Solving equation (1) in the Dirac–Pauli representation, where γ0=diag(1,−1)\gamma^0 = \mathrm{diag}(\mathbb 1,-\mathbb 1), gives the explicit two-block forms. The upper and lower two-component blocks are locked together by σ⋅p\boldsymbol\sigma\cdot\mathbf p, and a single prefactor N=E+mN = \sqrt{E+m} sets the scale.

u(s)(p)=E+m(χsσ⋅pE+m χs),v(s)(p)=E+m(σ⋅pE+m ηsηs)(2)u^{(s)}(p) = \sqrt{E+m}\begin{pmatrix} \chi^{s} \\[2pt] \dfrac{\boldsymbol\sigma\cdot\mathbf p}{E+m}\,\chi^{s}\end{pmatrix},\qquad v^{(s)}(p) = \sqrt{E+m}\begin{pmatrix} \dfrac{\boldsymbol\sigma\cdot\mathbf p}{E+m}\,\eta^{s} \\[2pt] \eta^{s}\end{pmatrix} \qquad (2)

where: χs,ηs\chi^{s},\eta^{s} — two-component Pauli spinors (s=1,2s = 1,2), e.g. χ1=(10)\chi^1 = \binom{1}{0}, χ2=(01)\chi^2 = \binom{0}{1}; σ\boldsymbol\sigma — Pauli matrices; the lower block of uu (and upper block of vv) is the "small component", suppressed by ∣p∣/(E+m)|\mathbf p|/(E+m). Reading: at rest (p=0\mathbf p = 0) the small components vanish and u=2m (χ0)u = \sqrt{2m}\,\binom{\chi}{0} is a pure Pauli spinor. The normalisation is the whole point of the article. The overall constant is fixed so that the Lorentz-scalar bilinear equals 2m2m; the Hermitian density then automatically equals 2E2E, matching the scalar-field convention of the previous article.

uˉ(r)(p) u(s)(p)=2m δrs,u(r)†(p) u(s)(p)=2E δrs(3)\bar u^{(r)}(p)\,u^{(s)}(p) = 2m\,\delta^{rs},\qquad u^{(r)\dagger}(p)\,u^{(s)}(p) = 2E\,\delta^{rs} \qquad (3)

where: uˉ=u†γ0\bar u = u^\dagger\gamma^0 — the Dirac adjoint; δrs\delta^{rs} — Kronecker delta enforcing orthogonality of the two spin states. Reading: uˉu\bar u u is a Lorentz scalar and so is frame-independent at 2m2m; u†u=j0u^\dagger u = j^0 is a density and so scales like the energy — the two facts are the spin-½ mirror of ρ=2E\rho = 2E for the Klein–Gordon field. The antiparticle spinors obey the same 2E2E density but a sign-flipped scalar, and particle and antiparticle spinors at the same momentum are mutually orthogonal — the four together form a complete basis of the spinor space.

vˉ(r)v(s)=−2m δrs,v(r)†v(s)=2E δrs,uˉ(r)v(s)=vˉ(r)u(s)=0(4)\bar v^{(r)} v^{(s)} = -2m\,\delta^{rs},\quad v^{(r)\dagger} v^{(s)} = 2E\,\delta^{rs},\quad \bar u^{(r)} v^{(s)} = \bar v^{(r)} u^{(s)} = 0 \qquad (4)

where: the minus sign on vˉv\bar v v is not optional — it is forced by (p̸+m)v=0(\slashed p + m)v = 0. Reading: uˉu>0\bar u u > 0 and vˉv<0\bar v v < 0 is how the scalar bilinear distinguishes matter from antimatter; the vanishing cross terms guarantee the four solutions are linearly independent. The completeness relations are the practical engine. Summing the outer product uuˉu\bar u over the two spins gives not the identity but a matrix built purely from the momentum — the positive-energy content of p̸+m\slashed p + m.

∑s=1,2u(s)(p) uˉ(s)(p)=p̸+m,∑s=1,2v(s)(p) vˉ(s)(p)=p̸−m(5)\sum_{s=1,2} u^{(s)}(p)\,\bar u^{(s)}(p) = \slashed p + m,\qquad \sum_{s=1,2} v^{(s)}(p)\,\bar v^{(s)}(p) = \slashed p - m \qquad (5)

where: the left side is a 4×44\times4 matrix (outer product uuˉu\bar u), unlike the scalar uˉu\bar u u of equation (3); mm multiplies the 4×44\times4 identity. Reading: these are the relations that make spin-averaged cross sections computable — a sum over unmeasured spins becomes a matrix you can take the trace of. Dividing the spin sums by 2m2m produces the energy projection operators, which single out the particle or antiparticle subspace and satisfy the defining properties of projectors.

Λ±=±p̸+m2m,Λ++Λ−=1,Λ±2=Λ±,Λ+Λ−=0(6)\Lambda_\pm = \frac{\pm\slashed p + m}{2m},\qquad \Lambda_+ + \Lambda_- = \mathbb 1,\qquad \Lambda_\pm^2 = \Lambda_\pm,\qquad \Lambda_+\Lambda_- = 0 \qquad (6)

where: Λ+=12m∑su(s)uˉ(s)\Lambda_+ = \tfrac{1}{2m}\sum_s u^{(s)}\bar u^{(s)} projects onto positive-energy (particle) states; Λ−=−12m∑sv(s)vˉ(s)\Lambda_- = -\tfrac{1}{2m}\sum_s v^{(s)}\bar v^{(s)} projects onto antiparticle states. Reading: that Λ++Λ−=1\Lambda_+ + \Lambda_- = \mathbb 1 is the true "resolution of the identity" — the individual spin sum is only half of it.

Derivation AFixing the spinor normalisation constant

The goal is to prove that the single prefactor N=E+mN = \sqrt{E+m} in the explicit spinors of equation (2) is exactly what delivers both uˉu=2m\bar u u = 2m and u†u=2Eu^\dagger u = 2E, with the second following from the first with no extra freedom. We compute the two bilinears directly from the block form, using only the algebraic identity (σ⋅p)2=p2(\boldsymbol\sigma\cdot\mathbf p)^2 = \mathbf p^2 and the on-shell relation p2=E2−m2\mathbf p^2 = E^2 - m^2, and show that the same NN serves the scalar and the density at once — which is why the two conventions are one convention.

Step 1. Compute the Hermitian density. With u=E+m (χ, σ⋅pE+mχ)Tu = \sqrt{E+m}\,\big(\chi,\ \tfrac{\boldsymbol\sigma\cdot\mathbf p}{E+m}\chi\big)^{\mathsf T} and χ†χ=1\chi^\dagger\chi = 1, form u†u=(E+m)[χ†χ+χ†(σ⋅p)2(E+m)2χ]=(E+m)[1+p2(E+m)2]u^\dagger u = (E+m)\big[\chi^\dagger\chi + \chi^\dagger\tfrac{(\boldsymbol\sigma\cdot\mathbf p)^2}{(E+m)^2}\chi\big] = (E+m)\big[1 + \tfrac{\mathbf p^2}{(E+m)^2}\big]. Why: u†uu^\dagger u sums the squared magnitudes of all four components, and the Pauli identity (σ⋅p)2=p2 1(\boldsymbol\sigma\cdot\mathbf p)^2 = \mathbf p^2\,\mathbb 1 turns the lower block into a pure number.

Step 2. Collapse using the mass shell. Combine over a common denominator: u†u=(E+m)2+p2E+mu^\dagger u = \tfrac{(E+m)^2 + \mathbf p^2}{E+m}, and substitute p2=E2−m2\mathbf p^2 = E^2 - m^2 so the numerator becomes (E+m)2+E2−m2=2E2+2Em=2E(E+m)(E+m)^2 + E^2 - m^2 = 2E^2 + 2Em = 2E(E+m). Hence u†u=2Eu^\dagger u = 2E. Why: the on-shell condition is what converts an expression full of p\mathbf p into the clean 2E2E — off shell there would be no such simplification.

Step 3. Compute the scalar. Insert γ0\gamma^0, which flips the sign of the lower block: uˉu=u†γ0u=(E+m)[1−p2(E+m)2]=(E+m)2−p2E+m\bar u u = u^\dagger\gamma^0 u = (E+m)\big[1 - \tfrac{\mathbf p^2}{(E+m)^2}\big] = \tfrac{(E+m)^2 - \mathbf p^2}{E+m}. With p2=E2−m2\mathbf p^2 = E^2 - m^2 the numerator is 2Em+2m2=2m(E+m)2Em + 2m^2 = 2m(E+m), giving uˉu=2m\bar u u = 2m. Why: the single sign flip from γ0\gamma^0 is the whole difference between the frame-dependent density 2E2E and the invariant scalar 2m2m — same spinor, two different bilinears.

Step 4. Read off the convention's logic. Both results came from one prefactor N=E+mN = \sqrt{E+m}; had we instead demanded uˉu=1\bar u u = 1 we would take N=(E+m)/2mN = \sqrt{(E+m)/2m} and be forced to reinsert factors of 2E2E by hand in every phase-space integral. Repeating Steps 1–3 for vv with the blocks swapped gives v†v=2Ev^\dagger v = 2E and vˉv=−2m\bar v v = -2m, the sign flowing from (p̸+m)v=0(\slashed p + m)v = 0. Why: choosing uˉu=2m\bar u u = 2m is precisely the choice that lets the spinor normalisation and the 2E2E scalar-field normalisation share one Lorentz-invariant phase space.

Go deeper · Derivation B · 1 minThe completeness relation as a sum over spins
Derivation BThe completeness relation as a sum over spins

The aim is to establish ∑su(s)(p) uˉ(s)(p)=p̸+m\sum_s u^{(s)}(p)\,\bar u^{(s)}(p) = \slashed p + m — the outer-product spin sum, a 4×44\times4 matrix, not to be confused with the scalar inner product of the previous derivation. We build the sum from the explicit block spinors, use the completeness of the Pauli two-spinors ∑sχsχs†=1\sum_s \chi^s\chi^{s\dagger} = \mathbb 1, and reassemble the four blocks into the covariant object p̸+m\slashed p + m. A trace check at the end confirms the normalisation is internally consistent with uˉu=2m\bar u u = 2m.

Step 1. Write the outer product block by block. Since uˉ=u†γ0=(E+m)(χ†, −χ†σ⋅pE+m)\bar u = u^\dagger\gamma^0 = (E+m)\big(\chi^\dagger,\ -\chi^\dagger\tfrac{\boldsymbol\sigma\cdot\mathbf p}{E+m}\big), the product u(s)uˉ(s)u^{(s)}\bar u^{(s)} is a 2×22\times2 array of two-component blocks, each carrying a factor χsχs†\chi^s\chi^{s\dagger}. Why: organising the 4×44\times4 matrix into upper-left, upper-right, lower-left, lower-right blocks lets the Pauli-spinor sum act on each block independently.

Step 2. Sum over spin. Using ∑s=1,2χsχs†=1\sum_{s=1,2}\chi^s\chi^{s\dagger} = \mathbb 1, the four blocks become: upper-left (E+m)1(E+m)\mathbb 1; upper-right −σ⋅p-\boldsymbol\sigma\cdot\mathbf p; lower-left σ⋅p\boldsymbol\sigma\cdot\mathbf p; lower-right −(σ⋅p)2E+m=−p2E+m=−(E−m)-\tfrac{(\boldsymbol\sigma\cdot\mathbf p)^2}{E+m} = -\tfrac{\mathbf p^2}{E+m} = -(E-m). Why: the two Pauli spinors are complete in the two-dimensional spin space, so their sum erases all spin labels and leaves only momentum-dependent blocks.

Step 3. Recognise the gamma-matrix structure. Assembling the blocks, ∑suuˉ=((E+m)1−σ⋅pσ⋅p−(E−m)1)\sum_s u\bar u = \begin{pmatrix} (E+m)\mathbb 1 & -\boldsymbol\sigma\cdot\mathbf p \\ \boldsymbol\sigma\cdot\mathbf p & -(E-m)\mathbb 1\end{pmatrix}. In the Dirac representation γ0=(100−1)\gamma^0 = \big(\begin{smallmatrix}\mathbb 1&0\\0&-\mathbb 1\end{smallmatrix}\big) and γi=(0σi−σi0)\gamma^i = \big(\begin{smallmatrix}0&\sigma^i\\-\sigma^i&0\end{smallmatrix}\big), so Eγ0−p⋅γ+m1=p̸+mE\gamma^0 - \mathbf p\cdot\boldsymbol\gamma + m\mathbb 1 = \slashed p + m reproduces exactly this matrix. Hence ∑su(s)uˉ(s)=p̸+m\sum_s u^{(s)}\bar u^{(s)} = \slashed p + m. Why: the block matrix is not an accident — it is p̸+m\slashed p + m written out, so the spin sum is the covariant projector onto positive-energy states.

Step 4. Cross-check by trace and repeat for vv. Take the trace: Tr(p̸+m)=Tr(p̸)+4m=0+4m=4m\mathrm{Tr}(\slashed p + m) = \mathrm{Tr}(\slashed p) + 4m = 0 + 4m = 4m, while Tr∑suuˉ=∑suˉ(s)u(s)=∑s2m=4m\mathrm{Tr}\sum_s u\bar u = \sum_s \bar u^{(s)} u^{(s)} = \sum_s 2m = 4m — they agree, tying completeness back to the uˉu=2m\bar u u = 2m normalisation. The identical computation with the blocks of vv swapped gives ∑sv(s)vˉ(s)=p̸−m\sum_s v^{(s)}\bar v^{(s)} = \slashed p - m. Why: the trace of an outer-product sum equals the sum of the inner products, so consistency between equations (3) and (5) is guaranteed by a single number, 4m4m.

Go deeper · Derivation C · 1 minEnergy projection operators and Casimir's trick
Derivation CEnergy projection operators and Casimir's trick

The purpose here is to show that the spin sums, divided by 2m2m, are genuine projection operators onto the particle and antiparticle subspaces, and then to demonstrate the single most important consequence — that completeness collapses a spin-averaged squared matrix element into a trace of gamma matrices, the technique known as Casimir's trick that underlies every unpolarised cross-section calculation. We verify the projector algebra directly from p̸ p̸=p2=m2\slashed p\,\slashed p = p^2 = m^2 and then walk the sum-over-spins through a generic amplitude.

Step 1. Verify idempotency. With Λ+=(p̸+m)/2m\Lambda_+ = (\slashed p + m)/2m, compute Λ+2=(p̸+m)24m2=p̸ p̸+2mp̸+m24m2\Lambda_+^2 = \tfrac{(\slashed p + m)^2}{4m^2} = \tfrac{\slashed p\,\slashed p + 2m\slashed p + m^2}{4m^2}; using p̸ p̸=pμpνγμγν=p2=m2\slashed p\,\slashed p = p_\mu p_\nu\gamma^\mu\gamma^\nu = p^2 = m^2 on shell, the numerator is 2m2+2mp̸=2m(p̸+m)2m^2 + 2m\slashed p = 2m(\slashed p + m), so Λ+2=(p̸+m)/2m=Λ+\Lambda_+^2 = (\slashed p + m)/2m = \Lambda_+. Why: p̸ p̸=p2\slashed p\,\slashed p = p^2 follows from the Clifford algebra {γμ,γν}=2gμν\{\gamma^\mu,\gamma^\nu\} = 2g^{\mu\nu}, and it is what makes Λ+\Lambda_+ idempotent — the hallmark of a projector.

Step 2. Confirm the pair resolves the identity. Add the two: Λ++Λ−=(p̸+m)+(−p̸+m)2m=2m2m=1\Lambda_+ + \Lambda_- = \tfrac{(\slashed p + m) + (-\slashed p + m)}{2m} = \tfrac{2m}{2m} = \mathbb 1, and multiply: Λ+Λ−=(p̸+m)(−p̸+m)4m2=m2−p̸ p̸4m2=0\Lambda_+\Lambda_- = \tfrac{(\slashed p+m)(-\slashed p+m)}{4m^2} = \tfrac{m^2 - \slashed p\,\slashed p}{4m^2} = 0. Why: complementary, orthogonal, idempotent — the three properties say the four spinors split the space cleanly into a particle half and an antiparticle half, which is the precise meaning of "completeness".

Step 3. Set up a matrix element. A typical amplitude is M=uˉ(p′) Γ u(p)\mathcal M = \bar u(p')\,\Gamma\,u(p) for some string of gammas Γ\Gamma; its modulus squared is ∣M∣2=[uˉ(p′)Γu(p)] [uˉ(p′)Γu(p)]∗=uˉ(p′)Γu(p) uˉ(p)Γˉu(p′)|\mathcal M|^2 = [\bar u(p')\Gamma u(p)]\,[\bar u(p')\Gamma u(p)]^{*} = \bar u(p')\Gamma u(p)\,\bar u(p)\bar\Gamma u(p'), where Γˉ=γ0Γ†γ0\bar\Gamma = \gamma^0\Gamma^\dagger\gamma^0. Why: writing the complex conjugate as a second Dirac sandwich is what lets the two u(p)u(p) factors sit next to each other, ready to be summed.

Step 4. Average and sum over spins. Summing over the final spin and averaging over the two initial spins, ∣M∣2‾=12∑s,s′uˉs′(p′)Γus(p) uˉs(p)Γˉus′(p′)\overline{|\mathcal M|^2} = \tfrac12\sum_{s,s'}\bar u^{s'}(p')\Gamma u^{s}(p)\,\bar u^{s}(p)\bar\Gamma u^{s'}(p'); replacing ∑susuˉs=p̸+m\sum_s u^s\bar u^s = \slashed p + m and ∑s′us′uˉs′≠p′+m\sum_{s'} u^{s'}\bar u^{s'} = \slashed p' + m turns the spinor sandwiches into ∣M∣2‾=12 Tr[(̸p′+m) Γ (p̸+m) Γˉ]\overline{|\mathcal M|^2} = \tfrac12\,\mathrm{Tr}\big[(\slashed p' + m)\,\Gamma\,(\slashed p + m)\,\bar\Gamma\big]. Why: this is Casimir's trick — the completeness relations remove every explicit spinor, leaving a trace of gamma matrices that gamma-algebra identities (or a computer) evaluate mechanically. Every unpolarised QED and QCD rate is computed this way.

Limiting cases

CaseResultMeaning
Rest frame (p=0\mathbf p = 0)u=2m (χ0)u = \sqrt{2m}\,\binom{\chi}{0}, small components vanishspinor is a pure Pauli two-spinor; uˉu=u†u=2m\bar u u = u^\dagger u = 2m coincide
Non-relativistic (∣p∣≪m\lvert\mathbf p\rvert \ll m)lower block ≈σ⋅p2m χ\approx \tfrac{\boldsymbol\sigma\cdot\mathbf p}{2m}\,\chi, order v/cv/c"large" and "small" components; the Pauli equation and its g=2g = 2 moment emerge
Ultra-relativistic (E≫mE \gg m)upper and lower blocks equal in magnitudehelicity →\to chirality; uˉu=2m\bar u u = 2m negligible beside u†u=2Eu^\dagger u = 2E
Massless (m→0m \to 0)∑su(s)uˉ(s)→p̸\sum_s u^{(s)}\bar u^{(s)} \to \slashed p; uˉu→0\bar u u \to 0Weyl spinors; Λ±\Lambda_\pm degenerate and chirality projectors 12(1±γ5)\tfrac12(1\pm\gamma^5) take over
Spin sum vs identity∑suuˉ=p̸+m≠1\sum_s u\bar u = \slashed p + m \ne \mathbb 1the sum is 2mΛ+2m\Lambda_+, a projector onto particles, not the full completeness Λ++Λ−\Lambda_+ + \Lambda_-

Worked examples

1 · An explicit electron spinor and its two norms (fully worked)

Build the spin-up positive-energy spinor for an electron (m=0.511m = 0.511 MeV) moving along z^\hat z with energy E=2mE = 2m, and verify u†u=2Eu^\dagger u = 2E and uˉu=2m\bar u u = 2m numerically.

Step 1. Fix the kinematics. With E=2m=1.022E = 2m = 1.022 MeV, the momentum is ∣p∣=E2−m2=(2m)2−m2=3 m=0.885|\mathbf p| = \sqrt{E^2 - m^2} = \sqrt{(2m)^2 - m^2} = \sqrt3\,m = \mathbf{0.885} MeV, so γ=E/m=2\gamma = E/m = 2 and β=∣p∣/E=3/2=0.866\beta = |\mathbf p|/E = \sqrt3/2 = 0.866. The prefactor is N=E+m=3m=1.533=1.238 MeV1/2N = \sqrt{E+m} = \sqrt{3m} = \sqrt{1.533} = \mathbf{1.238}\ \mathrm{MeV}^{1/2}, and the small-component ratio is ∣p∣E+m=3m3m=13=0.577\tfrac{|\mathbf p|}{E+m} = \tfrac{\sqrt3 m}{3m} = \tfrac{1}{\sqrt3} = \mathbf{0.577}.

Step 2. Write the spinor. For spin up, χ=(10)\chi = \binom{1}{0}, and σ⋅p=∣p∣ σz\boldsymbol\sigma\cdot\mathbf p = |\mathbf p|\,\sigma_z acts on χ\chi to give ∣p∣(10)|\mathbf p|\binom{1}{0}. Thus

u(1)=E+m(10∣p∣E+m0)=1.238(100.5770)=(1.23800.7150) MeV1/2.u^{(1)} = \sqrt{E+m}\begin{pmatrix} 1 \\ 0 \\ \tfrac{|\mathbf p|}{E+m} \\ 0\end{pmatrix} = 1.238\begin{pmatrix} 1 \\ 0 \\ 0.577 \\ 0\end{pmatrix} = \begin{pmatrix} 1.238 \\ 0 \\ 0.715 \\ 0\end{pmatrix}\ \mathrm{MeV}^{1/2}.

The third entry 0.715=1.238×0.5770.715 = 1.238\times0.577 is the "small component" — here not so small, because the electron is already relativistic.

Step 3. Verify both norms. The density is u†u=1.2382 (1+0.5772)=1.533×(1+0.333)=1.533×1.333=2.044u^\dagger u = 1.238^2\,(1 + 0.577^2) = 1.533\times(1 + 0.333) = 1.533\times1.333 = \mathbf{2.044} MeV, which equals 2E=2×1.022=2.0442E = 2\times1.022 = 2.044 MeV. The scalar inserts γ0\gamma^0, flipping the sign of the lower block: uˉu=1.533×(1−0.333)=1.533×0.667=1.022\bar u u = 1.533\times(1 - 0.333) = 1.533\times0.667 = \mathbf{1.022} MeV, which equals 2m=2×0.511=1.0222m = 2\times0.511 = 1.022 MeV. Both check exactly: u†u=2E=2.044u^\dagger u = 2E = 2.044 MeV (frame-dependent) and uˉu=2m=1.022\bar u u = 2m = 1.022 MeV (invariant), the two points read off the punchline figure at p=0.885p = 0.885 MeV.

2 · Completeness two ways: the spin-summed ∣uˉu∣2\lvert\bar u u\rvert^2

Evaluate S=∑r,s=1,2∣uˉ(r)(p) u(s)(p)∣2\displaystyle S = \sum_{r,s=1,2}\big|\bar u^{(r)}(p)\,u^{(s)}(p)\big|^2 first from orthonormality and then from the completeness relation, and confirm the two agree.

Step 1. Direct route via orthonormality. Equation (3) gives uˉ(r)u(s)=2m δrs\bar u^{(r)} u^{(s)} = 2m\,\delta^{rs}, so only the diagonal terms r=sr = s survive: S=∑r=s(2m)2=2×(2m)2=8m2S = \sum_{r=s}(2m)^2 = 2\times(2m)^2 = \mathbf{8m^2}. The factor of two counts the two spin states. This uses nothing but the normalisation.

Step 2. Trace route via completeness. Write ∣uˉrus∣2=uˉrus uˉsur|\bar u^{r} u^{s}|^2 = \bar u^{r} u^{s}\,\bar u^{s} u^{r} and sum: S=∑r,suˉrusuˉsur=Tr[(∑susuˉs)(∑ruruˉr)]=Tr[(p̸+m)2]S = \sum_{r,s}\bar u^{r} u^{s}\bar u^{s} u^{r} = \mathrm{Tr}\big[\big(\sum_s u^s\bar u^s\big)\big(\sum_r u^r\bar u^r\big)\big] = \mathrm{Tr}\big[(\slashed p + m)^2\big]. Expand: (p̸+m)2=p̸ p̸+2mp̸+m2=m2+2mp̸+m2=2m2+2mp̸(\slashed p + m)^2 = \slashed p\,\slashed p + 2m\slashed p + m^2 = m^2 + 2m\slashed p + m^2 = 2m^2 + 2m\slashed p. Take the trace, using Tr 1=4\mathrm{Tr}\,\mathbb 1 = 4 and Tr p̸=0\mathrm{Tr}\,\slashed p = 0: S=Tr(2m21)+2m Tr(p̸)=8m2+0=8m2S = \mathrm{Tr}(2m^2\mathbb 1) + 2m\,\mathrm{Tr}(\slashed p) = 8m^2 + 0 = \mathbf{8m^2}.

Step 3. Compare and read the lesson. Both routes give S=8m2S = 8m^2; for the electron, 8m2=8×0.5112=2.09 MeV28m^2 = 8\times0.511^2 = \mathbf{2.09\ MeV^2}. The agreement is the whole point of Casimir's trick in miniature: the completeness relation reproduces, as a mechanical gamma-matrix trace, an answer you could also get by hand-summing orthonormal spinors — and for any realistic operator sandwiched between the spinors, only the trace route stays tractable.

Physics of the portrait

The spinor, its four solutions, and the negative-energy pair that became antimatter are one person's legacy. The completeness relation ∑svvˉ=p̸−m\sum_s v\bar v = \slashed p - m is the modern, tidy descendant of an idea that once seemed like a disaster: an infinite sea of filled negative-energy states.

Origins · 1 minPaul Dirac and how this was discovered
PD
Paul Dirac
1902–1984 · the relativistic electron and the birth of antimatter

Paul Adrien Maurice Dirac was an English theoretical physicist of famously few words and unmatched formal instinct, who in 1928 set out to find a relativistic wave equation that was first order in time — repairing the negative-probability defect of the Klein–Gordon equation. Demanding that his operator square to E2=p2+m2E^2 = \mathbf p^2 + m^2 forced the coefficients to be anticommuting matrices, the gamma matrices, and forced the wavefunction to have four components: the Dirac spinor. Two of its four plane-wave solutions had positive energy and carried spin up and down — the electron, with the correct magnetic moment g=2g = 2 falling out for free, a triumph nobody had engineered. But two solutions stubbornly had negative energy. Rather than discard them, Dirac proposed in 1930 that the vacuum is a filled sea of negative-energy states, and that a hole in that sea would appear as a positive-energy particle of opposite charge. In 1931 he predicted this "anti-electron", and in 1932 Carl Anderson found the positron in cosmic-ray cloud-chamber tracks — one of the great predictive triumphs of theoretical physics, and the reason the vv-spinors of this article carry the label "antiparticle". The Dirac-sea picture is now superseded by quantum field theory, where the negative-energy solutions are simply antiparticle creation operators and the completeness relation ∑sv(s)vˉ(s)=p̸−m\sum_s v^{(s)}\bar v^{(s)} = \slashed p - m encodes their spin content cleanly, but the physics Dirac extracted from an apparent embarrassment — that every charged fermion has an antiparticle — remains exact. He shared the 1933 Nobel Prize with Schrödinger, and held Newton's Lucasian chair at Cambridge; his credo, that a physical law must possess mathematical beauty, is nowhere better vindicated than in the equation that bears his name.

History
1928
— Dirac writes his first-order relativistic wave equation and its four-component spinor solutions
1930
— Dirac proposes the filled negative-energy sea and the hole interpretation
1931
— Dirac predicts an anti-electron of equal mass and opposite charge
1932
— Anderson discovers the positron in cosmic-ray cloud-chamber tracks
1949
— Feynman and Dyson's covariant perturbation theory makes the spin sums u ū = p-slash + m and Casimir's trace trick routine

Going further

Texts. Thomson, Modern Particle Physics, §4.6–4.7 (solutions of the Dirac equation, spin sums and completeness) · Peskin & Schroeder, An Introduction to Quantum Field Theory, §3.3 (Dirac-field spinors, uˉu=2m\bar u u = 2m, ∑uuˉ=p̸+m\sum u\bar u = \slashed p + m) · Halzen & Martin, Quarks and Leptons, §5.3–5.4 (normalisation and projection operators) · Griffiths, Introduction to Elementary Particles, 2nd ed., §7.3 and Casimir's trick in §7.7 · Schwartz, Quantum Field Theory and the Standard Model, §11.1 · Bjorken & Drell, Relativistic Quantum Mechanics, §3.2 (energy projection operators). The key results to keep in mind: N=E+mN = \sqrt{E+m} gives u†u=2Eu^\dagger u = 2E and uˉu=2m\bar u u = 2m, with the electron numbers E=2m=1.022E = 2m = 1.022 MeV, ∣p∣=0.885|\mathbf p| = 0.885 MeV, small component 0.5770.577; the spin sum ∑suuˉ=p̸+m\sum_s u\bar u = \slashed p + m holds block by block with trace 4m4m; the projectors satisfy Λ±2=Λ±\Lambda_\pm^2 = \Lambda_\pm, Λ++Λ−=1\Lambda_+ + \Lambda_- = \mathbb 1, and Λ+Λ−=0\Lambda_+\Lambda_- = 0; and the spin-summed check gives ∑r,s∣uˉrus∣2=8m2=2.09 MeV2\sum_{r,s}|\bar u^r u^s|^2 = 8m^2 = 2.09\ \mathrm{MeV}^2, recoverable both directly from orthonormality and by the gamma-matrix trace route.

Prerequisites. Dirac equation · Gamma matrices and the Clifford algebra · Dirac solutions · Spin from the Dirac equation. Leads to. Spin and helicity states · Antiparticles · Matrix elements and spin averaging · Dirac adjoint and bilinear covariants.

Ask the tutor: Ask me to run Casimir's trick end to end on a real process — say electron–muon scattering — and show exactly where ∑suuˉ=p̸+m\sum_s u\bar u = \slashed p + m turns the spin sum into the trace that gives the cross section.

Why normalise to ū u = 2m rather than the tidier ū u = 1?

Because the spinor normalisation cannot be chosen in isolation — it has to mesh with the way the whole scattering machinery counts particles. The relativistic convention sets "2E2E particles per unit volume" for scalar fields, which is what makes Lorentz-invariant phase space and the flux factor come out clean, and u†u=2Eu^\dagger u = 2E is the spin-½ version of exactly that density. Its Lorentz-scalar partner is then forced to be uˉu=2m\bar u u = 2m. If you instead normalise to uˉu=1\bar u u = 1 (as some older texts do), nothing is wrong, but you must divide each external leg's phase-space factor by 2E2E by hand to compensate — you have only moved the 2E2E elsewhere. Keeping uˉu=2m\bar u u = 2m lets one invariant phase space serve particles of every spin, which is why Peskin, Thomson, Halzen–Martin, and Schwartz all adopt it.

What is the difference between ū u and u† u?

They are different bilinears of the same spinor and they transform differently. The combination u†uu^\dagger u is the time component of the probability current jμ=uˉγμuj^\mu = \bar u\gamma^\mu u, i.e. a density; it is not Lorentz invariant, and it equals 2E2E, growing with the particle's energy exactly as a density must under length contraction. The combination uˉu=u†γ0u\bar u u = u^\dagger\gamma^0 u inserts a γ0\gamma^0 and is a genuine Lorentz scalar; it equals 2m2m in every frame. The single factor of γ0\gamma^0 — which flips the sign of the lower two components — is the whole difference between the frame-dependent 2E2E and the invariant 2m2m. In figures and formulas, whenever a quantity must be the same to all observers (like a mass term in a Lagrangian) it is built from uˉu\bar u u; whenever it is a counting density it is built from u†uu^\dagger u.

Why does the spin sum give a matrix, but the normalisation gives a number?

Because they multiply the spinors in opposite orders. The normalisation uˉu\bar u u is an inner product: a row spinor uˉ\bar u times a column spinor uu contracts all four components into a single number, 2m2m. The completeness sum ∑su(s)uˉ(s)\sum_s u^{(s)}\bar u^{(s)} is an outer product: a column times a row produces a full 4×44\times4 matrix, and summing the two spin states gives p̸+m\slashed p + m. Both use the same spinors; one contracts them, the other spreads them into a matrix. The inner product answers "how big is this spinor?", the outer product answers "what subspace do these spinors span?" — which is why the first is a normalisation and the second is completeness.

What do the u and v spinors mean physically?

The uu spinors are the two positive-energy solutions and describe a spin-½ particle (say, an electron) with spin up or down; the vv spinors are the two solutions attached to the e+ip⋅xe^{+ip\cdot x} plane wave and describe its antiparticle (a positron). Historically the vv spinors are Dirac's negative-energy solutions reinterpreted: rather than a filled sea with a hole, modern quantum field theory treats them as antiparticle states outright. The tell-tale sign is the scalar bilinear — uˉu=+2m\bar u u = +2m for particles, vˉv=−2m\bar v v = -2m for antiparticles — so the sign of ψˉψ\bar\psi\psi literally distinguishes matter from antimatter. Both carry the same positive density u†u=v†v=2Eu^\dagger u = v^\dagger v = 2E, because both describe real, positive-energy quanta once the theory is quantised.

How exactly does completeness make a cross section computable?

Through Casimir's trick. A matrix element is a Dirac sandwich uˉ(p′)Γu(p)\bar u(p')\Gamma u(p); its square, summed over the spins nobody measured, contains ∑su(s)(p)uˉ(s)(p)\sum_s u^{(s)}(p)\bar u^{(s)}(p) and ∑s′u(s′)(p′)uˉ(s′)(p′)\sum_{s'} u^{(s')}(p')\bar u^{(s')}(p'). Replacing each sum by p̸+m\slashed p + m (or ̸p′+m\slashed p' + m) eliminates every explicit spinor and leaves ∣M∣2‾=14Tr[(̸p′+m)Γ(p̸+m)Γˉ]\overline{|\mathcal M|^2} = \tfrac14\mathrm{Tr}[(\slashed p' + m)\Gamma(\slashed p + m)\bar\Gamma] — a trace of gamma matrices. Traces of gamma-matrix strings are evaluated by a short list of Clifford-algebra identities (trace of an odd number vanishes; trace of two is 4gμν4g^{\mu\nu}; and so on), so what looked like an intractable sum over four or eight spin combinations becomes a mechanical algebraic exercise. Without completeness you would have to build every spinor explicitly and sum by brute force; with it, the spins are gone before you start.

Is the choice of the Pauli spinors χ¹, χ² physical?

No — it is a choice of basis, exactly like choosing the zz-axis to label ordinary spin. The Dirac equation determines a two-dimensional space of positive-energy solutions at each momentum, and any orthonormal pair of two-component spinors χ1,χ2\chi^1,\chi^2 picks a basis for it. The conventional choice χ1=(10)\chi^1 = \binom{1}{0}, χ2=(01)\chi^2 = \binom{0}{1} labels spin along zz, but a helicity basis (spin along p\mathbf p) or any rotated pair is equally valid; a general physical state is a linear combination. What is not a choice is the completeness relation: ∑su(s)uˉ(s)=p̸+m\sum_s u^{(s)}\bar u^{(s)} = \slashed p + m holds for any orthonormal basis, because the sum over a complete basis is basis-independent. So the individual spinors are conventional, but everything physical — cross sections, the projector p̸+m\slashed p + m — is not.

Why does the antiparticle scalar ū u come out negative, and is that a problem?

It is not a problem; it is a feature that labels antimatter. The sign is forced by the defining equation: uu satisfies (p̸−m)u=0(\slashed p - m)u = 0 while vv satisfies (p̸+m)v=0(\slashed p + m)v = 0, and running the normalisation computation with that opposite-sign mass term flips the scalar to vˉv=−2m\bar v v = -2m. Nothing negative happens to a probability — the density v†v=2Ev^\dagger v = 2E is perfectly positive, so there is always a positive number of antiparticles per unit volume. Only the Lorentz-scalar bilinear changes sign, and that sign is precisely what makes ψˉψ\bar\psi\psi a useful order parameter: a mass term mψˉψm\bar\psi\psi in the Lagrangian treats particles and antiparticles symmetrically because the sign flip is compensated by the way vv enters the field expansion. The negative vˉv\bar v v is the harmless bookkeeping trace of Dirac's negative-energy solutions, now understood as antiparticles.

Before this

Where this leads