---
title: Dark Matter Admixed Neutron Stars
url: https://www.emergentmind.com/topics/dark-matter-admixed-neutron-stars-dans
type: topic
---

# Dark Matter Admixed Neutron Stars

Searching arXiv for recent and foundational papers on dark matter admixed neutron stars to ground the article.
Dark matter admixed neutron stars (DANSs) are compact stars whose total mass-energy contains both ordinary baryonic or nuclear matter and a non-negligible dark matter component. In the literature summarized here, they are modeled primarily as two-fluid relativistic stars in which the visible neutron-star matter and the dark component are treated as separate perfect or degenerate fluids that interact only through gravity, sharing the same spacetime but generally retaining distinct pressures, energy densities, radii, and in some treatments distinct rotational states. This framework permits configurations ranging from modest dark cores embedded within otherwise ordinary neutron stars to dark-matter-dominated objects with a small normal-matter core inside an extended dark halo, and it has been used to study equilibrium structure, stability, oscillations, cooling thresholds, rotation, tidal response, binary initial data, X-ray observables, and observational bounds on dark matter content [1111.1787].

## 1. Foundational definition and two-fluid formalism

The foundational relativistic description of DANSs treats them as two gravitationally coupled fluids: ordinary nuclear matter (NM), or more generally baryonic matter (BM), and dark matter (DM). In the early general-relativistic treatment, the system is formulated through a master function
\[
\Lambda(n^2,p^2,x^2),
\]
with
\[
n^2=-n_\alpha n^\alpha,\qquad p^2=-p_\alpha p^\alpha,\qquad x^2=-n_\alpha p^\alpha,
\]
where the NM and DM currents are separately conserved. Under the assumption that NM and DM interact only through gravity, the master function is taken to be separable,
\[
\Lambda(n^2,p^2)=\Lambda_{\rm NM}(n^2)+\Lambda_{\rm DM}(p^2),
\]
with no dependence on \(x^2\), so entrainment-type coupling vanishes in that model [1111.1787].

In static spherical symmetry, the spacetime is written as
\[
ds^2 = -e^{\nu(r)}dt^2 + e^{\lambda(r)}dr^2 + r^2(d\theta^2+\sin^2\theta\,d\phi^2),
\]
and the coupled equilibrium equations generalize the Tolman–Oppenheimer–Volkoff system. In one formulation, the structure equations are
\[
A^0_0 p' + B^0_0 n' + \frac{1}{2}(Bn + Ap)\nu' = 0,
\]
\[
C^0_0 p' + A^0_0 n' + \frac{1}{2}(An + Cp)\nu' = 0,
\]
\[
\lambda' = \frac{1-e^\lambda}{r} - 8\pi r e^\lambda \Lambda,
\]
\[
\nu' = -\frac{1-e^\lambda}{r} + 8\pi r e^\lambda \Psi,
\]
while many later treatments adopt the equivalent two-fluid TOV form in which each fluid satisfies its own hydrostatic equation in the common gravitational potential [1111.1787]. A representative form is
\[
\frac{dp_N}{dr}= -\big[\mathcal{E}_N(r)+p_N(r)\big]\frac{d\nu}{dr},\qquad
\frac{dp_D}{dr}= -\big[\mathcal{E}_D(r)+p_D(r)\big]\frac{d\nu}{dr},
\]
with
\[
p(r)=p_N(r)+p_D(r),\qquad
M(r)=\int_0^r 4\pi \tilde r^2\big[\mathcal{E}_N(\tilde r)+\mathcal{E}_D(\tilde r)\big]\,d\tilde r,
\]
so that the fluids are dynamically distinct but gravitationally inseparable [1305.7354].

A persistent feature of this formalism is the existence of two characteristic radii. The baryonic or nuclear surface is defined where the visible fluid pressure vanishes, while the dark surface is defined by the vanishing of the dark pressure. Thus a DANS may be core-dominated, with \(R_D<R_B\), or halo-dominated, with \(R_D>R_B\). Because electromagnetic emission originates at the NM or BM surface, the observed radius may be the visible core radius rather than the full gravitational extent of the object [1305.7354].

## 2. Equations of state and dark-sector microphysics

