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.

It is quite conceivable that dark matter freeze-out occurred during an early period of matter domination, in which case the evolution and relic abundance differ from standard freeze-out calculations which assume a radiation dominated universe. Here we re-examine the classic models in which dark matter interactions with the Standard Model are mediated via either the Higgs or $Z$ boson in the context of matter dominated freeze-out. We highlight that while these classic models are largely excluded by searches in the radiation dominated case, matter dominated freeze-out can relax these limits and thus revive the Higgs and $Z$ portals. Additionally, we discuss the distinctions between matter dominated freeze-out and decoupling during the transition from matter domination to radiation domination, and we comment on the parameter regimes which lead to non-negligible dark matter production during this transition.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

We highlight that the anomalous orbits of Trans-Neptunian Objects (TNOs) and an excess in microlensing events in the 5-year OGLE dataset can be simultaneously explained by a new population of astrophysical bodies with mass several times that of Earth ($M_\oplus$). We take these objects to be primordial black holes (PBHs) and point out the orbits of TNOs would be altered if one of these PBHs was captured by the Solar System, inline with the Planet 9 hypothesis. Capture of a free floating planet is a leading explanation for the origin of Planet 9 and we show that the probability of capturing a PBH instead is comparable. The observational constraints on a PBH in the outer Solar System significantly differ from the case of a new ninth planet. This scenario could be confirmed through annihilation signals from the dark matter microhalo around the PBH.

A notable feature of UV freeze-in is that the relic density is strongly dependent on the highest temperatures of the thermal bath, and a common assumption is that the relevant 'highest temperature' should be the reheating temperature after inflation $T_\text{RH}$. However, the temperature of the thermal bath can be significantly higher in certain scenarios, reaching a value denoted T max , a fact which is only apparent away from the instantaneous decay approximation. Interestingly, it has been shown that if the operators are of sufficiently high mass dimension then the dark matter abundance can be enhanced by a 'boost factor' depending on ($T_\text{max}/T_\text{RH}$) relative to naive estimates assuming instantaneous reheating. We highlight here that in non-standard cosmological histories the critical mass dimension of the operator above at which the instantaneous decay approximation breaks down, and the exponent of the boost factor, depend on the equation of state $\omega$ prior to reheating. We highlight four examples in which the dark matter abundance receives a significant enhancement in the context of gravitino dark matter, the moduli portal, the Higgs portal, and the spin-2 portal (as might arise in bimetric gravity models). We comment on the transition from kination domination to radiation domination as a motivated example of non-standard cosmologies.

If the early universe is dominated by an energy density which evolves other than radiation-like the normal Hubble-temperature relation $H\propto T^2$ is broken and dark matter relic density calculations in this era can be significantly different. We first highlight that for a population of states $\phi$ sourcing an initial expansion rate of the form $H\propto T^{2+n/2}$ for $n\geq-4$, during the period of appreciable $\phi$ decays the evolution transitions to $H\propto T^4$. The decays of $\phi$ imply a source of entropy production in the thermal bath which alters the Boltzmann equations and impacts the dark matter relic abundance. We show that the form of the initial expansion rate leaves a lasting imprint on relic densities established while $H\propto T^4$ since the value of the exponent $n$ changes the temperature evolution of the thermal bath. In particular, a dark matter relic density set via freeze-in or non-thermal production is highly sensitive to the temperature dependance of the initial expansion rate. This work generalises earlier studies which assumed initial expansion rates due to matter or kination domination.

We highlight the general scenario of dark matter freeze-out whilst the energy density of the universe is dominated by a decoupled non-relativistic species. Decoupling during matter domination changes the freeze-out dynamics, since the Hubble rate is parametrically different for matter and radiation domination. Furthermore, for successful Big Bang Nucleosynthesis the state dominating the early universe energy density must decay, this dilutes (or repopulates) the dark matter. As a result, the masses and couplings required to reproduce the observed dark matter relic density can differ significantly from radiation dominated freeze-out.

