Hyperon polarization in isobaric Zr+Zr collisions at $\sqrt{s_{NN}}=200$ GeV: TRENRo3D + CLVisc with an initial longitudinal flow gradient

This paper presents a theoretical study using the TRENTo3D and CLVisc models with a novel initial longitudinal flow gradient to simultaneously describe STAR's global and azimuthal $\Lambda$ hyperon polarization measurements in isobaric Zr+Zr collisions, revealing that the azimuthal modulation $P_{y,\mathrm{c2}}$ is dominantly driven by shear effects while highlighting the challenges in achieving a unified description of all polarization observables.

The Gist

Imagine two atomic nuclei, specifically Zirconium-96, smashing into each other at nearly the speed of light. This isn't just a crash; it's a creation event. For a split second, the matter melts into a super-hot, super-dense soup called the Quark-Gluon Plasma (QGP). Scientists believe this soup behaves like a "perfect fluid," meaning it flows with almost zero friction, swirling and spinning with incredible force.

This paper is like a high-speed, 3D simulation of that crash, trying to understand how tiny particles inside this soup (called hyperons) get "spun up" or polarized, much like a spinning top.

Here is the breakdown of what the researchers did and found, using simple analogies:

1. The Setup: Building the "Perfect Storm"

To simulate this crash, the team used two main tools:

• TRENTo-3D: This is the "architect." It builds the initial shape of the crash. Imagine two soft, squishy balls (the nuclei) colliding. Usually, scientists assume the fluid flows straight out like a jet. But this team added a new twist: they allowed the fluid to have a longitudinal flow gradient.
  • Analogy: Think of a river. In the old model, the water flowed straight down the riverbed. In this new model, the water at the top of the river flows slightly faster or slower than the water at the bottom, creating a twisting motion (vorticity) right from the start.
• CLVisc: This is the "engine." It takes the shape built by TRENTo and simulates how the fluid expands, cools, and eventually freezes into particles we can detect.

2. The Mystery: Why Do Particles Spin?

When the nuclei collide off-center (like two cars grazing each other), they create a massive amount of orbital angular momentum. Think of it like a figure skater spinning with arms outstretched. The fluid created in the crash inherits this spin.

The researchers wanted to know: How does this spinning fluid make the tiny hyperon particles inside it spin?
They tested two main theories:

• The "Isothermal" Theory: Assumes the fluid freezes at a perfectly uniform temperature, like a block of ice forming evenly.
• The "Standard Thermal" Theory: Assumes the fluid has temperature gradients (hotter in the middle, cooler on the edges), like a cooling cup of coffee.

3. The Key Findings

A. The "Twist" Matters (The Longitudinal Flow)

The team discovered that the new "twist" they added to the initial flow (controlled by a parameter they call $f_v$) was essential.

• Analogy: If you try to spin a coin on a table, you need to flick it. Without that specific flick (the longitudinal flow gradient), the coin barely spins.
• Result: Without this new twist, their simulation predicted almost no polarization. With the twist set to the right amount ($f_v = 0.10$), their simulation matched the real-world data from the STAR experiment perfectly.

B. The Battle of Forces: Heat vs. Shear

The polarization of the particles comes from two competing sources:

1. Thermal Vorticity (The Spin): This comes from the fluid's rotation. It is strongest at lower speeds and gets weaker as particles move faster.
2. Shear (The Stretch): This comes from the fluid stretching and sliding past itself. It gets stronger as particles move faster.

• Result: At low speeds, the "Spin" wins. At high speeds, the "Stretch" takes over. The combination of these two forces explains why the polarization behaves the way it does across different speeds.

C. The Shape of the Nucleus Doesn't Matter Much

The researchers tested if the specific "shape" of the Zirconium nucleus (is it slightly squashed? does it have a weird bump?) changed the results.

• Analogy: Imagine trying to tell if a spinning top is made of wood or plastic just by watching how fast it spins.
• Result: It didn't matter. Whether they used the "standard" Zirconium shape or alternative shapes from the blind analysis, the polarization results were almost identical. The spin is driven more by the overall crash energy and flow than by the tiny details of the nuclear shape.