The DANS literature is distinguished less by the two-fluid gravitational framework, which is broadly shared, than by the diversity of equations of state assumed for the baryonic and dark sectors. For ordinary matter, the cited studies use realistic neutron-star EOSs including APR [1111.1787], relativistic mean-field models such as SLC, SLCd, and IU-FSU [1305.7354], microscopic chiral-interaction matter through the BL EOS derived in Brueckner-Hartree-Fock theory [2405.19251], SLy [2301.03568; 2102.08257], NL3\(\omega\rho\)L55 [2508.18434], DDME2 [2604.10011], unified NL3 from CompOSE [2606.27456], and covariant density functional families spanning PKA1, PKO1–3, PKDD, DD-ME2, DD-MEX, DD-LZ1, NL3, GM1, TM1e, and NL3\(\omega\rho\) [2512.23577].

The dark sector is treated in several distinct ways. A minimal and influential choice is a non-self-annihilating ideal relativistic Fermi gas of fermions with mass \(\sim 1\) GeV [1111.1787]. Closely related studies vary the fermion mass and introduce scalar attraction and vector repulsion through a relativistic mean-field-like Lagrangian,
\[
\mathcal{L}_D= \bar\psi_D\left[\gamma_\mu(i\partial^\mu-g_\mu V^\mu)-(M_D-g_s\phi)\right]\psi_D +\frac{1}{2}(\partial_\mu\phi\,\partial^\mu\phi-m_s^2\phi^2) -\frac{1}{4}D_{\mu\nu}D^{\mu\nu} +\frac{1}{2}m_v^2V_\mu V^\mu ,
\]
with effective potential
\[
V_{\rm eff}(r)=\frac{g_v^2}{4\pi}\frac{e^{-m_v r}}{r} -\frac{g_s^2}{4\pi}\frac{e^{-m_s r}}{r},
\]
so vector repulsion stiffens the DM EOS and scalar attraction softens it [1305.7354].

Other studies use repulsively self-interacting fermionic dark matter in the compact-star parametrization
\[
\varepsilon = \frac{1}{\pi^2}\int_0^{k_F} k^2\sqrt{m_D^2+k^2}\,dk + \left[\left(\frac{1}{3\pi^2}\right)^2 y^2 z^6\right],
\qquad
z=\frac{k_F}{m_D},\quad y=\frac{m_D}{m_I},
\]
with weakly interacting and strongly interacting regimes specified by \(m_I\sim 300\,\mathrm{GeV}\) and \(m_I\sim 0.1\,\mathrm{GeV}\), respectively [2604.10011]. A separate first-principles direction models DM as a confining QCD-like sector based on lattice \(G_2\)-QCD, with the dark matter candidate taken to be the stable three-dark-quark baryon and the EOS interpolated and extrapolated from lattice data [2503.19691].

Bosonic dark matter has also been explored. One rapidly rotating study treats DM as a self-interacting bosonic fluid derived from the complex-scalar Lagrangian
\[
\mathcal{L}=\frac{1}{2}\partial_\nu\phi^*\partial^\nu\phi-\frac{m^2}{2}\phi^*\phi-\frac{\lambda}{4}(\phi^*\phi)^2,
\]
leading in the zero-temperature Bose-condensed limit to
\[
n_{\rm DM} = \frac{\mu(\mu^2-m^2)}{\lambda},\qquad P_{\rm DM} = \frac{(\mu^2-m^2)^2}{4\lambda},
\]
\[
\varepsilon_{\rm DM}=3P_{\rm DM}+2m^2\sqrt{\frac{P_{\rm DM}}{\lambda}}, \qquad h_{\rm DM}=\log\!\left(\frac{\mu}{m}\right)
\]
[2502.17948]. Another study compares self-interacting bosonic and ideal-fermion DM and shows that both can generate halo-induced departures in the mass–radius–tidal response relation at fixed stellar mass [2312.17288].

This breadth of microphysical input is central to the field. It implies that DANS phenomenology is not controlled by a single dark parameter but by particle mass, self-interaction strength, annihilation properties, and, in practice, the central abundance or enthalpy ratio of the two fluids.

## 3. Internal morphology: dark cores, dark halos, and composition parameters

A defining structural distinction in DANS research is the difference between dark-core and dark-halo configurations. Early general-relativistic calculations established that sufficiently DM-rich equilibria can consist of a small NM core of radius of order a few km surrounded by a DM halo of total size \(\sim 10\) km [1111.1787]. In the strongly DM-dominated example with composition parameter
\[
\epsilon = \frac{M_{\rm DM}}{M_{\rm NM}+M_{\rm DM}},
\]
set to \(\epsilon=0.8\), the density profiles show NM concentrated at the center and DM forming the extended halo. For one representative configuration, the visible radius of the NM core is
\[
R_{\rm vis}\approx 0.56\,R,
\]
while the mass enclosed within the NM core is
\[
M_{\rm core}\approx 0.72\,M
\]
[1111.1787].

