ULYSSES the Third: An Odyssey Towards a Unified Python Toolkit for Leptogenesis

Abstract

We present the third release of $\texttt{ULYSSES}$, a Python package for the numerical evaluation of the baryon asymmetry generated through leptogenesis. This version includes code implementing state-of-the-art density matrix equations for low-scale leptogenesis with three quasi-degenerate right-handed neutrinos. We extend the validity of the code in this scenario beyond the 100 GeV right-handed neutrino mass scale, into the regime of resonant leptogenesis, by including neutrino production rates valid in both the relativistic and non-relativistic regimes. In addition, in the high-scale vanilla scenario, we provide routines for computing $ΔL = 1$ scattering processes, enabling full phase-space evolution of the right-handed neutrino and lepton asymmetry. A new $\texttt{--extended}$ parameter interface allows users to pass model-specific inputs beyond the standard leptogenesis runcard without modifying the core infrastructure and demonstrate its use with a toy module that simultaneously solves the vanilla leptogenesis equations and the freeze-in production of dark matter. On top of these improvements, we introduce an alternative parametrisation of the Casas-Ibarra matrix, update the default neutrino oscillation parameters and report cross-checks of the new low-scale leptogenesis module against published benchmarks and independent codes. $\texttt{ULYSSES}$ is publicly available on $\href{https://github.com/earlyuniverse/ulysses}{\texttt{GitHub}}$ and pip-installable from PyPI.