Pith. sign in

REVIEW 3 major objections 5 minor 4 cited by

Dark matter produced by late dark-sector decay must be heavier than about 0.1 GeV; combining Lyman-alpha forest and Big Bang nucleosynthesis bounds rules out lighter candidates.

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 · deepseek-v4-flash

2026-08-02 17:29 UTC pith:NBPKEFEC

load-bearing objection Careful, useful constraint on a dark-sector decay model; the thermal-WDM mapping is credible, but the paper should show validation across its actual parameter grid. the 3 major comments →

arxiv 2603.24331 v3 pith:NBPKEFEC submitted 2026-03-25 hep-ph astro-ph.CO

Lyman-α Forest Constraint on Dark Matter from Dark Sector Decay

Si-Yuan Zhao , Yi-Cheng Dai , Wei Liao , Yi-Song Lu This is my paper
classification hep-ph astro-ph.CO MSC 83F0585A40 PACS 95.35.+d98.80.-k
keywords dark matterLyman-alpha forestwarm dark matternon-thermal dark matterdark sector decaysub-GeV dark matterBig Bang nucleosynthesismatter power spectrum
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.

This paper argues that a specific non-thermal dark matter production mechanism—dark matter generated by the late-time decay of a heavier dark-sector particle—is almost fully closed off for sub-GeV candidates. By computing the momentum distribution of the decay-produced particles and mapping it to an equivalent thermal warm dark matter model, the authors translate existing Lyman-alpha forest limits into a direct constraint on the dark matter mass. Combined with Big Bang nucleosynthesis bounds on the parent particle's lifetime, they conclude that dark matter below roughly 0.1 GeV is excluded, and that even 0.01 GeV is ruled out under conservative assumptions. A reader should care because non-thermal production is a common way to evade WIMP search limits, and this paper shows small-scale structure can test it.

Core claim

The paper's central claim is that in the dark sector decay model, where a WIMP-like scalar φ freezes out and then decays through a right-handed neutrino into the dark matter particle χ, the resulting non-thermal χ particles are too hot to be the dominant dark matter if they are light. Using Boltzmann equations to evolve the phase-space distribution through the decay, the authors find that the distribution is single-peaked and can be matched to a thermal warm dark matter distribution by equating comoving number density and mean squared comoving momentum. The effective warm dark matter mass obtained this way is compared against the latest Lyman-alpha forest bound; the fiducial 5.7 keV limit ex

What carries the argument

The equivalent thermal warm dark matter (WDM) approximation is the load-bearing tool. The non-thermal decay distribution F(q) is replaced by a Fermi-Dirac distribution with the same comoving number density and the same mean squared comoving momentum; solving the two matching equations yields an effective temperature T_eff and effective mass m_eff. Because the Lyman-alpha forest has published lower bounds on thermal WDM mass, m_eff can be tested directly. A second, independent route fits the numerical transfer function to the standard WDM fitting form and inverts the mass–break-scale relation, giving consistent boundaries.

Load-bearing premise

The load-bearing premise is that matching only the number density and mean squared comoving momentum of the non-thermal distribution to a thermal warm dark matter distribution reproduces the small-scale power suppression everywhere; this was spot-checked for only a handful of parameter points.

What would settle it

Compute the exact linear matter power spectrum for the full non-thermal distribution at parameter points just inside and outside the exclusion boundary—for example mχ = 0.08 GeV, m_φ = 1000 GeV, m_N = 1 GeV, y_DS near 2×10^-12—and compare the half-mode scale (or the full T(k)) with the equivalent-WDM prediction. If the two disagree by more than the quoted ~2%, the exclusion boundary would shift and the central claim would need revision.

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

If this is right

  • Dark matter from this decay scenario cannot be the dominant component below ~0.1 GeV: the strongest Lyman-alpha bound excludes mχ ≲ 0.08 GeV, and even the weakest excludes mχ ≲ 0.01 GeV.
  • The joint Lyman-alpha plus BBN constraints require the decay coupling y_DS to be at least about 10^-12 for sub-GeV dark matter, because the parent must decay early enough to redshift the hot daughters.
  • Lyman-alpha data place an upper limit on the parent mass m_φ; for example, for mχ = 0.1 GeV the ceiling drops to about 500 GeV, and for mχ = 0.01 GeV the allowed window closes.
  • Highly degenerate parent–neutrino spectra (m_N/m_φ close to 1) are excluded by BBN for m_φ = 1000, 100, and 10 GeV once the ratio exceeds thresholds of 0.945, 0.810, and 0.001, respectively.
  • The mapping method is generic: any single-peaked non-thermal DM distribution can be translated into an effective WDM mass and constrained by the same published Lyman-alpha limits.

