Dynamics of metastable Standard Model particles from long-lived particle decays in the MeV primordial plasma

e investigate the cosmological impact of hypothetical unstable new physics particles that decay in the MeV-scale plasma of the Early Universe. Focusing on scenarios where the decays produce metastable species such as muons, pions, and kaons, we systematically analyze the dynamics of these particles using coupled Boltzmann equations governing their abundances. Our results demonstrate that the metastable species can efficiently annihilate or interact with nucleons, which often leads to their disappearance prior to decay. The suppression of decay significantly alters the properties of cosmic neutrinos, impacting cosmological observables like Big Bang nucleosynthesis and the Cosmic Microwave Background. To support further studies, we provide a public Mathematica code that traces the evolution of these metastable particles and apply it to several new physics models.