pith. machine review for the scientific record. sign in

arxiv: 2604.07618 · v1 · submitted 2026-04-08 · ✦ hep-ph

Recognition: no theorem link

Multi Component Dark Matter in a Minimal Model

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:52 UTC · model grok-4.3

classification ✦ hep-ph
keywords multi-component dark matterZ2 symmetrysinglet fermionssinglet scalarHiggs portalrelic densitydirect detectionneutrino floor
0
0 comments X

The pith

In a minimal Z2-symmetric model with two singlet fermions and one scalar, a viable parameter region matches the observed dark matter relic density while keeping direct detection signals below experimental limits.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs a simple extension of the Standard Model by adding two singlet fermions and one singlet scalar, all protected by a Z2 symmetry that prevents their decay into lighter particles. These three states interact with ordinary matter only through the Higgs boson, creating a multi-component dark matter scenario. Parameter space is scanned to find regions where the sum of their thermal relic densities equals the measured dark matter abundance. The scalar component couples at tree level to nucleons through the Higgs portal, yet can be tuned to supply only a small fraction of the total density so that its scattering rate stays below current bounds. The two fermions, whose scattering proceeds only through loops, naturally sit below the neutrino floor yet can carry most of the relic density.

Core claim

The model introduces two singlet fermions and a singlet scalar under a Z2 symmetry that communicate with SM particles through a scalar-Higgs portal. Regions in the parameter space are found where all three particles are kinematically stable, the total relic density matches the observed value, the scalar DM evades direct detection bounds with minimal relic contribution, and the loop-suppressed fermion DM cross sections lie below the neutrino floor while possessing a large fraction of the total relic density.

What carries the argument

The Z2 symmetry that kinematically stabilizes the two singlet fermions and the singlet scalar, together with the scalar-Higgs portal that sets both the thermal freeze-out rates and the tree-level versus loop-level direct detection cross sections.

If this is right

  • The scalar dark matter can occupy only a small share of the total relic density and still remain consistent with existing direct detection limits because its tree-level rate scales with that small abundance.
  • The fermion dark matter components can supply the dominant fraction of the observed density while their loop-induced scattering cross sections automatically fall below the neutrino floor.
  • Standard thermal freeze-out calculations suffice to determine the individual relic fractions once the masses and portal couplings are fixed.
  • The total relic density can be tuned to the observed value across a continuous range of mass hierarchies among the three particles.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Experiments that push sensitivity below the neutrino floor would directly test the fermion-dominated parts of this parameter space.
  • The same Z2-protected particles could be produced at colliders through Higgs decays, offering an independent search channel not explored in the paper.
  • Adding a small non-thermal production mechanism for the scalar would further enlarge the viable region without altering the loop suppression of the fermions.

Load-bearing premise

The three new particles remain stable because their masses forbid all decays allowed by the Z2 symmetry, and their abundances are set by ordinary thermal freeze-out without extra cosmological effects.

What would settle it

A direct detection experiment reporting a dark matter-nucleon scattering rate above the neutrino floor in the mass window allowed by the model would rule out the fermion-dominated relic region.

Figures

Figures reproduced from arXiv: 2604.07618 by Karim Ghorbani.

