pith. sign in

arxiv: 2604.14688 · v1 · submitted 2026-04-16 · ✦ hep-ph · astro-ph.CO

Exploring non-equilibrium effects in sequential freeze-in

Shiuli Chatterjee , Andrzej Hryczuk This is my paper

Pith reviewed 2026-05-10 11:18 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.CO
keywords dark matterfreeze-innon-equilibriumphase-spacerelic abundancemulti-component dark sectortwo-scalar modelnon-thermal evolution
0
0 comments X

The pith

In a two-scalar dark sector, non-equilibrium evolution shifts the dark matter relic abundance by up to an order of magnitude from standard calculations.

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

The paper examines freeze-in production in multi-component dark sectors where the dark sector itself may not reach local thermal equilibrium. In a minimal two-scalar model chosen for its potential observable signatures, the authors solve the full phase-space evolution instead of the usual number-density equations. They report that the resulting dark matter abundance can differ by as much as a factor of ten from the traditional approach in regions consistent with current constraints. This matters because relic density predictions directly affect expected rates in indirect detection and long-lived particle searches.

Core claim

Focusing on the phenomenologically viable regions, we analyse the impact of non-thermal evolution on the dark matter abundance, finding deviations of up to an order of magnitude between the full phase-space treatment and the traditional number-density approach.

What carries the argument

Numerical solution of the phase-space Boltzmann equations for the two-scalar model, tracking non-thermal momentum distributions during sequential freeze-in.

If this is right

  • Relic density calculations for multi-particle dark sectors must move beyond integrated number densities to full momentum distributions.
  • Parameter spaces derived from number-density freeze-in can mis-predict indirect detection fluxes by up to an order of magnitude.
  • Forward physics searches for long-lived particles will require adjusted production rates once non-equilibrium effects are included.
  • Dedicated phase-space codes become necessary for any dark sector with sequential production channels.

Where Pith is reading between the lines

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

  • Similar order-of-magnitude shifts may appear in other freeze-in scenarios involving multiple dark states, even if the paper only studies one model.
  • Existing cosmological bounds on dark matter from structure formation could tighten or loosen once non-thermal momentum distributions are properly accounted for.
  • The result motivates building general-purpose phase-space solvers that can handle arbitrary dark sector interactions beyond the two-scalar case.

Load-bearing premise

The chosen minimal two-scalar model and its numerical phase-space implementation faithfully capture the non-equilibrium dynamics without hidden fitting or uncontrolled approximations.

What would settle it

Re-running the phase-space solver for the same benchmark points and finding that the dark matter yield differs from the number-density result by less than a factor of a few.

read the original abstract

Freeze-in of multi-component dark sectors is governed not only by the interaction with the thermal plasma, but also by their internal dynamics. Full thermalisation within the dark sector is not guaranteed, raising the question of impact of departures from local thermal equilibrium onto the evolution and ultimately relic abundance and momentum distribution of dark matter. In this work we explore this question in a minimal two-scalar model, which can give rise to observable signatures in indirect detection and long-lived particle searches at forward physics experiments. Focusing on the phenomenologically viable regions, we analyse the impact of non-thermal evolution on the dark matter abundance, finding deviations of up to an order of magnitude between the full phase-space treatment and the traditional number-density approach. Our results highlight the importance of phase-space level computation for accurate freeze-in predictions and further motivate dedicated numerical tools for studying the evolution of multi-component dark sectors at the phase space level.

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

2 major / 2 minor

Summary. The paper studies non-equilibrium dynamics in sequential freeze-in of a multi-component dark sector using a minimal two-scalar model. It solves the full momentum-dependent Boltzmann equations for the phase-space distributions and compares the resulting dark matter relic abundance to the standard integrated number-density approach, reporting deviations of up to an order of magnitude in phenomenologically viable regions. The work argues that these differences arise from non-thermal evolution within the dark sector and calls for dedicated phase-space numerical tools.

Significance. If the reported deviations are confirmed to be physical rather than numerical artifacts, the result would demonstrate that the number-density approximation can fail at the order-of-magnitude level for sequential freeze-in scenarios, with direct consequences for indirect detection signals and long-lived particle searches. The choice of a concrete, minimal model and the quantitative comparison provide a useful benchmark for future phase-space studies in dark matter cosmology.

