We describe an algorithm for predicting the concentrations of dark matter halos via a random walk in energy space. Given a full merger tree for a halo, the total internal energy of each halo in that tree is determined by summing the internal and orbital energies of progenitor halos. For halos described by single-parameter density profiles (such as the NFW profile) the energy can be directly mapped to a scale radius, and so to a concentration. We show that this model can accurately reproduce the mean of the concentration mass relation measured in N-body simulations, and reproduces more of the scatter in that relation than previous models. However, our model underpredicts the kurtosis of the distribution of N-body concentrations. We test this model by examining both the autocorrelation of scale radii across time, and the correlations between halo concentration and spin, and comparing to results measured from cosmological N-body simulations. In both cases we find that our model closely matches the N-body results. Our model is implemented within the open source Galacticus toolkit.

Ultra-light axions (ULAs) are a promising dark-matter candidate. ULAs may have implications for small-scale challenges to the $\Lambda$CDM model, and arise in string scenarios. ULAs are already constrained by cosmic microwave background (CMB) experiments and large-scale structure surveys, and will be probed with much greater sensitivity by future efforts. It is challenging to compute observables in ULA scenarios with sufficient speed and accuracy for cosmological data analysis because the ULA field oscillates rapidly. In past work, an effective fluid approximation has been used to make these computations feasible. Here this approximation is tested against an exact solution of the ULA equations, comparing the induced error of CMB observables with the sensitivity of current and future experiments. In the most constrained mass range for a ULA dark matter component ($10^{-27}~{\rm eV}\leq m_{\rm ax}\leq 10^{-25}~{\rm eV}$), the induced bias on the allowed ULA fraction of dark matter from Planck data is less than $1\sigma$. In the cosmic-variance limit (including temperature and polarization data), the bias is $\lesssim 2\sigma$ for primary CMB anisotropies, with more severe biases (as high as $\sim 4\sigma$) resulting for less reliable versions of the effective fluid approximation. If all of the standard cosmological parameters are fixed by other measurements, the expected bias rises to $4-20\sigma$ (well beyond the validity of the Fisher approximation), though the required level of degeneracy breaking will not be achieved by any planned surveys.

We explore the viability of laboratory-scale mechanical resonators as detectors for ultralight scalar dark matter. The signal we investigate is an atomic strain due to modulation of the fine structure constant and the lepton mass at the Compton frequency of dark matter particles. The resulting stress can drive an elastic body with acoustic breathing modes, producing displacements that are accessible with opto- or electromechanical readout techniques. To address the unknown mass of dark matter particles (which determines their Compton frequency), we consider various resonator designs operating at kHz to MHz frequencies, corresponding to $10^{-12}-10^{-5}$ eV particle mass. Current resonant-mass gravitational wave detectors that have been repurposed as dark matter detectors weigh $\sim \! 10^3$ kg. We find that a large unexplored parameter space can be accessed with ultra-high-$Q$, cryogenically-cooled, cm-scale mechanical resonators possessing $\sim \! 10^7$ times smaller mass.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

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.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal).

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The fine-structure constant, $\alpha$, controls the strength of the electromagnetic interaction. There are extensions of the standard model in which $\alpha$ is dynamical on cosmological length and time-scales. The physics of the cosmic microwave background (CMB) depends on the value of $\alpha$. The effects of spatial variation in $\alpha$ on the CMB are similar to those produced by weak lensing: smoothing of the power spectrum, and generation of non-Gaussian features. These would induce a bias to estimates of the weak-lensing potential power spectrum of the CMB. Using this effect, Planck measurements of the temperature and polarization power spectrum, as well as estimates of CMB lensing, are used to place limits (95% C. L.) on the amplitude of a scale-invariant angular power spectrum of $\alpha$ fluctuations relative to the mean value ($C_L^{\alpha}=A^\alpha_{\rm SI}/[L(L+1)]$) of $A^\alpha_{\rm SI}\leq 1.6 \times 10^{-5}$. The limits depend on the assumed shape of the $\alpha$-fluctuation power spectrum. For example, for a white noise angular power spectrum ($C_L^{\alpha}=A^\alpha_{\rm WN}$), the limit is $A^\alpha_{\rm WN}\leq 2.3 \times 10^{-8}$. It is found that the response of the CMB to $\alpha$ fluctuations depends on a separate-universe approximation, such that theoretical predictions are only reliable for $\alpha$ multipoles with $L\lesssim 100$ . An optimal trispectrum estimator can be constructed and it is found that it is only marginally more sensitive than lensing techniques for Planck but significantly more sensitive when considering the next generation of experiments. For a future CMB experiment with cosmic-variance limited polarization sensitivity (e.g., CMB-S4), the optimal estimator could detect $\alpha$ fluctuations with $A^\alpha_{\rm SI}>1.9 \times 10^{-6}$ and $A^\alpha_{\rm WN} > 1.4 \times 10^{-9}$.

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.

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.

