Pith. sign in

REVIEW 2 major objections 1 minor 46 references

Two real scalar fields can form an asymmetric two-component FIMP dark matter model that reproduces the observed relic density.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

2026-06-29 21:33 UTC pith:GM6IMZIG

load-bearing objection This paper adds one concrete two-component asymmetric scalar FIMP example with quartic asymmetry transfer, but the relic density match is obtained by selecting benchmark parameters rather than predicting it. the 2 major comments →

arxiv 2605.25961 v1 pith:GM6IMZIG submitted 2026-05-25 hep-ph

Asymmetric Two-Component Scalar FIMP Dark Matter

S. Peyman Zakeri This is my paper
classification hep-ph
keywords dark matterFIMPasymmetric dark matterfreeze-in mechanismtwo-component dark matterscalar fieldsHiggs portalrelic density
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper constructs a model extending the Standard Model with two real scalar dark matter fields and a heavy scalar mediator, protected by a Z2 times Z2 prime symmetry. Dark matter is produced through the freeze-in mechanism via the Higgs portal, while the decay of the mediator creates an asymmetry in the first component that is partially shared with the second through quartic couplings. Numerical solution of the Boltzmann equations for benchmark parameters of 0.1 GeV and 0.5 GeV masses, a portal coupling of 7 times 10 to the minus 11, and eta equal to 0.01 yields the measured relic density, with the second component accounting for about 5 percent. The model keeps self-interactions far below observational limits and predicts an invisible Higgs decay rate well below current bounds.

Core claim

We propose a two-component asymmetric FIMP dark matter model in which both DM candidates are real scalar fields. The model is an extension of the Standard Model by two scalar DM components and a heavy scalar mediator, stabilized by a Z2 times Z2 prime symmetry. DM is produced via the freeze-in mechanism through the Higgs portal, while the out-of-equilibrium decay of the heavy mediator generates an asymmetry in the first component, which is partially transferred to the second component via quartic interactions. By solving the Boltzmann equations numerically, we compute the relic density and perform a detailed scan over the parameter space. The observed relic density is successfully reproduced

What carries the argument

The out-of-equilibrium decay of the heavy mediator generating asymmetry in the first DM component with partial transfer to the second via quartic interactions.

Load-bearing premise

The quartic interactions between the two dark matter components transfer the asymmetry generated by the mediator decay in a manner that permits the numerical solution of the Boltzmann equations to match the observed relic density.

What would settle it

Numerical integration of the coupled Boltzmann equations for the benchmark parameters m_phi1=0.1 GeV, m_phi2=0.5 GeV, lambda1H=7e-11, eta=0.01 that yields a total Omega_DM h^2 significantly different from 0.12 would falsify the reproduction claim.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The total relic density matches observations with the second component contributing only 5 percent.
  • The DM self-interaction cross section lies orders of magnitude below the Bullet cluster bound and the double radio relic limit.
  • The invisible Higgs decay branching ratio is about 4 times 10 to the minus 19, well below the LHC upper limit.
  • Direct detection prospects are negligible due to the small Higgs portal couplings.

Where Pith is reading between the lines

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

  • The small mass and coupling values suggest this mechanism could apply to other light scalar dark matter scenarios without conflicting with current data.
  • Adjusting the eta parameter might allow the second component to contribute varying fractions while still satisfying the relic density.
  • The stability provided by the separate Z2 symmetries could be tested by searching for signatures of the heavy mediator in collider data.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 1 minor

Summary. The paper proposes a two-component asymmetric FIMP dark matter model consisting of two real scalar fields stabilized by a Z₂ × Z₂' symmetry, along with a heavy scalar mediator. DM is produced via freeze-in through the Higgs portal, while the out-of-equilibrium decay of the mediator generates an asymmetry in the first component that is partially transferred to the second via quartic interactions. Numerical solution of the Boltzmann equations reproduces the observed relic density Ω_DM h² = 0.12 ± 0.001 for the benchmark parameters m_φ1=0.1 GeV, m_φ2=0.5 GeV, λ1H=7×10^{-11}, and η=0.01 (with the second component contributing ~5%), and the model satisfies self-interaction, invisible Higgs decay, and direct detection constraints.

