
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.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.
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.Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
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.
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
This experiment is kept only while this page is open. It does not change saved course progress.
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/ħ.
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.
- At momentum p the rest spinor's amplitudes on the two branches have ratio |b|/|a| = cp/(E + mc²), derived in this topic. Write u = cp/mc², so E = mc²√(1 + u²), and the condition is u/(√(1 + u²) + 1) = ½.
- Rearranging, 2u − 1 = √(1 + u²); squaring, 4u² − 4u + 1 = 1 + u², so 3u² = 4u and u = 4/3. Then E = (5/3)mc² and p = (4/3)mc = 0.681 MeV/c = 3.64 × 10⁻²² kg m s⁻¹.
- Negative-energy fraction: w₋ = (1 − mc²/E)/2 = (1 − 3/5)/2 = 0.20. Check from the ratio: (½)²/(1 + (½)²) = 0.25/1.25 = 0.20, the same.
- A Gaussian of rms width d has momentum spread σₚ = ħ/2d per axis, so (4/3)mc is about one standard deviation when ħ/2d ≈ (4/3)mc, that is d ≈ (3/8)ħ/mc = 145 fm. Limits: p → 0 gives ratio 0, the rest spinor itself; p → ∞ gives ratio 1 and w₋ → ½.
Answer. p = (4/3)mc = 0.681 MeV/c, where E = (5/3)mc² and 20% of that component's probability is negative energy; such momenta are typical in a packet of rms width about 145 fm, below ħ/mc.
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.
- E = mc²√(1 + 2.4²) = 2.6mc², so cp/E = 12/13 and mc²/E = 5/13. At fixed p and helicity, H = E n⋅τ with n = (12/13, 0, 5/13) and αz acting as τₓ, so the αz = +1 eigenstate is the Bloch vector along x.
- Its weight on +E is (1 + nₓ)/2 = (1 + 12/13)/2 = 25/26, so w₋ = 1/26 = 0.0385: at this momentum a velocity eigenstate is 96% positive energy.
- Drift: (w₊ − w₋)c²p/E = (24/26)(12/13)c = (12/13)²c = 0.852c. Velocity amplitude: 2√(w₊w₋)(mc²/E)c = (5/13)²c = 0.148c, so ⟨cαz⟩ swings between 1.000c at t = 0 and 0.704c.
- Angular frequency: 2E/ħ = 5.2mc²/ħ = 5.2 × 7.763 × 10²⁰ s⁻¹ = 4.04 × 10²¹ s⁻¹. Position amplitude: velocity amplitude over ω, (25/169)c × ħ/(5.2mc²) = (125/4394)ħ/mc = 0.0284 × 386.2 fm = 11.0 fm.
- Checks: the peak is exactly c, the eigenvalue, as it must be at t = 0, and ħ/mc is a length. A positive-energy packet at this momentum, kinetic energy 1.6mc² = 0.818 MeV, moves steadily at c²p/E = 0.923c, the relativistic speed that time-of-flight measurements confirm: Bertozzi's 1964 runs at 0.5 and 1.0 MeV, which bracket this energy, found electron speeds in line with c²p/E, quoted.
Answer. w₋ = 1/26; drift 0.852c; ⟨cαz⟩ oscillates between 1.000c and 0.704c at 2E/ħ = 4.04 × 10²¹ s⁻¹, an 11.0 fm tremble. A positive-energy packet would move steadily at 0.923c.
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.
- δ⋅δ = −(ħ²c²/4)H⁻¹η⋅H⁻¹η. Moving H⁻¹ through η flips its sign, ηH⁻¹ = −H⁻¹η, so δ⋅δ = (ħ²c²/4)H⁻²(η⋅η).
- η⋅η = Σ(αₖ − cpₖH⁻¹)² = 3 − cpₖ(αₖ, H⁻¹) + c²p²H⁻², and (αₖ, H⁻¹) = H⁻¹(αₖ, H)H⁻¹ = 2cpₖH⁻², so η⋅η = 3 − c²p²H⁻² and δ⋅δ = (ħ²c²/4)H⁻²(3 − c²p²H⁻²). It depends on p and H² only, so it is the same for both energy signs.
- For p ≪ mc, H⁻² → 1/m²c⁴ and ⟨δ²⟩ = (3/4)(ħ/mc)² = 0.75 × (3.862 × 10⁻¹³ m)² = 1.12 × 10⁻²⁵ m², shared equally by the axes since αᵢαⱼ + αⱼαᵢ = 2δᵢⱼ. Smearing: ⟨V(X + δ)⟩ ≈ V + (1/6)⟨δ²⟩∇²V = V + (ħ²/8m²c²)∇²V.
- For V = −e²/4πε₀r, ∇²V = (e²/ε₀)δ³(r), and |ψ₂ₛ(0)|² = 1/8πa₀³, quoted. So ΔE = (ħ²/8m²c²)(e²/ε₀)/(8πa₀³); with e²/ε₀ = 4παħc and a₀ = ħ/αmc this is α⁴mc²/16 = 9.06 × 10⁻⁵ eV, or 21.9 GHz.
- Quantum Mechanics II's kinetic and spin-orbit terms, quoted, leave 2S₁/₂ exactly α⁴mc²/16 below 2P₁/₂, so the smear closes the gap. The measured separation, 1.058 GHz with 2S₁/₂ above, is the Lamb shift: the size is right to 5% of the gap closed.
- Limits and status: dimensions of (ħ/mc)²∇²V are those of V; ⟨δ²⟩ falls as p grows. The step from V(x) to a classical average over δ is not justified, since δ's components do not commute with one another or with X, so this reproduces the Darwin term without deriving it.
Answer. ⟨δ²⟩ = (3/4)(ħ/mc)² = 1.12 × 10⁻²⁵ m², rms 334 fm; the smear gives (ħ²/8m²c²)∇²V and raises 2S₁/₂ by 21.9 GHz, closing its gap to 2P₁/₂ to within the measured 1.058 GHz Lamb shift: the right size, as a picture.
Quick check
Select an answer and explain it before checking.
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.