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

arXiv:2610.21966 · gr-qc, astro-ph.HE · Submitted 2026-10-15 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: We're looking at the paper: "Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations".

Jocelyn: The gist:

Vera: First, who's behind it and why it matters.

Paper summary: Vera: So, we're looking at this paper called "Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations." It's about how modeling viscous neutron stars using a certain theory of relativistic hydrodynamics—specifically the Bemfica– Disconzi–Noronha–Kovtun, or BDNK theory—can run into some problems.

Jocelyn: Basically, they argue that the idea of hydrodynamic universality breaks down in radial neutron star oscillations because things like causality and finite-size effects stop us from seeing that long-wavelength universal regime when the stars get sufficiently viscous.

Subrahmanyan: That means we can't always trust these simple hydrodynamic models for very thick, very viscous neutron stars. What they show is that this breakdown happens for a broad range of causal hydrodynamic frames in the BDNK theory when viscosity gets high enough.

Vera: They find that increasing the bulk viscosity pushes these causal frames away from what they call the BDNK regime of validity, eventually forcing all of them outside it when bulk viscosities are around zeta equal to one zero three zero or one zero three one g per centimeter per second.

Jocelyn: And this breakdown isn't just theoretical; it comes with real material changes in how the oscillation modes look. They find that non-hydrodynamic modes start getting damped on timescales that are comparable to the hydrodynamic modes, and the first overtone can actually become longer lived than the fundamental mode.

Subrahmanyan: That’s a significant point because it suggests that what we see in these oscillations isn't just simple fluid behavior anymore; you get these new effects appearing quickly. Plus, they mention that unphysical instabilities with sub-millisecond timescales can pop up.

Vera: Even though this universality breaks down, the paper points out something interesting about the gravitational collapse threshold. They find that the inviscid threshold of gravitational collapse remains unchanged in BDNK theory across all three viscosities and frames compared to the inviscid case.

Jocelyn: That's a bit like finding a stable anchor point even when everything else around it is shifting, which is important for understanding how these stars behave under extreme conditions.

Subrahmanyan: So, what this means for the bigger picture in cosmology and structure formation is that first-order theories of relativistic hydrodynamics might not be enough on their own when modeling these highly viscous neutron stars.

Vera: The authors suggest that second-order theory, like the one presented in reference eighty-five, would be better suited to study these cases because its first-order part is BDNK theory and it includes IS theory as a special case.

Jocelyn: To wrap up this paper on "Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations," the main idea is that causality and finite size effects restrict access to the long-wavelength universal hydrodynamic regime in sufficiently viscous stars, specifically when bulk viscosities are around one zero three zero to one zero three one g per centimeter per second.

Subrahmanyan: It really highlights that even with these specific first-order theories, we see these features—like the first overtone becoming the longest lived—even in a regime where universality has broken down.

Vera: The implication is that at least at this first order, we are seeing these complex behaviors in a region where universality has failed.

Jocelyn: We need to keep looking into whether this breakdown is just an artifact of the first-order theory, or if it points to something deeper about the physics of highly dissipative matter.

Subrahmanyan: And we need to see if those unphysical instabilities they mentioned are actually physically real, or if that's just a signal that we need a higher-order description.

Conclusion: Jocelyn: So we've been looking at how viscous neutron stars behave when you use these first-order hydrodynamics models, and this paper, "Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations," is really getting at the limits of those models.

Vera: Yeah, it tackles that idea of hydrodynamic universality—the general rule that says all these different simple equations should give you the same physical result in the long wavelength limit—and shows where it falls apart.

Subrahmanyan: It focuses on causal first-order theories, specifically using this Bemfica–Disconzi–Noronha–Kovtun or BDNK theory as a test case for how viscosity affects things.

Jocelyn: The authors are showing that causality and finite-size effects actually restrict access to that long wavelength universal regime when the stars get quite viscous.

