The excess of electron recoil events seen by the XENON1T experiment has been interpreted as a potential signal of axion-like particles (ALPs), either produced in the Sun, or constituting part of the dark matter halo of the Milky Way. It has also been explained as a consequence of trace amounts of tritium in the experiment. We consider the evidence for the solar and dark-matter ALP hypotheses from the combination of XENON1T data and multiple astrophysical probes, including horizontal branch stars, red giants, and white dwarfs. We briefly address the influence of ALP decays and supernova cooling. While the different datasets are in clear tension for the case of solar ALPs, all measurements can be simultaneously accommodated for the case of a sub-dominant fraction of dark-matter ALPs. Nevertheless, this solution requires the tuning of several a priori unknown parameters, such that for our choices of priors a Bayesian analysis shows no strong preference for the ALP interpretation of the XENON1T excess over the background hypothesis.

We use a modified version of the Peak-Patch excursion set formalism to compute the mass and size distribution of QCD axion miniclusters from a fully non-Gaussian initial density field obtained from numerical simulations of axion string decay. We find strong agreement with N-Body simulations at significantly lower computational cost. We employ a spherical collapse model, and provide fitting functions for the modified barrier in the radiation era. The halo mass function at $z=99$ has a power-law distribution $M^{-0.68}$ for small masses and $M^{-0.35}$ for large masses within the range $10^{-15}\lesssim M\lesssim 10^{-9}M_{\odot}$, with all masses scaling as $(m_a/50\mu\mathrm{eV})^{-0.5}$. We construct merger trees to estimate the collapse redshift and concentration mass relation, $C(M)$, which is well described using analytical results from the initial power spectrum and linear growth. Using the calibrated analytic results to extrapolate to $z=0$, our method predicts a mean concentration $C\sim \mathcal{O}(\text{few})\times10^4$. The low computational cost of our method makes future investigation of the statistics of rare, dense miniclusters easy to achieve.

We present a same-level comparison of the most prominent inversion methods for the reconstruction of the matter density field in the quasi-linear regime from the Ly$\alpha$ forest flux. Moreover, we present a pathway for refining the reconstruction in the framework of numerical optimization. We apply this approach to construct a novel hybrid method. The methods which are used so far for matter reconstructions are the Richardson-Lucy algorithm, an iterative Gauss-Newton method and a statistical approach assuming a one-to-one correspondence between matter and flux. We study these methods for high spectral resolutions such that thermal broadening becomes relevant. The inversion methods are compared on synthetic data (generated with the lognormal approach) with respect to their performance, accuracy, their stability against noise, and their robustness against systematic uncertainties. We conclude that the iterative Gauss-Newton method offers the most accurate reconstruction, in particular at small S/N, but has also the largest numerical complexity and requires the strongest assumptions. The other two algorithms are faster, comparably precise at small noise-levels, and, in the case of the statistical approach, more robust against inaccurate assumptions on the thermal history of the intergalactic medium (IGM). We use these results to refine the statistical approach using regularization. Our new approach has low numerical complexity and makes few assumptions about the history of the IGM, and is shown to be the most accurate reconstruction at small S/N, even if the thermal history of the IGM is not known. Our code will be made publicly available under https://github.com/hmuellergoe/reglyman.

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.

We study the impact of virtual axions on the polarization of photons inside a cavity during the interaction of high-power laser pulses. A novel detection scheme for measuring the axion-induced ellipticity signal during the Light-by-Light (LBL) scattering process is investigated. We show that a momentum exchange between photons in a probe laser beam and a high-intensity target beam may lead to a resonance at the physical mass of the axion. Consequently, the resonant enhancement of vacuum birefringence gives rise to a large ellipticity signal. This signal enhancement can be applied in order to discriminate between the axion contribution to LBL scattering and the standard model contribution due to electron-positron pairs. The sensitivity of the scheme is studied for experimentally feasible probe light sources and ultrahigh intensity laser backgrounds. It is shown that this technique has the potential to probe the QCD axion in the mass range $10^{-2} \textrm{eV} \lesssim m_{a} \lesssim 1 \textrm{eV}$. In this region the axion induced signal surpasses the standard model background.

The Sun provides the energy required to sustain life on Earth and drive our planet's atmospheric circulation. However, establishing a solid physical connection between solar and tropospheric variability has posed a considerable challenge across the spectrum of Earth-system science. The canon of solar variability, the solar fiducial clock, lies almost exclusively with the 400 years of human telescopic observations that demonstrates the waxing and waning number of sunspots, over an 11(ish) year period. Recent research has demonstrated the critical importance of the underlying 22-year magnetic polarity cycle in establishing the shorter sunspot cycle. Integral to the manifestation of the latter is the spatio-temporal overlapping and migration of oppositely polarized magnetic bands. The points when these bands emerge at high solar latitudes and cancel at the equator are separated by almost 20 years. Here we demonstrate the impact of these "termination" points on the Sun's radiative output and particulate shielding of our atmosphere through the dramatically rapid reconfiguration of solar magnetism. These events reset the Sun's fiducial clock and present a new portal to explore the Sun-Earth connection. Using direct observation and proxies of solar activity going back six decades we can, with high statistical significance, demonstrate an apparent correlation between the solar cycle terminations and the largest swings of Earth's oceanic indices---a previously overlooked correspondence. Forecasting the Sun's global behavior places the next solar termination in early 2020; should a major oceanic swing follow, our challenge becomes: when does correlation become causation and how does the process work?

An axion-like particle (ALP) with mass $m_\phi \sim 10^{-15}\text{ eV}$ oscillates with frequency $\sim$1 Hz and it appears naturally as a consequence of grand unification (GUT) in M-theory. However, with a GUT-scale decay constant such an ALP overcloses the Universe, and cannot solve the strong CP problem. In this paper, we present a two axion model in which the 1 Hz ALP constitutes the entirety of the dark matter (DM) while the QCD axion solves the strong CP problem but contributes negligibly to the DM relic density. The mechanism to achieve the correct relic densities relies on low-scale inflation ($m_\phi \lesssim H_{\rm inf}\lesssim 1 \text{ MeV}$) and stochastic drift of both axion fields, which also highly alleviates the cosmological moduli and gravitino problems. We present two variants of the model each compatible with cosmic microwave background anisotropies and big bang nucleosynthesis, one of which in addition satisfies the trans-Planckian censorship conjecture. The 1 Hz ALP has canonical GUT-scale couplings to both photons and nucleons, and may be detectable by ABRACADABRA and CASPEr. Furthermore the model predicts: (1) a value for the strong CP phase which oscillates around $\bar{\theta}_{\rm QCD}\sim 10^{-12}$, within reach of the proton storage ring electric dipole moment experiment, (2) that any intermediate mass black holes observed in gravitational waves from binary inspiral by LISA will have low dimensionless spin and reside on Regge trajectories.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

