Exploring the distribution and impact of bosonic dark matter in neutron stars

The presence of dark matter (DM) within neutron stars (NSs) can be introduced by different accumulation scenarios in which DM and baryonic matter (BM) may interact only through the gravitational force. In this work, we consider asymmetric self-interacting bosonic DM which can reside as a dense core inside the NS or form an extended halo around it. It is seen that depending on the boson mass ($m_{\chi}$), self-coupling constant ($\lambda$) and DM fraction ($F_{\chi}$), the maximum mass, radius and tidal deformability of NSs with DM admixture will be altered significantly. The impact of DM causes some modifications in the observable features induced solely by the BM component. Here, we focus on the widely used nuclear matter equation of state (EoS) called DD2 for describing NS matter. We show that by involving DM in NSs, the corresponding observational parameters will be changed to be consistent with the latest multi-messenger observations of NSs. It is seen that for $m_{\chi}\gtrsim200$ MeV and $\lambda\lesssim2\pi$, DM admixed NSs with $4\%\lesssim F_{\chi}\lesssim20\%$ are consistent with the maximum mass and tidal deformability constraints.


Introduction
Owing to the fact that DM constitutes the majority of matter in galaxies, several noteworthy ideas have been proposed about the presence of DM inside compact astrophysical objects [1][2][3][4][5][6][7][8][9].Among them, NSs, due to their content of high-density matter and extreme gravitational potential, provide an interesting astronomical environment where sizeable amounts of DM may be accumulated based on various scenarios [10][11][12][13][14].In this regard, the accretion of DM might occur during star evolution or over the lifetime of a NS.It is argued that high capture rates are more likely to happen towards the center of galaxies, where the density of DM is increasing [15][16][17].Moreover, DM production inside NSs and supernova explosions can lead to large values of the DM fraction [5,6,18,19].Among various models, a neutron decay anomaly by assuming a dark sector attracts considerable attention [20][21][22].It is notable to say that gravitationally stable objects so-called dark stars, may be formed entirely from DM [23][24][25][26], which can be considered as a capturing center for BM or even merge with NSs resulting in the formation of a mixed compact object [5,27,28].
In recent years, multi-messenger observations of NSs, have provided a unique opportunity to probe their internal structure and the possible existence of exotic configurations including DM [29,30].The impact of DM on NS properties has been investigated in various studies proposing a wide range of smoking guns [31][32][33][34][35][36][37].Generally, the presence of DM can be considered through two scenarios, i) a single-fluid object for which DM and BM interact also non-gravitationally via an EoS [18,19,[38][39][40][41][42].ii) Two-fluid DM admixed NSs where DM and BM are described by two EoSs and interact only by the gravitational force [10,[43][44][45][46][47][48][49].For the former approach, DM is distributed throughout the whole NSs while for the latter, DM can form a dense core inside NSs or an extended halo around it.Furthermore, in the single-fluid mixed object, DM usually softens the EoS, however, in the case of twofluid models, DM indicates both softening and stiffening behavior due to the formation of DM core or halo, respectively.Therefore, considering both approaches, it is seen that asymmetric DM could affect the maximum mass, radius and tidal deformability of NSs significantly [5,6,10,11,19].Moreover, for asymmetric DM, in which a particle-antiparticle asymmetry exists in the dark sector [50], and also self-annihilating (symmetric) DM, the cooling and heating of NSs could be altered which may help to interpret the thermal evolution of NSs [51][52][53][54][55].In addition, DM core or halo formations within NSs can influence the inspiral phase of their binary systems and also the corresponding gravitational wave (GW) signals especially in the post-merger stage [56][57][58][59][60][61].Thanks to the precise mass-radius measurements based on pulse-profile modeling by the NICER telescope [62][63][64], a novel approach has been introduced to probe DM halo around NSs [11,65], where the variation of the minima in pulse profiles can be severed as a notable signature of DM halo.
Among various DM candidates, bosonic DM models such as scalar fields, axions and sexaquarks are of great interest in various aspects of astrophysics and have been examined via their impacts on compact objects [19,66,67].In this research in order to model DM, we consider a self-repulsive complex scalar field with sub-GeV bosonic particles and coupling constant in the order of unity.Our DM EoS was first used in the pioneering work by Colpi et al. [68] to describe boson stars [69,70], and more recently revisited to model the two-fluid DM admixed NSs [10,11,33].For the BM component, we employ a widely used and wellknown nuclear matter EoS called DD2.Our utilized DD2 EoS is a specific parameterization for the density-dependent meson couplings in a generalized relativistic density functional (GRDF) for nuclear matter which consider σ, ω, and ρ mesons as exchange particles for describing the effective in-medium interaction of nucleons [71].Due to the unknown structure of DM admixed NSs, various nuclear matter EoSs including the ones with phase transition, strange quark matter, and hyperonic degrees of freedom have been proposed to be considered in the mixed compact objects [72][73][74][75][76][77][78][79].
Mass and radius measurements of NSs through observations of radio, optical and X-ray emissions provide reliable bounds for the maximum mass, M max ≳ 2M ⊙ [80][81][82], and the radius of 1.4M ⊙ NS, R 1.4 ≳ 11 km [29,83,84].Furthermore, in addition to constraints for mass and radius, the GW detection of binary NS mergers by the LIGO/Virgo collaboration [85], has introduced the tidal deformability (Λ) as a novel parameter to measure the deformation response of compact objects in a binary system [86].The tidal deformability, which is highly sensitive to the EoS of compact objects and their compactness, is inferred from the detected GW signal in GW170817 event and for a 1.4M ⊙ NS has an upper limit, Λ 1.4 ≲ 580 [87].These multi-messenger observations have been used extensively to constrain dense nuclear matter properties and to explore the possible existence of DM inside NSs by constraining the DM parameter space such as particle mass, coupling constant and its fraction inside the mixed object [10,11,13,14,46,88].Recently, by investigating the consistency of the DM admixed NS observable features with the above-mentioned astrophysical limits, an excluded region for DM parameters and a maximum possible fraction of DM have been obtained in [10,11].
DD2 EoS is considered one of the preeminent nuclear EoSs for which the density dependence of the couplings is adjusted to describe the properties of atomic nuclei around saturation density.However, since DD2 is a stiff EoS, it does not yield a value for Λ 1.4 that is consistent with astrophysical constraints.Considering the aim of the paper, we investigate the impact of bosonic DM on the maximum mass, radius, and particularly the tidal deformability derived from the DD2 EoS.Indeed, we are going to show how the presence of DM may modify the observable characteristics of NSs derived from nuclear EoSs.The allowed bosonic DM parameter space and fraction will be investigated for which all of the aforementioned parameters are in agreement with observational bounds.Our result indicates a modification of NS features obtained from the considered nuclear matter model thanks to the presence of DM.In light of current and upcoming promising astrophysical missions, we provide strong hints about the possibility of the existence of DM in NSs and its influence on various measurable parameters.
The paper is organized as follows.In Sec. 2 we present a two-fluid formalism for DM admixed NSs in addition to introducing the BM and DM EoSs.Various distributions of bosonic DM in NSs are considered in Sec. 3. The impact of DM on the tidal deformability of the mixed objects is examined in Sec. 4, in Sec. 5 the parameter space of DM model is considered in light of the latest multi-messenger constraints and finally the conclusion is given in 6.In this paper, we use units in which h = c = G = 1.

Two-fluid DM admixed NSs
Owing to the fact that in two-fluid DM admixed NSs, BM and DM interact only via the gravitational force, the energy-momentum tensors are conserved separately T µν = T µν DM + T µν BM .Thus, solving the Einstein equation leads to the two-fluid Tolman-Oppenheimer-Volkof (TOV) relations as follows [43,44,89] where m B (R B ) and m D (R D ) are the enclosed masses of BM and DM, respectively.If DM resides as core inside NS, R B > R D , and the outermost radius of the object is the BM radius.However, when DM forms a halo around a NS, R D > R B , and the radius of the object is defined by the DM fluid.It is worth mentioning that for both possible DM configurations, the visible radius of the star is R B .Moreover, the amount of DM inside the mixed object which is considered a model parameter in our analysis is defined as M T .In this research, the DM fluid is described with the following Lagrangian where a self-interacting complex scalar field, ϕ, has been considered as the bosonic DM with a repulsive interaction defined by V(ϕ) = λ 4 |ϕ| 4 potential for which λ is the dimensionless coupling constant and m χ is the particle mass [68].In the strong coupling regime, the system can be approximated as a perfect fluid and the anisotropy of pressure will be ignored, so one can reach the EoS of the self-repulsive bosonic DM as follows Where P and ρ are the pressure and density of the self-interacting bosonic DM (see the appendix of [10] for more details about the bosonic DM EoS in which a derivation has been proposed based on the mean-field approximation in a flat space-time).Here, the well-known DD2 model is utilized to describe the EoS of BM.The DD2 parameterization was obtained by fitting the properties of finite nuclei such as binding energy, charge and diffraction radii, surface thickness, and spin-orbit splitting [71].The predicted nuclear matter parameters at saturation point are consistent with recent constraints.In particular, a saturation density of 0.149065 fm −3 , a binding energy per nucleon of 16.02 MeV, and an incompressibility of 242.7 MeV are obtained for the symmetric nuclear matter within DD2 model [90].The symmetry energy at saturation and its slope parameter are given by J = 31.

DM distributions in NSs
In this section, we are going to investigate the distribution of self-repulsive bosonic DM inside NSs by considering the DD2 EoS as the BM component.Therefore, in Fig. (1) the variation of R B (solid line) and R D (dashed line) is depicted as a function of the DM fraction (F χ ) for various boson masses (m χ ) and coupling constants (λ).It should be mentioned that each colored line, both solid and dashed, is associated with a specific value of m χ (left panel) or λ (right panel), as labeled.Note that the radii for each given DM model parameter and fraction correspond to the maximum mass of the DM admixed NSs.For both plots, it is seen that in low F χ the DM radius (dashed lines) is smaller than the BM radius (solid lines) and a DM core is formed inside the NS.However, toward higher DM fractions, R D will increase and surpass the corresponding R B leading to a DM halo formation around the BM shell.Therefore, a transition from DM core to DM halo formation can be seen by enhancing F χ for each considered case.In the left panel, it is seen that heavier bosons (m χ ≳ 200 MeV), in comparison to lighter ones, reside mainly in the core inside NSs while for high DM fractions, they produce a small halo around the NS.Moreover, in the right panel where the impact of the coupling constant is considered for a fixed boson mass (m χ = 250 MeV), higher couplings result in larger DM halos for a given F χ .Obviously, the threshold DM fraction for which a DM core to DM halo transition occurs will increase by decreasing the coupling constant.Regarding the visible radius of the object (R B ), both of the plots show a decrease along higher DM fractions which for heavy bosons and lower couplings the rate of reduction is larger.For the sake of completeness, we performed a scan over the F χ − m χ parameter space for DM admixed NSs with M T = 1.4M ⊙ and fixed λ = 0.5π.The results are presented in Fig. (2), which shows the contours representing the values of R B /R D .It is seen that light bosons form an extended halo around NSs whilst heavier bosons make a dense DM core inside the star.The yellow line R B = R D = 1 depicts the region where DM is distributed inside the whole NS and a transition occurs from DM core to DM halo for any given parameter.Furthermore, for the regions where a halo is formed, higher values of F χ lead to a very large halo (R B /R D ≲ 0.1), however, for the part where a core is formed, lower values of DM fractions result in a very small DM core (R B /R D ≳ 3) in the BM shell.
Regarding the aforementioned results about DM core and DM halo formations, in the following sections, we are going to investigate the impacts of self-interacting bosonic DM on NSs observable features.

Tidal deformability of DM admixed NSs
The GW detection of a NS binary merger in the GW170817 event and measuring its tidal deformability provide a novel approach to probe high-density matter and the possibility to explore the existence of more exotic compact objects including DM.In this section, we investigate the tidal deformability (Λ) of DM admixed NSs for various bosonic DM model parameters and compare it with the pure NS which is modeled by the DD2 EoS.The dimensionless tidal deformability is defined as in which R and M are the outermost radius and total mass of the object, and k 2 is the tidal love number [91] calculated from the two-fluid TOV formalism [10].Notice that Λ for DM admixed NSs with 1.4M ⊙ is highly dependent on the outermost radius of the object which will be varied when DM core or DM halo is formed inside or around NSs.
In the following, Figs. 3, 4 and 5 indicate the tidal deformability of DM admixed NSs with respect to the mass of the object for various bosonic DM model parameters and fractions.In all of the plots, the curve related to pure BM NS (without DM) is shown by a black solid line and the magenta vertical line depicts the allowed observational range for the tidal deformability of a 1.4M ⊙ NS, Λ 1.4 = 190 +390 −120 reported by [87].Obviously, it is seen that the corresponding Λ 1.4 of DD2 EoS is higher than the upper allowed limit for tidal deformability, thus we will examine how bosonic DM will alter this parameter which may become compatible with the GW170817 constraint.
In Fig. (3), different m χ are considered for fixed F χ = 10% and λ = π.We see that light bosons increase Λ because of the rise in the outermost radius of the object and the formation of a large DM halo.However, for massive bosons (m χ ≥ 300 MeV) which make the star more compact, the tidal deformability will be reduced and is consistent with the allowed astrophysical limit.4) illustrates that for F χ ≲ 20%, the tidal deformability is improved, however, for higher fractions both maximum mass and Λ are not consistent with the observational constraints.Indeed, it is displayed that by raising F χ , the tidal deformability is enhanced which is due to a transition from a DM core to a DM halo formation by varying the amount of DM.Moreover in Fig. (5) we probe the impact of a changing coupling constant on the tidal deformability of DM admixed NSs for F χ = 10% and m χ = 200 MeV.It is seen that by decreasing the coupling constant, the tidal deformability will be reduced and becomes consistent with the observational range while the maximum mass is also in agreement with the 2M ⊙ limit.However, for λ ≳ 2π, we see an enhancement in the tidal deformability which is higher than the astrophysical bound because of the DM halo surrounding the NSs.
Therefore, concerning the aforementioned results, in the next section, we will probe the bosonic DM parameter space for which both M T max and Λ 1.4 are in agreement with the astrophysical limits and indicate the capability of bosonic DM in modifying the observable features related to the BM EoS.

