General relativistic study of $f$-mode oscillations in neutron stars with gravitationally bound dark matter

This paper presents a comprehensive general relativistic study of how gravitationally bound dark matter admixed with neutron stars modifies fundamental $f$-mode oscillation frequencies and damping times, establishing new asteroseismic universal relations and deriving multimessenger constraints from the GW170817 event.

The Gist

Imagine a neutron star as a cosmic super-ball, the densest object in the universe, squeezed so tight that a teaspoon of its material would weigh a billion tons on Earth. Now, imagine this super-ball isn't just made of normal matter (like protons and neutrons) but has a secret ingredient hidden inside: Dark Matter.

This paper is a detailed investigation into what happens to these "super-balls" when they are mixed with this invisible, mysterious dark stuff. The author, Pinku Routaray, uses the heavy math of Einstein's General Relativity to simulate how these stars behave, specifically looking at how they "vibrate" or "ring" like a bell after being hit.

Here is the breakdown of the study using simple analogies:

1. The Setup: The "Two-Fluid" Smoothie

Usually, scientists think of a neutron star as a single, uniform block of matter. However, this study treats the star like a smoothie made of two distinct ingredients:

• The Fruit (Normal Matter): The heavy, visible stuff (protons and neutrons).
• The Ice (Dark Matter): The invisible stuff that doesn't interact with light but has gravity.

The author uses a specific recipe called the "Higgs-portal" model. Think of this as a special blender setting that allows the "Ice" (Dark Matter) to interact with the "Fruit" (Normal Matter) just enough to stick together, but not so much that they turn into a new substance. They stay as two separate fluids swirling inside the same container.

2. The Secret Recipe: Where the Ice Sits

A key finding of this paper is that the Dark Matter doesn't spread out evenly like sugar in tea. Because the star's gravity is so incredibly strong, the Dark Matter gets sucked to the very center, forming a dense core, while the outer layers are mostly normal matter.

The author uses two "knobs" to control this mixture:

• Knob A (Amount): How much Dark Matter is in the mix.
• Knob B (Steepness): How sharply the Dark Matter piles up in the center versus the edges.

The Analogy: Imagine a crowd of people in a stadium. If the "steepness" knob is high, the people (Dark Matter) are all huddled tightly in the very center seats, leaving the outer seats empty. If the knob is low, they are spread out more evenly.

3. The Experiment: Ringing the Bell

The main goal of the study was to see how this hidden "Ice" changes the way the star vibrates. When a neutron star is disturbed (perhaps by a collision), it vibrates in a specific way called the f-mode (fundamental mode). Think of this like striking a bell:

• The Pitch (Frequency): How high or low the sound is.
• The Damping (How long it rings): How quickly the sound fades away.

What the study found:

• Higher Pitch: Adding Dark Matter makes the star "tighten up" (become more compact). Just like a tighter drum skin makes a higher sound, the star vibrates at a higher frequency.
• Faster Silence: The presence of Dark Matter also makes the vibrations fade away faster. The energy of the vibration leaks out as gravitational waves (ripples in space-time) more quickly.

4. The Universal Rules (The "Laws of Physics")

Scientists have discovered "Universal Relations" (URs). These are like rules of thumb that say, "If you know how heavy and compact a star is, you can predict exactly how it will ring, regardless of what it's made of."

The big question was: Does adding Dark Matter break these rules?

• The Result: No! The study found that even with the secret Dark Matter ingredient, the "Universal Relations" still hold true. The star still follows the same predictable patterns. This is great news because it means astronomers can still use these rules to figure out what stars are made of, even if they contain Dark Matter.

5. The "Fingerprint" of Dark Matter

The paper also looked at real-world data from a famous event called GW170817 (a collision of two neutron stars detected by gravitational wave detectors).

