⚛️ nuclear theory

Second-order evolution of spin shear tensor for massive spin-1/2 particles

Using a semi-classical kinetic theory within the relaxation time approximation, this paper derives a second-order relaxation-type evolution equation for the spin shear tensor in spin-polarized fluids of massive spin-1/2 particles, revealing new transport coefficients that govern its dynamics driven by thermodynamic forces and high-order gradient corrections.

Original authors: Pushpa Panday, Samapan Bhadury, Sadhana Dash

Published 2026-11-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Pushpa Panday, Samapan Bhadury, Sadhana Dash

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the extreme heat and density created when heavy atomic nuclei crash into one another at nearly the speed of light, matter behaves in ways that defy everyday experience. This state, known as a quark-gluon plasma, is a fluid so perfect that it flows with almost no resistance, and it is the only place in the universe where we can recreate the conditions that existed microseconds after the Big Bang. For years, physicists have been trying to understand how the tiny, intrinsic spins of the particles within this fluid align with the swirling motion of the collision itself. When two heavy ions collide off-center, they generate a massive amount of rotational motion, much like a figure skater pulling in their arms to spin faster. This rotation creates a kind of fluid turbulence that can force the particles inside to line up their internal magnetic axes, a phenomenon called polarization. While scientists have successfully measured how particles align with the overall rotation of the collision, a stubborn mystery has persisted regarding how they align along the direction of the beam itself. Theories that worked well for the overall rotation failed to explain this specific alignment, suggesting that the current models of how this fluid moves and spins were missing a crucial piece of the puzzle.

To solve this, a team of researchers has developed a more sophisticated mathematical framework that treats the spin of particles not just as a static property, but as a dynamic force that evolves over time. In their study, they focused on a specific, often overlooked aspect of the fluid's motion: the way the fluid shears, or slides past itself, and how that sliding motion interacts with the particles' spins. Previous attempts to describe this system were limited to simple, first-order descriptions that assumed the fluid adjusted instantly to changes. However, just as a heavy truck takes time to turn or stop, the spin of particles in this super-hot fluid takes a finite amount of time to adjust to the surrounding forces. The researchers realized that to understand the longitudinal polarization puzzle, they needed to look at the second-order effects—the delayed, more complex responses of the fluid as it evolves.

Using a semi-classical approach that bridges the gap between the quantum world of individual particles and the macroscopic world of fluid dynamics, the team derived a new set of equations to track the evolution of this spin shear. They modeled the massive spin-1/2 particles, which are the fundamental building blocks of matter like electrons and quarks, moving through this fluid. By applying a method that assumes the particles eventually relax back to a state of equilibrium after being disturbed, they calculated how the spin shear component changes over time. Their work revealed that the evolution of this spin shear is driven by a complex mix of thermodynamic forces, including the expansion of the fluid, its acceleration, and the gradients of the spin polarization itself. Crucially, they found that the fluid's vorticity, or its local swirling motion, influences the spin shear only through higher-order corrections, a detail that was previously missed in simpler models.

The most significant outcome of this research is the discovery of a new relaxation-type equation that governs how the spin shear tensor evolves. This equation introduces a fresh set of transport coefficients, which are essentially numbers that describe how efficiently the fluid transmits these spin-related forces. These new coefficients suggest that the dynamics of the dissipative spin-shear tensor are more intricate than previously thought, involving cross-couplings between different non-equilibrium currents. The authors emphasize that while their current model uses a simplified approximation for particle collisions, it successfully identifies new tensor structures that were not present in earlier theories. This finding is vital because it provides a more causal and stable theoretical foundation, meaning the equations do not predict impossible speeds or unstable behaviors that plagued earlier first-order attempts.

Despite these advances, the researchers are careful to note that their work is a step forward rather than a final solution. They point out that their current formulation relies on a simplified view of particle collisions that does not fully capture the complex interactions between the spin of the particles and the orbital motion of the medium. They suggest that future work must incorporate more realistic collision mechanisms and potentially allow for different relaxation times for different parts of the system, rather than assuming they all relax at the same rate. By providing this refined second-order theory, the study offers a clearer path for numerical simulations that could finally resolve the mystery of longitudinal spin polarization in heavy-ion collisions. The hope is that by plugging these new equations into computer models with realistic initial conditions, physicists will be able to match theoretical predictions with experimental observations, turning a long-standing ambiguity into a clear understanding of how spin behaves in the most extreme environments in the universe.

Technical Summary: Second-order evolution of spin shear tensor for massive spin-1/2 particles