Compensated isocurvature perturbations (CIPs) are primordial fluctuations that balance baryon and dark-matter isocurvature to leave the total matter density unperturbed. The effects of CIPs on the cosmic microwave background (CMB) anisotropies are similar to those produced by weak lensing of the CMB: smoothing of the power spectrum, and generation of non-Gaussian features. Previous work considered the CIP effects on the CMB power-spectrum but neglected to include the CIP effects on estimates of the lensing potential power spectrum (though its contribution to the non-Gaussian, connected, part of the CMB trispectrum). Here, the CIP contribution to the standard estimator for the lensing potential power-spectrum is derived, and along with the CIP contributions to the CMB power-spectrum, Planck data is used to place limits on the root-mean-square CIP fluctuations on CMB scales, $\Delta_{\rm rms}^2(R_{\rm CMB})$. The resulting constraint of $\Delta_{\rm rms}^2(R_{\rm CMB}) < 4.3 \times 10^{-3}$ using this new technique improves on past work by a factor of $\sim 3$. We find that for Planck data our constraints almost reach the sensitivity of the optimal CIP estimator. The method presented here is currently the most sensitive probe of the amplitude of a scale-invariant CIP power spectrum placing an upper limit of $A_{\rm CIP}< 0.017$ at 95% CL. Future measurements of the large-scale CMB lensing potential power spectrum could probe CIP amplitudes as low as $\Delta_{\rm rms}^2(R_{\rm CMB}) = 8 \times 10^{-5}$ ($A_{\rm CIP} = 3.2 \times 10^{-4}$).

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.

Inspired by the string axiverse idea, it has been suggested that the recent transition from decelerated to accelerated cosmic expansion is driven by an axion-like quintessence field with a sub-Planckian decay constant. The scenario requires that the axion field be rather near the maximum of its potential, but is less finely tuned than other explanations of cosmic acceleration. The model is parametrized by an axion decay constant $f$, the axion mass $m$, and an initial misalignment angle $|\theta_i|$ which is close to $\pi$. In order to determine the $m$ and $\theta_{i}$ values consistent with observations, these parameters are mapped onto observables: the Hubble parameter $H(z)$ at and angular diameter distance $d_{A}(z)$ to redshift $z= 0.57$, as well as the angular sound horizon of the cosmic microwave background (CMB). Measurements of the baryon acoustic oscillation (BAO) scale at $z\simeq 0.57$ by the BOSS survey and Planck measurements of CMB temperature anisotropies are then used to probe the $\left\{m,f,\theta_i\right\}$ parameter space. With current data, CMB constraints are the most powerful, allowing a fraction of only $\sim 0.2$ of the parameter-space volume. Measurements of the BAO scale made using the SPHEREx or SKA experiments could go further, observationally distinguishing all but $\sim 10^{-2}$ or $\sim 10^{-5}$ of the parameter-space volume (allowed by simple priors) from the $\Lambda$CDM model.

Compensated isocurvature perturbations (CIPs) are modes in which the baryon and dark matter density fluctuations cancel. They arise in the curvaton scenario as well as some models of baryogenesis. While they leave no observable effects on the cosmic microwave background (CMB) at linear order, they do spatially modulate two-point CMB statistics and can be reconstructed in a manner similar to gravitational lensing. Due to the similarity between the effects of CMB lensing and CIPs, lensing contributes nearly Gaussian random noise to the CIP estimator that approximately doubles the reconstruction noise power. Additionally, the cross correlation between lensing and the integrated Sachs-Wolfe (ISW) effect generates a correlation between the CIP estimator and the temperature field even in the absence of a correlated CIP signal. For cosmic-variance limited temperature measurements out to multipoles $l \leq 2500$, subtracting a fixed lensing bias degrades the detection threshold for CIPs by a factor of $1.3$, whether or not they are correlated with the adiabatic mode.

