It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

The Universe's early thermal history is poorly constrained, and it is possible that it underwent a period of early matter domination driven by a heavy particle or an oscillating scalar field that decayed into radiation before the onset of Big Bang nucleosynthesis. The entropy sourced by this particle's decay reduces the cross section required for thermal-relic dark matter to achieve the observed abundance. This degeneracy between dark matter properties and the thermal history vastly widens the field of viable dark matter candidates, undermining efforts to constrain dark matter's identity. Fortunately, an early matter-dominated era also amplifies density fluctuations at small scales and leads to early microhalo formation, boosting the dark matter annihilation rate and bringing smaller cross sections into the view of existing indirect-detection probes. Employing several recently developed models of microhalo formation and evolution, we develop a procedure to derive indirect-detection constraints on dark matter annihilation in cosmologies with early matter domination. This procedure properly accounts for the unique morphology of microhalo-dominated signals. While constraints depend on dark matter's free-streaming scale, the microhalos make it possible to obtain upper bounds as small as $\langle\sigma v\rangle \lesssim 10^{-32}$ cm$^3$s$^{-1}$ using Fermi-LAT observations of the isotropic gamma-ray background and the Draco dwarf galaxy.

In the cold dark matter scenario, the smallest dark matter halos may be earth mass or smaller. These microhalos would be the densest dark matter objects in the Universe, making their accurate characterization important for astrophysical dark matter detection efforts. Moreover, their properties are closely linked to the nature of dark matter and the physics of the early universe, making them valuable cosmological probes. Dark matter microhalos survive as subhalos within larger galactic halos, but due to their small size, they are susceptible to encounters with individual stars. We use a large number of $N$-body simulations to develop a framework that can predict the evolution of a microhalo's density profile due to stellar encounters. We find that there is a universal density profile for microhalos subjected to stellar encounters, and each encounter alters a microhalo's scale parameters in a way that is simply related to the energy the encounter injects. Our framework can rapidly and accurately characterize the impact of stellar encounters on whole ensembles of microhalos, making it a promising tool for understanding the populations of microhalos within galactic halos.

Dark matter halos grow by hierarchical clustering as they merge together to produce ever larger structures. During these merger processes, the smaller halo can potentially survive as a subhalo of the larger halo, so a galaxy-scale halo today likely possesses a rich abundance of substructure. This substructure can greatly boost the rate of dark matter annihilation within the host halo, but the precise magnitude of this boost is clouded by uncertainty about the survival prospects of these subhalos. In particular, tidal forces gradually strip material from the subhalos, reducing their annihilation signals and potentially destroying them. In this work, we use high-resolution idealized $N$-body simulations to develop and tune a model that can predict the impact of this tidal evolution on the annihilation rates within subhalos. This model predicts the time evolution of a subhalo's annihilation rate as a function of three physically motivated parameters of the host-subhalo system: the energy injected into subhalo particles per orbit about the host, the ratio of stretching to compressive tidal forces, and the radial distribution of tidal heating within the subhalo. Our model will improve the accuracy of predictions of the magnitude and morphology of annihilation signals from dark matter substructure. Additionally, our parametrization can describe the time evolution of other subhalo properties, so it has implications for understanding aspects of subhalo tidal evolution beyond the annihilation rate.

The first dark matter halos form by direct collapse from peaks in the matter density field, and evidence from numerical simulations and other analyses suggests that the dense inner regions of these objects largely persist today. These halos would be the densest dark matter structures in the Universe, and their abundance can probe processes that leave imprints on the primordial density field, such as inflation or an early matter-dominated era. They can also probe dark matter through its free-streaming scale. The first halos are qualitatively different from halos that form by hierarchical clustering, as evidenced by their $\rho\propto r^{-3/2}$ inner density profiles. In this work, we present and tune models that predict the density profiles of these halos from properties of the density peaks from which they collapsed. These models predict the coefficient $A$ of the $\rho=Ar^{-3/2}$ small-radius asymptote of the density profile along with the maximum circular velocity $v_\mathrm{max}$ and associated radius $r_\mathrm{max}$. These models are universal; they can be applied to any cosmology, and we confirm this by validating them using six $N$-body simulations carried out in wildly disparate cosmological scenarios. We find that these models can even predict the full population of halos with reasonable accuracy in scenarios with narrowly supported power spectra, although for broader power spectra, an understanding of the impact of halo mergers is needed. With their connection to the primordial density field established, the first dark matter halos will serve as probes of the early Universe and the nature of dark matter.

