Boost-invariant perfect Fermi-Dirac spin hydrodynamics

This paper demonstrates the feasibility of incorporating Fermi-Dirac statistics into boost-invariant perfect spin hydrodynamics simulations for spin-1/2 particles, revealing that while the resulting corrections are an order of magnitude smaller than spin feedback effects, the approach requires careful handling as numerical solutions can break down under extreme spin polarization.

The Gist

The Big Picture: A Spinning Soup

Imagine a giant, super-hot soup made of tiny particles (like protons and neutrons) created when heavy atoms smash into each other at nearly the speed of light. This soup is so hot and dense that it behaves like a fluid.

Now, imagine that every single particle in this soup is also a tiny spinning top. They aren't just moving around; they are all spinning in specific directions. This is what physicists call spin.

For a long time, scientists have been trying to write the "recipe" (the math) to describe how this spinning soup flows and evolves. This paper is about refining that recipe.

The Problem: The "Lazy" Approximation

In the past, when scientists tried to calculate how this soup behaves, they used a "lazy" shortcut. They assumed the particles were like a very thin, sparse gas where they rarely bump into each other. In physics, this is called the Boltzmann approximation.

Think of it like this: If you are trying to predict how a crowd of people moves in a huge, empty stadium, you can assume everyone walks in a straight line without bumping into anyone. That's the Boltzmann approach.

But in a heavy-ion collision, the "stadium" is packed shoulder-to-shoulder. The particles are crowded, and they interact constantly. In this crowded environment, the "lazy" shortcut isn't accurate enough. The particles follow Fermi–Dirac statistics, which is the rulebook for how crowded quantum particles behave (like how no two people can sit in the exact same seat).

The Goal of the Paper:
The authors wanted to swap the "lazy" shortcut for the "crowded" rulebook to see if it changes the results of their simulations. They asked: Does it matter if we treat these spinning particles as a dense crowd or a sparse gas?

The Experiment: Two Ways to Spin

To test this, the authors set up a simulation of the soup expanding. They looked at two specific ways the particles could be spinning:

1. The Longitudinal Spin: Imagine all the spinning tops are spinning along the direction the soup is flowing (like a train moving forward).
2. The Transverse Spin: Imagine the spinning tops are spinning sideways, perpendicular to the flow (like a wheel spinning as the train moves).

They ran the simulation twice for each scenario: once using the old "lazy" math (Boltzmann) and once using the new "crowded" math (Fermi–Dirac).

The Findings

1. The Differences are Small, but Real

When they compared the two results, they found that the "crowded" math didn't change the story of what happened. The soup still expanded and cooled down the same way.

However, the numbers were slightly different. The difference was about 8.5% in the worst-case scenarios.

• Analogy: Imagine you are baking a cake. The "lazy" recipe says you need 1 cup of sugar. The "crowded" recipe says you need 1.08 cups. The cake will still taste like a cake, but if you are a perfectionist, that extra sugar matters.

2. The "Exploding" Soup (The Singularity)

Here is the most dramatic part of the paper. When they used the "Longitudinal Spin" setup (spinning along the flow) with very strong initial spins, the simulation crashed.

In the math world, this is called a "singularity" or a "blow-up." It's like trying to divide by zero. The numbers got so big, so fast, that the computer couldn't handle them, and the simulation stopped working.

• Why did this happen? The authors realized it wasn't because the physics was "wrong" or because the particles were too crowded. It was a flaw in the specific mathematical recipe they were using for that specific direction of spin. It's like a bridge that holds up fine for cars, but if you drive a truck in a specific lane, the bridge collapses because of a design flaw, not because the truck is too heavy.

Interestingly, when they tried the "Transverse Spin" (spinning sideways), the simulation never crashed, even when they made the spins incredibly strong. The sideways spinning was much more stable.

The Solution: A New Tool

To make the "crowded" math work, the authors had to invent a new tool. The equations for the crowded particles involve some very complex, weird-looking integrals (mathematical sums) that don't exist in standard calculator software.

