Compensated isocurvature perturbations (CIPs) are opposite spatial fluctuations in the baryon and dark matter density. They can be generated for example in the curvaton model in the early Universe but are difficult to observe because their gravitational imprint nearly cancels. We therefore propose a new measurement method by searching for a spatial modulation of the baryon acoustic oscillation (BAO) scale that CIPs induce. We find that for a Euclid-like survey the sensitivity is marginally better than the WMAP cosmic microwave background (CMB) constraint, which exploits the CIP-induced modulation of the CMB sound horizon. For a cosmic-variance limited BAO survey using emission-line galaxies up to $z\sim7$ the sensitivity is between stage 3 and stage 4 CMB experiments. These results include using CIP-galaxy cross-correlations, which improves the sensitivity by a factor of $\sim2-3$ for correlated CIPs. The method could be further improved with an optimal estimator, similarly to the CMB, and could provide a useful cross-check of other CIP probes. Finally, if CIPs exist, they can bias cosmological measurements made assuming no CIPs. In particular, they can act as a super-sample fluctuation of the baryon density and bias measurements of the BAO scale. For modern BAO surveys, the largest 2$\sigma$ CIP fluctuation allowed by Planck's 95\% bound could bias BAO measurements of $H(z)$ by 2.2\%, partially reducing the tension with the local $H_0$ measurements from 3.1$\sigma$ to 2.3$\sigma$.

A common problem in cosmology is to integrate the product of two or more spherical Bessel functions (sBFs) with different configuration-space arguments against the power spectrum or its square, weighted by powers of wavenumber. Naively computing them scales as $N_{\rm g}^{p+1}$ with $p$ the number of configuration space arguments and $N_{\rm g}$ the grid size, and they cannot be done with Fast Fourier Transforms (FFTs). Here we show that by rewriting the sBFs as sums of products of sine and cosine and then using the product to sum identities, these integrals can then be performed using 1-D FFTs with $N_{\rm g} \log N_{\rm g}$ scaling. This "rotation" method has the potential to accelerate significantly a number of calculations in cosmology, such as perturbation theory predictions of loop integrals, higher order correlation functions, and analytic templates for correlation function covariance matrices. We implement this approach numerically both in a free-standing, publicly-available \textsc{Python} code and within the larger, publicly-available package \texttt{mcfit}. The rotation method evaluated with direct integrations already offers a factor of 6-10$\times$ speed-up over the naive approach in our test cases. Using FFTs, which the rotation method enables, then further improves this to a speed-up of $\sim$$1000-3000\times$ over the naive approach. The rotation method should be useful in light of upcoming large datasets such as DESI or LSST. In analysing these datasets recomputation of these integrals a substantial number of times, for instance to update perturbation theory predictions or covariance matrices as the input linear power spectrum is changed, will be one piece in a Monte Carlo Markov Chain cosmological parameter search: thus the overall savings from our method should be significant.

The forthcoming generation of wide-field galaxy surveys will probe larger volumes and galaxy densities, thus allowing for a much larger signal-to-noise ratio for higher-order clustering statistics, in particular the galaxy bispectrum. Extracting this information, however, is more challenging than using the power spectrum due to more complex theoretical modeling, as well as significant computational cost of evaluating the bispectrum signal and the error budget. To overcome these challenges, several proxy statistics have been proposed in the literature, which partially or fully capture the information in the bispectrum, while being computationally less expensive than the bispectrum. One such statistics are {\it weighted skew-spectra}, which are cross-spectra of the density field and appropriately weighted quadratic fields. Using Fisher forecasts, we show that the information in these skew-spectra is equivalent to that in the bispectrum for parameters that appear as amplitudes in the bispectrum model, such as galaxy bias parameters or the amplitude of primordial non-Gaussianity. We consider three shapes of the primordial bispectrum: local, equilateral and that due to massive particles with spin two during inflation. To obtain constraints that match those from a measurement of the full bispectrum, we find that it is crucial to account for the full covariance matrix of the skew-spectra.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

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.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

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

