Voids have emerged as a novel probe of cosmology and large-scale structure. These regions of extreme underdensity are sensitive to physics beyond the standard model of cosmology, and can potentially be used as a testing ground to constrain new physics. We present the first determination of the linear void bias measured in separate universe simulations. Our methods are validated by comparing the separate universe response bias with the clustering bias of voids. We find excellent agreement between the two methods for voids identified in the halo field and the down-sampled dark matter field. For voids traced by halos, we identify two different contributions to the bias. The first is due to the bias of the underlying halo field used to identify voids, while the second we attribute to the dynamical impact of long-wavelength density perturbations on void formation and expansion. By measuring these contributions individually, we demonstrate that their sum is consistent with the total void bias. We also measure the void profiles in our simulations, and determine their separate universe response. These can be interpreted as the sensitivity of the profiles to the background density of the Universe.

We study the implications of the recently proposed Trans-Planckian Censorship Conjecture (TCC) for early universe cosmology and in particular inflationary cosmology. The TCC leads to the conclusion that if we want inflationary cosmology to provide a successful scenario for cosmological structure formation, the energy scale of inflation has to be lower than $10^9$ GeV. Demanding the correct amplitude of the cosmological perturbations then forces the generalized slow-roll parameter $\epsilon$ of the model to be very small ($<10^{-31}$). This leads to the prediction of a negligible amplitude of primordial gravitational waves. For slow-roll inflation models, it also leads to severe fine tuning of initial conditions.

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.

In a universe with quintessence isocurvature, or perturbations in dark energy that are independent from the usual curvature perturbations, structure formation is changed qualitatively. The existence of two independent fields, curvature and isocurvature, causes the growth rate of matter perturbations to depend on their initial conditions. The quintessence perturbations cause their growth to depend on scale. We perform the first separate universe simulations for this cosmology. We demonstrate that the power spectrum response and the halo bias depend on scale and initial conditions and that the presence of the isocurvature mode changes the mapping from these quantities to the halo auto- and cross-power spectra, and the squeezed-limit bispectrum. We compare the bias to several models, finding reasonable agreement with both a power-spectrum-response model with one free parameter and a model that fits two independent bias parameters for curvature and isocurvature sourced fluctuations. We also verify that simulation responses to pure isocurvature and pure curvature modes can be linearly combined to reproduce responses with different ratios of isocurvature and curvature. This allows our results to be used to predict the halo power spectrum and stochasticity with arbitrary large-scale curvature and isocurvature power spectra. In an appendix, we study the generation of quintessence isocurvature during inflation and show that a modified kinetic term is typically required to produce observable isocurvature modes in a field with $w_Q\approx -1$.

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

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

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.

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.

Recent advances in cosmic observations have brought us to the verge of discovery of the absolute scale of neutrino masses. Nonzero neutrino masses are known evidence of new physics beyond the Standard Model. Our understanding of the clustering of matter in the presence of massive neutrinos has significantly improved over the past decade, yielding cosmological constraints that are tighter than any laboratory experiment, and which will improve significantly over the next decade, resulting in a guaranteed detection of the absolute neutrino mass scale.

Using $N$-body simulations with massive neutrino density perturbations, we detect the scale-dependent linear halo bias with high significance. This is the first time that this effect is detected in simulations containing neutrino density perturbations on all scales, confirming the same finding from separate universe simulations. The scale dependence is the result of the additional scale in the system, i.e. the massive neutrino free-streaming length, and it persists even if the bias is defined with respect to the cold dark matter plus baryon (instead of total matter) power spectrum. The separate universe approach provides a good model for the scale-dependent linear bias, and the effect is approximately $0.25f_\nu$ and $0.43f_\nu$ for halos with bias of 1.7 and 3.5, respectively. While the size of the effect is small, it is ${\it not}$ insignificant in terms of $f_\nu$ and should therefore be included to accurately constrain neutrino mass from clustering statistics of biased tracers. More importantly, this feature is a distinct signature of free-streaming particles and cannot be mimicked by other components of the standard cosmological model.

