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

# Dark Matter-Admixed Neutron Stars

Dark matter-admixed neutron stars, commonly abbreviated DMANSs or DANS, are compact stars in which ordinary neutron-star matter coexists with a dark component. In the literature, they are modeled either as genuinely two-fluid relativistic stars, with ordinary matter and dark matter coupled only through gravity, or as effectively single-fluid objects under stronger assumptions such as additive thermodynamics or chemical equilibrium between visible and dark sectors. Across these formulations, the presence of dark matter can modify the mass-radius relation, maximum mass, tidal deformability, moments of inertia, oscillation spectra, internal composition, and cooling thresholds, but the sign and magnitude of those changes depend strongly on the dark equation of state, particle mass, self-interaction strength, and spatial morphology of the dark component [1111.1787] [1208.3722] [2506.18892] [2508.19382].

## 1. Conceptual frameworks and terminology

Across the current literature, three recurrent modeling classes appear. In the most common class, ordinary matter and dark matter are treated as separate perfect fluids with separate hydrostatic equations and, in general, separate radii \(R_{\rm BM}\) and \(R_{\rm DM}\). A second class uses a one-fluid additive prescription, in which the total pressure and total energy density are taken as sums of baryonic and dark contributions. A third class treats the star as an effectively single fluid because the dark component is assumed to be in chemical equilibrium with neutrons and to follow the hadronic density profile [1111.1787] [1208.3722] [2506.18892].

| Modeling class | Defining assumption | Representative works |
|---|---|---|
| Two-fluid relativistic star | Separate BM and DM fluids, gravitational coupling only | [1111.1787], [2102.08257], [2508.19382] |
| Additive one-fluid model | \(P=P_N+P_\chi\), \(\mathcal{E}=\mathcal{E}_N+\mathcal{E}_\chi\) | [1208.3722], [2411.17105] |
| Chemically equilibrated mixed fluid | \(\mu_N=\mu_\chi\); no separate DM TOV sector | [2506.18892] |

A further standard distinction is morphological. In two-fluid models, one speaks of a **DM core** when \(R_{\rm DM}<R_{\rm BM}\), a **DM halo** when \(R_{\rm DM}>R_{\rm BM}\), and a **coextensive distribution** when \(R_{\rm DM}=R_{\rm BM}\). These are not merely descriptive labels: they determine which radius enters electromagnetic observables, which radius enters tidal response, and whether the dark component tends to increase or decrease compactness-related quantities [2508.19382] [2607.03840].

The early two-fluid study of Leung, Chu, and Lin identified a new class of stable compact stars consisting of a small normal-matter core with radius of a few km embedded in a ten-kilometer-sized dark-matter halo. In one representative DM-dominated configuration, the visible normal-matter radius was \(R_{\rm NM}=0.56\,R\), while the mass enclosed inside that visible core was \(0.72\,M\), implying a larger surface redshift than in an ordinary neutron star with the same total \(M\) and \(R\) [1111.1787].

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

The ordinary-matter sector has been modeled with a wide range of nuclear equations of state. These include APR and SLy-like descriptions, relativistic mean-field and quantum hadrodynamics models, the microscopic Brueckner-Hartree-Fock framework with hyperons, and a state-of-the-art BL equation of state for \(\beta\)-stable nuclear matter obtained in Brueckner-Hartree-Fock using two-body and three-body interactions derived from chiral effective field theory. Some recent work goes further and samples the hadronic EOS agnostically, anchoring it at low density by chiral EFT and at high density by perturbative QCD through a speed-of-sound construction [2405.19251] [2506.18892] [2607.03840].

The dark sector is even more heterogeneous. Representative choices include ideal relativistic fermion gases, self-interacting fermions parameterized by an interaction scale \(m_{\mathcal I}\), quartic self-interacting bosonic condensates, neutralino-like Higgs-portal fermions, effective dark equations of state inferred from galaxy rotation curves, and strongly interacting confining dark matter based on one-flavor \(G_2\)-QCD with lattice input [1208.3722] [2508.19382] [2411.17105] [1612.02804] [2503.19691].