Significance. If the numerical reproduction is robust, the work establishes a novel connection between asymmetric DM, freeze-in production, and multi-component FIMP scenarios. The explicit numerical Boltzmann solution and constraint checks (e.g., σ/m well below Bullet cluster and double radio relic limits, Br(h→inv) ~4×10^{-19}) are strengths that could motivate further model-building in this direction.

major comments (2)
  1. [Numerical results section] Numerical results section (and abstract): the observed relic density is reproduced for the stated benchmark values, but these parameters (m_φ1, m_φ2, λ1H, η) are selected to achieve the fit; the manuscript should report the scan ranges, priors, and fraction of points that reproduce Ω_DM h² within the quoted uncertainty to demonstrate that the match is not by construction.
  2. [Boltzmann equations section] Boltzmann equations section: the central mechanism relies on asymmetry generation via mediator decay and transfer via quartic interactions, but the explicit coupled equations for the number densities and chemical potentials of both components (including the transfer term proportional to η) are not provided, preventing verification of the ~5% contribution from the second component.
minor comments (1)
  1. [Abstract] Abstract: the invisible Higgs decay branching ratio is quoted without reference to the relevant equation or section where it is derived; a cross-reference would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and recommendation for minor revision. We address each major comment below and will incorporate the requested clarifications.

read point-by-point responses

  1. Referee: [Numerical results section] Numerical results section (and abstract): the observed relic density is reproduced for the stated benchmark values, but these parameters (m_φ1, m_φ2, λ1H, η) are selected to achieve the fit; the manuscript should report the scan ranges, priors, and fraction of points that reproduce Ω_DM h² within the quoted uncertainty to demonstrate that the match is not by construction.

    Authors: We agree that reporting the scan details would strengthen the presentation. In the revised manuscript we will add the parameter ranges scanned, the priors employed, and the fraction of points that reproduce Ω_DM h² = 0.12 ± 0.001, thereby showing that the benchmark is representative of viable parameter space rather than finely tuned. revision: yes

  2. Referee: [Boltzmann equations section] Boltzmann equations section: the central mechanism relies on asymmetry generation via mediator decay and transfer via quartic interactions, but the explicit coupled equations for the number densities and chemical potentials of both components (including the transfer term proportional to η) are not provided, preventing verification of the ~5% contribution from the second component.

    Authors: We acknowledge the value of explicit equations for verification. The revised manuscript will include the full set of coupled Boltzmann equations for the number densities and chemical potentials of both components, explicitly showing the mediator-decay asymmetry source and the η-proportional transfer term. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper defines a two-component scalar FIMP model stabilized by Z2 x Z2', writes Boltzmann equations incorporating freeze-in production, out-of-equilibrium mediator decay, and quartic-mediated asymmetry transfer, then numerically integrates those equations over a parameter scan. Matching the external Planck relic-density value for chosen benchmark points (m_φ1, m_φ2, λ1H, η) is a standard viability demonstration, not a reduction of the output to the input by construction. No self-definitional relations, no fitted quantities relabeled as predictions, and no load-bearing self-citations appear in the derivation chain. The numerical solution remains independently falsifiable against cosmological data.

Axiom & Free-Parameter Ledger

4 free parameters · 2 axioms · 2 invented entities

The model rests on several new fields and couplings introduced without independent evidence beyond fitting the relic density, plus standard assumptions of the freeze-in framework.

free parameters (4)
  • m_φ1 = 0.1 GeV
    Light DM mass chosen as benchmark to achieve correct relic density
  • m_φ2 = 0.5 GeV
    Second DM mass chosen as benchmark to achieve correct relic density
  • λ1H = 7e-11
    Higgs portal coupling chosen as benchmark to achieve correct relic density
  • η = 0.01
    Quartic coupling parameter chosen as benchmark to achieve correct relic density and asymmetry transfer
