Skip to main content

Advanced Quantum Mechanics course sequence · Interactive lesson

Dirac's sea and the radiative catastrophe

Unit 4 set the negative-energy solutions aside; this topic shows why they cannot stay there. We let a hydrogen electron radiate into them, find it would not last a microsecond, then follow Dirac's repair, filling every one, and count what the repair costs: an infinite vacuum and the end of a one-particle theory.

UAQM-08.03 · Lesson 3 of 9 in this unit. Position is not a completion score.
Starting model: Dirac's sea and the radiative catastropheLeft: hydrogen's 1s level just below +mc², and the negative-energy continuum below −mc², empty from −mc² to −2.0 mc² and filled (shaded) beneath. The arrow is a one-photon fall to −1.5 mc², emitting 1.277 MeV. Right: the lowest-order rate per unit landing energy, in α⁶mc²/ħ per mc², dashed for an empty sea and solid where states are still empty.sea empty from −mc² to −2.0 mc², shaded belowarrow: 1s → −1.5 mc², photon 1.277 MeVdashed: if the sea were empty1s ≈ +mc²−mc²rate per mc² of |E′||E′|/mc² from 1 to 6
Start with a picture. Change the model in the lab below.
1

Build the model

Empty, the sea drains hydrogen in nanoseconds; full, it gives the vacuum infinite charge.

The key idea

Dirac's sea is hole theory's vacuum: the state in which every negative-energy solution of the free Dirac equation, two spin states for each momentum, holds one electron, so that the Pauli principle forbids any further electron from radiating into the negative-energy continuum.

Example. Up to a momentum cutoff Λ = mc the sea holds Λ³/(3π²ħ³) = 5.9 × 10³⁵ electrons per cubic metre, a charge density of −9.4 × 10¹⁶ C m⁻³ and about 7 × 10⁶ times copper's conduction-electron density, all of it declared unobservable.

The complete physical model and assumptions

Minimal coupling hands the Dirac electron its photon vertex, H′ = −qcα⋅A, and the matrices αᵢ anticommute with β. At rest β is the sign of the energy, so the vertex has no matrix elements at all between resting positive-energy states, and at momentum p a fraction 1 − vₓ²/c² of what αₓ does lands in the negative-energy continuum. For a slow bound electron that is nearly everything, and nothing in the one-particle theory forbids the jump. Hydrogen's 1s electron can emit a photon of about 1.022 MeV or more and land below −mc², the nucleus absorbing the momentum the photon cannot balance; that costs α(Zα)⁵ rather than α, and the rate into the top band of the continuum alone is 2.2 × 10⁷ s⁻¹. Deeper bands give more, the sum over a continuum with no floor diverges, and every atom would collapse. Dirac's 1930 remedy fills every negative-energy state. Exclusion then blocks each of those transitions, the 1s level becomes the lowest an electron can reach, and the vacuum becomes a uniform sea that no one-electron measurement can see. The price is paid at once: the sea's charge density, −eΛ³/(3π²ħ³) below a momentum cutoff Λ, and its energy density, growing as Λ⁴, are infinite and are removed by declaring the filled sea the zero of charge and energy. The Dirac equation now describes one electron above infinitely many. That holds while the sea is left undisturbed and fails once it is not.

1The photon vertex points into the sea

Units are SI, and where a four-vector appears the metric is diag(1, −1, −1, −1); α with an index or a dot is a Dirac matrix, α alone the fine-structure constant. Minimal coupling gives H = cα⋅(p − qA) + βmc² + qφ, so a radiation field in Coulomb gauge acts through H′ = −qcα⋅A. Now use βαᵢβ = −αᵢ. At rest H₀ = βmc², so β is the sign of the energy, and for two positive-energy rest spinors ⟨+|αᵢ|+⟩ = ⟨+|βαᵢβ|+⟩ = −⟨+|αᵢ|+⟩ = 0: every matrix element of the vertex from a resting positive-energy electron ends at negative energy. At momentum p, sandwich αₓ between the projectors Λ± = (1 ± H₀/E)/2. Topic 8.01's H₀αₓ + αₓH₀ = 2cpₓ gives Λ₊αₓΛ₊ = (cpₓ/E)Λ₊, and since αₓ² = 1, a fraction vₓ²/c² of the weight αₓ carries out stays at positive energy while 1 − vₓ²/c² lands in the sea. Schrödinger's vertex p⋅A/m has no route to negative energy; for a slow electron, Dirac's is almost all route.

