Skip to main content

Advanced Quantum Mechanics course sequence · Interactive lesson

Dirac field with commutators: energy unbounded below

Here a quantisation rule is followed to the end and fails, and it fails where students do not expect. Commutators give back the Dirac equation and stay causal; what they cannot give is a vacuum that is lowest in energy with only positive norms above it. Choose which operator kills it, and either some states have negative norm or the energy has no floor.

UAQM-09.06 · Lesson 6 of 9 in this unit. Position is not a completion score.
Starting model: Dirac field with commutators: energy unbounded belowEnergy above the vacuum against the number n of quanta in one d-mode of momentum p = 0.75 mc, under the commutator rule [d, d†] = −1. Dashed, reading 1 (d|0⟩ = 0): (d†)ⁿ|0⟩ sits at +nEₚ but has norm (−1)ⁿn!, open marker when negative. Solid, reading 2 (d†|0⟩ = 0): dⁿ|0⟩ has norm n! but sits at −nEₚ = -1.916 MeV, and the line never stops.reading 1: d|0⟩ = 0, dashedE = 3 Eₚ, Q = −3 q, norm −6reading 2: d†|0⟩ = 0, solidE = −3 Eₚ, Q = 3 q, norm 6+10 mc²−10 mc²vacuumn quanta →
Start with a picture. Change the model in the lab below.
1

Build the model

The commutator forces [d, d†] = −1: negative norms, or an energy with no floor.

The key idea

Quantising the Dirac field with commutators imposes [ψₐ(x), ψ†b(y)] = δab δ³(x − y); the spinors' completeness then forces [d, d†] = −1, so either d†|0⟩ has norm −1 or the energy has no lower bound.

Example. In the reading with positive norms, three d-quanta in a mode with |p| = mc sit 3√2 mc² = 2.17 MeV below the vacuum with charge 3q = −4.81 × 10⁻¹⁹ C; a fourth takes the state to 2.89 MeV below, and nothing forbids a fifth.

The complete physical model and assumptions

Topic 9.05 left the Dirac field with H = Σ Eₚ(b†b − d d†) and Q = qΣ(b†b + d d†), and no rule yet for the operators. The rule that worked for every boson so far is the canonical commutator, and since the momentum conjugate to ψ is iħψ†, it reads [ψₐ(x), ψ†b(y)] = δab δ³(x − y) at equal times. This topic carries that rule to the end. The completeness of Unit 4's spinors turns it into one matrix condition per mode, and because the v-term of ψ carries d†, the condition fixes [b, b†] = +1 and [d, d†] = −1, derived. The minus sign is not a labelling slip: a commutator is antisymmetric, and unlike the scalar's momentum φ̇†/c², iħψ† has no time derivative to cancel it. What remains is to say which of d and d† annihilates the vacuum. If d does, its quanta have positive energy and the positron's charge but norm −1, so probabilities fail. If d† does, every norm is positive, but each quantum d creates carries −Eₚ and the electron's charge, and with no exclusion one mode holds any number of them, so H has no lower bound. The Heisenberg equation is the Dirac equation in both readings, free states are stationary, and the field's c-number bracket is causal: the dynamics is sound, and what fails is positivity. The measured positron, of positive energy, charge +e and ordinary probabilities, fits neither reading.

1Impose the commutator and let completeness choose the signs

SI units with ħ and c kept, metric (+, −, −, −). Topic 9.05's expansion, quoted, is ψ(x, t) = Σₚ,ₛ √(c/2EₚV) [b(p, s)u(p, s)e(−i(Eₚt − p⋅x)/ħ) + d†(p, s)v(p, s)e(i(Eₚt − p⋅x)/ħ)] in a box of volume V. Let [b, b†] = A and [d, d†] = B for equal labels, every other pair commuting. The d-term reaches [ψ, ψ†] as [d†, d] = −B, so after p → −p in that term the equal-time bracket is Σₚ (c/2EₚV) e(ip⋅(x − y)/ħ) [A Σₛ u u†(p) − B Σₛ v v†(−p)]. Unit 4's spin sums for u ū and v v̄, multiplied on the right by γ⁰, give Σₛ u u†(p) = (Eₚ + HD(p))/c and Σₛ v v†(−p) = (Eₚ − HD(p))/c, with HD(p) = cα⋅p + βmc². Since δab δ³(x − y) = (1/V)Σₚ e(ip⋅(x − y)/ħ) δab, every p needs (A − B)Eₚ + (A + B)HD(p) = 2Eₚ I. HD(p) has eigenvalues ±Eₚ and is no multiple of I, so A + B = 0 and A − B = 2: A = 1 and B = −1, derived.

