Neutron star evolution with the Bemfica-Disconzi-Noronha-Kovtun viscous hydrodynamics framework

This paper presents the first non-linear numerical simulation of spherically symmetric neutron stars using the causal and stable BDNK viscous hydrodynamics framework under the Cowling approximation, demonstrating stable evolutions and analyzing quasi-normal mode frequencies as a foundational step toward fully consistent astrophysical models.

The Gist

The Big Picture: Fixing the "Perfect Fluid" Problem

Imagine you are trying to predict how a giant, super-dense ball of dough (a neutron star) behaves when it gets shaken or squeezed. For decades, scientists have modeled these stars as "perfect fluids."

Think of a perfect fluid like honey that flows without any resistance. It's smooth, it doesn't stick to itself, and it has no internal friction. This is a great approximation for many things, but in reality, even honey has a little bit of stickiness (viscosity). When two neutron stars crash into each other, they get incredibly hot and chaotic. In this chaos, that "stickiness" (viscosity) might actually matter.

The problem is that the old math used to describe sticky fluids (called MIS theory) is broken. It's like trying to drive a car with a steering wheel that sometimes spins 360 degrees on its own. The math predicts things that can't happen in the real universe (like signals traveling faster than light), making the simulations unstable and unreliable.

The New Tool: The BDNK Framework

This paper introduces a new, better way to do the math, called the BDNK framework (named after the four scientists who built it: Bemfica, Disconzi, Noronha, and Kovtun).

The Analogy:
Imagine you are trying to navigate a ship through a storm.

• The Old Way (MIS): You have a map that works great in calm water, but as soon as the waves get high, the map tells you to sail into a cliff. The ship crashes.
• The New Way (BDNK): This is a brand-new GPS system. It was built from the ground up to handle the storm. It guarantees that your ship will stay on course, won't sail faster than light, and won't crash into mathematical cliffs.

The authors of this paper took this new GPS (BDNK) and used it to simulate a single, spinning neutron star for the first time.

What They Did: The "Cowling" Test Drive

To test this new system, they didn't try to simulate two stars crashing (which is like trying to drive a car while it's being hit by a tornado). Instead, they did a "test drive" in a controlled environment.

1. The Setup (Cowling Approximation): They kept the star's gravity fixed. Imagine the star is a heavy bowling ball sitting on a trampoline. They didn't let the trampoline bounce or warp; they just watched how the ball itself wobbled. This simplifies the math so they could focus on the fluid physics.
2. The Experiment: They took a stable neutron star and gave it a tiny "nudge" (using the tiny errors inherent in computer code). They then watched how the star wobbled and settled back down.
3. The Variables: They tested different "recipes" for the star's stickiness. Some stars were very sticky (high viscosity), some were less sticky, and some had different types of stickiness (shear vs. bulk).

The Results: What Did They Find?

1. The Math Works!
The most important result is that the simulation didn't crash. They were able to run the star for a very long time (8,000 "time units") without the numbers blowing up. This proves that the BDNK framework is stable and ready for real-world use.

2. The "Ring" Didn't Change Much
When you hit a bell, it rings at a specific pitch. When they "hit" the neutron star, it vibrated at a specific frequency (called a Quasi-Normal Mode).

• Surprise: They found that adding "stickiness" (viscosity) didn't change the pitch of the ring very much. The star still sang the same note, whether it was sticky or not.
• Why? The main vibration happens on a large scale (like the whole bell ringing). The stickiness mostly affects tiny, short-scale ripples. So, the big note stays the same.

3. The "Fade" Got Faster
While the pitch didn't change, the volume did.

• In a perfect fluid (no stickiness), the star would ring forever (or until the computer errors stopped it).
• In a sticky fluid, the energy gets eaten up by friction. The star's vibration died out faster.
• The Discovery: The more "sticky" the star was, the faster the vibrations faded away. This is the key signature scientists will look for in real gravitational waves. If we see a neutron star ring and then go silent very quickly, it tells us the star is very viscous.

Why This Matters for the Future

This paper is like building the engine for a new race car. They haven't driven it around the whole track yet (they haven't simulated two stars crashing with full gravity), but they proved the engine doesn't explode when you turn the key.