Later work sharpened the core–halo distinction. A DM halo forms when the DM distribution extends beyond the NM distribution, \(R_D>R_N\), and such halos are found to be more likely for low-mass DM, in the presence of repulsive DM self-interactions, and at lower central energy density [1305.7354]. In contrast, attractive self-interactions concentrate DM and suppress halo formation [1305.7354]. More recent analysis showed that the relative distribution is essentially determined by the ratio of the central enthalpy of the DM component to that of baryonic matter:
\[
h_D^c < h_B^c \Rightarrow \text{dark core},\qquad
h_D^c > h_B^c \Rightarrow \text{dark halo}
\]
[2204.05560]. At the critical configuration \(R_B=R_D\), the critical DM fraction satisfies the analytical estimate
\[
f_\chi^{\rm crit} \simeq 5\times10^{-7} \left(\frac{x}{0.1}\right)^3 \left(\frac{m_\chi}{100\,{\rm MeV}}\right)^4 \left(\frac{1.4}{M_T}\right) \left(\frac{R}{12\,{\rm km}}\right)^3,
\]
implying \(f_\chi^{\rm crit}\propto m_\chi^4\) at fixed self-interaction parameter and \(f_\chi^{\rm crit}\propto y^{-2}\) in the strong self-interaction regime [2204.05560].

Several papers provide useful classification schemes. One crust study defines the halo mass fraction
\[
\delta_M=\frac{M_h}{M_t}
\]
and labels configurations as dark core for \(\delta_M=0\), compact halo for \(0<\delta_M<0.6\), intermediate halo for \(0.6\le\delta_M<0.85\), and diffuse halo for \(0.85\le\delta_M<1\). In its fixed-\(f_\chi=0.05\) sequences, the transition proceeds from diffuse halo for \(m_\chi\lesssim 160\) MeV through compact halo for \(200\text{–}400\) MeV to dark core for \(m_\chi\gtrsim 400\) MeV [2606.27456].

A consistent trend across the literature is that heavier DM tends to be more centrally concentrated, while lighter or more repulsive DM tends to spread outward and form halos. This suggests that the core–halo dichotomy is not a secondary detail but one of the primary observables linking dark-sector microphysics to neutron-star phenomenology.

## 4. Mass–radius relations, compactness, and maximum-mass systematics

In ordinary neutron-star theory, a fixed EOS produces a one-parameter mass–radius sequence. A central result of DANS research is that this property is lost once dark matter is admitted. Because the equilibrium depends on two EOSs and on parameters such as DM fraction and DM self-interaction strength, DANSs populate a family of mass–radius curves rather than a unique sequence [1305.7354]. This “spread of mass-radius relationships” implies that a measured mass and radius cannot be mapped to a unique nuclear EOS without additional assumptions about DM content [1305.7354].

The earliest detailed DANS study found that increasing the DM fraction shifts the mass–radius relation away from that of ordinary neutron stars, decreases the maximum stable mass, and decreases the stellar radius. For the maximum stable-mass configuration, increasing \(\epsilon\) from 0 to 0.2 reduces \(M\) by about 35% and \(R\) by about 9% [1111.1787]. A later two-fluid TOV analysis confirmed that for fixed DM particle mass, increasing the DM fraction lowers the maximum stable mass, reduces the stellar radius, and increases the spread of possible mass–radius curves [1305.7354].

The effect of DM mass is more nuanced. In one study, heavy DM candidates with mass \(\gtrsim 2\) GeV have little effect on the mass–radius relation, whereas light DM of order \(500\) MeV produces major changes in both maximum mass and radius [1305.7354]. Another study using a realistic chiral nuclear EOS finds opposite trends depending on DM mass: for \(m_\chi=1\,\mathrm{GeV}\), DM reduces the maximum gravitational mass relative to the ordinary-star case, but for \(m_\chi=0.1\,\mathrm{GeV}\), DANSs form with a DM halo and the total gravitational maximum mass increases relative to the ordinary-star case [2405.19251]. The same work identifies a transition at \(m_\chi\simeq 0.188\,\mathrm{GeV}\): for \(m_\chi \lesssim 0.188\,\mathrm{GeV}\), the maximum mass remains compatible with observed \(\sim 2\,M_\odot\) neutron stars for any DM fraction, whereas for larger masses only restricted DM fractions are permitted [2405.19251].