In this paper we test the perturbative halo bias model at the field level. The advantage of this approach is that any analysis can be done without sample variance if the same initial conditions are used in simulations and perturbation theory calculations. We write the bias expansion in terms of modified bias operators in Eulerian space, designed such that the large bulk flows are automatically resummed and not treated perturbatively. Using these operators, the bias model accurately matches the Eulerian density of halos in N-body simulations. The mean-square model error is close to the Poisson shot noise for a wide range of halo masses and it is rather scale-independent, with scale-dependent corrections becoming relevant at the nonlinear scale. In contrast, for linear bias the mean-square model error can be higher than the Poisson prediction by factors of up to a few on large scales, and it becomes scale dependent already in the linear regime. We show that by weighting simulated halos by their mass, the mean-square error of the model can be further reduced by up to an order of magnitude, or by a factor of two when including $60\%$ mass scatter. We also test the Standard Eulerian bias model using the nonlinear matter field measured from simulations and show that it leads to a larger and more scale-dependent model error than the bias expansion based on perturbation theory. These results may be of particular relevance for cosmological inference methods that use a likelihood of the biased tracer at the field level, or for initial condition and BAO reconstruction that requires a precise estimate of the large-scale potential from the biased tracer density.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

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 Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

With low redshift probes reaching unprecedented precision, uncertainty of the CMB optical depth is expected to be the limiting factor for future cosmological neutrino mass constraints. In this paper, we discuss to what extent combinations of CMB lensing and galaxy surveys measurements at low redshifts $z\sim 0.5-5$ will be able to make competitive neutrino mass measurements without relying on any optical depth constraints. We find that the combination of LSST galaxies and CMB-S4 lensing should be able to achieve constraints on the neutrino mass sum of 25meV without optical depth information, an independent measurement that is competitive with or slightly better than the constraint of 30meV possible with CMB-S4 and present-day optical depth measurements. These constraints originate both in structure growth probed by cross-correlation tomography over a wide redshift range as well as, most importantly, the shape of the galaxy power spectrum measured over a large volume. We caution that possible complications such as higher-order biasing and systematic errors in the analysis of high redshift galaxy clustering are only briefly discussed and may be non-negligible. Nevertheless, our results show that new kinds of high-precision neutrino mass measurements at and beyond the present-day optical depth limit may be possible.

The lensing convergence measurable with future CMB surveys like CMB-S4 will be highly correlated with the clustering observed by deep photometric large scale structure (LSS) surveys such as the LSST, with cross-correlation coefficient as high as 95\%. This will enable use of sample variance cancellation techniques to determine cosmological parameters, and use of cross-correlation measurements to break parameter degeneracies. Assuming large sky overlap between CMB-S4 and LSST, we show that a joint analysis of CMB-S4 lensing and LSST clustering can yield very tight constraints on the matter amplitude $\sigma_8(z)$, halo bias, and $f_\mathrm{NL}$, competitive with the best stage IV experiment predictions, but using complementary methods, which may carry different and possibly lower systematics. Having no sky overlap between experiments degrades the precision of $\sigma_8(z)$ by a factor of 20, and that of $f_\mathrm{NL}$ by a factor of 1.5 to 2. Without CMB lensing, the precision always degrades by an order of magnitude or more, showing that a joint analysis is critical. Our results also suggest that CMB lensing in combination with LSS photometric surveys is a competitive probe of the evolution of structure in the redshift range $z\simeq 1-7$, probing a regime that is not well tested observationally. We explore predictions against other surveys and experiment configurations, finding that wide patches with maximal sky overlap between CMB and LSS surveys are most powerful for $\sigma_8(z)$ and $f_\mathrm{NL}$.

