An improved formula for Wigner function and spin polarization in a decoupling relativistic fluid at local thermodynamic equilibrium

This paper presents an improved formula for the Wigner function and spin polarization of fermions in a decoupling relativistic fluid by introducing a new expansion method that captures geometric dependencies in the spin-shear term, naturally justifies the isothermal condition, and extends the framework to particles with arbitrary spin.

The Gist

Imagine a massive, super-hot explosion, like the kind that happens when scientists smash heavy atoms together to recreate the conditions of the early universe. This explosion creates a tiny drop of "quark-gluon plasma" (QGP)—a soup of particles so hot and dense that they behave like a fluid. As this fluid expands and cools, it eventually reaches a point where the particles stop bumping into each other and fly off into space. Scientists call this moment "decoupling."

The paper you're asking about is like a new, upgraded instruction manual for predicting how these particles spin as they fly away.

Here is the breakdown of what the authors did, using simple analogies:

1. The Problem: The Old Map Was Too Rigid

Previously, scientists had a formula to predict the spin polarization (which way the particles are spinning) of these flying particles. However, that old formula relied on a very specific, idealized shape for the "edge" where the fluid stops.

• The Analogy: Imagine trying to predict how water splashes off a wall. The old formula only worked if the wall was perfectly flat and vertical. But in reality, the edge of this plasma fluid is wavy, curved, and irregular, like a crumpled piece of paper or the surface of a rolling wave. The old formula tried to force this complex shape into a flat box, which led to inaccuracies.

2. The Solution: A New Way to Look at the Edge

The authors developed a new mathematical method that works no matter what shape the edge of the fluid takes.

• The Analogy: Instead of forcing the fluid into a flat box, they invented a new way of "scanning" the edge. Imagine taking a photo of a curved, bumpy surface. The old method tried to flatten the photo before analyzing it. The new method analyzes the photo exactly as it is, respecting every curve and bump.
• The "Worldline" Trick: A key part of their new method involves looking at the path a particle takes (its "worldline"). They realized that to know how a particle spins at a specific point, you don't just look at that exact spot; you have to look at where that particle's path intersects the fluid's edge. Sometimes, a particle's path might cross the edge, go back in, and cross it again (like a boomerang path). Their formula accounts for all these crossing points, not just the first one.

3. The Big Discovery: Why "Isothermal" Matters

One of the most interesting findings is about temperature gradients (changes in temperature).

• The Old Confusion: In previous calculations, scientists had to manually assume that the temperature was the same everywhere along the edge of the fluid (an "isothermal" condition) to make the math work. It was like saying, "We'll just pretend the edge is all the same temperature because the math is too hard otherwise."
• The New Insight: The authors' new formula naturally shows that if the edge is indeed at a constant temperature, the messy temperature differences cancel themselves out automatically. You don't have to force the assumption; the math proves it happens on its own. It's like discovering that a complicated machine naturally balances itself without you needing to add a counterweight.

4. What They Found (The "Spin" Results)

Using this new, flexible formula, they updated the recipe for calculating spin. They found three main ingredients that determine how the particles spin:

1. Thermal Vorticity: Think of this as the "swirl" or "whirlpool" effect in the fluid. If the fluid is spinning like a tornado, the particles will spin with it.
2. Thermal Shear: This is like stretching or squeezing the fluid. If you pull the fluid in one direction and push it in another, it creates a different kind of spin. The new formula fixes how this stretching affects the spin, correcting errors from the old "flat wall" assumption.
3. Spin Hall Effect: This is a subtle quantum effect where particles drift sideways based on their spin, similar to how a car might drift on a wet road.

5. The "Ghost" Particles

The new math revealed some strange extra terms that seemed to suggest particles were coming from places they shouldn't be (like particles moving into the fluid from the outside). The authors suggest these are likely "ghosts" or mathematical artifacts caused by the way they modeled the fluid. They propose a simple fix: just ignore any paths where the particle is moving into the fluid, keeping only the ones flying out. This aligns with how other physicists have handled similar problems in the past.

Summary

In short, this paper provides a better, more flexible ruler for measuring how particles spin as they escape a hot, expanding fluid. It removes the need for unrealistic geometric assumptions, proves that temperature effects cancel out naturally under certain conditions, and offers a more accurate way to understand the quantum "spin" of matter in the most extreme environments in the universe.

Technical Summary

