Skip to content

Minimal Coupling in Wave Mechanics

Minimal coupling is the standard way a nonrelativistic charged particle couples to prescribed electromagnetic fields in wave mechanics. The ordinary momentum operator is replaced by a gauge-covariant momentum, and the scalar electric potential contributes to the energy.

Helpful background. The Time-Dependent Schrödinger Equation identifies the Hamiltonian’s dynamical role; Probability Current supplies the flow interpretation used later; Canonical Commutation Relations supports operator checks.

For a particle of mass mm and charge qq, the Hamiltonian is

H^=12m(−iℏ∇−qA)2+qΦ.\hat H = \frac{1}{2m} \left( -i\hbar\nabla-q\mathbf A \right)^2 +q\Phi .

Here Φ(r,t)\Phi(\mathbf r,t) is the scalar potential and A(r,t)\mathbf A(\mathbf r,t) is the vector potential. The time-dependent Schrödinger equation is

iℏ∂ψ∂t=[12m(−iℏ∇−qA)2+qΦ]ψ.i\hbar\frac{\partial\psi}{\partial t} = \left[ \frac{1}{2m} \left( -i\hbar\nabla-q\mathbf A \right)^2 +q\Phi \right]\psi .

This page treats Φ\Phi and A\mathbf A as classical background fields. Spin, radiation reaction, pair creation, and quantized electromagnetic fields require additional structure.

The electric and magnetic fields are obtained from the potentials by

E=−∇Φ−∂A∂t,B=∇×A.\mathbf E = -\nabla\Phi - \frac{\partial\mathbf A}{\partial t}, \qquad \mathbf B = \nabla\times\mathbf A.

The potentials are not unique. Many different pairs (Φ,A)(\Phi,\mathbf A) describe the same physical fields. This is not a defect of the formalism; it is gauge redundancy.

The Hamiltonian is nevertheless written in terms of potentials because quantum phases are sensitive to them. Magnetic fields enter local forces through B\mathbf B, but vector potentials enter the wave equation directly. The Aharonov–Bohm phase discussion shows where this distinction becomes unavoidable.

In position representation, the canonical momentum operator is

p^=−iℏ∇.\hat{\mathbf p}=-i\hbar\nabla.

Minimal coupling defines the kinetic, or mechanical, momentum operator

π^=p^−qA=−iℏ∇−qA.\hat{\boldsymbol\pi} = \hat{\mathbf p}-q\mathbf A = -i\hbar\nabla-q\mathbf A.

The Hamiltonian is then

H^=π^ 22m+qΦ.\hat H = \frac{\hat{\boldsymbol\pi}^{\,2}}{2m} +q\Phi.

The distinction matters. Canonical momentum is tied to translations and representation conventions. Kinetic momentum is tied to mechanical velocity:

v^=π^m.\hat{\mathbf v} = \frac{\hat{\boldsymbol\pi}}{m}.

In a magnetic field, the kinetic momentum components generally do not commute:

[π^i,π^j]=iℏq∑kϵijkBk.[\hat\pi_i,\hat\pi_j] = i\hbar q\sum_k\epsilon_{ijk}B_k.

This noncommutativity is one of the algebraic seeds of magnetic quantization and Landau levels.

The compact expression

(−iℏ∇−qA)2\left( -i\hbar\nabla-q\mathbf A \right)^2

must be read as an operator product. Since derivatives act on everything to their right, A(r,t)\mathbf A(\mathbf r,t) and ∇\nabla do not commute when A\mathbf A depends on position.

Acting on a wavefunction,

(p^−qA)2ψ=−ℏ2∇2ψ+iqℏ(∇⋅A)ψ+2iqℏ A⋅∇ψ+q2A2ψ.\begin{aligned} \left( \hat{\mathbf p}-q\mathbf A \right)^2\psi &= -\hbar^2\nabla^2\psi +iq\hbar(\nabla\cdot\mathbf A)\psi \\ &\quad +2iq\hbar\,\mathbf A\cdot\nabla\psi +q^2\mathbf A^2\psi . \end{aligned}