2Where the minus sign comes from, and why the scalar escaped it

The minus is an ordering sign. ψ multiplies v by d†, and ψ† multiplies v† by d, so [ψ, ψ†] meets the pair in the order [d†, d] = −[d, d†]; the canonical rule needs that term to enter with +1, so [d, d†] = −1. Topic 9.03's complex scalar had the same ordering, b† on its negative frequencies, and still gave [b, b†] = +1, because its rule pairs φ with π = (1/c²)∂φ†/∂t, and the time derivative gives the term carrying b the factor −iω where the term carrying a† gets +iω: a second minus sign, cancelling the first. The Dirac Lagrangian is first order in time, its momentum iħψ† carries no derivative, and nothing cancels. Nor can renaming help, since [d, d†] = −1 and [d†, d] = +1 are one statement. The rest frame shows the sign with the algebra stripped away: at p = 0, HD = βmc², the u's fill the upper two components and the v's the lower two, and the lower block reads −B × 2mc = 2mc. Dimensions agree: √(c/2EₚV) u carries m(−3/2), as ψ must, so both sides of the rule are in m⁻³.

3Reading 1: d|0⟩ = 0, and norms go negative

Take b|0⟩ = d|0⟩ = 0, as the names suggest. Then ⟨0|d d†|0⟩ = ⟨0|(d†d − 1)|0⟩ = −1, and since d(d†)ⁿ = (d†)ⁿd − n(d†)ⁿ⁻¹, induction gives ⟨0|dⁿ(d†)ⁿ|0⟩ = (−1)ⁿ n!, derived. The energies are fine. With H = Σ Eₚ(b†b − d d†), [H, d†] = −Eₚ[d, d†]d† = +Eₚd†, so d† raises the energy by Eₚ, and [Q, d†] = −qd† gives each quantum the charge −q. Positive energy, the electron's mass, the opposite charge: every one-particle property of a positron except a positive norm, and that one is fatal. The state 0.6|0⟩ + 0.8d†|0⟩ has norm 0.36 − 0.64 = −0.28, and (|0⟩ + d†|0⟩)/√2 has norm zero without being the zero vector. Born-rule ratios built on such norms come out negative or above 1, so no measurement statistics can be read from them. An electron, b†|0⟩, still has norm +1: the failure lives entirely in the antiparticle sector, and the one-particle electron limit sees nothing wrong.

4Reading 2: d†|0⟩ = 0, and the energy has no floor

Positivity returns if the roles swap: let d†|0⟩ = 0. Then c = d† obeys [c, c†] = +1, an ordinary oscillator, and ⟨0|(d†)ⁿdⁿ|0⟩ = n!, all positive. Now H = Σ Eₚ(b†b − c†c) and Q = qΣ(b†b + c†c), normal ordered already, with no constants to drop. Each quantum d creates carries energy −Eₚ and charge q, the electron's own: Dirac's negative-energy electrons, now bosons. One mode holds any number n of them at energy −nEₚ, and Eₚ grows without limit with |p|, so H has no lowest eigenvalue, derived. Hole theory's cure is gone: without exclusion no state fills a mode, so Unit 8's sea cannot be built. Nothing of charge −q exists in this reading, against the measured positron, of charge +e and rest energy +0.511 MeV, whose annihilation at rest gives the observed pair of 511 keV photons. And [b, b†] = +1 in both readings, so b†(p, s)²|0⟩ ≠ 0: electrons would be bosons, a further failure that topic 9.07 removes.

5The dynamics is innocent in both readings

The Heisenberg equation iħ∂ψₐ/∂t = [ψₐ, H] needs only [ψₐ(x), ψ†b(y)Xb(y)] with X = HDψ. With commutators it is [ψₐ, ψ†b]Xb + ψ†b[ψₐ, Xb] = δab δ³(x − y)Xb, because ψ commutes with ψ; with anticommutators the identity [A, BC] = {A, B}C − B{A, C} gives the same. Either way iħ∂ψ/∂t = −iħc α⋅∇ψ + βmc²ψ, derived. Mode by mode, [b, H] = Eₚb and [d, H] = −Eₚ[d, d d†] = +Eₚd, so b(t) = b e(−iEₚt/ħ) and d†(t) = d† e(iEₚt/ħ): the frequencies the expansion assumed. No reading enters, because the Heisenberg equation is an operator identity, blind to which state is called the vacuum. Extended to unequal times, the bracket [ψₐ(x), ψ†b(y)] is the same c-number function that topic 9.07 shows vanishes at spacelike separation, quoted, so causality is intact as well. Neither the field equation nor microcausality can reject this rule; only the norm or the spectrum can.

