Many models of high energy physics suggest that the cosmological dark sector consists of not just one, but a spectrum of ultralight scalar particles with logarithmically distributed masses. To study the potential signatures of low concentrations of ultralight axion (also known as fuzzy) dark matter, we modify Lagrangian perturbation theory (LPT) by distinguishing between trajectories of different dark matter species. We further adapt LPT to include the effects of a quantum pressure, which is necessary to generate correct initial conditions for ultralight axion simulations. Based on LPT, our modified scheme is extremely efficient on large scales and it can be extended to an arbitrary number of particle species at very little computational cost. This allows for computation of self-consistent initial conditions in mixed dark matter models. Additionally, we find that shell-crossing is delayed for ultralight particles and that the deformation tensor extracted from LPT can be used to identify the range of redshifts and scales for which the Madelung formalism of fuzzy dark matter is a reliable approximation.

Intensity mapping (IM) of spectral lines has the potential to revolutionize cosmology by increasing the total number of observed modes by several orders of magnitude compared to the cosmic microwave background (CMB) anisotropies. In this paper, we consider IM of neutral hydrogen (HI) in the redshift range $0 \lesssim z \lesssim 3$ employing a halo model approach where HI is assumed to follow the distribution of dark matter (DM) halos. If a portion of the DM is composed of ultralight axions then the abundance of halos is changed compared to cold dark matter below the axion Jeans mass. With fixed total HI density, $\Omega_{\rm HI}$, assumed to reside entirely in halos, this effect introduces a scale-independent increase in the HI power spectrum on scales above the axion Jeans scale, which our model predicts consistent with N-body simulations. Lighter axions introduce a scale-dependent feature even on linear scales due to its suppression of the matter power spectrum near the Jeans scale. We use the Fisher matrix formalism to forecast the ability of future HI surveys to constrain the axion fraction of DM and marginalize over astrophysical and model uncertainties. We find that a HIRAX-like survey is a very reliable IM survey configuration, being affected minimally by uncertainties due to non-linear scales, while the SKA1MID configuration is the most constraining as it is sensitive to non-linear scales. Including non-linear scales and combining a SKA1MID-like IM survey with the Simons Observatory CMB, the benchmark ''fuzzy DM'' model with $m_a = 10^{-22}\text{ eV}$ can be constrained at the 10% level. For lighter ULAs this limit improves below 1%, and allows the possibility to test the connection between axion models and the grand unification scale across a wide range of masses.

The light relic density affects the thermal and expansion history of the early Universe leaving a number of observable imprints. We focus on the primordial abundances of light elements produced during the process of Big Bang nucleosynthesis which are influenced by the light relic density. Primordial abundances can be used to infer the density of light relics and thereby serve as a probe of physics beyond the Standard Model. We calculate the observational uncertainty on primordial light element abundances and associated quantities that would be required in order for these measurements to achieve sensitivity to the light relic density comparable to that anticipated from upcoming cosmic microwave background surveys. We identify the nuclear reaction rates that need to be better measured to maximize the utility of future observations. We show that improved measurements of the primordial helium-4 abundance can improve constraints on light relics, while more precise measurements of the primordial deuterium abundance are unlikely to be competitive with cosmic microwave background measurements of the light relic density.

Cosmic string loops contain cusps which decay by emitting bursts of particles. A significant fraction of the released energy is in the form of photons. These photons are injected non-thermally and can hence cause spectral distortions of the Cosmic Microwave Background (CMB). Under the assumption that cusps are robust against gravitational back-reaction, we compute the fractional energy density released as photons in the redshift interval where such non-thermal photon injection causes CMB spectral distortions. Whereas current constraints on such spectral distortions are not strong enough to constrain the string tension, future missions such as the PIXIE experiment will be able to provide limits which rule out a range of string tensions between $G \mu \sim 10^{-15}$ and $G \mu \sim 10^{-12}$, thus ruling out particle physics models yielding these kind of intermediate-scale cosmic strings.