Nonlinear Causality in First Order Relativistic Viscous Hydrodynamics

1. Introduction: Relativistic Hydrodynamics and Causality

Relativistic hydrodynamics provides a framework for describing collective phenomena in systems moving at velocities comparable to the speed of light. It finds applications in diverse fields including high-energy nuclear physics, astrophysics, and cosmology. In such systems, causalitythe principle that influences cannot propagate faster than lightstands as a fundamental requirement of any physical theory.

The description of relativistic fluids centers on the energy-momentum tensor T, which encodes energy and momentum densities and their fluxes. Conservation of this tensor, expressed as:

provides the equations of motion for the fluid. For an ideal (non-dissipative) fluid, T = (+p)uu + pg, where is the energy density, p is the pressure, u is the fluid four-velocity, and g is the spacetime metric.

Real fluids, however, exhibit dissipative processes such as viscosity and heat conduction. When these effects are incorporated, the resulting theory must satisfy multiple physical requirements simultaneously: it must describe dissipative phenomena correctly, respect general covariance, maintain thermodynamic stability, and crucially, preserve causal signal propagation. As we will explore, maintaining causality in relativistic viscous hydrodynamics proves to be particularly challenging, especially when considering nonlinear effects.

2. First-Order Theories of Relativistic Viscous Hydrodynamics

The natural extension of ideal relativistic hydrodynamics to include dissipative effects is analogous to the non-relativistic Navier-Stokes equations. In what is termed first-order relativistic viscous hydrodynamics, dissipative fluxes are assumed to be linearly related to thermodynamic forces such as velocity gradients and expansion rates.

In this framework, the energy-momentum tensor takes the form:

where denotes the shear stress tensor, represents the bulk viscous pressure, and = g + uu is the projector orthogonal to the fluid four-velocity. The dissipative fluxes are given by constitutive relations:

where is shear viscosity, is bulk viscosity, is the shear tensor (symmetric traceless part of velocity gradients), and = u is the expansion scalar. These relations represent direct generalizations of the non-relativistic Navier-Stokes equations.

While this approach appears natural and consistent with the structure of non-relativistic viscous hydrodynamics, it faces profound problems in the relativistic context. Most significantly, the resulting equations can support signal propagation faster than lightviolating one of the fundamental tenets of relativity.

3. Origins of Acausality in First-Order Theories

The causality problem in first-order relativistic viscous hydrodynamics stems directly from the momentless (or algebraic) nature of the constitutive relations. Because the dissipative fluxes respond instantaneously to local changes in thermodynamic gradients, the system behaves analogously to heat-conducting solids or fluids described by parabolic equations rather than hyperbolic ones.

In non-relativistic physics, this instantaneous response doesn't pose a problem because there is no fundamental speed limit. However, in relativistic contexts, these instantaneously adjusting dissipative fluxes create effective signal velocities that can exceed the speed of light. Mathematically, this manifests as the presence of characteristic speeds in the linearized equations that grow without bound as the dissipative coefficients increase.

The problem becomes particularly acute when examining the characteristic velocities of the theory. For certain combinations of parameters, these velocities can exceed c = 1, indicating acausal behavior. Hiscock and Lindblom (1985) demonstrated this issue explicitly, showing that for first-order theories, there always exist physical conditions where signal velocities exceed the speed of light.

4. Exacerbation by Nonlinear Effects

While linear analysis already reveals causality violations, the problem becomes significantly more severe when nonlinear effects are incorporated. These nonlinearities arise from several sources in relativistic viscous hydrodynamics:

• Geometric nonlinearities from the projector and other tensor structures
• Nonlinear dependence of transport coefficients on thermodynamic variables
• Nonlinear coupling between different dissipative processes
• Nonlinearities inherent in the equations of state relating pressure, energy density, and other thermodynamic quantities

These nonlinear effects create additional pathways for causality violations. In regions where gradients become largea common situation in high-energy nuclear collisions, astrophysical explosions, or the early universethe nonlinear terms can dominate the dynamics, leading to characteristic velocities that not only exceed but grow arbitrarily large.