A related scaling appears in the self-gravitating pure-DM limit. For fermionic dark stars, one study quotes
\[
M_{\max} = 0.627\,M_\odot \left(\frac{1\,{\rm GeV}}{M_D}\right)^2,\qquad
R_D = 8.115\,{\rm km}\left(\frac{1\,{\rm GeV}}{M_D}\right)^2
\]
[1305.7354]. A 2026 extension to extremely light DM emphasizes the corresponding general scaling
\[
M_{\max} \approx 0.627\left(\frac{\mathrm{GeV}}{m_D}\right)^2 M_\odot,
\]
or, including self-interaction,
\[
M_D^{\max} = \left(0.269\,\frac{m_D}{m_I}+0.627\right) \left(\frac{\mathrm{GeV}}{m_D}\right)^2 M_\odot,
\]
and argues that for \(m_D<10^{-1}\) GeV the DANS becomes DM-dominated, with the NM core reduced to an embedded seed inside an enormous DM halo [2604.10011]. For \(m_D\sim 10^{-5}\,\mathrm{GeV}\), the paper reports
\[
M_D \sim 6.58\times 10^9\,M_\odot,\qquad
R_D \sim 7.86\times 10^{10}\,\mathrm{km},
\]
while the NM core remains near
\[
M_N \approx 2.48\,M_\odot,\qquad R_N \approx 11.7\,\mathrm{km}
\]
[2604.10011]. This suggests that the DANS concept can interpolate continuously between modestly admixed neutron stars and supermassive dark objects with embedded neutron-star cores.

Not all approaches use a strict two-fluid separation. Some mean-field studies incorporate DM directly into the core EOS through Higgs exchange or via fixed DM Fermi momentum. In those models, larger DM masses and larger DM Fermi momenta generally produce more compact stars with reduced radii, lower maximum masses, and smaller tidal deformabilities [2408.15312; 2411.17105]. A specific case study interprets XTE J1814-338, with inferred
\[
M = 1.21\pm 0.05\,M_\odot,\qquad R = 7.0\pm 0.4\ \text{km},
\]
as a DANS requiring a DM Fermi momentum \(k_f=0.08\) GeV, yielding \(R_{1.4}=7.10\) km and \(M_{\max}=1.55\,M_\odot\) in the nucleonic version of the model [2411.17105]. Because that work does not present a dedicated radial-stability analysis, its ultracompact configurations are best interpreted as phenomenological TOV solutions rather than fully established stable equilibria [2411.17105].

## 5. Stability, oscillations, and crustal seismology

The stability of DANSs has been examined through several complementary methods. The original relativistic two-fluid study solved the radial oscillation eigenvalue problem and found that for the DM-dominated sequence with \(\epsilon=0.8\), the fundamental mode satisfies \(\omega_0^2>0\) along the stable branch, crosses zero at the maximum-mass configuration, and becomes negative beyond it. Thus the onset of instability coincides with the maximum-mass point, preserving the standard one-fluid intuition in the two-fluid setting [1111.1787].

A later comprehensive analysis emphasized that in a two-fluid star the stability boundary is not a single point but a critical curve in the central-pressure plane \((p^c_{\rm om},p^c_{\rm dm})\). It established stability in two ways: by solving the harmonic radial perturbation problem with \(\xi_i(t,r)=\xi_i(r)e^{i\omega t}\), and by applying a turning-point criterion along directions tangent to contours of conserved particle numbers. Along the stability boundary there exists a direction \(\mathbf p\) such that
\[
\frac{dM}{d\mathbf p} = \frac{dN_{\rm om}}{d\mathbf p} = \frac{dN_{\rm dm}}{d\mathbf p} =0,
\]
and the two methods agree numerically [2102.08257]. A central conclusion is that there exist stable regions in parameter space where either the ordinary-matter or dark-matter central density exceeds its corresponding one-fluid critical value, so one-fluid stability criteria are insufficient for DANSs [2102.08257; 2405.19251].

