Examining the influence of anisotropy on the fundamental mode of nonradial oscillation in neutron stars on a complete general relativistic scheme

This paper investigates how local anisotropy significantly alters the fundamental oscillation frequency and dimensionless tidal deformability of neutron stars within a complete general relativistic framework, ultimately analyzing the correlation between these anisotropic effects and the GW170817 event.

The Gist

Imagine a neutron star as the ultimate cosmic weightlifter. It's a dead star, crushed so tightly by its own gravity that a single teaspoon of its material would weigh a billion tons on Earth. For decades, physicists have treated these cosmic giants as perfect, uniform spheres of fluid, like a giant, dense ball of water.

But this new paper asks a simple, yet profound question: What if the inside of these stars isn't perfectly uniform? What if the pressure pushing outwards is stronger in some directions than others? This is called anisotropy.

Think of a neutron star not as a smooth ball of water, but more like a giant, cosmic stress ball. If you squeeze a stress ball, it doesn't just squish evenly; it bulges out in specific directions. The authors of this paper wanted to see how this "squishing" (anisotropy) changes the way neutron stars behave, specifically how they "sing" and how they get squished by their neighbors.

Here is the breakdown of their findings using everyday analogies:

1. The "Singing" Star (Oscillations)

Neutron stars aren't silent. When they are disturbed—perhaps by a collision or a supernova explosion—they vibrate like a struck bell. The deepest, most fundamental note they can sing is called the "f-mode."

• The Old View: Scientists used to think the pitch of this "bell" depended only on how heavy and big the star was.
• The New Discovery: The authors found that the "stress" inside the star changes the pitch.
  • If the internal pressure is "stiffer" in certain directions (positive anisotropy), the star becomes slightly heavier and larger, and its "singing" note gets higher (faster frequency).
  • If the pressure is "looser" in certain directions (negative anisotropy), the star is smaller and lighter, and the note gets lower.
• The Analogy: Imagine a guitar string. If you tighten the string (increase tension/anisotropy), the note goes up. If you loosen it, the note goes down. The authors showed that the "tension" inside a neutron star is a hidden variable that changes the music the star makes.

2. The "Cosmic Squeeze" (Tidal Deformability)

When two neutron stars orbit each other, they pull on one another, stretching their shapes like taffy. This stretching ability is called tidal deformability.

• The Old View: We assumed this stretching depended only on the star's mass and size.
• The New Discovery: The internal "stress" of the star acts like a hidden spring.
  • If the star has high internal anisotropy (positive $\alpha$), it becomes stiffer. It resists being squished by its partner star. It's like a firm marshmallow vs. a soft one.
  • If the anisotropy is negative, the star is softer and squishes more easily.
• Why it matters: When the LIGO and Virgo detectors heard the "chirp" of two neutron stars colliding (the famous GW170817 event), they measured how much the stars squished. The authors show that by accounting for this internal "stress," their models fit the real-world data perfectly. It's like finally finding the right key to unlock the door to understanding what these stars are made of.

3. The "Recipe" (Equation of State)

To do this math, the authors had to create a recipe for the star's interior. They didn't just guess; they mixed two different types of physics:

1. Nuclear Physics: How atoms behave at normal (but still dense) levels.
2. Quantum Physics (QCD): How quarks behave at the extreme densities found in the core.

They stitched these two recipes together like a patchwork quilt. Then, they added the "anisotropy factor" (the stress) to see how the final dish tasted.

4. The Big Picture: Why Should We Care?

This paper is like finding a new dial on a radio.

