Pith. Sign in

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.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 27, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.34234 (version 1)
Pith IDC2BTGHGS, Pith's permanent record of this paper and its review
CitationDavood Rafiei Karkevandi, Gordon Baym, Hamid Reza Moshfegh, Mahboubeh Shahrbaf, Reza Karimi. "Probing In-Medium Neutron-to-Dark-Fermion Conversion with Multimessenger Neutron Star Observations." arXiv:2610.34234 (2026). https://pith.science/paper/C2BTGHGS
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

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

The paper studies neutron-to-dark-fermion conversion in neutron stars, where a single dark fermion χ with baryon number 1 is in chemical equilibrium with neutrons. The authors scan dark-ferion mass mχ from 939 to 1500 MeV and repulsive self-interaction strength Gχ from ǂ to 100 fm² across 21 relativistic mean-field (RMF) nucleonic equations of state. They confront the resulting stellar sequences with three multimessenger constraints: support for at least 2M⊙ neutron stars, GW170817 tidal deformability bounds, and NICER mass-radius inferences (both Illinois–Maryland and X‑PSI/Amsterdam analyses separately). The robust finding is that all 21 EOSs exclude the same region of parameter space: mχ up to about 1300 MeV and weak dark repulsion (Gχ ≲ 11 fm²). This means that any viable dark sector with a single fermion near the neutron mass must have sufficiently strong self-interactions to avoid excessive softening of the neutron star equation of state.

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-27: 1 request, 8,137 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.74 · human participants 0.01 · controlled comparison 0.42 · randomized assignment 0.01 · states quantitative result 0.97 · compares with prior results 0.58 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

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

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
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 kind of paper the classifier reports this is: simulation study 0.80 · evidence from analysis of existing data 0.58.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1. Purely nucleonic reference sequences for the 21
Figure 2. FIG. 2. Effect of neutron-to-dark-fermion conversion on the
Figure 3. FIG. 3. Particle fractions for the representative EOS exam
Figure 4. FIG. 4. Mass–radius relations (upper row) and dimensionless tidal deformability as a function of mass (lower row) for EOS 8
Figure 5. FIG. 5. Exclusion boundaries in the (
Figure 6. FIG. 6. Exclusion boundaries in the (
Figure 7. FIG. 7. EOS-by-EOS strict-consistency regions for the 13 main baseline EOSs under the strict X-PSI/Amsterdam criterion.
Figure 8. FIG. 8. EOS-overlap strict-consistency map for the 13 main
Figure 9. FIG. 9. Illustration of the mass–radius classification for the baseline EOS 20 using the Illinois–Maryland and X-PSI/Amsterdam
Figure 11. FIG. 11. Conservative exclusion structure for the six stiff
Figure 12. FIG. 12. Stiff-EOS not-excluded boundaries after adding
Figure 13. FIG. 13. Fully EOS-overlapped excluded regions for the

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.