Towards Uncovering Dark Matter Effects on Neutron Star Properties: A Machine Learning Approach

In recent years, researchers have become increasingly interested in understanding how dark matter affects neutron stars, helping them to better understand complex astrophysical phenomena. In this paper, we delve deeper into this problem by using advanced machine learning techniques to find potential connections between dark matter and various neutron star characteristics. We employ Random Forest classifiers to analyze neutron star (NS) properties and investigate whether these stars exhibit characteristics indicative of dark matter admixture. Our dataset includes 32,000 sequences of simulated NS properties, each described by mass, radius, and tidal deformability, inferred using recent observations and theoretical models. We explore a two-fluid model for the NS, incorporating separate equations of state for nucleonic and dark matter, with the latter considering a fermionic dark matter scenario. Our classifiers are trained and validated in a variety of feature sets, including the tidal deformability for various masses. Based on confusion matrices, these classifiers can identify NS with admixed dark matter with approximately 17% probability of misclassification. In particular, we find that additional tidal deformability data do not significantly improve the precision of our predictions. This article also delves into the potential of specific NS properties as indicators of the presence of dark matter. Radius measurements, especially at extreme mass values, emerge as particularly promising features. The insights gained from our study will guide future observational strategies and enhance dark matter detection capabilities. According to this study, neutron stars at 1.4 and 2.07 solar masses have radii that strongly suggest dark matter in neutron stars more likely than just hadronic composition, based on NICER data from pulsars PSR J0030+0451 and PSR J0740+6620.


