We study the cosmological evolution and phenomenological properties of scalar bosons in the keV to MeV range that have a tiny mixing with the Standard Model Higgs boson. The mixing determines both the abundance of light scalars produced via the freeze-in mechanism and their lifetime. Intriguingly, the parameters required for such scalars to account for all of the dark matter in the present Universe generically predict lifetimes comparable to the sensitivity of present and future indirect detection experiments. In order to accurately determine the relic abundance of light scalars, we calculate freeze-in yields including effects from finite temperatures and quantum statistics and develop a new approach for solving the Boltzmann equation for number-changing processes in the dark sector. We find that light scalars can potentially explain the anomalous x-ray emission at 3.5 keV, while evading constraints from structure formation and predicting potentially observable self-interaction cross sections.

We devise a new user-friendly tool interfaced with the Boltzmann code CLASS to deal with any kind of exotic electromagnetic energy injection in the universe and its impact on anisotropies of the Cosmic Microwave Background. It makes use of the results from standard electromagnetic cascade calculations develop in the context of WIMP annihilation, generalized to incorporate any injection history. We first validate it on a specific WIMP scenario, the Higgs Portal model, confirming that the standard effective on-the-spot treatment is accurate enough. We then analyze the more involved example of evaporating Primordial Black Holes (PBHs) with masses in the range $[3\times10^{13},5\times10^{16}]$g, for which the standard approximations break down. We derive robust CMB bounds on the relic density of evaporating PBHs, ruling out the possibility for PBHs with a monochromatic distribution of masses in the range $[3\times10^{13},2.5\times10^{16}]$g to represent all of the Dark Matter in our Universe. Remarkably, we confirm with an accurate study that the CMB bounds are several orders of magnitude stronger than those from the galactic gamma-ray background in the range $[3\times10^{13},3\times10^{14}]$g. A future CMB experiment like CORE+, or an experiment attempting at measuring the 21 cm signal from the Dark Ages could greatly improve the sensitivity to these models.