2A fall with no floor

Unit 7 puts the ground level at E₁ₛ = mc²√(1 − (Zα)²), and the negative-energy continuum runs from −mc² downwards without limit. A one-photon jump to a state of energy −|E′| releases ħω = E₁ₛ + |E′|: at least 1.0220 MeV less 13.6 eV for hydrogen, and still 0.890 MeV for Z = 92. A free electron cannot make the jump with one photon, because |p − p′| ≤ |p| + |p′| is less than (E + |E′|)/c = ħω/c, so the photon's momentum can never be matched; with two photons it can, so even a free electron would be unstable. In the atom the nucleus takes the excess, at least about mc, through its Coulomb field. Nothing stops the descent at the first landing: near the nucleus the electron at −|E′| can radiate again and land lower, and since |E′| has no maximum, neither does the energy one electron can give up. An atom in this theory has no ground state.

3Counting the powers of Zα in the rate

Count factors before doing any integral. The photon vertex costs e in amplitude, α in the rate. The momentum, at least mc, must come from one Coulomb exchange, Ṽ(Q) = −4πZαħ³c/Q², which at Q of order mc costs Zα in amplitude. The exchange happens where the electron is, which brings in the probability density at the nucleus, |ψ₁ₛ(0)|² = (Zαmc/ħ)³/π. Every other scale is mc, so Γ = C α(Zα)⁵ mc²/ħ, a scaling this lesson derives: a pure number times mc²/ħ, in s⁻¹, and vanishing as Z⁵ when Z → 0, as the free-electron argument above requires. The coefficient needs the full golden-rule sum, in which the Coulomb exchange may come before the photon or after it. Gauge invariance needs both orderings, and either alone gives a wrong answer. Summed numerically for this lesson, the two orderings reproduce the spectrum on the formula card, matching the quoted closed form f(w) to better than one part in 10⁴. Into the top band, |E′| from mc² to 2mc², the result is Γ = 0.192 α⁶mc²/ħ = 2.2 × 10⁷ s⁻¹ for hydrogen: a lifetime of 45 ns.

4Hydrogen is not falling

Set the numbers beside measurement. Hydrogen's 2p level decays in a measured 1.60 ns, as the golden rule's 1.596 ns predicts, so the top band alone would empty the ground state only 28 times more slowly than an excited level empties. Deeper bands make it worse: f(w) grows like w, so the band from mc² down to wmc² carries a rate that rises without limit, eventually as w²/π in units of α⁶mc²/ħ. The empty states down to −5.5mc² already drain the 1s level as fast as the 2p level decays, and the sum over the whole continuum diverges. Against that, hydrogen formed 13.8 billion years ago, 4.4 × 10¹⁷ s, is still most of the ordinary matter in the universe, and no 1 MeV glow rises from a flask of hydrogen gas. Even the top band alone misses by 25 orders of magnitude. Making the transitions slower would not rescue the theory; they have to be forbidden outright.

5Dirac's remedy: occupy every negative-energy state

Dirac's 1930 proposal makes the vacuum the state in which every negative-energy solution, two spin states for each momentum, already holds one electron. The Pauli principle, used here from Quantum Mechanics II, then blocks every transition counted above, because every final state is full. The 1s level becomes the lowest level an extra electron can reach, and a free positive-energy electron is stable for the same reason. The sea is invisible to one-electron experiments because it is perfectly uniform. Its momenta pair p with −p, so it carries no momentum and, since a negative-energy state moves with velocity −c²p/|E|, no current; its spins pair, so it has no spin; and a uniform background pushes nothing in any particular direction. A semiconductor shows the same structure in miniature: a filled valence band carries no current, and only excitations across the gap are seen. Silicon's gap is 1.12 eV; Dirac's is 2mc² = 1.022 MeV.

6The bill: infinite charge, energy and particle number