We propose an experiment for detecting Axion-Like Particles (ALPs) based on the axion-photon interaction in the presence of a non-uniform magnetic field. The impact of virtual ALPs on the polarization of the photons inside a cavity is studied and a detection scheme is proposed. We find that the cavity normal modes are dispersed differently owing to their coupling to the ALPs in the presence of a background magnetic field. This birefringence, in turn, can be observed as a phase difference between the cavity polarization modes. The signal is considerably enhanced close to the resonance frequencies of the cavity and further enhanced for a squeezed light source. We propose to scan the resonances with a variable frequency source. We argue that the amplified signal allows for exclusion of a broad range of axion mass $10^{-3}\text{eV} \lesssim m_{a} \lesssim 1 \text{eV}$ even at very small axion-photon coupling constant with the potential to reach sensitivity to the QCD axion. Our scheme allows for the exclusion of a range of axion masses that has not yet been covered by other experimental techniques.

If the symmetry breaking responsible for axion dark matter production occurs during the radiation-dominated epoch in the early Universe, then this produces large amplitude perturbations that collapse into dense objects known as axion miniclusters. The characteristic minicluster mass, $M_0$, is set by the mass inside the horizon when axion oscillations begin. For the QCD axion $M_0\sim 10^{-10}M_\odot$, however for an axion-like particle $M_0$ can approach $M_\odot$ or higher. Using the Press-Schechter formalism we compute the mass function of halos formed by hierarchical structure formation from these seeds. We compute the concentrations and collapse times of these halos and show that they can grow to be as massive as $10^6M_0$. Within the halos, miniclusters likely remain tightly bound, and we compute their gravitational microlensing signal taking the fraction of axion dark matter collapsed into miniclusters, $f_{\rm MC}$, as a free parameter. A large value of $f_{\rm MC}$ severely weakens constraints on axion scenarios from direct detection experiments. We take into account the non-Gaussian distribution of sizes of miniclusters and determine how this effects the number of microlensing events. We develop the tools to consider microlensing by an extended mass function of non-point-like objects, and use microlensing data to place the first observational constraints on $f_{\rm MC}$. This opens a new window for the potential discovery of the axion.

The fuzzy dark matter (FDM) model treats DM as a bosonic field with astrophysically large de Broglie wavelength. A striking feature of this model is $\mathcal{O}(1)$ fluctuations in the dark matter density on time scales which are shorter than the gravitational timescale. Including for the first time the effect of core oscillations, we demonstrate how such fluctuations lead to heating of star clusters, and thus an increase in their size over time. From the survival of the old star cluster in Eridanus II we infer $m_a\gtrsim 0.6\rightarrow 1\times 10^{-19}\text{ eV}$ within modelling uncertainty if FDM is to compose all of the DM, and derive constraints on the FDM fraction at lower masses. The subhalo mass function in the Milky Way implies $m_a\gtrsim 0.8\times 10^{-21}\text{ eV}$ to successfully form Eridanus II. The window between $10^{-21}\text{ eV}\lesssim m_a\lesssim 10^{-20}\text{ eV}$ is affected by narrow band resonances, and the limited applicability of the diffusion approximation. Some of this window may be consistent with observations of Eridanus II and more detailed investigations are required.

