Heavy Neutral Leptons (HNLs) are strongly motivated by theory due to their capability of simultaneously explaining the observed phenomena of dark matter, neutrino oscillations and the baryon asymmetry of the Universe. The existence of such particles would affect the expansion history of the Universe and the synthesis of primordial abundances of light elements. In this work we review, revise and extend the phenomenology of HNLs during the Big Bang Nucleosynthesis (BBN) epoch for masses up to 1 GeV. This is of great importance, as BBN is able to provide complementary bounds to those from upcoming and proposed laboratory experiments. To this end we have developed a high-precision Boltzmann code that simulates BBN in the presence of HNLs and takes into account all relevant HNL decay channels, as well as subsequent interactions of decay products (thermalization and decay showers), dilution due to QCD phase transition, active neutrino oscillations and corrections to the weak reaction rates. We present robust bounds on the lifetime and mixing angles of HNLs for masses $3\,\mathrm{MeV}\leq m_N \leq 1\,\mathrm{GeV}$ and show that BBN is able to constrain HNL lifetimes down to $0.03 - 0.05$ s, depending on the mixing pattern. Moreover, combining our results with current experimental searches, we can exclude HNLs that mix purely with electron neutrinos up to ${\sim}$450 MeV and those that mix purely with muon neutrinos up to ${\sim}$360 MeV, in both cases for lifetimes up to at least a few tens of seconds. Finally, we compare the BBN constraints with those obtained from Cosmic Microwave Background observations and explore how our results will be improved by a number of upcoming and proposed laboratory experiments.

Big Bang Nucleosynthesis (BBN) is very sensitive to the cosmological expansion rate. If the gravitational constant $G$ took a different value during the nucleosynthesis epoch than today, the primordial abundances of light elements would be affected. In this work, we improve the bounds on this variation using recent determinations of the primordial element abundances, updated nuclear and weak reaction rates and observations of the Cosmic Microwave Background (CMB). When combining the measured abundances and the baryon density from CMB observations by Planck, we find $G_\mathrm{BBN}/G_0 = 0.99^{+0.06}_{-0.05}$ at $2\sigma$ confidence level. If the variation of $G$ is linear in time, we find $\dot{G}/G_0 = 0.7^{+3.8}_{-4.3}\times 10^{-12} \, \mathrm{yr}^{-1}$, again at $2\sigma$. These bounds are significantly stronger than those from previous primordial nucleosynthesis studies, and are comparable and complementary to CMB, stellar, solar system, lunar laser ranging, pulsar timing and gravitational wave constraints.

New light states thermally coupled to the Standard Model plasma alter the expansion history of the Universe and impact the synthesis of the primordial light elements. In this work, we carry out an exhaustive and precise analysis of the implications of MeV-scale BSM particles in Big Bang Nucleosynthesis (BBN) and for Cosmic Microwave Background (CMB) observations. We find that, BBN observations set a lower bound on the thermal dark matter mass of $m_\chi > 0.4\,\text{MeV}$ at $2\sigma$. This bound is independent of the spin and number of internal degrees of freedom of the particle, of the annihilation being s-wave or p-wave, and of the annihilation final state. Furthermore, we show that current BBN plus CMB observations constrain purely electrophilic and neutrinophilic BSM species to have a mass, $m_\chi > 3.7\,\text{MeV}$ at $2\sigma$. We explore the reach of future BBN measurements and show that upcoming CMB missions should improve the bounds on light BSM thermal states to $m_\chi > (10-15)\,\text{MeV}$. Finally, we demonstrate that very light BSM species thermally coupled to the SM plasma are highly disfavoured by current cosmological observations.