Where Pith is reading between the lines

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

  • The same two-moment mapping could be applied to other late-decay production models (for instance, gravitino or axino superWIMPs), provided their distributions stay single-peaked; the paper's benchmark checks do not cover all those cases.
  • Because the paper fixes χ to account for all dark matter, a subdominant χ component would evade the mass bound; the exclusion should be viewed as applying to the dominant-DM interpretation.
  • Tighter future Lyman-alpha bounds (or 21-cm observations of the epoch of reionization) would push the lower mass limit upward approximately as m_eff ∝ ⟨v²⟩^{-3/8}, so the threshold scales slowly and substantial improvement requires significantly stronger limits.
  • A full test of the mapping near the boundary using a non-linear or N-body method, rather than the linear transfer function, would strengthen the claim; the five benchmark points do not exhaust parameter space.

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

3 major / 5 minor

Summary. The paper studies a non-thermal dark matter model in which the DM particle χ is produced by the late decay of a frozen-out dark scalar ϕ → N χ. The authors solve the Boltzmann equations for the phase-space distribution of χ, including a zero-momentum approximation for the decaying parent, and derive the resulting linear matter power spectrum. They then map the non-thermal χ distribution to an equivalent thermal WDM distribution by matching number density and mean squared comoving momentum (Eqs. 3.2–3.3), obtaining an effective WDM mass m_eff. This m_eff is compared to recent Lyman-α forest lower bounds on thermal WDM mass, and the results are combined with BBN lifetime limits. The main claim is that for the decay production scenario, DM masses m_χ ≲ 0.08 GeV are excluded, with even conservative limits excluding m_χ ≲ 0.01 GeV.

Significance. If the main result holds, the paper provides a strong and phenomenologically useful constraint on a well-motivated class of non-thermal dark matter models, showing that small-scale structure probes are complementary to BBN and can exclude sub-GeV DM from dark-sector decay. The paper has notable strengths: the Boltzmann collision terms are derived carefully, the comoving number conservation is checked explicitly, and the mapping to thermal WDM is cross-validated against CLASS transfer functions at several benchmark points. The two independent routes to the effective mass—velocity-dispersion matching and transfer-function fitting—agree at the displayed benchmarks, which is a good internal consistency check. The analysis is not circular: the Lyman-α and BBN limits are external observables. The main weakness is that the mapping is explicitly validated only for a handful of benchmark points, while the headline exclusion boundary lies outside those points; the asserted 2% accuracy across the parameter space is not exhibited.

major comments (3)
  1. [§3.2, Tables 1–2] The central claim rests on the equivalence mapping between the non-thermal χ distribution and a thermal WDM distribution. The CLASS validations are limited to m_ϕ = 1000 GeV, y_DS = 10^-12, m_N = 1 GeV, with m_χ = 0.01, 1, 10 GeV (Table 1) and m_χ = 0.62, 1.33 GeV (Table 2). The headline exclusion m_χ ≲ 0.08 GeV (Figs. 7–9) is not tested, and the text's statement that 'analogous calculations across a broader parameter space' remain within 2% is not supported by any displayed calculation. Since the boundary is precisely where the mapping error matters, please provide explicit α_χ vs α_WDM comparisons for points on the m_eff = 5.7 keV and 3.2 keV exclusion curves, including different m_ϕ and y_DS, or a dense scan of the relative error. Without this, the 0.08 GeV bound is an extrapolation.
  2. [Abstract and §5] The abstract states unconditionally that 'masses lighter than ~10^-1 GeV are excluded', but §5 correctly qualifies this as applying to the decay production scenario. The model also has a freeze-in production component, which becomes non-negligible for y_DS ≳ 2×10^-12 and is colder than the decay component. For larger y_DS the Lyα constraint is therefore weaker, and the full model may admit lighter m_χ. The abstract and conclusions should carry the y_DS ≲ 2×10^-12 qualification so that the headline does not overstate the scope of the bound.
  3. [§4, BBN limits] The BBN constraints are imposed as upper limits on the ϕ lifetime (τ ≤ 0.7 s or 0.07 s), but the decay chain is ϕ → N χ, and the energy that affects BBN is injected when N decays to SM particles through the seesaw Yukawa, not necessarily at the time of ϕ decay. The paper does not specify the N lifetime or the relevant Yukawa couplings, nor does it justify that the adopted τ limits are conservative with respect to the additional delay. Since the combined exclusion around m_χ ≈ 0.08 GeV uses these BBN lines, please clarify that the bounds are correctly applied or cite the analysis in Ref. [8] where this is established.
