Axial Oscillations of Viscous Neutron Stars

Using a causal and stable theory of relativistic hydrodynamics, this paper investigates the axial oscillation modes of viscous neutron stars, revealing new viscosity-driven mode families and long-lived spectra that lack perfect fluid counterparts.

The Gist

Imagine a neutron star not as a solid, silent ball of dead matter, but as a giant, cosmic drum. When two of these stars collide or when they are disturbed, they don't just sit there; they vibrate, ringing like a bell. These vibrations, called oscillation modes, carry secrets about what the star is made of and how it behaves under extreme pressure.

For decades, physicists have studied these "cosmic drums" assuming the star is a perfect fluid—a theoretical substance that flows without any friction or resistance, like a ghostly liquid that never slows down.

But real neutron stars are messy. They are made of super-dense matter where particles crash into each other, creating viscosity (think of it as internal friction or "stickiness," like honey vs. water). This paper asks a simple but crucial question: What happens to the sound of the cosmic drum if the fluid inside is sticky?

Here is the breakdown of their discovery, using some everyday analogies:

1. The Problem: The "Perfect" vs. The "Real"

Until now, scientists mostly ignored the "stickiness" (viscosity) because the math for sticky fluids in Einstein's gravity is incredibly messy and often breaks down (it predicts things moving faster than light or becoming unstable).

The authors used a brand-new, robust mathematical framework (called BDNK Hydrodynamics) that acts like a "traffic cop" for these equations. It ensures the physics stays logical, causal (nothing travels faster than light), and stable, even when the fluid is very sticky.

2. The Discovery: Two Families of Vibrations

When they turned on the "stickiness" in their simulations, they found two distinct types of vibrations:

• The "W-Modes" (The Spacetime Ripples):

  • Analogy: Imagine hitting a drum. The skin vibrates, but the air around it also ripples. These are the W-modes. They are vibrations of the star's very fabric of space and time.
  • The Result: The authors found that if the star is very sticky, these ripples slow down slightly and lose energy a bit faster. However, the change is relatively small (a few percent). It's like hitting a drum with a slightly damp cloth; the sound changes, but you still recognize the note.

• The "Eta-Modes" (The New Discovery):

  • Analogy: This is the big surprise. Imagine a slinky. If you shake it, it bounces. But if you fill that slinky with thick, sticky syrup, a new way of moving appears that only exists because of the stickiness. The syrup itself becomes the spring that pushes the motion back and forth.
  • The Result: The authors discovered a completely new family of vibrations called $\eta$-modes (Eta-modes). These modes do not exist in a perfect fluid. They are born entirely from the viscosity (the "stickiness") of the star's interior. They oscillate at high frequencies (thousands of times a second) and can last for milliseconds.

3. The "Dance of Avoidance"

Here is where it gets really cool. In physics, when two different types of vibrations get close to each other in frequency, they usually don't cross paths. Instead, they "repel" each other, like two magnets with the same pole facing each other.

• The Analogy: Imagine two dancers on a floor. As they move toward each other, instead of colliding, they suddenly step aside, changing their rhythm to avoid bumping into one another.
• The Result: As the star gets stickier, the new Eta-modes try to cross paths with the old W-modes. Instead of crossing, they "avoid" each other. This causes the W-modes to suddenly shift their frequency and damping rate in a chaotic way. This "mode avoidance" is a fingerprint that tells us the star is viscous.

4. Why Does This Matter?

• The Cosmic Detective: Future gravitational wave detectors (like the Einstein Telescope) will be able to "hear" the vibrations of neutron stars after they collide. If we hear these specific "Eta-modes" or see the "avoidance" dance, we can prove that the matter inside the star is incredibly sticky.
• The "Black Hole Mimicker" Test: Some objects look like black holes but might be ultra-dense stars. If they are sticky, their vibrations die out much faster than a perfect black hole would. This helps us distinguish between a real black hole and a weird, sticky star.
• The Limits of Reality: The authors had to use "super-sticky" values (stickier than we think is physically possible for normal matter) to see these effects clearly. However, if the core of a neutron star contains "strange matter" (exotic particles), it might be sticky enough to make these modes real.

Summary

Think of this paper as tuning a cosmic instrument.

• Before: We thought the instrument was made of frictionless glass.
• Now: We realized it's made of thick, sticky honey.
• The Change: The honey doesn't just muffle the sound; it creates new notes (Eta-modes) that we never knew existed and makes the old notes dance around each other (Mode Avoidance).

By listening for these new notes and dances, astronomers in the future might finally be able to taste the "texture" of the densest matter in the universe.

Technical Summary

1. Problem Statement

Neutron stars (NSs) are unique astrophysical laboratories where gravity, electromagnetism, and the weak and strong nuclear forces all play critical roles. While gravitational wave (GW) observations (e.g., from LIGO/Virgo and future third-generation detectors like Einstein Telescope) aim to probe the equation of state (EoS) of dense matter, current theoretical models often neglect dissipative effects such as shear and bulk viscosity.

• The Gap: Standard studies assume NSs are in chemical equilibrium or treat them as perfect fluids. However, violent events like binary mergers strongly perturb this equilibrium, making dissipative processes (mediated by reactions like modified Urca or hyperon interactions) non-negligible.
• The Challenge: Previous attempts to model viscosity in relativistic hydrodynamics often relied on the Israel-Stewart formalism (second-order theory), which introduces complex higher-order gradients and potential stability issues. There has been no systematic framework to study how viscosity affects non-radial oscillation modes (specifically axial modes) using a causal, stable, first-order theory.

2. Methodology

The authors employ a rigorous framework combining General Relativity with a modern formulation of relativistic hydrodynamics.