In the standard inflationary scenario, primordial perturbations are adiabatic. The amplitudes of most types of isocurvature perturbations are generally constrained by current data to be small. If, however, there is a baryon-density perturbation that is compensated by a dark-matter perturbation in such a way that the total matter density is unperturbed, then this compensated isocurvature perturbation (CIP) has no observable consequence in the cosmic microwave background (CMB) at linear order in the CIP amplitude. Here we search for the effects of CIPs on CMB power spectra to quadratic order in the CIP amplitude. An analysis of the Planck temperature data leads to an upper bound $\Delta_{\rm rms}^2 \leq 7.1\times 10^{-3}$, at the 68\% confidence level, to the variance $\Delta_{\rm rms}^2$ of the CIP amplitude. This is then strengthened to $\Delta_{\rm rms}^2\leq 5.0\times 10^{-3}$ if Planck small-angle polarization data are included. A cosmic-variance-limited CMB experiment could improve the $1\sigma$ sensitivity to CIPs to $\Delta^2_{\rm rms} \lesssim 9\times 10^{-4}$. It is also found that adding CIPs to the standard $\Lambda$CDM model can improve the fit of the observed smoothing of CMB acoustic peaks just as much as adding a non-standard lensing amplitude.

In the curvaton scenario, primordial curvature perturbations are produced by a second field that is sub-dominant during inflation. Depending on how the curvaton decays [possibly producing baryon number, lepton number, or cold dark matter (CDM)], mixtures of correlated isocurvature perturbations are produced, allowing the curvaton scenario to be tested using cosmic microwave background (CMB) data. Here, a full range of 27 curvaton-decay scenarios is compared with CMB data, placing limits on the curvaton fraction at decay, $r_D$, and the lepton asymmetry, $\xi_{\rm lep}$. If baryon number is generated by curvaton decay and CDM before (or vice-versa), these limits imply specific predictions for non-Gaussian signatures testable by future CMB experiments and upcoming large-scale-structure surveys.

Primordial fluctuations in the relative number densities of particles, or isocurvature perturbations, are generally well constrained by cosmic microwave background (CMB) data. A less probed mode is the compensated isocurvature perturbation (CIP), a fluctuation in the relative number densities of cold dark matter and baryons. In the curvaton model, a subdominant field during inflation later sets the primordial curvature fluctuation $\zeta$. In some curvaton-decay scenarios, the baryon and cold dark matter isocurvature fluctuations nearly cancel, leaving a large CIP correlated with $\zeta$. This correlation can be used to probe these CIPs more sensitively than the uncorrelated CIPs considered in past work, essentially by measuring the squeezed bispectrum of the CMB for triangles whose shortest side is limited by the sound horizon. Here, the sensitivity of existing and future CMB experiments to correlated CIPs is assessed, with an eye towards testing specific curvaton-decay scenarios. The planned CMB Stage 4 experiment could detect the largest CIPs attainable in curvaton scenarios with more than 3$\sigma$ significance. The significance could improve if small-scale CMB polarization foregrounds can be effectively subtracted. As a result, future CMB observations could discriminate between some curvaton-decay scenarios in which baryon number and dark matter are produced during different epochs relative to curvaton decay. Independent of the specific motivation for the origin of a correlated CIP perturbation, cross-correlation of CIP reconstructions with the primary CMB can improve the signal-to-noise ratio of a CIP detection. For fully correlated CIPs the improvement is a factor of $\sim$2$-$3.

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.

Spectral distortions of the cosmic microwave background (CMB) may become a powerful probe of primordial perturbations at small scales. Existing studies of spectral distortions focus almost exclusively on primordial scalar metric perturbations. Similarly, vector and tensor perturbations should source CMB spectral distortions. In this paper, we give general expressions for the effective heating rate caused by these types of perturbations, including previously neglected contributions from polarization states and higher multipoles. We then focus our discussion on the dissipation of tensors, showing that for nearly scale invariant tensor power spectra, the overall distortion is some six orders of magnitudes smaller than from the damping of adiabatic scalar modes. We find simple analytic expressions describing the effective heating rate from tensors using a quasi-tight coupling approximation. In contrast to adiabatic modes, tensors cause heating without additional photon diffusion and thus over a wider range of scales, as recently pointed out by Ota et. al 2014. Our results are in broad agreement with their conclusions, but we find that small-scale modes beyond k< 2x10^4 Mpc^{-1} cannot be neglected, leading to a larger distortion, especially for very blue tensor power spectra. At small scales, also the effect of neutrino damping on the tensor amplitude needs to be included.

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.