The axion is a hypothetical, well-motivated dark-matter particle whose existence would explain the lack of charge-parity violation in the strong interaction. In addition to this original motivation, an `axiverse' of ultra-light axions (ULAs) with masses $10^{-33}\,{\rm eV}\lesssim m_{\rm a}\lesssim 10^{-10}\,{\rm eV}$ also emerges from string theory. Depending on the mass, such a ULA contributes to the dark-matter density, or alternatively, behaves like dark energy. At these masses, ULAs' classical wave-like properties are astronomically manifested, potentially mitigating observational tensions within the $\Lambda$CDM paradigm on local-group scales. ULAs also provide signatures on small scales such as suppression of structure, interference patterns and solitons to distinguish them from heavier dark matter candidates. Through their gravitational imprint, ULAs in the presently allowed parameter space furnish a host of observational tests to target in the next decade, altering standard predictions for microwave background anisotropies, galaxy clustering, Lyman-$\alpha$ absorption by neutral hydrogen along quasar sightlines, pulsar timing, and the black-hole mass spectrum.

Axions and axion-like particles are compelling candidates for the missing dark matter of the universe. As they undergo gravitational collapse, they can form compact objects such as axion stars or even black holes. In this paper, we study the formation and distribution of such objects. First, we simulate the formation of compact axion stars using numerical relativity with aspherical initial conditions that could represent the final stages of axion dark matter structure formation. We show that the final states of such collapse closely follow the known relationship of initial mass and axion decay constant $f_a$. Second, we demonstrate with a toy model how this information can be used to scan a model density field to predict the number densities and masses of such compact objects. In addition to being detectable by the LIGO/VIRGO gravitational wave interferometer network for axion mass of $10^{-9} < m_a < 10^{-11}$ eV, we show using peak statistics that for $f_a < 0.2M_{pl}$, there exists a "mass gap" between the masses of axion stars and black holes formed from collapse.

Astrophysical observations of spinning BHs, which span $ 5M_\odot\lesssim M_{\rm BH}\lesssim 5\times 10^8 M_\odot$, can be used to exclude the existence of certain massive bosons via the superradiance phenomenon. In this work, we explore for the first time how these measurements can be used to constrain properties of statistical distributions for the masses of multiple bosonic fields. Quite generally, our methodology excludes $N_{\rm ax}\gtrsim 30$ scalar fields with a range of mass distribution widths and central values spanning many orders of magnitude. We demonstrate this for the specific example of axions in string theory and M-theory, where the mass distributions in certain cases take universal forms. We place upper bounds on $N_{\rm ax}$ for certain scenarios of interest realised approximately as mass distributions in M-theory, including the QCD axion, grand unified theories, and fuzzy dark matter.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The QCD axion was originally predicted as a dynamical solution to the strong CP problem. Axion like particles (ALPs) are also a generic prediction of many high energy physics models including string theory. Theoretical models for axions are reviewed, giving a generic multi-axion action with couplings to the standard model. The couplings and masses of these axions can span many orders of magnitude, and cosmology leads us to consider several distinct populations of axions behaving as coherent condensates, or relativistic particles. Light, stable axions are a mainstay dark matter candidate. Axion cosmology and calculation of the relic density are reviewed. A very brief survey is given of the phenomenology of axions arising from their direct couplings to the standard model, and their distinctive gravitational interactions.

If a significant fraction of the dark matter in the Universe is made of an ultra-light scalar field, named fuzzy dark matter (FDM) with a mass $m_a$ of the order of $10^{-22}-10^{-21}$ eV, then its de Broglie wavelength is large enough to impact the physics of large scale structure formation. In particular, the associated cut-off in the linear matter power spectrum modifies the structure of the intergalactic medium (IGM) at the scales probed by the Lyman-$\alpha$ forest of distant quasars. We study this effect by making use of dedicated cosmological simulations which take into account the hydrodynamics of the IGM. We explore heuristically the amplitude of quantum pressure for the FDM masses considered here and conclude that quantum effects should not modify significantly the non-linear evolution of matter density at the scales relevant to the measured Lyman-$\alpha$ flux power, and for $m_a \geq 10^{-22}$ eV. We derive a scaling law between $m_a$ and the mass of the well-studied thermal warm dark matter (WDM) model that is best adapted to the Lyman-$\alpha$ forest data, and differs significantly from the one infered by a simple linear extrapolation. By comparing FDM simulations with the Lyman-$\alpha$ flux power spectra determined from the BOSS survey, and marginalizing over relevant nuisance parameters, we exclude FDM masses in the range $10^{-22} \leq m_a < 2.3\times 10^{-21}$ eV at 95 % CL. Adding higher-resolution Lyman-$\alpha$ spectra extends the exclusion range up to $2.9\times 10^{-21}$ eV. This provides a significant constraint on FDM models tailored to solve the "small-scale problems" of $\Lambda$CDM.

It has been suggested that the internal dynamics of dwarf spheroidal galaxies (dSphs) can be used to test whether or not ultralight axions with $m_a\sim 10^{-22}\text{eV}$ are a preferred dark matter candidate. However, comparisons to theoretical predictions tend to be inconclusive for the simple reason that while most cosmological models consider only dark matter, one observes only baryons. Here we use realistic kinematic mock data catalogs of Milky Way dSph's to show that the "mass-anisotropy degeneracy" in the Jeans equations leads to biased bounds on the axion mass in galaxies with unknown dark matter halo profiles. In galaxies with multiple chemodynamical components this bias can be partly removed by modelling the mass enclosed within each subpopulation. However, analysis of the mock data reveals that the least-biased constraints on the axion mass result from fitting the luminosity-averaged velocity dispersion of the individual chemodynamical components directly. Applying our analysis to two dSph's with reported stellar subcomponents, Fornax and Sculptor, and assuming that the halo profile has not been acted on by baryons, yields core radii $r_{c}>1.5$ kpc and $r_c> 1.2$ kpc respectively, and $m_a<0.4\times 10^{-22}\text{eV}$ at 97.5\% confidence. These bounds are in tension with the number of observed satellites derived from simple (but conservative) estimates of the subhalo mass function in Milky Way-like galaxies. We discuss how baryonic feedback might affect our results, and the impact of such a small axion mass on the growth of structures in the Universe.

We report on a search for ultra-low-mass axion-like dark matter by analysing the ratio of the spin-precession frequencies of stored ultracold neutrons and $^{199}$Hg atoms for an axion-induced oscillating electric dipole moment of the neutron and an axion-wind spin-precession effect. No signal consistent with dark matter is observed for the axion mass range $10^{-24}~\textrm{eV} \le m_a \le 10^{-17}~\textrm{eV}$. Our null result sets the first laboratory constraints on the coupling of axion dark matter to gluons, which improve on astrophysical limits by up to 3 orders of magnitude, and also improves on previous laboratory constraints on the axion coupling to nucleons by up to a factor of 40.

The cosmic microwave background (CMB) places strong constraints on models of dark matter (DM) that deviate from standard cold DM (CDM), and on initial conditions beyond the scalar adiabatic mode. Here, the full \textit{Planck} data set (including temperature, $E$-mode polarisation, and lensing deflection) is used to test the possibility that some fraction of the DM is composed of ultralight axions (ULAs). This represents the first use of CMB lensing to test the ULA model. We find no evidence for a ULA component in the mass range $10^{-33}\leq m_a\leq 10^{-24}\text{ eV}$. We put percent-level constraints on the ULA contribution to the DM, improving by up to a factor of two compared to the case with temperature anisotropies alone. Axion DM also provides a low-energy window onto the high-energy physics of inflation through the interplay between the vacuum misalignment production of axions and isocurvature perturbations. We perform the first systematic investigation into the parameter space of ULA isocurvature, using an accurate isocurvature transfer function at all $m_{a}$ values. We precisely identify a "window of co-existence" for $10^{-25}\text{ eV}\leq m_a\leq10^{-24}\text{ eV}$ where the data allow, simultaneously, a $\sim10\%$ contribution of ULAs to the DM, and $\sim 1\%$ contributions of isocurvature and tensors to the CMB power. ULAs in this window (and \textit{all} lighter ULAs) are shown to be consistent with a large inflationary Hubble parameter, $H_I\sim 10^{14}\text{ GeV}$. The window of co-existence will be fully probed by proposed CMB-S4 observations with increased accuracy in the high-$\ell$ lensing power and low-$\ell$ $E$ and $B$-mode polarisation. If ULAs in the window exist, this could allow for two independent measurements of $H_I$ in the CMB using the axion DM content and isocurvature, and the tensor contribution to $B$-modes.

There are several issues to do with dwarf galaxy predictions in the standard $\Lambda$CDM cosmology that have suscitated much recent debate about the possible modification of the nature of dark matter as providing a solution. We explore a novel solution involving ultra-light axions that can potentially resolve the missing satellites problem, the cusp-core problem, and the `too big to fail' problem. We discuss approximations to non-linear structure formation in dark matter models containing a component of ultra-light axions across four orders of magnitude in mass, $10^{-24}\text{ eV}\lesssim m_a \lesssim 10^{-20}\text{ eV}$, a range too heavy to be well constrained by linear cosmological probes such as the CMB and matter power spectrum, and too light/non-interacting for other astrophysical or terrestrial axion searches. We find that an axion of mass $m_a\approx 10^{-21}\text{ eV}$ contributing approximately $85\%$ of the total dark matter can introduce a significant kpc scale core in a typical Milky Way satellite galaxy in sharp contrast to a thermal relic with a transfer function cut off at the same scale, while still allowing such galaxies to form in significant number. Therefore ultra-light axions do not suffer from the \emph{Catch 22} that applies to using a warm dark matter as a solution to the small scale problems of cold dark matter. Our model simultaneously allows formation of enough high redshift galaxies to allow reconciliation with observational constraints, and also reduces the maximum circular velocities of massive dwarfs so that baryonic feedback may more plausibly resolve the predicted overproduction of massive MWG dwarf satellites.

