⚛️ nuclear theory

A self-consistent Higgs-portal framework for dark matter--admixed neutron stars: Collider-motivated benchmarks meet multimessenger constraints

This paper proposes a self-consistent single-fluid relativistic mean-field framework for dark matter-admixed neutron stars using a Higgs-portal model with a massive $Z^\prime$ mediator, demonstrating how the resulting repulsive interaction dynamically couples baryonic and dark matter sectors to systematically soften the equation of state, reduce maximum stellar mass, and establish a direct link between collider-motivated WIMP models and multimessenger astrophysical constraints.

Original authors: Adamu Issifu

Published 2026-07-17

📖 3 min read🧠 Deep dive

CC BY 4.0

Original authors: Adamu Issifu

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

Imagine the universe as a giant, cosmic ocean. We can see the waves crashing on the shore—the stars, the planets, and us—but we know there is something else swimming beneath the surface. We can't see it, we can't touch it, but we know it's there because it pulls on the waves we can see. This invisible stuff is called "dark matter." It makes up about 27% of everything in the universe, while the stuff we are made of (like atoms and people) is a tiny, tiny fraction. Scientists have been trying to figure out what this dark matter is made of for decades, but it's been a ghostly puzzle. One of the most exciting places to look for clues isn't in a giant telescope or a particle collider, but inside the most extreme objects in the cosmos: neutron stars. These are the collapsed cores of dead stars, so dense that a teaspoon of their material would weigh a billion tons. If dark matter exists, it might get trapped inside these stars, hiding in plain sight and changing how they behave.

This paper is like a detective story where the authors build a new, super-smart simulation to see what happens if dark matter gets cozy inside a neutron star. They focus on a specific type of dark matter called a "WIMP" (Weakly Interacting Massive Particle), which is a heavy, slow-moving particle that might interact with normal matter through a "Higgs portal"—a theoretical doorway that connects our world to the dark one. The big problem with previous studies was that they treated the dark matter inside the star as a static guest, like a piece of furniture that was just placed there and didn't move. This paper argues that's wrong. Instead, they treat the dark matter as a living, breathing part of the star that reacts to the pressure and density of the star's core, just like the normal matter does.

The authors created a new mathematical framework to simulate this interaction. They found that when you let the dark matter adjust itself naturally inside the star, it acts like a heavy, invisible weight that squeezes the star from the inside out. This makes the star's internal "stiffness" (how hard it is to compress) drop significantly. In their simulations, adding even a tiny amount of these dark particles (about 0.2% of the total number of particles) caused the maximum possible mass of a neutron star to shrink by up to 29%. It's as if the star became a bit more like a squishy marshmallow and less like a hard rock.

They also discovered that the dark matter doesn't just sit there; it interacts with the normal matter through a "repulsive force" (mediated by a heavy particle they call a Z' boson), which pushes back against the star's gravity. This interaction is so strong that it changes the way sound waves travel through the star, making them slower. The study suggests that if we observe neutron stars that are smaller or lighter than we expect, it might be a sign that they are hiding a secret stash of dark matter. However, the authors are careful to note that these results come from computer simulations based on their specific model; they haven't found the dark matter yet, but they have provided a new, more realistic way to look for it. If future telescopes and gravitational wave detectors see stars that match these "squishier" predictions, it could be the first real evidence of what dark matter actually is.

Technical Summary: A Self-Consistent Higgs-Portal Framework for Dark Matter–Admixed Neutron Stars

Problem Statement
The theoretical modeling of Dark Matter (DM)-admixed neutron stars (DANSs) within a single-fluid framework has historically relied on phenomenological prescriptions that lack astrophysical self-consistency. Specifically, conventional approaches treat the local DM Fermi momentum ( $k^\chi_F$ ) or the local DM density ( $n_\chi$ ) as externally fixed parameters, decoupling the DM distribution from the stellar baryonic environment. This ad hoc freezing of the DM density neglects the physical processes of DM capture, thermalization, and gravitational settling within the star. Furthermore, while the Higgs-portal model is a widely used particle-physics realization for DM interactions, its implementation in neutron star (NS) studies often fails to dynamically couple the DM and baryonic matter (BM) sectors, leading to an incomplete description of the equation of state (EOS) and thermodynamic response.

Methodology
To address these limitations, the authors develop a self-consistent single-fluid relativistic mean-field (RMF) framework that extends the standard Higgs-portal model. The methodology involves the following key steps:

