Skip to main content

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.
Starting model: ψ̄ and the sixteen bilinear covariantsMinkowski diagram of an electron's bilinears V^μ = ψ̄ γ^μ ψ (solid) and A^μ = ψ̄ γ^μ γ⁵ ψ (dashed), in units of ψ̄ψ: time parts up, z parts across, light cone faint. At rapidity 0.80, v = 0.664c, V rides its hyperbola to ψ†ψ = 1.337 while ψ̄ψ = 1 stays at the vertex; with spin 40° from z, A⁰ = 0.680.η = 0.80: v = 0.664c, spin 40° from ztime partz partψ̄ψ = 1Vμ solid, Aμ dashedψ†ψ = V⁰ = 1.337A⁰ = 0.680A¹ = 0.643, untouchedVμAμ = 0.000
Start with a picture. Change the model in the lab below.
1

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.

2

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.
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

0.80
0 to 1.5; steps of 0.05.
40 °
0 to 180 °; steps of 20.

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.

Interactive Advanced Quantum Mechanics modelMinkowski diagram of an electron's bilinears V^μ = ψ̄ γ^μ ψ (solid) and A^μ = ψ̄ γ^μ γ⁵ ψ (dashed), in units of ψ̄ψ: time parts up, z parts across, light cone faint. At rapidity 0.80, v = 0.664c, V rides its hyperbola to ψ†ψ = 1.337 while ψ̄ψ = 1 stays at the vertex; with spin 40° from z, A⁰ = 0.680.η = 0.80: v = 0.664c, spin 40° from ztime partz partψ̄ψ = 1Vμ solid, Aμ dashedψ†ψ = V⁰ = 1.337A⁰ = 0.680A¹ = 0.643, untouchedVμAμ = 0.000
Minkowski diagram of an electron's bilinears Vμ = ψ̄ γμ ψ (solid) and Aμ = ψ̄ γμ γ⁵ ψ (dashed), in units of ψ̄ψ: time parts up, z parts across, light cone faint. At rapidity 0.80, v = 0.664c, V rides its hyperbola to ψ†ψ = 1.337 while ψ̄ψ = 1 stays at the vertex; with spin 40° from z, A⁰ = 0.680.

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 °

Baseline 0.80; choose 0 to 1.5, in steps of 0.05.
4. Compared with the saved baseline, this result will…

This experiment is kept only while this page is open. It does not change saved course progress.

4

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.

5

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.

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.

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.

6

Quick check

Select an answer and explain it before checking.

At one point a Dirac spinor has the value ψ = (1, 0, 0.5i, 0) in the Dirac representation. A student finds ψ̄γ⁵ψ = i and decides she has made an error, since a bilinear built from the Hermitian matrix γ⁵ ought to be real. Which statement is correct?
7

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.

Extra challenging questions

Connect several ideas, explore supplied data, and compare your reasoning with a worked solution.