The Future Goal:
Scientists want to use this new math to interpret the signals from future gravitational wave detectors (like the Einstein Telescope).

• The Dream: When we hear two neutron stars crash, we will hear a "chirp." By analyzing how that chirp fades out, we can use this BDNK math to figure out exactly how "sticky" the matter inside the star is.
• The Payoff: This will tell us about the fundamental laws of physics inside the densest objects in the universe, helping us understand how heavy elements (like gold and platinum) are created.

Summary in One Sentence

This paper successfully tested a new, mathematically stable way to simulate sticky neutron stars, proving that while the "pitch" of their vibrations stays the same, the "stickiness" makes them go silent much faster—a clue we can use to decode the secrets of the universe's most extreme matter.

Technical Summary

1. Problem Statement

Binary neutron star (NS) mergers are critical laboratories for testing fundamental physics, particularly the behavior of matter at extreme densities. While current numerical simulations of these events predominantly treat NS matter as an ideal fluid (neglecting viscosity), recent theoretical work suggests that viscous effects (dissipation) may play a significant role in the post-merger dynamics, potentially influencing gravitational wave signals and the restoration of beta equilibrium.

However, incorporating viscosity into relativistic hydrodynamics presents severe mathematical challenges:

• Ill-posedness: Classical first-order theories (Eckart, Landau-Lifshitz) are mathematically ill-posed, leading to acausal propagation and instability.
• Limitations of MIS: The widely used Israel-Stewart (MIS) second-order theories are often well-posed but can suffer from causality violations in realistic, highly dynamical simulations (e.g., early merger stages) where viscous terms are large.
• Need for New Frameworks: There is a pressing need for a first-order formulation that is causal, stable, and strongly hyperbolic (locally well-posed) to accurately model out-of-equilibrium effects in NS mergers.

The Bemfica-Disconzi-Noronha-Kovtun (BDNK) formulation was recently proposed as a solution, offering a mathematically sound first-order framework. This paper aims to test the viability of the BDNK framework for astrophysical applications by simulating neutron stars.

2. Methodology

Theoretical Framework (BDNK)

The authors utilize the BDNK formulation, which treats relativistic viscous hydrodynamics as an Effective Field Theory (EFT). Key features include:

• Field Redefinitions: It allows for arbitrary field redefinitions of thermodynamic variables up to first order in derivatives, defining a "hydrodynamic frame."
• Stress-Energy Tensor: The stress-energy tensor ($T^{\mu\nu}$) includes first-order derivative terms for energy density and velocity, parameterized by transport coefficients (shear viscosity $\eta$, bulk viscosity $\zeta$, and relaxation times $\tau_\epsilon, \tau_p, \tau_Q$).
• Well-Posedness: By imposing causality and stability constraints, the authors identify specific "open sets" of hydrodynamic frames where the initial value problem is well-posed. Crucially, this excludes the standard Landau and Eckart frames.

Numerical Setup

• Geometry: The simulations assume spherical symmetry and employ the Cowling approximation (fixed background spacetime, ignoring gravitational backreaction) to isolate hydrodynamic effects.
• Equation of State (EoS): A simplified, phenomenological EoS is constructed by combining a polytropic relation and an ideal gas law to ensure $p(\epsilon)$ is well-defined. This allows for a transition between cold and hot regimes.
• 3+1 Decomposition: The BDNK equations are decomposed into a system of balance laws. The authors introduce a first-order reduction in time, promoting time derivatives of primitive variables (energy density $\epsilon$ and velocity $v^i$) to dynamical fields ($\hat{\epsilon}, \hat{v}^i$).
• Primitive Recovery: A critical step involves mapping conservative variables (energy/momentum) back to primitive variables. Unlike ideal fluids (non-linear recovery), the BDNK formulation allows for an analytic linear recovery of the time-derivative variables, simplifying the numerical implementation.
• Numerical Scheme:
  • Time Integration: Third-order Strong Stability Preserving Runge-Kutta (SSP-RK).
  • Spatial Discretization: Third-order finite-volume scheme (equivalent to 4th-order finite difference with 3rd-order dissipation).
  • Stability: To handle the high characteristic speeds of the viscous system, the scheme uses the maximum characteristic velocity to determine numerical dissipation. A "staggered grid" is used at the origin ($r=0$) to avoid coordinate singularities.