The acoustic peaks in the angular power spectrum of cosmic microwave background (CMB) temperature and polarization anisotropies play an important role as a probe of the nature of new relativistic particles contributing to the radiation density in the early universe, parametrized by $\Delta N_{eff}$. The amplitude and phase of the acoustic oscillations provide information about whether the extra species are free-streaming particles, like neutrinos, or tightly-coupled, like the photons, during eras probed by the CMB. On the other hand, some extensions of the Standard Model produce new relativistic particles that decouple from their own non-gravitational interactions after neutrinos, but prior to photons. We study the signature of new relativistic species that decouple during this intermediate epoch. We argue that the decoupling species will cause a scale-dependent change in the amplitude and phase shift of the acoustic oscillations, different from the usual constant shifts on small scales. For intermediate decoupling times, the phase and amplitude shifts depend not only on $\Delta N_{eff}$ but the redshift $z_{dec,X}$ at which the new species decoupled. For $\Delta N_{eff} >0.334$, a Stage IV CMB experiment could determine $N_{eff}$ at the percent level and $z_{dec,X}$ at the $\sim 10\%$ level. For smaller values, $\Delta N_{eff}\sim 0.1$, constraints on $z_{dec,X}$ weaken but remain $\sim 20-50\%$ for $z_{dec,X} \sim \mathcal{O}(10^3-10^4)$. As an application, we study the contributions to $\Delta N_{eff}$ and determine the $z_{dec,X}$ values for simple implementations of the so-called $N$naturalness model.

Cosmic background neutrinos have a large velocity dispersion, which causes the evolution of long-wavelength density perturbations to depend on scale. This scale-dependent growth leads to the well-known suppression in the linear theory matter power spectrum that is used to probe neutrino mass. In this paper, we study the impact of long-wavelength density perturbations on small-scale structure formation. By performing separate universe simulations where the long-wavelength mode is absorbed into the local expansion, we measure the responses of the cold dark matter (CDM) power spectrum and halo mass function, which correspond to the squeezed-limit bispectrum and halo bias. We find that the scale-dependent evolution of the long-wavelength modes causes these quantities to depend on scale and provide simple expressions to model them in terms of scale and the amount of massive neutrinos. Importantly, this scale-dependent bias reduces the suppression in the linear halo power spectrum due to massive neutrinos by 13 and 26% for objects of bias $\bar{b}=2$ and $\bar{b} \gg1$, respectively. We demonstrate with high statistical significance that the scale-dependent halo bias ${\it cannot}$ be modeled by the CDM and neutrino density transfer functions at the time when the halos are identified. This reinforces the importance of the temporal nonlocality of structure formation, especially when the growth is scale dependent.

By absorbing fluctuations into a local background, separate universe simulations provide a powerful technique to characterize the response of small-scale observables to the long-wavelength density fluctuations, for example those of the power spectrum and halo mass function which lead to the squeezed-limit $n$-point function and halo bias, respectively. Using quintessence dark energy as the paradigmatic example, we extend these simulation techniques to cases where non-gravitational forces in other sectors establish a Jeans scale across which the growth of density fluctuations becomes scale dependent. By characterizing the separate universes with matching background expansion histories, we show that the power spectrum and mass function responses depend on whether the long-wavelength mode is above or below the Jeans scale. Correspondingly, the squeezed bispectrum and halo bias also become scale dependent. Models of bias that are effectively local in the density field at a single epoch, initial or observed, cannot describe this effect which highlights the importance of temporal nonlocality in structure formation. Validated by these quintessence tests, our techniques are applicable to a wide range of models where the complex dynamics of additional fields affect the clustering of matter in the linear regime and it would otherwise be difficult to simulate their impact in the nonlinear regime.

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.

Measuring the absolute scale of the neutrino masses is one of the most exciting opportunities available with near-term cosmological datasets. Two quantities that are sensitive to neutrino mass, scale-dependent halo bias $b(k)$ and the linear growth parameter $f(k)$ inferred from redshift-space distortions, can be measured without cosmic variance. Unlike the amplitude of the matter power spectrum, which always has a finite error, the error on $b(k)$ and $f(k)$ continues to decrease as the number density of tracers increases. This paper presents forecasts for statistics of galaxy and lensing fields that are sensitive to neutrino mass via $b(k)$ and $f(k)$. The constraints on neutrino mass from the auto- and cross-power spectra of spectroscopic and photometric galaxy samples are weakened by scale-dependent bias unless a very high density of tracers is available. In the high density limit, using multiple tracers allows cosmic-variance to be beaten and the forecasted errors on neutrino mass shrink dramatically. In practice, beating the cosmic variance errors on neutrino mass with $b(k)$ will be a challenge, but this signal is nevertheless a new probe of neutrino effects on structure formation that is interesting in its own right.

