We construct a model for the angular power spectrum of the instrumental noise in interferometer networks mapping gravitational wave backgrounds (GWBs) as a function of detector noise properties, network configuration and scan strategy. We use the model to calculate the noise power spectrum for current and future ground-based experiments, as well as for planned space missions. We present our results in a language similar to that used in cosmic microwave background and intensity mapping experiments, and connect the formalism with the sensitivity curves that are common lore in GWB analyses. Our formalism is implemented in a lightweight python module that we make publicly available at https://github.com/damonge/schNell.

We propose an observable for ultra-high energy cosmic ray (UHECR) physics: the harmonic-space cross-correlation power spectrum between the arrival directions of UHECRs and the large-scale cosmic structure mapped by galaxies. This cross-correlation has not yet been considered in the literature. We describe the expected form of the cross-correlation and show how, if the distribution of UHECR sources trace the large-scale cosmic structure, it could be easier to detect with current data than the UHECR auto-correlation. Moreover, the cross-correlation is more sensitive to UHECR anisotropies on smaller angular scales, more robust to systematic uncertainties, and it could be used to determine the redshift distribution of UHECR sources, as well as the UHECR chemical composition, making it a valuable tool in determining their origins and properties.

One of the most important checks for systematic errors in CMB studies is the cross correlation of maps made by independent experiments. In this paper we report on the cross correlation between maps from the Atacama B-mode Search (ABS) and Atacama Cosmology Telescope (ACT) experiments in both temperature and polarization. These completely different measurements have a clear correlation with each other and with the Planck satellite in both the EE and TE spectra at $\ell<400$ over the roughly $1100$ deg$^2$ common to all three. The TB, EB, and BB cross spectra are consistent with noise. Exploiting such cross-correlations will be important for future experiments operating in Chile that aim to probe the $30<\ell<8,000$ range.

The spatial and temporal discreteness of gravitational wave sources leads to shot noise that may, in some regimes, swamp any attempts at measuring the anisotropy of the gravitational wave background. Cross-correlating a gravitational wave background map with a sufficiently dense galaxy survey can alleviate this issue, and potentially recover some of the underlying properties of the gravitational wave background. We quantify the shot noise level and we explicitly show that cross-correlating the gravitational wave background and a galaxy catalog improves the chances of a first detection of the background anisotropy with a gravitational wave observatory operating in the frequency range (10Hz,100Hz), given sufficient sensitivity.

We cross-correlate maps of the thermal Sunyaev-Zeldovich (tSZ) Compton-$y$ parameter published by Planck with the projected distribution of galaxies in a set of low-redshift tomographic bins. We use the nearly full-sky 2MASS Photometric Redshift and WISE $\times$ SuperCOSMOS public catalogues, covering the redshift range $z\lesssim0.4$. Our measurements allow us to place constraints on the redshift dependence of the mass-observable relation for tSZ cluster count analyses in terms of the so-called 'hydrostatic mass bias' parameter $1-b_{\rm H}$. These results can also be interpreted as measurements of the bias-weighted average gas pressure $\langle bP_e\rangle$ as a function of redshift, a quantity that can be related to the thermodynamics of gas inside haloes and used to constrain energy injection processes. We measure $1-b_{\rm H}$ with $\sim13\%$ precision in 6 equispaced redshift bins, and find no evidence for a redshift-dependent mass bias parameter, in agreement with previous analyses. Our mean value of $1-b_{\rm H} = 0.59\pm0.03$ is also in good agreement with the one estimated by the joint analysis of Planck cluster counts and CMB anisotropies. Our measurements of $\langle bP_e\rangle$, at the level of $\sim10\%$ in each bin, are the most stringent constraints on the redshift dependence of this parameter to date, and agree well both with previous measurements and with theoretical expectations from shock-heating models.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

We analyze the clustering of galaxies in the first public data release of the HSC Subaru Strategic Program. Despite the relatively small footprints of the observed fields, the data are an excellent proxy for the deep photometric datasets that will be acquired by LSST, and are therefore an ideal test bed for the analysis methods being implemented by the LSST DESC. We select a magnitude limited sample with $i<24.5$ and analyze it in four redshift bins covering $0.15\lesssim z \lesssim1.5$. We carry out a Fourier-space analysis of the two-point clustering of this sample, including all auto- and cross-correlations. We demonstrate the use of map-level deprojection methods to account for fluctuations in the galaxy number density caused by observational systematics. Through an HOD analysis, we place constraints on the characteristic halo masses of this sample, finding a good fit up to scales $k_{\rm max}=1\,{\rm Mpc}^{-1}$, including both auto- and cross-correlations. Our results show monotonically decreasing average halo masses, which can be interpreted in terms of the drop-out of red galaxies at high redshifts for a flux-limited sample. In terms of photometric redshift systematics, we show that additional care is needed in order to marginalize over uncertainties in the redshift distribution in galaxy clustering, and that these uncertainties can be constrained by including cross-correlations. We are able to make a $\sim3\sigma$ detection of lensing magnification in the HSC data. Our results are stable to variations in $\sigma_8$ and $\Omega_c$ and we find constraints that agree well with measurements from Planck and low-redshift probes. Finally, we use our pipeline to study the clustering of galaxies as a function of limiting flux, and provide a simple fitting function for the linear galaxy bias for magnitude limited samples as a function of limiting magnitude and redshift. [abridged]

We investigate the sensitivity to the effects of lensing magnification on large-scale structure analyses combining photometric cosmic shear and galaxy clustering data (i.e. the now commonly called "3$\times$2-point" analysis). Using a Fisher matrix bias formalism, we disentangle the contribution to the bias on cosmological parameters caused by ignoring the effects of magnification in a theory fit from individual elements in the data vector, for Stage-III and Stage-IV surveys, assuming well-known redshift distributions, sample selection based on magnitude, and magnification strengths inferred from CFHTLenS and DES Y1 data. We show that the removal of elements of the data vectors that are dominated by magnification does not guarantee a reduction in the cosmological bias due to the magnification signal, but can instead increase the sensitivity to magnification. We find that the most sensitive elements of the data vector come from the shear-clustering cross-correlations, particularly between the highest-redshift shear bin and any lower-redshift lens sample, and that the parameters $\Omega_{M}$, $S_8=\sigma_8\sqrt{\Omega_M/0.3}$ and $w_0$ show the most significant biases for both survey models. Our forecasts predict that current analyses are not significantly biased by magnification, but this bias will become highly significant with the continued increase of statistical power in the near future. We therefore conclude that future surveys should measure and model the magnification as part of their flagship "3$\times$2-point" analysis.

