-body simulations of cosmologies with light massive relics
Phys. Rev. D 113, 103506 – Published 5 May, 2026
DOI: https://doi.org/10.1103/zmw5-935r
Abstract
The presence of additional relativistic particles at the time of recombination can be inferred through their contribution to . If these species have a finite but low mass (light massive relics—LiMRs), they act as a hot subcomponent of dark matter and impact late-time structure formation. Understanding these effects will be crucial to pin down the underlying particle physics properties of any future detection. While their impact has been well studied on linear scales, this work develops the framework for and presents results from the first set of cosmological N-body simulations that can track the effects of LiMRs, as a function of their mass and temperature, down to fully nonlinear scales. Importantly, our simulations model the impact of both the massive Standard Model neutrinos and LiMRs, which will be crucial in disentangling possible degeneracies. We systematically explore the effects of LiMR properties such as mass, temperature, and initial distribution, on various cosmological observables, including the total matter power spectrum, halo mass functions, mass-concentration relation, radial halo profiles, and weak lensing signals around massive clusters. The presence of LiMRs leads to a suppression of and in the total matter power spectrum for fermionic and bosonic species, respectively. This suppression is measured relative to a neutrino cosmology at for LiMRs fraction () . The framework and simulations developed here will enable detailed followup of the rich phenomenology of LiMR cosmologies.