The separate universe technique provides a means of establishing consistency relations between short wavelength observables and the long wavelength matter density fluctuations within which they evolve by absorbing the latter into the cosmological background. We extend it to cases where non-gravitational forces introduce a Jeans scale in other species like dynamical dark energy or massive neutrinos. The technique matches the synchronous gauge matter density fluctuations to the local expansion using the acceleration equation and accounts for the temporal nonlocality and scale dependence of the long wavelength response of small scale matter observables, e.g. the nonlinear power spectrum, halo abundance and the implied halo bias, and $N$-point correlation functions. Above the Jeans scale, the local Friedmann equation relates the expansion to real energy densities and a curvature that is constant in comoving coordinates. Below the Jeans scale, the curvature evolves and acts like a fake density component. In all cases, the matter evolution on small scales is correctly modeled as we illustrate using scalar field dark energy with adiabatic or isocurvature initial conditions across the Jeans scale set by its finite sound speed.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

The large-scale distribution of cold dark matter halos is generally assumed to trace the large-scale distribution of matter. In a universe with multiple types of matter fluctuations, as is the case with massive neutrinos, the relation between the halo field and the matter fluctuations may be more complicated. We develop a method for calculating the bias factor relating fluctuations in the halo number density to fluctuations in the mass density in the presence of multiple fluctuating components of the energy density. In the presence of massive neutrinos we find a small but pronounced feature in the halo bias near the neutrino free-streaming scale. The neutrino feature is a small step with amplitude that increases with halo mass and neutrino mass density. The scale-dependent halo bias lessens the suppression of the small-scale halo power spectrum and should therefore weaken constraints on neutrino mass from the galaxy auto-power spectrum and correlation function. On the other hand, the feature in the bias is itself a novel signature of massive neutrinos that can be studied independently.

The abundance of massive dark matter halos hosting galaxy clusters provides an important test of the masses of relic neutrino species. The dominant effect of neutrino mass is to lower the typical amplitude of density perturbations that eventually form halos, but for neutrino masses $> 0.4eV$ the threshold for halo formation can be changed significantly as well. We study the spherical collapse model for halo formation in cosmologies with neutrino masses in the range $m_{\nu i} =0.05eV$- $1eV$ and find that halo formation is differently sensitive to $\Omega_\nu$ and $m_\nu$. That is, different neutrino hierarchies with common $\Omega_\nu$ are in principle distinguishable. The added sensitivity to $m_{\nu}$ is small but potentially important for scenarios with heavier sterile neutrinos. Massive neutrinos cause the evolution of density perturbations to be scale-dependent at high redshift which complicates the usual mapping between the collapse threshold and halo abundance. We propose one way of handling this and compute the correction to the halo mass function within this framework. For $\sum m_{\nu i} < 0.3eV$, our prescription for the halo abundance is only $\sim< 15\%$ different than the standard calculation. However for larger neutrino masses the differences approach $50-10\%$ which, if verified by simulations, could alter neutrino mass constraints from cluster abundance.

Cold dark matter halos form within a smoothly distributed background of relic neutrinos -- at least some of which are massive and non-relativistic at late times. We calculate the accumulation of massive neutrinos around spherically collapsing cold dark matter halos in a cosmological background. We identify the physical extent of the "neutrino halo" in the spherical collapse model, which is large in comparison with the virial radius of the dark matter, and conditions under which neutrinos reaching the cold dark matter halo will remain bound to the halo at late times. We calculate the total neutrino mass and bound neutrino mass associated with isolated spherical halos for several neutrino mass hierarchies and provide fitting formulae for these quantities in terms of the cold dark matter halo mass and the masses of the individual neutrino species.

