Skip to main content

Advanced Quantum Mechanics course sequence · Interactive lesson

Dirac field with anticommutators: positrons without a sea

This is the topic where the Dirac sea stops being needed. Keep topic 9.05's modes, swap the commutator for an anticommutator and follow it through: positive energies, positrons as quanta, the Pauli principle and causal observables all come out, and the sea survives only as a constant that normal ordering drops.

UAQM-09.07 · Lesson 7 of 9 in this unit. Position is not a completion score.
Starting model: Dirac field with anticommutators: positrons without a seaOccupation of one mode of the anticommuting Dirac field against its quantum's energy Eₚ, at k_B T = 0.30 MeV and electron chemical potential μ = 0.20 MeV. Solid: electrons, 1/(eˣ + 1) with x = (Eₚ − μ)/k_B T; dashed: positrons, with Eₚ + μ. No state lies below mc² (shaded); no curve passes 1. At rest: 0.262 electrons and 0.085 positrons per mode.solid e⁻: x = (Eₚ − μ)/kB T, at rest 0.262dashed e⁺: x = (Eₚ + μ)/kB T, at rest 0.085⟨n⟩ = 1/(eˣ + 1) ≤ 1, as (b†)² = 010mc²123energy of the quantum Eₚ / MeV
Start with a picture. Change the model in the lab below.
1

Build the model

Change the bracket's symmetry and Dirac's sea shrinks to a constant you may drop.

The key idea

Quantising the Dirac field with anticommutators means imposing (ψₐ(t, x), ψ†b(t, y)) = δab δ³(x − y) instead of a commutator; it gives (b, b†) = (d, d†) = 1, so electrons and positrons are both quanta of positive energy and no mode can hold two.

Example. A positron at rest is d†|0⟩ with p = 0. Its norm is ⟨0|d d†|0⟩ = 1 − ⟨0|d†d|0⟩ = +1, its energy +mc² = +0.511 MeV and its charge +1.602 × 10⁻¹⁹ C, while d†d†|0⟩ = 0, so no second positron can join it in that mode.

The complete physical model and assumptions

Topic 9.05 fixed the Dirac field's dynamics and left its algebra open, and topic 9.06 showed that commutators fail: the term −Eₚ d d† in the Hamiltonian then leaves either negative norms or energies with no floor. The model here changes one thing. The canonical bracket keeps its right-hand side, iħδab δ³(x − y) against the conjugate momentum iħψ†, but becomes an anticommutator. Because (d, d†) = (d†, d), the operator that multiplies the negative-frequency wave can be declared a creator at no cost in norm, and the reordering d d† = 1 − d†d turns −Eₚ d d† into +Eₚ d†d plus a constant. Everything Unit 8 extracted from a filled sea then follows from operators acting on an empty vacuum, and more: positrons of positive energy and opposite charge, exclusion from (b†)² = 0, antisymmetric many-electron states and, through Unit 4's spin sums, an anticommutator that vanishes outside the light cone, so that bilinear observables commute there. The sea leaves a trace: its energy −Σ Eₚ and charge q Σ 1 are the infinite c-numbers that normal ordering drops. The anticommutator is an assumption justified by these consequences, not a deduction, and the model stops at free fields. Once the field couples to light the vacuum is no longer |0⟩, the modes Dirac pictured as filled respond to applied fields as vacuum polarisation, and handling that needs the renormalised field theory beyond this course.

1Keep the right-hand side, change the bracket

Units are SI, with metric (+, −, −, −), xμ = (ct, x), px = Eₚt − p⋅x for the four-vector product, and a periodic box of volume V so that momenta are discrete. Topic 9.05 fixed the dynamics and left the algebra open: π = iħψ† is the momentum conjugate to ψ, so the bracket of ψ with π must be iħδab δ³(x − y). Take it to be an anticommutator, {A, B} = AB + BA, and divide out iħ: (ψₐ(t, x), ψ†b(t, y)) = δab δ³(x − y), with (ψₐ, ψb) = 0. This is an assumption, judged by what follows. To read it in modes, multiply the expansion by e(ipx/ħ)u†(p, s) and integrate over the box. Only bₚₛ survives, because u†(p, s)u(p, s′) = (2Eₚ/c)δₛₛ′ and u†(p, s)v(−p, s′) = 0; v† isolates d†ₚₛ the same way. The rule then gives (bₚₛ, b†ₚ′ₛ′) = (dₚₛ, d†ₚ′ₛ′) = δₚₚ′δₛₛ′, every other pair anticommuting, derived. Dimensions agree: ψ†ψ is a number density, so ψ is in m(−3/2), matching δ³ in m⁻³, and b and d are dimensionless.

2A symmetric bracket lets d annihilate the vacuum