major comments (2)
  1. [§3.2] §3.2 (Numerical Implementation): No convergence tests are presented for the momentum-grid resolution, time-stepping accuracy, or collision-integral discretization in the phenomenologically viable parameter regions. The central claim of order-of-magnitude deviations between the full phase-space treatment and the number-density approach is therefore not yet shown to be robust against the numerical method, as the difference can be sensitive to these choices.
  2. [§4] §4 (Results): The comparison plots and tables do not include error estimates, sensitivity scans over numerical parameters, or explicit checks against known thermal limits. Without these, it remains unclear whether the reported deviations are driven by genuine non-equilibrium phase-space effects or by uncontrolled approximations in the solver.
minor comments (2)
  1. [§2] The notation for the interaction rates and distribution functions in §2 could be made more explicit to facilitate reproduction of the Boltzmann equations.
  2. A brief statement on the computational cost and typical runtime of the phase-space solver would help readers assess the practicality of the method.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive feedback on numerical robustness. We address the major comments point by point below. We agree that additional validation is needed and will revise the manuscript to incorporate the requested tests and checks.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (Numerical Implementation): No convergence tests are presented for the momentum-grid resolution, time-stepping accuracy, or collision-integral discretization in the phenomenologically viable parameter regions. The central claim of order-of-magnitude deviations between the full phase-space treatment and the number-density approach is therefore not yet shown to be robust against the numerical method, as the difference can be sensitive to these choices.

    Authors: We agree that explicit convergence tests were not included in the original submission and that they are necessary to substantiate the robustness of the results. In the revised manuscript we will add a new subsection (or appendix) that presents convergence studies with respect to momentum-grid resolution, time-stepping accuracy, and collision-integral discretization, performed specifically in the phenomenologically viable parameter regions. These tests will demonstrate that the reported order-of-magnitude deviations remain stable under refinement of the numerical parameters. revision: yes

  2. Referee: [§4] §4 (Results): The comparison plots and tables do not include error estimates, sensitivity scans over numerical parameters, or explicit checks against known thermal limits. Without these, it remains unclear whether the reported deviations are driven by genuine non-equilibrium phase-space effects or by uncontrolled approximations in the solver.

    Authors: We acknowledge that the original Section 4 lacks error estimates, sensitivity scans, and explicit thermal-limit checks. We will revise the results section to include (i) error estimates on the computed relic abundances, (ii) sensitivity scans over the principal numerical parameters, and (iii) direct comparisons to the expected thermal-equilibrium limits in appropriate regimes. These additions will confirm that the observed deviations originate from non-equilibrium phase-space evolution rather than numerical artifacts. revision: yes

Circularity Check

0 steps flagged

No circularity: numerical comparison of Boltzmann solvers in a minimal model

full rationale

The paper performs a numerical exploration of non-equilibrium dynamics in a two-scalar sequential freeze-in scenario by solving the full phase-space Boltzmann equations and comparing the resulting dark-matter abundance to the integrated number-density approximation. No equations, parameters, or results are shown to be defined in terms of themselves; the claimed order-of-magnitude deviations are outputs of the solver rather than inputs, and no self-citation chain or ansatz smuggling is present in the abstract or description. The derivation chain is therefore self-contained and independent of the target observable.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review prevents identification of specific free parameters, axioms, or invented entities; the central claim rests on an unstated numerical implementation of phase-space evolution whose assumptions are not visible.

