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.

A new component added to the standard model of cosmology that behaves like a cosmological constant at early times and then dilutes away as radiation or faster can resolve the Hubble tension. We show that a rolling axion coupled to a non-Abelian gauge group exhibits the behavior of such an extra component at the background level and can present a natural particle-physics model solution to the Hubble tension. We compare the contribution of this bottom-up model to the phenomenological fluid approximation and determine that CMB observables sensitive only to the background evolution of the Universe are expected to be similar in both cases, strengthening the case for this model to provide a viable solution to the Hubble tension.

The Hubble tension between the $\Lambda$CDM-model-dependent prediction of the current expansion rate $H_0$ using Planck data and direct, model-independent measurements in the local universe from the SH0ES collaboration disagree at $>3.5\sigma$. Moreover, there exists a milder $\sim 2\sigma$ tension between similar predictions for the amplitude $S_8$ of matter fluctuations and its measurement in the local universe. As explanations relying on unresolved systematics have not been found, theorists have been exploring explanations for these anomalies that modify the cosmological model, altering early-universe-based predictions for these parameters. However, new cosmological models that attempt to resolve one tension often worsen the other. In this paper, we investigate a decaying dark matter (DDM) model as a solution to both tensions simultaneously. Here, a fraction of dark matter density decays into dark radiation. The decay rate $\Gamma$ is proportional to the Hubble rate $H$ through the constant $\alpha_{\rm dr}$, the only additional parameter of this model. Then, this model deviates most from $\Lambda$CDM in the early universe, with $\alpha_{\rm dr}$ being positively correlated with $H_0$ and negatively with $S_8$. Hence, increasing $\alpha_{\rm dr}$ (and allowing dark matter to decay in this way) can then diminish both tensions simultaneously. When only considering Planck CMB data and the local SH0ES prior on $H_0$, $\sim 1$\% dark matter decays, decreasing the $S_8$ tension to $0.3\sigma$ and increasing the best-fit $H_0$ by $1.6$ km/s/Mpc. However, the addition of intermediate-redshift data (the JLA supernova dataset and baryon acoustic oscillation data) weakens the effectiveness of this model. Only $\sim 0.5$\% of the dark matter decays bringing the $S_8$ tension back up to $\sim 1.5 \sigma$ and the increase in the best-fit $H_0$ down to $0.4$ km/s/Mpc.

Early dark energy (EDE) that behaves like a cosmological constant at early times (redshifts $z\gtrsim3000$) and then dilutes away like radiation or faster at later times can solve the Hubble tension. In these models, the sound horizon at decoupling is reduced resulting in a larger value of the Hubble parameter $H_0$ inferred from the cosmic microwave background (CMB). We consider two physical models for this EDE, one involving an oscillating scalar field and another a slowly-rolling field. We perform a detailed calculation of the evolution of perturbations in these models. A Markov Chain Monte Carlo search of the parameter space for the EDE parameters, in conjunction with the standard cosmological parameters, identifies regions in which $H_0$ inferred from {\it Planck} CMB data agrees with the SH0ES local measurement. In these cosmologies, current baryon acoustic oscillation and supernova data are described as successfully as in \LCDM, while the fit to {\it Planck} data is slightly improved. Future CMB and large-scale-structure surveys will further probe this scenario.

Pulsars have been invoked to explain the origin of recently observed high-energy Galactic cosmic-ray positrons. Since the positron propagation distance decreases with energy, the number of pulsars that can contribute to the observed positrons decreases from $O(10^3)$ for positron energies $E\gtrsim10$ GeV to only a few for $E \gtrsim 500$ GeV. Thus, if pulsars explain these positrons, the positron energy spectrum should become increasingly bumpy at higher energies. Here we present a power-spectrum analysis that can be applied to seek such spectral features in the energy spectrum for cosmic-ray positrons and for the energy spectrum of the combined electron/positron flux. We account for uncertainties in the pulsar distribution by generating hundreds of simulated spectra from pulsar distributions consistent with current observational constraints. Although the current AMS-02 data do not exhibit evidence for spectral features, we find that such features would be detectable in $\simeq 10\%$ of our simulations, with 20 years of AMS-02 data or three years of DAMPE measurements on the electron-plus-positron flux.

