SciPost logo

SciPost Submission Page

Strongly Interacting Dark Matter admixed Neutron Stars

by Yannick Dengler, Suchita Kulkarni, Axel Maas, Kevin Radl

Submission summary

Authors (as registered SciPost users): Axel Maas
Submission information
Preprint Link: https://arxiv.org/abs/2503.19691v2  (pdf)
Date submitted: Aug. 22, 2025, 11:14 a.m.
Submitted by: Axel Maas
Submitted to: SciPost Physics Core
Ontological classification
Academic field: Physics
Specialties:
  • High-Energy Physics - Phenomenology
Approach: Phenomenological

Abstract

Dark matter may accumulate in neutron stars given its gravitational interaction and abundance. We investigate the modification of neutron star properties and confront them with the observations in the context of strongly-interacting dark matter scenario, specifically for a QCD-like theory with G$_2$ gauge group for which a first-principles equation-of-state from lattice calculations is available. We study the impact of various observational constraints and modeling of the QCD equation of state on the combined neutron stars. The results indicate that dark matter masses of a few hundred MeV to a few GeV are consistent with the latest observed neutron star properties.

Current status:
Awaiting resubmission

Reports on this Submission

Report #2 by Anonymous (Referee 2) on 2025-10-2 (Invited Report)

Strengths

There is a novel part of this work which is related to the construction of the EoS for the dark matter component of the star directly from lattice calculations. Based on that the authors can produce mass-radius relations and deformabilities relevant also for gravitational wave production.

Weaknesses

The model is at best contrived. The authors choose a very specific gauge group for dark matter namely G2 in order to my understanding to avoid having to deal with the notorious sign problem that arises in nonzero density studies of QCD-like theories. However there is not much theoretical justification and motivation for such a peculiar group. It is highly unlike that dark matter will have such a gauge interaction.

In addition contrary to what the authors claim, that this model produces unique mass-radius relations, this is hard to believe. One could easily mimic similar mass-radius relations by having dark matter particles exhibit some repulsive or attractive self-interactions via mediation of scalars or vectors with appropriate couplings. It is well known in neutron star physics that lots of EoS end up producing the same mass-radius relations and therefore it is hard to pinpoint which particular EoS is correct based on sole observations of mass and radius. The authors have to explain better in what sense they consider unique these EoS or they should simply tune down their statement.

Another relevant issue within the context of the studied model is how dark matter abundance is produced in the first place. The dark sector will have to go through a stage of confinement as it happened with QCD. The authors should comment on how their model produces the observed amount of dark matter abundance. Related to that is also the issue of dark matter self-interactions. If this component composes 100% of dark matter, they are constrained by limits from the bullet cluster and the ellipticity of galaxies. If not they should state more clearly that this is a small component of dark matter that must be less than 5% in order to avoid the aforementioned constraints.

Finally although the authors talk about accretion, there is practically not so much time to accrete such a huge component of mass in the lifetime of a neutron star. Therefore accretion is not the main mechanism of this admixed star formation. The authors should clarify this point and present or comment on a reasonable formation scenario.

I also see a couple of flaws not related to physics but rather to the way the paper is written. Firstly, there is a vast literature on dark matter effects on neutron stars ranging from cooling effects, black hole formation, admixed/dark stars etc that the authors ignore. They should cite relevant work on dark matter effects on neutron stars. Furthermore the abstract has to be rewritten so it reflects in a more precise way what has been done in this study. Accumulation as I mentioned does not seem to be valid. Moreover it should be clarified in the abstract if they refer to a 100% component of dark matter or a subdominant one, clarifying also that they choose a particular gauge group to model the dark matter interactions. At first read I thought that G2 is implemented to deal with the QCD baryonic part of the star.

Report

The authors have to answer/address all the issues mentioned.

Recommendation

Ask for major revision

  • validity: low
  • significance: low
  • originality: ok
  • clarity: ok
  • formatting: reasonable
  • grammar: excellent

Report #1 by Anonymous (Referee 1) on 2025-9-19 (Invited Report)

Disclosure of Generative AI use

The referee discloses that the following generative AI tools have been used in the preparation of this report:

I used the generative AI tool only to polish the language of my review. The scientific evaluation, analysis, and conclusions are fully my own.

Strengths