Nonlinear stability has also been verified dynamically. In the first dynamical evolution of DANSs, the equations of motion for multiple perfect fluids were cast in conservation form and evolved using finite volume and high-resolution shock-capturing methods. Stable static DANS remained stable nonlinearly, while unstable ones collapsed to black holes. The same simulations recovered radial oscillation frequencies by Fourier transforming central-pressure time series and showed excellent agreement with linear pulsation theory [2201.02274]. That work also demonstrated a dynamical formation scenario in which a preexisting neutron star surrounded by a Gaussian dark-matter cloud evolves into a stable mixed star, including a case where a dark-matter core forms by \(\bar t\sim 1250\) [2201.02274].

The crust has recently been incorporated explicitly. In DANS models with unified NL3 baryonic matter and repulsively self-interacting fermionic DM, the crust thickness
\[
\Delta R = R_b - R_{\rm cc}
\]
is always smaller than in the corresponding pure neutron star at fixed total mass and fixed DM fraction \(f_\chi=0.05\). For \(M_t=1.4\,M_\odot\), the maximum reduction is about 12% near \(m_\chi\sim 1\) GeV; for \(M_t=2.0\,M_\odot\), the maximum reduction is about 16% [2606.27456]. The effect is negligible for diffuse halos, where little DM lies inside the baryonic surface, and strongest for compact halos and dark cores [2606.27456]. The same paper derives the approximation
\[
\frac{\Delta R}{R_b}\approx \frac{\alpha(1-2C_b)}{C_b+\alpha(1-2C_b)},\qquad \alpha=\ln\!\left(\frac{\mu_{\rm cc}}{\mu_s}\right),
\]
with \(\alpha=0.0256\) for NL3, and finds that the rational form reproduces numerical crust thicknesses typically below 0.5% error and at worst about 1% [2606.27456].