The early universe could feature multiple reheating events, leading to jumps in the visible sector entropy density that dilute both particle asymmetries and the number density of frozen-out states. In fact, late time entropy jumps are usually required in models of Affleck-Dine baryogenesis, which typically produces an initial particle-antiparticle asymmetry that is much too large. An important consequence of late time dilution, is that a smaller dark matter annihilation cross section is needed to obtain the observed dark matter relic density. For cosmologies with high scale baryogenesis, followed by radiation-dominated dark matter freeze-out, we show that the perturbative unitarity mass bound on thermal relic dark matter is relaxed to $10^{10}$ GeV. We proceed to study superheavy asymmetric dark matter models, made possible by a sizable entropy injection after dark matter freeze-out, and identify how the Affleck-Dine mechanism would generate the baryon and dark asymmetries.

We show that Fermi repulsion can lead to cored density profiles in dwarf galaxies for sub-keV fermionic dark matter. We treat the dark matter as a quasi-degenerate self-gravitating Fermi gas and calculate its density profile assuming hydrostatic equilibrium. We find that suitable dwarf galaxy cores of larger than 130 pc can be achieved for fermion dark matter with mass in the range 70 eV - 400 eV. While in conventional dark matter scenarios, such sub-keV thermal dark matter would be excluded by free streaming bounds, the constraints are ameliorated in models with dark matter at lower temperature than conventional thermal scenarios, such as the Flooded Dark Matter model that we have previously considered. Modifying the arguments of Tremaine and Gunn we derive a conservative lower bound on the mass of fermionic dark matter of 70 eV and a stronger lower bound from Lyman-$\alpha$ clouds of about 470 eV, leading to slightly smaller cores than have been observed. We comment on this result and how the tension is relaxed in dark matter scenarios with non-thermal momentum distributions.

Most dark matter (DM) models set the DM relic density by some interaction with Standard Model particles. Such models generally assume the existence of Standard Model particles early on, with the DM relic density a later consequence of those interactions. Perhaps a more compelling assumption is that DM is not part of the Standard Model sector and a population of DM too is generated at the end of inflation. This democratic assumption does not necessarily provide a natural value for the DM relic density, and superficially leads to too much entropy in the dark sector. We address the latter issue by the late decay of heavy particles produced at early times, associating the DM relic density with the lifetime of a long-lived state. We ask what it would take for this scenario to be compatible with observations in what we call Flooded Dark Matter (FDM) and discuss several interesting consequences. One is that DM can be very light and furthermore, light DM is in some sense the most natural scenario in FDM as it is compatible with larger couplings of the decaying particle. Moreover, the decay of the field with the smallest coupling and hence the longest lifetime dominates the entropy and possibly the matter content of the Universe, a principle we refer to as 'Maximum Baroqueness'. We also show that the dark sector should be colder than the ordinary sector, relaxing the free-streaming constraints on light DM. We will discuss the implications for the core-cusp problem in a follow-up paper. FDM also has interesting baryogenesis implications. One possibility is that both DM and baryon asymmetries are simultaneously diluted by a late entropy dump. Alternatively, FDM is compatible with an elegant non-thermal leptogenesis implementation in which decays of a heavy RH neutrino lead to late time reheating of the Standard Model and provide suitable conditions for creation of a lepton asymmetry.

If dark matter is thermally decoupled from the visible sector, the observed relic density can potentially be obtained via freeze-in production of dark matter. Typically in such models it is assumed that the dark matter is connected to the thermal bath through feeble renormalisable interactions. Here, rather, we consider the case in which the hidden and visible sectors are coupled only via non-renormalisable operators. This is arguably a more generic realisation of the dark matter freeze-in scenario, as it does not require the introduction of diminutive renormalisable couplings. We examine general aspects of freeze-in via non-renormalisable operators in a number of toy models and present several motivated implementations in the context of Beyond the Standard Model physics. Specifically, we study models related to the Peccei-Quinn mechanism and Z' portals.

We derive the particle asymmetry due to inflationary baryogenesis involving a complex inflaton, obtaining a different result to that in the literature. While asymmetries are found to be significantly smaller than previously calculated, in certain parameter regions baryogenesis can still be achieved.