We investigate the impact of non-Gaussian lensing deflections on measurements of the CMB lensing power spectrum. We find that the false assumption of their Gaussianity significantly biases these measurements in current and future experiments at the percent level. The bias is detected by comparing CMB lensing reconstructions from simulated CMB data lensed with Gaussian deflection fields to reconstructions from simulations lensed with fully non-Gaussian deflection fields. The non-Gaussian deflections are produced by ray-tracing through snapshots of an N-body simulation and capture both the non-Gaussianity induced by non-linear structure formation and by multiple correlated deflections. We find that the amplitude of the measured bias is in agreement with analytical predictions by B\"ohm et al. 2016. The bias is largest in temperature-based measurements and we do not find evidence for it in measurements from a combination of polarization fields ($EB,EB$). We argue that the non-Gaussian bias should be even more important for measurements of cross-correlations of CMB lensing with low-redshift tracers of large-scale structure.

Future Baryon Acoustic Oscillation surveys aim at observing galaxy clustering over a wide range of redshift and galaxy populations at great precision, reaching tenths of a percent, in order to detect any deviation of dark energy from the $\LCDM$ model. We utilize a set of paired quasi-\Nb\, FastPM simulations that were designed to mitigate the sample variance effect on the BAO feature and evaluated the BAO systematics as precisely as $\sim 0.01\%$. We report anisotropic BAO scale shifts before and after density field reconstruction in the presence of redshift-space distortions over a wide range of redshift, galaxy/halo biases, and shot noise levels. We test different reconstruction schemes and different smoothing filter scales, and introduce physically-motivated BAO fitting models. For the first time, we derive a Galilean-invariant infrared resummed model for halos in real and redshift space. We test these models from the perspective of robust BAO measurements and non-BAO information such as growth rate and nonlinear bias. We find that pre-reconstruction BAO scale has moderate fitting-model dependence at the level of $0.1\%-0.2\%$ for matter while the dependence is substantially reduced to less than $0.07\%$ for halos. We find that post-reconstruction BAO shifts are generally reduced to below $0.1\%$ in the presence of galaxy/halo bias and show much smaller fitting model dependence. Different reconstruction conventions can potentially make a much larger difference on the line-of-sight BAO scale, upto $0.3\%$. Meanwhile, the precision (error) of the BAO measurements is quite consistent regardless of the choice of the fitting model or reconstruction convention.

As a major source of cosmological information, galaxy clustering is susceptible to long-wavelength density and tidal fluctuations. These long modes modulate the growth and expansion rate of local structures, shifting them in both amplitude and scale. These effects are often named the growth and dilation effects, respectively. In particular the dilation shifts the baryon acoustic oscillation (BAO) peak and breaks the assumption of the Alcock-Paczynski (AP) test. This cannot be removed with reconstruction techniques because the effect originates from long modes outside the survey. In redshift space, the long modes generate a large-scale radial peculiar velocity that affects the redshift-space distortion (RSD) signal. We compute the redshift-space response functions of the galaxy power spectrum to long density and tidal modes at leading order in perturbation theory, including both the growth and dilation terms. We validate these response functions against measurements from simulated galaxy mock catalogs. As one application, long density and tidal modes beyond the scale of a survey correlate various observables leading to an excess error known as the super-sample covariance, and thus weaken their constraining power. We quantify the super-sample effect on BAO, AP, and RSD measurements, and study its impact on current and future surveys.