Probing bosonic DM parameter space
In this section, we will examine the bosonic DM model parameters i.e. boson mass and coupling constant and its fraction, by considering observable quantities of the mixed compact object.The total maximum mass and tidal deformability of a 1.4M ⊙ DM admixed NS are taken into account to probe the DM parameter space which is consistent with both the M max ≥ 2M ⊙ and Λ 1.4 ≲ 580 constraints.It should be mentioned that DD2 is a stiff EoS which results in a large radius around 1.4M ⊙ , R 1.4 = 13.15 km, thus the R B 1.4 ≥ 11km constraint will be generally satisfied for this investigation, however, for more details about the variation of dark and visible radius of the star see [11].
Concerning DM core and DM halo formation, Fig. ( 6) illustrates the variation of M T max with respect to DM fraction.In both plots, for all of the considered DM parameters, the maximum mass of the object will decrease in low fractions which correspond to the DM core formation.However, by increasing F χ , the maximum mass of the star gradually rises due to the DM halo configuration.In the left panel, it is seen that for light bosons M T max increases sharply and is higher than the 2M ⊙ limit (shown by the gray dashed line) for the whole range of DM fractions.For more massive bosonic particles (m χ ≳ 200 MeV) there is a DM fraction limit where the maximum mass does not agree with the astrophysical bound and falls below 2M ⊙ .Regarding the right panel, it is indicated that for λ = 4π, M T max is always consistent with the observational constraint, however, for lower coupling constants, one can find a DM fraction beyond which M T max < 2M ⊙ .Furthermore, both plots demonstrate that for small amounts of DM fractions, the total maximum mass of DM admixed NSs is in agreement with M max ≥ 2M ⊙ constraint for all of the considered boson masses and coupling constants.Generally, in both of the plots, there are some regions where Λ 1.4 is modified by the bosonic DM and becomes less than the gray dashed line which shows the maximum confirmed observational value (Λ 1.4 =580).In the left panel, it is seen that for the light bosons, Λ 1.4 is higher than the astrophysical constraint for the whole considered DM fractions which is due to the DM halo formation.For m χ ≥ 200 MeV there are ranges of the DM fraction for which the Λ 1.4 ≲580 bound is satisfied thanks to the reduction in the outermost radius of the star by the DM core configuration.Moreover, in the right panel which shows the impact of the coupling constant, it is seen that lower values of λ are more compatible with the tidal deformability constraint for a wider range of DM fractions inside NSs.This behavior is caused by the fact that for low coupling constants DM mainly resides as a core which decreases the tidal deformability.Regarding both plots, it is that in low DM fractions (F χ ≲ 4%), Λ 1.4 is higher than 580 for all of the applied DM parameters, while it will be decreased by increasing F χ and become less than 580 for massive bosons and low coupling constants.In both Figures 6 and 7, we observe a decrease and increase in the values of the total maximum mass and tidal deformability by changing the DM fraction, this behavior is associated with the formation of a DM core or DM halo, respectively.Comparing the variation of M T max in Fig. (6) with the corresponding radius in Fig. (1), we note that for F χ < 10% where M T max decreases, the radius of the DM fluid is notably smaller than the BM fluid radius.Conversely, for higher F χ where M T max increases, the corresponding DM radius is noticeably larger than the BM one.Moreover, changes in the outermost radius due to the DM core and halo formation (see Fig. 2), significantly affect Λ 1.4 and cause the variation seen in Fig. (7).
In general, the softness/stiffness of DM EoS is varied with DM model parameters, while the inclusion of DM in NSs leads to the softening of the EoS of a DM admixed NS which results in a decrease of the maximum mass and tidal deformability compared to a pure NS.However, towards higher DM fractions a transition occurs from a softer EoS to a stiffer one which corresponds to a shift from DM core to DM halo configuration.Note that the DM core -halo transition is specified by the ratio of R B /R D while the increasing of M T max and Λ 1.4 , which is due to the DM halo formation, occurs at relatively higher DM fractions where the EoS is stiff enough.
Furthermore, as the main objective of this research, by comparing both of the Figures 6  and 7, it is shown that for 4% ≲ F χ ≲ 20%, massive bosons, m χ ≳ 200, and small coupling constants, λ ≲ 2π, present a bosonic DM parameter space for which both M T max and Λ 1.4 are consistent with the astrophysical bounds.This result indicates the modification of NS quantities resulting from DD2 EoS by involving self-interacting bosonic DM.