Axions arise in many theoretical extensions of the Standard Model of particle physics, in particular the "string axiverse". If the axion masses, $m_a$, and (effective) decay constants, $f_a$, lie in specific ranges, then axions contribute to the cosmological dark matter and dark energy densities. We compute the background cosmological (quasi-)observables for models with a large number of axion fields, $n_{\rm ax}\sim \mathcal{O}(10-100)$, with the masses and decay constants drawn from statistical distributions. This reduces the number of parameters from $2n_{\rm ax}$ to a small number of "hyperparameters". We consider a number of distributions, from those motivated purely by statistical considerations, to those where the structure is specified according to a class of M-theory models. Using Bayesian methods we are able to constrain the hyperparameters of the distributions. In some cases the hyperparameters can be related to string theory, e.g. constraining the number ratio of axions to moduli, or the typical decay constant scale needed to provide the correct relic densities. Our methodology incorporates the use of both random matrix theory and Bayesian networks.

The phase transition responsible for axion dark matter production can create large amplitude isocurvature perturbations which collapse into dense objects known as axion miniclusters. We use microlensing data from the EROS survey, and from recent observations with the Subaru Hyper Suprime Cam to place constraints on the minicluster scenario. We compute the microlensing event rate for miniclusters treating them as spatially extended objects with an extended mass function. Using the published bounds on the number of microlensing events we bound the fraction of DM collapsed into miniclusters, $f_{\rm MC}$. For an axion with temperature dependent mass consistent with the QCD axion we find $f_{\rm MC}<0.22(m_a/100\,\mu\text{eV})^{-0.57}$, which represents the first observational constraint on the minicluster fraction. We forecast that a high-efficiency observation of ten nights with Subaru would be sufficient to constrain $f_{\rm MC}\lesssim 0.1$ over the entire QCD axion mass range. We make various approximations to derive these constraints and dedicated analyses by the observing teams of EROS and Subaru are necessary to confirm our results. If accurate theoretical predictions for $f_{\rm MC}$ can be made in future then microlensing can be used to exclude, or discover, the QCD axion. Further details of our computations are presented in a companion paper.

We use state-of-art measurements of the galaxy luminosity function (LF) at z=6, 7 and 8 to derive constraints on warm dark matter (WDM), late-forming dark matter (LFDM) and ultra-light axion dark matter (ULADM) models alternative to the cold dark matter (CDM) paradigm. To this purpose we have run a suite of high-resolution N-body simulations to accurately characterise the low mass-end of the halo mass function and derive DM model predictions of the high-z luminosity function. In order to convert halo masses into UV-magnitudes we introduce an empirical approach based on halo abundance matching which allows us to model the LF in terms of the amplitude and scatter of the ensemble average star formation rate halo mass relation of each DM model, $\langle {\rm SFR}({\rm M_{ h}},z)\rangle$. We find that independent of the DM scenario the average SFR at fixed halo mass increases from z=6 to 8, while the scatter remains constant. At halo mass ${\rm M_{h}}\gtrsim 10^{12}\,{\rm M}_\odot$ h$^{-1}$ the average SFR as function of halo mass follows a double power law trend that is common to all models, while differences occur at smaller masses. In particular, we find that models with a suppressed low-mass halo abundance exhibit higher SFR compared to the CDM results. Using deviance statistics we obtain a lower limit on the WDM thermal relic particle mass, $m_{\rm WDM}\gtrsim 1.5$ keV at $2\sigma$. In the case of LFDM models, the phase transition redshift parameter is bounded to $z_t\gtrsim 8\cdot 10^5$ at $2\sigma$. We find ULADM best-fit models with axion mass $m_a\gtrsim 1.6\cdot 10^{-22}$ eV to be well within $2\sigma$ of the deviance statistics. We remark that measurements at $z=6$ slightly favour a flattening of the LF at faint UV-magnitudes. This tends to prefer some of the non-CDM models in our simulation suite, although not at a statistically significant level to distinguish them from CDM.

The classical equations of motion for an axion with potential $V(\phi)=m_a^2f_a^2 [1-\cos (\phi/f_a)]$ possess quasi-stable, localized, oscillating solutions, which we refer to as "axion stars". We study, for the first time, collapse of axion stars numerically using the full non-linear Einstein equations of general relativity and the full non-perturbative cosine potential. We map regions on an "axion star stability diagram", parameterized by the initial ADM mass, $M_{\rm ADM}$, and axion decay constant, $f_a$. We identify three regions of the parameter space: i) long-lived oscillating axion star solutions, with a base frequency, $m_a$, modulated by self-interactions, ii) collapse to a BH and iii) complete dispersal due to gravitational cooling and interactions. We locate the boundaries of these three regions and an approximate "triple point" $(M_{\rm TP},f_{\rm TP})\sim (2.4 M_{pl}^2/m_a,0.3 M_{pl})$. For $f_a$ below the triple point BH formation proceeds during winding (in the complex $U(1)$ picture) of the axion field near the dispersal phase. This could prevent astrophysical BH formation from axion stars with $f_a\ll M_{pl}$. For larger $f_a\gtrsim f_{\rm TP}$, BH formation occurs through the stable branch and we estimate the mass ratio of the BH to the stable state at the phase boundary to be $\mathcal{O}(1)$ within numerical uncertainty. We discuss the observational relevance of our findings for axion stars as BH seeds, which are supermassive in the case of ultralight axions. For the QCD axion, the typical BH mass formed from axion star collapse is $M_{\rm BH}\sim 3.4 (f_a/0.6 M_{pl})^{1.2} M_\odot$.

