Skip to main content

Advanced Quantum Mechanics course sequence · Interactive lesson

Negative-energy components and Zitterbewegung

Unit 4 left two loose ends: negative-energy solutions with no meaning, and a velocity operator whose only eigenvalues are ±c. This topic ties them together. Solve the free Heisenberg equations exactly and the trembling they predict turns out to be interference between the two energy branches, absent from any electron built from one.

UAQM-08.01 · Lesson 1 of 9 in this unit. Position is not a completion score.
Starting model: Negative-energy components and ZitterbewegungTop: velocity ⟨cα_z⟩/c against time for an electron of momentum p = 0.75 mc along z with weight 0.20 on the negative-energy branch, phased to start at its peak; dashed lines mark ±c and the group velocity c²p/E. Bottom: position minus drift, in units of ħ/2mc. Drift 0.36 c; tremble amplitude 99 fm.+c−c⟨cαz⟩ / cdashed: ±c and c²p/E = 0.60 cx minus drift, units of ħ/2mc = 193 fm+10t from 0 to 4π ħ/mc² = 1.62 × 10⁻²⁰ s
Start with a picture. Change the model in the lab below.
1

Build the model

The velocity eigenvalues ±c belong to states that mix both signs of energy.

The key idea

Zitterbewegung is the oscillation of a free Dirac electron's position at angular frequency 2E/ħ, with amplitude of order ħ/2mc, produced by interference between the positive- and negative-energy components of its spinor.

Example. An electron state at rest split equally between the two branches oscillates at 2mc²/ħ = 1.55 × 10²¹ s⁻¹, a period of 4.05 × 10⁻²¹ s, through ±ħ/2mc = ±193 fm; a packet of positive energy only does not oscillate at all.

The complete physical model and assumptions

Unit 4's negative-energy solutions cannot be kept out of a localised state. A positive-energy spinor at rest, squeezed into a packet of width d, contains plane waves with momenta up to about ħ/d, and at each of them it is no longer an eigenvector of cα⋅p + βmc²: it overlaps the negative-energy branch with amplitude cp/(E + mc²). Below a width of ħ/mc = 386 fm the negative branch's share climbs from about a tenth towards a half. Once both are present, the free Heisenberg equations, solved exactly with nothing but (H, α) = 2cp, split the position operator into a uniform drift at c²p/H and a term oscillating at 2E/ħ, never slower than 1.55 × 10²¹ s⁻¹, with amplitude about ħ/2mc. That term's operator anticommutes with H, so it connects +E only to −E: it is a beat between the branches, and in any state of one sign its expectation is exactly zero. This is Zitterbewegung, and it settles Unit 4's puzzle, because ±c are eigenvalues of an instantaneous velocity whose eigenstates mix both signs, while a positive-energy packet moves at c²p/E, below c. The model stops where it starts to matter. Resolving a 193 fm tremble needs probes energetic enough to make pairs, where a one-electron wavefunction no longer describes the experiment; and read as a smearing of the electron over ħ/mc, it gives the Darwin term's form and size as a picture, while Unit 6 supplied the derivation.

1Squeezing a positive-energy spinor brings in the other branch

Units are SI and the metric is (+, −, −, −). Take Unit 4's spin-up rest spinor ξ = (1, 0, 0, 0)ᵀ and confine it: ψ(x, 0) = G(x)ξ, with G a Gaussian of rms width d. At momentum p the eigenvectors of H(p) = cα⋅p + βmc² with energy +E have their lower pair equal to cσ⋅p/(E + mc²) times the upper, and those with −E have the upper pair −cσ⋅p/(E + mc²) times the lower, so ξ is an eigenvector only at p = 0. Projecting it gives a negative-energy amplitude cp/(E + mc²) times the positive one, derived, and a negative fraction w₋(p) = (1 − mc²/E)/2. A wide packet holds only small p, and averaging w₋ ≈ p²/4m²c² over the Gaussian gives a total negative weight (3/16)(ħ/mcd)², derived: 1.9 × 10⁻³ at d = 10ħ/mc = 3.9 pm. Summed numerically, it is 0.11 at d = ħ/mc and 0.22 at d = ħ/2mc, tending to ½ as d → 0. Below the reduced Compton wavelength, 386 fm, one sign of energy is no longer enough.

