⚛️ general relativity

Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations

This paper demonstrates that hydrodynamic universality breaks down in causal first-order viscous neutron star oscillations across a broad range of BDNK frames at high bulk viscosities, leading to significant alterations in oscillation mode structures and suggesting that such first-order theories are insufficient for modeling sufficiently viscous neutron stars.

Original authors: Lennox S. Keeble, Alejandro Cárdenas-Avendaño

Published 2026-10-16

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Lennox S. Keeble, Alejandro Cárdenas-Avendaño

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

Neutron stars are the dense, dead cores of massive stars that have exploded, packing more mass than our Sun into a sphere only about twenty kilometers wide. Inside these cosmic remnants, matter exists in states so extreme that it cannot be replicated in any laboratory on Earth. To understand how these stars behave, especially when they vibrate or collide, scientists rely on a set of rules called hydrodynamics. This is the study of how fluids flow, but applied to the relativistic speeds and crushing gravity found in space. For decades, physicists have used different mathematical frameworks to describe this flow, hoping that while the equations might look different, they would all predict the same physical reality when the fluid is moving slowly or over long distances. This expectation is known as "hydrodynamic universality." It is a comforting idea that nature has a single, consistent voice, regardless of which mathematical tool we use to listen to it. However, when these stars are subjected to extreme internal friction, or viscosity, that consistency might not hold.

A new study by researchers at Wake Forest University investigates what happens when a neutron star becomes so viscous that these different mathematical tools begin to disagree. The team focused on the radial oscillations of these stars, which are essentially the rhythmic pulsing of the star's entire volume, expanding and contracting like a beating heart. They used a modern, stable theory of fluid flow that respects the cosmic speed limit of light, known as the BDNK theory, to model these vibrations. By simulating stars with bulk viscosities ranging from $10^{29}$ to $10^{31}$ grams per centimeter per second, they discovered that the long-held belief in universality breaks down under these extreme conditions. When the internal friction reaches these high levels, the different ways of describing the fluid no longer agree on how the star should vibrate. The mathematical frameworks that were once interchangeable begin to produce different results, suggesting that the standard first-order theories used to model these stars are insufficient for describing matter with such high viscosity.

The researchers found that as the viscosity increases, the mathematical "window" where all theories agree shrinks and eventually disappears. In the simulations, this breakdown occurs when the bulk viscosity reaches values between $10^{30}$ and $10^{31}$ grams per centimeter per second. At this point, the different causal frames—different ways of defining the fluid's motion that all obey the laws of causality—stop predicting the same oscillation frequencies. Instead of a single, clear rhythm, the models begin to diverge. Some frames predict that the star's vibrations will die out quickly, while others suggest they will last much longer. In some cases, the models predict that the star could develop unphysical instabilities, where certain vibration modes grow wildly on timescales as short as a fraction of a millisecond. These are not real physical behaviors but rather artifacts indicating that the mathematical model has been pushed beyond its limits of validity.

One of the most striking findings is that this breakdown happens even in stars that are stable and not on the verge of collapsing. The study shows that the point at which a neutron star becomes unstable and collapses into a black hole remains the same regardless of the viscosity or the mathematical frame used. However, the journey to that point changes dramatically. As the viscosity rises, the usual order of the star's vibrations gets scrambled. The first overtone, which is a higher-pitched version of the main pulsation, can become longer-lived than the fundamental mode, a reversal of what is typically expected. Furthermore, modes of vibration that are usually invisible or quickly damped out become just as important as the main pulsations, complicating the picture of how the star moves.

The authors emphasize that these results are derived from numerical simulations of the equations governing the star's behavior. They do not claim to have observed this breakdown in a real neutron star, but rather that the current mathematical tools used to describe them may fail under such extreme conditions. The study suggests that if a neutron star does possess viscosity in this high range, which is within the bounds estimated from gravitational wave observations, then the first-order theories currently in use are not enough. To accurately model such an object, scientists may need to turn to more complex, second-order theories that can handle the far-from-equilibrium state of the fluid. Until then, the assumption that all causal hydrodynamic theories are equivalent in the long-wavelength limit cannot be trusted for sufficiently viscous neutron stars. This work serves as a crucial warning that our understanding of the most extreme matter in the universe may require a deeper, more robust mathematical foundation than previously thought.

Technical Summary: Breakdown of Hydrodynamic Universality in Causal First-Order Viscous Neutron Star Oscillations

Problem Statement
Hydrodynamic modeling of viscous neutron stars requires selecting a specific theory of relativistic hydrodynamics. While distinct formulations (e.g., different hydrodynamic frames) are expected to yield identical physical predictions in the long-wavelength, near-equilibrium limit—a property termed "hydrodynamic universality"—recent studies suggest this universality may fail in sufficiently viscous neutron stars. Specifically, causality constraints and finite-size effects can prevent the system from accessing the universal regime. Previous work by Denicol et al. demonstrated this breakdown in radial oscillations using Israel–Stewart (IS) theory versus acausal Navier–Stokes. However, it remained unclear whether this phenomenon is specific to second-order theories or extends to the broad class of causal first-order theories, such as the Bemfica–Disconzi–Noronha–Kovtun (BDNK) framework. This paper investigates whether hydrodynamic universality breaks down for a wide range of causal first-order frames in BDNK theory under realistic neutron star conditions.

