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

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.