Perfect Spin Hydrodynamics: Theory & Applications

• Perfect spin hydrodynamics is a framework that extends conventional fluid dynamics by incorporating intrinsic spin via a rank-3 tensor and modified conservation laws.
• It employs a generalized thermodynamic description and both kinetic and Wigner-function approaches to derive self-consistent constitutive relations for spin-polarized observables.
• The theory finds applications in modeling heavy-ion collisions and magnetized plasmas, providing predictions for spin polarization measurements under experimental conditions.

Perfect spin hydrodynamics is a generalization of relativistic perfect-fluid theory that systematically incorporates intrinsic spin degrees of freedom, allowing for a dynamical and non-dissipative evolution of spin polarization alongside energy, momentum, and conserved charges. By introducing a rank-3 spin tensor and promoting the spin-polarization tensor to an independent thermodynamic variable, this framework enables the self-consistent prediction of spin-polarized observables in relativistic multi-component media, notably relevant for heavy-ion collision phenomenology and strongly-interacting matter.

1. Fundamental Structure: Conservation Laws and Constitutive Relations

Perfect spin hydrodynamics extends the standard set of conservation laws to include the separate local conservation of the spin part of the total angular momentum. The key mathematical objects are the particle (baryon) current $N^\mu$, the energy–momentum tensor $T^{\mu\nu}$, and the spin tensor $S^{\lambda,\mu\nu}$, where the latter carries the intrinsic angular momentum density and is antisymmetric in its last two indices.

The general conservation equations are: $$\partial_\mu N^\mu = 0, \qquad \partial_\mu T^{\mu\nu} = 0, \qquad \partial_\lambda S^{\lambda,\mu\nu} = 0$$ when $T^{\mu\nu}$ is symmetric, as in the de Groot–van Leeuwen–van Weert (GLW) pseudo-gauge formalism (Singh, 2020). The total angular momentum conservation splits into conservation laws for orbital angular momentum (involving $T^{[\mu\nu]}$) and for spin ($S^{\lambda,\mu\nu}$).

The constitutive relations for locally equilibrated spin-1/2 fluids read: $$N^\mu = n u^\mu, \qquad T^{\mu\nu} = (\varepsilon + P) u^\mu u^\nu - P g^{\mu\nu},$$

$$S^{\lambda,\mu\nu} = {\cal C} \left[ n_{(0)}\,u^\lambda\,\omega^{\mu\nu} + S_\Delta^{\lambda,\mu\nu} \right]$$

with $u^\mu$ the fluid four-velocity ($u^2=1$), $n$ the charge density, $\varepsilon$ the energy density, $P$ the pressure, ${\cal C} = \cosh(\xi)$, $\xi = \mu / T$, and $S_\Delta^{\lambda,\mu\nu}$ capturing all symmetry-allowed orthogonal structures in spin space (Singh, 2020, Ryblewski et al., 2020).

2. Thermodynamic Potentials, Spin Tensors, and Generalized First Laws

Spin hydrodynamics requires an extended thermodynamic description. The spin-polarization tensor $\omega^{\mu\nu}$, antisymmetric and six-component, plays the role of a spin chemical potential—its conjugate is the spin density degrees of freedom encoded in $S^{\lambda,\mu\nu}$. The generalized thermodynamic identities are: $$\varepsilon + P = T s + \mu n + \frac{1}{2} \Omega_{\mu\nu} S^{\mu\nu}$$

$$d\varepsilon = T ds + \mu dn + \frac{1}{2} \Omega_{\mu\nu} dS^{\mu\nu}$$

where $\Omega_{\mu\nu}=T\,\omega_{\mu\nu}$ and $S^{\mu\nu}$ denotes the local spin density in the comoving frame (Florkowski et al., 2024, Bhadury et al., 2 Jul 2025).

A fully consistent framework emerges by deriving these identities from kinetic theory or Wigner-function approaches and constructing all currents as derivatives of a single scalar generating function ("Massieu potential") $\chi(\alpha, \beta_\mu, \tfrac12\Omega_{\mu\nu})$ with respect to the conjugate Lagrange multipliers, ensuring closure of the system and the correct entropy current (Florkowski et al., 2024, Abboud et al., 24 Jun 2025).

3. Microscopic Foundations: Classical and Quantum Descriptions

Two principal microscopic realizations underlie perfect spin hydrodynamics:

• Classical spin approach: Introduces classical spin vectors $s^\mu$ or spin tensors $s^{\mu\nu}$, extending phase space and leading to distribution functions of the form $$f^\pm(x, p, s) = \exp[\pm \xi - p\cdot\beta + \tfrac12\,\omega_{\mu\nu} s^{\mu\nu}]$$. All macroscopic conserved currents are computed as moments over $p$ and $s$ (Drogosz et al., 2 Jun 2025).
• Quantum (Wigner-function) approach: Describes spin-1/2 matter via matrix-valued phase-space densities and constructs the generating function as $$n(x) = 4 \int dP \cosh\xi\, e^{-p\cdot\beta} \cosh\sqrt{-a^2}$$, where $a_\mu \propto \widetilde\omega_{\mu\nu}p^\nu$ with $\widetilde\omega$ the dual of $\omega$ (Drogosz, 7 Sep 2025).

Both approaches yield formally identical tensorial structures for $N^\mu$, $T^{\mu\nu}$, and $S^{\lambda,\mu\nu}$, differing only in prefactor normalization at higher orders in $\omega_{\mu\nu}$. Agreement at $O(\omega^2)$ ensures universality for small spin polarization regimes (Drogosz, 7 Sep 2025, Drogosz et al., 2 Jun 2025).

4. Dynamical Equations and Solutions: Evolution of Spin and Polarization