6Nothing falls yet: what fails is the ground state

It is tempting to picture reading 2's vacuum collapsing into d-quanta. In a free theory it cannot: H commutes with every b†b and d d†, so every Fock state is stationary and no occupation ever changes. What fails is existence, not stability. H has no lowest eigenvalue, so no state can serve as the vacuum, and no thermal state exists: one mode's partition sum Σₙ e(nEₚ/kT) grows term by term by e(Eₚ/kT) = e(1.98 × 10⁷) at 300 K for p = 0, and diverges. Switch on any coupling to photons, as Unit 10 will, and conservation laws protect nothing: two d-quanta at rest, energy −2mc² and charge 2q, may become two at ±p plus two back-to-back photons of Eₚ − mc² each, allowed at every p, 1.10 MeV per photon at p = 3mc, so the energy released has no bound. Ordinary matter has a ground state and thermalises. For spin ½, commutators fail on positivity, the norm or the spectrum; topic 9.08 shows the scalar with anticommutators failing on causality instead.

2

Key equations

Read the units, assumptions and limits before substituting.

The rule under test: an equal-time commutator

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

One [x, p] = iħ per point and component, the rule that worked for the scalar field; this topic follows it through the modes to see whether it survives.

ψ in m(−3/2), δ³ in m⁻³; a, b = 1 to 4 label spinor components. Imposed, from π = iħψ† and [ψₐ, πb] = iħ δab δ³.

Completeness of u and v

Σₛ [u(p, s)u†(p, s) + v(−p, s)v†(−p, s)] = (2Eₚ/c) I

The u-part alone is (Eₚ + HD(p))/c and the v-part (Eₚ − HD(p))/c, so a mode's two brackets enter with I and with HD(p), and must match separately.

u†u = v†v = 2Eₚ/c in kg m s⁻¹; Eₚ = √(p²c² + m²c⁴) in J; I is the 4 × 4 unit matrix. Derived from Unit 4's spin sums.

What the commutator forces on the d-modes

[d(p, s), d†(p′, s′)] = −δₚₚ′ δₛₛ′

ψ carries d† on its v-term, so the field bracket meets [d†, d] = −[d, d†]. The scalar's momentum φ̇†/c² flips that sign back; iħψ† has no time derivative to do so.

b and d dimensionless; p, p′ box momenta, s, s′ spins. Derived, with [b, b†] = +δₚₚ′ δₛₛ′ from the same matching.

Reading 1, d|0⟩ = 0: norms alternate in sign

⟨0|dⁿ(d†)ⁿ|0⟩ = (−1)ⁿ n!

From [H, d†] = Eₚd† and [Q, d†] = −qd†, each quantum has energy +Eₚ and charge −q: a positron in every one-particle property but the norm, which is −1 for one quantum.

n = 0, 1, 2, … quanta in one mode, dimensionless. Derived by induction from d(d†)ⁿ = (d†)ⁿd − n(d†)ⁿ⁻¹.

Reading 2, d†|0⟩ = 0: an energy with no floor

E − Evac = Σₚ,ₛ Eₚ(nb − nd), with nb, nd = 0, 1, 2, …

Each d-quantum carries −Eₚ and the electron's charge q; one mode holds any number, and Eₚ grows with |p|, so H has no lowest eigenvalue to call the vacuum.

Eₚ in J; nb and nd are the eigenvalues of b†b and d d† in mode (p, s). Norms nb! nd!, all positive. Derived.

The Heisenberg equation does not choose

iħ ∂ψ/∂t = [ψ, H] = −iħc α⋅∇ψ + βmc²ψ

The Dirac equation returns in both readings and under topic 9.07's rule too, so no test of the dynamics rejects the commutator; only norms and the spectrum can.

Derived under either rule, since [ψₐ(x), ψ†b(y)Xb(y)] = δab δ³(x − y)Xb(y) for commutators and anticommutators alike.
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

3
0 to 6; steps of 1.
0.75 mc
0 to 2 mc; steps of 0.25.