Exclusion works only if the sea is really full, so count it. Each phase-space cell (2πħ)³ holds two electrons, so up to a momentum cutoff Λ the density is n = Λ³/(3π²ħ³): 5.9 × 10³⁵ m⁻³ for Λ = mc, a charge density of −9.4 × 10¹⁶ C m⁻³, and infinite as Λ → ∞. The energy density, −(1/π²ħ³)∫p²√(p²c² + m²c⁴) dp up to Λ, is −6.05 × 10²² J m⁻³ at Λ = mc and grows as −cΛ⁴/(4π²ħ³). Dirac removes both by definition: charge, energy and momentum are measured from the filled sea, Qobs = Q − Qvac, and Maxwell's equations are sourced by the difference. The one-electron equation now describes one electron moving above infinitely many, a many-body theory written as a one-body one. The sea stops being invisible once it is disturbed: its polarisation by a proton lowers hydrogen's 2s level by about 27 MHz of the measured 1058 MHz Lamb shift, quoted. Taking an electron out is topic 8.04; the full list of costs is topic 8.08.

2

Key equations

Read the units, assumptions and limits before substituting.

The photon vertex is odd under β

βαᵢβ = −αᵢ, so ⟨+, 0, s′|αᵢ|+, 0, s⟩ = 0

Derived. At rest β is the energy sign, so every matrix element the vertex has from a resting positive-energy electron ends at negative energy.

H′ = −qcα⋅A, q = −e = −1.602 × 10⁻¹⁹ C; |±, 0, s⟩ are rest spinors with β = ±1; αᵢ and β are dimensionless 4 × 4 matrices

Where αₓ sends a moving electron

Σ over s′ of |⟨−, p, s′|αₓ|+, p, s⟩|² = 1 − c²pₓ²/E²

Derived from Λ₊αₓΛ₊ = (cpₓ/E)Λ₊, itself from topic 8.01's H₀αₓ + αₓH₀ = 2cpₓ. At v = 0 it gives 1, the rest result: a slow electron's vertex points into the sea.

E = √(p²c² + m²c⁴) in J, pₓ in kg m s⁻¹, so cpₓ/E = vₓ/c; the remaining c²pₓ²/E² stays at positive energy

Energy the first photon carries off

ħω = E₁ₛ + |E′| ≥ mc²[1 + √(1 − (Zα)²)]

Derived from energy conservation. There is no upper limit, because |E′| has none: one electron can radiate an unbounded energy.

E₁ₛ from Unit 7, −|E′| ≤ −mc² the landing energy, ħω in J; hydrogen 1.0220 MeV less 13.6 eV, Z = 92 still 0.890 MeV

Rate per unit landing energy, empty sea

dΓ/d|E′| = (2/π) α(Zα)⁵ (mc²/ħ) f(w)/mc², w = |E′|/mc²

α(Zα)⁵ derived by counting; f quoted, and checked here by summing the golden rule numerically with both time orderings. f grows like w, so the total over the sea diverges.

f(w) = [√(w² − 1)(3w² + 2w + 4) − 3(w + 2) ln(w + √(w² − 1))]/(3(w + 1)²); Γ in s⁻¹; mc²/ħ = 7.763 × 10²⁰ s⁻¹

Sea density below a momentum cutoff Λ

n = 2 × (4πΛ³/3)/(2πħ)³ = Λ³/(3π²ħ³)

Derived. Λ = mc gives 5.9 × 10³⁵ m⁻³ and −9.4 × 10¹⁶ C m⁻³; as Λ → ∞ both diverge, so the sea's charge has to be subtracted by definition.

Two spin states per cell (2πħ)³; Λ in kg m s⁻¹ over ħ in J s leaves n in m⁻³; the sea's charge density is ρ = −en in C m⁻³

Sea energy density below Λ

ε = −(1/π²ħ³) ∫ p²√(p²c² + m²c⁴) dp, p from 0 to Λ

Derived. Dirac's subtraction sets Evac = 0 and measures every energy from the filled sea: the first infinity the theory removes by hand.

ε in J m⁻³; for Λ ≫ mc it tends to −cΛ⁴/(4π²ħ³); Λ = mc gives −6.05 × 10²² J m⁻³, 670 times water's rest-energy density
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

2.0 mc²
1 to 6 mc²; steps of 0.1.
1.5 mc²
1 to 6 mc²; steps of 0.1.

