Limits on dark matter existence in neutron stars from recent astrophysical observations and mass correlation analysis

TL;DR

This study constrains the maximum dark matter content in neutron stars using astrophysical observations, equations of state modeling, and a correlation with maximum neutron star mass, providing insights into dark matter's role in dense astrophysical objects.

Contribution

It establishes a strong linear correlation between dark matter mass fraction and neutron star maximum mass, enabling quantitative constraints on dark matter in neutron stars based on observational data.

Findings

Maximum dark matter mass fraction in neutron stars is approximately 0.15 solar masses.

A strong correlation (r=0.98) exists between dark matter fraction and neutron star maximum mass.

Constraints on dark matter content inform potential observational signatures like gravitational waves.

Abstract

Dark matter admixed neutron stars (DANSs) serve as a specific astrophysical laboratory for probing the features of dark matter (DM) and have emerged as a promising candidate for interpreting recent astrophysical observations (e.g., by NICER and LIGO/Virgo). Accurately constraining the internal DM content of DANSs is therefore of critical importance. In this work, we construct the equations of state (EoS) for DANS matter by employing twelve nuclear matter (NM) models within the covariant density functional (CDF) theory and a self-interacting fermionic model for DM. Using these EoSs as input, we solve the two-fluid Tolman-Oppenheimer-Volkov (TOV) equations to systematically investigate the influence of DM on the global properties of neutron stars (NSs). By incorporating recent observational constraints on NS properties, the maximum DM mass fraction $f_\chi^{\mathrm{max}}$ in DANSs is determined for each NM EoS model. Our analysis reveals a strong linear correlation (Pearson coefficient $r=0.98$) between $f_\chi^{\mathrm{max}}$ and the maximum mass of a pure NS, $M_{\rm{NS}}^{\mathrm{max}}$, described by $f_{\chi}^{\mathrm{max}} = 0.22 M_{\mathrm{NS}}^{\mathrm{max}} - 0.44$. Leveraging this correlation and the observed NS maximum mass distribution, $P(M_{\text{NS}}^{\max} \mid \text{EM})$, we derive the probability distribution function (PDF) for the maximum DM mass, $P(M_{\chi}^{\max} \mid \text{EM})$, in DANSs. We find that at the 68\% confidence level, $M_{\chi}^{\mathrm{max}}=0.150^{+0.070}_{-0.051}\ M_{\odot}$. This quantitative constraint on the DM mass provides a critical prior for interpreting potential observational signatures of DANSs, such as anomalous tidal deformabilities and distinctive gravitational-wave signals.