axioms (2)
  • domain assumption Dark matter is produced via the freeze-in mechanism through the Higgs portal
    Standard assumption invoked for FIMP production in the model setup
  • domain assumption The heavy mediator decays out of equilibrium to generate asymmetry
    Core mechanism assumption for asymmetry generation
invented entities (2)
  • heavy scalar mediator no independent evidence
    purpose: Generates asymmetry in first DM component via out-of-equilibrium decay
    New field postulated to produce the required asymmetry
  • Z2 x Z2' symmetry no independent evidence
    purpose: Stabilizes the two scalar DM components
    New discrete symmetry introduced to prevent decay into SM particles

pith-pipeline@v0.9.1-grok · 5846 in / 1752 out tokens · 34692 ms · 2026-06-29T21:33:23.776682+00:00 · methodology

0 comments
read the original abstract

We propose a two-component asymmetric FIMP (feebly interacting massive particle) dark matter (DM) model in which both DM candidates are real scalar fields. The model is an extension of the Standard Model (SM) by two scalar DM components and a heavy scalar mediator, stabilized by a $\mathbb{Z}_2 \times \mathbb{Z}_2'$ symmetry. DM is produced via the freeze-in mechanism through the Higgs portal, while the out-of-equilibrium decay of the heavy mediator generates an asymmetry in the first component, which is partially transferred to the second component via quartic interactions. By solving the Boltzmann equations numerically, we compute the relic density and perform a detailed scan over the parameter space. The observed relic density $\Omega_{\text{DM}} h^2 = 0.12 \pm 0.001$ is successfully reproduced for benchmark parameters $m_{\phi_1}=0.1$~GeV, $m_{\phi_2}=0.5$~GeV, $\lambda_{1H}=7\times10^{-11}$, and $\eta=0.01$, with the second component contributing only about $5\%$ to the total abundance. We also examine phenomenological constraints. The DM self-interaction cross section lies orders of magnitude below the Bullet cluster bound ($\sigma/m < 0.47$~cm$^2$/g) and the more stringent double radio relic limit ($\sigma/m < 0.22$~cm$^2$/g). The invisible Higgs decay branching ratio is $\sim 4\times10^{-19}$, well below the LHC upper limit, and direct detection prospects are negligible due to the small Higgs portal couplings. Our model establishes a novel connection between two-component DM, the freeze-in mechanism, and DM asymmetry.

Figures

Figures reproduced from arXiv: 2605.25961 by S. Peyman Zakeri.