Conclusion
The occurrence of self-repulsive bosonic DM has been taken into account to investigate the impact of DM on NS features for which the BM is described by the well-known DD2 model.We have seen that for low DM fractions, a dense core is mainly formed inside NSs while by increasing the fraction, bosonic particles distribute as a halo around the NS.It is shown that for massive bosons and low coupling constants, DM will reside as a core for wider ranges F χ , however, for light bosons and large values of λ, an extended DM halo will be formed for the majority of the considered DM fractions.Moreover, for each applied DM model parameter a transition occurs from DM core to DM halo formation by varying the parameter F χ .
Furthermore, we show that a DM core reduces the maximum mass, radius and tidal deformability, but a DM halo will enhance all of them.Therefore, these features demonstrate that including DM in NSs can modify the observable quantities induced from the considered nuclear matter model for the NS interior.Owing to the point that DD2 EoS provides a value for Λ 1.4 which is not consistent with the current observational bound, we examine the bosonic DM parameter space so as to find the regions that make the observable features (M T max and Λ 1.4 ) compatible with their astrophysical limits.We have indicated that for massive bosons (m χ ≳ 200) and small coupling constants (λ ≲ 2π), DM fractions within the range of 4% and 20%, are more favorable and improve the tidal deformability so that the results are consistent with both M max and Λ 1.4 constraints.There are no uncertainties given for the DD2 model and its parameterization which was adjusted at and below saturation density based on the properties of atomic nuclei.A modification of this parameterization to draw new conclusions is out of the scope of the present work and will be deferred to a further study since it requires a careful consideration of parameter correlations.
Finally, by considering the great possibility of accumulation of DM in NSs which alters its observable quantities, future work will account for a larger number of nuclear matter EoSs.Indeed, an overlap of allowed DM parameter spaces can be obtained for the collection of BM EoSs in which all of the observable features are compatible with their constraints.It is noteworthy to say that the analysis performed in this research shows the potential of DM in relaxing the uncertainties in the baryonic EoS space.Furthermore, regarding the advanced GW detectors and X-ray and radio telescopes, conducting similar investigations will help to shed more light on high-density matter in NSs and even the nature of DM.
B , p D , ϵ B and ϵ D are the corresponding pressure and energy density for each of the fluids, respectively, where the B and D indices stand for BM and DM components.Moreover, p = p B + p D denotes the summation of pressure for both BM and DM and M = r 0 4πr 2 ϵ B (r)dr + r 0 4πr 2 ϵ D (r)dr is the mass inside a radius r.Depending on the point where the pressure of one of the fluids gets zero first in the numerical calculations of the two-fluid TOV equations, mainly two possible DM distributions can be found within NSs, called DM core and DM halo configurations.For both cases, the total mass of the object is

