General relativistic study of $f$-mode oscillations in neutron stars with gravitationally bound dark matter
Pinku Routaray
TL;DR
This paper investigates nonradial $f$-mode oscillations in neutron stars admixed with gravitationally bound dark matter using full general relativity. The authors adopt a single-fluid, Higgs-portal DM model with a physically motivated non-uniform DM density described by two parameters, $\alpha M_χ$ and $β$, embedded into a relativistic mean-field hadronic EOS, and solve the TOV equations together with GR perturbations to obtain complex QNM frequencies via the Zerilli formalism. They derive analytic fits for the $f$-mode frequency versus compactness and tidal deformability ($f$-$C$-$τ$ and $f$-$Λ$-$τ$ relations), and examine correlations and universal relations (URs) in the presence of DM, including a GW170817-driven constraint on the canonical NS oscillation properties. The results show that DM softens the EOS, lowering the maximum mass and radius while increasing the $f$-mode frequency and decreasing the damping time; URs largely survive DM admixture, enabling DM–NS properties to be constrained by future multimessenger observations. Overall, the study provides a framework to connect microscopic DM parameters to macroscopic NS observables and GW signals, highlighting how asteroseismology can help bound Higgs-portal DM scenarios in light of GW data.
Abstract
A comprehensive investigation of nonradial oscillations in neutron star (NS) admixed with gravitationally bounded dark matter (DM) is carried out within the framework of full general relativity. The relativistic mean field (RMF) formalism is employed to illustrate the hadronic equation of state (EOS), while a physically motivated, gravitationally captured, non-uniform fermionic Higgs-portal DM component is incorporated to model DM-admixed NS. The DM distribution is characterized by two free parameters: $αM_χ$, an effective scaling factor that combines the DM concentration and the DM candidate mass, and $β$, a steepness index controlling the DM density distribution. The quasi normal mode (QNM) characteristics such as fundamental ($f$) mode frequency and its corresponding gravitational-wave (GW) damping time ($τ$) is calculated for DM-admixed NS by solving the general relativistic perturbed equations involving axial as well as polar modes. The study demonstrates how the inclusion of DM distribution modifies the $f$-mode frequency and enhances the damping rate, reflecting a stronger coupling between matter and spacetime perturbations. Considering DM effects, the correlation analysis among DM model parameters, NS observables and QNM characteristics also carried out. Analytic fits for the $f-C-τ$ and $f-Λ-τ$ relations are constructed and calibrated for DM-admixed NS models. Building upon asteroseismic universal relations (URs), multimessenger constraint from the GW170817 event is employed by mapping the tidal deformability $Λ_{1.4}$ into the $(f_{1.4},τ_{1.4})$ space, thereby providing observational bounds on the oscillation properties of canonical DM-admixed NS model.
