Thermal history of composite dark matter

We study the thermodynamic history of composite dark matter models. We start with classifying the models by means of the symmetries partially protecting the composite dark matter decays and constrain their lifetimes. For each model, we determine the impact of number-changing and number-conserving operators on its thermal history. We also develop the analytic formalism to calculate the asymptotic abundance of stable relics. We show how the relative strength between number- changing and number-conserving interactions together with the dark plasma lifetime affect the thermal fate of the various composite models. Additionally, we show that the final dark relic density of composite particles can be diluted due to an entropy increase stemming from dark plasma decay. Finally, we confront the models with experimental bounds. We find that indirect detection experiments are most promising in testing this large class of models.