⚛️ nuclear theory

NICER neutron stars with dark energy and dark matter: effects on the inferred equation of state

This study utilizes Bayesian inference on NICER mass-radius data to demonstrate that while dark energy cores significantly alter and can be constrained by neutron star observations, asymmetric dark matter cores remain observationally indistinguishable from purely baryonic stars for mass fractions up to 5%, implying that mass-radius measurements alone are insufficient to detect the latter.

Original authors: Nathan Rutherford, Chanda Prescod-Weinstein, Anna Watts

Published 2026-09-15

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Nathan Rutherford, Chanda Prescod-Weinstein, Anna Watts

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep in the cosmos, where gravity is so intense that it crushes matter into a state unlike anything found on Earth, lie neutron stars. These are the dense remnants of massive stars that have exploded, packing more mass than our Sun into a sphere only about the size of a city. Because they are so dense, they act as natural laboratories for physics, allowing scientists to test how matter behaves under extreme pressure. However, the universe is filled with mysteries that we cannot see directly: dark matter and dark energy. While ordinary matter makes up only a tiny fraction of the cosmos, dark matter and dark energy dominate the rest, yet we still do not know what they are made of. A central question in modern physics is whether these invisible substances can hide inside neutron stars, altering their size and weight in ways that our telescopes might detect.

A team of researchers recently used data from NASA's Neutron Star Interior Composition Explorer, known as NICER, to investigate this possibility. NICER is an X-ray telescope that measures the mass and radius of spinning neutron stars with remarkable precision. The scientists compared three different scenarios for what might be happening inside these stars. The first scenario assumes the stars are made entirely of ordinary, heavy matter. The second considers that the stars might contain a core of dark matter mixed in with the ordinary matter. The third, more exotic scenario, imagines that the very center of the star is a pocket of dark energy, surrounded by a shell of ordinary matter. By running complex computer simulations and comparing the results against the real measurements from NICER, the team sought to determine which of these hidden ingredients, if any, could be hiding in the cores of these cosmic giants.

The study focused on three specific neutron stars that NICER has observed in detail: PSR J0740+6620, PSR J0030+0451, and PSR J0437−4715. The researchers built mathematical models for each of the three scenarios and asked a simple question: which models produce stars that match the observed masses and sizes? When they looked at the possibility of dark matter cores, they found something surprising. Whether the dark matter was made of bosons or fermions—two different types of subatomic particles—the presence of a dark matter core did not change the star's appearance in a way that current telescopes could distinguish from a normal star. As long as the dark matter made up less than five percent of the star's total mass, the star looked exactly the same as a star made purely of ordinary matter. This suggests that simply measuring the size and weight of a neutron star is not enough to prove that dark matter is hiding inside it; we would need other methods to find it.

The results were very different when the researchers tested the dark energy scenario. In this model, the core of the star is not made of particles at all, but of a fluid with negative pressure, a property associated with the force driving the expansion of the universe. When the scientists included this dark energy core in their models, the predicted sizes and weights of the stars changed significantly. The models allowed for stars that could be much heavier and have different internal structures than the standard models. Crucially, the actual measurements from NICER were able to rule out many of the wilder possibilities for this dark energy fluid. The data narrowed down the range of allowed properties for this mysterious substance, effectively telling us what kind of dark energy could exist inside a neutron star and what kind cannot.

The researchers concluded that while current observations cannot confirm the presence of dark matter inside neutron stars, they can place very strict limits on the nature of dark energy cores. The data suggests that if dark energy exists in the center of these stars, it must follow very specific rules to be consistent with what we see. This finding is significant because it turns neutron stars into powerful tools for studying the dark sector of the universe. Even though we cannot see dark matter or dark energy directly, the way they would change the behavior of a neutron star leaves a fingerprint. For now, the evidence points to the idea that if dark matter is inside these stars, it is hiding in plain sight, indistinguishable from the ordinary matter around it. But if dark energy cores exist, the universe is giving us a way to measure them, tightening the noose on the theories that try to explain the invisible forces shaping our cosmos.

Technical Summary: NICER Neutron Stars with Dark Energy and Dark Matter

Problem Statement
The composition of the Universe is dominated by dark matter (DM) and dark energy (DE), yet their fundamental nature remains unknown. While the $\Lambda$ CDM model successfully explains large-scale cosmological observations, it faces challenges such as the Hubble tension and fine-tuning problems. Neutron stars (NSs) serve as extreme laboratories for probing these dark sectors. Previous literature suggests that DM can accumulate within NS interiors, altering their masses, radii, and tidal deformabilities. Similarly, if DE is distributed throughout the Universe, it could reside within NS cores, potentially forming hybrid stars with a dark energy core and a baryonic shell. However, the extent to which current observational data—specifically from NASA's Neutron Star Interior Composition Explorer (NICER)—can distinguish between purely baryonic stars, stars admixed with asymmetric dark matter (ADM), and stars with modified Chaplygin dark fluid (MCDF) cores remains unclear. This work investigates how accounting for these dark sector components affects the inferred equation of state (EoS) and neutron star properties.