minor comments (5)
  1. [Eq. (2.12)] The kinematic integration limits p_ϕ± are not written explicitly in the main text. Please state them in §2 or refer to Appendix A at first use, as the reader currently has to reconstruct them from the appendix.
  2. [§3.1 / Appendix B.2] The quantity ⟨v²⟩ is called the comoving pseudo-velocity dispersion but is defined as ⟨q²⟩_χ/m_χ². It would help to state explicitly that the physical velocity dispersion at scale factor a is a^{-2}⟨q²⟩_χ/m_χ² and that this is why comoving matching is sufficient.
  3. [§3.2, Eq. (3.5)] The fitting exponent µ = 1.12 is used without explanation. Please note that this is the standard value for a fermionic thermal relic and cite the original determination, as this matters for the interpretation of α and m_eff.
  4. [Table 3] Table 3 lists several observational limits (Milky Way satellite counts, stellar streams) that are not used in the analysis. It would improve clarity to state explicitly in the text that only the Lyman-α and UV luminosity function limits are used in the figures.
  5. [§3, Fig. 4] Figure 4 includes y_DS = 10^-11, which is outside the decay-production regime as defined in the paper. The caption notes this, but the main text could repeat it to avoid readers interpreting this benchmark as one that can be constrained by the presented BBN/Lyα bounds.

Circularity Check

0 steps flagged

No significant circularity: the Lyman-α and BBN constraints are external inputs, and the model-to-WDM mapping is an approximation validated by independent CLASS calculations, not an input recycled as a prediction.

full rationale

The paper's central claim—that sub-GeV dark matter from dark-sector decay is excluded by Lyman-α forest data—rests on a derivation chain that does not reduce to its own inputs. The model Lagrangian (eqs. 2.1–2.4) is defined independently; the phase-space distribution is computed by solving the Boltzmann equations (2.15–2.16); the mapping to an equivalent thermal WDM is a definitional matching of number density and mean-squared comoving momentum (eqs. 3.2–3.3), but the resulting effective WDM mass is then compared to externally derived Lyman-α lower bounds (e.g., m_WDM ≥ 5.7 keV from [16]). The mapping is not merely asserted: it is tested against the independent Boltzmann code CLASS in Figs. 5–6 and quantified in Tables 1–2. The alternative transfer-function method fits the model's transfer function and inverts the published analytical WDM formula (eq. 3.6), again referencing external results. BBN constraints are likewise imported from standard cosmology. The only notable self-citation is [8], the same group's earlier model paper, which supplies the model's starting assumptions and the freeze-in subdominance bound; these are not the target prediction, and the present paper recomputes the non-thermal distributions itself. The limited number of benchmark validations (five points, with the exclusion boundary near mχ ≈ 0.08 GeV not directly tested) is a legitimate robustness concern, but it is a matter of numerical verification, not circularity: no equation sets the predicted observable equal to a fitted parameter, and no part of the conclusion is assumed in the inputs. Therefore the analysis is self-contained against external observables and no circular step is present.

Axiom & Free-Parameter Ledger

5 free parameters · 6 axioms · 0 invented entities

No new particle or force is invented; the model is from ref. [8]. The main un-derived input is the equivalence-to-WDM ansatz plus the external cosmological limits. The four scan parameters are free model parameters rather than fitted constants; λHϕ is fixed by the relic-density condition.

free parameters (5)
  • mχ (dark matter mass) = scanned 0.01–100 GeV
    Central model parameter; the exclusion boundary is expressed as a lower bound on mχ.
  • mϕ (parent scalar mass) = scanned 10–1000 GeV
    Controls initial injection momentum and the BBN decay-width; one of the four scanned parameters.
  • mN (right-handed neutrino mass) = scanned 1–500 GeV, ratios mN/mϕ up to 0.99
    Affects decay kinematics and phase space; adopted as free in the phenomenology.
  • yDS (decay coupling) = scanned 10^-13–10^-11; decay dominance requires yDS≲2×10^-12
    Sets ϕ lifetime and hence the redshift at which χ is injected; central to the constraints.
  • λHϕ (scalar-Higgs portal) = not scanned; fixed by Ωχ=ΩDM via eq. (2.6)
    Listed for completeness; the paper removes it by requiring χ be all of the dark matter.