Global, symmetry-breaking phase transitions in the early universe can generate scaling seed networks which lead to metric perturbations. The acoustic waves in the photon-baryon plasma sourced by these metric perturbations, when Silk damped, generate spectral distortions of the cosmic microwave background (CMB). In this work, the chemical potential distortion ($\mu$) due to scaling seed networks is computed and the accompanying Compton $y$-type distortion is estimated. The specific model of choice is the $O(N)$ nonlinear $\sigma$-model for $N\gg 1$, but the results remain the same order of magnitude for other scaling seeds. If CMB anisotropy constraints to the $O(N)$ model are saturated, the resulting chemical potential distortion $\mu \lesssim 2\times 10^{-9}$.

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.

The damping of primordial perturbations at small scales gives rise to distortions of the cosmic microwave background (CMB). Here, the dependence of the distortion on the different types of cosmological initial conditions is explored, covering adiabatic, baryon/cold dark matter isocurvature, neutrino density/velocity isocurvature modes and some mixtures. The radiation transfer functions for each mode are determined and then used to compute the dissipative heating rates and spectral distortion signatures, utilizing both analytic estimates and numerical results from the thermalization code CosmoTherm. Along the way, the early-time super-horizon behavior for the resulting fluid modes is derived in conformal Newtonian gauge, and tight-coupling transfer function approximations are given. CMB spectral distortions caused by different perturbation modes can be estimated using simple k-space window functions which are provided here. Neutrinos carry away some fraction of the primordial perturbation power, introducing an overall efficiency factor that depends on the perturbation type. It is shown that future measurements of the CMB frequency spectrum have the potential to probe different perturbation modes at very small scales (corresponding to wavenumbers 1 Mpc^{-1} < k < few x 10^4 Mpc^{-1}). These constraints are complementary to those obtained at large scales and hence provide an exciting new window to early-universe physics.

Primordial isocurvature fluctuations between photons and either neutrinos or non-relativistic species such as baryons or dark matter are known to be sub-dominant to adiabatic fluctuations. Perturbations in the relative densities of baryons and dark matter (known as compensated isocurvature perturbations, or CIPs), however, are surprisingly poorly constrained. CIPs leave no imprint in the cosmic microwave background (CMB) on observable scales, at least at linear order in their amplitude and zeroth order in the amplitude of adiabatic perturbations. It is thus not yet empirically known if baryons trace dark matter at the surface of last scattering. If CIPs exist, they would spatially modulate the Silk damping scale and acoustic horizon, causing distinct fluctuations in the CMB temperature/polarization power spectra across the sky: this effect is first order in both the CIP and adiabatic mode amplitudes. Here, temperature data from the Wilkinson Microwave Anisotropy Probe (WMAP) are used to conduct the first CMB-based observational search for CIPs, using off-diagonal correlations and the CMB trispectrum. Reconstruction noise from weak lensing and point sources is shown to be negligible for this data set. No evidence for CIPs is observed, and a 95%-confidence upper limit of $1.1\times 10^{-2}$ is imposed to the amplitude of a scale-invariant CIP power spectrum. This limit agrees with CIP sensitivity forecasts for WMAP, and is competitive with smaller scale constraints from measurements of the baryon fraction in galaxy clusters. It is shown that the root-mean-squared CIP amplitude on 5-100 degree scales is smaller than 0.07-0.17 (depending on the scale) at the 95%-confidence level. Temperature data from the Planck satellite will provide an even more sensitive probe for the existence of CIPs, as will the upcoming ACTPol and SPTPol experiments on smaller angular scales.