2Solving the free Heisenberg equations exactly

Quoted from Quantum Mechanics II: dA/dt = (i/ħ)[H, A]. For H = cα⋅p + βmc², p commutes with H, and dx/dt = (i/ħ)[H, x] = cα, Unit 4's velocity operator. Its own derivative needs [H, α]. Because (αᵢ, αⱼ) = 2δᵢⱼ and (β, αᵢ) = 0, the anticommutator is simple: (H, α) = 2cp, derived. So αH = 2cp − Hα and dα/dt = (2i/ħ)(Hα − cp). Since H² = p²c² + m²c⁴ ≥ m²c⁴, H has an inverse that commutes with p, and η = α − cpH⁻¹ absorbs the constant: dη/dt = (2i/ħ)Hη, so η(t) = e(2iHt/ħ)η(0) and cα(t) = c²pH⁻¹ + c e(2iHt/ħ)η(0), derived. Two facts carry the rest. η anticommutes with H, since (H, η) = 2cp − 2cp = 0, so it turns an energy-E state into an energy −E state. And c²pH⁻¹ commutes with H, has commuting components, and on a plane wave of energy ±E takes the value ±c²p/E, smaller than c.

3Integrating once: a drift and a 2E/ħ tremble

Integrating cα(t) gives x(t) = x(0) + c²pH⁻¹t − (iħc/2)H⁻¹(e(2iHt/ħ) − 1)η(0), derived; with ηH = −Hη the last term is Hermitian, as a position must be. The second term is uniform motion. The third oscillates at 2E/ħ, with amplitude ħc/2E times the matrix elements of η, which are at most 1. Dimensions: ħc/E is J s × m s⁻¹ ÷ J = m, and 2E/ħ is s⁻¹. Limits: at rest the amplitude is ħ/2mc = 193 fm and the angular frequency 2mc²/ħ = 1.55 × 10²¹ s⁻¹; as p grows the beat quickens and ħc/2E shrinks, so fast electrons barely tremble. The combination X(t) = x(t) + (iħc/2)H⁻¹η(t) moves exactly uniformly, X(t) = X(0) + c²pH⁻¹t, derived: remove the tremble and a drift coordinate is left. No free electron's tremble has been observed, since 193 fm at 10²¹ s⁻¹ is beyond any probe that leaves the electron single; a 2010 trapped-ion simulation of the Dirac equation recorded it, quoted.

4Only a mixture of energy signs trembles

η anticommutes with H, so its matrix elements join +E only to −E. In a packet built from positive-energy spinors alone, ⟨ψ₊|H⁻¹(e(2iHt/ħ) − 1)η|ψ₊⟩ = 0 at every t, and ⟨x⟩ moves at ⟨c²p/E⟩ with no oscillation, derived. The tremble is interference. At one momentum and helicity the two branches form a two-level system, H = E n⋅τ with n = (cp/E, 0, mc²/E), α∥ acting as τₓ and β as τz: a positive-energy state is a Bloch vector along n and stands still, while any admixture precesses about n at 2E/ħ, which gives the two-branch formula. Atoms are almost pure. Projected on the free spinors, Unit 7's exact ground state has negative-energy weight 3.2 × 10⁻¹² in hydrogen and 2.9 × 10⁻³ in U⁹¹⁺, computed, mostly from momenta above 0.3mc, within about ħ/mc of the nucleus. Uranium's 13% small component belongs almost entirely to positive-energy spinors.

5Unit 4's velocity puzzle, resolved

Unit 4 found that cα has eigenvalues ±c along every axis and components that fail to commute, while no packet moved faster than light. The solution separates the two claims. The drift part c²pH⁻¹ is what a measured displacement over a time much longer than ħ/mc² = 1.29 × 10⁻²¹ s returns, and on positive energies it is the group velocity c²p/E: 0.707c at p = mc. The eigenvalues ±c belong to the instantaneous velocity, whose eigenstates mix both branches. In the two-level picture the α∥ = +1 eigenstate has w₋ = (1 − cp/E)/2: one half at rest, 1/26 at p = 2.4mc. To catch the electron moving at c you would have to resolve the 193 fm tremble within 10⁻²¹ s, which takes probes with momentum transfer above mc and energy transfer above 2mc² = 1.02 MeV, where pairs are created and the experiment is no longer about one electron. What remains, a velocity of ±c that only a pair-making probe could catch, is real in the one-particle theory and is dissolved, not solved, by the many-particle reading Units 8 and 9 build.