Vera: And they pinpoint this breakdown happening at bulk viscosities in the range of zeta, which is about one zero three zero to one zero three one g per centimeter per second.

Subrahmanyan: What this means for structure formation is that if you're modeling a neutron star with that much viscosity, these simple first-order theories might just not be enough to get the full picture.

Jocelyn: They also noted some weird changes in the oscillation modes themselves, like non-hydrodynamic modes getting damped on timescales similar to the hydrodynamic ones.

Vera: And they mentioned unphysical instabilities popping up with timescales less than one millisecond, which is something you don't want in a realistic model.

Subrahmanyan: But here's a thing—even when universality breaks down, the paper says the threshold for gravitational collapse remains the same across different viscosities and frames.

Jocelyn: That’s interesting because it suggests some fundamental stability is still there, even if our simple models are failing to describe it correctly.

Vera: So, in short, this paper uses these specific first-order equations to show that causality itself limits how far we can go with simple hydrodynamic approximations in very viscous stars.

Subrahmanyan: And what's next is figuring out if this breakdown is just a limitation of the theory or if it points toward needing a more complex second-order theory for these extreme objects.

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

Wake Forest University

gr-qc, astro-ph.HE

Submitted: 2026-10-15

Updated: 2026-10-15

Comments: 10 pages, 6 figures

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 85/100

The gist: The gist: Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations shows that causality and finite-size effects restrict access to the long-wavelength universal

Key concepts

Hydrodynamic Universality: This refers to the idea that different mathematical descriptions of relativistic fluid dynamics yield the same physical results when looking at very long wavelengths. The paper shows this universality fails in radial neutron star oscillations because causality and finite size limitations restrict access to this simple regime in viscous stars.
BDNK Theory: This is a specific causal first-order theory used to model the fluid dynamics of the neutron star. It is a covariant generalization of Navier-Stokes equations that ensures non-linear causality and linear stability in flat space, making it suitable for studying relativistic fluid behavior.
Causal Frame Deviation ($\Delta$): This metric measures how far a specific causal hydrodynamic frame deviates from the ideal universal regime. If this deviation is large (e.g., $\ge 10\%$), it indicates that the system is outside the universal hydrodynamic regime, signaling a breakdown of universality.
Bulk Viscosity ($\zeta$): This parameter quantifies the internal friction or resistance to flow within the neutron star fluid. The study finds that increasing this bulk viscosity pushes causal frames away from the BDNK validity region, eventually forcing them outside it at high values, leading to mode structure changes.

Summary

The gist: Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations shows that causality and finite-size effects restrict access to the long-wavelength universal hydrodynamic regime in sufficiently viscous stars, occurring at bulk viscosities around 1030–1031 g/cm/s. It has been recently shown that hydrodynamic universality can break down in radial neutron star oscillations because causality and finite-size effects restrict access to the long-wavelength universal hydrodynamic regime in sufficiently viscous stars. This breakdown occurs for a broad range of causal first-order hydrodynamic frames, and increasing viscosity collectively pushes this class of causal frames away from the BDNK regime of validity, eventually forcing all to lie outside it at bulk viscosities in the range ζ ∼ 1030–1031g/cm/s. This breakdown is accompanied by material changes in the oscillation mode structure, including non-hydrodynamic modes becoming damped on timescales comparable to those of hydrodynamic modes, and unphysical instabilities with sub-millisecond timescales can emerge.

Theoretical Framework and Causality

The study utilizes the Bemfica–Disconzi–Noronha–Kovtun (BDNK) theory as a representative causal first-order theory, which is covariant generalization of the Navier–Stokes equations that is nonlinearly causal, linearly stable in flat space, locally well-posed, and strongly hyperbolic. The validity of first-order hydrodynamics requires the relaxation times not to be too large, for otherwise the solution depends on the causal regulators τi. Causality conditions require sufficient conditions for the causal propagation of linear perturbations of the Einstein–BDNK system, which are given by specific inequalities involving various transport coefficients and viscosities. Increasing bulk viscosity pushes the causal region away from the origin, at which acausal Navier–Stokes is located.