D. The "Side-to-Side" vs. "Up-and-Down" Spin

The team looked at two types of polarization:

• Out-of-Plane ($P_y$): Spinning like a wheel rolling on the ground.
  • Result: The "Isothermal" model (uniform temperature) worked great here. It matched the data perfectly.
• Longitudinal ($P_z$): Spinning like a top standing up.
  • Result: This was tricky. The "Isothermal" model got the direction of the spin right (matching the real data), but it predicted the spin was too strong at high speeds. The "Standard Thermal" model (with temperature gradients) got the direction wrong (it predicted the opposite spin direction).
  • Conclusion: Neither model is perfect yet. The "Isothermal" model is better for the direction, but both struggle to explain why the spin isn't as strong as predicted at very high speeds.

4. What This Means

This paper is a major step forward because it successfully simulates a complex 3D collision and matches real experimental data for the first time in this specific setup.

• The Good News: They found that adding a specific "longitudinal flow" to the simulation is crucial to explain why particles spin. They also proved that the "Isothermal" (uniform temperature) approach is the better way to calculate the direction of the spin.
• The Open Question: They still can't fully explain why the spin is weaker than predicted at very high speeds. It suggests there are other physical forces (like bulk viscosity or electromagnetic fields) acting as a "brake" that their current model doesn't fully capture yet.

In short, the researchers built a better 3D map of the atomic crash, found the missing "twist" that makes the particles spin, and identified exactly where their current understanding of the physics needs a little more work.

Technical Summary

Technical Summary: Hyperon Polarization in Isobaric Zr+Zr Collisions

Problem Statement
The paper addresses the theoretical description of global and azimuthal-angle-dependent $\Lambda$ hyperon polarization in isobaric $^{96}_{40}\text{Zr}+^{96}_{40}\text{Zr}$ collisions at $\sqrt{s_{NN}} = 200$ GeV. While the global polarization ($-P_y$) in non-central heavy-ion collisions is established as a signature of the vortical nature of the Quark-Gluon Plasma (QGP), a comprehensive, self-consistent theoretical framework capable of simultaneously describing recent STAR measurements of global polarization, its azimuthal modulation coefficients ($P_{y,c0}$ and $P_{y,c2}$), and longitudinal polarization ($P_z$) has been lacking. Specific challenges include the treatment of temperature-gradient terms (the "sign puzzle" regarding $P_z$), the constraints on the three-dimensional initial fireball geometry, and the sensitivity of polarization to isobaric nuclear structure variations.

Methodology
The authors employ a hybrid theoretical framework combining:

1. TRENTo-3D Initial Conditions: A 3D initial condition model is used to generate the energy density profile. A novel extension is introduced: a tunable parameter $f_v$ controls an initial longitudinal flow velocity gradient ($v_{\eta_s}$), moving beyond the standard Bjorken flow approximation. This gradient is derived from the local center-of-mass rapidity $\eta_{cm}$ to provide a source of initial vorticity in the symmetric isobaric system. The model also incorporates a deformed Woods-Saxon nuclear distribution for $^{96}\text{Zr}$, including octupole deformation ($\beta_3$), and fragmentation regions controlled by a transverse momentum scale $k_T$.
2. CLVisc Hydrodynamics: The initial state evolves via a (3+1)-D viscous hydrodynamic model (CLVisc) with a specific shear viscosity $\eta/s = 0.08$. The study explores various bulk viscosity ($\zeta/s$) scenarios, including constant values and a temperature-dependent parametrization peaking near the QCD crossover.
3. Isothermal Polarization Formalism: The spin polarization is calculated using the isothermal equilibrium scenario, where the freeze-out hypersurface is approximated as isothermal ($T_H \approx 155$ MeV). In this limit, temperature gradients vanish, and the polarization decomposes into contributions from kinematic vorticity ($\omega_{\alpha\beta}$) and kinematic shear ($\Xi_{\alpha\beta}$). The authors also perform a comparative analysis using the "standard thermal" scenario, which retains full temperature gradients.

Key Contributions