The operator d† multiplies the negative-frequency wave v e(ipx/ħ), and nothing in the field says whether d or d† should annihilate the vacuum. For a commutator the choice is costly: [d, d†] = −1 (topic 9.06), and swapping the roles of d and d† makes the norm positive only by sending the energy the wrong way. An anticommutator is symmetric, (d, d†) = (d†, d), so relabelling d̃ = d† gives (d̃, d̃†) = 1 again: for fermions the particle-hole swap is a canonical transformation, and for bosons it never is. Choose bₚₛ|0⟩ = dₚₛ|0⟩ = 0. Then ⟨0|d d†|0⟩ = ⟨0|(1 − d†d)|0⟩ = +1, and every product of creators acting on |0⟩ has norm 1, or 0 when it fills a mode twice. One mode pair (p, s) has exactly four states, none, an electron, a positron or both, so a box holding N momenta, each with two spins, has a Fock space of 4(2N) dimensions whose norm is positive definite.

3Reorder d d† and every energy is positive

Topic 9.05 derived H = Σₚₛ Eₚ(b†b − d d†) before any algebra was chosen, the cross terms between b and d† having died by u†(p, s)v(−p, s′) = 0. Now reorder: d d† = 1 − d†d, so H = Σₚₛ Eₚ(b†b + d†d) − Σₚₛ Eₚ, derived. Because [d†d, d†] = d†(1 − d†d) = d† for fermions too, [H, d†ₚₛ] = +Eₚd†ₚₛ: a positron raises the energy by Eₚ ≥ mc², exactly as an electron does, and |0⟩ is the ground state. In topic 9.06, d d† = d†d − 1 left −Eₚd†d, and with it negative norms or no floor. Take one mode pair with p = 1.00 MeV/c, so Eₚ = 1.123 MeV. Before normal ordering its four states have energies Eₚ(nb + nd − 1): −1.123, 0, 0 and +1.123 MeV. After it they are 0, 1.123, 1.123 and 2.246 MeV, with the same spacings, because only a constant went. The momentum reorders the same way into P = Σₚₛ p(b†b + d†d), the constant Σ p vanishing by symmetry, so the positron carries +p although d† multiplies the wave of momentum −p.

4The constants normal ordering drops are Dirac's sea

The charge reorders too: Q = q Σₚₛ(b†b + d d†) = q Σₚₛ(b†b − d†d) + q Σₚₛ 1, derived, and since [Q, d†] = −q d†, each d-quantum carries −q = +e for the electron field, where q = −e. Now read the constants: −Σₚₛ Eₚ and q Σₚₛ 1 are the energy and charge of every negative-energy level filled once. Mode by mode, the unordered operators are hole theory's bookkeeping: the empty mode pair carries energy −Eₚ and charge q, one sea electron, and adding a positron brings both to zero, that electron removed. Normal ordering, which moves annihilators to the right with a minus sign for each fermion swap, :d d†: = −d†d, subtracts the sea. What Unit 8 lacked is now in hand: a vacuum that is one state of zero energy and charge, positrons as quanta in their own right, electrons and positrons treated alike, and no infinitely charged sea that must be declared unobservable. The measured check: b and d quanta share Eₚ, so their masses agree, and |m(e⁺) − m(e⁻)|/m is measured below 8 × 10⁻⁹; annihilation at rest is seen to yield two 511 keV photons.

5Exclusion and antisymmetry are consequences

Put both entries equal in (b†ₚₛ, b†ₚ′ₛ′) = 0 and (b†ₚₛ)² = 0: a second electron in a filled mode gives the zero vector, not a state, derived. The occupation n = b†b then obeys n² = b†(1 − b†b)b = n, so its only eigenvalues are 0 and 1. For two modes, moving ψ(x)ψ(y) through the creators gives ⟨0|ψ(x)ψ(y)b†₁b†₂|0⟩ = φ₂(x)φ₁(y) − φ₁(x)φ₂(y), a Slater determinant of the waves φₖ = √(c/2EₖV) uₖ e(−ipₖx/ħ). Quantum Mechanics II postulated that antisymmetry; here the anticommutator is the only assumption. Two limits check it: with one quantum, ⟨0|ψ(x)b†ₖ|0⟩ = φₖ(x), Unit 4's plane wave normalised to one electron in V, and for p ≪ mc the upper components of the b part become the non-relativistic electron field Unit 10 starts from. A thermal consequence follows at once. With K = H − μN and N = Σ(b†b − d†d), each mode is a two-state system, so ⟨b†b⟩ = 1/(eˣ + 1) with x = (Eₚ − μ)/kB T, and ⟨d†d⟩ takes Eₚ + μ instead: Fermi-Dirac occupations capped at 1, plotted in the figure.

6The fields anticommute, the observables commute