If the entire post-inflationary patch is large compared to our Hubble volume even a small level of non-Gaussianity can cause statistics of the primordial curvature field in our Hubble volume to be biased by mode-coupling. We explicitly compute the variation of locally measured statistics of the primordial curvature $\zeta$ from non-Gaussian mode coupling within a specific inflationary scenario: the curvaton model with a quadratic curvaton potential. This "super cosmic variance" is calculated in two ways: (i) as a super observer who has access to the curvature perturbation field across the entire post-inflationary patch and therefore sees local statistics as biased by mode coupling and (ii) as a local observer who sees the statistics of $\zeta$ determined by the local values of quantities in their Hubble patch. The two calculations agree and show that in the quadratic curvaton model patch-to-patch differences in statistics of $\zeta$ can be interpreted entirely as a shift in the value of the curvaton field at freeze out. Applying the same arguments to single-field slow-roll inflation gives a simple picture of how non-Gaussian mode-coupling between the curvature perturbations on very different physical scales must vanish in the attractor limit.

Local-type primordial non-Gaussianity couples statistics of the curvature perturbation \zeta on vastly different physical scales. Because of this coupling, statistics (i.e. the polyspectra) of \zeta in our Hubble volume may not be representative of those in the larger universe -- that is, they may be biased. The bias depends on the local background value of \zeta, which includes contributions from all modes with wavelength k ~< H_0 and is therefore enhanced if the entire post-inflationary patch is large compared with our Hubble volume. We study the bias to locally-measured statistics for general local-type non-Gaussianity. We consider three examples in detail: (i) the usual fNL, gNL model, (ii) a strongly non-Gaussian model with \zeta ~ \zeta_G^p, and (iii) two-field non-Gaussian initial conditions. In each scenario one may generate statistics in a Hubble-size patch that are weakly Gaussian and consistent with observations despite the fact that the statistics in the larger, post-inflationary patch look very different.

Large-scale clustering of highly biased tracers of large-scale structure has emerged as one of the best observational probes of primordial non-Gaussianity of the local type (i.e. f_{NL}^{local}). This type of non-Gaussianity can be generated in multifield models of inflation such as the curvaton model. Recently, Tseliakhovich, Hirata, and Slosar showed that the clustering statistics depend qualitatively on the ratio of inflaton to curvaton power \xi after reheating, a free parameter of the model. If \xi is significantly different from zero, so that the inflaton makes a non-negligible contribution to the primordial adiabatic curvature, then the peak-background split ansatz predicts that the halo bias will be stochastic on large scales. In this paper, we test this prediction in N-body simulations. We find that large-scale stochasticity is generated, in qualitative agreement with the prediction, but that the level of stochasticity is overpredicted by ~30%. Other predictions, such as \xi independence of the halo bias, are confirmed by the simulations. Surprisingly, even in the Gaussian case we do not find that halo model predictions for stochasticity agree consistently with simulations, suggesting that semi-analytic modeling of stochasticity is generally more difficult than modeling halo bias.

Models of inflation in which non-Gaussianity is generated outside the horizon, such as curvaton models, generate distinctive higher-order correlation functions in the CMB and other cosmological observables. Testing for violation of the Suyama-Yamaguchi inequality tauNL >= (6/5 fNL)^2, where fNL and tauNL denote the amplitude of the three-point and four-point functions in certain limits, has been proposed as a way to distinguish qualitative classes of models. This inequality has been proved for a wide range of models, but only weaker versions have been proved in general. In this paper, we give a proof that the Suyama-Yamaguchi inequality is always satisfied. We discuss scenarios in which the inequality may appear to be violated in an experiment such as Planck, and how this apparent violation should be interpreted. We analyze a specific example, the "ungaussiton" model, in which leading-order scaling relations suggest that the Suyama-Yamaguchi inequality is eventually violated, and show that the inequality always holds.

Cosmic inflation provides a mechanism for generating the early density perturbations that seeded the large-scale structures we see today. Primordial non-Gaussianity is among the most promising of few observational tests of physics at this epoch. At present non-Gaussianity is best constrained by the cosmic microwave background, but in the near term large-scale structure data may be competitive so long as the effects of primordial non-Gaussianity can be modeled through the non-linear process of structure formation. We discuss recent work modeling effects of a few types of primordial non-Gaussianity on the large-scale halo clustering and the halo mass function. More specifically, we compare analytic and N-body results for two variants of the curvaton model of inflation: (i) a "tauNL" scenario in which the curvaton and inflaton contribute equally to the primordial curvature perturbation and (ii) a "gNL" model where the usual quadratic fNL term in the potential cancels, but a large cubic term remains.