Introduction
Neutron stars (NS), because of their dense nature, are of great interest in the universe.Despite their importance in astrophysics, the insides of neutron stars remain a mystery.Neutron stars serve as natural laboratories for studying matter under extreme densities and pressures [1][2][3].Neutron star cores can reach densities approximately 4-5 times greater than normal nuclear saturation density, posing an intriguing astrophysical puzzle.The internal structure of neutron stars, particularly their cores, is one of the most mysterious topics in astrophysics.Despite advances in observational astronomy, current measurements of neutron star properties have not been able to definitively rule out the presence of exotic forms of matter, such as a deconfined quark phase or admixed dark matter.The extreme conditions within these stars make it possible for states of matter to exist that are not yet fully understood or observed under terrestrial conditions [4].Theoretical models suggest a range of scenarios for the core composition, from conventional nucleonic matter to more exotic forms, such as quark matter and strange matter, as discussed by Glendenning [1].
This uncertainty opens the door to considering the potential role of dark matter within the dense core of neutron stars [5][6][7][8].
Dark matter, an invisible yet omnipresent part of the universe, is a major focus in astrophysics.It is estimated to make up 27% of the mass-energy content of the universe, and its gravitational effects are seen in phenomena such as the rotational speeds of galaxies and the dynamics of galactic clusters [9][10][11].However, its mysterious nature and the difficulty in directly detecting it present a major challenge to our understanding of the cosmos.Neutron stars offer an exceptional and promising way to investigate the puzzle of dark matter.The intense gravitational fields of neutron stars make them a potential site for capturing dark matter particles [12][13][14], particularly Weakly Interacting Massive Particles (WIMPs), which are a leading candidate in dark matter models.The concept that dark matter could accumulate in the core of a neutron star and interact with its ordinary matter has generated considerable interest.Goldman and Nussinov (1989) first proposed this idea in their theoretical studies [15], and Kouvaris further explored the implications of dark matter annihilation within neutron stars, suggesting potential observable signs [16,17].The interaction between dark matter and neutron stars is not only about accumulation but also transformation.The dense environment within a neutron star could catalyze processes in which dark matter particles transition into different states [12][13][14].These transitions could affect the neutron star's temperature, it srotation, and even its life cycle.Moreover, the presence of dark matter within neutron stars could create innovative and exotic states of matter, which would challenge our current understanding of compact astrophysical objects.These theoretical possibilities have spurred a growing body of research aimed at using neutron stars as natural detectors for dark matter.Astrophysicists hope to detect signs indicative of dark matter interactions by analyzing the observational data of neutron stars, such as their mass-radius relationships, cooling rates, and magnetic-field structures.This effort not only increases our knowledge of neutron stars, but also potentially helps to solve the dark matter mystery that has long perplexed scientists [13][14][15][16][17][18][19][20][21][22][23][24][25].
The advanced LIGO and VIRGO collaborations have made a groundbreaking discovery of gravitational waves from the GW170817 event, which was a merger of two neutron stars.This, combined with multi-messenger observations of the binary, has opened up new possibilities for investigating the equation of state of dense matter and researching dark matter in compact objects [26,27].In recent years, the scientific community has made considerable progress in understanding the interaction between dark matter (DM) and neutron stars (NS).This research has mainly focused on the effects of DM on the characteristics and detectability of neutron stars.For example, studies have examined the mass-radius relationships and the stability of DM-admixed neutron stars, with a comprehensive analysis of stability aspects [5,28].Additionally, refs.[6,20] have demonstrated that a DM core can significantly reduce the maximum mass of a neutron star.Ref. [29] has extended these findings, discussing the formation of stable DM-admixed NSs and the potential variation in mass-radius relationships and the creation of a DM halo.This study has employed two equations of state for the neutron star's dense baryonic core, made up of piecewise generalized polytropes, and has investigated an asymmetric self-interacting fermionic dark matter component.The focus has been on different scenarios of admixed neutron stars under these specific conditions.Ref. [30] has taken this research further by simulating dark matter admixed neutron star binaries, showing that the presence of DM could shorten the merger remnant's lifespan and affect the brightness of electromagnetic signals.These collective studies have highlighted the potential impacts of DM in neutron stars, such as changes in maximum mass, tidal and surficial Love numbers, and mass-radius relationships.Research has also indicated that the inclusion of DM-bosonic or fermionic-might reduce the higher tidal deformability usually associated with a stiff nuclear equation of state (EOS).This phenomenon has been observed in both single-fluid and two-fluid approaches, as demonstrated in studies such as those by [20,[31][32][33].Furthermore, to better understand the EOS of neutron stars, statistical Bayesian methods have been applied to analyze astrophysical observations, with nuclear matter parameters playing a key role.
Although many studies have hypothesized the presence of DM in neutron stars, few have extensively explored the correlations between observable neutron star properties and dark matter parameters [34].
This research seeks to uncover the secrets of neutron stars and the impact of dark matter on their properties.Despite numerous studies, the universal correlations between neutron star properties in the presence and absence of dark matter remain unclear.To bridge this gap, our study explores these correlations, taking into account the uncertainties in nuclear EOS and dark matter EOS.We use machine learning classification algorithms that are trained on extensive datasets of neutron star properties, encompassing scenarios with and without dark matter, to investigate the feasibility of distinguishing neutron star characteristics based on the presence of dark matter.This machine learning endeavor has the potential to revolutionize our understanding of neutron stars.We also plan to determine the number of observations needed to conclusively classify the existence of dark matter within neutron stars.This could provide a new way to view and understand these enigmatic celestial bodies.Our research combines traditional astrophysical methods with cutting-edge machine learning techniques to venture into uncharted territories of neutron star and dark matter research.
The paper has the following organization.In Section 2, we introduce the basic formalism of the equation of state for nuclear matter and dark matter.In Section 3, we present and discuss the results of the current study.Finally, in Section 4, we provide concluding remarks.

The Equation of State (EOS)
To calculate the neutron star properties with admixed dark matter configurations, we use a two-fluid approach [5,6,32,33].In this approach, separate EOSs are required for two distinct fluids that interact only gravitationally, as described by two-fluid equations.For the nuclear matter component, the EOS is based on only nucleonic degrees of freedom, considering beta-equilibrated cold dense matter.For the dark matter component, we consider a fermionic particle with vector interactions only.Both EOSs are treated within the relativistic mean field (RMF) framework.It should be noted that our limited understanding of dark matter and its minimal constraints allow the validity of various types of modeling, including those involving bosonic particles [32,33,35].Despite this diversity, recent observational constraints do not definitively rule out any particular model.In the following section, we will explain the methodology used to construct both EOS.

Nuclear Matter EOS
We use the nuclear EOS dataset created in ref. [36], using the RMF approach, which only includes nucleonic degrees of freedom.Specifically, "Set 0" in the above-referenced article denotes the full dataset, which covers the entire range of prior.This approach is a mean-field theory that incorporates non-linear meson terms, both self-interactions and mixed terms.Here, a Bayesian inference method is employed to fine-tune the model parameters, guided by constraints such as nuclear saturation properties, the requirement that neutron stars exceed a maximum mass of 2 M ⊙ , and the low-density pure neutron matter EOS as determined by N3LO calculations in chiral effective field theory.This collection of 16,000 equations of state (EOS) is included in the dataset.The Lagrangian density that accounts for the baryonic components is expressed as follows: The Dirac spinor Ψ stands for the nucleon doublet (consisting of neutron and proton) with a bare mass m.The couplings of the nucleons to the meson fields σ, ω, and ϱ are denoted by g σ , g ω , and g ϱ respectively, with their respective masses being m σ , m ω , and m ϱ .The parameters b, c, ξ, and Λ ω , which control the strength of the non-linear terms, are determined in conjunction with the couplings g i (where i = σ, ω, ϱ) by imposing a set of conditions.

Dark Matter EOS
We consider a simple model for the dark matter equation of state that is similar to the Lagrangian of the nuclear model, as employed in ref. [34].This model has a single fermionic component (χ D ) and a dark vector meson V µ D that couples to the conserved DM current through g vd χD γ µ χ D V µ D .The Lagrangian and the corresponding EOS in the mean-field approximation are expressed as follows: (2) The ratio of g vd to m vd is denoted as c ω , and the bare mass of the fermionic dark matter is represented by m χ .These two factors, in addition to the dark matter Fermi momenta, are what determine the dark matter equation of state.The density of the fermionic dark matter is represented by ρ D , which is linked to the mean-field value of the 'dark vector' meson.The Fermi momenta of the dark matter are what determine the amount of dark matter density/mass that is accumulated inside neutron stars.The characteristics of neutron stars that have been mixed with dark matter depend on the dark matter equation of state and the mass fraction of the dark matter.
Our approach to analyzing the structure of NSs with an admixed dark matter configuration with a two-fluid approach involved using four different random nucleonic EOS with varying stiffness to account for uncertainties in the nuclear EOS domain.We also employed a range of dark matter EOS developed using the relativistic mean-field (RMF) approach.For each nuclear matter scenario, we examined a comprehensive set of 4000 unique combinations of dark matter parameters that gives NS maximum mass above 1.9 M ⊙ , which is a value within 3σ of the PSR J0348+0432 mass, 2.01 ± 0.04 M ⊙ .These parameters (c ω , m χ , and F χ ) were selected from uniform distributions within the specified ranges described in Table 1.Notably, the F χ is define a dark matter mass fraction represents the total accumulated dark matter gravitational mass within R χ , where the dark matter pressure reaches zero.In general, we consider 16,000 admixed configurations of the sequence of NS properties, such as mass, radius, and tidal deformability.This sequence represents the entire curve between mass, radius, and tidal deformability.We investigate dark matter with a mass range of 0.5-4.5 GeV [37], a self-interaction measure between 0.1 and 5 fm [29], and a mass fraction of up to 25% [38].For the nucleondark matter interaction within neutron stars, estimates suggest a dark matter capture rate around 10 25 GeV/s for dark matter mass near 1 GeV [39].Over the neutron star's lifetime, this implies a modest accumulation of dark matter, suggesting that alternative mechanisms might be necessary.For instance, dark matter production during supernovae events [40] or the conversion of neutrons into dark matter particles could contribute significantly to dark matter content in neutron stars, as explored in the context of dark matter admixed neutron stars [20].

The Datasets
Our study encompasses a comprehensive collection of 32,000 (32 K) neutron star (NS) properties, including mass-radius curves and mass-tidal deformability curves.Half of these (16,000, or 16 K) are derived from models with dark matter, and the remaining 16 K are from models without dark matter.To facilitate the machine learning classification task, we structured our datasets as follows.
We constructed a feature vector X comprising five randomly chosen masses, uniformly distributed between 1 solar mass and the maximum mass specified by each curve.Similarly to these five masses, we included five radii and five tidal deformabilities in the feature vector.Thus, each vector X consisted of 15 elements in total: 5 masses, 5 radii, and 5 tidal deformabilities.
For the target vector Y, we assigned a value of zero to the curves that represented nuclear-matter configurations only and a value of one to those that included dark matter configurations.Consequently, our target vector Y was a binary indicator, with zero signifying the absence of dark matter and one indicating its presence in the model.This binary classification system effectively distinguished between neutron star models with and without dark matter influence.

Results
In this section, we aim to examine the observational properties of neutron stars.Our analysis focuses on results obtained using a varied collection of nuclear matter EOS without dark matter, as well as a range of dark matter EOS.In the two fluid approach, four different nuclear EOS are considered, each with a different stiffness.
In Figure 1, we show the equation of state for both nuclear matter (highlighted in orange) and dark matter (highlighted in blue).The nuclear matter equation of state was obtained with minimal constraints on nuclear saturation properties, neutron star masses above two solar masses, and low density constraints from chiEFT theory [36].The dark matter equation of state was derived by varying the model parameters within the range given in Table 1.This range is wider than that for nuclear matter, as there are fewer constraints on dark matter equations of state.Additionally, the graph displays the domain of equations of state derived from GW170817, represented by the cross-hatched overlay.It should be noted that this domain is only applicable to nuclear matter equations of state, as it was derived using a single-fluid approach.Figure 2 presents the relationship between mass, radius, and tidal deformability for neutron stars, based on a set of 16,000 equations of state that describe nuclear matter without the influence of dark matter.In the mass-radius plot on the left, the gray areas represent the confidence intervals (CI) for the binary components observed in the GW170817 event, with solid lines indicating a 90% CI and dashed lines indicating a 50% CI.Additionally, we have incorporated data from NICER's X-ray observations of the pulsars PSR J0030+0451 and PSR J0740+6620.For PSR J0030+0451, the cyan and yellow shaded regions correspond to the 1σ (68%) confidence zone from its 2-D posterior distribution in the mass-radius space.The violet shaded region pertains to PSR J0740+6620.The error bars displayed-horizontal for radius and vertical for mass-represent the credible intervals 1σ from the marginalized posterior distributions of 1D derived from the same NICER data.In the right panel, which plots the mass against the tidal deformability, the blue bars mark the specific measurements for PSR J0740+6620: its radius at a mass of 2.08 M ⊙ and its tidal deformability at 1.36 M ⊙ .From the figure, it is evident that the nuclear matter EOS set we have used encompasses a wide array of neutron star properties.The maximum mass of these neutron stars ranges between 2 and 2.75 solar masses (M ⊙ ), indicating a wide spread in both the radius and tidal deformability values across the set, demonstrating the versatility of the EOS employed.Crucially, our results show good alignment with the various astrophysical constraints derived from recent observations.The figure referenced as Figure 3 illustrates the mass-radius (left panel) and mass-tidal deformability (right panel) relationships for neutron stars, based on 16,000 equations of state that include dark matter.These properties have been computed using a two-fluid approach, as previously described.The nuclear matter component in the two-fluid models is represented by four randomly selected EOS from the aforementioned set of 16,000 nuclear matter EOS.These four were specifically chosen for their varying stiffness, effectively spanning the current range of uncertainty within the nuclear EOS.It is important to note that the left panel of the figure presents the variation of the neutron star mass in relation to the total radius (R t ).This total radius (R t ) is defined as the point at which the pressure of both fluids, nuclear matter and dark matter, disappears.In this context, F χ represents the fraction of the total neutron star mass that is constituted by dark matter.An increase in F χ results in a corresponding decrease in both the mass and the radius of the neutron star, as shown in the left panel.Similarly, the right panel shows that dimensionless tidal deformability also decreases with an increase in F χ .The study cited as [34] concluded that F χ is strongly correlated with neutron star properties at any given mass, provided the nuclear matter EOS remains constant.However, when variations in nuclear matter EOS are introduced, effectively incorporating the uncertainty in the nuclear EOS, the correlations with F χ dissipate.This finding underscores the need to disentangle the effects of uncertainties stemming from both nuclear EOS and dark matter on the properties of neutron stars.Figure 4 examines the robustness of certain universal relations in the context of neutron star properties, particularly when dark matter is present.Previous research has suggested that some of these relations may not be valid when quark matter is taken into account.Therefore, our investigation focuses on the C − Λ relation, which has been found to be unaffected by the equation of state when only nucleonic degrees of freedom are considered, to determine the impact of admixed dark matter.The C-Love relationship is depicted in Figure 4 with light pink scatter plots for nuclear EOS sets without dark matter and light blue scatter plots for those with admixed fermionic dark matter.We used a simple polynomial curve, as described by the equation below, to fit the EOS data: The fitting process for the nuclear EOS data yielded coefficients a 0 = 0.36935858 ± 0.00001225, a 1 = −0.04011723± 0.00000397, and a 2 = 0.0010823 ± 0.0000003.In comparison, the admixed dark matter EOS set resulted in a 0 = 0.36054566 ± 0.00002011, a 1 = −0.0375908± 0.00000931, and a 2 = 0.00086283 ± 0.00000093.The lower panel in Figure 4 presents the absolute differences from these fits, revealing a deviation of approximately 1% that diminishes to 0.5% for higher values of Λ, with the nuclear matter EOS case consistently close to 0.5%.Our findings were juxtaposed with the results from Figure 4 in ref. [47], which reported on a constrained EOS with coefficients a 0 = 0.3617, a 1 = −0.03548,and a 2 = 0.0006194.The constrained EOS exhibited an absolute difference of less than 1%, whereas the unconstrained EOS showed around a 1% deviation.From this comparative analysis, we infer that dark matter does not significantly disrupt the universal C − Λ relation, suggesting that the relationship between compactness C and tidal deformability Λ remains stable even when dark matter is considered within the NS structure.The red line indicates the fit for the dark matter EOS, while the gray line shows the fit for the nuclear matter EOS, both fitted using Equation ( 5).The lower panel of the figure displays the residuals from these fittings.Additionally, the figure includes a comparison with the results from ref. [47], specifically focusing on the single-fluid Tolman-Oppenheimer-Volkoff (TOV) scenario without dark matter, which is highlighted in a pink-red color.
In our analysis, depicted in Figure 5, we explore the correlation between the radii of neutron stars with masses of 1.4 and 2.07 solar masses (M ⊙ ), comparing scenarios with and without dark matter.Inspired by a recent study [48], which presents a method to detect phase transitions in neutron stars by examining the radius correlations, including NICER's observations of PSR J0740+6620 and PSR 0030+0451, we extend this approach to investigate the presence of dark matter.We plotted the relationship between the radii of neutron stars with masses of 1.4 and 2.07 solar masses (M ⊙ ), comparing scenarios with and without dark matter.In the dark matter scenario, the total radius of the neutron star includes the dark matter component.We found that the Pearson correlation coefficient is 0.77 for the dark matter set and 0.85 for the nuclear matter set, indicating a slight weakening of the correlation with dark matter.Furthermore, we compared the constraints derived from NICER in blue, which were obtained by marginalizing the 1σ posteriors of the NICER measurements of these two pulsars over the NS mass.Then, the radius data of these two observations yielded this blue domain.This analysis considered only the data analyzed by the Riley et al. group.Significantly, we observed that the slope of the relationship between R 2.07 and R 1.4 undergoes a considerable change in the presence of dark matter.The slope is 1.85 in the absence of dark matter and 1.16 with admixed dark matter.The line fitted to each scenario is plotted along with the uncertainty band of 1σ.The presence of dark matter appears to reduce this slope.In particular, the slope for the NICER constraints is 0.85 ± 0.07, showing substantial overlap with the dark matter admixed set.This overlap suggests that the NICER measurements of the radii of PSR J0740+6620 and PSR 0030+0451 might indicate a higher likelihood of the presence of dark matter compared to solely hadronic matter.
The research conducted thus far has strongly encouraged us to explore the potential for dark matter to be revealed in the observational properties of neutron stars, such as mass, radius, and dimensionless tidal deformability, through the use of machine learning.The machine learning tools are becoming increasingly important in neutron star physics [49][50][51][52][53][54].We use the robust dataset of neutron star properties of both nuclear matter and admixed dark matter configurations that were previously discussed.Our goal is to create a classification system that can distinguish between the presence and absence of dark matter in neutron stars.We create two features, denoted as X1 and X2.Feature X1 consists of mass and radius values, while X2 includes an additional element-tidal deformability.These values are generated through random uniform sampling, with the range for mass and radius extending from 1 to their respective maximum values.The construction of X1 involves a total of 10 elements, each comprising 5 values of mass and radius.On the other hand, X2 has an added dimension of 5 tidal deformability measurements, bringing its total length to 15 elements.Section 2.3 discusses structuring the dataset for machine learning training.It should be noted that the target vector Y has only one element, which is set to zero for X with only nuclear matter and 1 for admixed dark matter.
In our study, we employ the Random Forest classifier [55], recognized as one of the most effective classification algorithms, to analyze our data.To begin, we combine two datasets, each containing 16,000 entries, ensuring a thorough shuffle to guarantee uniformity in the data distribution.Then, this merged dataset is divided into three distinct parts: 60% for training, 20% for validation, and the remaining 20% reserved for testing.
The training set, which constitutes 60% of the data, is used to teach the Random Forest model the patterns and relationships that are inherent to the data.The validation set, comprising 20% of the data, plays a critical role in fine-tuning the model.During this phase, we adjust the hyperparameters of the Random Forest classifier.Hyperparameters are the configuration settings used to structure the learning process and can significantly impact the performance of the model.By tweaking these parameters while observing the model's performance on the validation set, we aim to find the optimal configuration that yields the best results.This process of hyperparameter tuning is essential to prevent issues like overfitting, where the model performs well on the training data but poorly on unseen data.The final 20% of the data, which are discarded as a test dataset, are used to evaluate the performance of the fully trained model.This dataset is crucial, as it represents new unseen data for the model, providing a realistic assessment of how the model will perform in real-world scenarios.It is important to note that these test data are not involved in the training or validation process.By keeping these data separate, we ensure that our evaluation of the model's effectiveness is unbiased and indicative of its true predictive capabilities.
In Figure 6, two panels represent confusion matrices for a Random Forest classifier applied to train model with only the X1 feature.The left panel shows the prediction of the trained model on the test set and the right panel shows the entire dataset.A confusion matrix is a table that is often used to describe the performance of a classification model on a set of data for which the true values are known.It cross-tabulates the actual class labels with the predicted class labels, providing insight into the accuracy and types of errors made by the classifier.In the left panel, the confusion matrix for the test set reveals that the classifier accurately predicted NS with only nuclear matter (NM) 1872 times and with admixed dark matter (DM) 1644 times.However, there were instances of misclassification, indicated by off-diagonal numbers: 128 instances of NM were incorrectly classified as DM and 356 instances of DM were incorrectly classified as NM.These errors highlight the instances where the classifier was challenged to distinguish between the two classes.Moving to the right panel, the confusion matrix for the entire dataset shows that the classifier correctly identified NM 15,584 times and DM 15,458 times.Misclassifications are also present here, with NM being mistaken for DM 226 times and vice versa 542 times.The percentage of false positive for NS with only nuclear matter (NM) is 6.4% (1.98%) and NS with admixed dark matter (DM) is 17.8% (4.98%) for the set of tests (entire).It should be noted that the test set contains data separated from those not involved in training, but the entire set contains data involved in training.Figure 7 shows a confusion matrix for the model trained in the X2 feature set, which includes five additional features: tidal deformability for five different mass categories.The purpose of introducing these additional features was to determine whether the inclusion of tidal deformability constraints could improve the accuracy of the predictions.The results show that the rate of false positives, incorrectly predicting NM as dark matter DM, is 5.9% for the test set and 1.7% for the entire dataset.On the other hand, the rate of false positives for neutron stars with admixed dark matter is 17.1% for the test set and 4.69% for the entire set.These figures suggest that the inclusion of tidal deformability data does not significantly enhance the predictive power of the model.As Random Forest is a tree-based model, one can also extract the importance of different features (or the mean decrease in impurity) within a trained model in predicting the target variable.This is accomplished by monitoring the reduction in model precision when each feature is taken out of the model.To put it more simply, it shows the impact of each feature on the model's decision-making process, with higher numbers indicating a greater effect on the classification result.In Figure 8, we plot the importance of the feature that delineates the varying significance of the features employed by the Random Forest model in the context of classifying the existence of dark matter within neutron stars (NS).Radius measurements at various mass points ('R(M1)' through 'R(M5)') dominate the feature importance, with the plot revealing that 'R(M1)', the radius at a lower mass point, holds the greatest predictive power, followed by 'R(M4)', a higher mass point.This pattern suggests that radius measurements, especially at the extremes of the mass range considered, are critical for the model to discern the presence of dark matter in the NS.Tidal deformability features ('Λ(M1)' through 'Λ(M5)') are also included, although they exhibit less influence on the model's decisions.This insight underscores the necessity of precise radius measurements over a range of masses as a more determinant factor in revealing or classifying dark matter within NS.The comparative lower importance of the tidal deformability suggests that, while useful, they do not contribute as significantly to the model's classification ability as the radius measurements.

Conclusions
In this article, we explored the combination of dark matter and neutron stars through a two-fluid model, concentrating solely on a fermionic dark matter equation of state.To replicate this situation, we carefully employed EOS for nucleonic matter and dark matter within a relativistic mean-field (RMF) approach, which included a comprehensive collection of 16,000 different EOS.The nuclear EOS was the result of Bayesian inference with minimal constraints on nuclear saturation properties, an NS maximum mass of more than 2 solar masses, and a low-density pure neutron matter constraint.For the dark matter component, a straightforward model was adopted, featuring a singular fermionic entity interacting with a dark vector meson.The resulting dark matter EOS was contingent on the interaction strength with the dark vector meson and the intrinsic mass of the fermionic dark matter.The interplay between the dark matter EOS and the proportion of dark matter present critically determined the resultant properties of the neutron stars.
A comprehensive study was carried out on the combination of neutron stars and dark matter, using four different nucleonic equations of state and a range of dark matter equations of state.This research looked at 32,000 neutron stars properties curve, examining their mass, radius, and tidal deformability.The study also delved into the integrity of universal relations in the presence of dark matter, affirming the stability of the compactnesstidal deformability (C-Λ) relation.Furthermore, we investigate the correlation between the radii R 1.4 and R 2.07 for neutron stars with masses 1.4 and 2.07 M ⊙ , respectively.The two radii are closely related, but the inclusion of dark matter in neutron stars can alter the gradient.For the first time, our analysis has taken advantage of the most recent observational data from NICER, which strongly suggests that dark matter must be included in neutron star models.As our knowledge advances and our observational powers increase, the intricate part that dark matter plays in the universe will become more distinct, giving us new perspectives on the essential character of our cosmos.
We employed machine learning techniques, particularly Random Forest classifiers, to classify neutron stars (NS) in the presence of dark matter, based on their inherent properties.We evaluated the predictive accuracy of these classifiers, trained on various feature sets, using confusion matrices.Our findings suggest that the properties of NS with a mixed dark matter configuration can be identified with approximately 17% probability of misclassification as nuclear matter NS.Adding constraints related to tidal deformability does not significantly reduce this uncertainty.This figure reflects both the inherent uncertainty of the classifier model and the variability in NS properties.Future research should include a more focused analysis to distinguish between these sources of uncertainty, as well as a broader examination using a variety of machine learning algorithms.Additionally, our analysis sought to identify critical features that would be useful for future observations to detect dark matter within NS.The results suggest that measurements of radius at both lower and higher masses appear to be very promising indicators for the presence of dark matter.

Figure 1 .
Figure1.The 90% confidence interval of the dark matter (nuclear matter) equation of state is represented by the shaded blue (light orange) area.This is the pressure P χ (P B ) as a function of density ρ χ (ρ B ).The hatched grey band is the prediction from the GW170817 event.

Figure 2 .
Figure2.The mass-radius and mass-tidal deformability curves for the 16,000 nuclear matter equation of state, i.e., the properties of neutron stars without dark matter, are shown, respectively.The gray zones in the left panel indicate the 90% (solid) and 50% (dashed) confidence intervals for the binary components of the GW170817 event[41].The cyan and yellow zones in the same panel represent the 1σ (68%) credible area of the 2-D posterior distribution in the mass-radius domain from millisecond pulsar PSR J0030+0451[42,43], while the violet zone is from PSR J0740+6620[44,45], both derived from NICER X-ray data.The horizontal (radius) and vertical (mass) error bars reflect the 1σ credible interval derived from the same NICER data's 1-D marginalized posterior distribution.The blue bars depict the radius of PSR J0740+6620 at 2.08 M ⊙ (left panel) and its tidal deformability at 1.36 M ⊙ (right panel)[46].

Figure 3 .
Figure 3. Same as Figure 2, but the properties of neutron stars with a dark matter configuration of 16,000 EOS are depicted.The fraction of admixed dark matter over the total mass of the neutron star, F χ , is also shown.see text for details.The range of colors, starting from yellow and progressing to dark red, indicates the varying mass fraction of dark matter, ranging from zero to the maximum percentage shown.

Figure 4 .
Figure 4.The figure illustrates the C-Love universal relation for EOS involving only nuclear matter (depicted in light blue) and those that include dark matter (shown in light pink).The red line indicates the fit for the dark matter EOS, while the gray line shows the fit for the nuclear matter EOS, both fitted using Equation (5).The lower panel of the figure displays the residuals from these fittings.Additionally, the figure includes a comparison with the results from ref.[47], specifically focusing on the single-fluid Tolman-Oppenheimer-Volkoff (TOV) scenario without dark matter, which is highlighted in a pink-red color.

Figure 5 .
Figure 5.The figure displays the relationship between the radii of neutron stars with masses of 1.4 and 2.07 solar masses (M ⊙ ).The comparison is made between two different scenarios: one with dark matter and one without.The radius in the scenario that includes dark matter is the total radius of the neutron star, including the dark matter component.

Figure 6 .
Figure 6.The confusion matrices displayed are for the test set and the entire set, respectively (model trained with X1 feature, i.e., NS mass and radius).

Figure 7 .
Figure 7. Same as shown in Figure 6, but with the trained model on the X2 feature; it has five additional features compared to the previous, which are the tidal deformability for five different masses.

Figure 8 .
Figure 8. Feature importance for Random Forest classification from X2.

Table 1 .
The varied range of dark matter model parameters.