An improved estimator for the amplitude fnl of local-type non-Gaussianity from the cosmic microwave background (CMB) bispectrum is discussed. The standard estimator is constructed to be optimal in the zero-signal (i.e., Gaussian) limit. When applied to CMB maps which have a detectable level of non-Gaussianity the standard estimator is no longer optimal, possibly limiting the sensitivity of future observations to a non-Gaussian signal. Previous studies have proposed an improved estimator by using a realization-dependent normalization. Under the approximations of a flat sky and a vanishingly thin last-scattering surface, these studies showed that the variance of this improved estimator can be significantly smaller than the variance of the standard estimator when applied to non-Gaussian CMB maps. Here this technique is generalized to the full sky and to include the full radiation transfer function, yielding expressions for the improved estimator that can be directly applied to CMB maps. The ability of this estimator to reduce the variance as compared to the standard estimator in the face of a significant non-Gaussian signal is re-assessed using the full CMB transfer function. As a result of the late time integrated Sachs-Wolfe effect, the performance of the improved estimator is degraded. If CMB maps are first cleaned of the late-time ISW effect using a tracer of foreground structure, such as a galaxy survey or a measurement of CMB weak lensing, the new estimator does remove a majority of the excess variance, allowing a higher significance detection of fnl.

Measurements of cosmic microwave background (CMB) anisotropies constrain isocurvature fluctuations between photons and non-relativistic particles to be sub-dominant to adiabatic fluctuations. Perturbations in the relative number densities of baryons and dark matter, however, are surprisingly poorly constrained. In fact, baryon-density perturbations of fairly large amplitude may exist if they are compensated by dark-matter perturbations, so that the total density remains unchanged. These compensated isocurvature perturbations (CIPs) leave no imprint on the CMB at observable scales, at linear order in their amplitude. B modes in the CMB polarization are generated at reionization through the modulation of the optical depth by CIPs, but this induced polarization is small. The strongest known constraint $\lesssim 10%$ to the CIP amplitude comes from galaxy cluster baryon fractions. Here it is shown that modulation of the baryon density by the CIP at and before the decoupling of Thomson scattering at $z\sim 1100$ gives rise to CMB effects several orders of magnitude larger than those considered before. Polarization B modes are induced, as are correlations between temperature/polarization spherical-harmonic coefficients of different $lm$. It is shown that the CIP field at the surface of last scatter can be measured with these higher-order correlations. The sensitivity of ongoing and future experiments to these fluctuations is estimated. Data from the WMAP, ACT, SPT, and Spider experiments will be sensitive to fluctuations with amplitude $\sim 5-10%$. The Planck satellite and Polarbear experiment will be sensitive to fluctuations with amplitude $\sim 3%$. SPTPol, ACTPol, and future space-based polarization methods will probe amplitudes as low as $\sim 0.4%-0.6%$. In the cosmic variance limit, the lowest amplitude CIPs that could be detected with the CMB are of amplitude $\sim 0.05%$.

Baryon-density perturbations of large amplitude may exist if they are compensated by dark-matter perturbations so that the total density remains unchanged. Big-bang nucleosynthesis and galaxy clusters allow the amplitudes of these compensated isocurvature perturbations (CIPs) to be as large as $\sim10%$. CIPs will modulate the power spectrum of cosmic microwave background (CMB) fluctuations---those due to the usual adiabatic perturbations---as a function of position on the sky. This leads to correlations between different spherical-harmonic coefficients of the temperature/polarization map, and it induces B modes in the CMB polarization. Here, the magnitude of these effects is calculated and techniques to measure them are introduced. While a CIP of this amplitude can be probed on the largest scales with WMAP, forthcoming CMB experiments should improve the sensitivity to CIPs by at least an order of magnitude.

The calculation of a highly accurate cosmological recombination history has been the object of particular attention recently, as it constitutes the major theoretical uncertainty when predicting the angular power spectrum of Cosmic Microwave Background anisotropies. Lyman transitions, in particular the Lyman-alpha line, have long been recognized as one of the bottlenecks of recombination, due to their very low escape probabilities. The Sobolev approximation does not describe radiative transfer in the vicinity of Lyman lines to a sufficient degree of accuracy, and several corrections have already been computed in other works. In this paper, the impact of some previously ignored radiative transfer effects is calculated. First, the effect of Thomson scattering in the vicinity of the Lyman-alpha line is evaluated, using a full redistribution kernel incorporated into a radiative transfer code. The effect of feedback of distortions generated by the optically thick deuterium Lyman-alpha line blueward of the hydrogen line is investigated with an analytic approximation. It is shown that both effects are negligible during cosmological hydrogen recombination. Secondly, the importance of high-lying, non overlapping Lyman transitions is assessed. It is shown that escape from lines above Ly-gamma and frequency diffusion in Ly-beta and higher lines can be neglected without loss of accuracy. Thirdly, a formalism generalizing the Sobolev approximation is developed to account for the overlap of the high-lying Lyman lines, which is shown to lead to negligible changes to the recombination history. Finally, the possibility of a cosmological hydrogen recombination maser is investigated. It is shown that there is no such maser in the purely radiative treatment presented here.

