With several proposed MeV gamma-ray telescopes on the horizon, it is of paramount importance to perform accurate calculations of gamma-ray spectra expected from sub-GeV dark matter annihilation and decay. We present hazma, a python package for reliably computing these spectra, determining the resulting constraints from existing gamma-ray data, and prospects for upcoming telescopes. For high-level analyses, hazma comes with several built-in dark matter models where the interactions between dark matter and hadrons have been determined in detail using chiral perturbation theory. Additionally, hazma provides tools for computing spectra from individual final states with arbitrary numbers of light leptons and mesons, and for analyzing custom dark matter models. hazma can also produce electron and positron spectra from dark matter annihilation, enabling precise derivation of constraints from the cosmic microwave background.

The evaporation of primordial black holes with a mass in the $1\ {\rm gram}\lesssim M_{\rm PBH}\lesssim$1000 kg range can lead to the production of dark matter particles of almost any mass in the range $0.1\ {\rm MeV}\lesssim m_{\rm DM}\lesssim 10^{18}$ GeV with the right relic density at very early times, $\tau\lesssim 10^{-10}$ s. We calculate, as a function of the primordial black holes mass and initial abundance, the combination of dark matter particle masses and number of effective dark degrees of freedom leading to the right abundance of dark matter today, whether or not evaporation stops around the Planck scale. In addition, since black hole evaporation can also lead to the production of a baryon asymmetry, we calculate where dark matter production and baryogenesis can concurrently happen, under a variety of assumptions: baryogenesis via grand unification boson decay, via leptogenesis, or via asymmetric co-genesis of dark matter and ordinary matter. Finally, we comment on possible ways to test this scenario.