
Advanced Quantum Mechanics course sequence · Interactive lesson
ψ̄ and the sixteen bilinear covariants
Unit 4 used ψ†ψ as a probability density without asking whether a moving observer would agree. This is where you find out which densities a Dirac field can carry, why ψ̄ and not ψ† has to stand on the left, and how to tell at sight whether a bilinear is a scalar, a four-vector or a tensor before it is ever coupled to a field.
UAQM-05.05 · Lesson 5 of 9 in this unit. Position is not a completion score.Build the model
A spinor is never measured; the sixteen densities built from ψ̄ and ψ are.
The key idea
A bilinear covariant is a density ψ̄Γψ, with ψ̄ = ψ†γ⁰ and Γ one of sixteen matrices, 1, γ⁵, γμ (four), γμ γ⁵ (four) or σμν (six), chosen so that it transforms under Lorentz transformations as a scalar, pseudoscalar, four-vector, axial vector or antisymmetric tensor.
Example. An electron moving at 0.60c along z, with ψ̄ψ = 1 in units of its rest density, keeps ψ̄ψ = 1 while ψ̄ γμ ψ = (1.25, 0, 0, 0.75): the time component ψ†ψ has grown by γ = 1.25, and 1.25² − 0.75² = 1.00 has not changed.
The complete physical model and assumptions
A Dirac spinor is never measured; densities quadratic in it are, and covariance decides which of them mean the same thing in every frame. The obvious density fails: boosts act on ψ through matrices that are not unitary, so ψ†ψ changes from frame to frame. The repair is the adjoint ψ̄ = ψ†γ⁰. Because S⁻¹ = γ⁰S†γ⁰ (topic 5.03), ψ̄ transforms with S⁻¹, and every bilinear obeys ψ̄′Γψ′ = ψ̄(S⁻¹ΓS)ψ, so sorting densities becomes sorting 4 × 4 matrices under conjugation by S. Γ = 1 gives a scalar and Γ = γμ a four-vector, by the covariance condition alone: Unit 4's current is a four-vector, ψ†ψ is its time component, and conservation makes ∫ψ†ψ d³x, but not ∫ψ̄ψ d³x, the same in every frame. Products of gammas collapse to sixteen, 1 and γ⁵, four γμ and four γμ γ⁵, and six σμν; each squares to ±1, and any two different ones multiply to a phase times a member other than 1, which is traceless, so they are orthogonal under the trace and span every 4 × 4 matrix: there is no seventeenth kind of bilinear. Products transform one Λ per index, which makes σμν a tensor, and γ⁵ picks up det Λ, equal to 1 for every transformation reached from the identity. There the classification stops. Rotations and boosts treat 1 and γ⁵ identically, and γμ and γμ γ⁵ identically, so pseudoscalar and axial vector are, for now, names awaiting an improper transformation: parity, in topic 5.06.
1Why ψ† is the wrong partner and ψ̄ the right one
Units throughout: SI, with xμ = (ct, x), metric g = diag(1, −1, −1, −1) and ψ in m(−3/2), so every bilinear is a density in m⁻³. Quoted from topic 4.03 in one line each: the Clifford algebra, γμ γν + γν γμ = 2gμν; and the trace identities, tr 1 = 4, the trace of an odd number of gammas is zero, tr(γμ γν) = 4gμν and tr γ⁵ = 0. A density must be quadratic in ψ, and the first guess, ψ†Γψ, becomes ψ†(S†ΓS)ψ. For a rotation S is unitary, but a boost is a Hermitian stretch with S†S = cosh η + α⋅n̂ sinh η (topic 5.04), and S†ΓS bears no clean relation to Γ. Topic 5.03 found what replaces unitarity, γ⁰S†γ⁰ = S⁻¹. Put γ⁰ beside ψ†: ψ̄′ = ψ†S†γ⁰ = ψ†γ⁰(γ⁰S†γ⁰) = ψ̄S⁻¹, derived. Every density then obeys ψ̄′Γψ′ = ψ̄(S⁻¹ΓS)ψ, so the matrix alone decides how it transforms. Γ = 1 is a scalar at once, but not a probability: the rest spinor (1, 0, 0, 0) has ψ̄ψ = ψ†ψ, while the negative-energy rest spinor (0, 0, 1, 0) has ψ̄ψ = −ψ†ψ.
2The current is a four-vector, and its charge is invariant
Take Γ = γμ. The covariance condition S⁻¹ γμ S = Λμν γν turns the sandwich rule into ψ̄′ γμ ψ′ = Λμν ψ̄ γν ψ, so jμ = cψ̄ γμ ψ is a four-vector, derived. Since γ⁰γ⁰ = 1 and γ⁰γⁱ = αᵢ, its parts are j⁰ = cψ†ψ and j = cψ†αψ, Unit 4's density and current. So ψ†ψ is a time component and must grow under a boost: a plane wave at 0.60c has ψ†ψ = 1.25 ψ̄ψ, and |j|/j⁰ = 0.60 is its speed over c, so j = ρv holds in every frame, as it did non-relativistically. Now integrate. With ∂μ jμ = 0 (Unit 4), take the four-volume between the slice t = 0 of one frame and the slice t′ = 0 of another: the divergence integrates to zero, the sides at spatial infinity carry nothing, so equal flux crosses both slices, and because jμ is a four-vector the flux through t′ = 0 is ∫ψ′†ψ′ d³x′. So ∫ψ†ψ d³x is the same in every frame, derived: density up by γ, length down by γ. Measured: helium and H₂ are neutral to about 10⁻²⁰ e (King, 1960), while if charge followed ∫ψ̄ψ d³x, which for a bound electron falls below 1 by about its kinetic energy over mc², helium's two electrons, with 79.0 eV of kinetic energy between them, would by themselves come about 1.5 × 10⁻⁴ e short.
3Every product of gammas is one of sixteen
Any product of gammas can be reordered with γμ γν = −γν γμ for μ ≠ ν and shortened with (γμ)² = gμμ = ±1, no sum, until it is ±1 or ±i times a product of distinct gammas in increasing order. Such products are labelled by the subsets of the four indices 0, 1, 2, 3, so there are 2⁴ = 16: the empty product, four with one factor, six with two, four with three and one with four. Renamed, the two-factor products are −iσμν; the four-factor product is γ⁰γ¹γ²γ³ = −iγ⁵, and each three-factor product is γμ γ⁵ up to a phase, for example γ¹γ²γ³ = −iγ⁰γ⁵. Every one squares to ±1: (γ⁵)², (σ¹²)² and (γ¹γ⁵)² are 1, while (γ¹)², (σ⁰¹)² and (γ⁰γ⁵)² are −1. And every one except 1 is traceless by the quoted identities: γμ and γμ γ⁵ hold an odd number of gammas, σμν with μ ≠ ν has trace proportional to gμν = 0, and tr γ⁵ = 0. Splitting a two-factor product into its symmetric and antisymmetric halves gives γμ γν = gμν − iσμν: a number plus a tensor, with the i/2 in σμν chosen so that σⁱʲ = εᵢⱼₖΣₖ is Hermitian with eigenvalues ±1.
4Orthogonal under the trace, therefore a basis
The product of two different members is labelled by the symmetric difference of two different subsets, which is never empty, so ΓA ΓB is, whenever A ≠ B, a phase times a member other than 1, with zero trace. With ΓA = (ΓA)⁻¹ = ±ΓA this gives tr(ΓA ΓB) = 4δAB, derived. If Σ cA ΓA = 0, multiplying by ΓB and tracing gives 4cB = 0, so the sixteen are independent, and sixteen independent matrices in the sixteen-dimensional space of 4 × 4 matrices are a basis: M = ¼Σ tr(M ΓA) ΓA for every M, derived. Apply it to M = ψψ̄, whose trace against ΓA is ψ̄ ΓA ψ: ψψ̄ = ¼[ψ̄ψ + (ψ̄γ⁵ψ)γ⁵ + Vμ γμ − Aμ γμ γ⁵ + ½Tμν σμν]. The sixteen bilinears are the coordinates of one matrix. For (1, 0, 0, 0), ψψ̄ = diag(1, 0, 0, 0) = ¼(1 + γ⁰ + γ³γ⁵ + σ¹²), so ψ̄ψ = V⁰ = A³ = T¹² = 1 and the other twelve vanish. And a matrix commuting with all four γμ commutes with all sixteen, so with every matrix: it is a multiple of 1, the uniqueness topics 5.02 and 5.06 lean on.
5Products transform one Λ per index, and γ⁵ carries det Λ
Insert SS⁻¹ between neighbouring gammas: S⁻¹ γμ γν S = (S⁻¹ γμ S)(S⁻¹ γν S) = Λμα Λνβ γα γβ. A product of gammas thus transforms with one Λ per index, and σμν is an antisymmetric tensor, derived. Write γ⁵ = −(i/4!)εμνρσ γμ γν γρ γσ with ε₀₁₂₃ = −1: four factors of Λ contracted into ε give det Λ times ε, the definition of a determinant, so S⁻¹γ⁵S = (det Λ)γ⁵ for any Λ that has an S, derived. For rotations and boosts det Λ = 1, as a second route confirms: γ⁵ anticommutes with every γμ and so commutes with every σμν and every S. Hence iψ̄γ⁵ψ is invariant, ψ̄ γμ γ⁵ ψ a four-vector and ψ̄ σμν ψ a tensor. Quoting topic 5.06's SP = γ⁰ for ΛP = diag(1, −1, −1, −1), the same law makes iψ̄γ⁵ψ odd and sends the axial vector to (−A⁰, A) while the vector goes to (V⁰, −V). Reality fixes the last factors: (ψ̄Γψ)* = ψ̄(γ⁰Γ†γ⁰)ψ, and γ⁰Γ†γ⁰ = Γ for fifteen of the sixteen but −γ⁵ for γ⁵, so ψ̄γ⁵ψ alone is imaginary: hence the i in iψ̄γ⁵ψ.
6Reading the sixteen for a moving electron
For a free spinor with ūu = 2mc, all sixteen follow from pμ and the spin four-vector sμ: iūγ⁵u = 0; ū γμ u = 2pμ; ū γμ γ⁵ u = 2mc sμ; and ū σμν u = −2εμνρσ pρ sσ, each derived by checking it at rest, because two tensors equal in one frame are equal in all. Put the spin along z and boost along z to 0.60c, so γ = 1.25 and γβ = 0.75. In units of 2mc, V goes from (1, 0, 0, 0) to (1.25, 0, 0, 0.75); A goes from (0, 0, 0, 1) to (0.75, 0, 0, 1.25), V's mirror image in the light cone, its time part γβ now marking positive helicity; and the tensor, at rest only the spin density T¹² = 1, is untouched, as a boost along z leaves Bz untouched. A transverse spin fares the other way, keeping A and growing a σ⁰ⁱ part, as the second worked example shows. The invariants survive: VμVμ = 1.25² − 0.75² = 1, AμAμ = 0.75² − 1.25² = −1, VμAμ = 1.25 × 0.75 − 0.75 × 1.25 = 0 and ½Tμν Tμν = (T¹²)² = 1. For any spinor the Fierz identities, quoted, give VμVμ = (ψ̄ψ)² + (iψ̄γ⁵ψ)², AμAμ = −VμVμ and VμAμ = 0. What the laws cannot do is tell 1 from γ⁵, or γμ from γμ γ⁵: each pair transforms identically under every rotation and boost. Only a transformation with det Λ = −1 can, and parity, topic 5.06, is the one that does.
Key equations
Read the units, assumptions and limits before substituting.
The adjoint transforms with S⁻¹
ψ̄′(x′) = ψ̄(x) S⁻¹, where ψ̄ = ψ†γ⁰
Derived: ψ†S†γ⁰ = ψ†γ⁰(γ⁰S†γ⁰). Every density then obeys ψ̄′Γψ′ = ψ̄(S⁻¹ΓS)ψ, so sorting densities means sorting the matrices Γ under S⁻¹ΓS.
ψ in m(−3/2); ψ′(x′) = Sψ(x) with x′ = Λx, S a dimensionless 4 × 4 matrix; uses S⁻¹ = γ⁰S†γ⁰ from topic 5.03.Unit 4's current is a four-vector
j′μ(x′) = Λμν jν(x), with jμ = cψ̄ γμ ψ = (cψ†ψ, cψ†αψ)
Derived from S⁻¹ γμ S = Λμν γν. So ψ†ψ is a time component and grows by γ under a boost, while ψ̄ψ = ψ†ψ/γ for a free positive-energy plane wave.
jμ in m⁻² s⁻¹; Λμν dimensionless; xμ = (ct, x) with metric diag(1, −1, −1, −1).Total probability is frame-independent
Q = ∫ψ†ψ d³x = ∫ψ′†ψ′ d³x′
Derived: a conserved four-vector has one flux through every spacelike slice. ∫ψ̄ψ d³x is not invariant: for a free positive-energy packet it is Q/γ, because d³x contracts.
Q dimensionless; needs ∂μ jμ = 0 (Unit 4) and ψ → 0 at spatial infinity; each integral over a slice of constant time.Sixteen matrices, orthogonal under the trace
tr(ΓA ΓB) = 4δAB, with ΓA = (ΓA)⁻¹ = ±ΓA
Derived: for A ≠ B, ΓA ΓB is a phase times a member other than 1, hence traceless. So the sixteen are independent, and any 4 × 4 matrix is M = ¼Σ tr(M ΓA) ΓA.
ΓA: 1, γ⁵ = iγ⁰γ¹γ²γ³, four γμ and four γμ γ⁵, six σμν = (i/2)(γμ γν − γν γμ); all dimensionless.γ⁵ carries the determinant of Λ
S⁻¹γ⁵S = (det Λ) γ⁵
Derived from the covariance condition and the definition of a determinant. Under proper Λ, iψ̄γ⁵ψ is invariant and ψ̄ γμ γ⁵ ψ is a four-vector.
γ⁵ = −(i/4!)εμνρσ γμ γν γρ γσ with ε₀₁₂₃ = −1; det Λ = +1 for rotations and boosts, −1 for parity.The tensor bilinear
ψ̄′ σμν ψ′ = Λμα Λνβ ψ̄ σαβ ψ
Derived by inserting SS⁻¹ between the two gammas. A boost mixes σ⁰ⁱ with σⁱʲ as it mixes E with B, so a moving spin density acquires σ⁰ⁱ parts wherever the spin lies across the motion.
Six components, each a density in m⁻³; in the Dirac representation σⁱʲ = εᵢⱼₖΣₖ and σ⁰ⁱ = iαᵢ; Λ dimensionless.Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. Raise η and watch V climb its hyperbola: ψ†ψ grows as cosh η while the vertex, ψ̄ψ = 1, never moves. Then turn the spin from 0° to 180°: A slides along V's mirror image in the light cone, its time part changes sign with the helicity, and A¹ ignores the boost.
Live results
- DENSITY ψ†ψ ÷ ψ̄ψ
- 1.337
- SPEED v/c
- 0.664
- AXIAL TIME PART A⁰
- 0.680 × ψ̄ψ
- INVARIANT AμAμ
- -1.000 × (ψ̄ψ)²
Results follow this lesson’s model and assumptions. Displayed values are rounded.
Predict → change one variable → check
Plan a change before running it. Every other input stays at the saved baseline.
Saved baseline. Boost rapidity η along z: 0.80; Rest-frame spin angle from z: 40 °
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
Because ψ̄ψ is the Lorentz scalar, students conclude that the relativistic normalisation must be ∫ψ̄ψ d³x = 1, and that ∫ψ†ψ d³x = 1 is a non-relativistic leftover. It is the other way round. A scalar density integrated over a slice of constant time is not a scalar, because the volume of that slice contracts: for a free positive-energy packet, ψ̄ψ keeps its rest value while the packet shortens by γ, so ∫ψ̄ψ d³x = Q/γ, only 0.80 of its rest value at 0.60c. What is frame-independent is the flux of a conserved four-vector current, and that is ∫ψ†ψ d³x. The scalar fails as a probability a second way: ψ̄ψ = φ†φ − χ†χ goes negative whenever the lower components dominate, as in the negative-energy rest spinor. ∫ψ̄ψ d³x does mean something, just not probability: for Q = 1 in an energy eigenstate, the Hellmann–Feynman theorem makes it ∂E/∂(mc²), the response of the energy to the mass.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
An electron packet is a uniform cylinder of rest length 4.0 nm along z, with ψ†ψ = ψ̄ψ = ρ₀ inside it in its rest frame and ∫ψ†ψ d³x = 1. It moves along z with kinetic energy 0.400 MeV (mc² = 0.511 MeV). Find ψ†ψ, ψ̄ψ and the packet's length in the lab, then ∫ψ†ψ d³x and ∫ψ̄ψ d³x there.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- γ = 1 + 0.400/0.511 = 1.783, so β = √(1 − 1/γ²) = 0.828 and γβ = 1.476. At rest jμ = (cρ₀, 0, 0, 0), and the lower components are negligible, which is why ψ̄ψ = ψ†ψ there.
- jμ is a four-vector, so in the lab jμ = cρ₀(γ, 0, 0, γβ): ψ†ψ = j⁰/c = 1.783ρ₀, and the current travels at c j³/j⁰ = 0.828c, with the packet, as the non-relativistic j = ρv requires.
- ψ̄ψ is a scalar, so at each point of the packet it keeps its rest value ρ₀. The packet's length contracts to 4.0 nm/1.783 = 2.24 nm; its cross-section is unchanged.
- ∫ψ†ψ d³x = 1.783ρ₀ × (cross-section × 2.24 nm) = ρ₀ × cross-section × 4.0 nm = 1, as conservation demands. ∫ψ̄ψ d³x = ρ₀ × cross-section × 2.24 nm = 1/γ = 0.561.
- Checks: both integrals are pure numbers, since ψ carries m(−3/2); as the kinetic energy tends to 0, γ → 1 and both tend to 1; and 0.561 = mc²/E = 0.511/0.911, the free-particle ∂E/∂(mc²), which is what ∫ψ̄ψ d³x measures.
Answer. ψ†ψ = 1.78ρ₀ and ψ̄ψ = ρ₀ over a length of 2.24 nm; ∫ψ†ψ d³x = 1 in both frames, while ∫ψ̄ψ d³x falls to 1/γ = 0.561.
MediumExample 2
An electron at rest has spin along +x: ψ = (1, 1, 0, 0)/√2, so ψ̄ψ = 1. Expand ψψ̄ in the sixteen matrices by traces and read off all sixteen bilinears. Then boost along z to 0.80c, find the vector, axial-vector and tensor bilinears from the transformation laws alone, check one component against the boosted spinor, and check the invariants.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- ψ̄ = ψ†γ⁰ is the row (1, 1, 0, 0)/√2, since γ⁰ flips only lower entries, so ψψ̄ has ½ in each of its four upper-left entries and zero elsewhere. Its coefficient on ΓA is ¼tr(ψψ̄ ΓA) = ¼ψ̄ ΓA ψ: each trace is a bilinear.
- Only four traces survive: ψψ̄ = ¼(1 + γ⁰ + γ¹γ⁵ + σ²³), using Σₓ = σ²³ and γ⁰Σₓ = γ¹γ⁵. Matching ¼[ψ̄ψ + Vμ γμ − Aμ γμ γ⁵ + ½Tμν σμν + …] gives ψ̄ψ = 1, V⁰ = 1, A¹ = 1 (A₁ = −1 with the index down) and T²³ = 1; the other twelve, iψ̄γ⁵ψ among them, vanish.
- Boost: η = artanh 0.80 = ln 3, cosh η = 5/3, sinh η = 4/3. ψ̄ψ = 1 and iψ̄γ⁵ψ = 0 stay. Vμ = (5/3, 0, 0, 4/3); Aμ = (0, 1, 0, 0) is untouched, being transverse. The tensor law gives T²³ = cosh η = 5/3 and T⁰² = Λ⁰₃ Λ²₂ T³² = (4/3)(−1) = −4/3.
- Direct check: SB = cosh(η/2) + αz sinh(η/2) = (2 + αz)/√3, and αz sends (1, 1, 0, 0) to (0, 0, 1, −1), so ψ′ = (2, 2, 1, −1)/√6. Then ψ′†ψ′ = 10/6 = 5/3 = V⁰, ψ̄′ψ′ = (4 + 4 − 1 − 1)/6 = 1, and T⁰² = iψ′†γ²ψ′ = i(4i/3) = −4/3, as the law said.
- Invariants: VμVμ = 25/9 − 16/9 = 1 = (ψ̄ψ)², AμAμ = −1, VμAμ = 0, and ½Tμν Tμν = (T²³)² − (T⁰²)² = 1. Limit: as the speed tends to 0, T⁰² → 0 and only the rest-frame spin density T²³ = 1 remains.
Answer. At rest ψψ̄ = ¼(1 + γ⁰ + γ¹γ⁵ + σ²³). At 0.80c: ψ̄ψ = 1, V = (1.667, 0, 0, 1.333), A = (0, 1, 0, 0), T²³ = 1.667 and T⁰² = −1.333, with VμVμ = 1, AμAμ = −1 and VμAμ = 0.
HardExample 3
Hydrogen-like uranium, Z = 92, is in its ground state, normalised to ∫ψ†ψ d³x = 1. Taking the point-nucleus Dirac energy E = mc²√(1 − (Zα)²) as quoted from Unit 7, find ∫ψ̄ψ d³x and the probability in the lower components, compare with the non-relativistic estimate, and compare the Dirac binding energy with experiment, where the measured ground-state Lamb shift is 460.2 ± 4.6 eV.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- H = cα⋅p + βmc² − Zαħc/r, and only βmc² depends on m. The Hellmann–Feynman theorem, used and not re-derived, gives ∂E/∂(mc²) = ∫ψ†βψ d³x = ∫ψ̄ψ d³x, since β = γ⁰. With E proportional to mc², ∂E/∂(mc²) = E/mc² = √(1 − (Zα)²): derived, given the quoted energy.
- Zα = 92/137.036 = 0.67136 and (Zα)² = 0.45072, so ∫ψ̄ψ d³x = √0.54928 = 0.74113, while ∫ψ†ψ d³x = 1. The scalar integral is a quarter short; the charge is not.
- With ψ = (φ, χ), ψ†ψ = φ†φ + χ†χ and ψ̄ψ = φ†φ − χ†χ, so ∫χ†χ d³x = (1 − 0.74113)/2 = 0.1294. The non-relativistic χ ≈ (σ⋅p/2mc)φ gives ⟨p²⟩/(2mc)² = (Zα)²/4 = 0.1127, since ⟨p²⟩ = (Zαmc)² in the ground state: 13% low at this Z.
- Binding: mc²(1 − 0.74113) = 511.00 keV × 0.25887 = 132.28 keV. Experiment, Dirac's value less the measured Lamb shift, gives 131.82 keV: the binding is right to 0.35%, and E/mc² = 0.74204 against 0.74113, to 0.12%. Finite nuclear size and QED supply the difference (quoted).
- Checks: both integrals are pure numbers; for hydrogen, Z = 1, ∫ψ̄ψ d³x = 0.999973 and ∫χ†χ d³x = 1.331 × 10⁻⁵ = α²/4 to four figures, the non-relativistic limit; and as Zα → 1, ∫ψ̄ψ d³x → 0 while the charge stays e.
Answer. ∫ψ̄ψ d³x = 0.741 against ∫ψ†ψ d³x = 1; the lower components hold 0.129, against 0.113 non-relativistically; the Dirac binding of 132.28 keV is 0.35% above the measured 131.82 keV.
Quick check
Select an answer and explain it before checking.
Review & continue
A bilinear covariant is a density ψ̄Γψ, with ψ̄ = ψ†γ⁰ and Γ one of sixteen matrices, 1, γ⁵, γμ (four), γμ γ⁵ (four) or σμν (six), chosen so that it transforms under Lorentz transformations as a scalar, pseudoscalar, four-vector, axial vector or antisymmetric tensor.
- Explain the relationship without looking at the equation.
- Check a sign, a unit and a model assumption in your working.
- Return to the lab or worked example if a step is still unclear.