⚛️ nuclear theory

Relativistic viscous hydrodynamics in the Density Frame for heavy-ion collisions

This paper presents a robust numerical implementation of relativistic viscous hydrodynamics in the Density Frame, demonstrating its validity against kinetic theory and showing that, despite differences in bulk-viscous dynamics compared to traditional formulations, it yields final-state observables consistent with standard codes and is ready for Bayesian analysis of heavy-ion collision data.

Original authors: Jean-François Paquet, Mayank Singh, Derek Teaney

Published 2026-11-02

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jean-François Paquet, Mayank Singh, Derek Teaney

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the heart of the universe's most violent collisions, where atomic nuclei smash together at nearly the speed of light, matter behaves in ways that seem impossible under normal conditions. When protons or lead nuclei collide at facilities like the Large Hadron Collider, they create a fleeting, super-hot soup of subatomic particles known as the quark-gluon plasma. For a fraction of a second, this plasma flows like a fluid, expanding and cooling before freezing into the particles that detectors eventually record. To understand this process, physicists rely on a set of mathematical rules called hydrodynamics, which describe how fluids move and change. However, applying these rules to the extreme, relativistic speeds of these collisions has long been a challenge. The standard methods used to model this fluid often require adding extra, invisible variables to keep the math stable and prevent the equations from breaking down. These extra variables act like a safety net, but they introduce a layer of complexity that makes it difficult to know exactly which properties of the fluid are being measured and which are just artifacts of the mathematical model.

A team of researchers has now developed a new way to run these simulations, stripping away the unnecessary safety nets to create a cleaner, more direct model of the fluid's behavior. In a recent study, Jean-François Paquet, Mayank Singh, and Derek Teaney introduced a method they call the "Density Frame." Instead of tracking the fluid's flow and its internal stresses as separate, evolving things, this new approach treats the fluid's energy and momentum as the only things that change over time. Everything else, including the internal friction that resists the flow, is calculated instantly based on how the fluid is moving at that exact moment. This eliminates the need for extra, non-physical variables that have cluttered previous models. The researchers tested this new engine by comparing it against known solutions for how particles behave in a vacuum and by running it alongside the most widely used simulation software in the field, known as MUSIC. They found that while the new method smooths out tiny, sharp ripples in the fluid more quickly than the old one, the final results for the particles that emerge from the collision are remarkably similar. This suggests that the complex extra steps in the old models might not be necessary to understand the big picture of how the quark-gluon plasma flows.

The core of this work lies in how the researchers handle the concept of entropy, or disorder, within the fluid. In any real fluid, friction creates heat and increases disorder, a process that must happen continuously. The new method is built around a mathematical tool called a "dissipative function," which acts like a master ledger for this disorder. By designing their computer code so that every step of the simulation automatically increases this disorder in a consistent way, the researchers ensured that their model remains stable even when the fluid becomes very thin or turbulent. This stability allowed them to simulate the collision without needing the artificial "regulators" that other codes require to stop the math from crashing. In their tests, they started with a specific type of initial condition called "IP-Glasma," which mimics the chaotic state of the nuclei just before they collide. They ran the simulation forward in time, watching how the temperature and flow velocity evolved.

When they compared the new Density Frame results with the standard MUSIC code, they noticed a distinct difference in how the two handled the fluid's internal friction. In the standard code, the fluid's resistance to flow, known as shear stress, takes a little bit of time to adjust to changes in the flow, much like a heavy object that takes a moment to start moving or stop. In the new Density Frame, this resistance adjusts instantly. This difference was even more pronounced in the bulk viscosity, which is the fluid's resistance to being squeezed or stretched. In the standard model, this resistance relaxes slowly over time, especially near the point where the fluid transitions from a hot plasma to a cooler gas. In the new model, it responds immediately. Despite these differences in the middle of the process, the final outcome was surprisingly consistent. When the researchers looked at the particles produced at the end of the collision, such as the total number of particles created, their average speed, and how they spread out in different directions, the two methods produced nearly identical results.

The study also explored what happens when the fluid has no resistance to being squeezed, a scenario where the bulk viscosity is set to zero. In this case, the two methods agreed almost perfectly. When they added a temperature-dependent bulk viscosity, which mimics the complex behavior of the fluid near its transition point, small differences appeared, particularly in the outer edges of the collision zone. However, these differences were generally smaller than the effect of the bulk viscosity itself. The researchers concluded that the complex, non-hydrodynamic modes that the old models rely on to stay stable do not seem to play a major role in determining the final observable properties of the collision. This finding is significant because it suggests that the current uncertainty in measuring the fluid's properties might not be due to the physics of the fluid itself, but rather to the extra mathematical baggage carried by the simulation tools.

By removing these extra variables, the Density Frame offers a more direct path to understanding the fundamental properties of the quark-gluon plasma. The researchers have made their code public and documented, allowing other scientists to use it as a replacement for the hydrodynamic engines in existing simulation packages. This opens the door for a new kind of analysis where the only parameters being adjusted are the actual physical properties of the fluid, such as its viscosity and its equation of state, rather than the arbitrary time delays built into the mathematical models. The authors suggest that this cleaner approach will be particularly useful for studying smaller collision systems, where the fluid is less developed and the influence of these extra mathematical modes might be more pronounced. Ultimately, this work provides a more robust and transparent tool for physicists to extract the true transport properties of the quark-gluon plasma from the data collected at high-energy colliders, bringing us closer to a complete understanding of the universe's most extreme fluid.

Technical Summary: Relativistic Viscous Hydrodynamics in the Density Frame for Heavy-Ion Collisions