Calculations of cosmological hydrogen recombination are vital for the extraction of cosmological parameters from cosmic microwave background (CMB) observations, and for imposing constraints to inflation and re-ionization. The Planck} mission and future experiments will make high precision measurements of CMB anisotropies at angular scales as small as l~2500, necessitating a calculation of recombination with fractional accuracy of ~10^{-3}. Recent work on recombination includes two-photon transitions from high excitation states and many radiative transfer effects. Modern recombination calculations separately follow angular momentum sublevels of the hydrogen atom to accurately treat non-equilibrium effects at late times (z<900). The inclusion of extremely high-n (n>100) states of hydrogen is then computationally challenging, preventing until now a determination of the maximum n needed to predict CMB anisotropy spectra with sufficient accuracy for Planck. Here, results from a new multi-level-atom code (RecSparse) are presented. For the first time, `forbidden' quadrupole transitions of hydrogen are included, but shown to be negligible. RecSparse is designed to quickly calculate recombination histories including extremely high-n states in hydrogen. Histories for a sequence of values as high as n_max=250 are computed, keeping track of all angular momentum sublevels and energy shells of the hydrogen atom separately. Use of an insufficiently high n_max value (e.g., n_max=64) leads to errors (e.g., 1.8 sigma for Planck) in the predicted CMB power spectrum. Extrapolating errors, the resulting CMB anisotropy spectra are converged to 0.5 sigma at Fisher-matrix level for n_max=128, in the purely radiative case.

There is no direct evidence for radiation domination prior to big-bang nucleosynthesis, and so it is useful to consider how constraints to thermally-produced axions change in non-standard thermal histories. In the low-temperature-reheating scenario, radiation domination begins at temperatures as low as 1 MeV, and is preceded by significant entropy generation. Axion abundances are then suppressed, and cosmological limits to axions are significantly loosened. In a kination scenario, a more modest change to axion constraints occurs. Future possible constraints to axions and low-temperature reheating are discussed.

It is usually assumed that dark matter is produced during the radiation dominated era. There is, however, no direct evidence for radiation domination prior to big-bang nucleosynthesis. Two non-standard thermal histories are considered. In one, the low-temperature-reheating scenario, radiation domination begins as late as 1 MeV, and is preceded by significant entropy generation. Thermal axion relic abundances are then suppressed, and cosmological limits to axions are loosened. For reheating temperatures less than 35 MeV, the large-scale structure limit to the axion mass is lifted. The remaining constraint from the total density of matter is significantly relaxed. Constraints are also relaxed for higher reheating temperatures. In a kination scenario, a more modest change to cosmological axion constraints is obtained. Future possible constraints to axions and low-temperature reheating from the helium abundance and next-generation large-scale-structure surveys are discussed.

We consider a linearized, effective quantum theory of gravitation in which gravity weakens at energies higher than ~10^-3 eV in order to accommodate the apparent smallness of the cosmological constant. Such a theory predicts departures from the static Newtonian inverse-square force law on distances below ~0.05 mm. However, we show that such a modification also leads to changes in the long-range behavior of gravity and is inconsistent with observed gravitational lenses.

A search for optical line emission from the two-photon decay of relic axions was conducted in the galaxy clusters Abell 2667 and 2390, using spectra from the VIMOS (Visible Multi-Object Spectrograph) integral field unit at the Very Large Telescope. New upper limits to the two-photon coupling of the axion are derived, and are at least a factor of 3 more stringent than previous upper limits in this mass window. The improvement follows from larger collecting area, integration time, and spatial resolution, as well as from improvements in signal to noise and sky subtraction made possible by strong-lensing mass models of these clusters. The new limits either require that the two-photon coupling of the axion be extremely weak or that the axion mass window between 4.5 eV and 7.7 eV be closed. Implications for sterile-neutrino dark matter are discussed briefly also.