A strongly hyperbolic viscous relativistic hydrodynamics theory with first-order charge current

This paper extends the Bemfica-Disconzi-Noronha-Kovtun (BDNK) first-order dissipative relativistic hydrodynamics model to include a full first-order charge current with out-of-equilibrium corrections, demonstrating that this inclusion ensures a strongly hyperbolic, causal, and stable system coupled to Einstein's equations without requiring additional frame restrictions.

The Gist

Imagine you are trying to predict the weather. You have a giant, swirling ocean of air (the atmosphere) that moves, heats up, cools down, and carries moisture. In physics, we do something similar with fluids like water, gas, or even the super-hot soup of particles inside a neutron star. We use math to describe how these fluids flow.

For a long time, physicists had two main ways to describe these fluids:

1. The "Perfect" Fluid: Imagine a fluid with zero friction and no heat loss. It flows like a dream, but it's a fantasy. Real fluids get hot, they stick together (viscosity), and they conduct electricity.
2. The "Real" Fluid: This accounts for friction and heat. But for decades, the math used to describe real fluids was broken. It was like a weather forecast that said, "A storm will arrive in 5 minutes, but also 500 years from now, and it might travel faster than light." This made the math unstable and useless for computers.

The Problem: The "Charge" Glitch
Recently, a new, better way of doing this math was invented (called the BDNK theory). It fixed the broken math for fluids that just flow and get hot. But there was a catch: it didn't handle electric charge very well.

Think of a fluid as a busy highway.

• Energy is the cars.
• Heat is the traffic jams.
• Charge is the specific type of car (like red sports cars).

The old "perfect" math said: "The red cars just follow the traffic."
The new "real" math (BDNK) said: "The red cars can also drift sideways if the road is bumpy."

However, the previous version of this new math had a bug. When they tried to calculate how the "red cars" (charge) moved, the math created a "ghost lane." This ghost lane didn't actually move anything; it just sat there, confusing the computer. It made the system unstable, like a bridge that looks solid but collapses if you tap it.

The Solution: The "Ghost-Busting" Fix
The authors of this paper, Federico and Fernando, found a way to fix the "ghost lane."

They realized that the way the "red cars" (charge) move depends on how the whole highway (the fluid) is evolving. They added a new rule: "The movement of charge must be proportional to the movement of the fluid itself."

The Analogy of the Dance Floor
Imagine a crowded dance floor (the fluid).

• The Old Way: People were dancing, and someone shouted, "Charge!" The red-shirted dancers tried to move, but they didn't know how to move relative to the crowd. They froze or moved randomly, causing a panic (instability).
• The New Way: The authors realized that the red-shirted dancers should move in sync with the rhythm of the whole room. If the room spins left, the red shirts spin left. If the room speeds up, they speed up. By tying the charge's movement directly to the fluid's rhythm, the "ghost lane" disappears. Everyone moves smoothly, and the dance floor becomes stable.

Why This Matters

1. No Faster-Than-Light Signals: The math now guarantees that nothing travels faster than light. The "ghost" that was breaking the speed limit is gone.
2. Stability: The equations are now "strongly hyperbolic." In plain English, this means if you give a computer a starting picture of the fluid, it can predict the future without the numbers blowing up into infinity. It's like a weather model that actually works.
3. Real-World Applications: This is crucial for understanding:
   • Neutron Stars: These are dead stars so dense that a teaspoon weighs a billion tons. They spin fast and have huge magnetic fields (charge). To understand what happens when two of them crash (creating gravitational waves), we need this new math.
   • Black Holes: How gas swirls into a black hole.
   • The Big Bang: How the early universe expanded.

The "Frame" Choice
The paper also talks about choosing a "frame." Imagine you are watching a race.

• If you stand on the sidelines, the runners look fast.
• If you run alongside them, they look slow.

In physics, you can choose how you define "moving." The authors found a specific "viewpoint" (a frame) where the math works perfectly for almost any type of fluid, whether it's hot gas or dense star matter. They proved that if you pick this specific viewpoint, the fluid will always behave physically: it won't create energy out of nothing, and it will always obey the laws of thermodynamics (entropy always increases).