• The author used this data to set limits on how much Dark Matter could be hiding inside a typical neutron star.
• They found that if there is too much Dark Matter, or if it is too concentrated in the center, the star would become so small and heavy that it wouldn't match what we actually see in the sky.
• Conclusion: There is a "Goldilocks zone" for Dark Matter in neutron stars. It can be there, but not in huge amounts, or the star would look different than the ones we observe.

6. Can We Hear It?

Finally, the paper asks: "If a star with Dark Matter vibrates, can our current detectors hear it?"

• The Verdict: For stars very close to us (inside our galaxy), the vibrations might be loud enough for future, super-sensitive detectors (like the Einstein Telescope) to hear.
• However, for stars far away (in other galaxy clusters), the signal is too weak for our current equipment. The Dark Matter makes the star vibrate faster and quieter, which actually makes it harder to detect from far away, but potentially easier to spot if we have very sensitive ears nearby.

Summary

This paper is a theoretical simulation showing that if neutron stars are hiding Dark Matter in their cores, they will vibrate at a higher pitch and fade faster than normal stars. However, they still follow the same universal laws of physics. By comparing these predictions with real data from gravitational wave detectors, we can figure out how much Dark Matter is allowed to hide inside these cosmic giants without breaking the rules of the universe.

Technical Summary

Technical Summary: General Relativistic Study of f-Mode Oscillations in Neutron Stars with Gravitationally Bound Dark Matter

Problem Statement
Neutron stars (NSs) serve as natural laboratories for dense matter physics, yet their internal composition remains uncertain due to the potential presence of exotic components, including Dark Matter (DM). While previous studies have modeled DM-admixed NSs, many relied on simplified assumptions, such as constant DM density or the Cowling approximation (ignoring metric perturbations). A physically realistic model requires accounting for the non-uniform distribution of DM caused by gravitational capture within the star's deep potential well. Furthermore, understanding how such DM distributions alter the non-radial oscillation modes (specifically the fundamental f-mode) and their gravitational wave (GW) damping times is crucial for interpreting multimessenger observations. This study addresses the gap in modeling gravitationally bound, non-uniform DM distributions and their impact on the full general relativistic (GR) oscillation spectrum of NSs.

Methodology
The study employs a comprehensive framework within full general relativity to model NSs admixed with fermionic DM via a Higgs-portal interaction.

1. Equation of State (EOS) and DM Distribution:

   • Hadronic Matter: Modeled using the Relativistic Mean Field (RMF) formalism with the Hornick3 model, extended with a BPS crust.
   • Dark Matter: Modeled as a non-annihilating Weakly Interacting Massive Particle (WIMP) interacting with baryons via the Higgs portal. The system is treated as a single-fluid approximation, justified by strong coupling and thermalization timescales shorter than oscillation periods.
   • Density Profile: Unlike previous constant-density models, the DM density ($n_{DM}$) is defined as a function of the baryon density ($n_B$) using a phenomenological profile: $n_{DM}/n_0 = \alpha ((n_B - n_t)/n_0)^\beta$.
     • $\alpha M_\chi$: An effective control parameter combining the DM concentration scaling factor ($\alpha$) and the DM particle mass ($M_\chi$).
     • $\beta$: A steepness parameter controlling how sharply the DM density peaks toward the core.
   • Constraints: The model ensures beta equilibrium and charge neutrality. The scattering cross-section is constrained by direct detection experiments (XENON-1T, PandaX-II, LUX) and LHC limits.

2. Hydrostatic Equilibrium:

   • The Tolman-Oppenheimer-Volkoff (TOV) equations are solved to determine the mass-radius profiles and the dimensionless tidal deformability ($\Lambda$) for various $(\alpha M_\chi, \beta)$ configurations.

3. Oscillation Analysis:

   • Full GR Perturbation: The study solves the linearized Einstein equations for both polar (even-parity) and axial modes, treating fluid and metric perturbations consistently (avoiding the Cowling approximation).
   • Boundary Conditions: Regularity at the center, vanishing pressure perturbation at the surface, and purely outgoing gravitational waves at infinity (solved via the Zerilli equation).
   • Numerical Solution: The complex quasi-normal mode (QNM) frequencies ($\omega = 2\pi f + i/\tau$) are computed using an adaptive Runge-Kutta method to find the roots of the scattering amplitude function.

