Skip to content

Magnetic Moment

A magnetic moment is the coefficient of a particle’s linear spin response to a weak magnetic field. Minimal Dirac coupling fixes one contribution; a covariant Pauli interaction supplies an independent correction. This page matches that correction to a low-energy Hamiltonian and to electromagnetic form factors. The Pauli Equation owns the magnetic-moment convention, the minimal g=2g=2 result, Zeeman splitting, and spin precession.

Required background. The Pauli Equation defines the signed-charge convention, and Minimal Coupling defines the field tensor. Helpful background. The Gordon identity on Dirac Current relates the covariant vertex to its spin response.

Take a massive charged Dirac particle with signed charge q≠0q\ne0. Parameterize its moment by

μ=gq2mS,g=2(1+a),S=ℏ2σ.\boldsymbol\mu =g\frac{q}{2m}\mathbf S,\qquad g=2(1+a),\qquad \mathbf S=\frac{\hbar}{2}\boldsymbol\sigma.

For the charged lepton convention, gg is positive and the electron sign resides in q=−eq=-e. The parameter aa is the anomaly relative to minimal Dirac coupling. For general composites the coefficient can have either sign; no positive-gg assumption should be imposed on an independently matched moment.

The local Dirac–Pauli equation with a constant real coefficient aa is

[iℏcγμDμ−mc2−aqℏ4mσμνFμν]ψ=0,\left[ i\hbar c\gamma^\mu D_\mu-mc^2 -\frac{a q\hbar}{4m} \sigma^{\mu\nu}F_{\mu\nu} \right]\psi=0,

where σμν=i[γμ,γν]/2\sigma^{\mu\nu}=i[\gamma^\mu,\gamma^\nu]/2. Its field-strength interaction is gauge invariant: FμνF_{\mu\nu} is invariant and the spinor transforms by the same phase as in minimal coupling. It does not change the charge in DμD_\mu.

This is a specified effective interaction, not a derivation of aa from a wave equation. A measured coefficient, constituent model, or underlying quantum field theory must determine it. Taking a constant coefficient also omits resolved momentum dependence and further derivative operators.

With the tensor convention of Minimal Coupling,

F0i=Eic,Fij=−ϵijkBk,F_{0i}=\frac{E_i}{c},\qquad F_{ij}=-\epsilon_{ijk}B_k,

the contraction is

σμνFμν=2(icα⋅E−Σ⋅B).\sigma^{\mu\nu}F_{\mu\nu} =2\left( \frac{i}{c}\boldsymbol\alpha\cdot\mathbf E -\boldsymbol\Sigma\cdot\mathbf B \right).

Multiplication of the equation by β=γ0\beta=\gamma^0 therefore gives the additional Hamiltonian

δH=−aqℏ2mβΣ⋅B+iaqℏ2mcβα⋅E.\delta H =-\frac{a q\hbar}{2m}\beta\boldsymbol\Sigma\cdot\mathbf B +\frac{i a q\hbar}{2mc}\beta\boldsymbol\alpha\cdot\mathbf E.

Both terms are Hermitian for real fields: βΣ\beta\boldsymbol\Sigma is Hermitian and βα\beta\boldsymbol\alpha is anti-Hermitian. The electric term is odd in the energy-block decomposition. It is part of the Lorentz-covariant magnetic interaction, not a rest-frame electric dipole energy proportional to S⋅E\mathbf S\cdot\mathbf E.

At leading low-energy order the upper block yields

δHspin=−aqℏ2mσ⋅B.\delta H_{\rm spin} =-\frac{a q\hbar}{2m}\boldsymbol\sigma\cdot\mathbf B.

Adding the minimal Dirac contribution gives

Hspin=−(1+a)qℏ2mσ⋅B,H_{\rm spin} =-\frac{(1+a)q\hbar}{2m} \boldsymbol\sigma\cdot\mathbf B,

which matches g=2(1+a)g=2(1+a). This matching assumes weak, slowly varying fields and the positive-energy nonrelativistic regime. At higher order, the odd electric term and its commutators also contribute. Changing only the Zeeman coefficient in a relativistic Hamiltonian does not include all effects of the covariant Pauli operator.

For this paragraph set ℏ=c=1\hbar=c=1. For on-shell spinors of the same mass, define the momentum transfer k=p′−pk=p'-p and write the parity-even electromagnetic vertex as