pith-pipeline@v0.9.0 · 5445 in / 1137 out tokens · 21686 ms · 2026-05-10T11:18:28.286593+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. 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.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages · cited by 1 Pith paper · 3 internal anchors

  1. [1]

    Resonance enhancement of dark matter interactions: the case for early kinetic decoupling and velocity dependent resonance width

    M. Duch and B. Grzadkowski,Resonance enhancement of dark matter interactions: the case for early kinetic decoupling and velocity dependent resonance width,JHEP09(2017) 159 [1705.10777]

    work page Pith review arXiv 2017

  2. [2]

    Binder, T

    T. Binder, T. Bringmann, M. Gustafsson and A. Hryczuk,Early kinetic decoupling of dark matter: when the standard way of calculating the thermal relic density fails,Phys. Rev. D96 (2017) 115010 [1706.07433]

    work page arXiv 2017

  3. [3]

    Binder, T

    T. Binder, T. Bringmann, M. Gustafsson and A. Hryczuk,Dark matter relic abundance beyond kinetic equilibrium,Eur. Phys. J. C81(2021) 577 [2103.01944]

    work page arXiv 2021

  4. [4]

    Self-heating dark matter via semi-annihilation

    A. Kamada, H.J. Kim, H. Kim and T. Sekiguchi,Self-Heating Dark Matter via Semiannihilation,Phys. Rev. Lett.120(2018) 131802 [1707.09238]

    work page Pith review arXiv 2018

  5. [5]

    Cai and A

    Y. Cai and A. Spray,Low-Temperature Enhancement of Semi-annihilation and the AMS-02 Positron Anomaly,JHEP10(2018) 075 [1807.00832]

    work page arXiv 2018

  6. [6]

    Improved bounds on $\mathbb{Z}_{3}$ singlet dark matter

    A. Hektor, A. Hryczuk and K. Kannike,Improved bounds onZ 3 singlet dark matter,JHEP 03(2019) 204 [1901.08074]

    work page Pith review arXiv 2019

  7. [7]

    Hryczuk and M

    A. Hryczuk and M. Laletin,Impact of dark matter self-scattering on its relic abundance, Phys. Rev. D106(2022) 023007 [2204.07078]

    work page arXiv 2022

  8. [8]

    Chatterjee and A

    S. Chatterjee and A. Hryczuk,Conversions in two-component dark sectors: a phase space level analysis,JHEP07(2025) 279 [2502.08725]

    work page arXiv 2025

  9. [9]

    Dienes, F

    K.R. Dienes, F. Huang, J. Kost, S. Su and B. Thomas,Deciphering the archaeological record: Cosmological imprints of nonminimal dark sectors,Phys. Rev. D101(2020) 123511 [2001.02193]

    work page arXiv 2020

  10. [10]

    Beauchesne and C.-W

    H. Beauchesne and C.-W. Chiang,Impact of non-thermal phase-space distributions on dark matter abundance in secluded sectors,Eur. Phys. J. C84(2024) 950 [2401.03657]

    work page arXiv 2024

  11. [11]

    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,JCAP04(2022) 012 [2111.01267]

    work page arXiv 2022

  12. [12]

    Bringmann, P

    T. Bringmann, P.F. Depta, M. Hufnagel, J.T. Ruderman and K. Schmidt-Hoberg,Dark Matter from Exponential Growth,Phys. Rev. Lett.127(2021) 191802 [2103.16572]

    work page arXiv 2021

  13. [13]

    Hryczuk and M

    A. Hryczuk and M. Laletin,Dark matter freeze-in from semi-production,JHEP06(2021) 026 [2104.05684]. – 24 –

    work page arXiv 2021

  14. [14]

    Bhatia,A closer look at dark matter production in exponential growth scenarios,JHEP 04(2025) 185 [2308.09801]

    D. Bhatia,A closer look at dark matter production in exponential growth scenarios,JHEP 04(2025) 185 [2308.09801]

    work page arXiv 2025

  15. [15]

    Warm dark matter from freeze-in at stronger coupling

    D. Feiteira, O. Lebedev and V. Oliveira,Warm dark matter from freeze-in at stronger coupling,2602.20242

    work page internal anchor Pith review Pith/arXiv arXiv

  16. [16]

    Sakurai and W

    K. Sakurai and W. Yin,Stimulated emission of dark matter via thermal scattering: novel limits for freeze-in and eV cold dark matter,JHEP03(2025) 202 [2410.18968]

    work page arXiv 2025

  17. [17]

    Seasons of Dark Matter Freeze-In Shaped by the Weather of the Early Universe

    F. D’Eramo, A. Lenoci and T. Sassi,Seasons of dark matter freeze-in shaped by the weather of the early Universe,Phys. Rev. D113(2026) 083502 [2511.07511]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  18. [18]

    D’Eramo and A

    F. D’Eramo and A. Lenoci,Lower mass bounds on FIMP dark matter produced via freeze-in, JCAP10(2021) 045 [2012.01446]

    work page arXiv 2021

  19. [19]

    Decant, J

    Q. Decant, J. Heisig, D.C. Hooper and L. Lopez-Honorez,Lyman-αconstraints on freeze-in and superWIMPs,JCAP03(2022) 041 [2111.09321]

    work page arXiv 2022

  20. [20]

    B´ elanger, C

    G. B´ elanger, C. Delaunay, A. Pukhov and B. Zaldivar,Dark matter abundance from the sequential freeze-in mechanism,Phys. Rev. D102(2020) 035017 [2005.06294]

    work page arXiv 2020

  21. [21]

    Badziak and M

    M. Badziak and M. Laletin,Precise predictions for the QCD axion contribution to dark radiation with full phase-space evolution,JHEP02(2025) 108 [2410.18186]

    work page arXiv 2025

  22. [22]

    K.J. Bae, A. Kamada, S.P. Liew and K. Yanagi,Light axinos from freeze-in: production processes, phase space distributions, and Ly-αforest constraints,JCAP01(2018) 054 [1707.06418]

    work page arXiv 2018

  23. [23]

    Dark matter from dark photons: a taxonomy of dark matter production,

    T. Hambye, M.H.G. Tytgat, J. Vandecasteele and L. Vanderheyden,Dark matter from dark photons: a taxonomy of dark matter production,Phys. Rev. D100(2019) 095018 [1908.09864]

    work page arXiv 2019

  24. [24]

    Siqueira,Secluded Dark Matter in light of the Cherenkov Telescope Array (CTA),Phys

    C. Siqueira,Secluded Dark Matter in light of the Cherenkov Telescope Array (CTA),Phys. Lett. B797(2019) 134840 [1901.11055]

    work page arXiv 2019

  25. [25]

    Heikinheimo, T

    M. Heikinheimo, T. Tenkanen and K. Tuominen,Prospects for indirect detection of frozen-in dark matter,Phys. Rev. D97(2018) 063002 [1801.03089]

    work page arXiv 2018

  26. [26]

    Illuminating sequential freeze-in dark matter with dark photon signal at the CERN SHiP experiment

    X. Yin and S. Zheng,Illuminating sequential freeze-in dark matter with dark photon signal at the CERN SHiP experiment,Phys. Rev. D113(2026) 075013 [2512.05380]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  27. [27]

    Anchordoqui et al.,The forward physics facility: Physics opportunities and conceptual design,Nucl

    L.A. Anchordoqui et al.,The forward physics facility: Physics opportunities and conceptual design,Nucl. Phys. B1026(2026) 117398

    work page 2026

  28. [28]

    The Higgs portal to cosmology,

    O. Lebedev,The Higgs portal to cosmology,Prog. Part. Nucl. Phys.120(2021) 103881 [2104.03342]

    work page arXiv 2021

  29. [29]

    Freeze-in production of decaying dark matter in five steps

    S. Heeba, F. Kahlhoefer and P. St¨ ocker,Freeze-in production of decaying dark matter in five steps,JCAP11(2018) 048 [1809.04849]. [31]MAGIC, Fermi-LATcollaboration,Limits to Dark Matter Annihilation Cross-Section from a Combined Analysis of MAGIC and Fermi-LAT Observations of Dwarf Satellite Galaxies, JCAP02(2016) 039 [1601.06590]. [32]CTAcollaboration,S...

    work page Pith review arXiv 2018

  30. [30]

    M. Wood, J. Buckley, S. Digel, S. Funk, D. Nieto and M.A. Sanchez-Conde,Prospects for Indirect Detection of Dark Matter with CTA, inSnowmass 2013: Snowmass on the Mississippi, 5, 2013 [1305.0302]. – 25 –

    work page arXiv 2013

  31. [31]

    Minimal Extension of the Standard Model Scalar Sector

    D. O’Connell, M.J. Ramsey-Musolf and M.B. Wise,Minimal Extension of the Standard Model Scalar Sector,Phys. Rev.D75(2007) 037701 [hep-ph/0611014]. [35]Fermi-LAT, HA WC, H.E.S.S., MAGIC, VERITAScollaboration,Combined dark matter search towards dwarf spheroidal galaxies with Fermi-LAT, HAWC, H.E.S.S., MAGIC, and VERITAS,2508.20229. [36]CHARMcollaboration,Sea...

    work page Pith review arXiv 2007

  32. [32]

    Decay and Detection of a Light Scalar Boson Mixing with the Higgs

    M.W. Winkler,Decay and detection of a light scalar boson mixing with the Higgs boson, Phys. Rev. D99(2019) 015018 [1809.01876]. [41]F ASERcollaboration,FASER’s physics reach for long-lived particles,Phys. Rev. D99 (2019) 095011 [1811.12522]. [42]SHiPcollaboration,A facility to Search for Hidden Particles (SHiP) at the CERN SPS, 1504.04956. [43]MATHUSLAcol...

    work page Pith review arXiv 2019

  33. [33]

    Lanfranchi, M

    G. Lanfranchi, M. Pospelov and P. Schuster,The Search for Feebly-Interacting Particles, 2011.02157

    work page arXiv 2011

  34. [34]

    BBN for the LHC: constraints on lifetimes of the Higgs portal scalars

    A. Fradette and M. Pospelov,BBN for the LHC: constraints on lifetimes of the Higgs portal scalars,Phys. Rev.D96(2017) 075033 [1706.01920]

    work page Pith review arXiv 2017

  35. [35]

    Probing Light Thermal Dark-Matter With a Higgs Portal Mediator

    G. Krnjaic,Probing Light Thermal Dark-Matter With a Higgs Portal Mediator,Phys. Rev. D94(2016) 073009 [1512.04119]

    work page Pith review arXiv 2016

  36. [36]

    Cholis, Y.-M

    I. Cholis, Y.-M. Zhong, S.D. McDermott and J.P. Surdutovich,Return of the templates: Revisiting the Galactic Center excess with multimessenger observations,Phys. Rev. D105 (2022) 103023 [2112.09706]

    work page arXiv 2022

  37. [37]

    PPPC 4 DM ID: A Poor Particle Physicist Cookbook for Dark Matter Indirect Detection

    M. Cirelli, G. Corcella, A. Hektor, G. Hutsi, M. Kadastik, P. Panci et al.,PPPC 4 DM ID: A Poor Particle Physicist Cookbook for Dark Matter Indirect Detection,JCAP03(2011) 051 [1012.4515]

    work page Pith review arXiv 2011

  38. [38]

    Model-Independent Indirect Detection Constraints on Hidden Sector Dark Matter

    G. Elor, N.L. Rodd, T.R. Slatyer and W. Xue,Model-Independent Indirect Detection Constraints on Hidden Sector Dark Matter,JCAP06(2016) 024 [1511.08787]

    work page Pith review arXiv 2016

  39. [39]

    HDECAY: a Program for Higgs Boson Decays in the Standard Model and its Supersymmetric Extension

    A. Djouadi, J. Kalinowski and M. Spira,HDECAY: A Program for Higgs boson decays in the standard model and its supersymmetric extension,Comput. Phys. Commun.108(1998) 56 [hep-ph/9704448]

    work page Pith review arXiv 1998

  40. [40]

    Bringmann and S

    T. Bringmann and S. Hofmann,Thermal decoupling of WIMPs from first principles,JCAP 04(2007) 016 [hep-ph/0612238]. – 26 –

    work page arXiv 2007

  41. [41]

    Matter Power Spectrum in Hidden Neutrino Interacting Dark Matter Models: A Closer Look at the Collision Term

    T. Binder, L. Covi, A. Kamada, H. Murayama, T. Takahashi and N. Yoshida,Matter Power Spectrum in Hidden Neutrino Interacting Dark Matter Models: A Closer Look at the Collision Term,JCAP11(2016) 043 [1602.07624]

    work page Pith review arXiv 2016

  42. [42]

    Aboubrahim, M

    A. Aboubrahim, M. Klasen and L.P. Wiggering,Forbidden dark matter annihilation into leptons with full collision terms,JCAP08(2023) 075 [2306.07753]

    work page arXiv 2023

  43. [43]

    Yoon,Boltzmann Equation Solver for Thermalization,2603.28848

    J.-H. Yoon,Boltzmann Equation Solver for Thermalization,2603.28848. – 27 –

    work page arXiv