Recognition: unknown
Boost-invariant perfect Fermi-Dirac spin hydrodynamics
Pith reviewed 2026-05-10 06:00 UTC · model grok-4.3
The pith
Fermi-Dirac statistics produce hydrodynamic evolutions that differ from Boltzmann results by amounts smaller than spin feedback corrections in boost-invariant spin hydrodynamics.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In boost-invariant transversely homogeneous perfect spin hydrodynamics with second-order corrections to the baryon current and energy-momentum tensor in the spin polarization tensor ω and first-order spin tensor, the evolution of parameters under Fermi-Dirac statistics differs from the Boltzmann case by amounts about one order of magnitude smaller than the corrections arising from spin feedback, and the numerical approach is feasible through parametrization of the special functions involved.
What carries the argument
The ordering of corrections in powers of the spin polarization tensor ω, combined with parametrization of Fermi-Dirac integral functions in the hydrodynamic coefficients, under boost-invariant transversely homogeneous flow.
If this is right
- The numerical solutions can be obtained reliably for moderate values of spin polarization.
- Breakdown of solutions occurs at very large spin polarization in one of the geometric configurations.
- Spin feedback effects dominate the evolution over the choice of particle statistics.
- Previous works using Boltzmann approximation capture the main features accurately for the considered initial conditions.
Where Pith is reading between the lines
- This supports using Boltzmann statistics as a default in spin hydrodynamics simulations unless high precision on quantum effects is needed.
- Extensions to non-boost-invariant or inhomogeneous geometries could reveal larger statistical differences.
- Testing the breakdown condition experimentally in heavy-ion collisions might constrain the range of applicability of perfect spin hydrodynamics.
Load-bearing premise
The corrections to the baryon current and energy-momentum tensor can be treated as second order in the spin polarization while keeping the spin tensor first order, and the boost-invariant transversely homogeneous geometry holds during the entire evolution.
What would settle it
A direct numerical comparison showing that the evolution differences between Fermi-Dirac and Boltzmann statistics exceed the magnitude of spin feedback corrections for the initial conditions used in prior studies would falsify the main result.
Figures
read the original abstract
We analyze the effect of using the Fermi-Dirac statistics, rather than its Boltzmann approximation, in numerical simulations of perfect spin hydrodynamics of particles with spin 1/2. The system considered is boost invariant, transversely homogeneous, with corrections to the baryon current and the energy-momentum tensor that are second order in the spin polarization tensor $\omega$, and the spin tensor considered is first order in $\omega$. The study shows the feasibility of this approach, as the special functions defined by integrals that appear in the coefficients in the Fermi-Dirac case can be conveniently parametrized. For sets of initial conditions used in previous works, the differences in parameter evolution between the two underlying particle statistics are about one order of magnitude smaller than corrections coming from spin feedback. We also discuss when and why the numerical solutions of the equations of perfect spin hydrodynamics break down for very large values of spin polarization in one of the geometric configurations considered.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes the effect of Fermi-Dirac statistics versus the Boltzmann approximation in numerical simulations of perfect spin hydrodynamics for spin-1/2 particles. The setup is boost-invariant and transversely homogeneous, with baryon current and energy-momentum tensor corrections kept second-order in the spin polarization tensor ω while the spin tensor is first-order in ω. It demonstrates that the special functions in the Fermi-Dirac coefficients can be parametrized for numerical use, and reports that for initial conditions from prior works, the differences in parameter evolution between the statistics are about one order of magnitude smaller than spin-feedback corrections. The paper also discusses the conditions for numerical breakdown at very large spin polarization.
Significance. If the central claim holds, the result supports the practical sufficiency of the Boltzmann approximation in this class of spin-hydrodynamic simulations, since statistical differences are subdominant to spin effects. This could simplify modeling in applications such as heavy-ion collisions. The explicit parametrization of the Fermi-Dirac integrals is a concrete technical contribution that enables the reported numerics.
major comments (2)
- The central claim that Fermi-Dirac vs. Boltzmann differences are ~10× smaller than spin corrections for previous-work initial conditions rests on the maintained validity of the boost-invariant, transversely homogeneous ansatz and the strict second-order truncation in ω throughout the evolution. The manuscript notes numerical breakdown for large ω but does not provide explicit checks (e.g., time evolution of |ω| or higher-order residuals) confirming that the chosen initial conditions keep the system inside the regime where the ordering and geometry remain self-consistent; without this, the reported ordering of effects cannot be guaranteed inside the stated framework.
- The feasibility statement for the Fermi-Dirac case relies on parametrization of the special functions arising from the integrals. No quantitative error bounds, convergence tests against direct integration, or sensitivity analysis of the reported parameter differences to the parametrization accuracy are supplied; if the parametrization error is comparable to the ~10× smaller statistical difference, the comparison between statistics would be compromised.
minor comments (1)
- The abstract and introduction would benefit from a brief statement of the typical magnitude of ω (or range of initial conditions) for which the ordering is expected to hold, to make the domain of applicability clearer.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below and agree that the suggested additions will strengthen the presentation of our results.
read point-by-point responses
-
Referee: The central claim that Fermi-Dirac vs. Boltzmann differences are ~10× smaller than spin corrections for previous-work initial conditions rests on the maintained validity of the boost-invariant, transversely homogeneous ansatz and the strict second-order truncation in ω throughout the evolution. The manuscript notes numerical breakdown for large ω but does not provide explicit checks (e.g., time evolution of |ω| or higher-order residuals) confirming that the chosen initial conditions keep the system inside the regime where the ordering and geometry remain self-consistent; without this, the reported ordering of effects cannot be guaranteed inside the stated framework.
Authors: We agree that explicit verification of the regime of validity for the selected initial conditions is important to confirm self-consistency. In the revised manuscript we will add plots of the time evolution of |ω| for all initial conditions considered, together with estimates of the magnitude of higher-order terms in ω. These additions will demonstrate that the boost-invariant, transversely homogeneous ansatz and the second-order truncation remain valid throughout the evolution, thereby supporting the reported ordering of statistical versus spin-feedback effects. revision: yes
-
Referee: The feasibility statement for the Fermi-Dirac case relies on parametrization of the special functions arising from the integrals. No quantitative error bounds, convergence tests against direct integration, or sensitivity analysis of the reported parameter differences to the parametrization accuracy are supplied; if the parametrization error is comparable to the ~10× smaller statistical difference, the comparison between statistics would be compromised.
Authors: We acknowledge the need for quantitative validation of the parametrization. In the revised version we will include direct comparisons of the parametrized functions against numerical integration of the defining integrals, together with explicit error bounds and a sensitivity study showing the effect of parametrization inaccuracies on the reported differences between Fermi-Dirac and Boltzmann evolution. This will confirm that parametrization errors lie well below the scale of the statistical differences. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper performs numerical simulations of boost-invariant transversely homogeneous perfect spin hydrodynamics, parametrizing special functions arising in the Fermi-Dirac case to demonstrate feasibility. It then compares parameter evolution under Fermi-Dirac versus Boltzmann statistics for initial conditions drawn from prior literature, reporting that the observed differences are smaller than spin-feedback corrections. This comparison is an empirical outcome of integrating the stated hydrodynamic equations under fixed ordering (baryon current and EMT corrections second-order in ω, spin tensor first-order) and geometry; no step in the reported chain reduces by construction to a fitted parameter renamed as a prediction, a self-definitional relation, or a load-bearing self-citation that would force the result. The framework assumptions are explicitly declared and the breakdown for large ω is noted as a limitation, keeping the central claim self-contained as a numerical finding rather than a tautology.
Axiom & Free-Parameter Ledger
free parameters (1)
- parametrization coefficients for Fermi-Dirac integrals
axioms (2)
- domain assumption The system remains boost-invariant and transversely homogeneous throughout the evolution
- domain assumption Corrections to baryon current and energy-momentum tensor are second order in ω while spin tensor is first order
Forward citations
Cited by 1 Pith paper
-
Boost-invariant and cylindrically symmetric perfect spin hydrodynamics
In boost-invariant cylindrical spin hydrodynamics, azimuthal-longitudinal coupling in the spin tensor produces nonzero total polarization only via the longitudinal magnetic component coupled to the azimuthal electric ...
Reference graph
Works this paper leans on
-
[1]
Polarization phenomenon in heavy-ion collisions,
T. Niida and S. A. Voloshin, Int. J. Mod. Phys. E33, 2430010 (2024), arXiv:2404.11042 [nucl-ex]
-
[2]
Global $\Lambda$ hyperon polarization in nuclear collisions: evidence for the most vortical fluid
L. Adamczyket al.(STAR), Nature548, 62 (2017), arXiv:1701.06657 [nucl-ex]
work page Pith review arXiv 2017
-
[3]
J. Adamet al.(STAR), Phys. Rev. C98, 014910 (2018), arXiv:1805.04400 [nucl-ex]
-
[4]
J. Adamet al.(STAR), Phys. Rev. Lett.123, 132301 (2019), arXiv:1905.11917 [nucl-ex]
-
[5]
S. Acharyaet al.(ALICE), Phys. Rev. Lett.125, 012301 (2020), arXiv:1910.14408 [nucl-ex]
-
[6]
The ideal relativistic rotating gas as a perfect fluid with spin
F. Becattini and L. Tinti, Annals Phys.325, 1566 (2010), arXiv:0911.0864 [gr-qc]
work page Pith review arXiv 2010
-
[7]
Local Polarization and Isothermal Local Equilibrium in Relativistic Heavy Ion Collisions,
F. Becattini, M. Buzzegoli, G. Inghirami, I. Karpenko, and A. Palermo, Phys. Rev. Lett.127, 272302 (2021), arXiv:2103.14621 [nucl-th]
-
[8]
A. Palermo, E. Grossi, I. Karpenko, and F. Becattini, Eur. Phys. J. C 84, 920 (2024), arXiv:2404.14295 [nucl-th]
-
[9]
Relativistic fluid dynamics with spin
W. Florkowski, B. Friman, A. Jaiswal, and E. Speranza, Phys. Rev. C97, 041901 (2018), arXiv:1705.00587 [nucl-th] . REFERENCES21
work page Pith review arXiv 2018
- [10]
-
[11]
Relativistic first-order spin hydrodynamics via the Chapman-Enskog expansion,
J. Hu, Phys. Rev. D105, 076009 (2022), arXiv:2111.03571 [hep-ph]
-
[12]
Relativistic dissipative spin dynamics in the relaxation time approximation,
S. Bhadury, W. Florkowski, A. Jaiswal, A. Kumar, and R. Ryblewski, Phys. Lett. B814, 136096 (2021), arXiv:2002.03937 [hep-ph]
-
[13]
S. Bhadury, W. Florkowski, A. Jaiswal, A. Kumar, and R. Ryblewski, Phys. Rev. Lett.129, 192301 (2022), arXiv:2204.01357 [nucl-th]
-
[14]
Kinetic theory for massive spin-1/2 particles from the Wigner-function formalism,
N. Weickgenannt, X.-L. Sheng, E. Speranza, Q. Wang, and D. H. Rischke, Phys. Rev. D100, 056018 (2019), arXiv:1902.06513 [hep-ph]
-
[15]
Generating Spin Polarization from Vorticity through Nonlocal Collisions,
N. Weickgenannt, E. Speranza, X.-l. Sheng, Q. Wang, and D. H. Rischke, Phys. Rev. Lett.127, 052301 (2021), arXiv:2005.01506 [hep- ph]
-
[16]
N. Weickgenannt, E. Speranza, X.-l. Sheng, Q. Wang, and D. H. Rischke, Phys. Rev. D104, 016022 (2021), arXiv:2103.04896 [nucl-th]
-
[17]
Relativistic second-order dissipative spin hydrodynamics from the method of moments,
N. Weickgenannt, D. Wagner, E. Speranza, and D. H. Rischke, Phys. Rev. D106, 096014 (2022), arXiv:2203.04766 [nucl-th]
-
[18]
N. Weickgenannt and J.-P. Blaizot, Phys. Rev. D109, 056012 (2024), arXiv:2311.15817 [hep-ph]
- [19]
-
[20]
S. Banerjee, S. Bhadury, W. Florkowski, A. Jaiswal, and R. Ryblewski, Phys. Rev. C111, 064912 (2025), arXiv:2405.05089 [hep-ph]
-
[21]
Relativistic spin hydrodynamics with momentum and spin-dependent relaxation time,
S. Bhadury, Phys. Rev. C111, 034909 (2025), arXiv:2408.14462 [hep- ph]
- [22]
-
[23]
Fate of spin polarization in a relativistic fluid: An entropy-current analysis
K. Hattori, M. Hongo, X.-G. Huang, M. Matsuo, and H. Taya, Phys. Lett. B795, 100 (2019), arXiv:1901.06615 [hep-th]
work page Pith review arXiv 2019
-
[24]
K. Fukushima and S. Pu, Phys. Lett. B817, 136346 (2021), arXiv:2010.01608 [hep-th] . 22REFERENCES
-
[25]
Relativistic second-order spin hydrodynamics: An entropy-current analysis,
R. Biswas, A. Daher, A. Das, W. Florkowski, and R. Ryblewski, Phys. Rev. D108, 014024 (2023), arXiv:2304.01009 [nucl-th]
-
[26]
Causality and stability analysis for the minimal causal spin hydrodynamics,
X.-Q. Xie, D.-L. Wang, C. Yang, and S. Pu, Phys. Rev. D108, 094031 (2023), arXiv:2306.13880 [hep-ph]
- [27]
- [28]
-
[29]
Stability and causality of rest frame modes in second-order spin hydrodynamics,
A. Daher, W. Florkowski, R. Ryblewski, and F. Taghinavaz, Phys. Rev. D109, 114001 (2024), arXiv:2403.04711 [hep-ph]
- [30]
-
[31]
The ideal relativistic fluid limit for a medium with polarization
D. Montenegro, L. Tinti, and G. Torrieri, Phys. Rev. D96, 056012 (2017), [Addendum: Phys.Rev.D 96, 079901 (2017)], arXiv:1701.08263 [hep-th]
work page Pith review arXiv 2017
-
[32]
Linear response theory and effective action of relativistic hydrodynamics with spin,
D. Montenegro and G. Torrieri, Phys. Rev. D102, 036007 (2020), arXiv:2004.10195 [hep-th]
-
[33]
Perfect spinfluid: A divergence-type approach,
N. Abboud, L. Gavassino, R. Singh, and E. Speranza, Phys. Rev. D 112, 094043 (2025), arXiv:2506.19786 [nucl-th]
-
[34]
Nonlinear causality and stability of perfect spin hydrodynamics and its nonperturbative character,
S. Bhadury, Z. Drogosz, W. Florkowski, S. K. Kar, and V. Mykhaylova, (2025), arXiv:2511.19295 [hep-ph]
-
[35]
Relativistic hydrodynamics for spin-polarized fluids,
W. Florkowski, A. Kumar, and R. Ryblewski, Prog. Part. Nucl. Phys. 108, 103709 (2019), arXiv:1811.04409 [nucl-th]
-
[36]
Polarization and Vorticity in the Quark–Gluon Plasma,
F. Becattini and M. A. Lisa, Ann. Rev. Nucl. Part. Sci.70, 395 (2020), arXiv:2003.03640 [nucl-ex]
-
[37]
Huang,An introduction to relativistic spin hydrodynamics,Nucl
X.-G. Huang, Nucl. Sci. Tech.36, 208 (2025), arXiv:2411.11753 [nucl- th]
-
[38]
W. Florkowski, J. Subatomic Part. Cosmol.3, 100028 (2025), arXiv:2411.19673 [hep-ph]
-
[39]
W. Florkowski and M. Hontarenko, Phys. Rev. Lett.134, 082302 (2025), arXiv:2405.03263 [hep-ph] . REFERENCES23
-
[40]
Hybrid approach to perfect and dissipative spin hydrodynamics,
Z. Drogosz, W. Florkowski, and M. Hontarenko, Phys. Rev. D110, 096018 (2024), arXiv:2408.03106 [hep-ph]
-
[41]
Mathisson, Acta Phys
M. Mathisson, Acta Phys. Polon.6, 163 (1937)
1937
-
[42]
Mathisson, General Relativity and Gravitation42, 1011 (2010)
M. Mathisson, General Relativity and Gravitation42, 1011 (2010)
2010
-
[43]
Local equilibrium Wigner function for spin-1/2 particles,
S. Bhadury, Z. Drogosz, W. Florkowski, S. K. Kar, and V. Mykhaylova, (2025), arXiv:2505.02657 [hep-ph]
-
[44]
Drogosz, Physics7, 31 (2025), arXiv:2504.03396 [hep-ph]
Z. Drogosz, Physics7, 31 (2025), arXiv:2504.03396 [hep-ph]
- [45]
-
[46]
Perfect spin hydrodynamics at all orders in spin polarization,
Z. Drogosz, Phys. Lett. B873, 140205 (2026), arXiv:2509.06014 [hep- ph]
-
[47]
Application range of perfect spin hydrodynamics,
Z. Drogosz, W. Florkowski, and V. Mykhaylova, Phys. Rev. D112, L051901 (2025), arXiv:2506.01537 [hep-ph]
- [48]
-
[49]
Spin polarization evolution in a boost invariant hydrodynamical background
W. Florkowski, A. Kumar, R. Ryblewski, and R. Singh, Phys. Rev. C 99, 044910 (2019), arXiv:1901.09655 [hep-ph]
work page Pith review arXiv 2019
-
[50]
Z. Drogosz, W. Florkowski, N. Lygan, and R. Ryblewski, Phys. Rev. C111, 024909 (2025), arXiv:2411.06154 [hep-ph]
- [51]
- [52]
-
[53]
J. D. Bjorken, Phys. Rev. D27, 140 (1983)
1983
-
[54]
Florkowski,Phenomenology of Ultra-Relativistic Heavy-Ion Colli- sions(Singapore: World Scientific, 2010)
W. Florkowski,Phenomenology of Ultra-Relativistic Heavy-Ion Colli- sions(Singapore: World Scientific, 2010)
2010
-
[55]
Coleman,Lectures of Sidney Coleman on Quantum Field Theory, edited by B
S. Coleman,Lectures of Sidney Coleman on Quantum Field Theory, edited by B. G.-g. Chen, D. Derbes, D. Griffiths, B. Hill, R. Sohn, and Y.-S. Ting (WSP, Hackensack, 2018)
2018
- [56]
- [57]
-
[58]
Spin-dependent distribution functions for relativistic hydrodynamics of spin-1/2 particles
W. Florkowski, B. Friman, A. Jaiswal, R. Ryblewski, and E. Speranza, Phys. Rev. D97, 116017 (2018), arXiv:1712.07676 [nucl-th]
work page Pith review arXiv 2018
- [59]
- [60]
- [61]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.