Methodology
The authors employ a Bayesian inference framework to compare four distinct neutron star configurations against observational data:

Purely Baryonic Stars: Modeled using a piecewise polytropic (PP) EoS in the high-density regime, connected to chiral effective field theory ( $\chi$ EFT) results up to $1.5n_0$ (nuclear saturation density).
ADM-Admixed Stars: Neutron stars containing a core of either bosonic or fermionic asymmetric dark matter. These are modeled using the Nelson et al. ADM framework, where DM and baryons co-occupy the same spatial volume and interact primarily gravitationally. The structural properties are computed using the two-fluid Tolman-Oppenheimer-Volkoff (TOV) equations.
MCDF Hybrid Stars: Neutron stars featuring a dark energy core described by the Modified Chaplygin Dark Fluid (MCDF) EoS ( $P = A\epsilon - B/\epsilon$ ) surrounded by a baryonic shell. The transition between the non-mixing fluids is modeled using a sharp interface (Maxwell construction), and the structure is computed using standard single-fluid TOV equations.

The analysis utilizes mass-radius posterior distributions derived from NICER Pulse Profile Modeling (PPM) for three millisecond pulsars: PSR J0740+6620, PSR J0030+0451, and PSR J0437−4715. The inference is performed using the NEoST code (adapted for MCDF), employing nested sampling (MultiNest) to explore the parameter space of the EoS and dark sector parameters (e.g., ADM mass-fraction $F_\chi$ , MCDF parameters $A$ , $\alpha$ , and transition density $\epsilon^+_{dis}$ ).

Key Contributions and Results

Constraints on Dark Energy (MCDF):
- The presence of an MCDF core significantly broadens the posterior distributions for the total dense matter EoS, particularly at intermediate densities ( $10^{14.6} \text{--} 10^{15} \text{ g/cm}^3$ ).
- Unlike ADM models, MCDF cores predict substantially higher maximum masses ( $M_{TOV} \approx 2.29 M_\odot$ ) compared to the purely baryonic PP model ( $M_{TOV} \approx 2.16 M_\odot$ ).
- Crucially, the NICER mass-radius observations tightly constrain the MCDF parameter space. The posterior distributions for the barotropic pressure constant ( $A$ ) and the transition density ( $\epsilon^+_{dis}$ ) are significantly narrowed compared to their priors. Specifically, values of $A < 0.22$ and $\log_{10}(\epsilon^+_{dis}) > 15.0$ are strongly disfavored.
- However, the density jump parameter $\alpha$ remains unconstrained, as variations in $\alpha$ produce mass-radius changes smaller than current observational uncertainties.
Constraints on Dark Matter (ADM):
- Both bosonic and fermionic ADM cores broaden the total EoS prior at low densities ( $< 10^{14.65} \text{ g/cm}^3$ ) but converge to the purely baryonic PP model at higher densities.
- The inclusion of ADM cores reduces the inferred maximum mass slightly ( $M_{TOV} \approx 2.12 M_\odot$ ) but results in mass-radius posteriors that strongly coincide with the purely baryonic PP model.
- For ADM mass-fractions $F_\chi \leq 5\%$ , the posterior distributions for bosonic and fermionic ADM models are indistinguishable from each other and from the purely baryonic model.
- The analysis finds that current mass-radius measurements cannot constrain the ADM parameters or distinguish ADM cores from purely baryonic stars within this mass-fraction range.
Equation of State Inference:
- When isolating the baryonic EoS from the total model, the presence of an MCDF core biases the inferred baryonic EoS toward softer configurations at high densities to satisfy the pressure continuity at the core-shell interface.
- In contrast, ADM admixed models do not alter the posterior of the isolated baryonic EoS.
- All models predict a positive difference in radii between $2.0 M_\odot$ and $1.4 M_\odot$ stars ( $\Delta R = R_{2.0} - R_{1.4}$ ), with the MCDF model showing the most pronounced shift toward stiffening, though uncertainties remain too large for conclusive statements on high-density stiffness.

Significance and Claims
The paper demonstrates that while current NICER mass-radius measurements are insufficient to detect ADM cores inside neutron stars (even with improved uncertainties), they are highly effective at constraining the parameter space of dark energy cores. The study concludes that:

Dark Energy: Neutron stars can place strong constraints on the MCDF parameter space, ruling out broad regions of the $A$ and $\epsilon^+_{dis}$ parameter space. MCDF cores remain a viable explanation for the heaviest observed compact objects (e.g., PSR J0952−0607, GW190814) but require specific parameter values consistent with NICER data.
Dark Matter: Mass-radius measurements alone are insufficient to detect ADM cores for mass-fractions up to 5%. Independent astrophysical or particle physics probes are required to detect or constrain ADM accumulation in neutron stars.
Methodological: The work highlights the necessity of Bayesian inference to disentangle the interactions between the dark sector and neutron star observables, showing that different dark sector assumptions lead to distinct degeneracies in the inferred EoS.

The authors note that future work should address rotating configurations and the potential effects of dark matter halos, which were excluded from the current analysis due to computational complexity and the specific focus on core configurations.

More like this