Neutrino Direct Simulation Monte Carlo: Accurate modeling of out-of-equilibrium decaying cosmological relics
Abstract
We consider nonstandard cosmologies in which decaying relics inject energy at MeV temperatures, driving neutrinos out of equilibrium. The subsequent neutrino evolution shapes predictions for the cosmic radiation density and primordial light element abundances. Earlier work introduced Neutrino Direct Simulation Monte Carlo (νDSMC), which describes neutrino transport through Monte Carlo sampling of particle interactions. Here, we extend this method by jointly evolving the relic population, neutrino spectra, electromagnetic plasma, and cosmic expansion, and coupling this evolution to Big Bang nucleosynthesis. The framework includes direct and inverse decays of the relic, neutrino flavor conversion, quantum statistics, and pion production by energetic neutrinos. To reveal the physical impact of these effects and independently test νDSMC, we compare several approaches to neutrino evolution across a range of nonstandard cosmological scenarios. We find excellent agreement with independent solutions of quantum kinetic or quasi-classical Boltzmann equations. Compared with these calculations, νDSMC includes a broader range of physical processes and is substantially faster for strongly nonthermal neutrino populations. This makes broad, robust parameter scans for nonstandard cosmologies practical. We show that when departures from thermal equilibrium are large, approximations that average over neutrino momenta can give qualitatively incorrect predictions, including the wrong sign of the change in the effective number of neutrino species and substantially incorrect primordial nuclear abundances. We conclude that the extended νDSMC framework enables reliable cosmological tests of new physics by following neutrino thermalization and its observable consequences.
Related
- Cosmological preference for a positive neutrino mass at 2.7σ: A joint analysis of DESI DR2, DESY5, and DESY1 data
- Reanalyzing DESI DR1: IV. Percent-level cosmological constraints from combined probes and robust evidence for the normal neutrino mass hierarchy
- Cosmological Consequences of a Highly Decoupled Right-handed Neutrino
- Cosmological Evolution of Primordial Black Holes: UV/IR Decoupling and the KM3NeT 220 PeV Neutrino Anomaly
- Dark energy perturbations and the robustness of cosmological neutrino-mass constraints
BibTeX
@misc{akita2026neutrino,
title = {{Neutrino Direct Simulation Monte Carlo: Accurate modeling of out-of-equilibrium decaying cosmological relics}},
author = {Kensuke Akita and Miguel Escudero and Oleksii Ihnatenko and Maksym Ovchynnikov},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2609.31833},
eprint = {2609.31833},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.31833},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Neutrino Direct Simulation Monte Carlo: Accurate modeling of out-of-equilibrium decaying cosmological relics
AU - Kensuke Akita
AU - Miguel Escudero
AU - Oleksii Ihnatenko
AU - Maksym Ovchynnikov
PY - 2026
DA - 2026/10/28
PB - arXiv
DO - 10.48550/arxiv.2609.31833
UR - https://arxiv.org/abs/2609.31833
AB - We consider nonstandard cosmologies in which decaying relics inject energy at MeV temperatures, driving neutrinos out of equilibrium. The subsequent neutrino evolution shapes predictions for the cosmic radiation density and primordial light element abundances. Earlier work introduced Neutrino Direct Simulation Monte Carlo (νDSMC), which describes neutrino transport through Monte Carlo sampling of particle interactions. Here, we extend this method by jointly evolving the relic population, neutrino spectra, electromagnetic plasma, and cosmic expansion, and coupling this evolution to Big Bang nucleosynthesis. The framework includes direct and inverse decays of the relic, neutrino flavor conversion, quantum statistics, and pion production by energetic neutrinos. To reveal the physical impact of these effects and independently test νDSMC, we compare several approaches to neutrino evolution across a range of nonstandard cosmological scenarios. We find excellent agreement with independent solutions of quantum kinetic or quasi-classical Boltzmann equations. Compared with these calculations, νDSMC includes a broader range of physical processes and is substantially faster for strongly nonthermal neutrino populations. This makes broad, robust parameter scans for nonstandard cosmologies practical. We show that when departures from thermal equilibrium are large, approximations that average over neutrino momenta can give qualitatively incorrect predictions, including the wrong sign of the change in the effective number of neutrino species and substantially incorrect primordial nuclear abundances. We conclude that the extended νDSMC framework enables reliable cosmological tests of new physics by following neutrino thermalization and its observable consequences.
ER -