3D data compression techniques can be used to determine the natural basis of radial eigenmodes that encode the maximum amount of information in a tomographic large-scale structure survey. We explore the potential of the Karhunen-Lo\`eve decomposition in reducing the dimensionality of the data vector for cosmic shear measurements, and apply it to the final data from the \cfh survey. We find that practically all of the cosmological information can be encoded in one single radial eigenmode, from which we are able to reproduce compatible constraints with those found in the fiducial tomographic analysis (done with 7 redshift bins) with a factor of ~30 fewer datapoints. This simplifies the problem of computing the two-point function covariance matrix from mock catalogues by the same factor, or by a factor of ~800 for an analytical covariance. The resulting set of radial eigenfunctions is close to ell-independent, and therefore they can be used as redshift-dependent galaxy weights. This simplifies the application of the Karhunen-Lo\`eve decomposition to real-space and Fourier-space data, and allows one to explore the effective radial window function of the principal eigenmodes as well as the associated shear maps in order to identify potential systematics. We also apply the method to extended parameter spaces and verify that additional information may be gained by including a second mode to break parameter degeneracies. The data and analysis code are publicly available at https://github.com/emiliobellini/kl_sample.

The statistical power of Lyman-${\alpha}$ forest Baryon Acoustic Oscillation (BAO) measurements is set to increase significantly in the coming years as new instruments such as the Dark Energy Spectroscopic Instrument deliver progressively more constraining data. Generating mock datasets for such measurements will be important for validating analysis pipelines and evaluating the effects of systematics. With such studies in mind, we present LyaCoLoRe: a package for producing synthetic Lyman-${\alpha}$ forest survey datasets for BAO analyses. LyaCoLoRe transforms initial Gaussian random field skewers into skewers of transmitted flux fraction via a number of fast approximations. In this work we explain the methods of producing mock datasets used in LyaCoLoRe, and then measure correlation functions on a suite of realisations of such data. We demonstrate that we are able to recover the correct BAO signal, as well as large-scale bias parameters similar to literature values. Finally, we briefly describe methods to add further astrophysical effects to our skewers - high column density systems and metal absorbers - which act as potential complications for BAO analyses.

The disconnected part of the power spectrum covariance matrix (also known as the "Gaussian" covariance) is the dominant contribution on large scales for galaxy clustering and weak lensing datasets. The presence of a complicated sky mask causes non-trivial correlations between different Fourier/harmonic modes, which must be accurately characterized in order to obtain reliable cosmological constraints. This is particularly relevant for galaxy survey data. Unfortunately, an exact calculation of these correlations involves $O(\ell_{\rm max}^6)$ operations that become computationally impractical very quickly. We present an implementation of approximate methods to estimate the Gaussian covariance matrix of power spectra involving spin-0 and spin-2 flat- and curved-sky fields, expanding on existing algorithms. These methods achieve an $O(\ell_{\rm max}^3)$ scaling, which makes the computation of the covariance matrix as fast as the computation of the power spectrum itself. We quantify the accuracy of these methods on large-scale structure and weak lensing data, making use of a large number of Gaussian but otherwise realistic simulations. We show that, using the approximate covariance matrix, we are able to recover the true posterior distribution of cosmological parameters to high accuracy. We also quantify the shortcomings of these methods, which become unreliable on the very largest scales, as well as for covariance matrix elements involving cosmic shear $B$ modes. The algorithms presented here are implemented in the public code NaMaster (https://github.com/LSSTDESC/NaMaster).

This Science White Paper, prepared in response to the ESA Voyage 2050 call for long-term mission planning, aims to describe the various science possibilities that can be realized with an L-class space observatory that is dedicated to the study of the interactions of cosmic microwave background (CMB) photons with the cosmic web. Our aim is specifically to use the CMB as a backlight -- and survey the gas, total mass, and stellar content of the entire observable Universe by means of analyzing the spatial and spectral distortions imprinted on it. These distortions result from two major processes that impact on CMB photons: scattering by electrons (Sunyaev-Zeldovich effect in diverse forms, Rayleigh scattering, resonant scattering) and deflection by gravitational potential (lensing effect). Even though the list of topics collected in this White Paper is not exhaustive, it helps to illustrate the exceptional diversity of major scientific questions that can be addressed by a space mission that will reach an angular resolution of 1.5 arcmin (goal 1 arcmin), have an average sensitivity better than 1 uK-arcmin, and span the microwave frequency range from roughly 50 GHz to 1 THz. The current paper also highlights the synergy of our BACKLIGHT mission concept with several upcoming and proposed ground-based CMB experiments.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

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.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

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.

The $\alpha$-attractor inflationary models are nowadays favored by CMB Planck observations. Their similarity with canonical quintessence models motivates the exploration of a common framework that explains both inflation and dark energy. We study the expected constraints that next-generation cosmological experiments will be able to impose for the dark energy $\alpha$-attractor model. We systematically account for the constraining power of SNIa from WFIRST, BAO from DESI and WFIRST, galaxy clustering and shear from LSST and Stage-4 CMB experiments. We assume a tensor-to-scalar ratio, $10^{-3} < r < 10^{-2}$, which permits to explore the wide regime sufficiently close, but distinct, to a cosmological constant, without need of fine tunning the initial value of the field. We find that the combination S4CMB + LSST + SNIa will achieve the best results, improving the FoM by almost an order of magnitude; respect to the S4CMB + BAO + SNIa case. We find this is also true for the FoM of the $w_0 - w_a$ parameters. Therefore, future surveys will be uniquely able to probe models connecting early and late cosmic acceleration.

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the expected statistical uncertainty of upcoming experiments at separations under 10 Mpc. Successful extraction of information in such a regime thus requires, at the very least, unbiased models for the impact of galaxy formation on the matter distribution, and can benefit from complementary observational priors. This work reviews the current state of the art in the modelling of baryons for cosmology, from numerical methods to approximate analytical prescriptions, and makes recommendations for studies in the next decade, including a discussion of potential probe combinations that can help constrain the role of baryons in cosmological studies. We focus, in particular, on the modelling of the matter power spectrum, $P(k,z)$, as a function of scale and redshift, and of the observables derived from this quantity. This work is the result of a workshop held at the University of Oxford in November of 2018.

Upcoming observations from the Simons Observatory have been projected to constrain the tensor-to-scalar ratio, $r$, at the level of $\sigma(r)=$0.003. Here we describe one of the forecasting algorithms for the Simons Observatory in more detail, based on cleaning CMB polarization maps using a parametric model. We present a new code to perform this end-to-end forecast, and explore the assumptions in greater detail. If spatial uniformity of the spectral energy distribution of synchrotron radiation and thermal dust emission is assumed over the region planned for observations, covering almost a fifth of the sky, a bias of order 1--3$\sigma$ in $r$ is projected for foreground models consistent with current data. We find that by masking the most contaminated regions of sky, or by adopting more parameters to describe the spatial variation in spectral index for synchrotron and dust, such a bias can be mitigated for the foreground models we consider. We also explore strategies for testing whether the cleaned CMB polarization maps contain residual foreground contamination, including cross-correlating with maps tracing the foregrounds. This method also has applications for other CMB polarization experiments.

The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.

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.

Cosmological large-scale structure analyses based on two-point correlation functions often assume a Gaussian likelihood function with a fixed covariance matrix. We study the impact on cosmological parameter estimation of ignoring the parameter dependence of this covariance matrix, focusing on the particular case of joint weak-lensing and galaxy clustering analyses. Using a Fisher matrix formalism (calibrated against exact likelihood evaluation in particular simple cases), we quantify the effect of using a parameter dependent covariance matrix on both the bias and variance of the parameters. We confirm that the approximation of a parameter-independent covariance matrix is exceptionally good in all realistic scenarios. The information content in the covariance matrix (in comparison with the two point functions themselves) does not change with the fractional sky coverage. Therefore the increase in information due to the parameter dependent covariance matrix becomes negligible as the number of modes increases. Even for surveys covering less than $1\%$ of the sky, this effect only causes a bias of up to ${\cal O}(10\%)$ of the statistical uncertainties, with a misestimation of the parameter uncertainties at the same level or lower. The effect will only be smaller with future large-area surveys. Thus for most analyses the effect of a parameter-dependent covariance matrix can be ignored both in terms of the accuracy and precision of the recovered cosmological constraints.

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.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

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.

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.

The pseudo-$C_\ell$ is an algorithm for estimating the angular power and cross-power spectra that is very fast and, in realistic cases, also nearly optimal. The algorithm can be extended to deal with contaminant deprojection and $E/B$ purification, and can therefore be applied in a wide variety of scenarios of interest for current and future cosmological observations. This paper presents NaMaster, a public, validated, accurate and easy-to-use software package that, for the first time, provides a unified framework to compute angular cross-power spectra of any pair of spin-0 or spin-2 fields, contaminated by an arbitrary number of linear systematics and requiring $B$- or $E$-mode purification, both on the sphere or in the flat-sky approximation. We describe the mathematical background of the estimator, including all the features above, and its software implementation in NaMaster. We construct a validation suite that aims to resemble the types of observations that next-generation large-scale structure and ground-based CMB experiments will face, and use it to show that the code is able to recover the input power spectra in the most complex scenarios with no detectable bias. NaMaster can be found at https://github.com/LSSTDESC/NaMaster, and is provided with comprehensive documentation and a number of code examples.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

The Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC) will use five cosmological probes: galaxy clusters, large scale structure, supernovae, strong lensing, and weak lensing. This Science Requirements Document (SRD) quantifies the expected dark energy constraining power of these probes individually and together, with conservative assumptions about analysis methodology and follow-up observational resources based on our current understanding and the expected evolution within the field in the coming years. We then define requirements on analysis pipelines that will enable us to achieve our goal of carrying out a dark energy analysis consistent with the Dark Energy Task Force definition of a Stage IV dark energy experiment. This is achieved through a forecasting process that incorporates the flowdown to detailed requirements on multiple sources of systematic uncertainty. Future versions of this document will include evolution in our software capabilities and analysis plans along with updates to the LSST survey strategy.

We use full sky simulations, including the effects of foreground contamination and removal, to explore multi-tracer synergies between a SKA-like 21cm intensity mapping survey and a LSST-like photometric galaxy redshift survey. In particular we study ratios of auto and cross-correlations between the two tracers as estimators of the ratio of their biases, a quantity that should benefit considerably from the cosmic variance cancellation of the multi-tracer approach. We show how well we should be able to measure the bias ratio on very large scales (down to $\ell \sim 3$), which is crucial to measure primordial non-Gaussianity and general relativistic effects on large scale structure. We find that, in the absence of foregrounds but with realistic noise levels of such surveys, the multi-tracer estimators are able to improve on the sensitivity of a cosmic-variance contaminated measurement by a factor of $2-4$. When foregrounds are included, estimators using the 21cm auto-correlation become biased. However, we show that cross-correlation estimators are immune to this and do not incur in any significant penalty in terms of sensitivity from discarding the auto-correlation data. However, the loss of long-wavelength radial modes caused by foreground removal in combination with the low redshift resolution of photometric surveys, reduces the sensitivity of the multi-tracer estimator, albeit still better than the cosmic variance contaminated scenario even in the noise free case. Finally we explore different alternative avenues to avoid this problem.

Cosmological measurements over the next decade will enable us to shed light on the content and evolution of the Universe. Complementary measurements of the Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations are expected to allow an indirect determination of the sum of neutrino masses, within the framework of the flat $\Lambda$CDM model. However, possible deviations from $\Lambda$CDM such as a non-zero cosmological curvature or a dark energy equation of state with $w\neq -1$ would leave similar imprints on the expansion rate of the Universe and clustering of matter. We show how future CMB measurements can be combined with late-time measurements of galaxy clustering and cosmic shear from the Large Synoptic Survey Telescope to alleviate this degeneracy. Together, they are projected to reduce the uncertainty on the neutrino mass sum to 30 meV within this more general cosmological model. Achieving a 3$\sigma$ measurement of the minimal 60 meV mass (or 4$\sigma$ assuming $w=-1$) will require a five-fold improved measurement of the optical depth to reionization, obtainable through a large-scale CMB polarization measurement.

In this work, we present transit timing variations detected for the exoplanet TrES-5b. To obtain the necessary amount of photometric data for this exoplanet, we have organized an international campaign to search for exoplanets based on the Transit Timing Variation method (TTV) and as a result of this we collected 30 new light curves, 15 light curves from the Exoplanet Transit Database (ETD) and 8 light curves from the literature for the timing analysis of the exoplanet TrES-5b. We have detected timing variations with a semi-amplitude of A ~ 0.0016 days and a period of P ~ 99 days. We carried out the N-body modeling based on the three-body problem. The detected perturbation of TrES-5b may be caused by a second exoplanet in the TrES-5 system. We have calculated the possible mass and resonance of the object: M ~ 0.24MJup at a 1:2 Resonance.

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

We combine information from the clustering of HI galaxies in the 100% data release of the Arecibo Legacy Fast ALFA survey (ALFALFA), and from the HI content of optically-selected galaxy groups found in the Sloan Digital Sky Survey (SDSS) to constrain the relation between halo mass $M_h$ and its average total HI mass content $M_{\rm HI}$. We model the abundance and clustering of neutral hydrogen through a halo-model-based approach, parametrizing the $M_{\rm HI}(M_h)$ relation as a power law with an exponential mass cutoff. To break the degeneracy between the amplitude and low-mass cutoff of the $M_{\rm HI}(M_h)$ relation, we also include a recent measurement of the cosmic HI abundance from the $\alpha$.100 sample. We find that all datasets are consistent with a power-law index $\alpha=0.44\pm 0.08$ and a cutoff halo mass $\log_{10}M_{\rm min}/(h^{-1}M_\odot)=11.27^{+0.24}_{-0.30}$. We compare these results with predictions from state-of-the-art magneto-hydrodynamical simulations, and find both to be in good qualitative agreement, although the data favours a significantly larger cutoff mass that is consistent with the higher cosmic HI abundance found in simulations. Both data and simulations seem to predict a similar value for the HI bias ($b_{\rm HI}=0.875\pm0.022$) and shot-noise power ($P_{\rm SN}=92^{+20}_{-18}\,[h^{-1}{\rm Mpc}]^3$) at redshift $z=0$.

Time-dependent and systematic variations in band gain and central frequencies of instruments used to study the Cosmic Microwave Background are important factors in the data-to-map analysis pipeline. If not properly characterized, they could limit the ability of next-generation experiments to remove astrophysical foreground contamination. Uncertainties include the instrument detector band, which could systematically change across the focal plane, as well as the calibration of the instrument used to measure the bands. A potentially major effect is time-dependent gain and band uncertainties caused by atmospheric fluctuations. More specifically, changes in atmospheric conditions lead to frequency-dependent changes in the atmospheric transmission which, in turn, leads to variations in the effective gain and central frequency of the instrument's bandpass. Using atmospheric modeling software and ACTPol bandpasses, we simulate the expected variations in band gain and central frequency for 20, 40, 90, 150, and 240 GHz bands as a function of precipitable water vapor, observing angle, and ground temperature. Combining these effects enables us to set maximum and minimum limits on the expected uncertainties in band gain and central frequency over the course of a full observing season. We then introduce the uncertainties to parametric maximum-likelihood component separation methods on simulated CMB maps to forecast foreground removal performance and likelihoods on the tensor-to-scalar ratio r. We conclude that to measure a $\sigma(r=0)$ ~ $10^{-3}$ with a bias on the recovered $r$ under control, the limit on the uncertainty in the relative gain and central frequency of the bandpass must be <2% and <1%, respectively. We also comment on the possibility of self-calibrating bandpass uncertainties.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

Future surveys will access large volumes of space and hence very long wavelength fluctuations of the matter density and gravitational field. It has been argued that the set of secondary effects that affect the galaxy distribution, relativistic in nature, will bring new, complementary cosmological constraints. We study this claim in detail by focusing on a subset of wide-area future surveys: Stage-4 cosmic microwave background experiments and photometric redshift surveys. In particular, we look at the magnification lensing contribution to galaxy clustering and general relativistic corrections to all observables. We quantify the amount of information encoded in these effects in terms of the tightening of the final cosmological constraints as well as the potential bias in inferred parameters associated with neglecting them. We do so for a wide range of cosmological parameters, covering neutrino masses, standard dark-energy parametrizations and scalar-tensor gravity theories. Our results show that, while the effect of lensing magnification to number counts does not contain a significant amount of information when galaxy clustering is combined with cosmic shear measurements, this contribution does play a significant role in biasing estimates on a host of parameter families if unaccounted for. Since the amplitude of the magnification term is controlled by the slope of the source number counts with apparent magnitude, $s(z)$, we also estimate the accuracy to which this quantity must be known to avoid systematic parameter biases, finding that future surveys will need to determine $s(z)$ to the $\sim$5-10\% level. On the contrary, large-scale general-relativistic corrections are irrelevant both in terms of information content and parameter bias for most cosmological parameters, but significant for the level of primordial non-Gaussianity.

We cross-correlate the positions of damped Lyman-$\alpha$ systems (DLAs) and their parent quasar catalog with a convergence map derived from the Planck cosmic microwave background (CMB) temperature data. We make consistent detections of the lensing signal of both samples in both Fourier and configuration space. By interpreting the excess signal present in the DLA catalog with respect to the parent quasar catalog as caused by the large scale structure traced by DLAs, we are able to infer the bias of these objects: $b_{\rm DLA}=2.6\pm0.9$. These results are consistent with previous measurements made in cross-correlation with the Lyman-$\alpha$ forest, although the current noise in the lensing data and the low number density of DLAs limits the constraining power of this measurement. We discuss the robustness of the analysis with respect to a number different systematic effects and forecast prospects of carrying out this measurement with data from future experiments.

We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

We stack maps of the thermal Sunyaev-Zel'dovich effect produced by the Planck Collaboration around the centers of cosmic voids defined by the distribution of galaxies in the CMASS sample of the Baryon Oscillation Spectroscopic Survey, scaled by the void effective radii. We report a first detection of the associated cross-correlation at the $3.4\sigma$ level: voids are under-pressured relative to the cosmic mean. We compare the measured Compton-$y$ profile around voids with a model based solely on the spatial modulation of halo abundance with environmental density. The amplitude of the detected signal is marginally lower than predicted by an overall amplitude $\alpha_v=0.67\pm0.2$. We discuss the possible interpretations of this measurement in terms of modelling uncertainties, excess pressure in low-mass halos, or non-local heating mechanisms.

We study the correlations between parameters characterizing neutrino physics and the evolution of dark energy. Using a fluid approach, we show that time-varying dark energy models exhibit degeneracies with the cosmic neutrino background over extended periods of the cosmic history, leading to a degraded estimation of the total mass and number of species of neutrinos. We investigate how to break degeneracies and combine multiple probes across cosmic time to anchor the behaviour of the two components. We use Planck CMB data and BAO measurements from the BOSS, SDSS and 6dF surveys to present current limits on the model parameters, and then forecast the future reach from the CMB Stage-4 and DESI experiments. We show that a multi-probe analysis of current data provides only marginal improvement on the determination of the individual parameters and no reduction of the correlations. Future observations will better distinguish the neutrino mass and preserve the current sensitivity to the number of species even in case of a time-varying dark energy component.

Photometric redshift surveys map the distribution of matter in the Universe through the positions and shapes of galaxies with poorly resolved measurements of their radial coordinates. While a tomographic analysis can be used to recover some of the large-scale radial modes present in the data, this approach suffers from a number of practical shortcomings, and the criteria to decide on a particular binning scheme are commonly blind to the ultimate science goals. We present a method designed to separate and compress the data into a small number of uncorrelated radial modes, circumventing some of the problems of standard tomographic analyses. The method is based on the Karhunen-Lo\`{e}ve transform (KL), and is connected to other 3D data compression bases advocated in the literature, such as the Fourier-Bessel decomposition. We apply this method to both weak lensing and galaxy clustering. In the case of galaxy clustering, we show that the resulting optimal basis is closely associated with the Fourier-Bessel basis, and that for certain observables, such as the effects of magnification bias or primordial non-Gaussianity, the bulk of the signal can be compressed into a small number of modes. In the case of weak lensing we show that the method is able to compress the vast majority of the signal-to-noise into a single mode, and that optimal cosmological constraints can be obtained considering only three uncorrelated KL eigenmodes, considerably simplifying the analysis with respect to a traditional tomographic approach.

Imaging surveys of galaxies will have a high number density and angular resolution yet a poor redshift precision. Intensity maps of neutral hydrogen (HI) will have accurate redshift resolution yet will not resolve individual sources. Using this complementarity, we show how the clustering redshifts approach, proposed for spectroscopic surveys can also be used in combination with intensity mapping observations to calibrate the redshift distribution of galaxies in an imaging survey and, as a result, reduce uncertainties in photometric redshift measurements. We show how the intensity mapping surveys to be carried out with the MeerKAT, HIRAX and SKA instruments can improve photometric redshift uncertainties to well below the requirements of DES and LSST. The effectiveness of this method as a function of instrumental parameters, foreground subtraction and other potential systematic errors is discussed in detail.

An important indicator of modified gravity is the effect of the local environment on halo properties. This paper examines the influence of the local tidal structure on the halo mass function, the halo orientation, spin and the concentration-mass relation. We generalise the excursion set formalism to produce a halo mass function conditional on large-scale structure. Our model agrees well with simulations on large scales at which the density field is linear or weakly non-linear. Beyond this, our principal result is that $f(R)$ does affect halo abundances, the halo spin parameter and the concentration-mass relationship in an environment-independent way, whereas we find no appreciable deviation from $\Lambda$CDM for the mass function with fixed environment density, nor the alignment of the orientation and spin vectors of the halo to the eigenvectors of the local cosmic web. There is a general trend for greater deviation from $\Lambda$CDM in underdense environments and for high-mass haloes, as expected from chameleon screening.

The next generation of surveys will greatly improve our knowledge of cosmological gravity. In this paper we focus on how Stage IV photometric redshift surveys, including weak lensing and multiple tracers of the matter distribution and radio experiments combined with measurements of the cosmic microwave background will lead to precision constraints on deviations from General Relativity. We use a broad subclass of Horndeski scalar-tensor theories to forecast the accuracy with which we will be able to determine these deviations and their degeneracies with other cosmological parameters. Our analysis includes relativistic effects, does not rely on the quasi-static evolution and makes conservative assumptions about the effect of screening on small scales. We define a figure of merit for cosmological tests of gravity and show how the combination of different types of surveys, probing different length scales and redshifts, can be used to pin down constraints on the gravitational physics to better than a few percent, roughly an order of magnitude better than present probes. Future cosmological experiments will be able to constrain the Brans-Dicke parameter at a level comparable to Solar System and astrophysical tests

We present forecasts for cosmological parameters from future Cosmic Microwave Background (CMB) data measured by the Stage-4 (S4) generation of ground-based experiments in combination with large-scale anisotropy data from the PIXIE satellite. We demonstrate the complementarity of the two experiments and focus on science targets that benefit from their combination. We show that a cosmic-variance-limited measurement of the optical depth to reionization provided by PIXIE, with error $\sigma(\tau)=0.002$, is vital for enabling a 5$\sigma$ detection of the sum of the neutrino masses when combined with a CMB-S4 lensing measurement, and with lower-redshift constraints on the growth of structure and the distance-redshift relation. Parameters characterizing the epoch of reionization will also be tightly constrained; PIXIE's $\tau$ constraint converts into $\sigma(\rm{z_{re}})=0.2$ for the mean time of reionization, and a kinematic Sunyaev-Zel'dovich measurement from S4 gives $\sigma(\Delta \rm{z_{re}})=0.03$ for the duration of reionization. Both PIXIE and S4 will put strong constraints on primordial tensor fluctuations, vital for testing early-universe models, and will do so at distinct angular scales. We forecast $\sigma(r)\approx 5\times10^{-4}$ for a signal with a tensor-to-scalar ratio $r=10^{-3}$, after accounting for diffuse foreground removal and de-lensing. The wide and dense frequency coverage of PIXIE results in an expected foreground-degradation factor on $r$ of only $\approx$25%. By measuring large and small scales PIXIE and S4 will together better limit the energy injection at recombination from dark matter annihilation, with $p_{\rm ann}<0.09 \times 10^{-6} \ {\rm m^3/s/Kg}$ projected at 95% confidence. Cosmological parameters measured from the damping tail with S4 will be best constrained by polarization, which has the advantage of minimal contamination from extragalactic emission.

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.

CMB Stage-4 experiments will reduce the uncertainties on the gravitational lensing potential by an order of magnitude compared to current measurements, and will also produce a Sunyaev-Zel'dovich (SZ) cluster catalog containing $\sim10^{5}$ objects, two orders of magnitudes higher than what is currently available. In this paper we propose to combine these two observables and show that it is possible to calibrate the masses of the full Stage-4 cluster catalog internally owing to the high signal to noise measurement of the CMB lensing convergence field. We find that a CMB Stage-4 experiment will constrain the hydrostatic bias parameter to sub-percent accuracy. We also show constraints on a non parametric $Y-M$ relationship which could be used to study its evolution with mass and redshift. Finally we present a joint likelihood for thermal SZ (tSZ) flux and mass measurements, and show that it could lead to a $\sim5\sigma$ detection of the lower limit on the sum of the neutrino masses in the normal hierarchy ($\sum m_{\nu}=60 \textrm{meV}$) once combined with measurements of the primordial CMB and CMB lensing power spectra.

We quantitatively investigate the possibility of detecting baryonic acoustic oscillations (BAO) using single-dish 21cm intensity mapping observations in the post-reionization era. We show that the telescope beam smears out the isotropic BAO signature and, in the case of the Square Kilometer Array (SKA) instrument, makes it undetectable at redshifts $z\gtrsim1$. We however demonstrate that the BAO peak can still be detected in the radial 21cm power spectrum and describe a method to make this type of measurements. By means of numerical simulations, containing the 21cm cosmological signal as well as the most relevant Galactic and extra-Galactic foregrounds and basic instrumental effect, we quantify the precision with which the radial BAO scale can be measured in the 21cm power spectrum. We systematically investigate the signal-to-noise and the precision of the recovered BAO signal as a function of cosmic variance, instrumental noise, angular resolution and foreground contamination. We find that the expected noise levels of SKA would degrade the final BAO errors by $\sim5\%$ with respect to the cosmic-variance limited case at low redshifts, but that the effect grows up to $\sim65\%$ at $z\sim2-3$. Furthermore, we find that the radial BAO signature is robust against foreground systematics, and that the main effect is an increase of $\sim20\%$ in the final uncertainty on the standard ruler caused by the contribution of foreground residuals as well as the reduction in sky area needed to avoid high-foreground regions. We also find that it should be possible to detect the radial BAO signature with high significance in the full redshift range. We conclude that a 21cm experiment carried out by the SKA should be able to make direct measurements of the expansion rate $H(z)$ with competitive per-cent level precision on redshifts $z\lesssim2.5$.

We present a numerical code to simulate maps of Galactic emission in intensity and polarization at microwave frequencies, aiding in the design of Cosmic Microwave Background experiments. This Python code builds on existing efforts to simulate the sky by providing an easy-to-use interface and is based on publicly available data from the WMAP and Planck satellite missions. We simulate synchrotron, thermal dust, free-free, and anomalous microwave emission over the whole sky, in addition to the Cosmic Microwave Background, and include a set of alternative prescriptions for the frequency dependence of each component that are consistent with current data. We also present a prescription for adding small-scale realizations of these components at resolutions greater than current all-sky measurements. The code is available at https://github.com/bthorne93/PySM_public.

Detecting the imprint of inflationary gravitational waves on the $B$-mode polarization of the Cosmic Microwave Background (CMB) is one of the main science cases for current and next-generation CMB experiments. In this work we explore some of the challenges that ground-based facilities will have to face in order to carry out this measurement in the presence of Galactic foregrounds and correlated atmospheric noise. We present forecasts for Stage-3 (S3) and planned Stage-4 (S4) experiments based on the analysis of simulated sky maps using a map-based Bayesian foreground cleaning method. Our results thus consistently propagate the uncertainties on foreground parameters such as spatially-varying spectral indices, as well as the bias on the measured tensor-to-scalar ratio $r$ caused by an incorrect modelling of the foregrounds. We find that S3 and S4-like experiments should be able to put constraints on $r$ of the order $\sigma(r)=(0.5-1.0)\times10^{-2}$ and $\sigma(r)=(0.5-1.0)\times10^{-3}$ respectively, assuming instrumental systematic effects are under control. We further study deviations from the fiducial foreground model, finding that, while the effects of a second polarized dust component would be minimal on both S3 and S4, a 2\% polarized anomalous dust emission (AME) component would be clearly detectable by Stage-4 experiments.

In recent years there have been significant improvements in the sensitivity and the angular resolution of the instruments dedicated to the observation of the Cosmic Microwave Background (CMB). ACTPol is the first polarization receiver for the Atacama Cosmology Telescope (ACT) and is observing the CMB sky with arcmin resolution over about 2000 sq. deg. Its upgrade, Advanced ACTPol (AdvACT), will observe the CMB in five frequency bands and over a larger area of the sky. We describe the optimization and implementation of the ACTPol and AdvACT surveys. The selection of the observed fields is driven mainly by the science goals, that is, small angular scale CMB measurements, B-mode measurements and cross-correlation studies. For the ACTPol survey we have observed patches of the southern galactic sky with low galactic foreground emissions which were also chosen to maximize the overlap with several galaxy surveys to allow unique cross-correlation studies. A wider field in the northern galactic cap ensured significant additional overlap with the BOSS spectroscopic survey. The exact shapes and footprints of the fields were optimized to achieve uniform coverage and to obtain cross-linked maps by observing the fields with different scan directions. We have maximized the efficiency of the survey by implementing a close to 24 hour observing strategy, switching between daytime and nighttime observing plans and minimizing the telescope idle time. We describe the challenges represented by the survey optimization for the significantly wider area observed by AdvACT, which will observe roughly half of the low-foreground sky. The survey strategies described here may prove useful for planning future ground-based CMB surveys, such as the Simons Observatory and CMB Stage IV surveys.

Future ground-based CMB experiments will generate competitive large-scale structure datasets by precisely characterizing CMB secondary anisotropies over a large fraction of the sky. We describe a method for constraining the growth rate of structure to sub-1% precision out to $z\approx 1$, using a combination of galaxy cluster peculiar velocities measured using the kinetic Sunyaev-Zel'dovich (kSZ) effect, and the velocity field reconstructed from galaxy redshift surveys. We consider only thermal SZ-selected cluster samples, which will consist of $\mathcal{O}(10^4-10^5)$ sources for Stage 3 and 4 CMB experiments respectively. Three different methods for separating the kSZ effect from the primary CMB are compared, including a novel blind "constrained realization" method that improves signal-to-noise by a factor of $\sim 2$ over a commonly-used aperture photometry technique. Measurements of the integrated tSZ $y$-parameter are used to break the kSZ velocity-optical depth degeneracy, and the effects of including CMB polarization and SZ profile uncertainties are also considered. A combination of future Stage 4 experiments should be able to measure the product of the growth and expansion rates, $\alpha\equiv f H$, to better than 1% in bins of $\Delta z = 0.1$ out to $z \approx 1$ -- competitive with contemporary redshift-space distortion constraints from galaxy surveys.

We present a method to recover and study the projected gravitational tidal forces from a galaxy survey containing little or no redshift information. The method and the physical interpretation of the recovered tidal maps as a tracer of the cosmic web are described in detail. We first apply the method to a simulated galaxy survey and study the accuracy with which the cosmic web can be recovered in the presence of different observational effects, showing that the projected tidal field can be estimated with reasonable precision over large regions of the sky. We then apply our method to the 2MASS survey and present a publicly available full-sky map of the projected tidal forces in the local Universe. As an example of an application of these data we further study the distribution of galaxy luminosities across the different elements of the cosmic web, finding that, while more luminous objects are found preferentially in the most dense environments, there is no further segregation by tidal environment.

We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at $z\sim[3-5]$. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ~10 times higher (and down to ~4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.

Future surveys of large-scale structure will be able to measure perturbations on the scale of the cosmological horizon, and so could potentially probe a number of novel relativistic effects that are negligibly small on sub-horizon scales. These effects leave distinctive signatures in the power spectra of clustering observables and, if measurable, would open a new window on relativistic cosmology. We quantify the size and detectability of the effects for the most relevant future large-scale structure experiments: spectroscopic and photometric galaxy redshift surveys, intensity mapping surveys of neutral hydrogen, and radio continuum surveys. Our forecasts show that next-generation experiments, reaching out to redshifts $z\simeq 4$, will not be able to detect previously-undetected general-relativistic effects by using individual tracers of the density field, although the contribution of weak lensing magnification on large scales should be clearly detectable. We also perform a rigorous joint forecast for the detection of primordial non-Gaussianity through the excess power it produces in the clustering of biased tracers on large scales, finding that uncertainties of $\sigma(f_{\rm NL})\sim 1-2$ should be achievable. We study the level of degeneracy of these large-scale effects with several tracer-dependent nuisance parameters, quantifying the minimal priors on the latter that are needed for an optimal measurement of the former. Finally, we discuss the systematic effects that must be mitigated to achieve this level of sensitivity, and some alternative approaches that should help to improve the constraints. The computational tools developed to carry out this study, which requires the full-sky computation of the theoretical angular power spectra for $\mathcal{O}(100)$ redshift bins, as well as realistic models of the luminosity function, are publicly available.

Multiple tracers of the cosmic density field, with different bias, number and luminosity evolution, can be used to measure the large-scale properties of the Universe. We show how an optimal combination of tracers can be used to detect general-relativistic effects in the observed density of sources. We forecast for the detectability of these effects, as well as measurements of primordial non-Gaussianity and large-scale lensing magnification with current and upcoming large-scale structure experiments. In particular we quantify the significance of these detections in the short term with experiments such as the Dark Energy Survey (DES), and in the long term with the Large Synoptic Survey Telescope (LSST) and the Square Kilometre Array (SKA). We review the main observational challenges that must be overcome to carry out these measurements.

Using the 2MASS Photometric Redshift catalogue we perform a number of statistical tests aimed at detecting possible departures from statistical homogeneity and isotropy in the large-scale structure of the Universe. Making use of the angular homogeneity index, an observable proposed in a previous publication, as well as studying the scaling of the angular clustering and number counts with magnitude limit, we place constraints on the fractal nature of the galaxy distribution. We find that the statistical properties of our sample are in excellent agreement with the standard cosmological model, and that it reaches the homogeneous regime significantly faster than a class of fractal models with dimensions $D<2.75$. As part of our search for systematic effects, we also study the presence of hemispherical asymmetries in our data, finding no significant deviation beyond those allowed by the concordance model.

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

This chapter describes the assumed specifications and sensitivities for HI galaxy surveys with SKA1 and SKA2. It addresses the expected galaxy number densities based on available simulations as well as the clustering bias over the underlying dark matter. It is shown that a SKA1 HI galaxy survey should be able to find around $5\times 10^6$ galaxies over 5,000 deg$^2$ (up to $z\sim 0.8$), while SKA2 should find $\sim 10^9$ galaxies over 30,000 deg$^2$ (up to $z\sim 2.5$). The numbers presented here have been used throughout the cosmology chapters for forecasting.

21cm intensity mapping experiments aim to observe the diffuse neutral hydrogen (HI) distribution on large scales which traces the Cosmic structure. The Square Kilometre Array (SKA) will have the capacity to measure the 21cm signal over a large fraction of the sky. However, the redshifted 21cm signal in the respective frequencies is faint compared to the Galactic foregrounds produced by synchrotron and free-free electron emission. In this article, we review selected foreground subtraction methods suggested to effectively separate the 21cm signal from the foregrounds with intensity mapping simulations or data. We simulate an intensity mapping experiment feasible with SKA phase 1 including extragalactic and Galactic foregrounds. We give an example of the residuals of the foreground subtraction with a independent component analysis and show that the angular power spectrum is recovered within the statistical errors on most scales. Additionally, the scale of the Baryon Acoustic Oscillations is shown to be unaffected by foreground subtraction.

We investigate the dependence of the mass function of dark-matter haloes on their environment within the cosmic web of large-scale structure. A dependence of the halo mass function on large-scale mean density is a standard element of cosmological theory, allowing mass-dependent biasing to be understood via the peak-background split. On the assumption of a Gaussian density field, this analysis can be extended to ask how the mass function depends on the geometrical environment: clusters, filaments, sheets and voids, as classified via the tidal tensor (the Hessian matrix of the gravitational potential). In linear theory, the problem can be solved exactly, and the result is attractively simple: the conditional mass function has no explicit dependence on the local tidal field, and is a function only of the local density on the filtering scale used to define the tidal tensor. There is nevertheless a strong implicit predicted dependence on geometrical environment, because the local density couples statistically to the derivatives of the potential. We compute the predictions of this model and study the limits of their validity by comparing them to results deduced empirically from $N$-body simulations. We have verified that, to a good approximation, the abundance of haloes in different environments depends only on their densities, and not on their tidal structure. In this sense we find relative differences between halo abundances in different environments with the same density which are smaller than 13%. Furthermore, for sufficiently large filtering scales, the agreement with the theoretical prediction is good, although there are important deviations from the Gaussian prediction at small, non-linear scales. We discuss how to obtain improved predictions in this regime, using the 'effective-universe' approach.

We make use of a large set of fast simulations of an intensity mapping experiment with characteristics similar to those expected of the Square Kilometre Array (SKA) in order to study the viability and limits of blind foreground subtraction techniques. In particular, we consider different approaches: polynomial fitting, principal component analysis (PCA) and independent component analysis (ICA). We review the motivations and algorithms for the three methods, and show that they can all be described, using the same mathematical framework, as different approaches to the blind source separation problem. We study the efficiency of foreground subtraction both in the angular and radial (frequency) directions, as well as the dependence of this efficiency on different instrumental and modelling parameters. For well-behaved foregrounds and instrumental effects we find that foreground subtraction can be successful to a reasonable level on most scales of interest. We also quantify the effect that the cleaning has on the recovered signal and power spectra. Interestingly, we find that the three methods yield quantitatively similar results, with PCA and ICA being almost equivalent.

We present a code to generate mock observations of 21 cm intensity mapping experiments. The emphasis of the code is on reducing the computational cost of running a full-blown simulation, trading computational time for accuracy. The code can be used to generate independent realizations of the cosmological signal and foregrounds, which are necessary, for instance, in order to obtain realistic forecasts for future intensity mapping experiments. The code is able to reproduce the correct angular and radial clustering pattern for the cosmological signal, including redshift-space distortions, lightcone evolution and bias. Furthermore, it is possible to simulate a variety of foregrounds, including the potentially problematic polarized synchrotron emission.

We study the possibility of detecting the transition to homogeneity using photometric redshift catalogs. Our method is based on measuring the fractality of the projected galaxy distribution, using angular distances, and relies only on observable quantites. It thus provides a way to test the Cosmological Principle in a model-independent unbiased way. We have tested our method on different synthetic inhomogeneous catalogs, and shown that it is capable of discriminating some fractal models with relatively large fractal dimensions, in spite of the loss of information due to the radial projection. We have also studied the influence of the redshift bin width, photometric redshift errors, bias, non-linear clustering, and surveyed area, on the angular homogeneity index H2 ({\theta}) in a {\Lambda}CDM cosmology. The level to which an upcoming galaxy survey will be able to constrain the transition to homogeneity will depend mainly on the total surveyed area and the compactness of the surveyed region. In particular, a Dark Energy Survey (DES)-like survey should be able to easily discriminate certain fractal models with fractal dimensions as large as D2 = 2.95. We believe that this method will have relevant applications for upcoming large photometric redshift surveys, such as DES or the Large Synoptic Survey Telescope (LSST).

In this paper we present a new method to extract cosmological parameters using the radial scale of the Baryon Acoustic Oscillations as a standard ruler in deep galaxy surveys. The method consists in an empirical parametrization of the radial 2-point correlation function, which provides a robust and precise extraction of the sound horizon scale. Moreover, it uses data from galaxy surveys in a manner that is fully cosmology independent and therefore, unbiased. A study of the main systematic errors and the validation of the method in cosmological simulations are also presented, showing that the measurement is limited only by cosmic variance. We then study the full information contained in the Baryon Acoustic Oscillations, obtaining that the combination of the radial and angular determinations of this scale is a very sensitive probe of cosmological parameters, able to set strong constraints on the dark energy properties, even without combining it with any other probe.

In the advent of new large galaxy surveys, which will produce enormous datasets with hundreds of millions of objects, new computational techniques are necessary in order to extract from them any two-point statistic, the computational time of which grows with the square of the number of objects to be correlated. Fortunately technology now provides multiple means to massively parallelize this problem. Here we present a free-source code specifically designed for this kind of calculations. Two implementations are provided: one for execution on shared-memory machines using OpenMP and one that runs on graphical processing units (GPUs) using CUDA. The code is available at http://members.ift.uam-csic.es/dmonge/CUTE.html.

We study the non-linear gravitational collapse of dark matter into halos through numerical N-body simulations of Lemaitre-Tolman-Bondi void models. We extend the halo mass function formalism to these models in a consistent way. This extension not only compares well with the simulated data at all times and radii, but it also gives interesting clues about the impact of the background shear on the growth of perturbations. Our results give hints about the possibility of constraining the background shear via cluster number counts, which could then give rise to strong constraints on general inhomogeneous models, of any scale.

We propose a new method for cosmological parameters extraction using the baryon acoustic oscillation scale as a standard ruler in deep galaxy surveys with photometric determination of redshifts. The method consists in a simple empirical parametric fit to the angular 2-point correlation function w(theta). It is parametrized as a power law to describe the continuum plus a Gaussian to describe the BAO bump. The location of the Gaussian is used as the basis for the measurement of the sound horizon scale. This method, although simple, actually provides a robust estimation, since the inclusion of the power law and the use of the Gaussian removes the shifts which affect the local maximum. We discuss the effects of projection bias, non-linearities, redshift space distortions and photo-z precision, and apply our method to a mock catalog of the Dark Energy Survey, built upon a large N-body simulation provided by the MICE collaboration. We discuss the main systematic errors associated to our method and show that they are dominated by the photo-z uncertainty.

We perform numerical simulations of large scale structure evolution in an inhomogeneous Lemaitre-Tolman-Bondi (LTB) model of the Universe. We follow the gravitational collapse of a large underdense region (a void) in an otherwise flat matter-dominated Einstein-deSitter model. We observe how the (background) density contrast at the centre of the void grows to be of order one, and show that the density and velocity profiles follow the exact non-linear LTB solution to the full Einstein equations for all but the most extreme voids. This result seems to contradict previous claims that fully relativistic codes are needed to properly handle the non-linear evolution of large scale structures, and that local Newtonian dynamics with an explicit expansion term is not adequate. We also find that the (local) matter density contrast grows with the scale factor in a way analogous to that of an open universe with a value of the matter density OmegaM(r) corresponding to the appropriate location within the void.