Motivated by recent developments in perturbative calculations of the nonlinear evolution of large-scale structure, we present an iterative algorithm to reconstruct the initial conditions in a given volume starting from the dark matter distribution in real space. In our algorithm, objects are first moved back iteratively along estimated potential gradients, with a progressively reduced smoothing scale, until a nearly uniform catalog is obtained. The linear initial density is then estimated as the divergence of the cumulative displacement, with an optional second-order correction. This algorithm should undo nonlinear effects up to one-loop order, including the higher-order infrared resummation piece. We test the method using dark matter simulations in real space. At redshift $z=0$, we find that after eight iterations the reconstructed density is more than $95\%$ correlated with the initial density at $k\le 0.35\; h\mathrm{Mpc}^{-1}$. The reconstruction also reduces the power in the difference between reconstructed and initial fields by more than 2 orders of magnitude at $k\le 0.2\; h\mathrm{Mpc}^{-1}$, and it extends the range of scales where the full broadband shape of the power spectrum matches linear theory by a factor of 2-3. As a specific application, we consider measurements of the baryonic acoustic oscillation (BAO) scale that can be improved by reducing the degradation effects of large-scale flows. In our idealized dark matter simulations, the method improves the BAO signal-to-noise ratio by a factor of 2.7 at $z=0$ and by a factor of 2.5 at $z=0.6$, improving standard BAO reconstruction by $70\%$ at $z=0$ and $30\%$ at $z=0.6$, and matching the optimal BAO signal and signal-to-noise ratio of the linear density in the same volume. For BAO, the iterative nature of the reconstruction is the most important aspect.

We study the matter bispectrum of large-scale structure by comparing the predictions of different perturbative and phenomenological models with the full three-dimensional bispectrum from $N$-body simulations estimated using modal methods. We show that among the perturbative approaches, effective field theory succeeds in extending the range of validity furthest on intermediate scales, at the cost of free additional parameters. By studying the halo model, we show that although it is satisfactory in the deeply non-linear regime, it predicts a deficit of power on intermediate scales, worsening at redshifts $z>0$. By comparison with the $N$-body bispectrum on those scales, we show that there is a significant squeezed component underestimated in the halo model. On the basis of these results, we propose a new three-shape model, based on the tree-level, squeezed and constant bispectrum shapes we identified in the halo model; after calibration this fits the simulations on all scales and redshifts of interest. We extend this model further to primordial non-Gaussianity of the local and equilateral types by showing that the same shapes can be used to describe the additional non-Gaussian component in the matter bispectrum. This method provides a HALOFIT-like prototype of the bispectrum that could be used to describe and test parameter dependencies and should be relevant for the bispectrum of weak gravitational lensing and wider applications.

CMB and lensing reconstruction power spectra are powerful probes of cosmology. However they are correlated, since the CMB power spectra are lensed and the lensing reconstruction is constructed using CMB multipoles. We perform a full analysis of the auto- and cross-covariances, including polarization power spectra and minimum variance lensing estimators, and compare with simulations of idealized future CMB-S4 observations. Covariances sourced by fluctuations in the unlensed CMB and instrumental noise can largely be removed by using a realization-dependent subtraction of lensing reconstruction noise, leaving a relatively simple covariance model that is dominated by lensing-induced terms and well described by a small number of principal components. The correlations between the CMB and lensing power spectra will be detectable at the level of $\sim 5\sigma$ for a CMB-S4 mission, and neglecting those could underestimate some parameter error bars by several tens of percent. However we found that the inclusion of external priors or data sets to estimate parameter error bars can make the impact of the correlations almost negligible.

Modeling the large-scale structure of the universe on nonlinear scales has the potential to substantially increase the science return of upcoming surveys by increasing the number of modes available for model comparisons. One way to achieve this is to model nonlinear scales perturbatively. Unfortunately, this involves high-dimensional loop integrals that are cumbersome to evaluate. Trying to simplify this, we show how two-loop (next-to-next-to-leading order) corrections to the density power spectrum can be reduced to low-dimensional, radial integrals. Many of those can be evaluated with a one-dimensional Fast Fourier Transform, which is significantly faster than the five-dimensional Monte-Carlo integrals that are needed otherwise. The general idea of this FFT-PT method is to switch between Fourier and position space to avoid convolutions and integrate over orientations, leaving only radial integrals. This reformulation is independent of the underlying shape of the initial linear density power spectrum and should easily accommodate features such as those from baryonic acoustic oscillations. We also discuss how to account for halo bias and redshift space distortions.

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.

