Relativistic viscous hydrodynamics in the Density Frame for heavy-ion collisions
Abstract
We present a numerical study of relativistic viscous hydrodynamics in the Density Frame for boost-invariant flows in 2+1 dimensions. The Density Frame is a formulation of relativistic viscous hydrodynamics that is first order in time, requires no auxiliary fields, and has no non-hydrodynamic modes; the only parameters are the shear and bulk viscosities and the equation of state. We formulate the numerical method in terms of a discrete analog of the dissipative function, which guarantees that the discrete equations produce entropy on all length scales and leads to a remarkably robust algorithm. We test the method against the exact solution of kinetic theory in the relaxation time approximation with Gubser symmetry. In making this comparison, we clarify how the kinetic solution matches onto hydrodynamics: we organize the solution by the opacity and show how the Bjorken attractor passes onto the Gubser solution order by order in gradients. We then replace the hydrodynamic engine of the MUSIC code with the Density Frame, and compare the two for IP-Glasma initial conditions in Pb-Pb collisions at √(sNN) = 5.02 TeV. The Density Frame smooths sharp gradients more rapidly. The largest differences arise in the bulk-viscous sector, where the bulk pressure in the Density Frame responds instantaneously to the flow, while in the Müller-Israel-Stewart-based MUSIC it relaxes over a finite time. With shear viscosity alone, the multiplicity, mean transverse momentum, and elliptic flow of final-state hadrons from the two codes are very similar. The Density Frame is now ready for a complete Bayesian analysis of heavy-ion data.
Related
- Searching for the QCD critical point with relativistic viscous hydrodynamics
- Modeling Λ polarization in Au+Au collisions at √(sNN)=200 GeV using relativistic spin hydrodynamics
- Viscous hydrodynamics from correlated losses
- Levy-shaped femtoscopic sources from an exact solution of relativistic hydrodynamics
- Elliptic and triangular flow responses in O+O and Ne+Ne collisions from hydrodynamics and kinetic theory
BibTeX
@misc{paquet2026relativistic,
title = {{Relativistic viscous hydrodynamics in the Density Frame for heavy-ion collisions}},
author = {Jean-François Paquet and Mayank Singh and Derek Teaney},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.40909},
eprint = {2610.40909},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.40909},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Relativistic viscous hydrodynamics in the Density Frame for heavy-ion collisions
AU - Jean-François Paquet
AU - Mayank Singh
AU - Derek Teaney
PY - 2026
DA - 2026/10/29
PB - arXiv
DO - 10.48550/arxiv.2610.40909
UR - https://arxiv.org/abs/2610.40909
AB - We present a numerical study of relativistic viscous hydrodynamics in the Density Frame for boost-invariant flows in 2+1 dimensions. The Density Frame is a formulation of relativistic viscous hydrodynamics that is first order in time, requires no auxiliary fields, and has no non-hydrodynamic modes; the only parameters are the shear and bulk viscosities and the equation of state. We formulate the numerical method in terms of a discrete analog of the dissipative function, which guarantees that the discrete equations produce entropy on all length scales and leads to a remarkably robust algorithm. We test the method against the exact solution of kinetic theory in the relaxation time approximation with Gubser symmetry. In making this comparison, we clarify how the kinetic solution matches onto hydrodynamics: we organize the solution by the opacity and show how the Bjorken attractor passes onto the Gubser solution order by order in gradients. We then replace the hydrodynamic engine of the MUSIC code with the Density Frame, and compare the two for IP-Glasma initial conditions in Pb-Pb collisions at √(s_NN) = 5.02 TeV. The Density Frame smooths sharp gradients more rapidly. The largest differences arise in the bulk-viscous sector, where the bulk pressure in the Density Frame responds instantaneously to the flow, while in the Müller-Israel-Stewart-based MUSIC it relaxes over a finite time. With shear viscosity alone, the multiplicity, mean transverse momentum, and elliptic flow of final-state hadrons from the two codes are very similar. The Density Frame is now ready for a complete Bayesian analysis of heavy-ion data.
ER -