Belinfante-Rosenfeld Improvement
Statement
The canonical energy–momentum tensor \(T_c^{\mu\nu}\) obtained from translation invariance is conserved, \(\partial_\mu T_c^{\mu\nu}=0\), but is generally neither symmetric nor gauge invariant. Adding a divergence of a superpotential, \(T_B^{\mu\nu}=T_c^{\mu\nu}+\partial_\lambda B^{\lambda\mu\nu}\) with \(B^{\lambda\mu\nu}=-B^{\mu\lambda\nu}\), leaves the conserved charges and the conservation law untouched. Choosing the superpotential from the spin current \(S^{\lambda\mu\nu}\) fixed by Lorentz (rotation–boost) invariance, \(B^{\lambda\mu\nu}=\tfrac12\!\left(S^{\lambda\mu\nu}+S^{\mu\nu\lambda}+S^{\nu\mu\lambda}\right)\), yields the Belinfante–Rosenfeld tensor \(T_B^{\mu\nu}=T_B^{\nu\mu}\): symmetric, gauge invariant, and identical to the Hilbert tensor \(\tfrac{2}{\sqrt{-g}}\,\delta S/\delta g_{\mu\nu}\) that sources Einstein's equations.
Why it matters
Symmetry of the stress tensor is not optional book-keeping. It is the local statement of angular-momentum conservation, it is what makes \(T^{\mu\nu}\) a legitimate source in \(G^{\mu\nu}=8\pi G\,T^{\mu\nu}\) (the Einstein tensor is symmetric), and it is required for a gauge-invariant, physically meaningful momentum density. The canonical tensor for the electromagnetic field is famously asymmetric and gauge dependent; the Belinfante procedure is the systematic cure.
The construction shows that the "spin" carried by a field is not an independent bolt-on but is repackaged into the orbital moment \(x^{\mu}T_B^{\nu\lambda}-x^{\nu}T_B^{\mu\lambda}\) of a single symmetric tensor. Total angular momentum is unchanged; only its accounting is made covariant and geometric.
Assumptions
Derivation
Result
Reading. The intrinsic spin content \(S^{\lambda\mu\nu}\) is folded back into the stress tensor as a total divergence. The result is a single symmetric, conserved, gauge-invariant tensor whose orbital moment \(x^{\mu}T_B^{\nu\lambda}-x^{\nu}T_B^{\mu\lambda}\) already carries the full angular momentum — spin and orbital are unified. Because it is symmetric and gauge invariant, \(T_B^{\mu\nu}\) coincides with the Hilbert tensor obtained by varying the action with respect to the metric, and it is this tensor, not \(T_c\), that couples to gravity.
Units check. \(T_c^{\mu\nu}\) has dimension of energy density, \([\mathrm{J\,m^{-3}}]=[\mathrm{Pa}]\). The spin current \(S^{\lambda\mu\nu}\sim x\cdot T\) carries an extra length, \([\mathrm{J\,m^{-2}}]\); the derivative \(\partial_\lambda\) removes it, so \(\partial_\lambda B^{\lambda\mu\nu}\) is again \([\mathrm{J\,m^{-3}}]\), matching \(T_c^{\mu\nu}\). Dimensions are consistent.
Limiting cases
- Scalar field. \((\Sigma^{\mu\nu})=0\Rightarrow S^{\lambda\mu\nu}=0\), so \(T_B^{\mu\nu}=T_c^{\mu\nu}\): the canonical tensor is already symmetric and needs no improvement.
- Electromagnetism. \(S^{\lambda\mu\nu}=F^{\mu\lambda}A^{\nu}-F^{\nu\lambda}A^{\mu}\) gives \(\partial_\lambda B^{\lambda\mu\nu}=\partial_\lambda(F^{\mu\lambda}A^{\nu})\), converting the asymmetric, gauge-dependent \(T_c^{\mu\nu}\) into \(T_B^{\mu\nu}=-F^{\mu\alpha}F^{\nu}{}_{\alpha}+\tfrac14\eta^{\mu\nu}F_{\alpha\beta}F^{\alpha\beta}\), the symmetric gauge-invariant Maxwell tensor.
- Dirac field. The spin current \(S^{\lambda\mu\nu}=-\tfrac{i}{4}\bar\psi\{\gamma^{\lambda},\sigma^{\mu\nu}\}\psi\) is totally antisymmetric; Belinfante symmetrisation reproduces the standard \(T_B^{\mu\nu}=\tfrac{i}{4}\bar\psi(\gamma^{\mu}\!\overset{\leftrightarrow}{\partial^{\nu}}+\gamma^{\nu}\!\overset{\leftrightarrow}{\partial^{\mu}})\psi\).
- Flat-space / no gravity. Only \(P^{\nu}=\int T^{0\nu}d^3x\) and \(J^{\mu\nu}\) are observable; since improvements are pure divergences, \(T_c\) and \(T_B\) give identical total charges and are physically equivalent for global quantities.
Breaks when
- Higher-derivative or non-minimal Lagrangians. When \(\mathcal{L}\) depends on \(\partial\partial\phi\) or the symmetry is realised non-linearly, the spin current picks up additional terms and the three-term superpotential no longer symmetrises \(T_c\); the improvement must be rederived from scratch.
- Local (gauge) or conformal improvements confused with Belinfante. In gravitational or conformal contexts one may also need to add Weyl/curvature-improvement terms (e.g. \(\propto(\partial^{\mu}\partial^{\nu}-\eta^{\mu\nu}\Box)\phi^2\)) to make the trace vanish; the Belinfante procedure fixes only the antisymmetric part, not the trace, so a conformally-improved tensor requires a further, distinct step.
- Off shell / with sources. If the equations of motion fail (external currents, constraints not imposed), \(\partial_\lambda S^{\lambda\mu\nu}\neq T_c^{[\nu\mu]}\cdot 2\) and the antisymmetric parts no longer cancel, so \(T_B\) is neither symmetric nor conserved.
- Torsionful spacetime. If matter couples to an independent spin connection (Einstein–Cartan gravity), spin sources torsion and cannot be fully absorbed into a symmetric metric tensor; the Belinfante identification with the Hilbert tensor breaks and \(S^{\lambda\mu\nu}\) becomes an independent geometric source.
Failure modes
- Symmetrising by hand. Replacing \(T_c^{\mu\nu}\) with its symmetric part \(\tfrac12(T_c^{\mu\nu}+T_c^{\nu\mu})\). This changes the conserved charges and is generally not conserved — the improvement must be a specific divergence, not an ad hoc average.
- Wrong index symmetry of \(B\). Using a \(B^{\lambda\mu\nu}\) that is not antisymmetric in \(\lambda\mu\). Then \(\partial_\mu\partial_\lambda B^{\lambda\mu\nu}\neq0\) and conservation is destroyed.
- Dropping the on-shell condition. Trying to prove \(T_B^{[\mu\nu]}=0\) off shell; the cancellation in Step 8 relies on \(\partial_\lambda T_c^{\lambda\nu}=0\), which holds only on the equations of motion.
- Sign/ordering slips in the three cyclic \(S\) terms. Mis-ordering the indices of \(S^{\mu\nu\lambda}\) versus \(S^{\nu\mu\lambda}\); only the specific combination cancels the antisymmetric part — a swapped index gives a non-symmetric \(T_B\).
- Confusing spin current with spin density. Reading \(S^{\lambda\mu\nu}\) as "the spin" rather than a current whose divergence matters; the physical spin density is \(S^{0\mu\nu}\).
- Assuming EM \(T_c\) is already the physical tensor. Quoting \(T_c^{\mu\nu}=-F^{\mu\alpha}\partial^{\nu}A_{\alpha}+\ldots\) as the stress tensor — it is gauge dependent and asymmetric and gives a wrong momentum density.
Discussion
The deep content of the Belinfante–Rosenfeld construction is that a field theory does not hand you a unique local energy–momentum tensor. The canonical tensor is defined only up to the addition of an identically-conserved "improvement" divergence, an ambiguity that leaves all Poincaré charges invariant. Physics selects among these representatives by an external criterion: coupling to gravity. General relativity sources spacetime curvature with a symmetric tensor obtained by metric variation, \(T_H^{\mu\nu}=\tfrac{2}{\sqrt{-g}}\,\delta S/\delta g_{\mu\nu}\), and the theorem's punchline is that \(T_B^{\mu\nu}=T_H^{\mu\nu}\). The abstract requirement "be a good gravitational source" and the algebraic requirement "be symmetric" are the same requirement.
Symmetry is equivalent to local angular-momentum conservation. Once \(T_B\) is symmetric, the angular-momentum current is purely orbital, \(M_B^{\lambda\mu\nu}=x^{\mu}T_B^{\lambda\nu}-x^{\nu}T_B^{\lambda\mu}\), and \(\partial_\lambda M_B^{\lambda\mu\nu}=T_B^{\mu\nu}-T_B^{\nu\mu}=0\) follows automatically. The spin has not disappeared — it has been re-expressed as the moment of a redistributed momentum density. This is why the electromagnetic Poynting vector, the "correct" momentum density, emerges only after improvement: the canonical momentum density is gauge dependent and physically wrong at a point, even though its integral is right.
At the level of the covariant phase space, the improvement ambiguity is precisely the freedom to add exact forms to the Noether current — a cohomological statement. In the first-order (Palatini/tetrad) formulation the two notions genuinely split: the metric tensor \(\delta S/\delta e^{a}_{\mu}\) is Belinfante's symmetric object, while \(\delta S/\delta\omega^{ab}_{\mu}\) is the spin current, which sources torsion in Einstein–Cartan theory. There the Belinfante identity \(\partial_\lambda S^{\lambda\mu\nu}=T_c^{[\nu\mu]}\times 2\) is nothing but the on-shell content of the torsion equation, and the "improvement" that symmetrises \(T_c\) is exactly the geometric contribution of torsion integrated out. Belinfante–Rosenfeld is thus the flat-space shadow of the tetrad–spin-connection structure of gravity.
Common misconceptions. (1) "The improvement changes the physics." It does not change any measurable charge or the gravitational coupling; it changes only the unobservable local distribution, selecting the representative that is symmetric and gauge invariant. (2) "Belinfante makes the tensor traceless / conformal." No — it fixes the antisymmetric part only. Making the trace vanish (for a conformal theory) is an independent improvement. (3) "Spin is destroyed." Spin angular momentum is conserved throughout; it is merely re-bookkept from an explicit \(S^{\lambda\mu\nu}\) into the orbital moment of \(T_B\).
Worked examples
Example 1 — Electromagnetic plane wave: canonical vs symmetric momentum density.
Reading. Belinfante improvement replaces the gauge-dependent canonical momentum density with the Poynting result \(g=u/c\); the total momentum in any packet is unchanged because the correction integrates to zero.
Units check. \(\varepsilon_0 E^2\) has units \(\mathrm{(F\,m^{-1})(V\,m^{-1})^2}=\mathrm{J\,m^{-3}}\); dividing by \(c\ [\mathrm{m\,s^{-1}}]\) gives \(\mathrm{J\,s\,m^{-4}}=\mathrm{kg\,m^{-2}\,s^{-1}}\), a momentum density. Consistent.
Example 2 — Crossed static fields: field momentum from the symmetric tensor.
Reading. Although nothing moves, the crossed static fields carry a real momentum density along \(\hat{z}\). Only the symmetric, gauge-invariant Belinfante tensor exposes it; the canonical tensor hides it in a gauge artefact. This is the field momentum that balances the "hidden mechanical momentum" of the sources.
Units check. \(\varepsilon_0[\mathrm{F\,m^{-1}}]\cdot E[\mathrm{V\,m^{-1}}]\cdot B[\mathrm{T}]\); using \(\mathrm{F\,V=C}\), \(\mathrm{T=kg\,s^{-1}\,C^{-1}\,m^{0}}\)… combining gives \(\mathrm{kg\,m^{-2}\,s^{-1}}\), a momentum density. Consistent.
Problems
- Conservation of the improvement (grade: A). Show explicitly that for any \(B^{\lambda\mu\nu}=-B^{\mu\lambda\nu}\), the tensor \(T_B^{\mu\nu}=T_c^{\mu\nu}+\partial_\lambda B^{\lambda\mu\nu}\) satisfies \(\partial_\mu T_B^{\mu\nu}=\partial_\mu T_c^{\mu\nu}\).
Solution
\(\partial_\mu T_B^{\mu\nu}=\partial_\mu T_c^{\mu\nu}+\partial_\mu\partial_\lambda B^{\lambda\mu\nu}\). The operator \(\partial_\mu\partial_\lambda\) is symmetric under \(\mu\!\leftrightarrow\!\lambda\), while \(B^{\lambda\mu\nu}\) is antisymmetric under \(\lambda\!\leftrightarrow\!\mu\); the contraction of a symmetric with an antisymmetric pair vanishes. Hence \(\partial_\mu\partial_\lambda B^{\lambda\mu\nu}=0\) and \(\partial_\mu T_B^{\mu\nu}=\partial_\mu T_c^{\mu\nu}=0\). - Equal total charges (grade: A). Prove that \(P_B^{\nu}=\int d^3x\,T_B^{0\nu}\) equals \(P_c^{\nu}=\int d^3x\,T_c^{0\nu}\) for fields vanishing at spatial infinity.
Solution
\(P_B^{\nu}-P_c^{\nu}=\int d^3x\,\partial_\lambda B^{\lambda 0\nu}\). Split \(\lambda=0,i\). The \(\lambda=0\) term is \(\partial_0 B^{00\nu}=0\) because \(B^{00\nu}=-B^{00\nu}=0\) by antisymmetry in the first two indices. The remaining \(\int d^3x\,\partial_i B^{i0\nu}\) is a total spatial divergence, converting by Gauss's theorem to a surface integral at infinity, which vanishes for localised fields. Hence \(P_B^{\nu}=P_c^{\nu}\). - Scalar field (grade: B). For a real scalar with \(\mathcal{L}=\tfrac12(\partial\phi)^2-\tfrac12 m^2\phi^2\), compute \(T_c^{\mu\nu}\) and its spin current, and confirm no improvement is needed.
Solution
\(\partial\mathcal{L}/\partial(\partial_\mu\phi)=\partial^{\mu}\phi\), so \(T_c^{\mu\nu}=\partial^{\mu}\phi\,\partial^{\nu}\phi-\eta^{\mu\nu}\mathcal{L}\), which is manifestly symmetric in \(\mu\nu\). A scalar carries the trivial Lorentz representation, \((\Sigma^{\mu\nu})=0\), so \(S^{\lambda\mu\nu}=0\) and \(B^{\lambda\mu\nu}=0\). Therefore \(T_B^{\mu\nu}=T_c^{\mu\nu}\); the canonical tensor is already the Belinfante (and Hilbert) tensor. - Electromagnetic superpotential (grade: B). Using \(S^{\lambda\mu\nu}=F^{\mu\lambda}A^{\nu}-F^{\nu\lambda}A^{\mu}\), evaluate \(B^{\lambda\mu\nu}=\tfrac12(S^{\lambda\mu\nu}+S^{\mu\nu\lambda}+S^{\nu\mu\lambda})\) and show \(\partial_\lambda B^{\lambda\mu\nu}\) turns \(T_c^{\mu\nu}\) into the symmetric Maxwell tensor.
Solution
Substituting the three cyclic orderings, the terms proportional to \(A^{\lambda}\) cancel (\(F^{\mu\nu}+F^{\nu\mu}=0\)) and the survivors combine to \(B^{\lambda\mu\nu}=F^{\mu\lambda}A^{\nu}\), which is antisymmetric in \(\lambda\mu\) as required. Then \(\partial_\lambda B^{\lambda\mu\nu}=\partial_\lambda(F^{\mu\lambda}A^{\nu})=(\partial_\lambda F^{\mu\lambda})A^{\nu}+F^{\mu\lambda}\partial_\lambda A^{\nu}=F^{\mu\lambda}\partial_\lambda A^{\nu}\) on shell (\(\partial_\lambda F^{\mu\lambda}=0\)). Adding to \(T_c^{\mu\nu}=-F^{\mu\alpha}\partial^{\nu}A_{\alpha}+\tfrac14\eta^{\mu\nu}F^2\) and using \(F^{\mu\lambda}\partial_\lambda A^{\nu}-F^{\mu\alpha}\partial^{\nu}A_{\alpha}=-F^{\mu\alpha}(\partial^{\nu}A_{\alpha}-\partial_{\alpha}A^{\nu})=-F^{\mu\alpha}F^{\nu}{}_{\alpha}\) gives \(T_B^{\mu\nu}=-F^{\mu\alpha}F^{\nu}{}_{\alpha}+\tfrac14\eta^{\mu\nu}F_{\alpha\beta}F^{\alpha\beta}\), symmetric and gauge invariant. - Radiation-pressure numeric (grade: C). Sunlight at Earth delivers intensity \(I=1.36\times10^{3}\ \mathrm{W\,m^{-2}}\). Using the symmetric tensor, find the momentum-density and the radiation pressure on a perfectly absorbing surface, and state why the canonical tensor would give an ambiguous answer.
Solution
Energy density \(u=I/c=1.36\times10^{3}/3.0\times10^{8}=4.53\times10^{-6}\ \mathrm{J\,m^{-3}}\). Momentum density \(g=u/c=I/c^{2}=1.36\times10^{3}/(9.0\times10^{16})=1.51\times10^{-14}\ \mathrm{kg\,m^{-2}\,s^{-1}}\). Radiation pressure on an absorber equals the momentum flux \(T_B^{zz}=u\), so \(P_{\mathrm{rad}}=I/c=4.53\times10^{-6}\ \mathrm{Pa}\) (double for a perfect reflector, \(9.1\times10^{-6}\ \mathrm{Pa}\)). The canonical tensor \(T_c^{zz}\) is gauge dependent and asymmetric, so its momentum-flux component is not gauge invariant and cannot be equated with a measurable pressure; only \(T_B^{\mu\nu}\), which is symmetric and gauge invariant, yields the observed value.