Try this. Step n from 0 to 6: the dashed marker flips between filled and open as the norm changes sign, while the solid one only sinks, by {0.511*√(1+p²)|3} MeV per quantum. Raise p to steepen both lines; no setting gives positive norms and a floor together.

Live results

MODE ENERGY Eₚ
0.639 MeV
READING 1 NORM
-6
READING 2 ENERGY
-1.916 MeV
READING 2 CHARGE
-3 e

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

Interactive Advanced Quantum Mechanics modelEnergy above the vacuum against the number n of quanta in one d-mode of momentum p = 0.75 mc, under the commutator rule [d, d†] = −1. Dashed, reading 1 (d|0⟩ = 0): (d†)ⁿ|0⟩ sits at +nEₚ but has norm (−1)ⁿn!, open marker when negative. Solid, reading 2 (d†|0⟩ = 0): dⁿ|0⟩ has norm n! but sits at −nEₚ = -1.916 MeV, and the line never stops.reading 1: d|0⟩ = 0, dashedE = 3 Eₚ, Q = −3 q, norm −6reading 2: d†|0⟩ = 0, solidE = −3 Eₚ, Q = 3 q, norm 6+10 mc²−10 mc²vacuumn quanta →
Energy above the vacuum against the number n of quanta in one d-mode of momentum p = 0.75 mc, under the commutator rule [d, d†] = −1. Dashed, reading 1 (d|0⟩ = 0): (d†)ⁿ|0⟩ sits at +nEₚ but has norm (−1)ⁿn!, open marker when negative. Solid, reading 2 (d†|0⟩ = 0): dⁿ|0⟩ has norm n! but sits at −nEₚ = −1.916 MeV, and the line never stops.

Predict → change one variable → check

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

Saved baseline. Quanta n in the d-mode: 3; Mode momentum |p|: 0.75 mc

Baseline 3; choose 0 to 6, in steps of 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 relabel: the minus sign sits in [d, d†], so swap the names of d and d†, as topic 9.03 put b† on the scalar's negative frequencies. It cannot work. A commutator is antisymmetric, so [d, d†] = −1 and [d†, d] = +1 are the same statement, and renaming only decides which operator annihilates the vacuum, which is the choice between the two readings. Reading 1 keeps positive energies and loses positive norms; reading 2 keeps positive norms and loses the floor of the energy. Relabelling moves the failure from one to the other and never removes it. The scalar escaped because the time derivative in its momentum φ̇†/c² supplies a second minus sign; iħψ† has no derivative and supplies none. Only a symmetric bracket, with d d† + d†d = 1, lets d annihilate the vacuum with positive norm and positive energy at once, which is topic 9.07.

5

Practice & worked examples

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

EasyExample 1

Quantise the Dirac field with commutators and take the reading d†|0⟩ = 0. Find the energy relative to the vacuum, the charge and the norm of b†(0, ↑)d(p, ↑)³|0⟩: one electron at rest and three d-quanta in a mode with |p| = mc. Use mc² = 0.511 MeV and q = −e.

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

Take the reading d|0⟩ = 0 for one d-mode with |p| = 0.75 mc and the state |χ⟩ = 0.6|0⟩ + 0.8d†|0⟩. Find ⟨χ|χ⟩, each basis state's Born-rule share, and the energy ratio ⟨χ|H|χ⟩/⟨χ|χ⟩ measured from the vacuum, and compare it with the spectrum. Use mc² = 0.511 MeV.

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

Take p = 0.75 mc along z, so Eₚ = 1.25 mc², in the Dirac representation, and suppose [b, b†] = A and [d, d†] = B for this mode. Using Σₛ u u†(p) = (Eₚ + HD(p))/c and Σₛ v v†(−p) = (Eₚ − HD(p))/c, find the (1, 1) and (1, 3) elements of the mode's bracket matrix A Σₛ u u†(p) − B Σₛ v v†(−p), fix A and B from the canonical commutator, and say what B = +1 would give.

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

0 of 6 steps revealed.

6

Quick check

Select an answer and explain it before checking.

Quantise the Dirac field with commutators, so [d, d†] = −1, and take the reading in which d†|0⟩ = 0. For two quanta made in one mode, the state d(p, s)²|0⟩, what are its energy relative to the vacuum, its charge and its norm? Here q is the electron's charge.
7

Review & continue

Quantising the Dirac field with commutators imposes [ψₐ(x), ψ†b(y)] = δab δ³(x − y); the spinors' completeness then forces [d, d†] = −1, so either d†|0⟩ has norm −1 or the energy has no lower bound.

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