At unequal times, insert the expansion and Unit 4's spin sums Σₛ uū = p̸ + mc and Σₛ vv̄ = p̸ − mc: (ψ(x), ψ̄(y)) = Σₚ (c/2EₚV)[(p̸ + mc)e(−ipz/ħ) + (p̸ − mc)e(ipz/ħ)], with z = x − y. Since iħγμ∂μ acting on e(∓ipz/ħ) gives ±p̸, this is (iħγμ∂μ + mc)[D(z) − D(−z)], with D the mode sum of topic 9.04's scalar amplitude up to a constant factor, derived. Topic 9.04 showed that D(z) = D(−z) for spacelike z in the continuum limit, so the bracket and its derivatives vanish there. At equal times the bracket returns γ⁰δ³(x − y), the canonical rule again. Without the d† term it would be (iħγμ∂μ + mc)D(z) alone, which falls off over a reduced Compton wavelength, ħ/mc = 386 fm, but is not zero. The fields anticommute rather than commute at spacelike separation, which is harmless, because they change fermion number and are not observables. Observables are bilinears, and [WX, YZ] = W{X, Y}Z − WY{X, Z} + {W, Y}ZX − Y{W, Z}X writes their commutator through anticommutators, all zero there, so observables commute at spacelike separation, derived.

2

Key equations

Read the units, assumptions and limits before substituting.

Canonical anticommutator (assumed)

(ψₐ(t, x), ψ†b(t, y)) = δab δ³(x − y)

The right-hand side is the one topic 9.06's commutator had; only the bracket's symmetry changes, and that is the assumption every later result tests.

ψ in m(−3/2), δ³ in m⁻³; a, b = 1 to 4 are spinor indices; it is the bracket with π = iħψ† divided by iħ

Mode expansion in a box of volume V

ψ(x) = Σₚₛ √(c/2EₚV) [bₚₛ u(p, s) e(−ipx/ħ) + d†ₚₛ v(p, s) e(ipx/ħ)]

Projecting with u† and v†, using u†(p, s)v(−p, s′) = 0, turns the field rule into (bₚₛ, b†ₚ′ₛ′) = (dₚₛ, d†ₚ′ₛ′) = δₚₚ′δₛₛ′, all other pairs anticommuting. Derived.

px = Eₚt − p⋅x; Eₚ = √(p²c² + m²c⁴) in J; V in m³; u†u = v†v = 2Eₚ/c; bₚₛ and dₚₛ are dimensionless operators

Hamiltonian: every quantum costs energy

H = Σₚₛ Eₚ (b†ₚₛbₚₛ + d†ₚₛdₚₛ) − Σₚₛ Eₚ

Derived from topic 9.05's Σ Eₚ(b†b − d d†) with d d† = 1 − d†d. Normal ordering, :d d†: = −d†d, removes only the constant; the sign flip is the anticommutator's work.

Eₚ in J; b†b and d†d count quanta and take the values 0 and 1; the c-number −Σₚₛ Eₚ is the filled sea's energy

Charge: d-quanta carry −q

Q = q Σₚₛ (b†ₚₛbₚₛ − d†ₚₛdₚₛ) + q Σₚₛ 1

Derived from q Σ(b†b + d d†). Since [Q, d†] = −q d†, each d-quantum carries +e: the positron, with the charge hole theory gave it, and no sea behind it.

q = −e = −1.602 × 10⁻¹⁹ C for the electron field; Σₚₛ 1 counts modes, so the constant is the sea's charge

Exclusion from the algebra

(b†ₚₛ)² = ½(b†ₚₛ, b†ₚₛ) = 0, so n² = n for n = b†ₚₛbₚₛ

No mode holds two electrons, and b†₁b†₂|0⟩ = −b†₂b†₁|0⟩: the antisymmetry Quantum Mechanics II postulated is derived here from the field's algebra.

n is a dimensionless occupation number with eigenvalues 0 and 1 only; the same holds for d†ₚₛdₚₛ

Anticommutator at any separation

(ψ(x), ψ̄(y)) = (iħγμ∂μ + mc)[D(x − y) − D(y − x)]

Derived; with topic 9.04's D(z) = D(−z) at spacelike z it vanishes there, though ⟨0|ψ(x)ψ̄(y)|0⟩ does not, so bilinear observables commute there.

D(z) = Σₚ (c/2EₚV) e(−ipz/ħ), in m⁻³ (kg m s⁻¹)⁻¹; uses Unit 4's spin sums Σₛ uū = p̸ + mc and Σₛ vv̄ = p̸ − mc
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

0.30 MeV
0.05 to 1.5 MeV; steps of 0.05.
0.20 MeV
−1.5 to 1.5 MeV; steps of 0.05.