Try this. Slide the first control to 1 mc², closing the empty band: the solid curve vanishes and the rate at every landing energy reads zero, which is Dirac's filled sea. Then open it to 6 mc² and drag the landing point down: the rate per unit energy keeps rising, so an empty sea gives a divergent total.

Live results

PHOTON ENERGY ħω
1.277 MeV
RATE PER mc² OF |E′|
2.10 10⁷ s⁻¹ / mc²
EMPTY SEA STATES
30.5 10³⁵ m⁻³
MOMENTUM TO NUCLEUS ≥
1.38 mc

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

Interactive Advanced Quantum Mechanics modelLeft: hydrogen's 1s level just below +mc², and the negative-energy continuum below −mc², empty from −mc² to −2.0 mc² and filled (shaded) beneath. The arrow is a one-photon fall to −1.5 mc², emitting 1.277 MeV. Right: the lowest-order rate per unit landing energy, in α⁶mc²/ħ per mc², dashed for an empty sea and solid where states are still empty.sea empty from −mc² to −2.0 mc², shaded belowarrow: 1s → −1.5 mc², photon 1.277 MeVdashed: if the sea were empty1s ≈ +mc²−mc²rate per mc² of |E′||E′|/mc² from 1 to 6
Left: hydrogen's 1s level just below +mc², and the negative-energy continuum below −mc², empty from −mc² to −2.0 mc² and filled (shaded) beneath. The arrow is a one-photon fall to −1.5 mc², emitting 1.277 MeV. Right: the lowest-order rate per unit landing energy, in α⁶mc²/ħ per mc², dashed for an empty sea and solid where states are still empty.

Predict → change one variable → check

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

Saved baseline. Sea empty down to |E| (filled below): 2.0 mc²; Landing energy |E′| of the electron: 1.5 mc²

Baseline 2.0; choose 1 to 6, in steps of 0.1.
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 tempting repair is to throw the negative-energy solutions away and keep the Dirac equation as a one-particle theory on its positive-energy half. That cannot be done consistently. The positive-energy states are not closed under interaction: α⋅A is odd, so acting on a positive-energy state it puts weight into the sea, nearly all of it for a slow electron, and projecting that weight away changes the predictions rather than tidying them. The change is measurable. At second order a Dirac electron scatters low-energy light only through negative-energy intermediate states; drop them and the predicted cross-section falls to zero instead of Thomson's (8π/3)rₑ² = 6.65 × 10⁻²⁹ m², the low-energy limit that measured Compton cross-sections approach (quoted). Dirac's sea keeps every negative-energy state in the theory; it only decides who may occupy them.

5

Practice & worked examples

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

EasyExample 1

An electron moves at 0.80c along x. Of the unit weight the vertex matrix αₓ carries out of its positive-energy spinor, what fraction lands in negative-energy states of the same momentum, and what fraction does αᵧ send there? Then estimate the αₓ fraction for hydrogen's 1s electron.

Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.

0 of 4 steps revealed.

MediumExample 2

In the empty-sea theory hydrogen's 1s electron radiates into negative-energy states at a rate per unit landing energy of g(w) α⁶mc²/ħ per mc², with w = |E′|/mc² and g = 0, 0.1788 and 0.4442 at w = 1, 1.5 and 2. Estimate the rate into the band from −mc² to −2mc², its lifetime and photon energies, and compare with the measured 2p lifetime of 1.60 ns.

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

Give Dirac's sea a momentum cutoff Λ = mc. Find the number of sea electrons inside a sphere of radius a₀ = 5.292 × 10⁻¹¹ m, their charge, and their energy, using ∫ x²√(1 + x²) dx = 0.4202 for x from 0 to 1. By what factors do the number and the energy grow if Λ = 10mc, where the same integral up to 10 is 2525? Why does hydrogen show none of it?

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.

In the empty-sea Dirac theory, the rate at which a hydrogen-like 1s electron radiates one photon into the top of the negative-energy continuum scales as α(Zα)⁵mc²/ħ, although the vertex αₓ sends a slow electron into the sea with weight close to 1. Where do the five powers of Zα come from?
7

Review & continue

Dirac's sea is hole theory's vacuum: the state in which every negative-energy solution of the free Dirac equation, two spin states for each momentum, holds one electron, so that the Pauli principle forbids any further electron from radiating into the negative-energy continuum.

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