Skip to main content

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.
Starting model: Klein-Gordon minimal coupling and gauge invarianceEnergy against canonical momentum ħk for a spinless charge q, in units of mc² and mc. Dashed: the field-free shell. Solid: the plateau's shell, lifted by V₀ = 0.2 mc² and shifted sideways by qA₀ = 0.4 mc. At E = 1.45 mc² the plateau wave has P = 0.750 mc and ħκ = 0.000 mc; the arrow, from the gap's centre to E, is ρ/|φ|² = 1.25.dashed: field-free shell, E² = m²c⁴ + c²p²solid: plateau, V₀ = 0.2 mc², qA₀ = 0.4 mcarrow: ρ/|φ|² = (E − V₀)/mc² = 1.25EE/mc²0canonical ħk/mc, from −3.4 to +3.4
Start with a picture. Change the model in the lab below.
1

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.

2

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 = 1

Gauge-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 momentum

Local 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 = 1
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

0.2 mc²
0 to 3.5 mc²; steps of 0.1.
1.45 mc²
1.05 to 1.95 mc²; steps of 0.1.
0.4 mc
−1 to 1 mc; steps of 0.1.

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.

Interactive Advanced Quantum Mechanics modelEnergy against canonical momentum ħk for a spinless charge q, in units of mc² and mc. Dashed: the field-free shell. Solid: the plateau's shell, lifted by V₀ = 0.2 mc² and shifted sideways by qA₀ = 0.4 mc. At E = 1.45 mc² the plateau wave has P = 0.750 mc and ħκ = 0.000 mc; the arrow, from the gap's centre to E, is ρ/|φ|² = 1.25.dashed: field-free shell, E² = m²c⁴ + c²p²solid: plateau, V₀ = 0.2 mc², qA₀ = 0.4 mcarrow: ρ/|φ|² = (E − V₀)/mc² = 1.25EE/mc²0canonical ħk/mc, from −3.4 to +3.4
Energy against canonical momentum ħk for a spinless charge q, in units of mc² and mc. Dashed: the field-free shell. Solid: the plateau's shell, lifted by V₀ = 0.2 mc² and shifted sideways by qA₀ = 0.4 mc. At E = 1.45 mc² the plateau wave has P = 0.750 mc and ħκ = 0.000 mc; the arrow, from the gap's centre to E, is ρ/|φ|² = 1.25.

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

Baseline 0.2; choose 0 to 3.5, in steps of 0.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

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.

5

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.

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.

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.

6

Quick check

Select an answer and explain it before checking.

A π⁺ (mc² = 139.6 MeV) is in a stationary state of energy E = 300 MeV inside a long, field-free drift tube held at Φ₀ = 50 MV, so V = eΦ₀ = 50 MeV there. The gauge function Λ = Φ₀t sets the potential inside the tube to zero. What does the new description give?
7

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.

Extra challenging questions

Connect several ideas, explore supplied data, and compare your reasoning with a worked solution.