Forming a three dimensional view of the Universe is a long-standing goal of astronomical observations, and one that becomes increasingly difficult at high redshift. In this paper we discuss how tomography of the intergalactic medium (IGM) at $z\simeq 2.5$ can be used to estimate the redshifts of massive galaxies in a large volume of the Universe based on spectra of galaxies in their background. Our method is based on the fact that hierarchical structure formation leads to a strong dependence of the halo density on large-scale environment. A map of the latter can thus be used to refine our knowledge of the redshifts of halos and the galaxies and AGN which they host. We show that tomographic maps of the IGM at a resolution of $2.5\,h^{-1}$Mpc can determine the redshifts of more than 90 per cent of massive galaxies with redshift uncertainty $\Delta z/(1+z)=0.01$. Higher resolution maps allow such redshift estimation for lower mass galaxies and halos.

Given the importance of future large scale structure surveys for delivering new cosmological information, it is crucial to reliably predict their observables. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a manifestly convergent perturbative scheme to compute the clustering of dark matter in the weakly nonlinear regime in an expansion in $k/k_{\rm NL}$, where $k$ is the wavenumber of interest and $k_{\rm NL}$ is the wavenumber associated to the nonlinear scale. It has been recently shown that the EFTofLSS matches to $1\%$ level the dark matter power spectrum at redshift zero up to $k\simeq 0.3 h\,$Mpc$^{-1}$ and $k\simeq 0.6 h\,$Mpc$^{-1}$ at one and two loops respectively, using only one counterterm that is fit to data. Similar results have been obtained for the momentum power spectrum at one loop. This is a remarkable improvement with respect to former analytical techniques. Here we study the prediction for the equal-time dark matter bispectrum at one loop. We find that at this order it is sufficient to consider the same counterterm that was measured in the power spectrum. Without any remaining free parameter, and in a cosmology for which $k_{\rm NL}$ is smaller than in the previously considered cases ($\sigma_8=0.9$), we find that the prediction from the EFTofLSS agrees very well with $N$-body simulations up to $k\simeq 0.25 h\,$Mpc$^{-1}$, given the accuracy of the measurements, which is of order a few percent at the highest $k$'s of interest. While the fit is very good on average up to $k\simeq 0.25 h\,$Mpc$^{-1}$, the fit performs slightly worse on equilateral configurations, in agreement with expectations that for a given maximum $k$, equilateral triangles are the most nonlinear.

The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional Fast Fourier Transforms. This is exact and a few orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. The method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.

The rapidly improving precision of measurements of gravitational lensing of the Cosmic Microwave Background (CMB) also requires a corresponding increase in the precision of theoretical modeling. A commonly made approximation is to model the CMB deflection angle or lensing potential as a Gaussian random field. In this paper, however, we analytically quantify the influence of the non-Gaussianity of large-scale structure lenses, arising from nonlinear structure formation, on CMB lensing measurements. In particular, evaluating the impact of the non-zero bispectrum of large-scale structure on the relevant CMB four-point correlation functions, we find that there is a bias to estimates of the CMB lensing power spectrum. For temperature-based lensing reconstruction with CMB Stage-III and Stage-IV experiments, we find that this lensing power spectrum bias is negative and is of order one percent of the signal. This corresponds to a shift of multiple standard deviations for these upcoming experiments. We caution, however, that our numerical calculation only evaluates two of the largest bias terms and thus only provides an approximate estimate of the full bias. We conclude that further investigation into lensing biases from nonlinear structure formation is required and that these biases should be accounted for in future lensing analyses.

