Nonlinear analysis of causality for heat flow in heavy-ion collisions: constraints from equation of state

This study analyzes the causal constraints of the Mueller-Israel-Stewart theory for heat flow in heavy-ion collisions, revealing that realistic conditions lead to causality-violating heat fluxes that suggest either an overestimation of transport coefficients or a breakdown of the fluid approximation.

The Gist

Imagine you are trying to predict the weather in a storm. You have a set of rules (equations) that tell you how wind, rain, and temperature move. Usually, these rules work great. But what if the storm gets so violent that your rules start predicting that a raindrop could travel faster than the speed of light? That would be a disaster for physics, because nothing can beat light speed.

This paper is about checking the "rules of the road" for a very specific, extreme kind of weather: the hot, dense soup of particles created when heavy atoms smash into each other in giant particle accelerators (like the RHIC at Brookhaven National Lab). This soup is called a quark-gluon plasma.

Here is a breakdown of the paper's story, using simple analogies:

1. The Problem: The "Too Fast" Heat

When these atoms smash together, they create a fireball that is incredibly hot in the center and cooler on the edges. Heat naturally wants to rush from the hot center to the cool edges. This rushing heat is called heat flux.

Physicists use a set of advanced math rules (called the Mueller-Israel-Stewart theory) to describe how this heat moves. These rules are designed to make sure nothing moves faster than light (causality). However, the author of this paper asked: "What happens if the heat is rushing so fast that our math breaks down?"

2. The Detective Work: Checking the "Speed Limit"

The author treated the math like a speed limit sign.

• The Goal: Find the "safe zone" where the heat moves slower than light.
• The Variables: The author looked at two main things that change the speed limit:
  1. The Equation of State (EoS): Think of this as the "personality" of the fluid. Is it stiff like a rock or squishy like jelly? The paper found that if the fluid is "stiffer," the heat can move faster without breaking the rules.
  2. Relaxation Time: Imagine you are driving a car. If you slam on the brakes, the car doesn't stop instantly; it takes a moment to react. This delay is "relaxation time." The paper found that if the fluid takes a longer time to react to temperature changes, it can handle much more heat flow without breaking the speed limit.

3. The Shocking Discovery: The "Impossible" Heat

The author then tried to calculate how much heat is actually rushing around in a real heavy-ion collision using the best estimates we have for how well this "soup" conducts heat (thermal conductivity).

The Result was terrifying:
The math predicted that the heat is rushing at a speed roughly 300 to 800 times faster than the energy density of the soup itself.

• Analogy: Imagine a tiny ant (the energy) trying to push a massive boulder (the heat). The math says the ant is pushing the boulder so hard that the boulder is flying away at 800 times the speed of light.

This is physically impossible. It's like trying to run a marathon at the speed of a bullet.

4. Why is this happening? (The Two Suspects)

The paper suggests two possible reasons for this impossible result:

• Suspect A: We are overestimating the "conductivity."
  We don't know exactly how well this particle soup conducts heat. We are guessing based on models. The author suspects our guess is way too high. If the soup is actually a terrible conductor of heat (like a thick wool sweater), the heat wouldn't rush so fast, and the math would work.
• Suspect B: The "Fluid" idea is broken.
  Maybe, under these extreme conditions, we can't even treat the particles as a smooth fluid anymore. It's like trying to describe a crowd of people as a single river of water. If the crowd is too chaotic, the "river" math stops working entirely.

5. The "Pressure" Correction

The author tried to fix the math by adding a "pressure gradient" term (think of it as the crowd pushing back against the heat). This helped a little—it reduced the crazy speed by about 15%—but it didn't fix the problem. The heat was still moving way too fast.

6. The Conclusion: We Need Better Data

The paper ends with a call to action. The current math says the heat flow in these collisions is "unrealistically large." This means either:

1. Our models for how heat moves through this soup are wrong (we need better data from supercomputers called Lattice QCD).
2. Or, the fluid description we are using is completely breaking down in these extreme moments.

In a nutshell:
This paper is a "safety check" for the physics of particle collisions. It found that if we use our current best guesses, the heat in these collisions is moving so fast that it violates the laws of the universe. This tells us we need to either find better ways to measure how heat moves in these tiny explosions, or admit that our current way of describing them is too simple for such extreme conditions.

Technical Summary

1. Problem Statement

The study addresses the issue of causality and hyperbolicity in relativistic dissipative fluid dynamics, specifically within the Mueller-Israel-Stewart (MIS) second-order theory for heat-conducting fluids.

• Context: In ultra-relativistic heavy-ion collisions (e.g., at RHIC and LHC), the system is often modeled as a fluid. However, standard first-order theories (Navier-Stokes) are acausal and unstable. While MIS theory resolves these issues in the linear regime, its behavior in the non-linear regime (specifically with large heat fluxes) remains less explored.
• Specific Challenge: The paper investigates whether the MIS framework remains causal (i.e., signal propagation speeds do not exceed the speed of light) when subjected to large heat fluxes ($q$) relative to energy density ($\varepsilon$).
• Frame Choice: The analysis focuses on the Eckart frame, where the fluid four-velocity is defined by the conserved baryon current. This frame is particularly relevant for heavy-ion collisions at lower beam energies (finite baryon density) where heat conduction is a distinct dissipative process, unlike in the Landau frame where it is often absorbed into diffusion currents.

2. Methodology

The author employs a combination of analytical formulation and numerical analysis to map the causal parameter space.