For scalar-field dark matter at zero temperature in a Bose-Einstein-condensed phase, one representative EOS used in gravitational-wave inference is
\[
p= \frac{m_\chi^4}{9\lambda_\chi} \left( \sqrt{1+\frac{3\lambda_\chi}{m_\chi^4}\epsilon}-1 \right)^2,
\]
where \(m_\chi\) is the particle mass and \(\lambda_\chi\) the quartic self-coupling. For ideal fermionic dark matter in a realistic two-fluid study, the EOS is written parametrically as
\[
\varepsilon(x)=\alpha(m_\chi,\nu)\,\chi(x),\qquad
p(x)=\alpha(m_\chi,\nu)\,\phi(x),\qquad
n(x)=n_0(m_\chi,\nu)\,x^3,
\]
with \(x=\hbar k_F/(m_\chi c)\) and \(\nu=2\) for spin-\(\tfrac12\) particles [2508.19382] [2405.19251].

Despite this diversity, one modeling assumption recurs: in most two-fluid studies, ordinary matter and dark matter are taken to interact only gravitationally. Notable exceptions exist. In the Higgs-portal single-fluid models used for some ultracompact-star and radial-oscillation studies, the dark fermion is heavy, the coupling is mediated by the Standard Model Higgs, and the DM content is parameterized phenomenologically by a dark Fermi momentum \(k\) or \(k_f^{\rm DM}\) rather than derived from a capture history [2411.17105] [2211.12808]. In the neutron-dark-decay scenario, by contrast, the dark fermion abundance is fixed by chemical equilibrium \(\mu_N=\mu_\chi\), and the dark density profile follows the hadronic one closely enough that the star is treated as an effectively single fluid [2506.18892].

## 3. Hydrostatic structure, mass fractions, and core-halo morphology

In the two-fluid formalism, the total stress tensor is the sum of the ordinary and dark perfect-fluid tensors, and each sector is separately conserved. A standard form of the hydrostatic equations is
\[
\frac{dp_I}{dr} = -(p_I+\epsilon_I)\, \frac{m(r)+4\pi r^3 p}{r[r-2m(r)]},
\qquad
\frac{dm}{dr}=4\pi r^2 \epsilon,
\]
with
\[
p=p_{\rm DM}+p_{\rm BM},\qquad \epsilon=\epsilon_{\rm DM}+\epsilon_{\rm BM},\qquad
m(r)=m_{\rm DM}(r)+m_{\rm BM}(r).
\]
Integration starts from chosen central energy densities for the two fluids and proceeds outward until one pressure vanishes, then until the second vanishes, thereby defining \(R_{\rm BM}\) and \(R_{\rm DM}\) separately [2508.19382].

Different papers adopt different DM-fraction variables. In gravitational-wave parameter estimation one often uses
\[
F_\chi=\frac{M_{\rm DM}}{M_T},
\]
with \(M_T=M_{\rm BM}+M_{\rm DM}\). In the chiral-EOS and rapid-rotation literature the notation
\[
f_\chi=\frac{M_{\rm DM}}{M_{\rm tot}}
\]
is common, whereas in the chemically equilibrated neutron-dark-decay model the global fraction is defined energetically,
\[
f_{\rm DM} = \frac{\int_V \epsilon_\chi\, dV}{\int_V \epsilon\, dV}.
\]
These are not interchangeable without care, because they refer to different model classes and different effective fluids [2508.19382] [2405.19251] [2506.18892] [2502.17948].

Morphology is controlled primarily by the dark EOS stiffness and the dark particle mass. In the agnostic two-fluid framework, light fermionic DM tends to form extended halos, while heavy DM tends to form compact cores; specifically, \(m_D=0.2~\mathrm{GeV}\) yields only halo configurations and \(m_D=1.1~\mathrm{GeV}\) yields only core configurations in the sampled model space. The same core-halo dichotomy appears in particle-specific studies: \(m_\chi=1\,\mathrm{GeV}\) typically produces a compact core, whereas \(m_\chi=0.1\,\mathrm{GeV}\) produces a halo [2607.03840] [2405.19251].