Γμ(p′,p)=q[F1(k2)γμ+iF2(k2)2mσμνkν].\Gamma^\mu(p',p) =q\left[ F_1(k^2)\gamma^\mu +\frac{iF_2(k^2)}{2m} \sigma^{\mu\nu}k_\nu \right].

The physical charge has been factored out, so F1(0)=1F_1(0)=1. With this normalization,

a=F2(0),g=2[F1(0)+F2(0)].a=F_2(0),\qquad g=2\bigl[F_1(0)+F_2(0)\bigr].

One way to see why both form factors enter the spin coefficient is the Gordon identity:

uˉ′γμu=uˉ′[(p′+p)μ2m+iσμνkν2m]u.\bar u'\gamma^\mu u =\bar u'\left[ \frac{(p'+p)^\mu}{2m} +\frac{i\sigma^{\mu\nu}k_\nu}{2m} \right]u.

After this substitution the spin-dependent part of the vertex carries F1+F2F_1+F_2, whereas the convection part carries F1F_1. The static, long-wavelength magnetic response selects the former combination. Grozin (2005) develops this form-factor interpretation and its perturbative calculation.

The local Pauli coefficient matches the zero-transfer value. At finite transfer, form factors describe additional structure and radiative response; a constant anomalous moment cannot reproduce that information. Schwinger’s leading QED correction is discussed on the Pauli owner page; its loop calculation requires quantized fields and is outside this one-particle matching argument.

Intrinsic spin response and orbital energy

Section titled “Intrinsic spin response and orbital energy”

Return to SI units. For a uniform field, choose the symmetric gauge A=B×r/2\mathbf A=\mathbf B\times\mathbf r/2. The leading low-energy orbital term expands as

(p−qA)22m=p22m−q2mB⋅L+q28m∣B×r∣2.\frac{(\mathbf p-q\mathbf A)^2}{2m} =\frac{\mathbf p^2}{2m} -\frac{q}{2m}\mathbf B\cdot\mathbf L +\frac{q^2}{8m}|\mathbf B\times\mathbf r|^2.

Thus the energy linear in a weak field contains both orbital and intrinsic contributions,

Hlinear=−q2m(L+gS)⋅B.H_{\rm linear} =-\frac{q}{2m} (\mathbf L+g\mathbf S)\cdot\mathbf B.

The quadratic term is the diamagnetic contribution. These are terms in a Hamiltonian with a specified orbital state and boundary conditions. The derivative −∂E/∂B-\partial E/\partial B of a complete energy level can include orbital motion, state mixing, and relativistic binding corrections. It need not equal a free particle’s intrinsic spin moment.

For an isolated, differentiable nondegenerate eigenvalue, the Hellmann–Feynman theorem relates that derivative to −⟨∂H/∂B⟩-\langle\partial H/\partial B\rangle. At a degeneracy, first diagonalize the perturbation in the degenerate subspace. Atomic Landé factors and material-dependent effective spin Hamiltonians therefore require the appropriate state and environment; they are not alternative measurements of the same free Dirac coefficient.

Neutral particles and independently signed moments

Section titled “Neutral particles and independently signed moments”

The parameterization aq/(2m)a q/(2m) is inconvenient when q=0q=0. Instead introduce a real coefficient μ∗\mu_* with units of magnetic moment and use

[iℏcγμDμ−mc2−μ∗2σμνFμν]ψ=0.\left[ i\hbar c\gamma^\mu D_\mu-mc^2 -\frac{\mu_*}{2}\sigma^{\mu\nu}F_{\mu\nu} \right]\psi=0.

For a neutral Dirac particle, Dμ=∂μD_\mu=\partial_\mu, but the Pauli term can remain nonzero. Its low-energy interaction is

Hspin=−μ∗σ⋅B,μ=2μ∗ℏS.H_{\rm spin}=-\mu_*\boldsymbol\sigma\cdot\mathbf B, \qquad \boldsymbol\mu=\frac{2\mu_*}{\hbar}\mathbf S.

For the additional moment of a charged particle, μ∗=aqℏ/(2m)\mu_*=a q\hbar/(2m) recovers the previous equations. For a neutral particle it is an independent signed coefficient. This observation does not assert that every neutral particle permits such a diagonal operator. For an identical anticommuting Majorana field, the tensor bilinear ΨˉMσμνΨM\bar\Psi_{\rm M}\sigma^{\mu\nu}\Psi_{\rm M} vanishes, excluding this diagonal Pauli moment (Dreiner, Haber, and Martin, 2010). This is a field-statistics statement. Majorana Spinors explains why identities for anticommuting fields must be distinguished from those for commuting classical spinor components.

Checking Hermiticity. Show that iβαii\beta\alpha_i is Hermitian. Why does the explicit ii in the electric interaction not make the Hamiltonian non-Hermitian?

Solution

Since αi†=αi\alpha_i^\dagger=\alpha_i, β†=β\beta^\dagger=\beta, and {αi,β}=0\{\alpha_i,\beta\}=0,

(iβαi)†=−iαiβ=iβαi.(i\beta\alpha_i)^\dagger =-i\alpha_i\beta=i\beta\alpha_i.

The matrix that multiplies ii was anti-Hermitian. For real multiplicative EiE_i, the product is Hermitian on the appropriate common domain.

Orbital contamination. In a weak field Bz^B\hat{\mathbf z}, suppose the perturbation is diagonal in a state with Lz=ℏmℓL_z=\hbar m_\ell and Sz=ℏmsS_z=\hbar m_s. Find its first-order slope.

Solution

The linear shift is

ΔE=−qℏB2m(mℓ+gms).\Delta E =-\frac{q\hbar B}{2m}(m_\ell+g m_s).

The slope measures the total combination. Only after determining or removing the orbital contribution can it isolate the intrinsic gg.

Neutral matching. A neutral spin-half particle has μ=γS\boldsymbol\mu=\gamma\mathbf S with signed gyromagnetic ratio γ\gamma. Determine μ∗\mu_* and explain why assigning it a finite charge-normalized anomaly is unsuitable.

Solution

Matching gives μ∗=γℏ/2\mu_*=\gamma\hbar/2. The product aqℏ/(2m)a q\hbar/(2m) vanishes for q=0q=0 at fixed finite aa, so it cannot represent a nonzero neutral moment. The independently signed coefficient avoids a singular normalization.