Mode Spectrum Analysis

The analysis in flat space reveals that imposing sufficient conditions for causal propagation leads to a minimum frame deviation ∆min scaling as ζ2. Defining the fractional frame deviation ∆ being at least 10% for all causal frames, i.e., ∆min = 0.1, the breakdown of hydrodynamic universality is diagnosed. For the background star considered here, Eq. (19) gives ζ∗ ∼ 3 × 1030 g/cm/s, which is close to the value of ζ∗ ∼ 2 × 1030 g/cm/s from Fig. 3 below. The effect of finite size manifests clearly in Eqs. (18–19), and if the system can access the universal regime k → 0, ∆ can be made arbitrarily small and ζ∗ arbitrarily large.

Neutron Star Oscillations and Instabilities

For the viscosities considered, it is possible that at the upper range there are order-unity corrections to the energy-momentum tensor which would necessarily place the theory outside the EFT regime. The deviation ∆GE scales linearly with ζ, and for the background star with ϵc = 4 × 1015 g/cm3 considered in this work, ranges from 0.0038 to 0.38 across ζc = 1029–1031 g/cm/s. The deviation ∆EFT assesses first-order frame robustness, a lack of which indicates that the system is not in the universal hydrodynamic regime. For ζc = 1031 g/cm/s, all the causal frames lie outside the EFT regime (i.e., ∆EFT ≥ 0.85). Furthermore, for ζc = 1029 g/cm/s, frames outside the EFT regime still give a fundamental mode close to that of acausal Navier–Stokes (i.e., ∆NS ≤ 0.016).

Conclusion on Model Sufficiency

The results provide evidence that first-order theories may be insufficient for hydrodynamic modeling of sufficiently viscous neutron stars. The breakdown of universality occurs for a broad class of causal first-order frames, providing evidence that first-order theories may be insufficient for hydrodynamic modeling of sufficiently viscous neutron stars. The central density at which the fundamental mode becomes unstable, i.e., the threshold of gravitational collapse, remains unchanged for all three viscosities and frames relative to the inviscid case. These results provide evidence that, at least at first order, such features occur in a regime in which universality has broken down. The second-order theory presented in Ref. [85] would be well suited for such a study, since its first-order truncation is BDNK theory and it includes IS theory as a special case.

Validation and Comparison

The frequency-domain methods employed allow us to more accurately investigate the mode structure near this threshold by solving directly for the eigenvalues and eigenvectors of the BDNK system, and employing physical boundary conditions at the surface which are able to preserve the location of the zero-frequency mode which, as argued above, is shared with the inviscid spectrum. The time-domain simulations reported a slight shift in the inviscid threshold of gravitational collapse within BDNK theory. The frequency-domain results cure the numerical viscosity-driven plateau at small viscosities. These results provide evidence that, although linearly stable in flat space [53, 59] and in zero-diffusion radial oscillations in a small-viscosity expansion [37], BDNK theory can admit unphysical instabilities, in this case occurring as hydrodynamic universality breaks down at sufficiently large viscosities. The breakdown of universality in sufficiently viscous, gravitationally stable neutron stars can lead to causal BDNK frames developing unphysical modes that grow on sub-millisecond timescales. These results provide evidence that, at least at first order, such features occur in a regime in which universality has broken down. The systematic second-order causal-frame analysis would determine whether universality is restored at these viscosities and such features are physical, or if they also occur outside the second-order EFT regime. The focus in this work on causal first-order oscillations is motivated by the fact that causality is an essential ingredient for the general stability properties of dissipative field equations [59]. This breakdown of universality occurs due to the fundamental ingredients of causality and finite size, extending beyond the specific case of radial neutron star oscillations considered here. The properties of the first-order mode spectrum—the first overtone becoming the longest lived and the emergence of dynamically important non-hydrodynamic modes—are consistent with the IS mode spectrum studied in Ref. [41]. The final sentence states that even in lower-viscosity regimes, a detailed comparison of time-domain evolutions and long-term dynamics between different first-order theories is necessary to establish whether a more general, physically-realizable universal hydrodynamic regime exists.