Parameter Space Exploration

The authors systematically explore the BDNK parameter space by fixing the hydrodynamic frame parameters ($\hat{s}, \hat{a}, \hat{q}$) to satisfy well-posedness and causality conditions. They then vary the transport coefficients to create four distinct test cases:

1. smallSB-F2: Small shear and bulk viscosity.
2. medS-F2: Medium shear viscosity, zero bulk viscosity.
3. highB-F9: High bulk viscosity, low shear viscosity.
4. medSB-F9: Moderate shear and bulk viscosity.

3. Key Contributions

1. First Non-Linear Simulation: This is the first work to perform fully non-linear numerical simulations of neutron stars using the BDNK formulation. Previous works were limited to linear perturbations or ideal fluids.
2. Validation of Stability: The authors demonstrate that stable, long-term evolutions (up to $t \approx 8000 M_\odot$) are achievable within the BDNK framework for a wide range of viscosity parameters, provided the hydrodynamic frame is chosen correctly.
3. Primitive Recovery Implementation: They successfully implemented and validated the analytic primitive variable recovery scheme specific to the BDNK stress-energy tensor, proving its feasibility for complex viscous systems.
4. Separation of Numerical vs. Physical Viscosity: The study develops a robust methodology to disentangle physical viscous decay from numerical dissipation by performing convergence tests and extrapolating decay rates to the continuum limit.

4. Results

Stability and Evolution

• Stable Evolutions: All four parameter sets resulted in stable neutron star configurations. The stars remained in equilibrium, with only minor deviations near the center and surface attributed to numerical truncation errors and dissipation.
• EFT Regime: The simulations confirmed that the system remains within the validity regime of first-order hydrodynamics (dimensionless viscous terms $\ll 1$).

Quasi-Normal Modes (QNMs)

• Frequency Spectrum: The authors analyzed the power spectral density of central energy density oscillations.
  • Fundamental Mode (f-mode): The frequency of the f-mode ($\approx 2.69$ kHz) was found to be insensitive to the viscosity and the choice of hydrodynamic frame. It matched the Perfect Fluid (PF) prediction and literature values for the Cowling approximation.
  • Overtones: Higher-order modes (H1, H2) showed slight dependence on viscosity, likely because viscous terms affect short-length-scale physics more significantly.
• Decay Rates:
  • Viscous cases exhibited significantly faster decay rates compared to the Perfect Fluid (where decay is purely numerical).
  • Bulk vs. Shear: While bulk viscosity generally dominated the decay, the study found that shear viscosity also contributed non-negligibly to the decay rate even in spherical symmetry (a result that challenges the assumption that shear is irrelevant in spherical models).
  • Extrapolation: By extrapolating decay rates to zero grid spacing ($\Delta r \to 0$), the authors isolated the physical decay rate, confirming that the BDNK framework captures physical dissipation distinct from numerical artifacts.

5. Significance and Future Outlook

• Proof of Concept: This work establishes the BDNK formulation as a viable, robust alternative to MIS for modeling neutron star mergers, particularly in regimes where causality violations might plague other theories.
• Astrophysical Implications: The insensitivity of the f-mode frequency to viscosity suggests that gravitational wave frequency measurements alone may not easily constrain viscosity; however, the decay rate (damping) of the post-merger signal is a sensitive probe.
• Path Forward: The authors outline necessary steps for future research to make these models fully realistic:
  • Relaxing the Cowling approximation to include dynamical gravity (General Relativity).
  • Removing spherical symmetry to study non-radial modes and gravitational wave emission.
  • Implementing independent baryon number conservation to allow for more complex, realistic Equations of State.
  • Incorporating microphysics (e.g., weak interaction rates) to derive transport coefficients from first principles rather than phenomenological fitting.

In conclusion, this paper successfully bridges the gap between the mathematical rigor of the BDNK formulation and practical astrophysical simulations, paving the way for more accurate modeling of viscous effects in neutron star mergers and the interpretation of future gravitational wave data.