We study the matter bispectrum of the large-scale structure by comparing different perturbative and phenomenological models with measurements from $N$-body simulations obtained with a modal bispectrum estimator. Using shape and amplitude correlators, we directly compare simulated data with theoretical models over the full three-dimensional domain of the bispectrum, for different redshifts and scales. We review and investigate the main perturbative methods in the literature that predict the one-loop bispectrum: standard perturbation theory, effective field theory, resummed Lagrangian and renormalised perturbation theory, calculating the latter also at two loops for some triangle configurations. We find that effective field theory (EFT) succeeds in extending the range of validity furthest into the mildly nonlinear regime, albeit at the price of free extra parameters requiring calibration on simulations. For the more phenomenological halo model, we confirm that despite its validity in the deeply nonlinear regime it has a deficit of power on intermediate scales, which worsens at higher redshifts; this issue is ameliorated, but not solved, by combined halo-perturbative models. We show from simulations that in this transition region there is a strong squeezed bispectrum component that is significantly underestimated in the halo model at earlier redshifts. We thus propose a phenomenological method for alleviating this deficit, which we develop into a simple phenomenological "three-shape" benchmark model based on the three fundamental shapes we have obtained from studying the halo model. When calibrated on the simulations, this three-shape benchmark model accurately describes the bispectrum on all scales and redshifts considered, providing a prototype bispectrum HALOFIT-like methodology that could be used to describe and test parameter dependencies.

As galaxy surveys begin to measure the imprint of baryonic acoustic oscillations (BAO) on large-scale structure at the sub-percent level, reconstruction techniques that reduce the contamination from nonlinear clustering become increasingly important. Inverting the nonlinear continuity equation, we propose an Eulerian growth-shift reconstruction algorithm that does not require the displacement of any objects, which is needed for the standard Lagrangian BAO reconstruction algorithm. In real-space DM-only simulations the algorithm yields 95% of the BAO signal-to-noise obtained from standard reconstruction. The reconstructed power spectrum is obtained by adding specific simple 3- and 4-point statistics to the pre-reconstruction power spectrum, making it very transparent how additional BAO information from higher-point statistics is included in the power spectrum through the reconstruction process. Analytical models of the reconstructed density for the two algorithms agree at second order. Based on similar modeling efforts, we introduce four additional reconstruction algorithms and discuss their performance.

Clustering of large-scale structure provides significant cosmological information through the power spectrum of density perturbations. Additional information can be gained from higher-order statistics like the bispectrum, especially to break the degeneracy between the linear halo bias $b_1$ and the amplitude of fluctuations $\sigma_8$. We propose new simple, computationally inexpensive bispectrum statistics that are near optimal for the specific applications like bias determination. Corresponding to the Legendre decomposition of nonlinear halo bias and gravitational coupling at second order, these statistics are given by the cross-spectra of the density with three quadratic fields: the squared density, a tidal term, and a shift term. For halos and galaxies the first two have associated nonlinear bias terms $b_2$ and $b_{s^2}$, respectively, while the shift term has none in the absence of velocity bias (valid in the $k \rightarrow 0$ limit). Thus the linear bias $b_1$ is best determined by the shift cross-spectrum, while the squared density and tidal cross-spectra mostly tighten constraints on $b_2$ and $b_{s^2}$ once $b_1$ is known. Since the form of the cross-spectra is derived from optimal maximum-likelihood estimation, they contain the full bispectrum information on bias parameters. Perturbative analytical predictions for their expectation values and covariances agree with simulations on large scales, $k\lesssim 0.09h/\mathrm{Mpc}$ at $z=0.55$ with Gaussian $R=20h^{-1}\mathrm{Mpc}$ smoothing, for matter-matter-matter, and matter-matter-halo combinations. For halo-halo-halo cross-spectra the model also needs to include corrections to the Poisson stochasticity.