Causality and Frame Dependence

The analysis is confined to the causal region of parameter space where sufficient conditions for the causal propagation of linear perturbations are met. The causality conditions require specific relationships between relaxation times and viscosities. Increasing bulk viscosity pushes the causal region away from the origin, at which acausal Navier–Stokes is located. The analysis shows that increasing viscosity pushes these causal frames away from the origin and the EFT regime.

Breakdown Metrics

The breakdown of universality is monitored using several indicators. First, the fractional frame deviation ∆ being at least 10% for all causal frames is used as a diagnostic. Second, a proxy for the validity of first-order truncation is the condition ∆GE ≡ 1 − ωNS/ωPF ≪ 1, which scales linearly with ζ. Third, the deviation ∆EFT assesses first-order frame robustness, indicating a lack of universality.

Improvements for AI systems

This AI can perform systematic second-order causal-frame analysis to determine if hydrodynamic universality is restored at high viscosities, as suggested by The second-order theory presented in Ref. [85] would be well suited for such a study.
The improved system can differentiate between the effects of first-order frame transformations on mode frequencies by calculating the fractional deviations, specifically identifying when the change in the fundamental mode frequency upon transforming to a causal frame from acausal Navier–Stokes should be small relative to the change due to including viscosity (referencing Eq. 22).
The AI can predict and flag potential unphysical instabilities in numerical simulations by monitoring mode crossings and avoided crossings, where too large a step between successive viscosities in the walk can cause the tracked mode to jump to a distinct mode, which is indicated by checking if the overlap O(x, y) = ⟨x, y⟩/(∥x∥ ∥y∥) in ξ between neighboring viscosities in the walk is less than 0.9.
The system can provide a quantitative assessment of the breakdown of universality by calculating metrics such as ∆GE and ∆EFT, specifically identifying that for sufficiently large viscosity, all causal frames lie outside the EFT regime (i.e., ∆EFT ≥ 0.85) for stars with central bulk viscosities ζc = 1031 g/cm/s.
The AI can perform a sensitivity analysis on the threshold of gravitational collapse by confirming that the central density at which the fundamental mode becomes unstable, i.e., the threshold of gravitational collapse, remains unchanged across different causal frames and viscosities relative to the inviscid case.

Abstract

Hydrodynamic modeling of viscous neutron stars requires a choice of theory of relativistic hydrodynamics. That distinct formulations contain the same physical content in the long-wavelength limit is referred to as ``hydrodynamic universality.'' It has been recently shown that hydrodynamic universality can break down in radial neutron star oscillations: causality and finite-size effects restrict access to the long-wavelength universal hydrodynamic regime in sufficiently viscous stars. We show that this breakdown occurs for a broad range of causal hydrodynamic frames in Bemfica--Disconzi--Noronha--Kovtun (BDNK) theory. Increasing viscosity collectively pushes this class of causal frames away from the BDNK regime of validity, eventually forcing all to lie outside it at bulk viscosities in the range ζ about10 30 -- 10 31 g/cm/s. This breakdown of universality is accompanied by material changes in the oscillation mode structure: non-hydrodynamic modes become damped on timescales comparable to those of hydrodynamic modes, the first overtone can become longer lived than the fundamental mode, and unphysical instabilities with sub-millisecond timescales can emerge. Even in regimes where universality has broken down, we find that the inviscid threshold of gravitational collapse remains unchanged in BDNK theory. These results provide evidence that first-order theories of relativistic hydrodynamics may be insufficient for modeling sufficiently viscous neutron stars.