In models of freeze-in the dark matter (DM) is decoupled from the visible sector and initially has a depleted number density. The hidden and visible sectors are connected only via a feeble portal interaction by which DM can be produced. Asymmetric freeze-in (AFI) combines this scenario with ideas from asymmetric DM and provides a potential cogenesis mechanism. However, it has been argued that existing AFI models do not produce suitably large particle asymmetries due to cancellations which arises because the mediator state remains in thermal equilibrium. We examine AFI via an out-of-equilibrium mediator and using a simple scalar model show that in this case sizeable asymmetries may be generated.

The apparent excess of gamma rays in an extended region in the direction of the galactic center has a spatial distribution and amplitude that are suggestive of dark matter annihilations. If this excess is indeed due to dark matter annihilations, it would indicate the presence of both dark matter and an additional particle beyond the Standard Model that mediates the interactions between the dark matter and Standard Model states. We introduce reference models describing dark matter annihilation to pairs of these new mediators, which decouples the SM-mediator coupling from the thermal annihilation cross section and easily explains the lack of direct detection signals. We determine the parameter regions that give good descriptions of the gamma ray excess for several motivated choices of mediator couplings to the SM. We find fermion dark matter with mass 7-26 GeV and a dark vector mediator, or scalar dark matter in the 10-50 GeV range (Higgs portal mediator) or 10-65 GeV range (gluophilic mediator) can provide a comparable or improved fit, compared to the case of direct annihilation. We demonstrate that these models can easily satisfy all constraints from collider experiments, direct detection, and cosmology.

In the simplest models of asymmetric dark matter (ADM) annihilation signals are not expected, since the DM is non-self-conjugate and the relic density of anti-DM is negligible. We investigate a new class of models in which a symmetric DM component, in the `low-mass' 1-10 GeV regime favoured for linking the DM and baryon asymmetries, is repopulated through decays. We find that, in models without significant velocity dependence of the annihilation cross section, observational constraints generally force these decays to be (cosmologically) slow. These late decays can give rise to gamma-ray signal morphologies differing from usual annihilation profiles. A distinctive feature of such models is that signals may be absent from dwarf spheroidal galaxies.

If R-parity is only mildly violated then the lightest supersymmetric particle (LSP) can be stable over cosmologically time-scales and still account for the dark matter relic density. We examine the possibility of generating detectable X-ray lines from R-parity violating decays of keV-scale LSP dark matter to neutrino-photon pairs. Specifically, we consider scenarios in which the LSP is a light gravitino, bino, or hidden sector photino. Potential signals are discussed in the context of recent claims of an unidentified 3.5 keV X-ray line in studies of stacked galaxy clusters. We comment on the difficulties in obtaining the observed relic density for keV scale bino or hidden photino dark matter and some possible resolutions.

We propose a new framework for explaining the proximity of the baryon and dark matter relic densities \Omega_{DM} \approx 5\Omega_B. The scenario assumes that the number density of the observed dark matter states is generated due to decays from a second hidden sector which simultaneously generates the baryon asymmetry. In contrast to asymmetric dark matter models, the dark matter can be a real scalar or Majorana fermion and thus presents distinct phenomenology. We discuss aspects of model building and general constraints in this framework. We present a simple supersymmetric implementation of this mechanism and show that it can be used to obtain the correct dark matter relic density for a bino LSP.

It is argued that experimental constraints on theories of asymmetric dark matter (ADM) almost certainly require that the DM be part of a richer hidden sector of interacting states of comparable mass or lighter. A general requisite of models of ADM is that the vast majority of the symmetric component of the DM number density must be removed in order to explain the observed relationship $\Omega_B\sim\Omega_{DM}$ via the DM asymmetry. Demanding the efficient annihilation of the symmetric component leads to a tension with experimental limits if the annihilation is directly to Standard Model (SM) degrees of freedom. A comprehensive effective operator analysis of the model independent constraints on ADM from direct detection experiments and LHC monojet searches is presented. Notably, the limits obtained essentially exclude models of ADM with mass 1GeV$\lesssim m_{DM} \lesssim$ 100GeV annihilating to SM quarks via heavy mediator states. This motivates the study of portal interactions between the dark and SM sectors mediated by light states. Resonances and threshold effects involving the new light states are shown to be important for determining the exclusion limits.