The equations of motion project into scalar and vector components along $u^\mu$ and transverse to it. For the baryon current and energy-momentum conservation, the projected equations are: $$D n + n\,\theta = 0, \quad D \varepsilon + (\varepsilon + P)\theta = 0, \quad (\varepsilon + P) D u^\alpha = \nabla^\alpha P$$ where $D = u^\mu \partial_\mu$ and $\theta = \partial_\mu u^\mu$.

The evolution of the spin-polarization tensor follows from

$$u\cdot\partial \omega^{\mu\nu} + \omega^{\mu\nu}\theta + 2 u^{[\mu} \omega^{\nu]\alpha} D u_\alpha = 0$$

so fluid acceleration and expansion drive the dynamics of local spin polarization (Singh, 2020).

Under boost-invariant (Bjorken) symmetry, all fields depend on proper time $\tau$ only, and the conservation equations reduce to a set of decoupled $1/\tau$-damped ODEs for each spin degree of freedom: $$\frac{dC}{d\tau} + \frac{C}{\tau} = 0 \implies C(\tau) = C_0 \frac{\tau_0}{\tau}$$ where $C$ are coefficients in a basis decomposition of $\omega^{\mu\nu}$ (Singh, 2020, Singh, 2020).

At freeze-out, the mean polarization vector measured in experiments is computed from moments of the hydrodynamic fields and $\omega^{\mu\nu}$ on the freeze-out hypersurface: $$E_p \frac{d\Pi_\mu(p)}{d^3p} = -\frac{\cosh \xi}{(2\pi)^3 m} \int_{\Sigma_f} d\Sigma_\lambda\, p^\lambda e^{-\beta\cdot p} \widetilde\omega_{\mu\nu} p^\nu$$ with the invariant momentum spectrum in the denominator (Singh, 2020).

5. Applicability, Validity Regime, and Causality

The regime of validity is set by the requirement that local-equilibrium integrals remain convergent and the kinetic (or Wigner) distribution is positive-definite. Defining the “electric-like” and “magnetic-like” components of $\omega^{\mu\nu}$ as $e^i, b^i$, one finds for a spin-1/2 fluid in its local rest frame: $$\text{Classical:}\quad s \sqrt{ e'^2 + b'^2 + 2|e' \times b'| } < \frac{m}{T}$$

$$\text{Quantum:}\quad \frac{1}{2} \sqrt{ e'^2 + b'^2 + 2|e' \times b'| } < \frac{m}{T}$$

with $s$ the (classical) spin length (Drogosz et al., 2 Jun 2025). Violation signals a breakdown of the perfect spin fluid description and the need for inclusion of spin-dissipative corrections.

All properly constructed formulations (classical or quantum, Boltzmann or Fermi-Dirac statistics) can be cast as divergence-type (symmetric-hyperbolic) theories (Bhadury et al., 24 Nov 2025, Abboud et al., 24 Jun 2025). For all such systems with the exact, nonperturbative generating function, nonlinear causality (signal speeds $|v|\leq 1$) and linear stability are rigorously guaranteed. Truncations at finite order in $\omega$ can break this property, making the nonperturbative structure essential at large spin polarization (Bhadury et al., 24 Nov 2025).

6. Extensions, Feedback Effects, and Applied Scenarios

The theory admits systematic second-order corrections in $\omega$ that feedback on the evolution of the fluid background. The energy–momentum tensor and current acquire terms quadratic in spin, leading to additional contributions in the hydrodynamic equations and modified scaling relations for the background fields. These feedback effects remain quantitatively small as long as $|\omega^{\mu\nu}| \ll 1$ but introduce constraints on permissible spin configurations under Bjorken symmetry and other symmetries (Drogosz et al., 2024).

Applications include:

• Modeling the time evolution of spin polarization in heavy-ion collisions and predicting $\Lambda$-hyperon polarization observables. Successful phenomenology requires delayed initialization of spin hydrodynamics relative to the background fluid, reflecting the timescale for suppression of spin-orbit dissipation and the onset of (quasi)ideal angular momentum conservation in the later stages of collisions (Singh et al., 2024).
• Extension to magnetized plasmas (guiding-center spin hydrodynamics), where motion perpendicular to strong magnetic fields is suppressed and the field structure modifies the spin dynamics and polarization evolution (Singh, 29 Mar 2025).
• Inclusion of electromagnetic field coupling, where the joint evolution of spin and electromagnetic fields yields observable effects such as polarization splitting in $\Lambda$–$\bar{\Lambda}$ measurements (Singh et al., 2022).

7. Synthesis and Outlook

Perfect spin hydrodynamics provides a symmetry-driven, thermodynamically consistent, and computationally tractable relativistic field theory for spinful matter. Its fundamental structure as a divergence-type theory ensures causal, stable evolution as long as local spin polarization is moderate. The explicit connection between microscopic (kinetic/Wigner) and macroscopic (hydrodynamic) formulations is now well-established, with complete closure relations and all-orders generating functions constructed (Florkowski et al., 2024, Drogosz, 7 Sep 2025, Abboud et al., 24 Jun 2025).

Future developments and open directions include:

• Systematic construction and implementation of dissipative (spin-relaxing) extensions suitable for the early time, high-polarization regimes of heavy-ion collisions and neutron star interiors.
• Detailed numerical simulations of spin hydrodynamics with realistic equations of state, freeze-out procedures, and feedback on global observables.
• Cross-coupling with electromagnetic fields, vorticity, chiral effects, and generalized anomaly-driven transport in quantum relativistic media.

The mathematical and physical infrastructure now enables precision modeling of spin-polarized observables in a variety of settings, with experimental applications and further formal generalizations anticipated.