REVIEW
Probing In-Medium Neutron-to-Dark-Fermion Conversion with Multimessenger Neutron Star Observations
Nuclear Theory · Submitted 25 Oct 2026 · 136 references
Pith's verdict: The most systematic multimessenger constraint on the single-χ conversion model to date, well executed and transparent.
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
Neutron stars rule out weak dark repulsion for dark fermions near neutron mass.
The reading
What carries the argument
The chemical‑equilibrium condition µχ=µn ensures that neutron‑to‑dark‑fermion conversion is treated thermodynamically. The dark fermion carries baryon number 1 and its repulsive self-action is modelled via an effective vector coupling Gχ. The nucleonic backgrounds are an ensemble of 21 unified RMF equations of state that vary in high‑density stiffness and symmetry‑energy slope. Stellar sequences are computed and tested against a mass‑radius constraint, a tidal‑deformability constraint, and the NICER mass‑radius credible regions. The EOS‑overlapped common excluded region identifies parameter points excluded by every one of the 21 EOSs.
Extended reading notes
Core claim
Across a 21-EOS ensemble constructed within a non-linear σ‑ω‑ρ relativistic mean-field framework, and using the combined requirements of a 2M⊙ maximum mass, GW170817 tidal deformability (70≤Λ≤580), and the NICER mass-radius inferences, the authors find that the low‑Gχ part of the dark-sector parameter space is robustly excluded for all EOSs. For the Illinois–Maryland analysis the common excluded region covers mχ=939–1275 MeV and Gχ=1–10 fm²; for the X‑PSI/Amsterdam analysis it covers mχ=939–1300 MeV and Gχ=1–11 fm². Weak dark self-repulsion is thus disfavoured across a broad dark‑fermion mass range, while the remaining not‑excluded parameter space is highly dependent on the nucleonic EOS and the choice of NICER analysis.
Load-bearing premise
The 21 nucleonic RMF EOSs, drawn from a single non‑linear σ‑ω‑ρ family, adequately span the range of possible dense‑matter behavior, and the dark sector contains only a single χ fermion with baryon number 1 and repulsive vector self‑interactions, with no other interactions with baryons.
Editorial extensions
If this is right
- Dark‑sector models that contain a single fermion with baryon number 1 and weak repulsive self‑interactions (Gχ≲10 fm²) near the neutron mass are robustly excluded by current neutron star observations, irrespective of which of the 21 RMF EOSs is adopted.
- The common not‑excluded region is much smaller and more EOS‑dependent, implying that future improvements in the neutron‑star equation of state (e.g., through gravitational‑wave detections or better X‑ray timing) will likely narrow the allowed dark‑sector parameter space further.
- The method demonstrates that neutron stars can probe GeV‑scale dark sectors through in‑medium conversion channels that are inaccessible to terrestrial experiments, even when the dark fermion is heavier than the neutron.
- Repulsive dark self‑interactions strong enough to compensate for the EOS softing caused by conversion are required to keep the maximum mass above 2M⊙, making such interactions a necessary ingredient for models that try to explain the neutron‑life anomaly via dark decay.
- The exclusions are robust across two independent NICER mass‑radius inference sets (Illinois–Maryland and X‑PSI/Amsterdam), strengthening confidence in the result despite current differences between the two analyses.
Reading between the lines
- If the dark sector contains interactions beyond the single repulsive vector, such as scalar attraction or coupling between dark and baryonic sectors, the exclusions could be modified; exploring those scenarios would test the robustness of the present result.
- The absence of a common not‑excluded region for the full 21OS ensemble implies that any claim of a 'favored' dark‑matter parameter space from neutron star data must be conditioned on the chosen nucleonic EOS; the present analysis provides a more honest representation of the current uncertainty.
- The method could be extended to multi‑particle dark decay channels (e.g., n→χχχ) where the dark particles carry fractional baryon number, possibly avoiding the need for strong self‑interactions and testing a different region of parameter space.
- Combining the exlusion with future observations of even more massive or more compact neutron stars will likely push the bounds to larger mχ and higher Gχ, potentially ruling out all but the most repulsive dark‑sector scenarios.
Original abstract
Neutron stars provide a natural environment for probing in-medium neutron-to-dark-fermion conversion and its impact on the dense-matter equation of state (EOS). We use multimessenger neutron star observations to map excluded and not-excluded regions of dark-sector parameter space for an EOS ensemble within a non-linear $\sigma$--$\omega$--$\rho$ relativistic mean-field (RMF) framework and determine their overlap across the ensemble. We include a fermionic dark component with repulsive self-interactions of strength $G_\chi$. We consider masses $m_\chi$ near the neutron mass and heavier states accessible in neutron star cores at sufficiently large neutron chemical potentials. To reduce sensitivity to the nucleonic EOS within this framework, we scan $m_\chi=939$--$1500$ MeV and $G_\chi=1$--$100$ fm$^2$ across 21 nucleonic EOSs: a main baseline set of 13 EOSs and soft and stiff limiting cases. Stellar sequences are tested against the requirement that the EOS support at least $2M_\odot$ neutron stars, GW170817 tidal-deformability bounds, and NICER mass--radius constraints, considering the Illinois--Maryland and X-PSI/Amsterdam inference sets separately. For the full ensemble of 21 EOSs, both NICER analyses identify similar regions excluded by all EOSs, with approximate minimum--maximum extents of $m_\chi=939$--$1300$ MeV and $G_\chi=1$--$11$ fm$^2$. Unlike the full 21-EOS ensemble, the main baseline set has a common not-excluded region whose extent depends on the NICER analysis. The soft and stiff limiting cases show how these regions shift when the nucleonic EOS predicts a maximum neutron star mass near $2M_\odot$ or large tidal deformabilities. Across these 21 calibrated nucleonic RMF EOSs, weak dark-sector repulsion is robustly disfavored across a broad dark-fermion mass range, while the extent of the remaining not-excluded region depends on the NICER analysis and baryonic EOS stiffness.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
Neutron stars provide a natural environment for probing in-medium neutron-to-dark-fermion conversion and its impact on the dense-matter equation of state (EOS). We use multimessenger neutron star observations to map excluded and not-excluded regions of dark-sector parameter space for an EOS ensemble within a non-linear $\sigma$--$\omega$--$\rho$ relativistic mean-field (RMF) framework and determine their overlap across the ensemble. We include a fermionic dark component with repulsive self-interactions of strength $G_\chi$. We consider masses $m_\chi$ near the neutron mass and heavier states accessible in neutron star cores at sufficiently large neutron chemical potentials. To reduce sensitivity to the nucleonic EOS within this framework, we scan $m_\chi=939$--$1500$ MeV and $G_\chi=1$--$100$ fm$^2$ across 21 nucleonic EOSs: a main baseline set of 13 EOSs and soft and stiff limiting cases. Stellar sequences are tested against the requirement that the EOS support at least $2M_\odot$ neutron stars, GW170817 tidal-deformability bounds, and NICER mass--radius constraints, considering the Illinois--Maryland and X-PSI/Amsterdam inference sets separately. For the full ensemble of 21 EOSs, both NICER analyses identify similar regions excluded by all EOSs, with approximate minimum--maximum extents of $m_\chi=939$--$1300$ MeV and $G_\chi=1$--$11$ fm$^2$. Unlike the full 21-EOS ensemble, the main baseline set has a common not-excluded region whose extent depends on the NICER analysis. The soft and stiff limiting cases show how these regions shift when the nucleonic EOS predicts a maximum neutron star mass near $2M_\odot$ or large tidal deformabilities. Across these 21 calibrated nucleonic RMF EOSs, weak dark-sector repulsion is robustly disfavored across a broad dark-fermion mass range, while the extent of the remaining not-excluded region depends on the NICER analysis and baryonic EOS stiffness.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | Neutron stars provide a natural environment for probing in-medium neutron-to-dark-fermion conversion and its impact on the dense-matter equation of state (EOS). | background 1.00 | 0.08 | 0.04 | asserted 0.84 |
| s1 | We use multimessenger neutron star observations to map excluded and not-excluded regions of dark-sector parameter space for an EOS ensemble within a non-linear $\sigma$--$\omega$--$\rho$ relativistic mean-field (RMF) framework and determine their overlap across the ensemble. | contribution 0.52 | 0.96 | 0.07 | asserted 0.99 |
| s2 | We include a fermionic dark component with repulsive self-interactions of strength $G_\chi$. | method 0.92 | 0.81 | 0.02 | asserted 0.91 |
| s3 | We consider masses $m_\chi$ near the neutron mass and heavier states accessible in neutron star cores at sufficiently large neutron chemical potentials. | method 0.86 | 0.63 | 0.02 | asserted 0.80 |
| s4 | To reduce sensitivity to the nucleonic EOS within this framework, we scan $m_\chi=939$--$1500$ MeV and $G_\chi=1$--$100$ fm$^2$ across 21 nucleonic EOSs: a main baseline set of 13 EOSs and soft and stiff limiting cases. | method 0.99 | 0.95 | 0.03 | asserted 0.97 |
| s5 | Stellar sequences are tested against the requirement that the EOS support at least $2M_\odot$ neutron stars, GW170817 tidal-deformability bounds, and NICER mass--radius constraints, considering the Illinois--Maryland and X-PSI/Amsterdam inference sets separately. | method 0.99 | 0.89 | 0.04 | asserted 0.93 |
| s6 | For the full ensemble of 21 EOSs, both NICER analyses identify similar regions excluded by all EOSs, with approximate minimum--maximum extents of $m_\chi=939$--$1300$ MeV and $G_\chi=1$--$11$ fm$^2$. | result 1.00 | 0.94 | 0.22 | asserted 1.00 |
| s7 | Unlike the full 21-EOS ensemble, the main baseline set has a common not-excluded region whose extent depends on the NICER analysis. | result 1.00 | 0.93 | 0.09 | asserted 1.00 |
| s8 | The soft and stiff limiting cases show how these regions shift when the nucleonic EOS predicts a maximum neutron star mass near $2M_\odot$ or large tidal deformabilities. | result 0.60 | 0.88 | 0.05 | asserted 0.99 |
| s9 | Across these 21 calibrated nucleonic RMF EOSs, weak dark-sector repulsion is robustly disfavored across a broad dark-fermion mass range, while the extent of the remaining not-excluded region depends on the NICER analysis and baryonic EOS stiffness. | result 0.75 | 0.93 | 0.87 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Showing 12 of the paper's 13 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2610.34234 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.