The one-fluid additive literature uses a different language but often reaches an analogous phenomenology. When heavy, weakly pressurized dark matter is added to the total EOS, the resulting object can become so compact that the authors themselves remark that it looks “more like DM-stars rather than NSs.” Conversely, when very light or strongly pressure-supported dark matter dominates the structure, ordinary distinctions between nucleonic and hyperonic baryonic sectors can become negligible in the mass-radius sequence [1208.3722].

## 4. Static properties: mass-radius relations, maximum mass, and composition

No single monotonic rule covers all DMANS models. In the additive one-fluid analysis of non-self-annihilating self-interacting fermions, dark matter generally softens the EOS more strongly than hyperons and can reduce the maximum mass severely. For the interpretation of PSR J1614-2230 as a hyperon star within that framework, compatibility with the measured \(1.97\pm0.04\,M_\odot\) requires
\[
m_\chi \lesssim 0.64~\mathrm{GeV}
\]
for strongly interacting DM and
\[
m_\chi \lesssim 0.16~\mathrm{GeV}
\]
for weakly interacting DM. The same work also showed that an ambient DM halo around the star would contribute only about \(10^{-24}\,M_\odot\) to the inferred mass, so the measured mass must be attributed to the compact object itself rather than to ordinary Galactic halo material [1208.3722].

In the realistic two-fluid study based on the BL chiral-BHF EOS, the sign of the effect changes with \(m_\chi\). For \(m_\chi=1\,\mathrm{GeV}\), the stellar gravitational maximum mass decreases as the DM fraction increases, whereas for \(m_\chi=0.1\,\mathrm{GeV}\), \(M_{\max}\) increases relative to ordinary neutron stars. For the \(1\,\mathrm{GeV}\) case, the paper gives
\[
M_{\max}=2.07\,M_\odot \quad (f_\chi=0),
\qquad
1.99\,M_\odot \quad (f_\chi=3.10\%),
\qquad
1.94\,M_\odot \quad (f_\chi=5.41\%),
\]
and finds that the direct-URCA threshold mass drops from
\[
M_{tot}^{durca}=0.96\,M_\odot
\]
to \(0.87\,M_\odot\) and \(0.82\,M_\odot\) at the same two fractions. In addition, for a \(1.4\,M_\odot\) star, the observable radius \(R_{OM}\) decreases by about \(0.6\,\mathrm{km}\), or about \(5\%\), at \(f_\chi=6\%\) [2405.19251].

The agnostic two-fluid study reaches an important negative conclusion: dark matter does not generically compactify the star. Light halo-forming DM raises the effective radius relevant for tides and therefore increases \(\Lambda\), whereas heavy core-forming DM lowers it. This shifts the dominant observational constraint from gravitational-wave tidal deformability in the light-halo regime to NICER mass-radius measurements in the heavy-core regime [2607.03840].

The same theme appears in candidate-object modeling. In the heavy, weakly interacting neutralino-like one-fluid construction proposed for XTE J1814-338, reproducing
\[
M = 1.21 \pm 0.05\,M_\odot,\qquad R = 7.0 \pm 0.4~{\rm km}
\]
requires a DM Fermi momentum \(k=0.08~\mathrm{GeV}\), yielding
\[
I_{1.2}\approx (0.56\!-\!0.58)\times10^{45}\ {\rm g\,cm^2},
\qquad
\Lambda_{1.2}\approx 50\!-\!57.
\]
The same paper stresses, however, that such unusually small-radius stars are not unique evidence for dark matter admixture, since hybrid stars, strange quark stars, and other exotic interpretations remain open [2411.17105].

## 5. Oscillations and stability

Radial stability has been studied most rigorously in the two-fluid framework. A broad analysis using a multi-fluid TOV system with either free-fermion or mirror dark matter showed that the stability boundary can be located both by solving the harmonic radial perturbation problem and by a static turning-point/tangency criterion; the two methods agree. One noteworthy result is that stable two-fluid stars can exist in regions where both central pressures exceed the corresponding one-fluid critical pressures, so the stability landscape is not reducible to ordinary neutron-star intuition [2102.08257].