If the dark matter is part of a hidden sector with only very feeble couplings to the Standard Model, the lightest particle in the hidden sector will generically be long-lived and could come to dominate the energy density of the universe prior to the onset of nucleosynthesis. During this early matter-dominated era, density perturbations will grow more quickly than otherwise predicted, leading to a large abundance of sub-earth-mass dark matter microhalos. Since the dark matter does not couple directly to the Standard Model, the minimum halo mass is much smaller than expected for weakly interacting dark matter, and the smallest halos could form during the radiation-dominated era. In this paper, we calculate the evolution of density perturbations within the context of such hidden sector models and use a series of $N$-body simulations to determine the outcome of nonlinear collapse during radiation domination. The resulting microhalos are extremely dense, which leads to very high rates of dark matter annihilation and to large indirect detection signals that resemble those ordinarily predicted for decaying dark matter. We find that the Fermi Collaboration's measurement of the high-latitude gamma-ray background rules out a wide range of parameter space within this class of models. The scenarios that are most difficult to constrain are those that feature a very long early matter-dominated era; if microhalos form prior to the decay of the unstable hidden sector matter, the destruction of these microhalos effectively heats the dark matter, suppressing the later formation of microhalos.

Enhanced density fluctuations on small scales would lead to the formation of numerous dark matter minihalos, so limits on the minihalo abundance can place upper bounds on the small-scale primordial power spectrum. In particular, the ultracompact minihalo (UCMH), a dark matter structure hypothesized to possess a $\rho\propto r^{-9/4}$ density profile due to its formation at $z\geq 1000$, has been used to establish an upper bound on the primordial power spectrum at scales smaller than 2 Mpc. The extreme slope of this density profile amplifies the observational signals of UCMHs. However, we recently showed via N-body simulations that the $\rho\propto r^{-9/4}$ density profile does not develop in realistic formation scenarios, throwing UCMH-derived power spectrum constraints into question. Instead, minihalos develop shallower inner profiles with power-law indices between $-3/2$ and $-1$. In this paper, we expand on that result and discuss its implications. Using a model that is calibrated to simulation results and predicts halo structures in spiked power spectra based on their formation times, we calculate new upper bounds on the primordial power spectrum based on limits on the dark matter annihilation rate within the Galaxy. We find that despite assuming shallower profiles, this minihalo model actually yields stronger constraints than the previous UCMH picture owing to its inclusion of all minihalos instead of only the earliest-forming ones.

Ultracompact minihalos (UCMHs) have emerged as a valuable probe of the primordial power spectrum of density fluctuations at small scales. UCMHs are expected to form at early times in regions with ${\delta\rho/\rho \gtrsim 10^{-3}}$, and they are theorized to possess an extremely compact ${\rho\propto r^{-9/4}}$ radial density profile, which enhances their observable signatures. Nonobservation of UCMHs can thus constrain the primordial power spectrum. Using $N$-body simulations to study the collapse of extreme density peaks at ${z \simeq 1000}$, we show that UCMHs forming under realistic conditions do not develop the ${\rho\propto r^{-9/4}}$ profile and instead develop either ${\rho\propto r^{-3/2}}$ or ${\rho\propto r^{-1}}$ inner density profiles depending on the shape of the power spectrum. We also demonstrate via idealized simulations that self-similarity---the absence of a scale length---is necessary to produce a halo with the ${\rho\propto r^{-9/4}}$ profile, and we argue that this implies such halos cannot form from a Gaussian primordial density field. Prior constraints derived from UCMH nonobservation must be reworked in light of this discovery. Although the shallower density profile reduces UCMH visibility, our findings reduce their signal by as little as $\mathcal O(10^{-2})$ while allowing later-forming halos to be considered, which suggests that new constraints could be significantly stronger.