Model Extension: The authors introduce a massive neutral vector mediator, $X_\mu$ (identified as a dark $Z'$ boson), associated with a hidden $U(1)_X$ gauge symmetry. This mediator generates a density-dependent repulsive vector mean-field that dynamically couples the BM and DM sectors by modifying their chemical potentials. This approach is analogous to the role of the $\omega$ -meson in standard RMF theory but applied to the DM sector.
Self-Consistent Coupling: Instead of prescribing a fixed $k^\chi_F$ , the framework imposes a global particle fraction, $F_\chi = N_\chi/N_B$ , as the sole free astrophysical parameter. The local DM density is derived self-consistently as $n_\chi = F_\chi n_B$ , ensuring that the DM distribution evolves naturally with the baryonic density profile.
Particle Physics Benchmarks: The model parameters are constrained by collider-motivated benchmarks (ATLAS) and direct detection limits. The DM mass is set to $m_\chi = 200$ GeV (WIMP paradigm), the mediator mass to $m_X = 1.5$ TeV, and couplings are fixed to $g_{qX} = 0.25$ and $g_{\chi X} = 1$ . The Higgs-DM coupling is set to $g_h = 0.07$ .
Equation of State (EOS) Derivation: The authors derive exact analytical expressions for the total energy density and pressure, incorporating contributions from the scalar Higgs portal and the vector $Z'$ interaction. Crucially, they derive rigorous analytical expressions for the adiabatic speed of sound ( $c_s^2$ ), its density derivative, and the adiabatic index ( $\Gamma$ ). These derivations explicitly account for the density dependence of $n_\chi$ and $k^\chi_F$ , distinguishing the framework from models with fixed $k^\chi_F$ .
Numerical Implementation: The framework is applied using three distinct RMF EOSs for baryonic matter: NL3 $\omega\rho$ , FSU2R, and DD2. The Tolman–Oppenheimer–Volkoff (TOV) equations are solved to determine the macroscopic structure (mass, radius, binding energy) of the stars for varying values of $F_\chi$ (0 to 0.2%).

Key Contributions

Self-Consistent Formalism: The paper establishes a framework where the local DM density is not an external input but a self-consistent outcome of the global particle fraction and local baryonic density ( $n_\chi = F_\chi n_B$ ). This eliminates the need for arbitrary prescriptions of the DM Fermi momentum.
Analytical Thermodynamic Benchmarks: The authors provide exact analytical expressions for $c_s^2$ , $dc_s^2/dn_B$ , and $\Gamma$ in the context of coupled BM-DM systems. These expressions identify the specific microscopic contributions (kinetic, scalar, vector, and interaction terms) to the thermodynamic response, offering a rigorous benchmark for assessing causality and stability.
Interaction Decomposition: The study quantitatively decomposes the BM-DM interaction energy into scalar (Higgs) and vector ( $Z'$ ) components. It demonstrates that for the adopted collider-motivated benchmarks, the repulsive vector channel overwhelmingly dominates the scalar channel (by a factor of $>20$ ), providing a microscopic explanation for the observed macroscopic effects.
Micro-Macro Connection: The work bridges particle physics constraints (collider limits on $Z'$ mass and couplings) with astrophysical observables (mass-radius relations, tidal deformability), creating a unified framework for multimessenger constraints.

Results

EOS Softening: Increasing the DM fraction $F_\chi$ systematically softens the nuclear EOS. This occurs because the heavy DM particles ( $m_\chi \gg m_N$ ) contribute significantly to the energy density but provide negligible pressure support. Consequently, the maximum NS mass ( $M_{\text{max}}$ ) decreases by approximately 25–29% as $F_\chi$ increases from 0 to 0.2%, with the reduction being most pronounced for softer EOSs (FSU2R).
Stellar Compactness: The presence of DM leads to more compact stellar configurations. For a canonical 1.4 $M_\odot$ star, the radius decreases by 2–3 km, and the central baryon density increases by up to ~20%.
Thermodynamic Response: The adiabatic speed of sound ( $c_s^2$ ) is systematically suppressed, particularly at high densities, due to the dynamical contribution of the DM sector to thermodynamic derivatives. The density derivative of $c_s^2$ also shows a monotonic decrease with increasing $F_\chi$ , indicating a more gradual thermodynamic response.
Binding Energy Anomalies: The study reports that for high $F_\chi$ , the conventional gravitational binding energy ( $E_b = M_B - M$ ) can become negative. The authors clarify that this does not imply hydrostatic instability (which is determined by $dM/d\varepsilon_c > 0$ ) but rather reflects the decoupling of the baryon mass budget ( $M_B$ ) from the total rest mass ( $M_0$ ), as the latter includes the substantial rest mass of the heavy WIMPs.
Microscopic Properties: The DM remains deeply non-relativistic throughout the stellar interior ( $k^\chi_F \lesssim 100$ MeV, $\mu_\chi \approx m_\chi$ ). The chemical potential shift is driven primarily by the density-dependent vector interaction, while the scalar contribution remains subdominant.

Significance
The paper claims that enforcing self-consistency in the treatment of the dark sector is essential for obtaining physically meaningful single-fluid EOSs. By deriving exact analytical expressions for thermodynamic quantities, the framework provides a general benchmark for assessing the stability and causality of DM-admixed matter. The work establishes a direct, quantitative link between collider-motivated WIMP models and multimessenger astrophysics. It demonstrates that compact stars can serve as complementary laboratories to terrestrial experiments, where observational constraints from gravitational waves (e.g., GW170817) and X-ray timing (e.g., NICER) can be translated into limits on the amount of DM accumulated inside neutron stars for a fixed particle-physics scenario. The results suggest that massive pulsars provide a stringent upper bound on the DM content, and that the specific signatures of heavy WIMP admixture (such as EOS softening and modified $c_s^2$ profiles) can be distinguished from uncertainties in the nuclear EOS.

Technical Summary: A Self-Consistent Higgs-Portal Framework for Dark Matter–Admixed Neutron Stars

More like this