In the single-fluid Higgs-portal model, dark matter softens the EOS but raises radial oscillation frequencies because the resulting stars are more compact. For a canonical \(1.4\,M_\odot\) star, the fundamental radial frequency increases from
\[
2.8450~{\rm kHz}
\]
at \(k_f^{\rm DM}=0\) to
\[
3.1841~{\rm kHz}
\]
at \(k_f^{\rm DM}=0.03~{\rm GeV}\) and
\[
3.7404~{\rm kHz}
\]
at \(k_f^{\rm DM}=0.05~{\rm GeV}\). In the same sequence, the critical central density shifts from
\[
2.0022 \times 10^{15}~{\rm g/cm^3}
\]
to \(2.2022 \times 10^{15}~{\rm g/cm^3}\) and \(2.9177 \times 10^{15}~{\rm g/cm^3}\), while the usual criterion \(\partial M/\partial \rho_c>0\) remains valid along the stable branch [2211.12808].

A different oscillation signature arises in the chemically equilibrated neutron-dark-decay scenario, where the star is treated as compositionally stratified single-fluid matter. There the buoyancy enhancement produced by the dark component raises the \(g\)-mode spectrum, and the shift depends predominantly on the global dark fraction rather than on hadronic effective-mass variations. The resulting quasi-universal fits are
\[
\Delta \nu_{g_1} = 14.23(\pm 0.02)\times f_{\rm DM}[\%],
\qquad
\Delta \nu_{g_2} = 10.18(\pm 0.02)\times f_{\rm DM}[\%].
\]
The same study argues that observing \(g_1\) frequencies as high as \(\sim 500\) Hz for \(M\sim 1\,M_\odot\) would be a “tell-tale signature” of dark matter, because ordinary hadronic stars or threshold-driven hyperon/quark phases would not generally reach such high low-mass frequencies [2506.18892].

Taken together, these results suggest that DMANS stability and seismology are highly sensitive to which degree of freedom is dynamical: two-fluid radial modes probe coupled hydrostatics, whereas the neutron-dark-decay scenario probes composition-gradient buoyancy in an effectively single fluid.

## 6. Rotation, tides, binaries, and waveform modeling

Rotation exposes a conceptual asymmetry between electromagnetic and gravitational observables. In a slow-rotation two-fluid treatment, the star possesses separate ordinary-matter and dark-matter angular momenta and moments of inertia, while the total exterior spacetime and quadrupole are sourced by both fluids. Because electromagnetic observations directly probe ordinary matter rather than dark matter, the standard single-fluid I-Love-Q relations can deviate significantly in DM-admixed stars once the sector actually measured is distinguished from the sector that sources the external field [2311.07714].

In the rapid-rotation regime, a fully relativistic two-fluid extension of the RNS framework shows that independent BM and DM rotation changes the equilibrium problem qualitatively. For a \(1\,\mathrm{GeV}\) bosonic DM core, allowing the DM itself to rotate raises the maximum mass by about \(4.5\%\) relative to the nonrotating-DM case, whereas in a nominal halo case the corresponding change is a decrease of about \(2.5\%\). The same study finds much stronger deformation for counter-rotating DM than for co-rotating DM: in representative \(M_T=2\,M_\odot\), \(f_\chi=5\%\) core configurations,
\[
\frac{R_{p,\mathrm{DM}}}{R_{e,\mathrm{DM}}}\approx 0.51
\]
for co-rotation and
\[
\frac{R_{p,\mathrm{DM}}}{R_{e,\mathrm{DM}}}\approx 0.34
\]
for counter-rotation. Rapid rotation can also blur the static core-halo distinction, producing mixed geometries in which BM and DM alternate in radial dominance between pole and equator [2502.17948].

Binary dynamics amplify these structural differences. The first consistent quasi-equilibrium initial data for DM-admixed binaries showed that DM-core stars are less deformable than ordinary neutron stars of the same total rest mass, whereas DM-halo stars can undergo large halo deformation and approach mass shedding much earlier. In a representative halo case, the DM surfaces touched and mass shedding was inferred to begin between \(28\) and \(29\,M_\odot\) separation, while in a representative DM-core case the baryonic matter approached mass shedding only at much smaller separations [2301.03568].