• The Analogy: Imagine you need to calculate the weight of a cloud, but you don't have a scale. You have to build a custom scale out of wood and string.
• What they did: They created a "lookup table" (a pre-calculated map) for these complex math functions. This allows computers to quickly find the answers without doing the heavy lifting every time. They proved this method works and is fast enough for real-world simulations.

The Takeaway

1. It matters, but not a lot: Using the correct "crowded" math (Fermi–Dirac) instead of the "lazy" math (Boltzmann) gives slightly more accurate results (up to ~8.5% difference), which is important for precision physics but doesn't change the big picture.
2. Math can break: Even if a theory is stable and makes sense physically, the specific equations used to simulate it can sometimes "break" (crash) if you push the conditions too hard, especially in certain directions.
3. We are ready for the next step: The authors have built the necessary tools (the lookup tables) to run these more accurate simulations. This brings us one step closer to perfectly understanding the "spinning soup" created in particle colliders, which helps us understand the very early universe.

In short: They upgraded the software for simulating spinning particle soup. The new version is more accurate, and they discovered that while the soup is stable in some directions, it can cause the math to explode in others—a crucial clue for future scientists.

Technical Summary

1. Problem Statement

Relativistic spin hydrodynamics is a theoretical framework developed to describe spin-polarization phenomena (such as $\Lambda$ hyperon polarization) observed in heavy-ion collisions. Previous numerical simulations of spin-1/2 particles in this framework relied on the Boltzmann approximation for particle statistics. This approximation assumes a dilute system where quantum statistical effects (Fermi–Dirac or Bose–Einstein) are negligible.

The authors address two primary gaps:

1. Statistical Accuracy: The validity of the Boltzmann approximation in realistic heavy-ion collision scenarios, particularly at moderate energies where the baryon chemical potential ($\mu$) is high and the system is not extremely dilute.
2. Numerical Stability: Previous works vaguely noted that solutions for perfect spin hydrodynamics could break down for large spin polarization values. The specific conditions and mechanisms for this breakdown in different geometric configurations were not rigorously investigated.

2. Methodology

The authors performed a numerical study of perfect spin hydrodynamics for spin-1/2 particles obeying Fermi–Dirac statistics, comparing the results directly with the Boltzmann approximation.

• Theoretical Framework:

  • The system is modeled using the hybrid approach (combining kinetic theory for spin-1/2 particles with classical spin and the Israel–Stewart entropy-based method for dissipative effects, though this study focuses on the perfect fluid limit).
  • The hydrodynamic tensors (baryon current $N^\mu$, energy-momentum tensor $T^{\mu\nu}$, and spin tensor $S^{\lambda\mu\nu}$) are expanded up to second order in the dimensionless spin polarization tensor $\omega_{\mu\nu}$.
  • The spin tensor is treated at first order in $\omega$.
  • The study utilizes the Bjorken expansion geometry (boost-invariant, transversely homogeneous).

• Geometric Configurations:
  To make the system of equations exactly determined (reducing the overdetermined system of 11 equations for 8 unknowns), two specific symmetry configurations were analyzed:

  1. Longitudinal Configuration: Spin vectors $C_k$ and $C_\omega$ align with the beam axis ($z$-axis).
  2. Transverse Configuration: Spin vectors lie in the transverse ($x-y$) plane and are parallel to each other.

• Implementation of Fermi–Dirac Statistics:

  • Unlike the Boltzmann case, which yields coefficients expressed via modified Bessel functions ($K_n$), the Fermi–Dirac case introduces special functions defined by integrals ($J_{mnk}$).
  • Since these functions are not standard in numerical packages, the authors constructed interpolating functions based on tabulated values of the logarithms of these integrals over a wide parameter range ($\xi = \mu/T$ from -100 to 100, and $z = m/T$ from 0.1 to 100).