Problem Statement
Relativistic viscous hydrodynamics is the standard framework for describing collective flow in heavy-ion collisions at RHIC and the LHC. However, conventional formulations face conceptual and practical challenges. The standard Landau or Eckart frame Navier-Stokes theory is unstable and acausal. The widely used Müller-Israel-Stewart (MIS) and BDNK formulations resolve causality by introducing non-hydrodynamic modes (dynamical fields like shear stress $\pi^{\mu\nu}$ and bulk pressure $\Pi$ ) that relax over finite timescales ( $\tau_\pi, \tau_\Pi$ ). These extra modes introduce non-universal parameters (e.g., relaxation times) that lie outside the realm of universal hydrodynamics, potentially obscuring the extraction of fundamental transport coefficients (shear and bulk viscosities) from experimental data, particularly in small systems.

Methodology
The authors present a numerical implementation of the Density Frame, a formulation of relativistic viscous hydrodynamics that is first-order in time, requires no auxiliary dynamical fields, and possesses no non-hydrodynamic modes. In this frame, the conserved energy ( $E = T^{00}$ ) and momentum ( $M^i = T^{0i}$ ) densities on a spatial slice are the sole slow variables. All other quantities, including the spatial stress tensor $T^{ij}$ , are determined instantaneously via a gradient expansion.

The numerical method employs a discrete analog of the dissipative function $R$ , a positive semi-definite quadratic form in the strains that governs entropy production.

Discretization: The viscous stresses are implemented as functional derivatives of a discrete dissipative function $R$ with respect to the hydrodynamic variables at lattice sites.
Algorithm: The authors utilize an Implicit-Explicit (IMEX) time-stepping scheme. The ideal (non-linear) evolution is treated explicitly using a Kurganov-Tadmor scheme, while the viscous evolution is treated implicitly.
Entropy Guarantee: By constructing the viscous step as the derivative of a discrete $R$ , the algorithm guarantees that entropy is produced on all length scales in the semi-discrete limit, leading to a remarkably robust numerical stability. This allows the simulation of dilute plasma edges with only a minimal energy-density pedestal, avoiding the complex regulators often required in MIS solvers.

Key Contributions and Benchmarks

2+1D Extension: The authors extend their previous 1+1D work to 2+1 dimensions in Bjorken coordinates, establishing the Density Frame as a fully developed tool for heavy-ion collision simulations.
Kinetic Theory Benchmark (Gubser Flow): The method is tested against the exact solution of kinetic theory in the Relaxation Time Approximation (RTA) with Gubser symmetry.
- The authors clarify the matching between kinetic theory and hydrodynamics by organizing the solution via an opacity expansion ( $\chi_1$ ).
- They demonstrate how the Bjorken attractor passes smoothly onto the Gubser solution order-by-order in gradients.
- The Density Frame solution, evolved using Bjorken time $\tau$ (where the fluid flows radially), agrees with the Gubser hydrodynamic solution (where the fluid is at rest in de Sitter time $\rho$ ) up to controlled gradient corrections, validating the frame's approximate foliation independence and Lorentz invariance.
Comparison with MUSIC (MIS Framework): The Density Frame engine replaces the hydrodynamic solver in the MUSIC code and is compared against the standard MIS-based MUSIC for IP-Glasma initial conditions in Pb-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

Results

Dissipation of Gradients: The Density Frame smooths sharp gradients and small-scale structures more rapidly than MIS. In the Density Frame, viscous stresses respond instantaneously to flow gradients, whereas in MIS, they relax over a finite time $\tau_\pi$ .
Bulk Viscous Sector: The most significant difference arises in the bulk sector. In the Density Frame, bulk pressure responds instantaneously to expansion ( $\Pi = -\zeta \partial_\mu u^\mu$ ). In MIS, bulk pressure relaxes over $\tau_\Pi$ , which can be long near the QCD crossover, causing a lag in the response.
Final State Observables: Despite differences in the evolution of local fields:
- With shear viscosity alone, the final-state hadron multiplicity, mean transverse momentum ( $\langle p_T \rangle$ ), and elliptic flow ( $v_2$ ) from both codes are very similar. The non-hydrodynamic modes of MIS have little effect on these observables.
- With bulk viscosity, differences are more pronounced (e.g., $v_2$ in peripheral collisions is $\sim 10\%$ lower in the Density Frame), but these differences remain smaller than the overall effect of including bulk viscosity itself.
- The authors note that the Density Frame yields slightly lower multiplicity ( $\sim 5\%$ ) and $\langle p_T \rangle$ ( $\sim 2-3\%$ ) compared to MUSIC when bulk viscosity is included, but these shifts are correlated and could be compensated by re-calibrating model parameters.

Significance
The paper claims that the Density Frame resolves the conceptual issues of non-universal parameters and non-hydrodynamic modes inherent in second-order hydrodynamics. By eliminating relaxation times and auxiliary fields, the framework reduces the parameter space to only the shear and bulk viscosities and the equation of state. The authors assert that the method is now ready for complete Bayesian analyses of heavy-ion data. This approach promises to constrain the transport coefficients of QCD more directly, without the theoretical uncertainties introduced by the choice of relaxation times or frame parameters. Furthermore, the framework is positioned to investigate the sensitivity to non-hydrodynamic modes in small systems (e.g., p+Pb, O+O) and to facilitate the study of stochastic hydrodynamics via Metropolis steps.

Technical Summary: Relativistic Viscous Hydrodynamics in the Density Frame for Heavy-Ion Collisions

More like this