For gravitational waves, the central observable is the dimensionless tidal deformability,
\[
\Lambda=\frac{2}{3}k_2 C^{-5},
\qquad C=\frac{M}{R}.
\]
Recent infrastructure work has translated this stellar physics into waveform-generation tools. The proposed package `Darksuite` constructs a bank of two-fluid TOV solutions, tabulates \((C,f,k_2)\), interpolates them with linear splines, and injects the resulting \(\Lambda\) values into `LALSuite` tidal phase corrections such as the NRTidal family, thereby enabling waveform generation for binaries containing one or more DM-admixed stars [2508.19382] [2507.22415].

## 7. Observational constraints, candidate objects, and unresolved issues

The first direct attempt to constrain DMANSs with real gravitational-wave data reanalyzed GW230529, GW200115, GW200105, and GW190814 as NS-BH systems. For GW200105, the strongest bound among the standard events was obtained:
\[
F_\chi < 6.1\% \quad (95\%~{\rm CL}),\qquad m_\chi>435~{\rm MeV}
\]
with SLy4, and
\[
F_\chi < 9.52\%,\qquad m_\chi>403~{\rm MeV}
\]
with APR4. GW230529 and GW200115 admitted much weaker bounds but still favored DM-core configurations. GW190814 was exceptional: if one insists on soft SLy4 or APR4 hadronic EOSs, the secondary can be interpreted as a DM-halo star with
\[
F_\chi>78.7\%,\qquad m_\chi=216^{+30}_{-20}\ {\rm MeV}
\]
or
\[
F_\chi>75.6\%,\qquad m_\chi=214^{+40}_{-20}\ {\rm MeV},
\]
whereas with the stiffer MPA1 EOS the same event moves to a DM-core-like regime with
\[
F_\chi<2.2\%,\qquad m_\chi>366\ {\rm MeV}.
\]
This is a clear EOS-DM degeneracy rather than evidence for a unique halo interpretation [2508.19382].

A more model-independent multimessenger analysis with agnostic hadronic and dark EOSs finds that current data at \(1\sigma\) constrain the DM fraction to
\[
f_{\mathrm{DM}} \lesssim 0.11
\]
for light dark matter. In that framework, GW170817 tidal deformability dominates the light, halo-dominated regime, whereas NICER mass-radius data dominate the heavy, core-dominated regime. The same study argues that neutron stars with similar masses but very different tidal deformabilities could be a smoking-gun signature of dark matter in neutron stars [2607.03840].

Candidate-object interpretations remain conditional. The heavy neutralino-like model for XTE J1814-338 reproduces the reported ultracompact \(M\)-\(R\) point only with a relatively large internal DM Fermi momentum, while HESS J1731-347 is modeled with a smaller admixture. Yet that same work states explicitly that small-radius compact stars are not definitive evidence for dark matter admixture because bosonic DM cores, hybrid stars with strong hadron-to-quark phase transitions, strange quark stars, and observational systematic effects remain plausible alternatives [2411.17105].

Several misconceptions can therefore be excluded. Dark matter is not a generic cure for the hyperon puzzle; in the additive fermion model it usually worsens the softening problem [1208.3722]. Dark matter does not generically compactify neutron stars; light halos can instead raise \(\Lambda\) [2607.03840]. And conclusions from one model class do not automatically transfer to another: a two-fluid halo/core system, a one-fluid additive EOS, and a chemically equilibrated neutron-dark-decay star probe distinct physical regimes [2506.18892].

The principal open issue is model dependence. Many studies parameterize the dark content phenomenologically rather than deriving it from capture, formation, or chemical kinetics; many assume gravity-only BM-DM coupling; many ignore rotation, temperature, magnetic fields, or nongravitational interactions; and observational bounds remain entangled with the uncertainty of the ordinary high-density EOS. A plausible implication is that future progress will depend less on any single “DMANS signature” than on combining multimessenger observables—mass, radius, tidal deformability, moment of inertia, and oscillation spectra—within frameworks that track dark morphology and visible-sector EOS uncertainty simultaneously.

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