
Advanced Quantum Mechanics course sequence · Interactive lesson
Why four components: the Dirac representation constructed
This is the step where ψ stops being one function. Two lines of trace algebra and one fact about Pauli matrices force four components; then we build the two 4 × 4 representations every later calculation uses, Dirac's for slow particles and Weyl's for fast ones, and see why the choice changes nothing measurable.
UAQM-04.02 · Lesson 2 of 9 in this unit. Position is not a completion score.Build the model
Four anticommuting matrices need four components; two run out of room.
The key idea
The Dirac representation is the choice β = [[1, 0], [0, −1]] and αᵢ = [[0, σᵢ], [σᵢ, 0]] in 2 × 2 blocks: four Hermitian 4 × 4 matrices that square to 1 and anticommute in pairs, the smallest matrices that can.
Example. At rest the Dirac Hamiltonian is βmc², already diagonal: for an electron its entries are +0.511, +0.511, −0.511 and −0.511 MeV, so the trace is zero and the upper pair of components sits 1.022 MeV above the lower pair.
The complete physical model and assumptions
The previous topic turned Dirac's requirements into algebra: four Hermitian matrices αₓ, αy, αz and β, each squaring to 1, anticommuting in pairs, with nothing yet saying how large they are. Two lines decide. Multiply αᵢβ = −βαᵢ by β to get βαᵢβ = −αᵢ, and the trace, which ignores cyclic reordering, then says Tr αᵢ = −Tr αᵢ, so Tr αᵢ = 0; the same trick empties the trace of β. A Hermitian matrix that squares to 1 has eigenvalues ±1, so a zero trace means as many eigenvalues +1 as −1, and the dimension is even. The even candidate 2 then fails by geometry. A traceless Hermitian 2 × 2 matrix is a⋅σ for a real vector a, two such matrices anticommute exactly when their vectors are perpendicular, and three-dimensional space holds three perpendicular unit vectors, not four. So four is the minimum, and β = [[1, 0], [0, −1]] with αᵢ = [[0, σᵢ], [σᵢ, 0]] reaches it: the doubled space supplies the missing direction. The Weyl representation, β off the diagonal, is the same answer seen through a unitary matrix, and Pauli's theorem, quoted, says every 4 × 4 answer is. The argument stops at counting. It fixes how many components ψ has, not what they mean: the mass term is what demands the fourth matrix, the vanishing trace of H forces two negative-energy states at every momentum in every representation, and their meaning waits for Unit 8.
1One anticommuting partner makes a matrix traceless
The previous topic turned Dirac's requirements into algebra: αₓ, αy, αz and β are Hermitian, each squares to the N × N identity, and any two anticommute. Nothing yet fixes N. Take αᵢ and multiply αᵢβ = −βαᵢ on the left by β; since β² = 1, βαᵢβ = −αᵢ. Now take the trace. A trace is unchanged when the last factor of a product is moved to the front, so Tr(βαᵢβ) = Tr(β²αᵢ) = Tr αᵢ, while βαᵢβ = −αᵢ says the same trace is −Tr αᵢ. A number equal to its own negative is zero: Tr αᵢ = 0, derived. Running the same steps with any αᵢ as the partner gives αᵢβαᵢ = −β and Tr β = 0. Only two facts were used, that the partner squares to 1 and anticommutes, and every one of the four matrices has three such partners. The trace of a finite matrix is the whole engine; it is the same argument that shows [x, p] = iħ has no finite matrices, since the left side would be traceless and the right side would have trace iħN.
2Eigenvalues ±1 in equal numbers make the dimension even
A Hermitian matrix is diagonalised by a unitary and has real eigenvalues, and M² = 1 makes each one obey λ² = 1, so λ = +1 or −1. With n₊ of the first and n₋ of the second, the basis-independent trace is Tr M = n₊ − n₋, which the previous step set to zero: n₊ = n₋ and N = 2n₊, derived. One component is out, since a traceless 1 × 1 matrix is the number 0, and so is three, whose trace can only be ±1 or ±3. The same count applies to the Hamiltonian. Work in SI units, H = cα⋅p + βmc², keeping the metric diag(1, −1, −1, −1) for the four-vectors that arrive next topic; both terms carry joules, cp as (m s⁻¹)(kg m s⁻¹) and mc² as kg m² s⁻², because α and β are pure numbers. H is a sum of traceless matrices, so Tr H = 0, and the Klein-Gordon requirement gave H² = (p²c² + m²c⁴)1. So H/Eₚ, with Eₚ = √(p²c² + m²c⁴), is Hermitian, squares to 1 and is traceless: half the energies at every momentum are +Eₚ and half are −Eₚ, in any representation. No choice of matrices removes the negative branch; Unit 8 interprets it.
3Two components run out of directions
Evenness still allows N = 2, so test it. The identity and the three Pauli matrices span every 2 × 2 Hermitian matrix with real coefficients, a₀1 + a⋅σ, and tracelessness removes a₀. Quantum Mechanics I's identity σᵢσⱼ = δᵢⱼ1 + iεᵢⱼₖσₖ gives (a⋅σ)(b⋅σ) = (a⋅b)1 + i(a × b)⋅σ, and the cross-product term flips sign when the order is reversed, so the anticommutator is 2(a⋅b)1. Dirac's conditions become geometry: squaring to 1 means a unit vector, anticommuting means perpendicular vectors. Three mutually perpendicular unit vectors fit in three dimensions, which is why σₓ, σy and σz exist, but a fourth would have to be perpendicular to all three, hence zero, and the zero vector squares to 0, not 1. Two components hold at most three of Dirac's four matrices; with N even and N ≠ 2, N ≥ 4, derived. The limiting case shows what the fourth matrix pays for. At m = 0, β leaves H, three matrices suffice, the σᵢ serve, and H = ±cσ⋅p is a two-component equation, Weyl's. It is the mass term that costs two extra components.
4Doubling the space supplies the missing direction
Write 4 × 4 matrices in 2 × 2 blocks and try β = [[1, 0], [0, −1]], αᵢ = [[0, σᵢ], [σᵢ, 0]]. Each is Hermitian because 1 and σᵢ are. Block multiplication checks the algebra: αᵢαⱼ = [[σᵢσⱼ, 0], [0, σᵢσⱼ]], so αᵢαⱼ + αⱼαᵢ = 2δᵢⱼ times the 4 × 4 identity; αᵢβ = [[0, −σᵢ], [σᵢ, 0]] and βαᵢ = [[0, σᵢ], [−σᵢ, 0]] sum to zero; β² = 1 on sight, and all four traces vanish as they must. Why it works shows best if the block index gets Pauli matrices τ of its own: αᵢ = τₓ ⊗ σᵢ and β = τz ⊗ 1. The σ factor makes the α's anticommute among themselves, and τz, a Pauli matrix of the block index whose vector is perpendicular to that of τₓ, is the fourth direction the old space lacked. β is diagonal, so at rest H = βmc² = diag(mc², mc², −mc², −mc²). For the electron, with measured mₑc² = 0.51100 MeV, the upper pair sits at +0.511 MeV and the lower pair at −0.511 MeV, 1.022 MeV apart. That diagonal rest energy is why this is the representation for slow particles, and the one Unit 6 uses for the non-relativistic limit.
5The Weyl representation is one unitary away
Nothing in the algebra prefers a diagonal β. If U†U = 1, then UαᵢU† and UβU† are still Hermitian, still square to 1 and still anticommute, since U†U = 1 lets products pass through. Take U = (1/√2)[[1, −1], [1, 1]] in blocks. Multiplying out, UβU† = [[0, 1], [1, 0]] and UαᵢU† = [[−σᵢ, 0], [0, σᵢ]], derived: the Weyl, or chiral, representation. Wavefunctions follow, ψW = UψD, so Dirac halves φ and χ become ψL = (φ − χ)/√2 and ψR = (φ + χ)/√2, with ψW†ψW = ψD†ψD. Here H = [[−cσ⋅p, mc²], [mc², cσ⋅p]], and the mass is the only thing linking the halves. At m = 0 the equation falls apart into two Weyl equations, the two-component limit found above, which is why this representation suits fast or massless particles. The measured record fits that limit, quoted: in 1958 Goldhaber, Grodzins and Sunyar found the neutrino from electron capture in ¹⁵²ᵐEu to be left-handed, as one Weyl half alone describes, and tritium-decay data from KATRIN now bound the neutrino mass below 0.45 eV, so any mass linking its halves is more than a million times smaller than the electron's.
6Pauli's theorem: every 4 × 4 answer is the same answer
Could some other 4 × 4 set be genuinely different? Pauli's fundamental theorem, quoted here and not proved, says no: any two sets of 4 × 4 matrices obeying the algebra are related by M′ = SMS⁻¹ with one invertible S, unitary when both sets are Hermitian, and larger solutions are built from copies of the 4 × 4 one, so N is a multiple of four (also quoted). Measurement cannot tell the representations apart: a similarity transformation keeps every eigenvalue of H, and with ψ′ = Uψ every density ψ†ψ and expectation value ψ†Aψ is unchanged. Components are not physics. An electron at rest with spin up is (1, 0, 0, 0) in the Dirac representation and (1, 0, 1, 0)/√2 in the Weyl one, half its probability in each half. The freedom remains even with the α's held fixed: β = τy ⊗ 1 also squares to 1 and anticommutes with every τₓ ⊗ σᵢ, and the unitary e(iπτₓ/4) ⊗ 1, a quarter turn of the block space about τₓ, carries τz ⊗ 1 into it while leaving each αᵢ alone. Pick the representation that shortens the calculation; nothing measured depends on the choice.
Key equations
Read the units, assumptions and limits before substituting.
One anticommuting partner empties the trace
Tr M = Tr(MB²) = Tr(BMB) = −Tr M, so Tr M = 0
Derived from the cyclic property of the trace alone. Every one of the four matrices has such a partner, so all four are traceless.
M, B any two of αₓ, αy, αz, β: dimensionless N × N matrices with B² = 1 and MB = −BMEqual numbers of +1 and −1
Tr M = n₊ − n₋ = 0, so N = n₊ + n₋ = 2n₊
Derived: the dimension is even. That rules out 1 and 3 but not 2, which needs the geometry of the next card.
n₊, n₋ count eigenvalues +1 and −1, the only ones since M† = M and M² = 1; N is the number of components of ψTwo-by-two matrices are vectors
(a⋅σ)(b⋅σ) + (b⋅σ)(a⋅σ) = 2(a⋅b)1
Derived: squaring to 1 means |a| = 1 and anticommuting means a ⊥ b. Space has room for three such vectors, Dirac needs four, so N ≥ 4.
a, b real dimensionless 3-vectors; every traceless Hermitian 2 × 2 matrix is a⋅σ; from σᵢσⱼ = δᵢⱼ1 + iεᵢⱼₖσₖThe Dirac representation
β = [[1, 0], [0, −1]], αᵢ = [[0, σᵢ], [σᵢ, 0]]
Built and checked by block multiplication. β is diagonal, so at rest H = βmc² puts the upper pair at +mc² and the lower pair at −mc².
2 × 2 blocks built from the identity 1 and the Pauli matrices σᵢ, so each matrix is 4 × 4, Hermitian and dimensionlessThe Weyl representation, one unitary away
U = (1/√2)[[1, −1], [1, 1]]: UβU† = [[0, 1], [1, 0]], UαᵢU† = [[−σᵢ, 0], [0, σᵢ]]
Derived by block multiplication. Here H = [[−cσ⋅p, mc²], [mc², cσ⋅p]], so at m = 0 it splits into two-component Weyl equations.
U is 4 × 4 in 2 × 2 blocks with U†U = 1; ψW = UψD turns the Dirac halves φ, χ into (φ − χ)/√2 and (φ + χ)/√2Negative energies at every momentum
Tr H = 0 and H² = Eₚ²1, so H has eigenvalues +Eₚ and −Eₚ, N/2 of each
Derived: H/Eₚ is Hermitian, squares to 1 and is traceless. This holds in every representation; Unit 8 says what the negative branch means.
H = cα⋅p + βmc² and Eₚ = √(p²c² + m²c⁴) in J; p in kg m s⁻¹; m in kg; c in m s⁻¹Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. Swing θ and φ: the three anticommutator bars trade height, but their squares always sum to |b|², so one stays at least |b|/√3. Shrinking |b| lowers them only by pushing β² − 1 towards −1, and M bottoms out at 0.75 near |b| = 0.71, never at zero.
Live results
- LARGEST |bᵢ|
- 0.579
- SUM OF bᵢ²
- 1.000
- β² − 1
- 0.000
- MISMATCH M
- 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. Polar angle θ of b: 55 °; Azimuth φ of b: 45 °; Length |b|: 1.00
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
The four components are not four particles or four spin states, and the upper pair is not "the electron" with the lower pair "the positron". The number four comes from needing four mutually anticommuting matrices; which components count as upper is a choice of representation. An electron at rest does sit entirely in the upper pair of the Dirac representation, but the same state has exactly half its probability in each half of the Weyl representation, and even in the Dirac representation a moving electron has lower components: at pc = 0.800 MeV they carry 23% of the probability. What every representation shares is the spectrum, two eigenvalues at +Eₚ and two at −Eₚ at each momentum. That is a statement about solutions, not about components, and the meaning of the negative pair is Unit 8's to give, not something read off the layout of the matrices.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
A classmate proposes three-component Dirac matrices and, failing that, a two-component set with αᵢ = σᵢ and β = (σₓ + σy + σz)/√3. Use the trace argument to rule out N = 3, then show that the proposed β squares to 1 and is traceless yet fails, by computing βσₓ + σₓβ.
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.
- Every Dirac matrix M is Hermitian with M² = 1, so its eigenvalues are ±1, and βαᵢβ = −αᵢ with the cyclic trace makes Tr M = 0. For N = 3, Tr M = n₊ − n₋ with n₊ + n₋ = 3 can only be ±1 or ±3, never 0, so no 3 × 3 set exists.
- The proposed β is b⋅σ with b = (1, 1, 1)/√3, so |b|² = 3 × (1/3) = 1. By (a⋅σ)(b⋅σ) = (a⋅b)1 + i(a × b)⋅σ, β² = |b|²1 = 1, and Tr β = 0 because every σᵢ is traceless: it passes both tests the trace argument sets.
- The anticommutator with αₓ = σₓ, whose vector is x̂, is βσₓ + σₓβ = 2(b⋅x̂)1 = (2/√3)1 = 1.155 × 1, not 0. By symmetry σy and σz give the same 1.155.
- Check by components: βσₓ = (1/√3)(1 − iσz + iσy) and σₓβ = (1/√3)(1 + iσz − iσy), using σyσₓ = −iσz and σzσₓ = iσy. The imaginary parts cancel in the sum, leaving (2/√3)1.
- No better choice exists: the three residuals 2bₓ, 2by, 2bz have squares summing to 4|b|² = 4, so the most even split leaves each at 2/√3 = 1.155, which this β already achieves.
Answer. N = 3 is impossible because a traceless matrix with eigenvalues ±1 needs an even dimension; the proposed β squares to 1 and is traceless, but each anticommutator βσᵢ + σᵢβ is 1.155 times the identity instead of 0.
MediumExample 2
An electron moves along z with pc = 0.800 MeV. In the Dirac representation, confirm by block multiplication that αz and β anticommute, write H = cpαz + βmc² as a 4 × 4 matrix, and find its eigenvalues using mₑc² = 0.511 MeV. Check the trace and the limit p → 0.
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.
- αzβ = [[0, σz], [σz, 0]][[1, 0], [0, −1]] = [[0, −σz], [σz, 0]] and βαz = [[0, σz], [−σz, 0]]. Their sum is the zero matrix, and αz² = [[σz², 0], [0, σz²]] = 1.
- H = [[mc², cpσz], [cpσz, −mc²]]. Since σz = diag(1, −1), H links only component 1 with 3 and 2 with 4: in MeV the blocks are [[0.511, 0.800], [0.800, −0.511]] on (ψ₁, ψ₃) and [[0.511, −0.800], [−0.800, −0.511]] on (ψ₂, ψ₄).
- Each block is traceless with determinant −(0.511² + 0.800²) = −(0.2611 + 0.6400) = −0.9011 MeV², so its eigenvalues are ±√0.9011 MeV = ±0.9493 MeV.
- The eigenvalues of H are +0.949, +0.949, −0.949 and −0.949 MeV. They sum to 0 = Tr H = mc² + mc² − mc² − mc², as tracelessness requires: two positive, two negative.
- Checks: every entry of each block is an energy in MeV. As p → 0 the eigenvalues go to ±0.511 MeV, the rest-energy pair, and when pc ≫ mc² they approach ±pc.
Answer. E = ±0.949 MeV, each twice: the positive pair is √(p²c² + m²c⁴), and the traceless H forces the matching negative pair.
HardExample 3
For the same electron (pc = 0.800 MeV along z, mₑc² = 0.511 MeV), take the positive-energy eigenvector of H in the (ψ₁, ψ₃) block of the Dirac representation. Transform it to the Weyl representation with U = (1/√2)[[1, −1], [1, 1]], find the probability in each half in both representations, and show that ⟨β⟩ is the same in both. Check the limits p → 0 and p → ∞.
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.
- In the block [[mc², pc], [pc, −mc²]] the eigenvector for E = +0.9493 MeV is proportional to (E + mc², pc) = (1.4603, 0.8000). With r = pc/(E + mc²) = 0.5478, the normalised spinor is (1, 0, r, 0)/√(1 + r²), and 1 + r² = 1.3001.
- Dirac weights: the upper half carries 1/(1 + r²) = 0.769 and the lower half r²/(1 + r²) = 0.231.
- Weyl halves: ψL = (φ − χ)/√2 = (1 − r, 0)/√(2(1 + r²)) and ψR = (φ + χ)/√2 = (1 + r, 0)/√(2(1 + r²)). Their weights are 0.2044/2.6003 = 0.079 and 2.3958/2.6003 = 0.921, summing to 1.000 as U†U = 1 demands.
- Dirac: ⟨β⟩ = 0.769 − 0.231 = 0.538. Weyl, with β = [[0, 1], [1, 0]]: ⟨β⟩ = 2 Re(ψL†ψR) = (1 − r²)/(1 + r²) = 0.6999/1.3001 = 0.538. Algebra turns (1 − r²)/(1 + r²) into mc²/E = 0.511/0.9493 = 0.538, a pure number as it must be.
- Limits: as p → 0, r → 0, and the Dirac weights go to 1 and 0 while the Weyl weights go to 0.5 and 0.5. As p → ∞, r → 1, and the Dirac weights go to 0.5 each while the Weyl weight collects in ψR, a right-handed state, since spin and p both point along +z. Each representation gathers the state into one half, in opposite limits.
Answer. Dirac halves 0.769 and 0.231; Weyl halves 0.079 (ψL) and 0.921 (ψR); ⟨β⟩ = mc²/E = 0.538 in both, so the components depend on the representation and the expectation value does not.
Quick check
Select an answer and explain it before checking.
Review & continue
The Dirac representation is the choice β = [[1, 0], [0, −1]] and αᵢ = [[0, σᵢ], [σᵢ, 0]] in 2 × 2 blocks: four Hermitian 4 × 4 matrices that square to 1 and anticommute in pairs, the smallest matrices that can.
- 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.