Recent measurements of cosmic-ray electron and positron spectra at energies from a GeV to 5 TeV, as well as radio, X-ray and a wide range of gamma-ray observations of pulsar-wind nebulae, indicate that pulsars are significant sources of high-energy cosmic-ray electrons and positrons. Here we calculate the local cosmic-ray $e^\pm$ energy spectra from pulsars taking into account models for (a) the distribution of the pulsars spin-down properties; (b) the cosmic-ray source spectra; and (c) the physics of cosmic-ray propagation. We then use the measured cosmic-ray fluxes from AMS-02, CALET and DAMPE to constrain the space of pulsar and cosmic-ray-propagation models and in particular, local cosmic-ray diffusion and energy losses, the pulsars' energy-loss time-dependence, and the injected $e^{\pm}$ spectra. We find that the lower estimates for the local $e^{\pm}$ energy losses are inconsistent with the data. We also find that pulsar braking indexes of 2.5 or less for sources with ages more than 10 kyr are strongly disfavored. Moreover, the cosmic-ray data are consistent with a wide range of assumptions on the $e^{\pm}$ injection spectral properties and on the distribution of initial spin-down powers. Above a TeV in energy, we find that pulsars can easily explain the observed change in the $e^{+} + e^{-}$ spectral slope. These conclusions are valid as long as pulsars contribute $\gtrsim10\%$ of the observed cosmic-ray $e^\pm$ at energies $\gtrsim100$ GeV.

Cosmological observations are used to test for imprints of an ultra-light axion-like field (ULA), with a range of potentials $V(\phi)\propto[1-\cos(\phi/f)]^n$ set by the axion-field value $\phi$ and decay constant $f$. Scalar field dynamics dictate that the field is initially frozen and then begins to oscillate around its minimum when the Hubble parameter drops below some critical value. For $n\!=\!1$, once dynamical, the axion energy density dilutes as matter; for $n\!=\!2$ it dilutes as radiation and for $n\!=\!3$ it dilutes faster than radiation. Both the homogeneous evolution of the ULA and the dynamics of its linear perturbations are included, using an effective fluid approximation generalized from the usual $n=1$ case. ULA models are parameterized by the redshift $z_c$ when the field becomes dynamical, the fractional energy density $f_{z_c} \equiv \Omega_a(z_c)/\Omega_{\rm tot}(z_c)$ in the axion field at $z_c$, and the effective sound speed $c_s^2$. Using Planck, BAO and JLA data, constraints on $f_{z_c}$ are obtained. ULAs are degenerate with dark energy for all three potentials if $1+z_c \lesssim 10$. When $3\times10^4 \gtrsim 1+z_c \gtrsim 10 $, $f_{z_c}$ is constrained to be $ \lesssim 0.004 $ for $n=1$ and $f_{z_c} \lesssim 0.02 $ for the other two potentials. The constraints then relax with increasing $z_c$. These results strongly constrain ULAs as a resolution to cosmological tensions, such as discrepant measurements of the Hubble constant, or the EDGES measurement of the global 21 cm signal.

Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess.

Precise measurements of the cosmic microwave background (CMB) power spectrum are in excellent agreement with the predictions of the standard $\Lambda$CDM cosmological model. However, there is some tension between the value of the Hubble parameter $H_0$ inferred from the CMB and that inferred from observations of the Universe at lower redshifts, and the unusually small value of the dark-energy density is a puzzling ingredient of the model. In this paper, we explore a scenario with a new exotic energy density that behaves like a cosmological constant at early times and then decays quickly at some critical redshift $z_c$. An exotic energy density like this is motivated by some string-axiverse-inspired scenarios for dark energy. By increasing the expansion rate at early times, the very precisely determined angular scale of the sound horizon at decoupling can be preserved with a larger Hubble constant. We find, however, that the Planck temperature power spectrum tightly constrains the magnitude of the early dark-energy density and thus any shift in the Hubble constant obtained from the CMB. If the reionization optical depth is required to be smaller than the Planck 2016 $2\sigma$ upper bound $\tau\lesssim 0.0774$, then early dark energy allows a Hubble-parameter shift of at most 1.6 km~s$^{-1}$~Mpc$^{-1}$ (at $z_c\simeq 1585$), too small to fully alleviate the Hubble-parameter tension. Only if $\tau$ is increased by more than $5\sigma$ can the CMB Hubble parameter be brought into agreement with that from local measurements. In the process, we derive strong constraints to the contribution of early dark energy at the time of recombination---it can never exceed $\sim2\%$ of the radiation/matter density for $10 \lesssim z_c \lesssim 10^5$.