• Theoretical Framework:
  • The study utilizes the MIS equations for a one-dimensional fluid flow with planar symmetry in the Eckart frame.
  • Shear and bulk viscosities are neglected to isolate the causal structure associated specifically with heat conduction.
  • The system is described by four variables: baryon density ($n$), energy density ($\varepsilon$), fluid rapidity ($\rho$), and heat flux magnitude ($q$).
• Characteristic Analysis:
  • The governing equations are cast into a matrix form $A^\mu_\nu \partial_\mu Y^\nu + B = 0$.
  • Causality is determined by the hyperbolicity of the system. The system is causal if the eigenvalues of the characteristic matrix (representing wave propagation speeds, $v$) are real and subluminal ($|v| \leq 1$).
  • The characteristic equation (a quartic polynomial in $v$) is derived and solved numerically.
• Equation of State (EoS) Variations:
  • Constant $c_s^2$: The speed of sound squared ($c_s^2$) is varied (0.05, 0.10, 0.33, 0.50) to test sensitivity.
  • Lattice QCD EoS: A realistic, temperature-dependent EoS derived from Lattice QCD (Wuppertal-Budapest collaboration) is used to account for the non-monotonic behavior of $c_s^2$ near the QCD crossover temperature ($T_c \approx 150$ MeV).
• Relaxation Parameter ($\lambda$):
  • The second-order transport coefficient $\beta_1$ is parameterized as $\beta_1 = \lambda \frac{5}{4p}$. The study varies $\lambda$ (0.5, 0.75, 1.0, 10.0, 20.0) to assess the impact of the relaxation time on causality.
• Heat Flux Estimation:
  • To assess physical relevance, the author estimates the magnitude of heat flux ($|q|/\varepsilon$) in RHIC conditions using the Navier-Stokes approximation (limit $\tau_q \to 0$).
  • Thermal conductivity ($\kappa$) is estimated using kinetic theory models (BGK/quasiparticle models) and the NJL model, as first-principle lattice QCD results for $\kappa$ are currently unavailable.

3. Key Contributions

1. Mapping Nonlinear Causality: The paper provides a comprehensive map of the causal parameter space for the MIS theory in the Eckart frame, identifying regions of strong hyperbolicity, weak hyperbolicity, and non-hyperbolicity (acausality) as functions of $q/\varepsilon$, $c_s^2$, and $\lambda$.
2. EoS Dependence: It demonstrates that the causal domain is highly sensitive to the Equation of State. Specifically, a realistic Lattice QCD EoS (with a dip in $c_s^2$ near $T_c$) significantly narrows the causal window compared to constant-speed-of-sound models.
3. Quantitative Constraints on Transport Coefficients: The study quantifies how the relaxation time parameter ($\lambda$) influences the maximum allowable heat flux before causality is violated.
4. Physical Reality Check: By estimating actual heat flux values in heavy-ion collisions, the paper highlights a potential breakdown of the fluid approximation or an overestimation of transport coefficients in current models.

4. Key Results

• Causal Domain Expansion:
  • Stiffer EoS: Larger values of $c_s^2$ (stiffer equations of state) expand the range of $|q/\varepsilon|$ for which the system remains causal.
  • Relaxation Time: Increasing the relaxation parameter $\lambda$ (longer relaxation time) significantly enlarges the causal domain. For example, for $c_s^2 = 1/3$, the causal limit increases from $|q/\varepsilon| \lesssim 0.18$ ($\lambda=1$) to $|q/\varepsilon| \lesssim 0.25$ ($\lambda=10$).
• Lattice QCD Impact: When using the realistic Lattice QCD EoS, the causal region is more constrained than with constant $c_s^2$. The dip in the speed of sound around $T \approx 130$ MeV leads to a narrowing of the causal window for intermediate energy densities.
• Heat Flux Estimates (The "Crisis"):
  • Using kinetic theory estimates for thermal conductivity ($\kappa$) and typical RHIC parameters ($T_0 = 200$ MeV, $R=5$ fm), the calculated heat flux ratio is unrealistically large: $|q|/\varepsilon \approx 330$ (for $T=200$ MeV) and $\approx 811$ (for $T=150$ MeV).
  • These values are orders of magnitude larger than the causal limits ($|q/\varepsilon| \lesssim 0.25$).
  • Including pressure gradient corrections reduces the heat flow by only $\approx 15\%$, which is insufficient to resolve the causality violation.
• Energy Condition Violation: The estimated heat fluxes violate the weak energy condition ($|q|/(\varepsilon+p) > 1/2$). This implies that for some observers, the energy density would be negative, indicating a fundamental breakdown of the fluid description itself, independent of causality.

5. Significance and Implications

• Limitations of Current Models: The results suggest that either the current kinetic theory estimates for thermal conductivity ($\kappa$) in QCD matter are significantly overestimated, or the hydrodynamic approximation breaks down under the extreme temperature gradients found in heavy-ion collisions.
• Need for First-Principles Data: The paper strongly argues that first-principles calculations of thermal conductivity from Lattice QCD are essential to accurately constrain the causal parameter space and validate hydrodynamic models for heavy-ion collisions.
• Frame Independence of Causality: While the analysis is performed in the Eckart frame, the author notes that causality is a property of the constitutive structure, not the frame choice. Therefore, these constraints likely manifest in Landau-frame formulations at finite baryon density as well.
• Future Directions: The findings underscore the need to investigate higher-dimensional flows, more realistic EoS at high baryon chemical potentials, and the transition from non-equilibrium to hydrodynamic regimes to fully understand causal limits in relativistic heavy-ion collisions.

In summary, the paper reveals that under realistic heavy-ion collision conditions, standard estimates of heat flow lead to severe causality violations, suggesting a critical gap in our understanding of transport coefficients or the applicability of fluid dynamics in these extreme environments.