Figure 1 .
Figure 1.Variation for the radii of BM (R B ) and DM (R D ) are shown with respect to the DM fraction (F χ ) by solid and dashed lines, respectively.All the radii are related to the maximum mass of DM admixed NSs for various DM model parameters.Note that each colored solid and dashed line corresponds to a specific boson mass (left) or coupling constant (right), as indicated in the legends.In the left panel, different bosonic particle masses are considered as labeled for λ = π, while the right panel is for coupling constants varied from 0.5π to 4π and m χ = 250 MeV.

Figure 2 .
Figure 2. The contour plot shows R B /R D ratio for a scan over the F χ − m χ parameter space considering a fixed coupling constant λ = 0.5π for 1.4M ⊙ DM admixed NSs.The yellow line indicates where R B = R D and the DM-core to DM-halo transition occurs.

Figure 3 .
Figure 3.The dimensionless tidal deformability, denoted as Λ, is presented as a function of the stellar mass for different bosonic particle masses as labeled, considering fixed F χ = 10% and λ = π.The solid black line represents the Λ − M graph for a pure NS without DM.The magenta vertical line signifies the Λ 1.4 constraint derived from the low-spin prior analysis of GW170817, as reported in [87].Taking various DM fractions into account for m χ = 200 MeV and λ = π, Fig. (4) illustrates that for F χ ≲ 20%, the tidal deformability is improved, however, for higher fractions both maximum mass and Λ are not consistent with the observational constraints.Indeed, it is displayed that by raising F χ , the tidal deformability is enhanced which is due to a transition from a DM core to a DM halo formation by varying the amount of DM.Moreover in Fig.(5) we probe the impact of a changing coupling constant on the tidal deformability of DM admixed NSs for F χ = 10% and m χ = 200 MeV.It is seen that by decreasing the coupling constant, the tidal deformability will be reduced and becomes consistent with the observational range while the maximum mass is also in agreement with the 2M ⊙ limit.However, for λ ≳ 2π, we see an enhancement in the tidal deformability which is higher than the astrophysical bound because of the DM halo surrounding the NSs.Therefore, concerning the aforementioned results, in the next section, we will probe the bosonic DM parameter space for which both M T max and Λ 1.4 are in agreement with the astrophysical limits and indicate the capability of bosonic DM in modifying the observable features related to the BM EoS.