1-The work is timely and conceptually valuable.
2-Introducing a lattice-derived EOS for dark matter into two-fluid neutron star modelling is a novelty and useful blueprint for future studies.
3-The main qualitative conclusions, light $m_c$ tends to halos while heavy $m_c$ tends to cores; sub-percent DM fractions yield tidal properties close to ordinary neutron stars, are plausible and well illustrated.
4-The authors perform an extended investigation by varying the central pressures of both ordinary and dark matter components and systematically scanning the parameter space according to their chosen criterion [Eq. (4.1)].
5-By considering a broad range of observables, the study provides a useful framework that future works can employ for comparison and further refinement.
6-The study employs three different ordinary matter EOSs with varying stiffness, allowing for an exploration of model dependence in the ordinary matter sector.
7-For the dark matter sector, two distinct EOSs are considered, which broadens the scope of the analysis and strengthens the robustness of the results.
8-The authors carry out a careful and detailed analysis of the ordinary and dark components separately, specifically the masses (M_O,M_D) and radii (R_O,R_D), which is particularly valuable in the context of two-fluid dark-matter–admixed neutron stars.

Weaknesses

1-Observational constraints (mass, radius, tidal deformability) not applied tightly enough or not shown under stricter cuts. 2-Tidal deformability analysis is too limited (only the softest EOS and the DM mass corresponding to the core configuration have been presented). 3-Figures are not easily readable, and the corresponding texts are not clarifying (the connection to M–R phenomenology is unclear.) 4-While the broad range of observables is valuable, it does not establish the consistency of the selected parameters with current astrophysical measurements. As a result, the dark matter parameters are not effectively constrained by the presented calculations. Consequently, the main conclusion stated in the abstract, that dark matter masses of a few hundred MeV to a few GeV are consistent with the latest neutron star observations, appears overstated, since the latest observational constraints were not actually incorporated into the analysis. 5-Lack of standard fixed central-pressure ratio or fixed dark matter fraction sequences for a full M-R line, which are widely used in neutron star literature. 6-Unclear treatment of dark matter self-interaction (implicit in lattice EOS but not clarified). 7- Lack of any report about the speed of sound. 8-In one instance, the authors cite around 30 references immediately after the sentence “The amount of dark matter may vary depending on the age, history, and location of the neutron star, but most estimates conclude that it is very unlikely to exceed 1% of the total neutron-star mass.” This is unusual and should be reconsidered.

Report

The authors need to clarify some questions and remarks before publication.
With this clarification, the manuscript will offer a clear and robust reference for dark matter admixed neutron stars using first-principles dark-sector input. The additions do not change the scope but improve clarity and completeness.

Requested changes

1- First of all, I would like to emphasize that the authors should compare their results directly with the most recent observational constraints on neutron stars if they intend to propose their parameter space for dark matter (even if it is only the dark matter particle mass) as a consistent range. In particular, they need to present both mass–radius and tidal deformability relations in a way that clearly illustrates the behavior of their models relative to the data. They may clarify this with some plots with exact observational constraints included, or by providing tables or texts which explain explicitly which range of parameters and which EOSs are consistent with the observational data (For example, $\Lambda_{1.4}$<580). This is the main missing part of the paper. 2-The authors should elaborate on the role of dark matter self-interactions in their framework. 3-The authors should at least provide some comments on the speed of sound in the employed EoSs. 4-The authors state that they interpolate piecewise polytropes between NChPT at low densities and pQCD at high densities. However, the maximum energy density reached in their models is well below the regime where pQCD is valid. This statement should therefore be revised to accurately reflect the interpolation procedure actually used. 5-On page 7, the authors claim that the crust would hardly affect the radius. This is not correct, as the crust has a non-negligible impact on neutron star radii. The statement about not including a crust should be modified accordingly. 6-For clarity and consistency, I recommend revising the use of terminology and abbreviations. Some common terms (e.g., "equation of state," "neutron stars") are repeated many times without abbreviation, while certain abbreviations (e.g., WIMPs) appear without definition, or are defined only after their first occurrence (e.g., NChPT). Defining each abbreviation upon its first use and employing it consistently thereafter would improve readability and flow. 7-Correct the thermodynamic identity to $n\mu = \varepsilon + p$ (instead of $n\mu = p\varepsilon$) before Eq. (6.1). 8-Correct the expression for the chemical potential to $\mu = c + \tfrac{\Gamma}{\Gamma - 1} K n^{\Gamma - 1}$, not $\Gamma/(\Gamma+1)$. 9-Review the manuscript for typographical errors such as “we use choose …” or “is a a system …” and correct them throughout.

Recommendation

Ask for minor revision

  • validity: good
  • significance: good
  • originality: good
  • clarity: low
  • formatting: reasonable
  • grammar: good

Login to report or comment