
Advanced Quantum Mechanics course sequence · Interactive lesson
Klein-Gordon minimal coupling and gauge invariance
A charged particle feels a field in exactly one way here: through its four-momentum. This topic builds that coupling into the Klein-Gordon equation, proves the phase freedom that comes with it, and reads off two things the course needs: a density that turns negative under a high potential, and magnetism that is orbital only.
UAQM-03.04 · Lesson 4 of 9 in this unit. Position is not a completion score.Build the model
Charge couples through the kinetic momentum, and the potentials become a choice of gauge.
The key idea
Minimal coupling puts a charge q into the Klein-Gordon equation by replacing every four-momentum pμ with the kinetic four-momentum pμ − qAμ (iħ∂μ − qAμ as an operator), which leaves the equation unchanged when the potentials are gauge-transformed and φ gains the phase e(iqΛ/ħ).
Example. A π⁺ of energy 200 MeV meets a plateau where V = qΦ = 120 MeV: there ρ = (80/139.6)|φ|² = 0.57|φ|² and the wave decays over 1.7 fm. Raise V to 280 MeV and ρ = −0.57|φ|², with the same decay, although E is still +200 MeV.
The complete physical model and assumptions
Of a charge's two momenta, quantum mechanics turns the canonical p into −iħ∇, but the mass shell constrains the kinetic P = p − qA = γmv, paired with the energy E − qΦ. So the Klein-Gordon equation acquires its fields by pμ → pμ − qAμ and becomes (iħ∂/∂t − V)²φ = [c²(−iħ∇ − qA)² + m²c⁴]φ with V = qΦ. The potentials are not unique, and the equation does not care which ones are used: changing Aμ by −∂μ Λ is undone by the local phase e(iqΛ/ħ), derived here. So the energy E and the canonical momentum ħk are bookkeeping, while E − V, P, the density and the current are physics. A stationary state has density ((E − V)/mc²)|φ|², negative wherever V exceeds E: the first sight of the Klein step that Unit 8 resolves. In the non-relativistic limit the equation becomes Quantum Mechanics II's (p − qA)²/2m + V, which holds orbital magnetism with g = 1 and diamagnetism but nothing for spin, the baseline against which Unit 6's Dirac reduction finds g = 2. The model stops where its assumptions do: a point charge with no structure (the pion's measured 0.66 fm charge radius and its polarisability need extra gauge-invariant terms), an external field that is not itself quantised, and a one-particle reading that fails once V climbs past E.
1Couple through the kinetic four-momentum
Units are SI, the metric is diag(1, −1, −1, −1), xμ = (ct, r), Aμ = (Φ/c, A) and pμ = iħ∂μ (topic 3.01). Classically a charge has two momenta: the canonical p = γmv + qA, conjugate to r and quantised as −iħ∇, and the kinetic P = p − qA = γmv. The mass shell constrains the kinetic pair, (E − qΦ)²/c² − P² = m²c², quoted from relativistic Lagrangian mechanics. So the field enters by pμ → pμ − qAμ (equivalently ∂μ → Dμ = ∂μ + (iq/ħ)Aμ), and (pμ − qAμ)(pμ − qAμ)φ = m²c²φ becomes, derived, (iħ∂/∂t − V)²φ = [c²(−iħ∇ − qA)² + m²c⁴]φ with V = qΦ. Checks: qA has units C × T m = kg m s⁻¹, a momentum, and V is in joules; with Aμ = 0 it is topic 3.01's free equation. V enters squared: for constant potentials (E − V)² = c²P² + m²c⁴, so V lifts both energy branches together, which is how a potential can push a positive-energy wave onto the lower one.
2Gauge invariance, derived, and what it leaves physical
Change the potentials by any real Λ(r, t): A → A + ∇Λ and Φ → Φ − ∂Λ/∂t, which is Aμ → Aμ − ∂μ Λ and leaves the electric and magnetic fields alone. Let φ → e(iqΛ/ħ)φ. Acting on the phase, iħ∂μ produces −q ∂μ Λ, which cancels the +q ∂μ Λ coming from the −qAμ term, so Dμ φ → e(iqΛ/ħ)Dμ φ and the equation, ρ and j are unchanged. Read backwards, demanding a free local phase forces ∂μ to become Dμ, although it also permits extra gauge-invariant terms that minimal coupling leaves out. A constant Φ₀ is pure gauge: Λ = Φ₀t removes it and lowers E by qΦ₀, so the zero of energy is a choice and only differences of V between regions are physical. Around a closed loop the phases add to q∮A⋅dl/ħ = qΦB/ħ, which no Λ removes: superconducting rings trap flux in steps of h/2e = 2.07 × 10⁻¹⁵ Wb, measured in 1961, and electron holography confirmed the Aharonov–Bohm phase in 1986.
3The density and the current carry E − V
Replace ∂ by D in topic 3.02's current: jμ = (iħ/2m)[φ*Dμ φ − φ(Dμ φ)*]. It is conserved because Dμ Dμ in the equation cancels between φ and φ* just as ∂μ ∂μ did, and gauge invariant because it is built from φ* times D φ. Its parts are ρ = (iħ/2mc²)(φ*∂φ/∂t − φ∂φ*/∂t) − (V/mc²)|φ|² and j = (1/m)Re[φ*(−iħ∇ − qA)φ]. For a stationary state e(−iEt/ħ)u(r), derived, ρ = [(E − V)/mc²]|u|²: the energy that counts is the gauge-invariant E − V. Checks: the factor is dimensionless; at V = 0 it is topic 3.02's E/mc²; in the non-relativistic limit it tends to 1. For a plane wave in constant V and A, j = (P/m)|φ|², so j/ρ = c²P/(E − V), the group velocity. For a pion, atomic and nuclear fields never threaten the sign: the Coulomb potential at the centre of a lead nucleus, about 25 MeV for a π⁺, leaves a pion with 30 MeV of kinetic energy at ρ = (169.6 − 24.9)/139.6 = 1.04|φ|².
4Where V exceeds E: three regions and the Klein step
On a plateau the shell (E − V)² = c²P² + m²c⁴ sorts every wave into three regions, each derived from it. For V ≤ E − mc² it oscillates, with positive density, moving along P. For E − mc² ≤ V ≤ E + mc² it decays with ħκc = √(m²c⁴ − (E − V)²), and its density changes sign at V = E. For V ≥ E + mc² it oscillates again, with negative density and a velocity opposite to P. Complex conjugation explains the third region: φ* solves the equation with q → −q and frequency −E/ħ, so a region-three wave of charge q is a wave of charge −q with E − V = V − E ≥ mc², for which the plateau is a well. Reading that as an antiparticle waits for Unit 9, and the reflection and transmission at such a step for Unit 8. A wave arriving in region one reaches region three only by tunnelling across the gap, 2mc² of potential wide, so V must rise by 2mc² within about ħ/mc: 279 MeV over 1.41 fm for a pion, a field near 2 × 10²³ V m⁻¹, against 2.3 × 10²¹ V m⁻¹ at the surface of a lead nucleus.
5The non-relativistic limit keeps the coupling intact
Factor out the rest-energy oscillation, φ = e(−imc²t/ħ)ψ. Then (iħ∂/∂t − V)φ = e(−imc²t/ħ)(mc² + W)ψ with W = iħ∂/∂t − V, and squaring, the m²c⁴ cancels: 2mc²Wψ + W²ψ = c²(p − qA)²ψ. Where W, the local kinetic energy, is small against 2mc², drop W² to find, derived, iħ∂ψ/∂t = [(p − qA)²/2m + V]ψ: Quantum Mechanics II's Hamiltonian, with the same gauge rule, since Λ never touches the factor e(−imc²t/ħ). The conditions are |E − mc² − V| ≪ mc² and |p − qA| ≪ mc wherever ψ is appreciable, with fields that vary slowly over times of order ħ/mc². Kept to the next order, W² returns topic 3.03's correction with p replaced by p − qA. The size of what is dropped: in hydrogen's ground state, kinetic energy over 2mc² is 13.6 eV/1.022 MeV = 1.3 × 10⁻⁵, which is α²/4.
6Orbital magnetism only: the spinless baseline
In a uniform B take A = ½B × r, for which ∇⋅A = 0, and expand: (p − qA)²/2m = p²/2m − (q/2m)B⋅L + (q²/8m)(B × r)². The linear term is an orbital moment μ = (q/2m)L with g = 1; the quadratic one raises the energy, which is diamagnetism. That is all, derived: a spinless particle in an s state has no linear Zeeman shift. Singlet atomic lines test exactly this, because S = 0 silences spin: the cadmium 643.8 nm line splits into Zeeman's normal triplet, its components μB B/h apart, 14.0 GHz at 1 T, as measured. An unpaired electron does not fit. Unit 6 reduces the Dirac equation and finds (σ⋅(p − qA))²/2m = (p − qA)²/2m − (qħ/2m)σ⋅B, an extra −(q/m)S⋅B with g = 2, quoted here; the measured g = 2.00231930436 differs from 2 only by a radiative correction beyond this course. The Klein-Gordon result is the baseline that makes g = 2 a prediction rather than an input.
Key equations
Read the units, assumptions and limits before substituting.
Minimally coupled Klein-Gordon equation
(iħ ∂/∂t − V)²φ = [c²(−iħ∇ − qA)² + m²c⁴]φ
Derived from the mass shell of the kinetic four-momentum. V enters squared, so it raises both energy branches together; with Aμ = 0 it is topic 3.01's free equation.
V = qΦ in J, A in T m, q in C; covariantly (pμ − qAμ)(pμ − qAμ)φ = m²c²φGauge transformation
A → A + ∇Λ, Φ → Φ − ∂Λ/∂t, φ → e(iqΛ/ħ)φ
Derived: the phase absorbs the change of potentials. E and the canonical ħk shift; E − V, the kinetic P = p − qA, ρ and j do not.
Λ(r, t) any real function, in T m² = V s; covariantly Aμ → Aμ − ∂μ Λ, and Dμ φ → e(iqΛ/ħ)Dμ φDensity of a stationary state
ρ = [(E − V)/mc²]|φ|²
Derived from the gauge-covariant current. Negative wherever V exceeds E, even with E positive; tends to |φ|² in the non-relativistic limit.
E the total energy including mc², in J; V = qΦ static; ρ in m⁻³ when ∫ρ d³r = 1Gauge-invariant current
j = (1/m)Re[φ*(−iħ∇ − qA)φ]
Derived. Divided by ρ it gives v = c²P/(E − V): along P on the upper branch, against P where E − V ≤ −mc².
j in m⁻² s⁻¹; for a plane wave j = (P/m)|φ|², with P = ħk − qA the kinetic momentumLocal mass shell on a plateau
(E − V)² = c²P² + m²c⁴
Derived for constant V and A. The wave oscillates for V ≤ E − mc² and again for V ≥ E + mc²; between them it decays, and its density changes sign at V = E.
P in kg m s⁻¹; where |E − V| ≤ mc² the wave decays with ħκc = √(m²c⁴ − (E − V)²)Non-relativistic limit
iħ ∂ψ/∂t = [(p − qA)²/2m + V]ψ, with φ = e(−imc²t/ħ)ψ
Derived. In uniform B with A = ½B × r it is p²/2m − (q/2m)B⋅L + (q²/8m)(B × r)²: orbital g = 1 and diamagnetism, no spin term. Dirac's −(q/m)S⋅B, g = 2, is quoted.
valid for |E − mc² − V| ≪ mc² and |p − qA| ≪ mc where ψ lives; ∫|ψ|² d³r = 1Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. With E = 1.45 mc², raise V₀: the crossings close once V₀ passes 0.45 mc², the arrow turns down past 1.45, and beyond 2.45 the line cuts the lower branch, oscillating again with negative density and v against P. Slide qA₀: ħk moves, but ρ, P and v do not.
Live results
- DENSITY ρ/|φ|² ON PLATEAU
- 1.25
- KINETIC MOMENTUM P
- 0.750 mc
- CANONICAL ħk = P + qA₀
- 1.150 mc
- VELOCITY v = c²P/(E − V₀)
- 0.600 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. Plateau potential energy V₀ = qΦ₀: 0.2 mc²; Particle energy E: 1.45 mc²; Vector potential qA₀ (pure gauge): 0.4 mc
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
A positive energy does not make a positive density. Carrying topic 3.02's free result ρ = (E/mc²)|φ|² into a potential, students conclude that every solution with E above zero is an ordinary particle. But E is not gauge invariant: Λ = Φ₀t shifts it by qΦ₀ and changes nothing measurable, so no physical statement can rest on its sign. The density is [(E − V)/mc²]|φ|², and E − V is what every gauge agrees on. Under a potential that rises above E, a wave that arrived from the field-free side with positive E carries negative density. The sign flips at V = E, not at the Klein threshold V = E + mc², so even the decaying tail between them is negative. Oscillation is decided by |E − V| against mc², the sign by E − V against zero: two different tests.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
A π⁺ (charge e, mc² = 139.57 MeV) moves along x in a uniform electric field ℰ = 20 MV m⁻¹ along +x. Gauge 1 takes Φ = −ℰx and A = 0; gauge 2 takes Φ = 0 and A = −ℰt x̂. Find the Λ that turns gauge 1 into gauge 2. For the gauge-1 stationary state of energy E = 250 MeV, find ρ/|φ|² and the local kinetic momentum at x = 5.0 m in each gauge.
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.
- Gauge 2 needs A + ∇Λ = −ℰt x̂ and Φ − ∂Λ/∂t = 0. Λ = −ℰxt does both: ∇Λ = −ℰt x̂, and ∂Λ/∂t = −ℰx cancels Φ = −ℰx. Both gauges give the same field, −∇Φ − ∂A/∂t = ℰx̂. The wavefunction becomes φ′ = e(−ieℰxt/ħ)φ.
- Gauge 1: V = eΦ = −eℰx = −(20 MeV m⁻¹)(5.0 m) = −100 MeV, so E − V = 350 MeV and ρ/|φ|² = (E − V)/mc² = 350/139.57 = 2.51. The pion has gained 100 MeV of kinetic energy from the field.
- Gauge 2: Φ′ = 0, so ρ′ = (iħ/2mc²)(φ′*∂φ′/∂t − φ′∂φ′*/∂t) = (ħω′/mc²)|φ′|², with local frequency ħω′ = E + eℰx = 350 MeV at x = 5.0 m. The state is no longer stationary, since its frequency grows with x, yet ρ′/|φ′|² = 2.51 again, and |φ′| = |φ|.
- Kinetic momentum: in gauge 1, A = 0 and the local ħk is P, with cP = √(350² − 139.57²) MeV = 321 MeV. In gauge 2 the canonical ħk′ = ħk − eℰt drifts with time, but P = ħk′ − qA′ = ħk − eℰt + eℰt = ħk: the same 321 MeV/c, and v = c²P/(E − V) = 0.917c in both.
- Checks: (E − V)/mc² is dimensionless and eℰxt/ħ is a pure phase. As ℰ → 0, Λ → 0 and the ratio returns to the free E/mc² = 1.79 = γ. What differs between the gauges, the energy label and the canonical momentum, is exactly what cannot be measured.
Answer. Λ = −ℰxt, so φ′ = e(−ieℰxt/ħ)φ. At x = 5.0 m both gauges give ρ/|φ|² = 2.51, cP = 321 MeV and v = 0.917c; only the energy label and the canonical momentum differ.
MediumExample 2
A π⁺ of total energy E = 250 MeV leaves a field-free region and crosses onto a plateau at electrostatic potential Φ₀. For Φ₀ = 60 MV, 200 MV and 450 MV, find on the plateau ρ/|φ|², whether the wave oscillates or decays, its kinetic momentum or decay length, and its velocity. Then redo the 450 MV plateau for a π⁻ of the same energy.
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.
- V = eΦ₀ for the π⁺, and on the plateau (E − V)² = c²P² + m²c⁴ with mc² = 139.57 MeV. Oscillation needs |E − V| ≥ mc², that is V ≤ 110.4 MeV or V ≥ 389.6 MeV; in between the wave decays.
- 60 MV: E − V = 190 MeV, so ρ = (190/139.57)|φ|² = 1.36|φ|², cP = √(190² − 139.57²) MeV = 129 MeV and v = c²P/(E − V) = 0.679c along P: an ordinary pion, slowed from 0.830c.
- 200 MV: E − V = 50 MeV lies in the gap. ħκc = √(139.57² − 50²) MeV = 130.3 MeV, so the decay length is ħ/κ = (197.33 MeV fm)/(130.3 MeV) = 1.51 fm, and ρ = 0.358|φ|², still positive because V is below E.
- 450 MV: E − V = −200 MeV, past the threshold 389.6 MeV. The wave oscillates with cP = √(200² − 139.57²) MeV = 143 MeV, but ρ = −1.43|φ|² and v = c²P/(E − V) = −0.716c: a wave with P along +x moves towards −x, so a wave leaving the step must carry P pointing back at it.
- π⁻ on the same plateau: V = −450 MeV, E − V = 700 MeV, ρ = 5.02|φ|², cP = 686 MeV and v = 0.980c, because for the opposite charge the plateau is a well. Conjugating the π⁺ wave of the previous step also gives a charge −e wave on this well, with E − V = +200 MeV and ρ = +1.43|φ|²: the reading Unit 9 makes precise.
- Checks: every ρ/|φ|² is a ratio of energies; as Φ₀ → 0 every case returns to the free ρ = 1.79|φ|², cP = 207 MeV and v = 0.830c; and the sign of ρ changes at V = E = 250 MeV, between the second and third plateaus.
Answer. 60 MV: ρ = 1.36|φ|², cP = 129 MeV, v = 0.679c. 200 MV: decays over 1.51 fm with ρ = 0.358|φ|². 450 MV: cP = 143 MeV, ρ = −1.43|φ|², v = −0.716c. π⁻ at 450 MV: ρ = 5.02|φ|², cP = 686 MeV, v = 0.980c.
HardExample 3
Helium's ground state has L = 0 and S = 0, and precise variational calculations give ⟨r₁² + r₂²⟩ = 2.387a₀² (quoted). Treat each electron with the non-relativistic limit of the minimally coupled equation. Find helium's energy shift in B = 2.00 T and its molar susceptibility, and compare with the measured −1.88 × 10⁻⁶ cm³ mol⁻¹ (CGS). Then predict hydrogen's 1s shift at 2.00 T the same way, beside the measured splitting gμB B with g = 2.0023.
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.
- Symmetric gauge A = ½B × r, so ∇⋅A = 0. With q = −e, each electron's (p + eA)²/2m = p²/2m + (e/2m)B⋅Lᵢ + (e²/8m)(B × rᵢ)². The linear terms add to (e/2m)B⋅L, zero on L = 0, and since Lz annihilates the state they give nothing at second order either. The spinless Hamiltonian has no spin term, and with S = 0 none would act anyway.
- For B along z, (B × r)² = B²(x² + y²), with spherical average (2/3)B²r². So ΔE = (e²B²/12m)⟨r₁² + r₂²⟩ = (1.602 × 10⁻¹⁹ C)²(2.00 T)²(2.387)(5.292 × 10⁻¹¹ m)²/(12 × 9.109 × 10⁻³¹ kg) = 6.28 × 10⁻²⁹ J = 3.92 × 10⁻¹⁰ eV. Units: C² T² m² kg⁻¹ = kg m² s⁻² = J.
- The energy rises as B², so the induced moment opposes the field: μ = −∂ΔE/∂B = −(e²B/6m)⟨r₁² + r₂²⟩ = −6.28 × 10⁻²⁹ J T⁻¹ at 2.00 T. Per mole, χ = μ₀NAμ/B = −μ₀NA e²⟨r₁² + r₂²⟩/(6m) = −2.38 × 10⁻¹¹ m³ mol⁻¹, independent of B as a linear response must be.
- Divide by 4π × 10⁻⁶ to reach CGS molar units: χ = −1.89 × 10⁻⁶ cm³ mol⁻¹, against the measured −1.88 × 10⁻⁶, a 0.6% difference. For a state with S = 0 and L = 0 the spinless baseline is the whole magnetic response.
- Hydrogen's 1s also has L = 0, with ⟨r²⟩ = 3a₀², so the same Hamiltonian predicts only ΔE = e²B²(3a₀²)/(12m) = e²B²a₀²/(4m) = 7.89 × 10⁻²⁹ J = 4.93 × 10⁻¹⁰ eV at 2.00 T: one level, shifted and unsplit.
- Measured: in 2.00 T the 1s level splits in two, hyperfine structure aside, as electron-spin resonance shows, by gμB B = 2.0023 × 5.788 × 10⁻⁵ eV T⁻¹ × 2.00 T = 2.32 × 10⁻⁴ eV (56.0 GHz), 4.7 × 10⁵ times the spinless shift. The missing term is −(q/m)S⋅B with g = 2, which Unit 6 derives from the Dirac equation; helium hid it by pairing its spins.
Answer. Helium: ΔE = 3.92 × 10⁻¹⁰ eV at 2.00 T and χ = −2.38 × 10⁻¹¹ m³ mol⁻¹ = −1.89 × 10⁻⁶ cm³ mol⁻¹, matching the measured −1.88. Hydrogen 1s: a 4.93 × 10⁻¹⁰ eV shift with no splitting, against a measured 2.32 × 10⁻⁴ eV splitting.
Quick check
Select an answer and explain it before checking.
Review & continue
Minimal coupling puts a charge q into the Klein-Gordon equation by replacing every four-momentum pμ with the kinetic four-momentum pμ − qAμ (iħ∂μ − qAμ as an operator), which leaves the equation unchanged when the potentials are gauge-transformed and φ gains the phase e(iqΛ/ħ).
- 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.