Try this. Set μ = 0 and drag kB T down to 0.05 MeV: both curves sink to zero, an empty vacuum with no sea. Then raise μ to 1.00 MeV: electrons fill the modes below μ up to 1 and no further, the Pauli step, while positrons, whose x uses Eₚ + μ, disappear.

Live results

e⁻ PER MODE AT REST ⟨b†b⟩
0.262
e⁺ PER MODE AT REST ⟨d†d⟩
0.085
HOLE THEORY: SEA AT −mc² FULL
0.915
e⁺ TO e⁻ RATIO AT REST
0.327

Results follow this lesson’s model and assumptions. Displayed values are rounded.

Interactive Advanced Quantum Mechanics modelOccupation of one mode of the anticommuting Dirac field against its quantum's energy Eₚ, at k_B T = 0.30 MeV and electron chemical potential μ = 0.20 MeV. Solid: electrons, 1/(eˣ + 1) with x = (Eₚ − μ)/k_B T; dashed: positrons, with Eₚ + μ. No state lies below mc² (shaded); no curve passes 1. At rest: 0.262 electrons and 0.085 positrons per mode.solid e⁻: x = (Eₚ − μ)/kB T, at rest 0.262dashed e⁺: x = (Eₚ + μ)/kB T, at rest 0.085⟨n⟩ = 1/(eˣ + 1) ≤ 1, as (b†)² = 010mc²123energy of the quantum Eₚ / MeV
Occupation of one mode of the anticommuting Dirac field against its quantum's energy Eₚ, at kB T = 0.30 MeV and electron chemical potential μ = 0.20 MeV. Solid: electrons, 1/(eˣ + 1) with x = (Eₚ − μ)/kB T; dashed: positrons, with Eₚ + μ. No state lies below mc² (shaded); no curve passes 1. At rest: 0.262 electrons and 0.085 positrons per mode.

Predict → change one variable → check

Plan a change before running it. Every other input stays at the saved baseline.

Saved baseline. Temperature kB T: 0.30 MeV; Electron chemical potential μ: 0.20 MeV

Baseline 0.30; choose 0.05 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

The usual slip is to credit normal ordering with the positive energy. Normal ordering subtracts a c-number; it cannot change the sign of an operator. Apply it to topic 9.06's commutator theory and the operator part survives: with d annihilating the vacuum, :H: = Σₚₛ Eₚ(b†b − d†d) and d†|0⟩ has norm −1; with d† annihilating it, :H: = Σₚₛ Eₚ(b†b − d d†), each quantum d creates lowers the energy by Eₚ, and there is no ground state. The sign flip happens one step earlier, when d d† is reordered: an anticommutator gives d d† = 1 − d†d, a commutator gives d d† = d†d − 1. Only after that step is there a constant to drop, and the constant is not a nuisance term but the energy and charge of Dirac's filled sea. The two operations do different jobs. The anticommutator makes each positron cost +Eₚ; normal ordering discards the sea's infinite and unobservable energy −Σ Eₚ and charge q Σ 1, leaving every energy difference as it was.

5

Practice & worked examples

3 authored examples. Make an attempt, then reveal the reasoning one step at a time.

EasyExample 1

An electron of momentum 0.800 MeV/c along +z with spin up and a positron of momentum 0.800 MeV/c along −z with spin down are created from the vacuum: |Ψ⟩ = b†ₚ↑d†₋ₚ↓|0⟩. Using only the anticommutators, find the norm of |Ψ⟩ and its energy, momentum and charge from the normal-ordered operators. What changes if the positron is replaced by a second spin-up electron of the same momentum?

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

Before normal ordering the Dirac Hamiltonian and charge carry the constants −Σₚₛ Eₚ and q Σₚₛ 1. Keep only modes with |p| ≤ mc, a cutoff at the electron's own mass, and find the number of (p, s) modes per cubic metre and the charge and energy densities that normal ordering removes. Compare the first with the electron density of copper, 2.46 × 10³⁰ m⁻³.

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

Two spin-up electrons occupy box modes with momenta p₁ = +0.0100 MeV/c and p₂ = −0.0100 MeV/c along z: |Ψ⟩ = b†₁b†₂|0⟩. Using only the anticommutators, find the equal-time amplitude Ψab(x, y) = ⟨0|ψₐ(x)ψb(y)|Ψ⟩. Then find the pair density at a separation of 5.00 pm along z, relative to that of two distinguishable particles, and where it first peaks.

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.

With anticommutators, ψ(x) and ψ̄(y) anticommute when x − y is spacelike, but they do not commute there. A student calls this a violation of causality. Which reply is correct?
7

Review & continue

Quantising the Dirac field with anticommutators means imposing (ψₐ(t, x), ψ†b(t, y)) = δab δ³(x − y) instead of a commutator; it gives (b, b†) = (d, d†) = 1, so electrons and positrons are both quanta of positive energy and no mode can hold two.

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