Problem Statement
Relativistic heavy-ion collisions create a quark-gluon plasma (QGP) that evolves through local thermodynamic equilibrium (LTE) before decoupling into hadrons. Calculating physical observables, particularly spin polarization, at the decoupling hypersurface ($\Sigma_D$) requires a rigorous framework based on quantum statistical field theory. Previous calculations of the Wigner function and spin polarization at LTE, specifically those in references [1] and [2], relied on significant geometric approximations. These included assuming the decoupling hypersurface was a hyperplane with infinitely distant time-like boundaries or that it was strictly perpendicular to the fluid four-velocity field. However, analytical models and numerical simulations indicate that the actual decoupling hypersurface in heavy-ion collisions possesses arbitrary geometry and does not satisfy these restrictive conditions. Furthermore, previous derivations did not naturally exclude contributions from space-time gradients in the normal direction of the hypersurface, requiring the isothermal condition ($T=$ const) to be imposed a priori to eliminate temperature gradient contributions.

Methodology
The authors propose a new expansion method to evaluate observables at LTE on decoupling hypersurfaces of arbitrary geometry. The core of this approach involves inverting the order of integration in the linear response theory calculation of the Wigner function.

1. Linear Response Framework: Starting from the Zubarev density operator, the mean value of a local operator is expanded around global equilibrium (GE) using the difference in the four-temperature field, $\Delta\beta_\nu(y, x) = \beta_\nu(y) - \beta_\nu(x)$.
2. Integration Order Exchange: Unlike previous works that integrated over the hypersurface $\Sigma_D$ first, this method performs the momentum integration ($q$) first. This is facilitated by the hydrodynamic limit, where $\Delta\beta_\nu$ is a slowly varying function, allowing for a small-$q$ expansion of the integrand.
3. Geometric Handling: The method utilizes a decomposition of the hypersurface into single-valued branches and employs the Leibniz rule to handle derivatives. A key operator, $D_y$, is introduced to act on the intersection points $\bar{y}(x, p)$ between the particle world-line and the hypersurface. This formalism naturally incorporates gradients of the normal vector $n^\mu$ to the hypersurface.
4. Spin Polarization Calculation: The derived Wigner function is applied to calculate the spin polarization of spin-1/2 fermions. The authors also extend the formalism to include finite chemical potential by replacing the stress-energy tensor with the conserved current operator.

Key Contributions and Results
The paper derives an upgraded, compact formula for the Wigner function and the resulting spin polarization vector $S^\mu(p)$.

• Generalized Geometry: The new formula is applicable to hypersurfaces with arbitrary geometry, removing the need for the rectangular or perpendicular-velocity approximations used in prior literature.
• Shear-Induced Polarization: The expression for the shear-induced polarization is upgraded. The geometric factor previously dependent on specific assumptions (e.g., $u^\mu/p \cdot u$ or $\hat{t}^\mu/p \cdot \hat{t}$) is replaced by the general term $n^\mu / |p \cdot n|$, where $n^\mu$ is the unit normal vector to the hypersurface. This resolves ambiguities in the shear-induced polarization formula found in previous studies.
• Natural Exclusion of Normal Gradients: A significant theoretical result is that the new method naturally excludes contributions from space-time gradients in the normal direction of the hypersurface. Consequently, if the decoupling hypersurface is isothermal, the temperature gradient contributions to spin polarization vanish naturally without being imposed as an external condition.
• Thermal Vorticity and Shear: The leading-order spin polarization includes contributions from thermal vorticity ($\varpi_{\mu\nu}$), thermal shear ($\xi_{\mu\nu}$), and the spin Hall effect. The thermal vorticity term acquires a modification via the sign function $\text{sgn}(p \cdot n)$, which is necessary when the hypersurface is not entirely space-like and future-oriented.
• Spurious Intersections: The derivation reveals contributions from non-trivial intersections where the particle world-line crosses the hypersurface at points $\bar{y} \neq x$. The authors suggest these may be spurious effects arising from the use of free fields in the Wigner operator for the in-plasma region and propose discarding them using a cutoff $\theta(p \cdot n)$, analogous to the Cooper-Frye formula.
• Curvature Independence: The authors demonstrate that, at the leading order, the curvature of the decoupling hypersurface does not contribute to the spin polarization.

Significance
The paper claims to provide a theoretically robust framework for calculating spin polarization in relativistic heavy-ion collisions that is consistent with the arbitrary geometry of the decoupling hypersurface observed in numerical simulations. By deriving the isothermal condition naturally from the formalism rather than imposing it, the work offers a deeper theoretical justification for previous assumptions. The upgraded formula for shear-induced polarization resolves long-standing ambiguities in the literature. The authors state that this framework can be extended to particles with arbitrary spin and merits verification in numerical simulations of relativistic nuclear collisions.