In Summary
This paper is like a mechanic fixing a high-performance race car engine. The engine (the math) was powerful but had a glitch when you added a turbocharger (electric charge). The mechanic (the authors) rewired the fuel injection system so that the turbo works perfectly with the engine. Now, the car can race at the speed of light (metaphorically) without falling apart, allowing scientists to simulate the most violent events in the universe with confidence.

Technical Summary

Here is a detailed technical summary of the paper "A strongly hyperbolic viscous relativistic hydrodynamics theory with first-order charge current" by Schianchi and Abalos.

1. Problem Statement

Relativistic hydrodynamics is essential for modeling high-energy physics phenomena (e.g., heavy-ion collisions, neutron star mergers). While ideal fluid models are well-posed, dissipative models face significant challenges:

• Causality and Stability: Classical first-order theories (Eckart, Landau-Lifshitz) are parabolic, leading to acausal signal propagation and instability. Second-order theories (Müller-Israel-Stewart, MIS) fix this but introduce extra dynamical variables and stiff relaxation equations, complicating numerical simulations.
• The BDNK Framework: The Bemfica-Disconzi-Noronha-Kovtun (BDNK) theory offers a first-order dissipative model that is causal, stable, and strongly hyperbolic without extra variables. However, previous formulations (specifically Ref. [20]) treated the charge current ($J^\mu$) as that of an ideal fluid.
• The Specific Deficit: In the ideal-charge-current approximation of Ref. [20], the system becomes weakly hyperbolic (ill-posed) due to a degenerate null eigenvalue in the principal matrix. Furthermore, the resulting evolution equation for charge conservation is not in a conservative form, making it unsuitable for handling discontinuities (shocks) in numerical relativity.

Goal: Extend the BDNK framework to include a fully first-order charge current with non-equilibrium corrections proportional to the ideal fluid evolution equations. The goal is to achieve a strongly hyperbolic, causal, stable, and entropy-positive system that remains in a conservative form suitable for numerical simulation.

2. Methodology

A. Theoretical Formulation

The authors propose a new constitutive relation for the charge current $J^\mu$ and stress-energy tensor $T^{\mu\nu}$.

• Variables: The system evolves the energy density ($\epsilon$), charge number density ($n$), and fluid four-velocity ($u^\mu$).
• Constitutive Relations: The dissipative terms ($A, \Pi, Q^\mu, N, J^\mu$) are expanded to first order in gradients. Crucially, the non-equilibrium terms are constructed to be proportional to the ideal fluid evolution equations.
  • The charge current includes terms proportional to $\nabla_\lambda (\mu/T)$ (diffusion) and terms proportional to the acceleration and pressure gradients (advection).
  • This specific structure ensures that the evolution equations for $\epsilon, n, u^\mu$ become second-order partial differential equations (PDEs) that can be cast in a conservative form (flux-conservative).
• Frame Choice: The theory allows for freedom in choosing the "frame" (definitions of $u^\mu, T, \mu$). The authors define a specific class of frames parameterized by dimensionless coefficients ($\hat{\tau}, \hat{\sigma}, \tilde{\sigma}$, etc.) to satisfy physical constraints.

B. Mathematical Analysis Tools

• Matrix Pencil Theory: Instead of performing an explicit first-order reduction (which increases algebraic complexity), the authors utilize a newly developed technique (Ref. [42]) to analyze the strong hyperbolicity of second-order PDE systems directly. This involves analyzing the principal symbol (matrix pencil) $N_{\mu\nu} l^\mu l^\nu$.
• Hyperbolicity Conditions: They derive conditions ensuring:
  1. Real Characteristic Velocities: All signal speeds are real.
  2. Multiplicity: The dimension of the kernel matches the algebraic multiplicity of the eigenvalues.
  3. Uniformity: Eigenvectors depend continuously on the wave vector.