This is why it is safer to keep the gauge-covariant form until a specific gauge has been chosen. Expanding too early is a common source of sign and ordering errors.

A gauge transformation is specified by a scalar function χ(r,t)\chi(\mathbf r,t):

A′=A+∇χ,Φ′=Φ−∂χ∂t.\mathbf A' = \mathbf A+\nabla\chi, \qquad \Phi' = \Phi-\frac{\partial\chi}{\partial t}.

These changes leave the physical fields E\mathbf E and B\mathbf B unchanged. The wavefunction must transform at the same time:

ψ′=exp⁡(iqχℏ)ψ.\psi' = \exp\left( \frac{iq\chi}{\hbar} \right)\psi .

This equation is often written with the exponent displayed as iqχ/ℏiq\chi/\hbar; the exponential is a phase because χ\chi is real.

The covariant momentum transforms consistently:

(−iℏ∇−qA′)ψ′=eiqχ/ℏ(−iℏ∇−qA)ψ.\left( -i\hbar\nabla-q\mathbf A' \right)\psi' = e^{iq\chi/\hbar} \left( -i\hbar\nabla-q\mathbf A \right)\psi.

Thus the transformed Schrödinger equation has the same physical content. Gauge transformations change the description, not the measured probabilities.

The density remains

ρ=∣ψ∣2.\rho=\lvert\psi\rvert^2.

For the minimally coupled Hamiltonian, the conserved probability current is

j=1mRe⁡[ψ∗(−iℏ∇−qA)ψ].\mathbf j = \frac{1}{m} \operatorname{Re} \left[ \psi^* \left( -i\hbar\nabla-q\mathbf A \right)\psi \right].

Equivalently,

j=ℏ2mi(ψ∗∇ψ−ψ∇ψ∗)−qmA∣ψ∣2.\mathbf j = \frac{\hbar}{2mi} \left( \psi^*\nabla\psi - \psi\nabla\psi^* \right) - \frac{q}{m}\mathbf A\lvert\psi\rvert^2.

With this definition,

∂ρ∂t+∇⋅j=0.\frac{\partial\rho}{\partial t} +\nabla\cdot\mathbf j =0.

The free-particle current formula is recovered when A=0\mathbf A=0. Keeping the free-particle expression after introducing a vector potential gives a non-gauge-covariant current and generally the wrong local conservation law. The derivation and flux interpretation are developed in Probability Current and Continuity Equation.

Write

ψ(r,t)=R(r,t)eiS(r,t)/ℏ,R≥0.\psi(\mathbf r,t) = R(\mathbf r,t)e^{iS(\mathbf r,t)/\hbar}, \qquad R\ge0.

Then the density is ρ=R2\rho=R^2, and the current becomes

j=ρm(∇S−qA).\mathbf j = \frac{\rho}{m} \left( \nabla S-q\mathbf A \right).

The gauge-invariant local velocity field is therefore controlled by

∇S−qA,\nabla S-q\mathbf A,

not by ∇S\nabla S alone. Under a gauge transformation,

S′=S+qχ,A′=A+∇χ,S'=S+q\chi, \qquad \mathbf A'=\mathbf A+\nabla\chi,

so

∇S′−qA′=∇S−qA.\nabla S'-q\mathbf A' = \nabla S-q\mathbf A.

This is the simplest way to see why the vector potential term in the current is not optional.

If A=0\mathbf A=0 and Φ\Phi is time independent, minimal coupling reduces to the ordinary scalar-potential Hamiltonian

H^=p^22m+qΦ(r).\hat H = \frac{\hat{\mathbf p}^2}{2m} +q\Phi(\mathbf r).

The potential energy is V=qΦV=q\Phi.

For a uniform magnetic field

B=Bz^,\mathbf B=B\hat{\mathbf z},

two common choices are the Landau gauge

A=Bx y^\mathbf A=Bx\,\hat{\mathbf y}

and the symmetric gauge

A=12B×r.\mathbf A = \frac{1}{2}\mathbf B\times\mathbf r.

Both give the same magnetic field:

∇×A=Bz^.\nabla\times\mathbf A = B\hat{\mathbf z}.

The wavefunctions and conserved-looking quantum numbers differ between gauges, but gauge-invariant energies and measurable densities agree. This is why magnetic-field problems should be phrased in terms of fields, fluxes, currents, and kinetic momenta, not only in terms of one chosen vector potential.

  • Writing qΦq\Phi but forgetting the −qA-q\mathbf A shift in momentum.
  • Treating p^\hat{\mathbf p} and π^\hat{\boldsymbol\pi} as the same observable in a magnetic field.
  • Expanding (p^−qA)2\left(\hat{\mathbf p}-q\mathbf A\right)^2 as if ∇\nabla did not act on A\mathbf A.
  • Calling a gauge-dependent wavefunction phase unphysical without checking gauge-invariant interference effects.
  • Using the free-particle probability current after turning on A\mathbf A.
  • Thinking that a particular gauge choice is the physical magnetic field.
  1. Verify the gauge covariance of the kinetic momentum operator.
Solution

Let

ψ′=eiqχ/ℏψ,A′=A+∇χ.\psi'=e^{iq\chi/\hbar}\psi, \qquad \mathbf A'=\mathbf A+\nabla\chi.

Then

−iℏ∇ψ′=eiqχ/ℏ(−iℏ∇ψ+q∇χ ψ).-i\hbar\nabla\psi' = e^{iq\chi/\hbar} \left( -i\hbar\nabla\psi +q\nabla\chi\,\psi \right).

Therefore

(−iℏ∇−qA′)ψ′=eiqχ/ℏ(−iℏ∇ψ+q∇χ ψ−qAψ−q∇χ ψ)=eiqχ/ℏ(−iℏ∇−qA)ψ.\begin{aligned} \left( -i\hbar\nabla-q\mathbf A' \right)\psi' &= e^{iq\chi/\hbar} \left( -i\hbar\nabla\psi +q\nabla\chi\,\psi -q\mathbf A\psi -q\nabla\chi\,\psi \right) \\ &= e^{iq\chi/\hbar} \left( -i\hbar\nabla-q\mathbf A \right)\psi . \end{aligned}
  1. Starting from ψ=ReiS/ℏ\psi=Re^{iS/\hbar}, derive the current j=ρ(∇S−qA)/m\mathbf j=\rho(\nabla S-q\mathbf A)/m.
Solution

The covariant momentum acting on ψ\psi is

(−iℏ∇−qA)ψ=eiS/ℏ[−iℏ∇R+R∇S−qRA].\left( -i\hbar\nabla-q\mathbf A \right)\psi = e^{iS/\hbar} \left[ -i\hbar\nabla R +R\nabla S -qR\mathbf A \right].

Multiplying by ψ∗=Re−iS/ℏ\psi^*=Re^{-iS/\hbar} gives

ψ∗(−iℏ∇−qA)ψ=−iℏR∇R+R2(∇S−qA).\psi^* \left( -i\hbar\nabla-q\mathbf A \right)\psi = -i\hbar R\nabla R +R^2 \left( \nabla S-q\mathbf A \right).

The first term is imaginary, so the real part is

ρ(∇S−qA).\rho \left( \nabla S-q\mathbf A \right).

Dividing by mm gives the stated current.

  1. Check that the Landau gauge A=Bx y^\mathbf A=Bx\,\hat{\mathbf y} gives B=Bz^\mathbf B=B\hat{\mathbf z}.
Solution

In Cartesian components,

Ax=0,Ay=Bx,Az=0.A_x=0, \qquad A_y=Bx, \qquad A_z=0.

The curl has

Bz=∂Ay∂x−∂Ax∂y=B−0=B.B_z = \frac{\partial A_y}{\partial x} - \frac{\partial A_x}{\partial y} =B-0=B.

The other components vanish:

Bx=0,By=0.B_x=0, \qquad B_y=0.

Thus

∇×A=Bz^.\nabla\times\mathbf A = B\hat{\mathbf z}.

These pages develop the current, gauge symmetry, magnetic models, and applications:

  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.