The highly populated, low-energy excitations of a scalar field of mass $m_a\sim 10^{-22}\,\textrm{eV}$ can represent the full dark matter content of the universe and alleviate some tensions in the standard cosmological scenario on small-scales. This {\it fuzzy dark matter} component is commonly assumed to arise as the consequence of a new axion-like particle in the matter sector, yet for simplicity it is usually modeled in terms of a simple free quadratic field. In this paper we consider how the cosmological constraints are modified when the effects of an instanton potential and temperature-dependent mass are included. Current isocurvature and tensor bounds confirms that this particle should be formed before the end of a low-scale inflation period with Hubble parameter $H_I\lesssim 2.5\times 10^{12}\,\textrm{GeV}$, in accordance with previous free-field analysis. The axion decay constant, $f_a$, which fixes axion couplings, appears in the instanton potential and determines the relic abundance, the stability of galaxy cores to axion emission, and the direct searches of fuzzy dark matter. If the axion mass is $T$-independent, we find that $f_a\gtrsim 10^{16}\,\text{GeV}$ is required in order to reach the observed relic density without fine tuning the initial conditions, while for a $T$-dependent case this bound can be lowered by an order of magnitude. However, the anharmonicities in the instanton potential, and mainly a $T$-dependent mass, can delay the onset of field oscillations, leading to larger physical suppression scales in the matter power spectrum for a fixed zero-temperature axion mass. This may favor a string or accidental axion over one emerging from a strongly coupled gauge sector if this model is required to provide large galactic halo cores while simultaneously satisfying observational constraints from cosmic structure formation.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

Measurements of cosmic microwave background (CMB) anisotropies provide strong evidence for the existence of dark matter and dark energy. They can also test its composition, probing the energy density and particle mass of different dark-matter and dark-energy components. CMB data have already shown that ultra-light axions (ULAs) with mass in the range $10^{-32}~{\rm eV} \to 10^{-26}~{\rm eV}$ compose a fraction $< 0.01$ of the cosmological critical density. Here, the sensitivity of a proposed CMB-Stage IV (CMB-S4) experiment (assuming a 1 arcmin beam and $< 1~\mu K{\rm-arcmin}$ noise levels over a sky fraction of 0.4) to the density of ULAs and other dark-sector components is assessed. CMB-S4 data should be $\sim 10$ times more sensitive to the ULA energy-density than Planck data alone, across a wide range of ULA masses $10^{-32}< m_{a}< 10^{-23}~{\rm eV}$, and will probe axion decay constants of $f_{a}\approx 10^{16}~{\rm GeV}$, at the grand unified scale. CMB-S4 could improve the CMB lower bound on the ULA mass from $\sim 10^{-25}~{\rm eV}$ to $10^{-23}~{\rm eV}$, nearing the mass range probed by dwarf galaxy abundances and dark-matter halo density profiles. These improvements will allow for a multi-$\sigma$ detection of percent-level departures from CDM over a wide range of masses. Much of this improvement is driven by the effects of weak gravitational lensing on the CMB, which breaks degeneracies between ULAs and neutrinos. We also find that the addition of ULA parameters does not significantly degrade the sensitivity of the CMB to neutrino masses. These results were obtained using the axionCAMB code (a modification to the CAMB Boltzmann code), presented here for public use.