• Stability Analysis: Linear stability around thermal equilibrium is analyzed using the Routh-Hurwitz theorem applied to the characteristic polynomials of the shear and sound channels.
• Entropy Production: The divergence of the entropy current is calculated to ensure the Second Law of Thermodynamics ($\nabla_\mu S^\mu \geq 0$) holds up to second-order gradients.

C. Coupling to Gravity

The authors analyze the coupling of this hydrodynamic system to the Einstein Field Equations (EFE) in the generalized harmonic gauge. They prove that the coupled system remains strongly hyperbolic provided the fluid characteristic speeds do not degenerate with the gravitational light cone.

3. Key Contributions

1. Resolution of Weak Hyperbolicity: The paper demonstrates that the previous BDNK formulation (Ref. [20]) is only weakly hyperbolic for charged fluids due to a specific degeneracy ($D=0$ in the characteristic polynomial). By including first-order terms in the charge current, this degeneracy is lifted, restoring strong hyperbolicity without imposing restrictive extra conditions on the transport coefficients.
2. Conservative Formulation: The new formulation yields a set of second-order evolution equations that can be written in a flux-conservative form. This is critical for numerical relativity, as it allows the use of high-resolution shock-capturing (HRSC) schemes to handle discontinuities correctly.
3. Frame Parametrization: The authors identify a specific family of frames (Proposition 7) where the system is simultaneously:
   • Strongly Hyperbolic: Well-posed initial value problem.
   • Causal: Signal speeds $< c$.
   • Linearly Stable: Perturbations decay over time.
   • Thermodynamically Consistent: Positive entropy production.
4. Generalized Equation of State (EoS) Compatibility: The analysis is not limited to ideal gases. The authors show that the conditions hold for piecewise polytropic EoS with thermal contributions, which are standard in neutron star merger simulations.

4. Key Results

• Hyperbolicity Conditions: The system is strongly hyperbolic if the discriminant of the sound channel polynomial is positive ($\Delta > 0$) and specific inequalities involving transport coefficients ($\tau, \beta, \lambda$) are met.
• Causality: All characteristic velocities (shear and sound modes) are proven to be less than the speed of light ($c=1$) under the derived frame constraints.
• Entropy: The entropy production is non-negative up to $O(\nabla^3)$ provided the bulk viscosity ($\zeta$), shear viscosity ($\eta$), and effective thermal conductivity ($\sigma + \sigma_0$) are non-negative.
• Numerical Feasibility: Through numerical sampling of the parameter space (using a provided Python script), the authors confirm that there exists a large region of parameters (specifically large relaxation times $\hat{\tau}$) where all causality, stability, and hyperbolicity conditions are satisfied for realistic EoS parameters found in binary neutron star mergers.
• Coupling to GR: The coupled EFE-BDNK system preserves strong hyperbolicity as long as the fluid's effective light cones remain strictly inside the spacetime light cone.

5. Significance

• Numerical Relativity: This work provides the mathematical foundation for implementing charged, dissipative relativistic hydrodynamics in numerical codes. Previous attempts were hindered by ill-posedness or non-conservative forms. This formulation enables robust simulations of charged neutron star mergers and accretion disks where charge diffusion and viscosity are relevant.
• Theoretical Advancement: It resolves a long-standing issue in the BDNK framework regarding charged fluids, proving that a first-order theory can be fully well-posed without the need for the auxiliary variables of MIS theory.
• Practical Application: The identification of a specific "safe" frame (Proposition 7) gives practitioners a concrete set of parameters to use in simulations, bridging the gap between abstract theoretical constraints and practical computational implementation.
• Future Work: The authors outline the next steps: implementing the formulation in a flux-conservative first-order reduction (to handle shocks explicitly) and testing it in full 3D general relativistic simulations of binary neutron star mergers.

In summary, this paper successfully extends the BDNK theory to include charge diffusion in a way that is mathematically rigorous (strongly hyperbolic), physically consistent (causal, stable, entropic), and numerically viable (conservative form), opening the door for realistic simulations of charged relativistic fluids in astrophysics.