• Simulation Parameters:

  • Initial proper time $\tau_0 = 1$ fm, final $\tau_f = 10$ fm.
  • Particle mass: $\Lambda$ hyperon ($1.116$ GeV).
  • Initial conditions: $T_0 = 155$ MeV, $\mu_0 = 800$ MeV (typical for moderate-energy collisions).
  • Initial spin polarization coefficients ($C_{k}, C_{\omega}$) varied from 0.25 to 0.75 (and higher for stability tests).

3. Key Contributions

• First Numerical Implementation of FD Spin Hydro: This is the first study to numerically simulate perfect spin hydrodynamics using exact Fermi–Dirac statistics rather than the Boltzmann approximation.
• Parametrization of Special Functions: The authors successfully developed a method to handle the non-standard special functions required for the Fermi–Dirac case, demonstrating that the approach is computationally feasible.
• Singularity Analysis: A rigorous investigation into the "breakdown" of numerical solutions was conducted, identifying the specific mathematical mechanism causing finite-time blow-ups in the longitudinal configuration.
• Configuration Comparison: The study revealed a stark difference in stability between longitudinal and transverse spin configurations.

4. Results

A. Comparison: Fermi–Dirac vs. Boltzmann

• Evolution Trends: The general evolution of the spin polarization coefficients ($C_k$ and $C_\omega$) remains qualitatively similar between the two statistics (monotonic increase for $C_k$, slow growth or decrease for $C_\omega$).
• Quantitative Differences: The relative differences between the Fermi–Dirac and Boltzmann results are small but non-negligible.
  • Differences reach approximately 8.5% in specific cases (e.g., transverse configuration with high initial coefficients).
  • The authors conclude that while the Boltzmann approximation is a reasonable first step, the Fermi–Dirac corrections are significant enough to warrant inclusion in precision studies, especially at high baryon densities.

B. Numerical Stability and Singularities

• Longitudinal Configuration (Unstable):
  • Solutions develop finite-time singularities (blow-ups) when the initial magnitude of the spin polarization coefficients exceeds a critical threshold (approx. $C_0 \approx 0.87$ for the tested parameters).
  • Mechanism: The breakdown is caused by the denominator of the explicit ordinary differential equation for temperature ($\dot{T}$) vanishing and changing sign. This occurs when the algebraic constraint derived from baryon conservation and the energy equation become incompatible for the given initial conditions.
  • This failure happens regardless of whether Boltzmann or Fermi–Dirac statistics are used, indicating it is a structural issue of the second-order expansion in this geometry.
• Transverse Configuration (Stable):
  • Solutions remained stable and converged up to $\tau_f = 10$ fm for all tested initial values, including extremely large coefficients (up to $10^6$).
  • Mechanism: The transverse equations involve damping terms (due to the specific dependence on $\tau$ and the coefficients $A$ and $A_1$) that prevent the denominator from vanishing, unlike the longitudinal case.

C. Applicability Criteria

The study clarifies that the standard convergence criterion for the generating function (which limits the magnitude of $\omega$) is irrelevant for theories expanded to a finite order in $\omega$. In a finite-order expansion, the exponential term becomes a polynomial, ensuring integrals always converge. Therefore, the observed singularities are not due to the breakdown of the statistical definition but are artifacts of the specific hydrodynamic equations and geometry.

5. Significance and Outlook

• Feasibility: The paper proves that incorporating exact Fermi–Dirac statistics into spin hydrodynamics simulations is technically feasible and computationally manageable via interpolation.
• Precision: The ~8.5% deviation highlights that the Boltzmann approximation may introduce systematic errors in quantitative comparisons with experimental data, particularly in high-baryon-density regimes.
• Theoretical Limits: The identification of singularities in the longitudinal configuration serves as a warning that "causality and stability" (properties of the differential equations) do not guarantee "global existence" of solutions for all initial data.
• Future Directions: The authors suggest future work should include:
  • Higher-order expansions in $\omega$.
  • Inclusion of dissipative effects.
  • Simulations in more complex geometries (e.g., 3D cylindrically symmetric).

In summary, this work refines the theoretical foundation of spin hydrodynamics by moving beyond the dilute gas approximation and provides critical insights into the numerical stability limits of the theory under different geometric constraints.