Dark Cores and Halos: How Dark Matter Changes Neutron Stars
Dark matter admixed neutron stars are studied via a two-fluid formalism using agnostic hadronic and dark matter equations of state. Dark matter is modeled as a Fermi gas defined only by its low-density equation of state and mass. Hadronic matter is anchored at low densities by chiral effective field theory and at high densities by perturbative QCD; a sound-speed parametrization covers intermediate densities for hadrons and high densities for dark matter, ensuring thermodynamic consistency without bias. Light dark matter forms extended halos that increase tidal deformability, while heavy dark matter forms compact cores that decrease it. Observational constraints thus shift from gravitational-wave tidal deformability for light dark matter to NICER mass–radius data for heavy dark matter. At 1σ, the dark matter fraction is bounded to f_DM ≲ 0.11 for light dark matter. This agnostic framework yields conservative, broadly applicable bounds, and neutron stars with similar masses but very different tidal deformabilities could be a smoking-gun signature of dark matter.
Context
Neutron stars—the remnants of supernovae—are among the densest objects in the universe. Their central density exceeds nuclear density many times over, and their gravity is second only to black holes. That's why neutron stars serve as ideal laboratories for testing fundamental physics: from the behavior of matter at supranuclear densities to the properties of dark matter. Dark matter, predicted by Vera Rubin and born after the Big Bang, makes up a quarter of the universe's energy, yet remains undetected directly. Neutron stars can capture dark matter particles, altering their observable parameters—mass, radius, and tidal deformability, measured through gravitational waves. Modern data from pulsars, discovered by Jocelyn Bell Burnell and studied via X-ray spectroscopy (the NICER mission), and gravitational-wave detectors (LIGO, Rainer Weiss and others) let us constrain the dark matter fraction in these objects, without relying on specific models beyond the Standard Model.
Methods
To avoid bias, the scientists took an agnostic approach to the equations of state for both nuclear and dark matter. The nuclear part is anchored by chiral effective field theory calculations at low densities and perturbative quantum chromodynamics at high densities, with the intermediate region parameterized by random sound speed profiles (c²_s), limited only by causality (the sound speed cannot exceed the speed of light) and thermodynamic stability. For dark matter, a similar parameterization was used but without a high-density anchor: at densities below 0.1 of nuclear saturation, it was described as a degenerate Fermi gas with particle mass m_D; above that, by sound speed interpolation. The star was treated as a two-fluid system: the two components interact only through gravity, and each satisfies its own hydrostatic equilibrium equation (a modified Tolman–Oppenheimer–Volkoff equation). Ensembles of ~10⁵ mass–radius sequences were generated for various m_D (0.2–1.1 GeV) and dark matter fractions f_DM (0.01–0.15), then compared with observational data.
Results
The analysis revealed that dark matter's influence depends dramatically on its particle mass. Light dark matter (m_D ~ 0.2–0.3 GeV) forms a diffuse halo around the star, increasing its effective radius. This boosts the tidal deformability Λ (proportional to R⁵), and such configurations are tightly constrained by data from the gravitational-wave event GW170817. In contrast, heavy dark matter (m_D > 0.7 GeV) concentrates in a compact core, shrinking the star and lowering Λ. Here, the main constraints come from NICER's mass and radius measurements. At intermediate particle masses (0.5–0.7 GeV), both scenarios are possible, along with transitional "core-halo" structures. A statistical analysis using all available data (NICER, GW170817, and a maximum mass requirement of ≥ 2.01 M⊙) found that for light dark matter, the allowed fraction is no more than 11%. For heavy dark matter, constraints are stronger at higher fractions, but overall the parameter space remains much broader. An intriguing prediction is that spotting two neutron stars with similar masses but wildly different tidal deformabilities could be a "smoking gun"—a hint of dark matter in one of them.
Implications
These results are model-independent and apply to a broad class of dark matter theories. They show that modern multimessenger observations can significantly constrain the dark matter content in neutron stars even without knowing its microscopic nature. The specific constraint—tidal deformability or mass-radius—serves as a diagnostic tool: the dominance of gravitational-wave constraints points to halo structures, while X-ray constraints dominating suggests core dark matter configurations.
Future development
In the future, new NICER measurements for other pulsars, plus data from next-generation gravitational-wave detectors (including space-based ones), should significantly sharpen these limits. Further development of the agnostic method may incorporate bosonic dark matter and account for stellar rotation, bringing models closer to reality. Studying the dependence on the chosen anchor density for the dark sector is also promising.
Impact
The results will impact dark matter physics, neutron star astrophysics, and gravitational-wave astronomy, offering a new tool for searching for light dark matter particle candidates.
Next steps
Next steps include analyzing higher-precision data from future missions like STROBE-X or Athena, and redoing the study for bosonic dark matter.
Key open problems
The research is directly tied to the unsolved problem of dark matter's nature and the question of neutron star composition at extreme densities, where exotic matter might emerge.
🎯 A teaspoon of neutron star material would weigh about a billion tons—roughly as much as Mount Everest.
Key numbers
- maximum neutron star mass: 2.01 M⊙
- dark matter particle mass range: 0.2–1.1 GeV
- allowed dark matter fraction (light): ≲0.11
- tidal deformability constraint: 70 ≤ Λ₁.₄ ≤ 580
- central density: up to several n₀ (~0.16 fm⁻³)