Figure 1
Figure 1. Figure 1 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: All the necessary interactions to compute the relic densit [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: shown is the viable region in the mass parameter space, res [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The ratio ξφ/ξ1, the relative size of the relic density of the scalar DM and the fermion DM, is shown in terms of the scalar mass. The corresponding fermion mass is also shown in a vertical indicator. Bound from invisible Higgs decay is not imposed. the scalar mass and the mass of a fermion. The results in [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Shown are the Feynman diagrams for the fermionic DM scat [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The elastic DM-nucleon cross section for three DM compon [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
read the original abstract

We study a minimal model with an imposed $\mathbb{Z}_2$ symmetry incorporating two singlet fermions and a singlet scalar which communicate with the SM particles through a scalar-Higgs portal. We probe regions in the parameter space where stability of the three new particles are guaranteed kinematically, and thus introducing a multi component dark matter (DM) scenario. The regions in the parameter space with predicted total relic density in accordance with the observed DM abundance are found, and the contribution of each species to the total relic density is determined. The elastic DM-nucleon scattering cross section of the two fermions DM are loop suppressed, while that of the scalar DM starts at tree level and thus presumably being dominant. It is found that there exist a viable region in the parameter space that the scalar DM can evade the current direct detection (DD) experimental bounds while it has a minimal contribution to the observed relic density. The DD cross section of the fermion DM being loop suppressed, resides below the lower limit from the {\it neutrino floor}, possessing a large fraction of the total relic density.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 1 minor

Summary. The manuscript introduces a minimal Z2-symmetric model containing two singlet fermions and one singlet scalar as multi-component dark matter, with interactions to the SM mediated solely by a scalar-Higgs portal. Parameter regions are identified in which all three new states are kinematically stable, the sum of their thermal relic densities matches the observed value, the scalar contributes only a small fraction (allowing its tree-level DD cross section to evade current bounds), and the fermions supply the dominant fraction with loop-suppressed DD cross sections lying below the neutrino floor.

Significance. If the relic-density partitioning and DD calculations are robust, the work would illustrate how a minimal multi-component setup can simultaneously satisfy the observed abundance and current direct-detection limits by assigning the tree-level portal coupling to the sub-dominant component. This is a standard but still useful demonstration for DM model-building. No machine-checked proofs, public code, or parameter-free derivations are provided, so the result remains dependent on the numerical implementation.

major comments (3)
  1. [Relic density calculation] The relic-density analysis determines the individual contributions of the scalar and two fermions by applying the standard single-species Boltzmann equation to each component separately and then summing. Because the scalar couples to the Higgs and, to communicate with the fermions, must possess Yukawa-type interactions with them, number-changing processes such as ψψ ↔ SS and ψS ↔ ψS are present. These channels couple the three Boltzmann equations; omitting them means the reported minimal scalar fraction and dominant fermion fraction are not guaranteed to hold, directly affecting the claimed viable region that evades DD bounds.
  2. [Direct detection cross sections] The fermion DD cross sections are asserted to be loop-suppressed and to lie below the neutrino floor while still accounting for a large fraction of the total relic density. No explicit one-loop diagrams, suppression factors, or numerical values evaluated at the selected benchmark points are supplied, preventing verification that the loop suppression is sufficient across the reported viable parameter space.
  3. [Parameter space scan] Viable points are obtained by scanning parameters so that the total relic density equals the observed value and then computing the individual contributions and DD cross sections from those same points. This procedure renders the evasion of bounds a post-selection outcome rather than an independent prediction, weakening the claim that the model naturally accommodates both relic density and DD constraints.
minor comments (1)
  1. [Abstract] The abstract contains the vague qualifier 'presumably being dominant' for the scalar DD cross section; this should be replaced by a quantitative statement once the cross sections are computed.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and valuable feedback on our manuscript. We address each of the major comments point by point below. Where appropriate, we have revised the manuscript to incorporate the referee's suggestions and strengthen the presentation of our results.

read point-by-point responses
  1. Referee: [Relic density calculation] The relic-density analysis determines the individual contributions of the scalar and two fermions by applying the standard single-species Boltzmann equation to each component separately and then summing. Because the scalar couples to the Higgs and, to communicate with the fermions, must possess Yukawa-type interactions with them, number-changing processes such as ψψ ↔ SS and ψS ↔ ψS are present. These channels couple the three Boltzmann equations; omitting them means the reported minimal scalar fraction and dominant fermion fraction are not guaranteed to hold, directly affecting the claimed viable region that evades DD bounds.

    Authors: We agree that the Boltzmann equations are in principle coupled through the conversion processes enabled by the Yukawa interactions. However, in the parameter regions we consider, the masses of the fermions and scalar are such that the conversion processes are kinematically suppressed or occur at rates much smaller than the Hubble expansion at the time of freeze-out. We have added a new subsection in the revised manuscript discussing the validity of the decoupled approximation, including estimates of the relevant reaction rates, to support our results. A fully coupled numerical solution is beyond the scope of this work but would be an interesting extension. revision: partial

  2. Referee: [Direct detection cross sections] The fermion DD cross sections are asserted to be loop-suppressed and to lie below the neutrino floor while still accounting for a large fraction of the total relic density. No explicit one-loop diagrams, suppression factors, or numerical values evaluated at the selected benchmark points are supplied, preventing verification that the loop suppression is sufficient across the reported viable parameter space.

    Authors: We apologize for the omission. In the revised manuscript, we have included the relevant one-loop Feynman diagrams for the fermion-nucleon scattering, along with the analytical expressions for the suppression factors. Additionally, we provide numerical values of the DD cross sections for the benchmark points listed in the paper, confirming that they lie below the neutrino floor. revision: yes

  3. Referee: [Parameter space scan] Viable points are obtained by scanning parameters so that the total relic density equals the observed value and then computing the individual contributions and DD cross sections from those same points. This procedure renders the evasion of bounds a post-selection outcome rather than an independent prediction, weakening the claim that the model naturally accommodates both relic density and DD constraints.

    Authors: The scan is used to delineate the viable parameter space where the total relic density matches the observed value. Within this space, we identify sub-regions where the scalar's contribution is minimal, allowing it to evade DD bounds, while the fermions dominate the density with loop-suppressed cross sections. This demonstrates that the model can accommodate the constraints through the different interaction structures of the components. We have clarified the language in the manuscript to emphasize that we are showing the existence of such viable regions rather than claiming an independent prediction without selection. revision: no

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper identifies regions of parameter space where the total relic density matches the observed value via standard thermal freeze-out for each Z2-odd species, then computes direct-detection cross sections (tree-level for the scalar, loop-suppressed for the fermions) from those same points. This is a conventional scan procedure in which relic density supplies a constraint and DD observables are derived quantities for the surviving points; the two are independent observables and the existence claim does not reduce to a tautology or self-definition. No equations or self-citations in the abstract or described text exhibit any of the enumerated circular patterns, and the derivation remains self-contained against external relic-density and DD benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The central claim rests on introducing new particles stabilized by an ad-hoc symmetry and fitting multiple free parameters to match the observed relic density while deriving DD rates; since only the abstract is available, exact parameter values and full derivation details are unknown.

free parameters (2)
  • Masses of the two singlet fermions and the singlet scalar
    Chosen to ensure kinematic stability of all three particles and to reproduce the observed total relic density.
  • Higgs portal coupling strength
    Determines the interaction rates that control both relic density contributions and direct detection cross sections.
axioms (2)
  • domain assumption Imposition of Z2 symmetry to stabilize the new particles
    The Z2 is imposed by hand on the new fields to prevent their decay into SM particles or each other.
  • standard math Standard thermal freeze-out in Lambda-CDM cosmology for relic density
    Assumes the usual Boltzmann equation solution in standard Big Bang cosmology without non-standard effects.
invented entities (1)
  • Two singlet fermions and one singlet scalar no independent evidence
    purpose: To serve as stable multi-component dark matter candidates
    New particles are postulated in the model; no independent evidence is provided beyond the assumption of kinematic stability.

pith-pipeline@v0.9.0 · 5471 in / 1708 out tokens · 67991 ms · 2026-05-10T16:52:39.454339+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Topological solitons of two-field scalar theories in rotationally symmetric backgrounds

    hep-th 2026-05 unverdicted novelty 6.0

    Explicit radial dependence in the scalar potential stabilizes topological solitons in higher-dimensional rotationally symmetric backgrounds, enabling exact solutions via shared target-space orbits and a geometry-encod...

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Bergstr¨ om,Nonbaryonic dark matter: Observational evidence and detection methods,Rept

    L. Bergstr¨ om, “Nonbaryonic dark matter: Observational ev idence and detection methods,” Rept. Prog. Phys. 63 (2000) 793, arXiv:hep-ph/0002126. 9

  2. [2]

    Feng,Dark Matter Candidates from Particle Physics and Methods of Detection,Ann

    J. L. Feng, “Dark Matter Candidates from Particle Physics and M ethods of Detection,” Ann. Rev. Astron. Astrophys. 48 (2010) 495–545, arXiv:1003.0904 [astro-ph.CO]

  3. [3]

    The Waning of the WIMP? A Review of Models, Searches, and Constraints

    G. Arcadi, M. Dutra, P. Ghosh, M. Lindner, Y. Mambrini, M. Pierre , S. Profumo, and F. S. Queiroz, “The waning of the WIMP? A review of models, searche s, and constraints,” Eur. Phys. J. C 78 no. 3, (2018) 203, arXiv:1703.07364 [hep-ph]

  4. [4]

    Precise Relic WIMP Abundance and its Impact on Searches for Dark Matter Annihilation

    G. Steigman, B. Dasgupta, and J. F. Beacom, “Precise Relic WIMP Abundance and its Impact on Searches for Dark Matter Annihilation,” Phys. Rev. D 86 (2012) 023506, arXiv:1204.3622 [hep-ph]

  5. [5]

    McDonald, Phys

    J. McDonald, “Gauge singlet scalars as cold dark matter,” Phys. Rev. D 50 (1994) 3637–3649, arXiv:hep-ph/0702143

  6. [6]

    Singlet fermionic dark matter

    Y. G. Kim, K. Y. Lee, and S. Shin, “Singlet fermionic dark matter,” JHEP 05 (2008) 100, arXiv:0803.2932 [hep-ph]

  7. [7]

    Fermionic dark matter with pseudo-scalar Yukawa interaction,

    K. Ghorbani, “Fermionic dark matter with pseudo-scalar Yukawa interaction,” JCAP 1501 (2015) 015, arXiv:1408.4929 [hep-ph]

  8. [8]

    Pseudoscalar P ortal Dark Matter and New Signatures of Vector-like Fermions,

    J. Fan, S. M. Koushiappas, and G. Landsberg, “Pseudoscalar P ortal Dark Matter and New Signatures of Vector-like Fermions,” JHEP 01 (2016) 111, arXiv:1507.06993 [hep-ph]

  9. [9]

    Fermionic Dark Matter through a Light Pseudoscala r Portal: Hints from the DAMA Results,

    K.-C. Yang, “Fermionic Dark Matter through a Light Pseudoscala r Portal: Hints from the DAMA Results,” Phys. Rev. D 94 no. 3, (2016) 035028, arXiv:1604.04979 [hep-ph]

  10. [10]

    Minimal renormalizable simplified dark ma tter model with a pseudoscalar mediator,

    S. Baek, P. Ko, and J. Li, “Minimal renormalizable simplified dark ma tter model with a pseudoscalar mediator,” Phys. Rev. D 95 no. 7, (2017) 075011, arXiv:1701.04131 [hep-ph]

  11. [11]

    Electroweak Baryogenesis and Dark Matter v ia a Pseudoscalar vs. Scalar,

    P. H. Ghorbani, “Electroweak Baryogenesis and Dark Matter v ia a Pseudoscalar vs. Scalar,” JHEP 08 (2017) 058, arXiv:1703.06506 [hep-ph]

  12. [12]

    Leading Loop Effects in Pseud oscalar-Higgs Portal Dark Matter,

    K. Ghorbani and P. H. Ghorbani, “Leading Loop Effects in Pseud oscalar-Higgs Portal Dark Matter,” JHEP 05 (2019) 096, arXiv:1812.04092 [hep-ph]

  13. [13]

    Díaz Sáez, P

    B. D ´ ıaz S´ aez, P. Escalona, S. Norero, and A. R. Zerwekh, “Fermion singlet dark matter in a pseudoscalar dark matter portal,” JHEP 10 (2021) 233, arXiv:2105.04255 [hep-ph]

  14. [14]

    Light thermal dark matter beyond p-wave annihilation in minimal Higgs portal model,

    Y.-T. Chen, S. Matsumoto, T.-P. Tang, Y.-L. S. Tsai, and L. Wu, “Light thermal dark matter beyond p-wave annihilation in minimal Higgs portal model,” JHEP 05 (2024) 281, arXiv:2403.02721 [hep-ph] . 10

  15. [15]

    Bhattacharya, A

    S. Bhattacharya, A. Drozd, B. Grzadkowski, and J. Wudka, “ Two-Component Dark Matter,” JHEP 10 (2013) 158, arXiv:1309.2986 [hep-ph]

  16. [16]

    Search for invisible decays of a Higgs boson produced through vector boson fusion in proton-proton c ollisions at √ s = 13 TeV,

    CMS Collaboration, A. M. Sirunyan et al. , “Search for invisible decays of a Higgs boson produced through vector boson fusion in proton-proton c ollisions at √ s = 13 TeV,” Phys. Lett. B 793 (2019) 520–551, arXiv:1809.05937 [hep-ex]

  17. [17]

    Alguero, G

    G. Alguero, G. Belanger, F. Boudjema, S. Chakraborti, A. Gou delis, S. Kraml, A. Mjallal, and A. Pukhov, “micrOMEGAs 6.0: N-component dark matt er,” Comput. Phys. Commun. 299 (2024) 109133, arXiv:2312.14894 [hep-ph]

  18. [18]

    Planck 2018 results. VI. Cosmological parameters

    Planck Collaboration, N. Aghanim et al. , “Planck 2018 results. VI. Cosmological parameters,” Astron. Astrophys. 641 (2020) A6, arXiv:1807.06209 [astro-ph.CO] . [Erratum: Astron.Astrophys. 652, C4 (2021)]

  19. [19]

    Alguero, G

    G. Alguero, G. Belanger, S. Kraml, and A. Pukhov, “Co-scatte ring in micrOMEGAs: A case study for the singlet-triplet dark matter model,” SciPost Phys. 13 (2022) 124, arXiv:2207.10536 [hep-ph]

  20. [20]

    Dark Matter Search Results from a One Ton-Year Exposure of XENON1T,

    XENON Collaboration, E. Aprile et al. , “Dark Matter Search Results from a One Ton-Year Exposure of XENON1T,” Phys. Rev. Lett. 121 no. 11, (2018) 111302, arXiv:1805.12562 [astro-ph.CO]

  21. [21]

    First Dark Matter Search with Nuclear Recoils from the XENONnT Experiment,

    XENON Collaboration, E. Aprile et al. , “First Dark Matter Search with Nuclear Recoils from the XENONnT Experiment,” Phys. Rev. Lett. 131 no. 4, (2023) 041003, arXiv:2303.14729 [hep-ex]

  22. [22]

    Billardet al., Rept

    J. Billard et al. , “Direct Detection of Dark Matter – APPEC Committee Report,” arXiv:2104.07634 [hep-ex] . 11