Methodology
The authors employ a perturbative framework to study radial oscillations of viscous neutron stars within the BDNK theory.

Background and Constitutive Relations: The study utilizes a spherically symmetric spacetime coupled to a barotropic BDNK fluid with zero chemical potential. The constitutive relations include bulk viscosity ( $\zeta$ ), shear viscosity ( $\eta$ ), and relaxation times ( $\tau_{e,p,q}$ ). The background star is modeled using the Tolman–Oppenheimer–Volkoff (TOV) equations with a polytropic equation of state.
Numerical Approach: The radial perturbation equations are solved in the frequency domain. The authors discretize the system using second-order-accurate centered finite differences and a spatial Crank–Nicolson approach for the first-order Einstein equation to avoid unphysical grid modes. This results in a quadratic eigenvalue problem $Q(\mu)U = 0$ , where $\mu = \omega/\omega_{PF}$ is the frequency normalized by the perfect fluid eigenvalue.
Solvers: Two methods are used to solve the eigensystem: a Newton iteration (for tracking modes across viscosity) and the QZ algorithm (for non-hydrodynamic modes).
Causality and Parameter Space: The analysis focuses on the "causal region" of the parameter space defined by sufficient conditions for causal propagation. The authors parameterize the relaxation times and viscosities to define "rays" in the three-dimensional parameter space $(\tau_e, \tau_p, \tau_q)$ . They specifically analyze the "minimally causal ray," which represents the frames closest to the origin (where $\tau \to 0$ ).
Diagnostic Metrics: To quantify the breakdown of universality, the authors define fractional deviations:
- $\Delta_{EFT} = |(\omega_{BDNK} - \omega_{NS}) / (\omega_{NS} - \omega_{PF})|$ : Measures the deviation of the causal BDNK frequency from the acausal Navier–Stokes (NS) frequency relative to the viscous shift. $\Delta_{EFT} \ll 1$ indicates the system is in the universal EFT regime.
- $\Delta_{NS} = |(\omega_{BDNK} - \omega_{NS}) / \omega_{NS}|$ : Measures the shift from NS to BDNK.
- The breakdown of universality is diagnosed when $\Delta_{EFT} \geq 0.1$ for all causal frames.

Key Contributions and Results

Breakdown of Universality in First-Order Theory: The authors demonstrate that hydrodynamic universality breaks down for a broad class of causal first-order frames in BDNK theory. As bulk viscosity increases, the causal region of parameter space is pushed away from the origin (where the EFT regime resides) toward larger relaxation times.
Critical Viscosity Threshold: The breakdown occurs at bulk viscosities $\zeta \sim 10^{30}–10^{31} \text{ g/cm/s}$ . At these values, all causal frames lie outside the first-order EFT regime ( $\Delta_{EFT} \gtrsim 0.85$ for $\zeta_c = 10^{31} \text{ g/cm/s}$ ). This range aligns with previous estimates for neutron star bulk viscosity and the upper bound derived from GW170817.
Material Changes in Oscillation Modes: The breakdown of universality is accompanied by significant changes in the oscillation mode structure:
- Non-hydrodynamic modes become damped on timescales comparable to hydrodynamic modes.
- The first overtone can become longer-lived than the fundamental mode.
- Unphysical Instabilities: For sufficiently high viscosities ( $\zeta_c = 10^{31} \text{ g/cm/s}$ ), unphysical instabilities with sub-millisecond timescales emerge in the causal BDNK frames, even though the theory is linearly stable in flat space.
Invariance of Gravitational Collapse Threshold: Despite the breakdown of universality and the emergence of unphysical modes, the authors find that the inviscid threshold of gravitational collapse remains unchanged in BDNK theory. The zero-frequency mode (indicating the onset of instability) persists at the same central density as in the perfect fluid case, provided the boundary conditions are handled correctly in the frequency domain.
Frame Dependence: The study highlights that while the EFT regime is lost, the physical predictions (eigenvalues) can still vary significantly depending on the choice of causal frame, particularly as viscosity increases.

Significance and Claims
The paper provides evidence that first-order theories of relativistic hydrodynamics, including the causal BDNK framework, may be insufficient for modeling sufficiently viscous neutron stars. The breakdown of hydrodynamic universality implies that different causal formulations will no longer agree on physical predictions, and the theory may enter a regime where it admits unphysical instabilities.

The authors emphasize that while BDNK theory is designed to be causal and stable, these properties are guaranteed only within its regime of validity (the EFT regime). When viscosity is high enough to force the system outside this regime, the theory's predictive power regarding mode structure and stability is compromised. The results suggest that for astrophysical scenarios involving high bulk viscosity (potentially relevant to binary neutron star inspirals and post-merger phases), second-order theories or far-from-equilibrium prescriptions may be necessary to capture the correct physics. The authors explicitly state that their findings do not confirm that these unphysical instabilities will manifest in nature but indicate a limitation in the applicability of first-order causal hydrodynamics for such extreme conditions.

Technical Summary: Breakdown of Hydrodynamic Universality in Causal First-Order Viscous Neutron Star Oscillations

More like this