How new physics affects primordial neutrinos decoupling: Direct Simulation Monte Carlo approach

Abstract

Cosmological observations from Big Bang Nucleosynthesis and the Cosmic Microwave Background (CMB) offer crucial insights into the Early Universe, enabling us to trace its evolution back to lifetimes as short as 0.01 seconds. Upcoming CMB spectrum measurements, such as those underway at the Simons Observatory, will achieve unprecedented precision, allowing for more accurate extraction of information about the properties of the primordial plasma and, in particular, primordial neutrinos. This provides an opportunity to test whether these properties align with the predictions of the standard cosmological model or indicate the presence of new physics that influenced the evolution of the MeV-temperature plasma. A key component in understanding how new physics may have affected primordial neutrinos is solving the neutrino Boltzmann equation. In this paper, we present a novel approach to solving this equation that offers model independence, transparency, and computational efficiency - features that current state-of-the-art methods lack. We demonstrate a proof-of-concept implementation and apply it to several toy scenarios, showcasing key aspects of the primordial plasma's evolution in the presence of new physics.