Key Contributions and Results

• Impact of Non-Uniform DM on Structure:

  • The inclusion of non-uniform DM softens the effective EOS, leading to a reduction in the maximum stable mass and the radius of the NS.
  • Higher values of $\alpha M_\chi$ and $\beta$ result in more compact configurations. For steep profiles ($\beta=4$), the mass-radius relation deviates significantly from observational constraints (e.g., PSR J0740+0432, GW170817) at high DM concentrations.
  • The maximum allowed DM mass fraction ($f_{\chi, \max}$) is found to be sensitive to $\beta$; while it increases with $\alpha M_\chi$ for moderate $\beta$, it decreases for very steep profiles ($\beta=4$), suggesting that excessively concentrated DM suppresses the total allowable DM fraction to maintain stability.

• Modification of f-Mode Oscillations:

  • Frequency ($f$): The presence of DM, particularly with higher $\alpha M_\chi$ and $\beta$, shifts the f-mode frequency to higher values compared to pure hadronic stars of the same mass. This is attributed to the increased central concentration and stronger gravitational potential.
  • Damping Time ($\tau$): Conversely, the GW damping time decreases in the presence of DM. The study finds a stronger coupling between matter and spacetime perturbations, leading to enhanced energy dissipation.
  • Correlations: A strong correlation exists between the DM parameters ($\alpha M_\chi$) and the DM mass fraction. The f-mode frequency and $\tau$ show strong correlations with stellar compactness and tidal deformability, even in the presence of DM.

• Universal Relations (URs):

  • The study investigates whether the presence of non-uniform DM breaks the established "universal relations" between oscillation properties and global NS parameters.
  • $f-C-\tau$ and $f-\Lambda-\tau$: Analytic fits are constructed for the relations between f-mode frequency, damping time, compactness ($C$), and tidal deformability ($\Lambda$).
  • Robustness: The results indicate that the inclusion of gravitationally bound DM does not break these universal relations. The relations remain robust, allowing for the construction of fits that include both nucleonic and DM degrees of freedom.
  • Constraints: By mapping the tidal deformability constraint from GW170817 ($\Lambda_{1.4} = 190^{+390}_{-120}$) into the $(f, \tau)$ space, the study provides observational bounds for a canonical DM-admixed NS: $f_{1.4} = 2.110^{+0.445}_{-0.445}$ kHz and $\tau_{1.4} = 0.164^{+0.093}_{-0.044}$ s.

• Detectability:

  • The study estimates the GW energy ($E_{GW}$) required for detection by Advanced LIGO/Virgo and third-generation detectors (Einstein Telescope).
  • While DM-induced softening increases the f-mode frequency (potentially moving it into a more detectable range), the required energy for detection at extragalactic distances (e.g., Virgo cluster) remains far above expected emission levels for current detectors. Detection is deemed plausible only for galactic sources or with significantly improved third-generation sensitivity.

Significance and Claims
The paper claims that modeling DM as a gravitationally captured, non-uniform component within a full general relativistic framework is essential for accurate predictions of NS oscillations. The primary significance lies in demonstrating that:

1. Physical Realism: The non-uniform density profile (governed by $\beta$) significantly alters the structural and oscillation properties compared to constant-density models.
2. Preservation of Universality: Despite the complex microphysics introduced by DM, the fundamental universal relations between oscillation modes and macroscopic observables (compactness, tidal deformability) remain intact. This suggests that URs can still be used as reliable tools to constrain NS properties and potentially infer DM effects from future GW data.
3. Observational Bounds: The study provides specific, calibrated constraints on the f-mode frequency and damping time for canonical NSs in the presence of DM, derived from the GW170817 event, offering a pathway for future multimessenger asteroseismology to probe the nature of dark matter in dense stellar environments.