Pseudo-gauge transformations are redefinitions of local energy–momentum and spin tensors using total divergences that preserve global conservation laws while altering local distributions.
They play a critical role in distinguishing the split between orbital and spin angular momentum and ensuring thermodynamic consistency in relativistic spin hydrodynamics.
Recent advances use pseudo-gauge invariant local equilibrium operators to remove ambiguities in spin-polarization observables, enabling more robust hydrodynamic modeling in high-energy collisions.
Pseudo-gauge transformations are redefinitions of the local energy–momentum and spin tensors by total divergences or superpotentials that preserve the conservation laws for energy, linear momentum, and angular momentum while changing the local redistribution of these quantities. In relativistic field theory this arbitrariness is described as pseudogauge freedom or symmetry, and it becomes operationally important whenever observables are computed from local densities rather than global charges. Recent work has shown that pseudo-gauge freedom controls the split between orbital and spin angular momentum, the form of local-equilibrium density operators, and the thermodynamic consistency of relativistic spin hydrodynamics (Dey et al., 2023, Becattini et al., 12 Jul 2025).
1. Formal definition and conservation structure
A commonly used formulation starts from a stress–energy tensor Tμν and a spin tensor Sλ,μν=Sλ,[μν] obeying
∂μTμν=0,∂λSλ,μν=−2T[μν].
The associated total angular-momentum current is
Jλ,μν=xμTλν−xνTλμ+Sλ,μν,
whose conservation follows from the Ward identities. A pseudo-gauge transformation is generated by an arbitrary operator-valued superpotential Φα,μν=Φα,[μν] through
S′λ,μν=Sλ,μν−Φλ,μν,
T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).
These transformed tensors satisfy the same Ward identities, and the total charges
A more general formulation allows an additional rank-4 superpotential Zμν,λρ with
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,
so that
Sλ,μν=Sλ,[μν]0
Under standard boundary conditions, Sλ,μν=Sλ,[μν]1 and Sλ,μν=Sλ,[μν]2 vanish after reduction to surface integrals at spatial infinity (Dey et al., 2023).
The physical content of the transformation is therefore local rather than global. What changes are the local densities and currents; what does not change are the exactly conserved global generators.
2. Canonical, Belinfante, GLW, HW, and related realizations
For microscopic field theories, the canonical Noether currents provide the starting point. In flat Sλ,μν=Sλ,[μν]3 spacetime, a Lagrangian density Sλ,μν=Sλ,[μν]4 yields the canonical energy–momentum tensor
Sλ,μν=Sλ,[μν]5
and canonical spin current
Sλ,μν=Sλ,[μν]6
while the Belinfante tensor is
Sλ,μν=Sλ,[μν]7
which is symmetric on the equations of motion. For free Dirac and Proca fields, explicit expressions for canonical, Belinfante, Hilgevoord–Wouthuysen, de Groot–van Leeuwen–van Weert, and Klein–Gordon-type pseudo-gauges are available, and the same framework extends to interacting theories and electromagnetic backgrounds (Armas et al., 20 Jan 2026, Weickgenannt et al., 2022).
Pseudogauge
Defining potentials
Resulting property
Canonical (Noether)
Sλ,μν=Sλ,[μν]8
Original Noether Sλ,μν=Sλ,[μν]9 and ∂μTμν=0,∂λSλ,μν=−2T[μν].0
For free fields, GLW and HW provide symmetric energy–momentum tensors and conserved spin tensors. In interacting or out-of-equilibrium settings, this simplification generally fails: interactions can spoil exact spin conservation, and one has Jλ,μν=xμTλν−xνTλμ+Sλ,μν,3 (Weickgenannt et al., 2022). In electromagnetic backgrounds, a mixed KG pseudo-gauge for matter combined with a Belinfante shift for the electromagnetic field yields a manifestly gauge-invariant splitting
Jλ,μν=xμTλν−xνTλμ+Sλ,μν,4
with Lorentz-force and spin-precession equations following for the matter sector (Weickgenannt et al., 2022).
These constructions exhibit the main structural feature of pseudo-gauge freedom: several locally inequivalent tensor decompositions can encode the same total Jλ,μν=xμTλν−xνTλμ+Sλ,μν,5 and Jλ,μν=xμTλν−xνTλμ+Sλ,μν,6.
3. Pseudo-gauge-invariant local equilibrium
A central recent development is the construction of a local thermodynamic equilibrium density operator that is itself pseudo-gauge invariant. Seeking a density operator of the form
Jλ,μν=xμTλν−xνTλμ+Sλ,μν,7
linear in Jλ,μν=xμTλν−xνTλμ+Sλ,μν,8 and Jλ,μν=xμTλν−xνTλμ+Sλ,μν,9, one considers the ansatz
Φα,μν=Φα,[μν]0
Requiring invariance under arbitrary Φα,μν=Φα,[μν]1 fixes
Φα,μν=Φα,[μν]2
Identifying Φα,μν=Φα,[μν]3, with Φα,μν=Φα,[μν]4 the inverse-temperature four-vector, and defining
Φα,μν=Φα,[μν]5
one obtains the unique pseudo-gauge-invariant local-equilibrium operator
Φα,μν=Φα,[μν]6
At global equilibrium, Φα,μν=Φα,[μν]7 is a Killing field, Φα,μν=Φα,[μν]8, and the operator reduces to
Φα,μν=Φα,[μν]9
A naïve addition of a conserved charge term,
S′λ,μν=Sλ,μν−Φλ,μν,0
is not in general pseudo-gauge invariant unless S′λ,μν=Sλ,μν−Φλ,μν,1 is constant or S′λ,μν=Sλ,μν−Φλ,μν,2 has a suitable improvement that drops out of the hypersurface integral (Becattini et al., 12 Jul 2025).
The principal consequence is that expectation values
S′λ,μν=Sλ,μν−Φλ,μν,3
become independent of the chosen split between orbital and spin angular momentum. Ambiguities in local energy flow and spin density drop out, and for all local operators built from fields and their derivatives one obtains exactly the same result as in the Belinfante pseudo-gauge. In the Belinfante frame, where S′λ,μν=Sλ,μν−Φλ,μν,4 and S′λ,μν=Sλ,μν−Φλ,μν,5 is symmetric,
S′λ,μν=Sλ,μν−Φλ,μν,6
yet the correlators coincide with those computed from the full pseudo-gauge-invariant operator (Becattini et al., 12 Jul 2025).
This construction is directly motivated by high-energy nuclear collisions. After a short time S′λ,μν=Sλ,μν−Φλ,μν,7, the system is often assumed to have locally thermalized on the hyperboloidS′λ,μν=Sλ,μν−Φλ,μν,8. Replacing the true pseudo-gauge-independent quantum state by S′λ,μν=Sλ,μν−Φλ,μν,9 on T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).0 removes the unphysical dependence of spin-polarization observables on the choice of spin tensor. As a corollary, mean spin-polarization calculations performed in the Belinfante pseudo-gauge are identified as the unique physically correct result (Becattini et al., 12 Jul 2025).
4. Thermodynamic pseudo-gauges and second-order hydrodynamics
In relativistic hydrodynamics with spin, one introduces the spin chemical potential T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).1 and the free-energy current
T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).2
Ideal spin hydrodynamics posits
T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).3
with constitutive relations
T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).4
and thermodynamic identities
T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).5
However, if one computes improved currents T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).6 in an arbitrary pseudo-gauge and then defines the rest-frame densities from T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).7, the standard first-law and Gibbs–Duhem relations need not hold. This leads to the notion of a smaller family of “thermodynamic” pseudo-gauges in which the ideal-fluid form is retained (Armas et al., 20 Jan 2026).
At quadratic order in spin, these thermodynamic pseudo-gauges are parameterized by an arbitrary function T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).8. Under such a shift, the pressure transforms as
T′μν=Tμν+21∂α(Φα,μν+Φμ,να+Φν,μα).9
Thus the spin equation of state retains a residual one-parameter ambiguity even within the thermodynamic class; in conformal theories this ambiguity is fixed by scale invariance. By contrast, combinations such as
Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν0
are pseudo-gauge invariant. Explicit free-field examples were worked out for free massless Dirac fermions and scalar fields (Armas et al., 20 Jan 2026).
A complementary hydrodynamic result is that extending hydrodynamics by a spin variable is equivalent to modifying conventional symmetric hydrodynamics by a set of non-dissipative, second-order terms. In a parity-even sector one may write
Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν1
Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν2
Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν3
The second law requires five entropy-production coefficients Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν4 to vanish, which reduces the seven transport functions Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν5 to two independent ones. These may be taken as Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν6 and Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν7, equivalently Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν8 and the spin susceptibility Pν=∫ΣdΣμTμν,Jμν=∫ΣdΣλJλ,μν9. In the same framework one obtains the nondissipative heat-current relation
Zμν,λρ0
which is described as a vorticity-driven thermal Hall effect in the no-drag frame (Li et al., 2020).
5. Dynamical constraints and symmetry-restricted residual freedom
Classical pseudo-gauge transformations in hydrodynamic models can be analyzed by decomposing the superpotential Zμν,λρ1 into Lorentz-invariant pieces relative to the fluid velocity Zμν,λρ2. Writing
Zμν,λρ3
and separating a fully transverse part
Zμν,λρ4
one arrives at a decomposition with Zμν,λρ5 independent components, matching a general Zμν,λρ6 (Drogosz et al., 2024).
If both the original and transformed energy–momentum tensors are required to be symmetric, then the superpotential must satisfy
Zμν,λρ7
the STS condition. When Zμν,λρ8 is constructed only from hydrodynamic fields Zμν,λρ9 and without higher gradients, the analysis implies that essentially only
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,0
survives. Even then, the six independent components of the STS condition become six equations for one scalar function Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,1, which generically has no solution unless further symmetries are present (Drogosz et al., 2024).
A special case is Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,2-dimensional boost-invariant Bjorken flow, where Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,3 holds identically for any function of proper time. The residual pseudo-gauge transformation is
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,4
For the hydrodynamic tensor
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,5
this residual transformation shifts the viscosities according to
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,6
while preserving the combination
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,7
Hence the bulk and shear viscosity coefficients are pseudo-gauge dependent in this symmetry class, whereas the linear combination entering the Bjorken equation of motion is invariant (Drogosz et al., 2024).
This suggests that pseudo-gauge transformations can alter local transport parametrizations without modifying the hydrodynamic structure that actually governs the evolution.
6. Operator algebra, kinetic theory, and broader usage
Pseudo-gauge freedom is not innocuous at the level of local spin operators. For equal-time spin generators
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,8
the canonical spin density yields
Zμν,λρ=−Zνμ,λρ=−Zμν,ρλ,9
By contrast, in the GLW and HW gauges the transformed spin operators acquire improvement terms whose commutators do not vanish in general, so the resulting operators fail to satisfy the Sλ,μν=Sλ,[μν]00 algebra. This is the basis for the conclusion that only the canonical spin tensor defines bona fide spin-angular-momentum operators. At the same time, polarization vectors constructed through the Pauli–Lubanski vector can still be pseudogauge-independent, and GLW/HW remain useful in classical or semiclassical settings or whenever one prefers a conserved spin tensor (Dey et al., 2023).
In Wigner-function kinetic theory, different pseudo-gauges correspond to different reorganizations of the microscopic Clifford components Sλ,μν=Sλ,[μν]01, Sλ,μν=Sλ,[μν]02, Sλ,μν=Sλ,[μν]03, and related quantities into energy–momentum and spin densities. For interacting Dirac and Proca fields, the gradient-expanded Boltzmann equation
Sλ,μν=Sλ,[μν]04
shows that nonlocal collisions mix orbital and spin degrees of freedom and spoil naïve spin conservation in the canonical pseudo-gauge, while HW, GLW, or KG-type choices absorb the same transfer differently. In global equilibrium one recovers a thermal-vorticity-driven polarization with
The term “pseudo-gauge” also appears in a mathematically distinct setting: time-dependent diffeomorphisms Sλ,μν=Sλ,[μν]07 acting on autonomous ordinary differential equations. In that usage, the transformed system is
Sλ,μν=Sλ,[μν]08
and the linear gauge transformation Sλ,μν=Sλ,[μν]09 is the special case Sλ,μν=Sλ,[μν]10. Solution correspondence, symmetry transfer, and transformation of first integrals persist, but the identification problem becomes a PDE for Sλ,μν=Sλ,[μν]11 rather than the matrix ODESλ,μν=Sλ,[μν]12 of the linear case (Gaeta et al., 8 Jun 2025).
Across these contexts, pseudo-gauge transformations encode the same abstract pattern: they preserve the underlying dynamical or conserved content while reshuffling its local representation. In relativistic spin physics, the contemporary emphasis has shifted from the mere existence of this freedom to the identification of invariant quantities, admissible thermodynamic gauges, and operator constructions that remain physically meaningful under it.