axioms (6)
  • domain assumption χ is never in thermal equilibrium and freeze-in production is negligible for yDS ≲ 2×10^-12; late decay of ϕ dominates.
    §2 after eq. (2.16): freeze-in contribution <10% and is 'substantially colder' than decay component, so neglected.
  • domain assumption The parent ϕ thermalizes with the SM bath and freezes out; its subsequent decay determines the χ relic abundance.
    §2, eq. (2.6) fixes λHϕ by Ωχ=ΩDM.
  • domain assumption A thermal WDM model with the same n and ⟨q²⟩ reproduces the small-scale power suppression of the non-thermal decay distribution.
    §3.1 eqs. (3.2)-(3.3); validated only at selected benchmarks in Tables 1-2.
  • domain assumption The BBN constraints on ϕ lifetime (τ≤0.7 s and τ≤0.07 s) are correct for this model.
    §4: adopted from ref. [8]; not rederived here.
  • domain assumption The transfer-function parameterization (3.5) with μ=1.12 and the WDM mass–α relation (3.6) are accurate in the mass range used.
    Eqs. (3.5)–(3.6) taken from refs. [13,17]; claimed <1.5% in 3–5 keV range.
  • domain assumption The zero-momentum approximation for parent decay is accurate.
    Appendix B.1: residual momentum at decay time gives ~0.5% correction; compared in Fig. 14.

pith-pipeline@v1.3.0-alltime-deepseek · 17849 in / 15951 out tokens · 163794 ms · 2026-08-02T17:29:47.593902+00:00 · methodology

0 comments
read the original abstract

By exploiting small-scale structure formation probed by Lyman-$\alpha$ forest observations, we study constraints on a model of dark matter from dark sector decay. We compute the phase space distribution of the dark matter and the linear matter power spectrum. We map the non-thermal dark matter distribution in this dark matter model to an approximate thermal warm dark matter distribution, and use this approximation to obtain a constraint from the Lyman-$\alpha$ forest observation. We combine the latest Lyman-$\alpha$ forest bounds with the constraint from the Big Bang Nucleosynthesis. As these two probes offer highly complementary constraints, we impose strong limits on sub-GeV dark matter. Consequently, masses lighter than $\sim 10^{-1}$ GeV are excluded, thereby significantly limiting the allowed parameter space. More broadly, our findings demonstrate the utility of small-scale structure observations in testing non-thermal dark matter paradigms, offering valuable insights for exploring a wider class of late-time decay models.

discussion (0)

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

Forward citations

Cited by 4 Pith papers

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

  1. Cosmological evolution with decaying dark matter: an integral-equation approach

    astro-ph.CO 2026-07 accept novelty 6.0

    CLASSIER-DDM evaluates generic two-body decaying dark matter perturbations via iterative integral equations at O(0.1%) accuracy and O(1 min) cost without a Boltzmann hierarchy or fluid approximation.

  2. KineticXGPU: A Tensorized Collision Operator for Dark-Sector Self-Scattering

    hep-ph 2026-07 unverdicted novelty 6.0

    Presents a tensorized GPU implementation of the 2-to-2 elastic self-collision operator for dark-sector particles and applies it to a two-source freeze-in scenario where self-interactions erase bimodal features.

  3. Lyman-$\alpha$ forest constraints on pure and mixed fuzzy dark matter

    astro-ph.CO 2026-06 unverdicted novelty 5.0

    Lyman-alpha forest data yield m_FDM > 1.9e-21 eV (95% CL) for pure FDM and f_FDM upper limits of 0.07-0.65 for mixed FDM at log10(m_FDM/eV) = -23 to -21.

  4. Machine Learning Does It and Does It Better: Unearthing Primordial Dark-Matter Velocities from the Matter Power Spectrum

    astro-ph.CO 2026-06 unverdicted novelty 4.0

    A 1D convolutional neural network reconstructs the dark-matter phase-space distribution from the matter power spectrum with greater accuracy and broader applicability than an earlier empirical formula.

Reference graph

Works this paper leans on