Primordial non-Gaussianity has emerged as one of the most promising probes of the inflationary epoch. While the cosmic microwave background and large-scale halo bias currently provide the most stringent constraints on the non-Gaussian parameter fNL, the abundance of dark matter halos is a complementary probe which may allow tests of Gaussianity which are independent of the precise form of non-Gaussian initial conditions. We study the halo mass function in N-body simulations with a range of non-Gaussian initial conditions. In addition to the usual fNL model, we consider gNL Phi^3-type non-Gaussianity and models where the 4-point amplitude tauNL is an independent parameter. We introduce a new analytic form for the halo mass function in the presence of primordial non-Gaussianity, the "log-Edgeworth" mass function, and find good agreement with the N-body simulations. The log-Edgeworth mass function introduces no free parameters and can be constructed from first principles for any model of primordial non-Gaussianity.

A wide range of multifield inflationary models generate non-Gaussian initial conditions in which the initial adiabatic fluctuation is of the form (zeta_G + g_{NL} zeta_G^3). We study halo clustering in these models using two different analytic methods: the peak-background split framework, and brute force calculation in a barrier crossing model, obtaining agreement between the two. We find a simple, theoretically motivated expression for halo bias which agrees with N-body simulations and can be used to constrain g_{NL} from observations. We discuss practical caveats to constraining g_{NL} using only observable properties of a tracer population, and argue that constraints obtained from populations whose observed bias is <~ 2.5 are generally not robust to uncertainties in modeling the halo occupation distribution of the population.

One of the primary scientific targets of current and future CMB polarization experiments is the search for a stochastic background of gravity waves in the early universe. As instrumental sensitivity improves, the limiting factor will eventually be B-mode power generated by gravitational lensing, which can be removed through use of so-called delensing algorithms. We forecast prospects for delensing using lensing maps which are obtained externally to CMB polarization: either from large-scale structure observations, or from high-resolution maps of CMB temperature. We conclude that the forecasts in either case are not encouraging, and that significantly delensing large-scale CMB polarization requires high-resolution polarization maps with sufficient sensitivity to measure the lensing B-mode. We also present a simple formalism for including delensing in CMB forecasts which is computationally fast and agrees well with Monte Carlos.

In Lyman-alpha forest measurements it is generally assumed that quasars are mere background light sources which are uncorrelated with the forest. Gravitational lensing of the quasars violates this assumption. This effect leads to a measurement bias, but more interestingly it provides a valuable signal. The lensing signal can be extracted by correlating quasar magnitudes with the flux power spectrum and with the flux decrement. These correlations will be challenging to measure but their detection provides a direct measure of how features in the Lyman-alpha forest trace the underlying mass density field. Observing them will test the fundamental hypothesis that fluctuations in the forest are predominantly driven by fluctuations in mass, rather than in the ionizing background, helium reionization or winds. We discuss ways to disentangle the lensing signal from other sources of such correlations, including dust, continuum and background residuals. The lensing-induced measurement bias arises from sample selection: one preferentially collects spectra of magnified quasars which are behind overdense regions. This measurement bias is ~0.1-1% for the flux power spectrum, optical depth and the flux probability distribution. Since the effect is systematic, quantities such as the amplitude of the flux power spectrum averaged across scales should be interpreted with care.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Gravitational lensing of the cosmic microwave background by large-scale structure in the late universe is both a source of cosmological information and a potential contaminant of primordial gravity waves. Because lensing imprints growth of structure in the late universe on the CMB, measurements of CMB lensing will constrain parameters to which the CMB would not otherwise be sensitive, such as neutrino mass. If the instrumental noise is sufficiently small (<~ 5 uK-arcmin), the gravitational lensing contribution to the large-scale B-mode will be the limiting source of contamination when constraining a stochastic background of gravity waves in the early universe, one of the most exciting prospects for future CMB polarization experiments. High-sensitivity measurements of small-scale B-modes can reduce this contamination through a lens reconstruction technique that separates the lensing and primordial contributions to the B-mode on large scales. A fundamental design decision for a future CMB polarization experiment such as CMBpol is whether to have coarse angular resolution so that only the large-scale B-mode (and the large-scale E-mode from reionization) is measured, or high resolution to additionally measure CMB lensing. The purpose of this white paper is to evaluate the science case for CMB lensing in polarization: constraints on cosmological parameters, increased sensitivity to the gravity wave B-mode via lens reconstruction, expected level of contamination from non-CMB foregrounds, and required control of beam systematics.

