
Advanced Quantum Mechanics course sequence · Interactive lesson
Canonical quantisation of the real Klein-Gordon field
In this lesson Unit 3's negative energies stop being a problem. We stop reading φ as a wavefunction, quantise it as a mechanical system with one canonical pair per point of space, and push the algebra until the Hamiltonian is Quantum Mechanics I's oscillator, once per mode. Every later field in the course is quantised this way.
UAQM-09.02 · Lesson 2 of 9 in this unit. Position is not a completion score.Build the model
One [x, p] = iħ per point of space turns a classical field into a set of oscillators.
The key idea
Canonical quantisation of a real field imposes [φ(x), π(y)] = iħ δ³(x − y) at equal times, which makes each normal mode of the field an independent quantum oscillator whose quanta are the field's particles.
Example. For the π⁰ field, mc² = 134.98 MeV, the mode with k = 1.00 fm⁻¹ is an oscillator whose quantum carries ħωₖ = √(197.33² + 134.98²) MeV = 239.1 MeV and momentum 197.3 MeV/c; two quanta in it carry 478.2 MeV.
The complete physical model and assumptions
A real Klein-Gordon field is a mechanical system with one coordinate, φ(x), at every point of space, and quantising it repeats Quantum Mechanics I's step for a particle once per point. Topic 9.01's density makes the momentum conjugate to φ equal to π = (1/c²)∂φ/∂t, and [φ(x, t), π(y, t)] = iħ δ³(x − y) is imposed, the limit of one [x, p] = iħ per lattice cell. The rest is derived, quoting only Quantum Mechanics I's oscillator and topic 9.01's Noether momentum. The field equation is linear, so φ expands in the plane waves of a periodic box; φ is Hermitian, so the operator on each negative-frequency wave must be the adjoint of the one on the matching positive-frequency wave. Inverting the expansion turns the field commutator into [aₖ, aₖ'†] = δₖₖ', and in the Hamiltonian the terms that would make or destroy pairs at ±k cancel by the dispersion relation, leaving H = Σₖ ħωₖ(aₖ†aₖ + ½): one oscillator per mode, its quanta carrying energy √(ħ²k²c² + m²c⁴) and momentum ħk. Unit 3's negative frequencies survive, but now multiply creation operators, which raise the energy. Creation operators commute, so the quanta are bosons by construction. The model stops in three places: the vacuum's Σ½ħωₖ is infinite, and normal ordering only sets it aside for Unit 11; a real field has no charge, so its quantum is its own antiparticle; and without interactions no quantum is ever made or destroyed.
1One canonical pair per point of space
SI units throughout, with xμ = (ct, x) and the metric (+, −, −, −) of Unit 3. Topic 9.01's density, ℒ = ½[(1/c²)(∂φ/∂t)² − (∇φ)² − (mc/ħ)²φ²] in J m⁻³, gives φ the unit √(J m⁻¹) and π = ∂ℒ/∂(∂φ/∂t) = (1/c²)∂φ/∂t. To see what to impose, cut space into cells of side ℓ. The Lagrangian becomes Σⱼ ℓ³ℒⱼ, so cell j has coordinate φⱼ and momentum pⱼ = ℓ³πⱼ, and Quantum Mechanics I's rule for each pair, [φᵢ, pⱼ] = iħδᵢⱼ, reads [φᵢ, πⱼ] = iħδᵢⱼ/ℓ³. As ℓ → 0, δᵢⱼ/ℓ³ becomes δ³(x − y). That is a postulate, and an equal-time one: it says nothing about φ and π at different times. Dimensions agree, φπ and iħδ³ both in J s m⁻³. A consistency check follows: with H = ∫ℋ d³x and ℋ = ½[c²π² + (∇φ)² + (mc/ħ)²φ²], the Heisenberg equations give ∂φ/∂t = c²π and ∂π/∂t = ∇²φ − (mc/ħ)²φ, so the operator φ obeys Unit 3's Klein-Gordon equation, derived.
2Expand in a box, and let Hermiticity place a†
The field equation is linear, so its operator solution is a sum of plane waves. In a periodic cube of side L the allowed wavevectors are k = 2πn/L, n a vector of integers, and e(i(k⋅x − ωt)) is a solution only for ω = ±ωₖ with ωₖ = c√(k² + (mc/ħ)²). Each k therefore brings a positive-frequency wave with some operator Aₖ and a negative-frequency wave e(−i(k⋅x − ωₖt)) with an operator Bₖ. In Unit 3 the two were independent, and the second was a negative-energy state. Here φ is Hermitian, and taking the adjoint of the sum swaps the two waves, so Bₖ = Aₖ† for every k: one operator per mode, not two. Writing Aₖ = √(ħc²/(2ωₖV)) aₖ, a factor in √(J m⁻¹) chosen so that aₖ is dimensionless with unit commutator, gives the expansion on the card. For the π⁰ field in a cube of side 20.0 fm the spacing is 2π/L = 0.314 fm⁻¹: the k = 0 mode has ħωₖ = 134.98 MeV, the six with |n| = 1 have 148.5 MeV, and 147 modes have |k| ≤ 1 fm⁻¹, against the continuum estimate VK³/(6π²) = 135.
3Invert the expansion to find the ladder algebra
At t = 0, project the field onto one wave: ∫d³x e(−ik⋅x)φ = √(ħc²V/(2ωₖ))(aₖ + a₋ₖ†), and ∫d³x e(−ik⋅x)π = −i√(ħωₖV/(2c²))(aₖ − a₋ₖ†). Adding the two with the right weights isolates aₖ = (2ħV)(−1/2)∫d³x e(−ik⋅x)[(√ωₖ/c)φ + i(c/√ωₖ)π]. Now feed in the postulate. In [aₖ, aₖ'†] only the φπ and πφ cross terms survive, and each contributes ½ through ∫d³x e(−i(k − k')⋅x) = Vδₖₖ', so [aₖ, aₖ'†] = δₖₖ', while [φ, φ] = [π, π] = 0 give [aₖ, aₖ'] = 0: derived. Once k' = k, the √ωₖ in each term's φ weight cancels the 1/√ωₖ in its π weight, which is why the mass never enters the algebra. Dimensions: both sides are pure numbers. The one-mode limit is Quantum Mechanics I's [a, a†] = 1, and the Heisenberg equation with the Hamiltonian of the next step gives aₖ(t) = aₖe(−iωₖt): the time dependence written into the expansion is what the operators do of their own accord.
4The Hamiltonian falls apart into oscillators
Substitute the expansion into H = ∫ℋ d³x. The c²π² term gives (ħωₖ/4)(aₖaₖ† + aₖ†aₖ − aₖa₋ₖe(−2iωₖt) − h.c.) for each k, and the gradient and mass terms together give (ħc²(k² + m²c²/ħ²)/(4ωₖ))(aₖaₖ† + aₖ†aₖ + aₖa₋ₖe(−2iωₖt) + h.c.). The second prefactor equals ħωₖ/4 only because ωₖ² = c²(k² + m²c²/ħ²), so the terms that would make or destroy pairs at ±k cancel, and H = Σₖ ħωₖ(aₖ†aₖ + ½), derived. Quantum Mechanics I's norm argument, quoted, makes each nₖ = 0, 1, 2, …, and topic 9.01's Noether momentum, P = −∫π∇φ d³x, becomes Σₖ ħk aₖ†aₖ. So aₖ†|0⟩ has energy ħωₖ = √(ħ²k²c² + m²c⁴) and momentum ħk: a free relativistic particle. Limits: m = 0 gives ħωₖ = ħc|k|; for k ≪ mc/ħ, ħωₖ ≈ mc² + ½ħ²k²/m, and at k = 0.1 fm⁻¹ for the π⁰ the two give 136.412 and 136.419 MeV. Measured: the π⁰, chargeless and its own antiparticle as a real field's quantum must be, has mc² = 134.977 MeV.
5Where Unit 3's negative energies went
The algebra gives [H, aₖ†] = ħωₖaₖ† and [P, aₖ†] = ħk aₖ†: acting on any state, aₖ† adds energy +ħωₖ and momentum +ħk, although in φ it multiplies e(−i(k⋅x − ωₖt)), which Unit 3 would have read as energy −ħωₖ and momentum −ħk. The sign of a frequency in the field is not the energy of any state. It records whether the operator beside it removes a quantum or adds one, and aₖ†(t) = aₖ†e(+iωₖt) is how every raising operator evolves in the Heisenberg picture. Every state has energy at least E₀, so there is nothing below to cascade into. The one-quantum wavefunction, ⟨0|φ(x, t)|1ₖ⟩ = √(ħc²/(2ωₖV)) e(i(k⋅x − ωₖt)), has positive frequency only: the negative-frequency solutions have been given a job, not thrown away. What a real field cannot supply is a second species. With one operator per mode the quantum aₖ† adds is the one aₖ removes, chargeless and its own antiparticle; distinct antiparticles need a complex field, topic 9.03.
6Normal ordering, and bosons by construction
The vacuum, aₖ|0⟩ = 0 for every k, has energy E₀ = Σₖ½ħωₖ. With V/(2π)³ modes per unit k-space volume and |k| ≤ K, E₀/V = (ħc/(4π²))∫₀K k²√(k² + (mc/ħ)²) dk ≈ ħcK⁴/(16π²): for K = 10 fm⁻¹, 12.5 GeV fm⁻³ or 2.0 × 10³⁶ J m⁻³. Normal ordering, every aₖ† moved left of every aₖ, replaces H by :H: = Σₖħωₖaₖ†aₖ, which differs from H only by the constant E₀ and so leaves every energy difference unchanged. Only gravity responds to an absolute energy, and the dark-energy density cosmology measures, about 5 × 10⁻¹⁰ J m⁻³, is some 4 × 10⁴⁵ times smaller: a mismatch this course leaves open. When the constant's change with boundaries is measurable is Unit 11's question. Symmetry needs no repair. Creation operators commute, so aₖ₁†aₖ₂†|0⟩ = aₖ₂†aₖ₁†|0⟩, derived: the two-quantum wavefunction ⟨0|φ(x₁)φ(x₂)|1ₖ₁1ₖ₂⟩ is symmetric in x₁ and x₂, and (aₖ†)²|0⟩ = √2|2ₖ⟩ is allowed. The quanta are bosons, and the excess of identical-pion pairs at small relative momentum, measured in collisions since 1960, is that symmetry seen.
Key equations
Read the units, assumptions and limits before substituting.
Equal-time canonical commutator, imposed
[φ(x, t), π(y, t)] = iħ δ³(x − y)
A postulate, not a derivation: the ℓ → 0 limit of Quantum Mechanics I's [φⱼ, pⱼ] = iħ for each lattice cell of side ℓ, where pⱼ = ℓ³πⱼ.
π = (1/c²)∂φ/∂t; φ in √(J m⁻¹), so φπ and iħδ³ both carry J s m⁻³. Also [φ, φ] = [π, π] = 0.Mode expansion of a Hermitian field
φ(x, t) = Σₖ √(ħc²/(2ωₖV)) [aₖ e(i(k⋅x − ωₖt)) + aₖ† e(−i(k⋅x − ωₖt))]
Derived: φ = φ† forces the operator on each negative-frequency wave to be aₖ†, the adjoint of the one on the positive-frequency wave, so a real field has one operator per mode.
Periodic box, V = L³; k = 2πn/L, n a vector of integers; ωₖ = c√(k² + (mc/ħ)²) in s⁻¹; aₖ dimensionless.Ladder algebra of each mode
[aₖ, aₖ'†] = δₖₖ', [aₖ, aₖ'] = 0
Derived by inverting the expansion, aₖ = (2ħV)(−1/2) ∫d³x e(−ik⋅x)[(√ωₖ/c)φ + i(c/√ωₖ)π] at t = 0: Quantum Mechanics I's [a, a†] = 1, once per mode.
Pure numbers; δₖₖ' is a Kronecker delta because the box makes k discrete.Hamiltonian of the free real field
H = Σₖ ħωₖ (aₖ†aₖ + ½)
Derived: the aₖa₋ₖ terms cancel because ωₖ² = c²(k² + m²c²/ħ²). With P = Σₖ ħk aₖ†aₖ a quantum carries ħωₖ and ħk, so E² = p²c² + m²c⁴.
H in J; from ℋ = ½[c²π² + (∇φ)² + (mc/ħ)²φ²], topic 9.01's T⁰⁰; ħωₖ = √(ħ²k²c² + m²c⁴).Bose symmetry from the algebra
aₖ₁† aₖ₂† |0⟩ = aₖ₂† aₖ₁† |0⟩
Derived: exchange symmetry and unlimited occupation follow from the commutator, so Quantum Mechanics II's symmetrisation postulate becomes a theorem for these quanta.
From [aₖ₁†, aₖ₂†] = 0. Normalised states (aₖ†)ⁿ|0⟩/√(n!) exist for every n = 0, 1, 2, …Zero-point energy density, cutoff K
E₀/V = (1/V) Σₖ ½ħωₖ ≈ ħcK⁴/(16π²) for K ≫ mc/ħ
Derived, and divergent: K = 10 fm⁻¹ gives 2.0 × 10³⁶ J m⁻³. Normal ordering, :H: = Σₖ ħωₖaₖ†aₖ, drops it; no energy difference ever contained it.
E₀/V in J m⁻³; K, the largest |k| kept, in m⁻¹; a box holds V/(2π)³ modes per unit volume of k-space.Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. Set m = 0 for the massless limit, where the branches become the straight lines ±ħc|k|. Restore m = 135 for the π⁰ and step n from 0 to 3: the dot climbs in equal steps of ħωₖ, one oscillator's ladder. Slide k through 0: the arrow always ends at positive energy.
Live results
- ONE QUANTUM ħωₖ
- 167.2 MeV
- DROPPED ZERO POINT ½ħωₖ
- 83.6 MeV
- n QUANTA, E − E₀
- 167.2 MeV
- GROUP SPEED v/c
- 0.590
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. Mode wavenumber k: 0.50 fm⁻¹; Mass of the quantum m: 135 MeV/c²; Quanta in mode k, n: 1
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
A first reading of the expansion takes aₖ†e(−i(k⋅x − ωₖt)) to be Unit 3's negative-energy solution with an operator attached, so that aₖ† makes a quantum of energy −ħωₖ, or else an antiparticle. Neither is right. The energy of a state is an eigenvalue of H, not the sign of the frequency of the c-number function beside an operator: [H, aₖ†] = +ħωₖaₖ†, so aₖ† adds ħωₖ and, by [P, aₖ†] = ħk aₖ†, momentum +ħk, not the −ħk that e(−ik⋅x) would carry as a wavefunction. The factor e(+iωₖt) is simply how a raising operator evolves in the Heisenberg picture. And a Hermitian field has one operator per mode, so what aₖ† adds is exactly what aₖ removes: the quantum is its own antiparticle, as the chargeless π⁰ is. A distinct antiparticle needs a second operator, which only a complex field supplies, topic 9.03.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
Treat the π⁰ (mc² = 134.977 MeV) as the quantum of a free real Klein-Gordon field in a periodic cube of side L = 10.0 fm. Let a₀ belong to the k = 0 mode and a₁ to k = (2π/L, 0, 0). For the state (a₁†)²a₀†|0⟩, find the energy above the vacuum, the momentum and the norm, and check the invariant mass.
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.
- Mode 1 has k = 2π/(10.0 fm) = 0.62832 fm⁻¹, so ħkc = 197.327 MeV fm × 0.62832 fm⁻¹ = 123.984 MeV and ħω₁ = √(123.984² + 134.977²) MeV = 183.278 MeV. Mode 0 has ħω₀ = mc² = 134.977 MeV.
- With :H: = Σₖ ħωₖaₖ†aₖ the state has n₁ = 2 and n₀ = 1, so E − E₀ = (2 × 183.278 + 134.977) MeV = 501.53 MeV. Normal ordering removed E₀ from both terms of this difference, so it changes nothing here.
- P = Σₖ ħk aₖ†aₖ gives 2ħk₁ + 0 = 247.97 MeV/c along x. The k = 0 quantum carries energy mc² and no momentum.
- Norm: moving annihilators to the right with [a, a†] = 1 gives ⟨0|a₀a₁²(a₁†)²a₀†|0⟩ = 2! × 1! = 2, so the normalised state is (a₁†)²a₀†|0⟩/√2 = |n₁ = 2, n₀ = 1⟩. Two quanta share mode 1 with no objection: there is no exclusion.
- Check: the invariant mass √(501.53² − 247.97²) MeV = 435.94 MeV exceeds 3mc² = 404.93 MeV, as it must for three quanta that are not all moving together.
Answer. E − E₀ = 501.5 MeV and P = 248.0 MeV/c along x; the squared norm is 2; the invariant mass, 435.9 MeV, exceeds 3mc² = 404.9 MeV.
MediumExample 2
On a ring of circumference L = 20.0 fm, in one space dimension, a real Klein-Gordon field is φ(x, t) = Σₖ √(ħc²/(2ωₖL)) [aₖe(i(kx − ωₖt)) + aₖ†e(−i(kx − ωₖt))] with k = 2πj/L. Keep only |j| ≤ J. Find the equal-time [φ(x), π(y)], its value at x = y and its first zero for J = 50, and the J that brings the first zero within 0.05 fm.
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.
- π = (1/c²)∂φ/∂t has coefficients −i√(ħωₖ/(2c²L)), and the product of the two square roots is ħ/(2L): ωₖ, c and the mass all cancel, so the answer is the same for every mass.
- With [aₖ, aₖ'†] = δₖₖ', [φ(x), π(y)] = (iħ/(2L))Σₖ(e(iks) + e(−iks)) = (iħ/L)[1 + 2 × (sum of cos(2πjs/L) from j = 1 to J)], with s = x − y. Summed, it is (iħ/L) sin((2J + 1)πs/L)/sin(πs/L).
- At s = 0 it equals iħ(2J + 1)/L, which for J = 50 is iħ × 101/(20.0 fm) = iħ × 5.05 fm⁻¹. The first zero is at s = L/(2J + 1) = 20.0 fm/101 = 0.198 fm.
- Integrated over the ring only the j = 0 term survives, so ∫[φ(x), π(y)] dy = iħ for every J, exactly as for iħδ(x − y), whose unit here is J s m⁻¹. The truncated sum is a δ of the right weight but finite width.
- A first zero within 0.05 fm needs 2J + 1 ≥ 20.0/0.05 = 400, so J = 200: 401 modes, first zero at 20.0 fm/401 = 0.0499 fm. One [x, p] = iħ per point is a statement about all the modes together.
Answer. [φ(x), π(y)] = (iħ/L) sin((2J + 1)πs/L)/sin(πs/L), independent of m. For J = 50 it peaks at iħ × 5.05 fm⁻¹ with its first zero at 0.198 fm; J = 200, 401 modes, brings the zero to 0.0499 fm.
HardExample 3
A massless real Klein-Gordon field fills a large box in thermal equilibrium at temperature T. Starting from :H: = Σₖ ħωₖaₖ†aₖ, derive its energy density u(T), check it against the blackbody law for photons, and evaluate it at the cosmic microwave background temperature, 2.7255 K. Say why the zero-point energy plays no part.
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.
- The modes are independent oscillators, so the canonical density operator e(−:H:/kBT)/Z factorises mode by mode. For one mode, with x = ħωₖ/kBT, ⟨nₖ⟩ is the sum of n e(−nx) over the sum of e(−nx), which is 1/(eˣ − 1), every n = 0, 1, 2, … being allowed because the quanta are bosons.
- A box holds V/(2π)³ modes per unit k-space volume and ωₖ = c|k| at m = 0, so u = ∫d³k/(2π)³ ħck/(e(ħck/kBT) − 1) = ((kBT)⁴/(2π²(ħc)³)) ∫₀^∞ x³ dx/(eˣ − 1).
- The integral is π⁴/15 = 6.494, so u = π²(kBT)⁴/(30(ħc)³). Dimensions: J⁴/(J m)³ = J m⁻³.
- Photons have two polarisations, each a set of oscillators of these frequencies (Unit 11), so they carry 2u = π²(kBT)⁴/(15(ħc)³) = 4σT⁴/c with σ = π²kB⁴/(60ħ³c²) = 5.670 × 10⁻⁸ W m⁻² K⁻⁴, the measured Stefan-Boltzmann constant.
- At 2.7255 K, kBT = 3.763 × 10⁻²³ J and ħc = 3.1615 × 10⁻²⁶ J m, so u = 2.09 × 10⁻¹⁴ J m⁻³, half the 4.17 × 10⁻¹⁴ J m⁻³ in the background photons, whose spectrum FIRAS measured to follow Planck's form within 50 parts per million of its peak.
- E₀ = Σₖ½ħωₖ is the same at every temperature, so it cancels from u, the energy above the vacuum: normal ordering removes nothing that depends on T.
Answer. u = π²(kBT)⁴/(30(ħc)³) = 2σT⁴/c; at 2.7255 K, u = 2.09 × 10⁻¹⁴ J m⁻³, half the 4.17 × 10⁻¹⁴ J m⁻³ carried by the background photons.
Quick check
Select an answer and explain it before checking.
Review & continue
Canonical quantisation of a real field imposes [φ(x), π(y)] = iħ δ³(x − y) at equal times, which makes each normal mode of the field an independent quantum oscillator whose quanta are the field's particles.
- 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.