Nuggets---very large stable bound objects arising in the presence of a sufficiently attractive and long-range force and in the absence of a dark Coulomb force---are a smoking gun signature for Asymmetric Dark Matter (ADM). The cosmology of ADM nuggets is both generic and unique: nuggets feature highly exothermic fusion processes, which can impact the shape of the core in galaxies, as well as give rise to rare dark star formation. We find, considering the properties of nuggets in a generic extended nuclear model with both attractive and repulsive forces, that self-interaction constraints place an upper bound on nugget masses at the freeze-out of synthesis in the ballpark of $M_{\rm fo} \lesssim 10^{16}$ GeV. We also show that indirect detection strongly constrains models where the scalar mediator binding the nuggets mixes with the Higgs.

We compute the mass function of bound states of Asymmetric Dark Matter--nuggets--synthesized in the early Universe. We apply our results for the nugget density and binding energy computed from a nuclear model to obtain analytic estimates of the typical nugget size exiting synthesis. We numerically solve the Boltzmann equation for synthesis including two-to-two fusion reactions, estimating the impact of bottlenecks on the mass function exiting synthesis. These results provide the basis for studying the late Universe cosmology of nuggets in a future companion paper.

Motivated by gauge coupling unification and dark matter, we present an extension to the Standard Model where both are achieved by adding an extra new matter multiplet. Such considerations lead to a Grand Unified Theory with very heavy WIMPzilla dark matter, which has mass greater than ~10^7 GeV and must be produced before reheating ends. Naturally, we refer to this scenario as GUTzilla dark matter. Here we present a minimal GUTzilla model, adding a vector-like quark multiplet to the Standard Model. Proton decay constraints require the new multiplet to be both color and electroweak charged, which prompts us to include a new confining SU(3) gauge group that binds the multiplet into a neutral composite dark matter candidate. Current direct detection constraints are evaded due to the large dark matter mass; meanwhile, next-generation direct detection and proton decay experiments will probe much of the parameter space. The relic abundance is strongly dependent on the dynamics of the hidden confining sector, and we show that dark matter production during the epoch of reheating can give the right abundance.