19 extracted references · 13 linked inside Pith · cited by 4 Pith papers

  1. [7]

    J.L. Feng, A. Rajaraman and F. Takayama,Superweakly interacting massive particles,Phys. Rev. Lett.91(2003) 011302 [hep-ph/0302215]

    Pith/arXiv arXiv 2003

  2. [8]

    Cheng and W

    Y. Cheng and W. Liao,Light dark matter from dark sector decay,Phys. Lett. B815(2021) 136118

    2021

  3. [9]

    Cheng, W

    Y. Cheng, W. Liao and Q.-S. Yan,Collider search of light dark matter model with dark sector decay,Chin. Phys. C46(2022) 063103. – 22 –

    2022

  4. [10]

    Liu, F.-L

    A. Liu, F.-L. Shao, Z.-L. Han, Y. Jin and H. Li,Constraints on sterile neutrino portal dark matter from cosmology and collider,Phys. Rev. D108(2023) 115028 [2212.10043]

    Pith/arXiv arXiv 2023

  5. [11]

    Ghosh, P

    D.K. Ghosh, P. Ghosh and S. Jeesun,CMB signature of non-thermal Dark Matter produced from self-interacting dark sector,JCAP01(2024) 023 [2301.13754]

    Pith/arXiv arXiv 2024

  6. [12]

    Liu, F.-L

    A. Liu, F.-L. Shao, Z.-L. Han, Y. Jin and H. Li,Sterile neutrino portal dark matter in νTHDM,Eur. Phys. J. C83(2023) 415

    2023

  7. [13]

    Non-cold

    R. Murgia, A. Merle, M. Viel, M. Totzauer and A. Schneider,“Non-cold” dark matter at small scales: a general approach,JCAP11(2017) 046 [1704.07838]

    Pith/arXiv arXiv 2017

  8. [14]

    Doroshkevich and M.Y

    A.G. Doroshkevich and M.Y. Khlopov,Formation of structure in the Universe with unstable neutrinos,Mon. Not. Roy. Astron. Soc.211(1984) 279

    1984

  9. [15]

    M. Viel, J. Lesgourgues, M.G. Haehnelt, S. Matarrese and A. Riotto,Constraining warm dark matter candidates including sterile neutrinos and light gravitinos with WMAP and the Lyman-alpha forest,Phys. Rev. D71(2005) 063534 [astro-ph/0501562]

    Pith/arXiv arXiv 2005

  10. [16]

    Iršič et al.,Unveiling dark matter free-streaming at the smallest scales with the high redshift Lyman-αforest,Phys

    V. Iršič et al.,Unveiling dark matter free-streaming at the smallest scales with the high redshift Lyman-αforest,Phys. Rev. D109(2024) 043511 [2309.04533]

    Pith/arXiv arXiv 2024

  11. [17]

    Decant, J

    Q. Decant, J. Heisig, M.E. Krauss and T. Ma,Lyman-αconstraints on freeze-in and superWIMPs,JCAP03(2022) 041 [2111.09321]

    Pith/arXiv arXiv 2022

  12. [18]

    Y. Du, F. Huang, H.-L. Li, Y.-Z. Li and J.-H. Yu,Revisiting Dark Matter Freeze-in and Freeze-out through Phase-Space Distribution,JHEP04(2022) 012 [2111.01267]

    Pith/arXiv arXiv 2022

  13. [19]

    Dodelson and F

    S. Dodelson and F. Schmidt,Modern Cosmology, Academic Press, Amsterdam, 2nd ed. (2020)

    2020

  14. [20]

    D. Blas, J. Lesgourgues and T. Tram,The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes,JCAP07(2011) 034 [1104.2933]

    Pith/arXiv arXiv 2011

  15. [21]

    Lesgourgues,The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview, 1104.2932

    J. Lesgourgues,The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview, 1104.2932

    Pith/arXiv arXiv

  16. [22]

    Nadler, R

    E.O. Nadler, R. An, V. Gluscevic, A. Benson and X. Du,COZMIC. I. Cosmological Zoom-in Simulations with Initial Conditions Beyond Cold Dark Matter,Astrophys. J.986(2025) 127 [2410.03635]

    Pith/arXiv arXiv 2025

  17. [23]

    J. Liu, Y. Gong and K. Liao,Joint Constraints on Fuzzy and Warm Dark Matter from Satellite Populations of the Milky Way and Andromeda,Astrophys. J.(2026) [2512.01361]

    arXiv 2026

  18. [24]

    Banik, J

    N. Banik, J. Bovy, G. Bertone, D. Erkal and T.J.L. de Boer,Novel constraints on the particle nature of dark matter from stellar streams,Mon. Not. Roy. Astron. Soc.502(2021) 2364 [1911.02663]

    Pith/arXiv arXiv 2021

  19. [25]

    B. Liu, H. Shan and J. Zhang,New galaxy UV luminosity constraints on warm dark matter from JWST,Astrophys. J.968(2024) 79 [2404.13596]. – 23 –

    Pith/arXiv arXiv 2024