• Before: We could tune the radio (the star's mass and size) to get a signal, but the sound was sometimes fuzzy.
• Now: We have a new dial (anisotropy). By turning this dial, we can get a much clearer signal that matches what our detectors actually hear.

The authors conclude that anisotropy is a crucial ingredient. If we ignore it, we might misunderstand the size, weight, and composition of these mysterious objects. By including this "internal stress," their models predict that these stars could be detected by our current telescopes if they are within our own galaxy, and they fit perfectly with the data we already have from the 2017 collision.

In short: Neutron stars aren't just uniform balls of dough. They are complex, stressed-out objects that sing different songs and squish differently depending on how their internal pressure is arranged. This paper gives us the sheet music to finally understand their song.

Technical Summary

1. Problem Statement

Neutron stars (NS) are extreme astrophysical objects where matter exists at densities far exceeding nuclear saturation. While standard models often assume the fluid inside these stars is isotropic (pressure is equal in all directions), theoretical evidence suggests that extreme conditions (e.g., strong magnetic fields, phase transitions, or superfluidity) can induce anisotropic pressure (where radial pressure $p_r$ differs from tangential pressure $p_t$).

The paper addresses the gap in understanding how this anisotropy affects two critical observables:

1. The f-mode (fundamental mode) oscillation frequency, which is a key signature in gravitational wave (GW) asteroseismology.
2. The dimensionless tidal deformability ($\Lambda$), a parameter constrained by GW events like GW170817.

The authors aim to derive the complete general relativistic (GR) equations for nonradial oscillations in anisotropic matter and quantify how anisotropy alters the star's equilibrium structure, oscillation frequencies, and detectability.

2. Methodology

Theoretical Framework

• General Relativity: The study employs the full Einstein field equations ($G_{\mu\nu} = 8\pi T_{\mu\nu}$) with a stress-energy tensor modified to include an anisotropic factor $\sigma$.
• Equation of State (EOS): The authors use a sophisticated EOS constructed by matching microscopic nuclear calculations (Chiral Effective Field Theory) at low densities with perturbative Quantum Chromodynamics (pQCD) at high densities. These are connected via a piecewise polytropic interpolation scheme (specifically "EOS III" from reference [58]), ensuring the model supports $\sim 2 M_\odot$ stars.
• Anisotropic Profile: A local anisotropy model is adopted:
  $$\sigma = \alpha p_r \left(1 - \frac{1}{g_{11}}\right)^n$$
  where $\alpha$ is a dimensionless anisotropy constant (ranging from $-1$ to $1$),$g_{11}$ is the metric component, and $n=2$ is chosen based on previous successful applications to quasinormal modes. The model ensures $\sigma=0$ at both the center and the surface.

Derivation of Equations

• Static Equilibrium: The Tolman-Oppenheimer-Volkoff (TOV) equations are modified to include the anisotropic term $\sigma$.
• Nonradial Perturbations: The authors derive the linearized perturbation equations for even-parity modes within the full GR framework.
  • They decompose the metric perturbation $h_{\mu\nu}$ using Regge-Wheeler gauge.
  • Crucially, they modify the first law of thermodynamics to account for the work done by anisotropic pressure varying with direction.
  • They derive a coupled system of differential equations for perturbation variables ($\tilde{H}, \tilde{H}_1, \tilde{K}, \tilde{W}, \tilde{V}$) and introduce an auxiliary function $\tilde{X}$ to close the system.
• Boundary Conditions:
  • Center ($r=0$): Regularity conditions are derived via Taylor series expansions to ensure finite solutions.
  • Surface ($r=R$): The Lagrangian pressure perturbation $\Delta p_r$ is set to zero.
  • Exterior: The equations are matched to the Zerilli equation to determine quasi-normal modes.

Numerical Implementation

• Oscillation Frequencies: A "shooting method" is used. Solutions are integrated from the center and surface to a midpoint, then matched to find the complex eigenfrequencies ($\omega = \omega_R + i\omega_I$).
• Tidal Deformability: The Riccati-type differential equation for the tidal Love number $y(r)$ is integrated alongside the structure equations.
• Detectability: Signal-to-noise ratios (SNR) are estimated for Advanced LIGO/Virgo and the Einstein Telescope (ET) based on the f-mode frequency and damping time.

3. Key Contributions

1. Derivation of Complete GR Equations: The paper provides a rigorous derivation of nonradial oscillation equations for anisotropic fluids, correcting inconsistencies found in previous literature (specifically differing from Refs [36, 37] regarding the thermodynamic work term).
2. Systematic Parameter Study: It systematically explores the impact of the anisotropy parameter $\alpha$ (from $-1.0$ to $+1.0$) on stellar structure and dynamics, rather than just assuming a fixed value.
3. Correlation with Observables: The study directly links theoretical anisotropy parameters to observable GW quantities (f-mode frequency and tidal deformability), offering a pathway to constrain anisotropy using future GW data.

4. Key Results

Equilibrium Configurations

• Mass and Radius: Anisotropy significantly alters the mass-radius relationship.
  • Positive $\alpha$ ($\alpha > 0$): Leads to higher maximum masses and slightly larger radii. The anisotropic pressure provides additional support against gravitational collapse.
  • Negative $\alpha$ ($\alpha < 0$): Results in lower maximum masses and smaller radii.
• Observational Consistency: The model configurations with various $\alpha$ values remain consistent with NICER observations of PSR J0030+0451 and PSR J0740+6620.

Oscillation Frequencies (f-mode)

• Frequency Shift: The f-mode frequency ($f_f$) is highly sensitive to anisotropy.
  • As $\alpha$ increases from $-1.0$ to $+1.0$, the f-mode frequency increases.
  • For a central density $\rho_c = 600 \text{ MeV/fm}^3$, changing $\alpha$ from $-1$ to $+1$ increases the f-mode frequency by approximately 23%.
• Scaling: The normalized frequency $\omega_f \sqrt{R^3/M}$ shows an almost linear relationship with mass up to the maximum mass limit, but the slope depends heavily on $\alpha$.

Tidal Deformability ($\Lambda$)

• Inverse Correlation: There is a clear inverse relationship between anisotropy and tidal deformability.
  • Positive $\alpha$: Reduces $\Lambda$ (stars become "stiffer" or more compact).
  • Negative $\alpha$: Increases $\Lambda$ (stars become "softer").
• GW170817 Constraints: All calculated curves for $\Lambda$ vs. Mass fall within the 90% confidence interval reported by LIGO-Virgo for the GW170817 event ($\Lambda_{1.4} = 190^{+390}_{-120}$). This suggests that anisotropic models are viable candidates for explaining the event.

Detectability

• Energy Requirements: The study calculates the minimum GW energy ($E_{gw}$) required to detect the f-mode signal with SNR > 5.
  • For a source at 10 kpc, Einstein Telescope (ET) can detect signals with energies as low as $\sim 10^{-10} M_\odot$, whereas Advanced LIGO/Virgo requires $\sim 10^{-7} M_\odot$.
  • Anisotropy affects the damping time ($\tau$) and frequency, thereby shifting the detection thresholds. Higher $\alpha$ generally leads to shorter damping times and higher frequencies.

5. Significance and Conclusion

This work demonstrates that anisotropy is a critical parameter in neutron star physics that cannot be ignored when interpreting gravitational wave data.

• Asteroseismology: The f-mode frequency serves as a potential "diagnostic" for the internal pressure anisotropy. Future high-precision GW detections could constrain the value of $\alpha$, distinguishing between different microphysical models of dense matter.
• Equation of State: The results show that anisotropy can mimic or mask the effects of the Equation of State. Ignoring anisotropy could lead to incorrect inferences about the stiffness of nuclear matter.
• Future Observations: The findings suggest that next-generation detectors (like the Einstein Telescope) will be capable of detecting f-mode oscillations from galactic sources, providing a new window to test the validity of anisotropic fluid models in General Relativity.

In summary, the paper establishes a robust theoretical framework showing that anisotropy significantly modifies the structural and oscillatory properties of neutron stars, with measurable consequences for current and future gravitational wave astronomy.