1. Introduction 2. Models: the QCD axion; the strong CP problem; PQWW, KSVZ, DFSZ; anomalies, instantons and the potential; couplings; axions in string theory 3. Production and I.C.'s: SSB and non-perturbative physics; the axion field during inflation and PQ SSB; cosmological populations - decay of parent, topological defects, thermal production, vacuum realignment 4. The Cosmological Field: action; background evolution; misalignment for QCD axion and ALPs; cosmological perturbation theory - i.c.'s, early time treatment, axion sound speed and Jeans scale, transfer functions and WDM; the Schrodinger picture; simualting axions; BEC 5. CMB and LSS: Primary anisotropies; matter power; combined constraints; Isocurvature and inflation 6. Galaxy Formation; halo mass function; high-z and the EOR; density profiles; the CDM small-scale crises 7. Accelerated expansion: the c.c. problem; axion inflation (natural and monodromy) 8. Gravitational interactions with black holes and pulsars 9. Non-gravitational interactions: stellar astrophysics; LSW; vacuum birefringence; axion forces; direct detection with ADMX and CASPEr; Axion decays; dark radiation; astrophysical magnetic fields; cosmological birefringence 10. Conclusions A Theta vacua of gauge theories B EFT for cosmologists C Friedmann equations D Cosmological fluids E Bayes Theorem and priors F Degeneracies and sampling G Sheth-Tormen HMF

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Cold dark matter (CDM) is a well established paradigm to describe cosmological structure formation, and works extraordinarily well on large, linear, scales. Progressing further in dark matter physics requires being able to understand structure formation in the non-linear regime, both for CDM and its alternatives. This short note describes a calculation, and accompanying code, WarmAndFuzzy, incorporating the popular models of warm and fuzzy dark matter (WDM and FDM) into the standard halo model to compute the non-linear matter power spectrum. The FDM halo model power spectrum has not been computed before. The FDM implementation models ultralight axions and other scalar fields with $m_a\approx 10^{-22}\text{ eV}$. The WDM implementation models thermal WDM with mass $m_X\approx 1\text{ keV}$. The halo model shows that differences between WDM, FDM, and CDM survive at low redshifts in the quasi-linear and fully non-linear regimes. The code uses analytic transfer functions for the linear power spectrum, modified collapse barriers in the halo mass function, and a modified concentration-mass relationship for the halo density profiles. Modified halo density profiles (for example, cores) are not included, but are under development. Cores are expected to have very minor effects on the power spectrum on observable scales. Applications of this code to the Lyman-$\alpha$ forest flux power spectrum and the cosmic microwave background lensing power spectrum will be discussed in companion papers. \textsc{WarmAndFuzzy} is available online at \url{https://github.com/DoddyPhysics/HMcode}, where collaboration in development is welcomed.

The particle nature of dark matter remains a mystery. In this paper, we consider two dark matter models---Late Forming Dark Matter (LFDM) and Ultra-Light Axion (ULA) models---where the matter power spectra show novel effects on small scales. The high redshift universe offers a powerful probe of their parameters. In particular, we study two cosmological observables: the neutral hydrogen (HI) redshifted 21-cm signal from the epoch of reionization, and the evolution of the collapsed fraction of HI in the redshift range $2 < z < 5$. We model the theoretical predictions of the models using CDM-like N-body simulations with modified initial conditions, and generate reionization fields using an excursion-set model. The N-body approximation is valid on the length and halo mass scales studied. We show that LFDM and ULA models predict an increase in the HI power spectrum from the epoch of reionization by a factor between 2--10 for a range of scales $0.1<k<4 \, \rm Mpc^{-1}$. Assuming a fiducial model where a neutral hydrogen fraction $\bar{x}_{HI}=0.5$ must be achieved by $z=8$, the reionization process allows us to put approximate bounds on the redshift of dark matter formation $z_f > 4 \times 10^5$ (for LFDM) and the axion mass $m_a > 2.6 \times 10^{-23} \, \rm eV$ (for ULA). The comparison of the collapsed mass fraction inferred from damped Lyman-$\alpha$ observations to the theoretical predictions of our models lead to the weaker bounds: $z_f > 2 \times 10^5$ and $m_a > 10^{-23} \, \rm eV$. These bounds are consistent with other constraints in the literature using different observables; we briefly discuss how these bounds compare with possible constraints from the observation of luminosity function of galaxies at high redshifts. In the case of ULAs, these constraints are also consistent with a solution to the cusp-core problem of CDM.

Using a fundamental discrete symmetry, $\mathbb{Z}_N$, we construct a two-axion model with the QCD axion solving the strong-$CP$ problem, and an ultralight axion (ULA) with $m_{\rm ULA}\approx 10^{-22}\text{ eV}$ providing the dominant form of dark matter (DM). The ULA is light enough to be detectable in cosmology from its imprints on structure formation, and may resolve the small-scale problems of cold DM. The necessary relative DM abundances occur without fine tuning in constructions with decay constants $f_{\rm ULA}\sim 10^{17}\text{ GeV}$, and $f_{\rm QCD}\sim 10^{11}\text{ GeV}$. An example model achieving this has $N=24$, and we construct a range of other possibilities. We compute the ULA couplings to the Standard Model, and discuss prospects for direct detection. The QCD axion may be detectable in standard experiments through the $\vec{E}\cdot\vec{B}$ and $G\tilde{G}$ couplings. In the simplest models, however, the ULA has identically zero coupling to both $G\tilde{G}$ of QCD and $\vec{E}\cdot\vec{B}$ of electromagnetism due to vanishing electromagnetic and color anomalies. The ULA couples to fermions with strength $g\propto 1/f_{\rm ULA}$. This coupling causes spin precession of nucleons and electrons with respect to the DM wind with period $t\sim$months. Current limits do not exclude the predicted coupling strength, and our model is within reach of the CASPEr-Wind experiment, using nuclear magnetic resonance.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Self-gravitating bosonic fields can support stable and localised field configurations. For real fields, these solutions oscillate in time and are known as oscillatons. The density profile is static, and is soliton. Such solitons should be ubiquitous in models of axion dark matter, with the soliton characteristic mass and size depending on some inverse power of the axion mass. Stable configurations of non-relativistic axions are studied numerically using the Schr\"{o}dinger-Poisson system. This method, and the resulting soliton density profiles, are reviewed. Using a scaling symmetry and the uncertainty principle, the core size of the soliton can be related to the central density and axion mass, $m_a$, in a universal way. Solitons have a constant central density due to pressure-support, unlike the cuspy profile of cold dark matter (CDM). One consequence of this fact is that solitons composed of ultra-light axions (ULAs) may resolve the `cusp-core' problem of CDM. In DM halos, thermodynamics will lead to a CDM-like Navarro-Frenk-White profile at large radii, with a central soliton core at small radii. Using Monte-Carlo techniques to explore the possible density profiles of this form, a fit to stellar-kinematical data of dwarf spheroidal galaxies is performed. In order for ULAs to resolve the cusp-core problem (without recourse to baryon feedback or other astrophysical effects) the axion mass must satisfy $m_a<1.1\times 10^{-22}\text{ eV}$ at 95\% C.L. On the other hand, ULAs with $m_a\lesssim 1\times 10^{-22}\text{ eV}$ are in some tension with cosmological structure formation. An axion solution to the cusp-core problem thus makes novel predictions for future measurements of the epoch of reionisation. On the other hand, this can be seen as evidence that structure formation could soon impose a \emph{Catch 22} on axion/scalar field DM, similar to the case of warm DM.

The non-linear hydrodynamic equations for axion/scalar field dark matter (DM) in the non-relativistic Madelung-Shcr\"{o}dinger form are derived in a simple manner, including the effects of universal expansion and Hubble drag. The hydrodynamic equations are used to investigate the relative velocity between axion DM and baryons, and the moving-background perturbation theory (MBPT) derived. Axions massive enough to be all of the DM do not affect the coherence length of the relative velocity, but the MBPT equations are modified by the inclusion of the axion effective sound speed. These MBPT equations are necessary for accurately modelling the effects of axion DM on the formation of the first cosmic structures, and suggest that the 21cm power spectrum could improve constraints on axion mass by up to four orders of magnitude with respect to the current best constraints. A further application of these results uses the "quantum force" analogy to model scalar field gradient energy in a smoothed-particle hydrodynamics model of axion DM. Such a model can treat axion DM in the non-linear regime and could be incorporated into existing N-body codes.

If the dark matter (DM) were composed of axions, then structure formation in the Universe would be suppressed below the axion Jeans scale. Using an analytic model for the halo mass function of a mixed DM model with axions and cold dark matter, combined with the abundance-matching technique, we construct the UV-luminosity function. Axions suppress high-$z$ galaxy formation and the UV-luminosity function is truncated at a faintest limiting magnitude. From the UV-luminosity function, we predict the reionization history of the universe and find that axion DM causes reionization to occur at lower redshift. We search for evidence of axions using the Hubble Ultra Deep Field UV-luminosity function in the redshift range $z=6$-$10$, and the optical depth to reionization, $\tau$, as measured from cosmic microwave background polarization. All probes we consider consistently exclude $m_a\lesssim 10^{-23}\text{ eV}$ from contributing more than half of the DM, with our strongest constraint ruling this model out at more than $8\sigma$ significance. In conservative models of reionization a dominant component of DM with $m_a=10^{-22}\text{ eV}$ is in $3\sigma$ tension with the measured value of $\tau$, putting pressure on an axion solution to the cusp-core problem. Tension is reduced to $2\sigma$ for the axion contributing only half of the DM. A future measurement of the UV-luminosity function in the range $z=10$-$13$ by JWST would provide further evidence for or against $m_a=10^{-22}\text{ eV}$. Probing still higher masses of $m_a=10^{-21}\text{ eV}$ will be possible using future measurements of the kinetic Sunyaev-Zel'dovich effect by Advanced ACTPol to constrain the time and duration of reionization.

Ultra-light axions (ULAs) with masses in the range 10^{-33} eV <m <10^{-20} eV are motivated by string theory and might contribute to either the dark-matter or dark-energy density of the Universe. ULAs could suppress the growth of structure on small scales, or lead to an enhanced integrated Sachs-Wolfe effect on large-scale cosmic microwave-background (CMB) anisotropies. In this work, cosmological observables over the full ULA mass range are computed, and then used to search for evidence of ULAs using CMB data from the Wilkinson Microwave Anisotropy Probe (WMAP), Planck satellite, Atacama Cosmology Telescope, and South Pole Telescope, as well as galaxy clustering data from the WiggleZ galaxy-redshift survey. In the mass range 10^{-32} eV < m <10^{-25.5} eV, the axion relic-density \Omega_{a} (relative to the total dark-matter relic density \Omega_{d}) must obey the constraints \Omega_{a}/\Omega_{d} < 0.05 and \Omega_{a}h^{2} < 0.006 at 95%-confidence. For m> 10^{-24} eV, ULAs are indistinguishable from standard cold dark matter on the length scales probed, and are thus allowed by these data. For m < 10^{-32} eV, ULAs are allowed to compose a significant fraction of the dark energy.

The archetypal theory of dark energy is quintessence: a minimally coupled scalar field with a canonical kinetic energy and potential. By studying random potentials we show that quintessence imposes a restricted set of priors on the equation of state of dark energy. Focusing on the commonly-used parametrisation, $w(a)\approx w_0+w_a(1-a)$, we show that there is a natural scale and direction in the $(w_0, w_a)$ plane that distinguishes quintessence as a general framework. We calculate the expected information gain for a given survey and show that, because of the non-trivial prior information, it is a function of more than just the figure of merit. This allows us to make a quantitative case for novel survey strategies. We show that the scale of the prior sets target observational requirements for gaining significant information. This corresponds to a figure of merit FOM$\gtrsim 200$, a requirement that future galaxy redshift surveys will meet.

A model of high scale inflation is presented where the radial part of the Peccei-Quinn (PQ) field with a non-minimal coupling to gravity plays the role of the inflaton, and the QCD axion is the dark matter. A quantum fluctuation of $\mathcal{O}(H/2\pi)$ in the axion field will result in a smaller angular fluctuation if the PQ field is sitting at a larger radius during inflation than in the vacuum. This changes the effective axion decay constant, $f_a$, during inflation and dramatically reduces the production of isocurvature modes. This mechanism opens up a new window in parameter space where an axion decay constant in the range $10^{12}\text{ GeV}\lesssim f_a\lesssim 10^{15}\text{ GeV}$ is compatible with observably large $r$. The exact range allowed for $f_a$ depends on the efficiency of reheating. This model also predicts a minimum possible value of $r=10^{-3}$. The new window can be explored by a measurement of $r$ possible with \textsc{Spider} and the proposed CASPEr experiment search for high $f_a$ axions.

The recent detection of B-modes by BICEP2 has non-trivial implications for axion dark matter implied by combining the tensor interpretation with isocurvature constraints from Planck. In this paper the measurement is taken as fact, and its implications considered, though further experimental verification is required. In the simplest inflation models $r=0.2$ implies $H_I=1.1\times 10^{14}\text{ GeV}$. If the axion decay constant $f_a<H_I/2\pi$ constraints on the dark matter (DM) abundance alone rule out the QCD axion as DM for $m_a \lesssim 52\chi^{6/7}\,\mu\text{eV}$ (where $\chi>1$ accounts for theoretical uncertainty). If $f_a>H_I/2\pi$ then vacuum fluctuations of the axion field place conflicting demands on axion DM: isocurvature constraints require a DM abundance which is too small to be reached when the back reaction of fluctuations is included. High $f_a$ QCD axions are thus ruled out. Constraints on axion-like particles, as a function of their mass and DM fraction, are also considered. For heavy axions with $m_a\gtrsim 10^{-22}\text{ eV}$ we find $\Omega_a/\Omega_d\lesssim 10^{-3}$, with stronger constraints on heavier axions. Lighter axions, however, are allowed and (inflationary) model-independent constraints from the CMB temperature power spectrum and large scale structure are stronger than those implied by tensor modes.

If low-energy SUSY exists, LHC data favors a high mass scale for scalar superpartners (above a TeV), while sfermions and the dark matter can be parametrically lighter -- leading to a so-called split-spectrum. When combining this fact with the motivation from fundamental theory for shift-symmetric scalars (moduli) prior to SUSY breaking, this leads to a non-thermal history for the early universe. Such a history implies different expectations for the microscopic properties of dark matter, as well as the possibility of dark radiation and a cosmic axion background. In this paper we examine how correlated and mixed isocurvature perturbations are generated in such models, as well as the connection to dark radiation. WMAP constraints on multiple correlated isocurvature modes allow up to half of the primordial perturbations to be isocurvature, contrary to the case of a single isocurvature mode where perturbations must be dominantly adiabatic. However, such bounds are strongly prior dependent, and have not been investigated with the latest Planck data. In this paper we use the example of a SUSY non-thermal history to establish theoretical priors on cosmological parameters. Of particular interest, we find that priors on dark radiation are degenerate with those on the total amount of isocurvature -- they are inversely correlated. Dark radiation is tightly constrained in the early universe and has been used recently to place stringent constraints on string-based approaches to beyond the standard model. Our results suggest such constraints can require more input from theory. Data from cosmological probes -- such as Planck -- can help guide model building in models of the MSSM, split-SUSY, and beyond. Our model also demonstrates the utility of UV models in constructing cosmological priors. (abridged)

Dark sectors, consisting of new, light, weakly-coupled particles that do not interact with the known strong, weak, or electromagnetic forces, are a particularly compelling possibility for new physics. Nature may contain numerous dark sectors, each with their own beautiful structure, distinct particles, and forces. This review summarizes the physics motivation for dark sectors and the exciting opportunities for experimental exploration. It is the summary of the Intensity Frontier subgroup "New, Light, Weakly-coupled Particles" of the Community Summer Study 2013 (Snowmass). We discuss axions, which solve the strong CP problem and are an excellent dark matter candidate, and their generalization to axion-like particles. We also review dark photons and other dark-sector particles, including sub-GeV dark matter, which are theoretically natural, provide for dark matter candidates or new dark matter interactions, and could resolve outstanding puzzles in particle and astro-particle physics. In many cases, the exploration of dark sectors can proceed with existing facilities and comparatively modest experiments. A rich, diverse, and low-cost experimental program has been identified that has the potential for one or more game-changing discoveries. These physics opportunities should be vigorously pursued in the US and elsewhere.

Ultra-light axions ($m_a\lesssim 10^{-18}$eV), motivated by string theory, can be a powerful probe of the energy scale of inflation. In contrast to heavier axions the isocurvature modes in the ultra-light axions can coexist with observable gravitational waves. Here it is shown that large scale structure constraints severely limit the parameter space for axion mass, density fraction and isocurvature amplitude. It is also shown that radically different CMB observables for the ultra-light axion isocurvature mode additionally reduce this space. The results of a new, accurate and efficient method to calculate this isocurvature power spectrum are presented, and can be used to constrain ultra-light axions and inflation.

Measuring spectral distortions of the cosmic microwave background (CMB) is attracting considerable attention as a probe of high energy particle physics in the cosmological context, since PIXIE and PRISM have recently been proposed. In this paper, CMB distortions due to resonant conversions between CMB photons and light axion like particles (ALPs) are investigated, motivated by the string axiverse scenario which suggests the presence of a plenitude of light axion particles. Since these resonant conversions depend on the strength of primordial magnetic fields, constraints on CMB distortions can provide an upper limit on the product of the photon-ALP coupling constant g and the comoving strength of primordial magnetic fields B. Potentially feasible constraints from PIXIE/PRISM can set a limit g B < 10^{-16} GeV^{-1} nG for ALP mass, m_\phi < 10^{-14} eV. Although this result is not a direct constraint on g and B, it is significantly tighter than the product of the current upper limits on g and B.

This paper is concerned with string cosmology and the dynamics of multiple scalar fields in potentials that can become negative, and their features as (Early) Dark Energy models. Our point of departure is the "String Axiverse", a scenario that motivates the existence of cosmologically light axion fields as a generic consequence of string theory. We couple such an axion to its corresponding modulus. We give a detailed presentation of the rich cosmology of such a model, ranging from the setting of initial conditions on the fields during inflation, to the asymptotic future. We present some simplifying assumptions based on the fixing of the axion decay constant $f_a$, and on the effective field theory when the modulus trajectory is adiabatic, and find the conditions under which these assumptions break down. As a by-product of our analysis, we find that relaxing the assumption of fixed $f_a$ leads to the appearance of a new meta-stable de-Sitter region for the modulus without the need for uplifting by an additional constant. A dynamical systems analysis reveals the existence of many fixed point attractors, repellers and saddle points, which we analyse in detail. We also provide geometric interpretations of the phase space. The fixed points can be used to bound the couplings in the model. A systematic scan of certain regions of parameter space reveals that the future evolution of the universe in this model can be rich, containing multiple epochs of accelerated expansion.

A generic prediction of string theory is the existence of many axion fields. It has recently been argued that many of these fields should be light and, like the well known QCD axion, lead to observable cosmological consequences. In this paper we study in detail the effect of the so-called string axiverse on large scale structure, focusing on the morphology and evolution of density perturbations, anisotropies in the cosmic microwave background and weak gravitational lensing of distant galaxies. We quantify specific effects that will arise from the presence of the axionic fields and highlight possible degeneracies that may arise in the presence of massive neutrinos. We take particular care understanding the different physical effects and scales that come into play. We then forecast how the string axiverse may be constrained and show that with a combination of different observations, it should be possible to detect a fraction of ultralight axions to dark matter of a few percent.

A model is presented in the philosophy of the "String Axiverse" of Arvanitaki et al (arXiv:0905.4720v2 [hep-th]) that incorporates a coupling of ultralight axions to their corresponding moduli through the mass term. The light fields roll in their potentials at late times and contribute to the dark sector energy densities in the cosmological expansion. The addition of a coupling and extra field greatly enrich the possible phenomenology of the axiverse. There are a number of interesting phases where the axion and modulus components behave as Dark Matter or Dark Energy and can have considerable and distinct effects on the expansion history of the universe by modifying the equation of state in the past or causing possible future collapse of the universe. In future such a coupling may help to alleviate fine tuning problems for cosmological axions. We motivate and present the model, and briefly explore its cosmological consequences numerically.

Ultra-light scalar fields, with masses of between m=10^{-33} eV and m=10^{-22} eV, can affect the growth of structure in the Universe. We identify the different regimes in the evolution of ultra-light scalar fields, how they affect the expansion rate of the universe and how they affect the growth rate of cosmological perturbations. We find a number of interesting effects, discuss how they might arise in realistic scenarios of the early universe and comment on how they might be observed.