Radial Oscillations of Viscous Stars

TL;DR

This paper investigates how viscosity affects the radial oscillations of neutron stars, revealing damping effects, frequency shifts, and implications for gravitational collapse, which are crucial for gravitational wave astronomy.

Contribution

It provides a comparative analysis of viscous effects using two hydrodynamic frameworks, advancing understanding of neutron star oscillations and stability under viscosity.

Findings

Viscosity damps radial modes on millisecond timescales.

Frequency shifts increase with star compactness and viscosity.

Viscosity cannot prevent gravitational collapse in the studied models.

Abstract

Oscillation modes of neutron stars, a key target for third-generation gravitational wave detectors, encode key information about their constituent nuclear matter. In this work, we study the effect of viscosity on oscillations of cold, polytropic, spherically symmetric neutron stars. We focus on purely radial oscillations and work perturbatively to linear order within two hydrodynamic frameworks: the acausal covariant generalization of the Navier-Stokes equations proposed by Eckart, and the causal generalization formulated by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). We find that viscosity damps the radial modes on millisecond timescales and induces fractional shifts in the oscillation frequency which increase both with the compactness and viscosity of the star, reaching up to the percent level for the fundamental mode with bulk viscosities $\zeta\sim10^{30}\mathrm{g}/\mathrm{cm}/\mathrm{s}$. For more viscous stars, the oscillation frequency decreases, becoming zero (i.e., an overdamped mode) for $\zeta\gtrsim10^{31}\mathrm{g}/\mathrm{cm}/\mathrm{s}$. We also study the linear threshold of gravitational collapse. Consistent with recent analytic results in the small viscosity regime, we find that viscosity in Eckart theory cannot stabilize an unstable inviscid star. We provide numerical evidence that viscosity in BDNK theory is similarly unable to prevent gravitational collapse, but it slightly modifies the threshold of collapse. Overall, our results advance our understanding of the impact of viscosity on the oscillation modes of neutron stars, a key component of viscous asteroseismology with next-generation gravitational wave detectors.