Problem Statement
Recent experimental observations of spin polarization in non-central heavy-ion collisions at RHIC and LHC have highlighted the necessity of understanding the dynamics of spin in relativistic fluids, specifically the Quantum Chromodynamics (QCD) medium. While global spin polarization is well-explained by thermal vorticity, the longitudinal polarization of $\Lambda$ -hyperons remains ambiguous. Previous first-order theories of spin hydrodynamics suffer from acausality and instability, similar to early relativistic hydrodynamics. Furthermore, existing frameworks often rely on assumptions of instantaneous thermalization of spin degrees of freedom, which contradicts findings that spin relaxation times are significantly longer than other diffusive modes. To address the causality issues and the complex dynamics of longitudinal polarization, a second-order theory of spin hydrodynamics is required. Specifically, recent studies suggest that the evolution of the spin shear component of the spin tensor is crucial for a successful description of experimental data.

Methodology
The authors employ a semi-classical kinetic theory framework for massive spin-1/2 particles to derive the second-order evolution equations. The methodology involves the following steps:

Framework: The study utilizes the de Groot–Leeuwen–Weert (GLW) pseudogauge to define the energy-momentum tensor and the spin tensor. The spin degrees of freedom are treated classically within an extended phase-space distribution function, $f_s(x, p, s)$ .
Conservation Laws: The authors start from the conservation of the energy-momentum tensor ( $T^{\mu\nu}$ ) and the spin tensor ( $S^{\lambda,\mu\nu}$ ). The spin tensor is decomposed into various components, including spin density, spin pressure, spin shear ( $\Sigma^{\langle\mu\nu\rangle}_{(s)}$ ), spin vorticity, and spin diffusion.
Kinetic Theory: The relativistic Boltzmann equation is solved under the Relaxation Time Approximation (RTA). The distribution function is expanded using the Chapman-Enskog method into equilibrium ( $f_{0s}$ ) and non-equilibrium ( $\delta f_s$ ) parts.
Derivation:
- First-order corrections to the distribution function are calculated to determine dissipative transport coefficients.
- Second-order evolution equations are derived by iterating the Boltzmann equation. Specifically, the authors integrate the time-derivative of the distribution function weighted by $u^\mu s_\mu \langle \nu p_\lambda \rangle$ to isolate the evolution of the spin shear component.
- The formulation is limited to local collisions and Boltzmann statistics.

Key Contributions and Results
The primary result of the paper is the derivation of a relaxation-type evolution equation for the spin shear tensor ( $\Sigma^{\langle\lambda\nu\rangle}_{(s)}$ ) up to second-order gradients in spacetime. The derived equation (Eq. 20) takes the form:
$\dot{\Sigma}^{\langle\lambda\nu\rangle}_{(s)} + \frac{\Sigma^{\langle\lambda\nu\rangle}_{(s)}}{\tau_R} = \text{Source Terms}$
Key findings include:

Driving Forces: The evolution of the spin shear is driven by thermodynamic forces including scalar expansion ( $\theta$ ), fluid acceleration ( $\dot{u}^\alpha$ ), shear stress ( $\sigma^{\alpha\beta}$ ), and the spacelike gradient of the spin polarization tensor ( $\nabla_\gamma \omega_{\alpha\beta}$ ).
Second-Order Corrections: The equation includes higher-order gradient corrections and cross-coupling terms involving the fluid vorticity tensor ( $\Omega_{\alpha\beta}$ ) and gradients of non-equilibrium currents (such as bulk pressure $\Pi$ , shear stress $\pi^{\mu\nu}$ , and spin diffusion $\delta S^{\rho,\alpha\beta}$ ).
New Transport Coefficients: The derivation yields a new set of transport coefficients (denoted as $B$ coefficients in the text) that govern the dynamics of the dissipative spin-shear tensor. These coefficients depend on thermodynamic integrals derived from the equilibrium distribution.
Coupling Mechanism: The presence of the term $(\nabla_\gamma \delta S_{\rho,\alpha\beta})$ explicitly demonstrates how the spin shear couples with other components of the spin tensor, a feature captured at the second-order level.
Relaxation Time: Within the RTA framework used, the relaxation time for the spin shear evolution is found to be identical to the relaxation time $\tau_R$ used for other hydrodynamic modes. The authors note this is an artifact of the approximation and that phenomenological applications may require distinct relaxation times.

Significance and Claims
The paper claims to provide a causal and stable second-order formulation for the evolution of the spin shear component, addressing the limitations of first-order theories. By constructing this equation using kinetic theory, the authors establish a foundation for investigating the dynamics of spin polarization without assuming immediate thermalization of spin degrees of freedom.

The authors emphasize that this work is a necessary step toward resolving the puzzle of longitudinal spin polarization in heavy-ion collisions. They note that while the current formulation is restricted to local collisions and the RTA, the identification of new tensor structures and transport coefficients suggests that the spin shear component plays a non-trivial role in spin dynamics. The paper concludes that numerical implementation of these equations with appropriate initial conditions could help clarify the mechanisms behind longitudinal polarization, though it explicitly states that further investigation is needed to understand the impact of the newly identified transport coefficients and to incorporate non-local collision kernels that account for interactions between particle spin and the medium's orbital angular momentum.

More like this