6The Darwin term as a smeared potential: a picture

δ = x − X = −(iħc/2)H⁻¹η equals ε[ε, x]/2 with ε = H/|H|: it is the odd part of x, the part joining the two branches, and X is the even part. Its mean vanishes in a one-sign state, but its square does not: δ⋅δ = (ħ²c²/4)H⁻²(3 − c²p²H⁻²), derived, so a slow electron of either sign has ⟨δ²⟩ = (3/4)(ħ/mc)² = 1.12 × 10⁻²⁵ m², an rms 334 fm, shared equally among the axes. Treat δ as a random displacement and a slowly varying potential is felt as ⟨V(X + δ)⟩ ≈ V + (1/6)⟨δ²⟩∇²V = V + (ħ²/8m²c²)∇²V, which is Unit 6's Darwin term with its coefficient. It stays a picture: δ's components do not commute with one another or with X, and only the Foldy-Wouthuysen reduction derives the term. The measurement agrees with the coefficient: in hydrogen it raises 2S₁/₂ by α⁴mc²/16 = 21.9 GHz, closing the gap to 2P₁/₂ to within the measured Lamb shift, 1.058 GHz.

2

Key equations

Read the units, assumptions and limits before substituting.

A rest spinor seen at momentum p

|b(p)| / |a(p)| = cp / (E + mc²)

The price of localisation: the negative-energy fraction at each p is w₋ = (1 − mc²/E)/2, so a packet narrower than ħ/mc is up to half negative energy.

a, b: amplitudes on the +E and −E branches, dimensionless; p in kg m s⁻¹; E = √(p²c² + m²c⁴) in J. Derived.

The anticommutator that drives the motion

(H, α) = 2cp

Turns dα/dt = (i/ħ)[H, α] into the linear equation dα/dt = (2i/ħ)(Hα − cp), which integrates in closed form.

H = cα⋅p + βmc² in J; α dimensionless; from (αᵢ, αⱼ) = 2δᵢⱼ and (β, αᵢ) = 0. Derived.

Velocity operator at time t

cα(t) = c²pH⁻¹ + c e(2iHt/ħ) η(0), with η = α − cpH⁻¹

The constant part is the drift velocity, with eigenvalues ±c²p/E inside ±c; the part that beats at 2E/ħ is where the eigenvalues ±c of cα come from.

η is dimensionless and anticommutes with H; H⁻¹ = H/(p²c² + m²c⁴) in J⁻¹. Derived from the free Heisenberg equation.

Position operator: drift plus tremble

x(t) = x(0) + c²pH⁻¹t − (iħc/2) H⁻¹ (e(2iHt/ħ) − 1) η(0)

Uniform motion at c²p/H plus an oscillation of amplitude about ħc/2E: Zitterbewegung, with the operator η deciding which states show it.

x in m, t in s; ħc/2E is a length, 193 fm at E = mc². Derived by integrating cα(t).

Frequency of the tremble

ωZB = 2E/ħ ≥ 2mc²/ħ = 1.55 × 10²¹ s⁻¹

A beat between +E and −E, not the phase of either: twice mc²/ħ at rest, a period of 4.05 × 10⁻²¹ s, and faster for a moving electron.

E is the magnitude of each branch's energy, in J; ħ in J s. Derived as the beat of e(−iEt/ħ) against e(+iEt/ħ).

Two-branch state at one momentum

⟨cα∥⟩(t) = (w₊ − w₋)c²p/E + 2√(w₊w₋)(mc²/E) c cos(2Et/ħ + φ)

Zero weight on either branch kills the beat. The peak never exceeds c and reaches it only at w₋ = (1 − cp/E)/2, an eigenstate of α∥.

w₊, w₋: branch weights summing to 1; φ their relative phase; α∥ = α⋅p̂, one helicity. Derived.
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

0.75 mc
0 to 3 mc; steps of 0.25.
0.20
0 to 1; steps of 0.05.