As confusion with lensing B-modes begins to limit experiments that search for primordial B-mode polarization, robust methods for delensing the CMB polarization sky are becoming increasingly important. We investigate in detail the possibility of delensing the CMB with the cosmic infrared background (CIB), emission from dusty star-forming galaxies that is an excellent tracer of the CMB lensing signal, in order to improve constraints on the tensor-to-scalar ratio $r$. We find that the maps of the CIB, such as current Planck satellite maps at 545 GHz, can be used to remove more than half of the lensing B-mode power. Calculating optimal combinations of different large-scale-structure tracers for delensing, we find that co-adding CIB data and external arcminute-resolution CMB lensing reconstruction can lead to significant additional improvements in delensing performance. We investigate whether measurement uncertainty in the CIB spectra will degrade the delensing performance if no model of the CIB spectra is assumed, and instead the CIB spectra are marginalized over, when constraining $r$. We find that such uncertainty does not significantly affect B-mode surveys smaller than a few thousand degrees. Even for larger surveys it causes only a moderate reduction in CIB delensing performance, especially if the surveys have high (arcminute) resolution, which allows self-calibration of the delensing procedure. Though further work on the impact of foreground residuals is required, our overall conclusions for delensing with current CIB data are optimistic: this delensing method can tighten constraints on $r$ by a factor up to $\approx2.2$, and by a factor up to $\approx4$ when combined with external $\approx 3 \mu$K-arcmin lensing reconstruction, without requiring the modeling of CIB properties. CIB delensing is thus a promising method for the upcoming generation of CMB polarization surveys.

Gravitational lensing provides a significant source of cosmological information in modern CMB parameter analyses. It is measured in both the power spectrum and trispectrum of the temperature fluctuations. These observables are often treated as independent, although as they are both determined from the same map this is impossible. In this paper, we perform a rigorous analysis of the covariance between lensing power spectrum and trispectrum analyses. We find two dominant contributions coming from: (i) correlations between the disconnected noise bias in the trispectrum measurement and sample variance in the temperature power spectrum; and (ii) sample variance of the lenses themselves. The former is naturally removed when the dominant N0 Gaussian bias in the reconstructed deflection spectrum is dealt with via a partially data-dependent correction, as advocated elsewhere for other reasons. The remaining lens-cosmic-variance contribution is easily modeled but can safely be ignored for a Planck-like experiment, justifying treating the two observable spectra as independent. We also test simple likelihood approximations for the deflection power spectrum, finding that a Gaussian with a parameter-independent covariance performs well.

We present the implementation of a fast estimator for the full dark matter bispectrum of a three-dimensional particle distribution relying on a separable modal expansion of the bispectrum. The computational cost of accurate bispectrum estimation is negligible relative to simulation evolution, so the isotropic bispectrum can be used as a standard diagnostic whenever the power spectrum is evaluated. As an application we measure the evolution of gravitational and primordial dark matter bispectra in $N$-body simulations with Gaussian and non-Gaussian initial conditions of the local, equilateral, orthogonal and flattened shape. The results are compared to theoretical models using a 3D visualisation, 3D shape correlations and the cumulative bispectrum signal-to-noise, all of which can be evaluated extremely quickly. Our measured bispectra are determined by $\mathcal{O}(50)$ coefficients, which can be used as fitting formulae in the nonlinear regime and for non-Gaussian initial conditions. In the nonlinear regime with $k<2h\,\mathrm{Mpc}^{-1}$, we find an excellent correlation between the measured dark matter bispectrum and a simple model based on a `constant' bispectrum plus a (nonlinear) tree-level gravitational bispectrum. In the same range for non-Gaussian simulations, we find an excellent correlation between the measured additional bispectrum and a constant model plus a (nonlinear) tree-level primordial bispectrum. We demonstrate that the constant contribution to the non-Gaussian bispectrum can be understood as a time-shift of the constant mode in the gravitational bispectrum, which is motivated by the one-halo model. The final amplitude of this extra non-Gaussian constant contribution is directly related to the initial amplitude of the constant mode in the primordial bispectrum. We also comment on the effects of regular grid and glass initial conditions on the bispectrum.

In this paper we present the implementation of an efficient formalism for the generation of arbitrary non-Gaussian initial conditions for use in N-body simulations. The methodology involves the use of a separable modal approach for decomposing a primordial bispectrum or trispectrum. This approach allows for the far more efficient generation of the non-Gaussian initial conditions already described in the literature, as well as the generation for the first time of non-separable bispectra and the special class of diagonal-free trispectra. The modal approach also allows for the reconstruction of the spectra from given realisations, a fact which is exploited to provide an accurate consistency check of the simulations.