• Introduction of Longitudinal Flow Gradient: The paper introduces the parameter $f_v$ into the TRENTo-3D initialization for the first time, providing a mechanism to generate initial vorticity in symmetric isobaric collisions where participant thickness imbalances average to zero.
• Simultaneous Description of Observables: The framework achieves a simultaneous description of STAR data for global polarization ($-P_y$) across centrality, transverse momentum ($p_T$), and pseudorapidity ($\eta$), as well as the azimuthal modulation coefficients $P_{y,c0}$ and $P_{y,c2}$.
• Decomposition of Polarization Mechanisms: The study explicitly separates the contributions of thermal vorticity and shear stress. It demonstrates that $P_{y,c0}$ is dominated by thermal vorticity, while $P_{y,c2}$ is primarily shear-driven, offering a clean probe for shear-induced polarization.
• Systematic Parameter Scans: The authors systematically scan $f_v$, $k_T$, nuclear structure configurations (comparing $^{96}\text{Zr}$ with alternative $^{96}\text{Ru}$ parameter sets), and bulk viscosity to constrain the initial state and assess model sensitivity.

Results

• Global Polarization ($-P_y$): The model reproduces the centrality dependence and the rising trend of $-P_y$ with $p_T$. The $p_T$ dependence is shown to result from a competition: the thermal vorticity term decreases with $p_T$, while the shear term increases, eventually dominating at higher $p_T$.
• Parameter Sensitivity:
  • $f_v$ primarily controls the overall magnitude of $-P_y$ and the relative weight of the thermal vorticity contribution. A value of $f_v = 0.10$ is favored to match data.
  • $k_T$ affects the low-$p_T$ level and the $\eta$ profile; increasing $k_T$ enhances the shear contribution and creates a convex $\eta$ profile.
  • Nuclear Structure: Calculations using five different nuclear structure configurations (including variations in $\beta_2$ and $\beta_3$) yield nearly indistinguishable polarization results. The polarization is insensitive to nuclear deformations at the $\beta \sim 0.1$ level, suggesting that global polarization is not a sensitive probe for distinguishing isobaric nuclear structures within current experimental precision.
• Azimuthal Modulation: The model successfully describes the Fourier coefficients $P_{y,c0}$ and $P_{y,c2}$, confirming the distinct physical origins of the constant and quadrupole terms.
• Longitudinal Polarization ($P_z$): The isothermal scenario correctly reproduces the sign and amplitude of the azimuthal modulation of $P_z$ (negative sinusoidal pattern). However, the model significantly overpredicts the $p_T$-dependent modulation amplitude $\langle P_z \sin[2(\phi-\Psi_2)] \rangle$ at high $p_T$ ($>1.2$ GeV).
• Isothermal vs. Standard Thermal: A comparison reveals that while the standard thermal scenario offers a modest improvement in describing the $p_T$ dependence of $-P_y$, it fails qualitatively for $P_z$, predicting the wrong sign for the azimuthal modulation. The isothermal scenario correctly captures the $P_z$ phase but struggles with the high-$p_T$ magnitude.

Significance and Claims
The paper claims that the TRENTo-3D + CLVisc framework with the isothermal polarization scenario provides a solid baseline for understanding hyperon polarization in isobaric collisions. The introduction of the longitudinal flow gradient $f_v$ is identified as essential for generating vorticity in symmetric systems. The work highlights that while the isothermal approximation captures the essential physics of the azimuthal phase of $P_z$ and the global polarization magnitude, it does not provide a unified description of all observables, particularly the high-$p_T$ longitudinal polarization.

The authors modestly conclude that the persistent discrepancy in high-$p_T$ $P_z$ suggests the need for additional physical mechanisms beyond the current model, such as bulk viscosity effects (which were found insufficient to resolve the tension), electromagnetic fields, or late-stage hadronic rescattering. The study establishes that global polarization is a powerful probe of the initial flow velocity field and fireball geometry but is not sensitive to nuclear deformations at the level probed by the Zr/Ru isobaric system. The decomposition into Fourier coefficients is proposed as a valuable tool for isolating shear-induced polarization in future measurements.