Figure 6 .
Figure 6.The change of total maximum mass of DM admixed NSs as a function of F χ are indicated for different DM model parameters.The left panel corresponds to various boson masses and λ = π, while the right panel is related to several coupling constants as labeled and m χ = 250 MeV.The gray dashed line depicts the 2M ⊙ limit for the mass of NSs.Considering the fact that Λ 1.4 of DD2 EoS is not consistent with the observational constraint, Fig. (7) depicts the change of tidal deformability for 1.4M ⊙ DM admixed NSs as a function of F χ .Various boson masses and a fixed coupling constant λ = π are considered in the left panel, however, for the right panel λ is varied for the boson mass m χ = 250 MeV.Generally, in both of the plots, there are some regions where Λ 1.4 is modified by the bosonic DM and becomes less than the gray dashed line which shows the maximum confirmed observational value (Λ 1.4 =580).In the left panel, it is seen that for the light bosons, Λ 1.4 is

Figure 7 .
Figure 7. Variation of Λ 1.4 is shown in terms of DM fraction for different boson masses and coupling constants.In the left panel λ = π is fixed while in the right panel m χ = 250 MeV for all of the cases.The maximum observational limit for the tidal deformability of 1.4M ⊙ NS (580) is indicated by the gray dashed line.
7MeV and L = 55.0MeV.Further, the observable properties of NSs obtained from the DD2 EoS are M max = 2.4M ⊙ , R 1.4 = 13.15 km and Λ 1.4 = 681.61.The large maximum mass, radius and tidal deformability indicate that DD2 EoS is rather stiff.Even though M max and R 1.4 are consistent with astrophysical observations, Λ 1.4 is larger than the observational constraint i.e.Λ 1.4 ≲ 580 from GW170817.In the following sections, we are going to check how the presence of DM could alter the observable features and modify the tidal deformability.