This crustal compression raises torsional mode frequencies. In the relativistic Cowling approximation, the torsional eigenvalue problem is
\[
Y_\ell'' +\left(\frac{4}{r}+\Phi'-\Lambda'+\frac{\mu'}{\mu}\right)Y_\ell' +e^{2\Lambda}\left[ \frac{\omega^2}{v_s^2}e^{-2\Phi} -\frac{(\ell+2)(\ell-1)}{r^2} \right]Y_\ell=0,
\]
with boundary conditions at the crust-core interface and surface [2606.27456]. For \(M_t=1.4\,M_\odot\), the maximum frequency shifts relative to pure neutron stars are about 6.4% for the quadrupole fundamental \({}_2 t_0\) and about 16.4% for the first overtone \({}_2 t_1\) as \(m_\chi\to 1\) GeV [2606.27456]. Because electron screening softens the crust and lowers frequencies, while DM thins the crust and raises them, the two effects can be degenerate; the paper finds that this degeneracy can be broken in some regions of parameter space, especially for massive stars and the first overtone [2606.27456].

## 6. Rotation, tidal response, binaries, and multimessenger observables

Rotation exposes the multi-fluid nature of DANSs particularly clearly. In the slow-rotation Hartle formalism generalized to multiple fluids, each fluid has its own angular velocity \(\Omega_x\), four-velocity
\[
u_x^\mu = u_x^t(1,0,0,\Omega_x),
\]
and moment of inertia \(I_x=J_x/\Omega_x\), while the total angular momentum is \(J=\sum_x J_x\) [2311.07714]. The first-order frame-dragging equation is sourced by the sum of all fluids, but each fluid contributes separately to the angular-momentum budget [2311.07714]. A key observational point is that electromagnetic measurements probe ordinary matter alone, whereas gravitational observables such as the tidal deformability depend on the total mass distribution. Accordingly, the observationally relevant dimensionless quantities are
\[
\overline{I}_{\rm om}=\frac{I_{\rm om}}{M_*^3},\qquad
\overline{\lambda}=\frac{\lambda}{M_*^5},\qquad
\overline{Q}_{\rm om}=\frac{QM_*}{J_{\rm om}^2},
\]
not their total-fluid analogues [2311.07714]. Using these definitions, the paper shows that DANS can deviate significantly from standard single-fluid I–Love–Q relations [2311.07714].

Rapid rotation has been studied in two distinct regimes. One axisymmetric two-fluid calculation extends the RNS code to rapidly rotating DANS in which BM is spun up by accretion and rotates rigidly, while DM remains torque-free and differentially rotates through frame dragging:
\[
v_D=0 \implies \Omega_D=\omega
\]
[2508.18434]. For three representative halo types with non-rotating DM fraction \(f_{\chi,\mathrm{TOV}}=0.05\), corresponding at \(1.4\,M_\odot\) to \(R_D\approx 19.4\) km, \(48.5\) km, and \(530\) km, rapid BM rotation reduces DM halo sizes when central densities are held fixed. At \(\bar r_{B,\mathrm{ratio}}=0.9\), the baryonic equatorial radius increases by about 5% and the polar radius decreases by about 4%, while the diffuse halo radius can shrink by about 80% and the cloud mass outside the baryonic surface by about 97%; the corresponding reductions are about 20% and 39% for the intermediate halo and about 5% and 37% for the compact halo [2508.18434]. The same study emphasizes that local TOV-like component masses and global Komar-type component masses differ and should not be interpreted as separately measurable gravitational masses [2508.18434].

A complementary rapid-rotation study allows BM and DM to rotate independently, including co-rotating and counter-rotating states, in models with self-interacting bosonic DM [2502.17948]. There the defining feature is independent fluid angular velocities,
\[
\Omega_{\rm BM}\neq \Omega_{\rm DM},
\]
with even opposite signs permitted [2502.17948]. In the baseline case where BM rotates rapidly but DM is non-rotating, pure baryonic stars show the expected rotational increase of maximum mass by about 20% and radius by about 40%, while DM core configurations reduce both maximum mass and moment of inertia relative to no-DM stars [2502.17948]. Turning on DM rotation reveals strong qualitative differences between co-rotation and counter-rotation. In one example the ratio of DM polar to equatorial radius drops from about 0.51 in co-rotation to about 0.34 in counter-rotation, and counter-rotating sequences become numerically difficult before the last stable configuration because of extreme deformation [2502.17948].

Tidal response has likewise emerged as a major diagnostic. One study proposes a new criterion: at fixed stellar mass, DANS can exhibit a negative correlation between normal-matter radius \(R_N\) and tidal deformability \(\Lambda\), so that a star with smaller visible radius has larger \(\Lambda\). This does not occur for ordinary neutron stars in that analysis and arises because \(\Lambda\) depends on the larger of \(R_N\) and \(R_D\), while the observed radius tracks only \(R_N\) [2312.17288]. For APR3 with self-interacting bosonic DM of mass \(m_{\chi,b}=200\) MeV, the paper reports that \(\Lambda_{1.4}\) can reach about 8000 for large DM fraction \(F_X\sim 0.8\), compared with \(\Lambda_{1.4}\sim 400\) in ordinary APR3 stars [2312.17288].

Binary DANS systems have also begun to be modeled self-consistently. Using an adapted SGRID code, the first constraint-solved quasi-equilibrium initial data for equal-mass binary neutron stars each containing dark matter were constructed in both dark-core and dark-halo scenarios [2301.03568]. The stars are treated as two non-interacting perfect fluids minimally coupled to gravity, with irrotational flow \(w^{(s)i}=0\), and the data solve the coupled XCTS equations plus fluid equilibrium conditions [2301.03568]. DM modifies tidal response in a morphology-dependent way: halo configurations are much more easily distorted, and the DM component in the halo case reaches a mass-shedding diagnostic \(\chi\approx 0.9\) already at relatively large separation; mass shedding occurs at separations between about \(28 M_\odot\) and \(29 M_\odot\), whereas in the dark-core case the baryonic fluid reaches mass shedding only around \(16 M_\odot\) in the representative sequence [2301.03568].

## 7. Observational signatures, constraints, and open interpretations

The principal observational motivation for DANSs is that dark matter can decouple the gravitational and electromagnetic appearance of a compact star. A DANS may have a total mass in the neutron-star range and a total radius of order 10 km, yet only a smaller NM-emitting core contributes to thermal radiation. The original DANS paper therefore emphasized that such objects may be observed as extraordinarily small neutron stars incompatible with realistic nuclear matter models [1111.1787]. A later two-fluid study illustrated this with a comparison at fixed total mass \(1.35M_\odot\): the DANS had a smaller NM core radius of about 9 km inside a total radius \(12.76\) km and a larger redshift, roughly \(z=0.340\) versus \(z=0.206\) for the ordinary neutron star [1305.7354].

Moment of inertia offers an additional discriminator. In the non-interacting two-fluid limit of the original relativistic model, one can write
\[
I_{\rm NM} = \frac{J_{\rm NM}}{\Omega_{\rm NM}},\qquad
I_{\rm DM} = \frac{J_{\rm DM}}{\Omega_{\rm DM}},\qquad
I=I_{\rm NM}+I_{\rm DM},
\]
and define the scaled quantity
\[
\tilde I = \frac{I}{MR^2}
\]
[1111.1787]. Ordinary neutron stars lie near the Bejger–Haensel universal fit in terms of compactness \(z=(M/M_\odot)(\mathrm{km}/R)\), but DANS do not; for DM-dominated stars with \(\epsilon=0.8\), \(\tilde I\) is significantly smaller than expected for an ordinary neutron star of the same compactness [1111.1787]. This line of argument is strengthened in the slow-rotation I–Love–Q analysis, which shows that once observationally relevant BM-only moments of inertia are used, the deviations from single-fluid relations are larger still [2311.07714].

X-ray pulse-profile modeling is particularly sensitive to dark halos because the metric exterior to the baryonic surface is not vacuum if DM extends beyond it. A focused study of halo/core formation and pulse profiles found that the peak flux deviation is controlled by the compactness-like ratio \(M_{\rm halo}/R_D\), with
\[
\log\left(\frac{|\Delta F_{\rm peak}|}{F_{\rm peak}}\right) = 0.998\log\left(\frac{M_{\rm halo}}{R_D}\right)+0.854,
\]
and \(R^2=0.995\) [2204.05560]. For a compact halo with \(M_{\rm halo}/R_D\sim 0.01\), the peak flux deviation can reach \(\sim 10\%\), while nearly halo-free systems with \(M_{\rm halo}/R_D\le 10^{-3}\) should distort the flux by less than \(\sim 1\%\) [2204.05560]. Bayesian analysis using NICER posteriors for PSR J0030+0451 and PSR J0740+6620 then yielded a 90% upper limit \(m_\chi\lesssim 1.5\,\mathrm{GeV}\) in the nearly halo-free case, while in the pure dark-core case the posterior peaks around \(m_\chi\sim 0.6\,\mathrm{GeV}\); the self-interaction strength \(y\) remained unconstrained [2204.05560].

Several papers convert recent multimessenger observations into bounds on dark content. With twelve covariant density functional EOSs and self-interacting fermionic DM fixed at \(m_\chi=1000\) MeV and \(y=0.1\), one study finds a strong linear correlation
\[
f_\chi^{\max}=0.22\,M_{\mathrm{NS}^{\max}}-0.44
\]
between the maximum allowed DM fraction and the maximum mass of the corresponding pure neutron star, with Pearson coefficient \(r=0.98\) [2512.23577]. Propagating the observational maximum-mass distribution then yields
\[
M_{\chi}^{\max}=0.150^{+0.070}_{-0.051}\,M_\odot
\]
at 68% confidence [2512.23577]. A more targeted study of PSR J0740+6620 argues that the simultaneous mass–radius measurement reduces the uncertainty in the DM central density ratio \(\varepsilon_D^c/\varepsilon_N^c\) by more than 50% compared with using mass or radius separately, and concludes that the total DM fraction should be smaller than 2% when constrained by the observed neutron-star maximum mass alone and smaller than 0.3% with the simultaneous mass–radius measurement [2509.07767].

Observational interpretation remains model dependent. Some studies invoke DANSs to explain unusually compact objects such as EXO 1745-248, 4U 1608-52, 4U 1820-30 [1612.02804], HESS J1731-347 [2411.17105], or XTE J1814-338 [2411.17105]. In each case the mechanism is the same: dark matter softens or compresses the visible star, allowing smaller apparent radii than standard hadronic EOSs typically admit. A plausible implication is that DANSs are most difficult to distinguish from ordinary neutron stars when the dark fraction is small and core-like, and easiest to distinguish when the star carries a sizable halo, an anomalously small visible radius, or a moment of inertia or tidal response incompatible with one-fluid expectations.

The main unresolved issue across this literature is not whether mathematically consistent DANS equilibria exist—they do in many formulations—but how much dark matter real neutron stars can accumulate and which microphysical model is realized in nature. The review literature therefore treats DANSs not as a single astrophysical object class with settled parameters, but as a broad two-fluid framework connecting dark-sector microphysics to neutron-star structure, stability, and multimessenger phenomenology [2412.09381].

Source: https://www.emergentmind.com/topics/dark-matter-admixed-neutron-stars-dans