We develop a systematic derivation for the Limber approximation to the angular cross-power spectrum of two random fields, as a series expansion in 1/(\ell+1/2). This extended Limber approximation can be used to test the accuracy of the Limber approximation and to improve the rate of convergence at large \ell's. We show that the error in ordinary Limber approximation is O(1/\ell^2). We also provide a simple expression for the second order correction to the Limber formula, which improves the accuracy to O(1/\ell^4). This correction can be especially useful for narrow redshift bins, or samples with small redshift overlap, for which the zeroth order Limber formula has a large error. We also point out that using \ell instead of (\ell+1/2), as is often done in the literature, spoils the accuracy of the approximation to O(1/\ell).

The detection of primordial non-Gaussianity could provide a powerful means to test various inflationary scenarios. Although scale-invariant non-Gaussianity (often described by the $f_{NL}$ formalism) is currently best constrained by the CMB, single-field models with changing sound speed can have strongly scale-dependent non-Gaussianity. Such models could evade the CMB constraints but still have important effects at scales responsible for the formation of cosmological objects such as clusters and galaxies. We compute the effect of scale-dependent primordial non-Gaussianity on cluster number counts as a function of redshift, using a simple ansatz to model scale-dependent features. We forecast constraints on these models achievable with forthcoming data sets. We also examine consequences for the galaxy bispectrum. Our results are relevant for the Dirac-Born-Infeld model of brane inflation, where the scale-dependence of the non-Gaussianity is directly related to the geometry of the extra dimensions.

In paper I of this series we discuss how magnification bias distorts the 3D correlation function by enhancing the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. This lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification {\it and} magnification-magnification correlations. Here we extend the discussion to the power spectrum and also to redshift space. In real space, pairs oriented close to the LOS direction are not protected against nonlinearity even if the pair separation is large; this is because nonlinear fluctuations can enter through gravitational lensing at a small transverse separation (or i.e. impact parameter). The situation in Fourier space is different: by focusing on a small wavenumber $k$, as is usually done, linearity is guaranteed because both the LOS and transverse wavenumbers must be small. This is why magnification distortion of the galaxy correlation appears less severe in Fourier space. Nonetheless, the effect is non-negligible, especially for the transverse Fourier modes, and should be taken into account in interpreting precision measurements of the galaxy power spectrum, for instance those that focus on the baryon oscillations. The lensing induced anisotropy of the power spectrum has a shape that is distinct from the more well known redshift space anisotropies due to peculiar motions and the Alcock-Paczynski effect. The lensing anisotropy is highly localized in Fourier space while redshift space distortions are more spread out. This means that one could separate the magnification bias component in real observations, implying that potentially it is possible to perform a gravitational lensing measurement without measuring galaxy shapes.

It is well known that magnification bias, the modulation of galaxy or quasar source counts by gravitational lensing, can change the observed angular correlation function. We investigate magnification-induced changes to the shape of the observed correlation function w(\theta) and the angular power spectrum C_{\ell}, paying special attention to the matter-radiation equality peak and the baryon wiggles. Lensing mixes the correlation function of the source galaxies with the matter correlation at the lower redshifts of the lenses. Since the lenses probe structure nearer to the observer, the angular scale dependence of the lensing terms is different from that of the sources, thus the observed correlation function is distorted. We quantify how the lensing corrections depend on the width of the selection function, the galaxy bias b, and the number count slope s. The correction increases with redshift and larger corrections are present for sources with steep number count slopes and/or broad redshift distributions. The most drastic changes to C_{\ell} occur for measurements at z >~1.5 and \ell <~ 100. For the source distributions we consider, magnification bias can shift the matter-radiation equality scale by 1-6% at z ~ 1.5 and by z ~ 3.5 the shift can be as large as 30%. The baryon bump in \theta^2w(\theta) is shifted by <~ 1% and the width is typically increased by ~10%. Shifts of >~ 0.5% and broadening of >~ 20% occur only for very broad selection functions and/or galaxies with (5s-2)/b>~2. However, near the baryon bump the magnification correction is not constant but a gently varying function which depends on the source population. Depending on how the w(\theta) data is fitted, this correction may need to be accounted for when using the baryon acoustic scale for precision cosmology.