Try this. Set w₋ to 0 or 1 and the tremble vanishes: one sign of energy never jitters. At the default, p = 0.75mc and w₋ = 0.20, the state starts as the αz = +1 eigenvector and touches +c; raise p and the beat quickens as 2E/ħ while the velocity swing shrinks as mc²/E and the position tremble as (mc²/E)².

Live results

ANGULAR FREQUENCY 2E/ħ
1.94 10²¹ s⁻¹
TREMBLE AMPLITUDE
99 fm
DRIFT VELOCITY
0.36 c
PEAK VELOCITY
1.00 c

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

Interactive Advanced Quantum Mechanics modelTop: velocity ⟨cα_z⟩/c against time for an electron of momentum p = 0.75 mc along z with weight 0.20 on the negative-energy branch, phased to start at its peak; dashed lines mark ±c and the group velocity c²p/E. Bottom: position minus drift, in units of ħ/2mc. Drift 0.36 c; tremble amplitude 99 fm.+c−c⟨cαz⟩ / cdashed: ±c and c²p/E = 0.60 cx minus drift, units of ħ/2mc = 193 fm+10t from 0 to 4π ħ/mc² = 1.62 × 10⁻²⁰ s
Top: velocity ⟨cαz⟩/c against time for an electron of momentum p = 0.75 mc along z with weight 0.20 on the negative-energy branch, phased to start at its peak; dashed lines mark ±c and the group velocity c²p/E. Bottom: position minus drift, in units of ħ/2mc. Drift 0.36 c; tremble amplitude 99 fm.

Predict → change one variable → check

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

Saved baseline. Momentum p: 0.75 mc; Negative-energy weight w₋: 0.20

Baseline 0.75; choose 0 to 3, in steps of 0.25.
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

A common reading is that every Dirac electron trembles: that an electron in a beam or an atom jitters through ħ/2mc at 2mc²/ħ, and that this jitter is what the eigenvalues ±c of cα describe. It does not. The oscillating part of x(t) is built from η = α − cpH⁻¹, which anticommutes with H and so has matrix elements only between +E and −E. In any state of one energy sign its expectation vanishes at every time, and ⟨x⟩ moves uniformly at ⟨c²p/E⟩. Zitterbewegung is an interference term, present only when a state mixes the branches, as it must once squeezed below ħ/mc. What a positive-energy electron does keep is a spread: the odd part δ = x − X has ⟨δ⟩ = 0 but ⟨δ²⟩ = (3/4)(ħ/mc)² at rest, the size the Darwin picture uses. A spread is not a motion, and nothing in a one-sign state oscillates at 2E/ħ.

5

Practice & worked examples

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

EasyExample 1

A packet is built by multiplying Unit 4's spin-up rest spinor by a narrow Gaussian. At what momentum does a plane-wave component of this packet have a negative-energy amplitude half its positive-energy amplitude? What fraction of that component's probability is negative energy, and roughly how narrow must the packet be for such momenta to be typical?

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

An electron of momentum p = 2.4mc along z is prepared at t = 0 in the eigenstate of αz with eigenvalue +1 and definite helicity. Find its weights on the two energy branches, its drift velocity, and the angular frequency, velocity range and position amplitude of its Zitterbewegung. Compare with a packet of positive energy only.

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

Write the Dirac position as x = X + δ, with δ = −(iħc/2)H⁻¹η and η = α − cpH⁻¹. Show that δ⋅δ = (ħ²c²/4)H⁻²(3 − c²p²H⁻²), evaluate ⟨δ²⟩ for a slow electron, and smear the hydrogen Coulomb potential over δ. What shift does the picture give 2S₁/₂, and how does it compare with the measured 2S₁/₂ to 2P₁/₂ separation of 1.058 GHz?

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.

A broad electron packet has a sharply defined momentum p = (4/3)mc, with 10% of its probability on the negative-energy branch at that momentum and helicity and 90% on the positive. What are the angular frequency of its Zitterbewegung and its drift velocity?
7

Review & continue

Zitterbewegung is the oscillation of a free Dirac electron's position at angular frequency 2E/ħ, with amplitude of order ħ/2mc, produced by interference between the positive- and negative-energy components of its spinor.

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