Figure 1
Figure 1. Figure 1: FIG. 1: Total DM abundance [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Total DM abundance [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Relic density Ω [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Relic density Ω [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Self-interaction parameter space for both components ( [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

46 extracted references · 5 canonical work pages · 3 internal anchors

  1. [1]

    Dependence onλ 1H In Fig. 1, we plot the total DM abundanceY DM =Y total ϕ1 +Y total ϕ2 as a function of the temperatureTfor four different values of the Higgs portal couplingλ 1H = 5×10 −11,7× 10−11,1×10 −10,1.5×10 −10, while keepingη= 0.01 fixed. The temperature range is restricted to 2–11 GeV to better visualize the differences between the curves. As e...

  2. [2]

    2, we show the temperature evolution ofY DM for four different values of the asymmetry parameterη= 0,0.01,0.05,0.1, while fixingλ 1H = 7×10 −11

    Dependence onη In Fig. 2, we show the temperature evolution ofY DM for four different values of the asymmetry parameterη= 0,0.01,0.05,0.1, while fixingλ 1H = 7×10 −11. The effect of the asymmetry is clearly visible: largerηleads to a larger overall abundance. This is because the total abundance ofϕ 1 is given by (1+η)Y sym ϕ1 , while the symmetricY sym ϕ1...

  3. [3]

    3, we plot the relic density as a function ofλ 1H for four different values of the asymmetry parameterη= 0,0.01,0.05,0.1

    Dependence onλ 1H In Fig. 3, we plot the relic density as a function ofλ 1H for four different values of the asymmetry parameterη= 0,0.01,0.05,0.1. As expected from the freeze-in mechanism, the relic density scales as Ω DMh2 ∝λ 2 1H. The figure also clearly shows the effect of the asymmetry: for a fixedλ 1H, largerηleads to a larger relic density. This is...

  4. [4]

    4, we show the relic density as a function of the asymmetry parameterηfor four different values ofλ 1H = 5×10 −11,7×10 −11,1×10 −10,1.5×10 −10

    Dependence onη In Fig. 4, we show the relic density as a function of the asymmetry parameterηfor four different values ofλ 1H = 5×10 −11,7×10 −11,1×10 −10,1.5×10 −10. The results exhibit a clear linear dependence, Ω DMh2 ∝(1 +η), which follows directly fromY total ϕ1 = (1 +η)Y sym ϕ1 . The slope of each curve is determined by the symmetric partY sym ϕ1 , ...

  5. [5]

    Invisible Higgs Decay The Higgs boson can decay invisibly into pairs of DM particles through the portal cou- plingsλ 1H andλ 2H. The partial decay width forh→ϕ iϕi is given by Γh→ϕiϕi = λ2 iHv2 32πmh s 1− 4m2 ϕi m2 h .(16) For our benchmark parameters, withλ 1H = 7×10 −11,λ 2H = 7×10 −12, andm ϕi ≪m h, the decay widths are Γ h→ϕ1ϕ1 ≈2.36×10 −20 GeV and Γ ...

  6. [6]

    In our model, the spin-independent scattering cross section is proportional toλ 2 iH and is further suppressed by the small portal couplings

    Direct Detection Direct detection experiments search for elastic scattering of DM particles off nuclei. In our model, the spin-independent scattering cross section is proportional toλ 2 iH and is further suppressed by the small portal couplings. For the freeze-in regime withλ iH ∼10 −11–10−12, the predicted cross section lies far below the current sensiti...

  7. [7]

    V. C. Rubin and W. K. Ford, Jr., Astrophys. J.159, 379 (1970)

    1970

  8. [8]

    Clowe et al., Astrophys

    D. Clowe et al., Astrophys. J.648, L109 (2006)

    2006

  9. [9]

    Astrophys.641, A6 (2020)

    Planck Collaboration, Astron. Astrophys.641, A6 (2020)

    2020

  10. [10]

    Ya. B. Zel’dovich, Zh. Eksp. Teor. Fiz.48, 986 (1965); Ya. B. Zel’dovich, L. B. Okun, and S. B. Pikelner, Usp. Fiz. Nauk84, 113 (1965)

    1965

  11. [11]

    H. -Y. Chiu, Phys. Rev. Lett.17, 712 (1966)

    1966

  12. [12]

    Aalbers et al

    J. Aalbers et al. (LZ Collaboration), Phys. Rev. Lett.131, 041002 (2023)

    2023

  13. [13]

    Aprile et al

    E. Aprile et al. (XENON Collaboration), Phys. Rev. Lett.131, 041003 (2023)

    2023

  14. [14]

    L. J. Hall, K. Jedamzik, J. March-Russell, and S. M. West, JHEP1003, 080 (2010)

    2010

  15. [15]

    Bernal, M

    N. Bernal, M. Heikinheimo, T. Tenkanen, K. Tuominen, and V. Vaskonen, Int. J. Mod. Phys. A32, 1730023 (2017)

    2017

  16. [16]

    S. Y. Ayazi, S. M. Firouzabadi, and S. P. Zakeri, J. Phys. G43, 095006 (2016)

    2016

  17. [17]

    B´ elanger, S

    G. B´ elanger, S. Chakraborti, and A. Pukhov, Eur. Phys. J. Spec. Top.233, 2135 (2024)

    2024

  18. [18]

    Profumo, K

    S. Profumo, K. Sigurdson, and L. Ubaldi, JCAP0912, 016 (2009)

    2009

  19. [19]

    Bhattacharya, A

    S. Bhattacharya, A. Drozd, B. Grzadkowski, and J. Wudka, JHEP1310, 158 (2013)

    2013

  20. [20]

    S. Esch, M. Klasen, and C. E. Yaguna, JHEP1409, 108 (2014)

    2014

  21. [21]

    J. P. Carvalho-Corrˆ ea, I. M. Pereira, B. L. S´ anchez-Vega, and A. C. D. Viglioni, Eur. Phys. J. C85, 1353 (2025)

    2025

  22. [22]

    S. P. Zakeri, S. M. Moosavi Nejad, M. Zakeri, and S. Y. Ayazi,Chinese Physics C42, 073101 (2018)

    2018

  23. [23]

    B´ elanger, A

    G. B´ elanger, A. Pukhov, C. E. Yaguna, and ´O. Zapata, JHEP09, 030 (2020)

    2020

  24. [24]

    Bhattacharya, P

    S. Bhattacharya, P. Ghosh, A. K. Saha, and A. Sil, JHEP2003, 090 (2020)

    2020

  25. [25]

    Chatterjee and A

    S. Chatterjee and A. Das, Phys. Rev. D111, 015023 (2025). 16

    2025

  26. [26]

    Kumar and G

    N. Kumar and G. Lekshmi, Phys. Rev. D110, 123015 (2024)

    2024

  27. [27]

    Tulin and H.-B

    S. Tulin and H.-B. Yu, Phys. Rept.730, 1 (2018)

    2018

  28. [28]

    Yang, Y.-L

    D. Yang, Y.-L. S. Tsai, and Y.-Z. Fan, arXiv:2504.02303 (2025)

    work page arXiv 2025

  29. [29]

    Bhattacharya, D

    S. Bhattacharya, D. Mahanta, N. Mondal, and D. Pradhan, JCAP09, 032 (2025)

    2025

  30. [30]

    Nussinov, Phys

    S. Nussinov, Phys. Lett. B165, 55 (1985)

    1985

  31. [31]

    Petraki and R

    K. Petraki and R. R. Volkas, Int. J. Mod. Phys. A28, 1330028 (2013)

    2013

  32. [32]

    Becker and W.-C

    M. Becker and W.-C. Huang, arXiv:1911.06788 (2019)

    work page arXiv 1911

  33. [33]

    Kapustin, Phys

    A. Kapustin, Phys. Rev. D59, 035006 (1999)

    1999

  34. [34]

    Dutta and J

    B. Dutta and J. Kumar, Phys. Lett. B699, 364 (2011)

    2011

  35. [35]

    A. C. Ritter and R. R. Volkas, Phys. Rev. D110, 015032 (2024)

    2024

  36. [36]

    Bodas, M

    A. Bodas, M. A. Buen-Abad, A. Hook, and R. Sundrum, JHEP06, 070 (2024)

    2024

  37. [37]

    Falkowski, J

    A. Falkowski, J. T. Ruderman, and T. Volansky, JHEP1105, 106 (2011)

    2011

  38. [38]

    Neutron Portal and Dark Matter-Baryon Coincidence: from UV Completion to Phenomenology

    S. Girmohanta, R. S. Gupta, Y. Shigekami, and A. Thapa, arXiv:2604.21168 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  39. [39]

    Unwin, JHEP10, 190 (2014)

    J. Unwin, JHEP10, 190 (2014)

    2014

  40. [40]

    Goudelis, D

    A. Goudelis, D. Karamitros, P. Papachristou, and V. C. Spanos, Eur. Phys. J. C82, 921 (2022)

    2022

  41. [41]

    K. Asai, S. Kanemura, and K. Matsuda, Phys. Lett. B836, 137627 (2023)

    2023

  42. [42]

    S. P. Zakeri, arXiv:2602.20570 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  43. [43]

    Harvey, R

    D. Harvey, R. Massey, T. Kitching, A. Taylor, and E. Tittley, Science347, 1462 (2015)

    2015

  44. [44]

    M. J. Jee, I. Khabibullin, and A. Vikhlinin, arXiv:2605.00093 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  45. [45]

    ATLAS Collaboration, Nature628, 8009 (2024)

    2024

  46. [46]

    CMS Collaboration, Phys. Rev. D109, 012011 (2024). 17

    2024