It is well known gravitational lensing, mainly via magnification bias, modifies the observed galaxy/quasar clustering. Such discussions have largely focused on the 2D angular correlation. Here and in a companion paper (Paper II) we explore how magnification bias distorts the 3D correlation function and power spectrum, as first considered by Matsubara. The interesting point is: the distortion is anisotropic. Magnification bias preferentially enhances the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. For example at LOS separation of ~100 Mpc/h, where the intrinsic galaxy-galaxy correlation is rather weak, the observed correlation can be enhanced by lensing by a factor of a few, even at a modest redshift of z ~ 0.35. The opportunity: this lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification and magnification-magnification correlations, without measuring galaxy shapes. The anisotropy is distinguishable from the well known distortion due to peculiar motions, as will be discussed in Paper II. The challenge: the magnification distortion of the galaxy correlation must be accounted for in interpreting data as precision improves. For instance, the ~100 Mpc/h baryon acoustic oscillation scale in the correlation function is shifted by up to ~3% in the LOS orientation, and up to ~0.6% in the monopole, depending on the galaxy bias, redshift and number count slope. The corresponding shifts in the inferred Hubble parameter and angular diameter distance, if ignored, could significantly bias measurements of the dark energy equation of state. Lastly, magnification distortion offers a plausible explanation for the well known excess correlations seen in pencil beam surveys.

Integrated Sachs-Wolfe (ISW) measurements, which involve cross-correlating the CMB with the foreground large-scale structure (e.g. galaxies/quasars), have proven to be an interesting probe of dark energy. We show that magnification bias, which is the inevitable modulation of the foreground number counts by gravitational lensing, alters both the shape and amplitude of the observed ISW signal. This is true especially at high redshifts because (1) the intrinsic galaxy-temperature signal diminishes greatly back in the matter dominated era, (2) the lensing efficiency increases with redshift and (3) the number count slope generally steepens with redshift in a magnitude limited sample. At z >~ 2, the magnification-temperature correlation dominates over the intrinsic galaxy-temperature correlation and causes the observed ISW signal to increase with z, despite dark energy subdominance -- a result of the fact that magnification probes structures between the observer and the sources. Ignoring magnification bias can then lead to erroneous conclusions about dark energy. While the lensing modulation opens up an interesting high z window for ISW measurements, high z measurements are not expected to add much new information to low z ones if dark energy is the cosmological constant. This is because lensing introduces significant covariance across redshifts. The most compelling reason to pursue high z ISW measurements is to look for a potential surprise such as early dark energy domination or the signature of modified gravity. We conclude with a discussion of existing measurements, the highest z of which is at the margin of being sensitive to magnification bias. We also develop a formalism which might be of general interest: to predict biases in estimating parameters when certain physical effects are ignored in interpreting data.

The Advanced Liquid-Mirror Probe of Asteroids, Cosmology and Astrophysics (ALPACA) is a proposed 8-meter liquid mirror telescope surveying 1000 square degree of the southern-hemisphere sky. It will be a remarkably simple and inexpensive telescope, that nonetheless will deliver a powerful sample of optical data for studying dark energy. The bulk of the cosmological data consists of nightly, high signal-to-noise, multiband light curves of SN Ia. At the end of the three-years run ALPACA is expected to collect 100,000 SNe Ia up to z~1. This will allow us to reduce present systematic uncertainties affecting the standard-candle relation. The survey will also provide several other datasets such as the detection of baryon acoustic oscillations in the matter power spectrum and shear weak lensing measurements. In this preliminary analysis we forecast constraints on dark energy parameters from SN Ia and baryon acoustic oscillations. The combination of these two datasets will provide competitive constraints on the dark energy parameters under minimal prior assumptions. Further studies are needed to address the accuracy of weak lensing measurements.