A. Theoretical Framework: BDNK Hydrodynamics

Instead of the traditional Israel-Stewart approach, the authors utilize BDNK Hydrodynamics (Bemfica-Disconzi-Noronha-Kovtun).

• Key Feature: This is a causal and stable first-order theory. It achieves well-posedness by allowing for frame transformations of the hydrodynamic variables (density $\rho$ and 4-velocity $u^\mu$) rather than introducing second-order gradient terms.
• Stress-Energy Tensor: The tensor includes first-order corrections involving shear viscosity ($\eta$) and bulk viscosity ($\zeta$). The authors set $\zeta=0$ for axial perturbations as it does not couple to this sector.
• Parametrizations: Two distinct parametrizations (labeled A and B) for the transport coefficients are tested to ensure results are independent of the specific choice of hydrodynamic frame. These involve dimensionless constants $\hat{\eta}$ and $\hat{\tau}$ and a length scale $L_0$.

B. Perturbation Analysis

• Background: Spherically symmetric equilibrium stars described by the Tolman-Oppenheimer-Volkoff (TOV) equations with polytropic and constant-density equations of state.
• Perturbations: The study focuses on axial (odd-parity) perturbations. The metric and fluid velocity are perturbed, leading to a coupled system of two Ordinary Differential Equations (ODEs) for the metric perturbation function $\psi$ and the fluid velocity perturbation $Z$.
• Governing Equations:
  1. A modified Regge-Wheeler equation for $\psi$ (spacetime modes).
  2. A second equation for $Z$ involving a "viscous propagation speed" (second sound) $c_\eta^2 = \eta / [\tau(p+\rho)]$.
• Boundary Conditions:
  • Center: Regularity ($\psi, Z \sim r^{\ell+1}$).
  • Surface: A specific regularity condition is derived to handle divergences caused by vanishing viscosity at the surface.
  • Infinity: Outgoing gravitational waves.

C. Numerical Implementation

• Quasinormal Modes (QNMs): The authors solve for complex frequencies $\omega = 2\pi f - i/\tau$ where solutions are regular at the origin and outgoing at infinity.
• Technique: A "shooting method" combined with a continued fraction expansion (Leaver's method) is used.
  • In-solution: Integrated from the center to a matching radius $a$.
  • Up-solution: Integrated from infinity (using a continued fraction ansatz) to $a$.
  • Matching: The Wronskian of the two solutions is minimized to find the roots (QNMs).

3. Key Contributions

1. First Systematic Study of Axial Viscous Modes: This is the first work to systematically characterize how shear viscosity affects the axial oscillation spectrum of neutron stars using a causal first-order relativistic hydrodynamic framework.
2. Discovery of $\eta$-modes: The authors identify a new family of oscillation modes driven purely by shear viscosity, dubbed $\eta$-modes. These modes have no counterpart in perfect fluid stars; their restoring force is provided by the shear viscosity and the regulator parameter $\tau$.
3. Mode Avoidance Phenomenon: The study demonstrates that as shear viscosity increases, the traditional spacetime w-modes and the new $\eta$-modes approach each other in the complex frequency plane but undergo avoided crossing (they repel rather than cross).
4. Frame Robustness: The authors prove that the frequencies of the standard w-modes are largely independent of the specific hydrodynamic frame parametrization (differences $<1\%$), validating the robustness of the BDNK approach for gravitational modes.

4. Key Results

• Impact on w-modes (Spacetime Modes):
  • For realistic to high viscosities ($\eta_c \sim 10^{29} - 10^{31}$ g/cm/s), the w-mode frequencies shift by up to 10% and damping times change significantly.
  • The shift is more pronounced in less compact stars.
  • The dependence on viscosity is approximately linear.
• Characterization of $\eta$-modes:
  • These modes oscillate at kHz frequencies (comparable to w-modes) and have damping times in the order of milliseconds.
  • They are highly sensitive to the choice of hydrodynamic frame and transport coefficients.
  • In the limit of zero viscosity, these modes become undamped ($\text{Im}(\omega) \to 0$).
• Mode Avoidance:
  • As viscosity increases, the $\eta$-modes approach the w-modes. Instead of crossing, they repel, causing a sudden change in the frequency and damping time of the w-modes. This suggests that in specific parameter regimes, the observed spectrum could be highly sensitive to the microscopic transport coefficients.
• Ultracompact Objects:
  • For ultracompact stars (with stable light rings), viscosity drastically alters the imaginary part of the frequency (damping time).
  • While perfect fluid ultracompact stars support long-lived modes (trapped by the light ring), shear viscosity damps these modes efficiently, potentially quenching nonlinear instabilities associated with stable light rings.

5. Significance and Implications

• Gravitational Wave Astronomy: The results imply that future high-precision GW asteroseismology (e.g., from post-merger signals) must account for viscous effects to accurately infer the Equation of State. Ignoring viscosity could lead to misinterpretation of the oscillation spectrum.
• New Physics Probes: The existence of $\eta$-modes and the phenomenon of mode avoidance offer a potential new channel to probe the viscous transport coefficients of dense matter, which are otherwise difficult to constrain.
• Theoretical Validation: The work validates the BDNK framework as a viable and computationally efficient tool for studying dissipative relativistic stars, offering a simpler alternative to second-order theories while maintaining causality and stability.
• Future Directions: The authors note that this study is limited to the axial sector. The polar sector (containing f-, p-, and g-modes) and rotating stars (r-modes) are expected to exhibit even richer structures involving dissipation, which are left for future investigation.

In summary, this paper establishes that viscosity is not merely a damping mechanism but a fundamental driver of new oscillation modes in neutron stars, fundamentally altering the interpretation of their dynamical spectrum in the era of multi-messenger astronomy.