An important concept in understanding nonlinear causality violations is that of the causality horizon. This is a surface in spacetime beyond which the theory ceases to provide causal evolution. In first-order relativistic viscous hydrodynamics, the position of this horizon depends on the local values of gradients and transport coefficients. As nonlinearities amplify in high-gradient regions, the horizon can shift unpredictably, potentially encompassing entire domains of interest.

Mathematically, the causality horizon can be understood by examining the characteristic determinant of the full nonlinear system. The condition for acausal behavior becomes det(g - vv) < 0, where v are the characteristic velocities. When v becomes spacelike (v² > 1), the determinant becomes negative, indicating acausal propagation. In nonlinear regimes, v can develop complex spatial dependence, leading to irregular causality horizons that are difficult to track in numerical simulations.

5. The Israel-Stewart Second-Order Formalism

The resolution to these causality problems lies in extending first-order theories to include additional termsspecifically, relaxation terms that give dissipative fluxes inertia. This approach was pioneered through the independent work of Israel (1976) and Stewart (1979), leading to what is now known as the Israel-Stewart formalism.

In the Israel-Stewart framework, the dissipative fluxes obey evolution equations rather than algebraic relations. For shear viscosity and bulk viscosity, these equations typically take the form:

where and are relaxation times for shear and bulk viscous processes, respectively. These equations indicate that dissipative fluxes respond to changes in thermodynamic forces not instantaneously but with a characteristic relaxation time, making the evolution hyperbolic rather than parabolic.

The inclusion of second-order termssuch as or ²adds further nonlinear corrections to the theory. While these terms are often neglected in practical applications, they can be important in regimes where dissipative fluxes become large. Importantly, their inclusion does not reintroduce causality violations provided the relaxation terms remain present.

6. Modern Extensions and Applications

Recent developments in relativistic viscous hydrodynamics have extended the Israel-Stewart formalism in several important directions. Significant advances include:

• Extended Israel-Stewart theories: Recent work by Denicol, Niemi, and others has identified additional second-order terms that were omitted in the original formulation. These terms contribute to the accuracy of the theory, particularly in far-from-equilibrium situations encountered in high-energy nuclear collisions.
• Effective field theory approaches: These methods treat hydrodynamics as a low-energy effective theory of a consistent underlying microscopic quantum field theory, systematically controlling approximations and ensuring that fundamental symmetries are respected.
• Holographic insights: Through the AdS/CFT correspondence in string theory, strongly coupled quantum field theories can be studied via dual gravitational systems. This holographic approach has yielded important insights into the microscopic origins of transport coefficients and the nature of causality constraints.
• Kinetic theory foundations: Studies of the Boltzmann equation and other kinetic frameworks provide microscopic derivations of hydrodynamic equations, helping to identify which terms are necessary for a consistent theory and how transport coefficients relate to microscopic dynamics.

In high-energy nuclear physics, where relativistic hydrodynamics is routinely used to model the quark-gluon plasma formed in heavy-ion collisions, the Israel-Stewart formalism has become standard practice. These simulations must handle extreme conditionstemperatures far exceeding nuclear densities, rapid expansion (expansion rates reaching 0.1-1 fm/c), and significant velocity gradientsmaking a fully consistent causal formulation essential.

7. Conclusion

Nonlinear causality violations in first-order relativistic viscous hydrodynamics represent a fundamental theoretical limitation that necessitated the development of more sophisticated frameworks. The transition from first-order to second-order theories, particularly through the Israel-Stewart formalism and its modern extensions, has resolved these causality problems while providing a more accurate description of relativistic dissipative systems.

As precision measurements in high-energy nuclear physics, astrophysics, and cosmology continue to improve, maintaining causal and consistent hydrodynamic descriptions will remain essential for interpreting observations and advancing our understanding of collective phenomena in relativistic systems. Future developments in this field will likely focus on refining our understanding of second-order transport coefficients, extending hydrodynamic theories to even more far-from-equilibrium situations, and strengthening connections between hydrodynamics and microscopic theories.