Are the stellar-mass merging binary black holes, recently detected by their gravitational wave signal, of stellar or primordial origin? Answering this question will have profound implications for our understanding of the Universe, including the nature of dark matter, the early Universe and stellar evolution. We build on the idea that the clustering properties of merging binary black holes can provide information about binary formation mechanisms and origin. The cross-correlation of galaxy with gravitational wave catalogues carries information about whether black hole mergers trace more closely the distribution of dark matter -- indicative of primordial origin -- or that of stars harboured in luminous and massive galaxies -- indicative of a stellar origin. We forecast the detectability of such signal for several forthcoming and future gravitational wave interferometers and galaxy surveys, including, for the first time in such analyses, an accurate modelling for the different merger rates, lensing magnification and other general relativistic effects. Our results show that forthcoming experiments could allow us to test most of the parameter space of the still viable models investigated, and shed more light on the issue of binary black hole origin and evolution.

In the era of precision cosmology, establishing the correct magnitude of statistical errors in cosmological parameters is of crucial importance. However, widely used approximations in galaxy surveys analyses can lead to parameter uncertainties that are grossly mis-estimated. These approximations can be introduced at three different levels: in the form of the likelihood, in the theoretical modelling of the observable and in the numerical computation of the observable. Their consequences are important both in data analysis through e.g., Markov Chain Monte Carlo parameter inference, and when survey instrument and strategy are designed and their constraining power on cosmological parameters is forecasted, for instance using Fisher matrix analyses. In this work, considering the galaxy angular power spectrum as the target observable, we report one example for each of such three categories and derive their impact on the estimated parameters uncertainties and degeneracies. We show that there are cases where these commonly used approximations lead, perhaps counter-intuitively, to unacceptably large mis-estimates of parameters errors and correlations. Furthermore, we stress how these approximations might even spoil the benefits of the nascent multi-tracer and multi-messenger cosmology. Hence we recommend that the type of analysis presented here should be repeated for every approximation adopted in survey design or data analysis, to quantify how it may affect the results. To this aim, we have developed \texttt{Multi\_CLASS}, a new extension of \texttt{CLASS} that includes the angular power spectrum for multiple (galaxy and other tracers such as gravitational waves) populations. The public release of \texttt{Multi\_CLASS} is associated with this paper.

Cosmological parameter estimation from forthcoming experiments promise to reach much greater precision than current constraints. As statistical errors shrink, the required control over systematic errors increases. Therefore, models or approximations that were sufficiently accurate so far, may introduce significant systematic biases in the parameter best-fit values and jeopardize the robustness of cosmological analyses. We present a general expression to estimate a priori the systematic error introduced in parameter inference due to the use of insufficiently good approximations in the computation of the observable of interest or the assumption of an incorrect underlying model. Although this methodology can be applied to measurements of any scientific field, we illustrate its power by studying the effect of modeling the angular galaxy power spectrum incorrectly. We also introduce Multi_CLASS, a new, public modification of the Boltzmann code CLASS, which includes the possibility to compute angular cross-power spectra for two different tracers. We find that significant biases in most of the cosmological parameters are introduced if one assumes the Limber approximation or neglects lensing magnification in modern galaxy survey analyses, and the effect is in general larger for the multi-tracer case, especially for the parameter controlling primordial non-Gaussianity of the local type, $f_{\rm NL}$.

The detection and characterization of the Stochastic Gravitational Wave Background (SGWB) is one of the main goals of Gravitational Wave (GW) experiments. The observed SGWB will be the combination of GWs from cosmological (as predicted by many models describing the physics of the early Universe) and astrophysical origins, which will arise from the superposition of GWs from unresolved sources whose signal is too faint to be detected. Therefore, it is important to have a proper modeling of the astrophysical SGWB (ASGWB) in order to disentangle the two signals; moreover, this will provide additional information on astrophysical properties of compact objects. Applying the Cosmic Rulers formalism, we compute the observed ASGWB angular power spectrum, hence using gauge invariant quantities, accounting for all effects intervening between the source and the observer. These are the so-called projection effects, which include Kaiser, Doppler and gravitational potentials effect. Our results show that these projection effects are the most important at the largest scales, and they contribute to up to tens of percent of the angular power spectrum amplitude, with the Kaiser term being the largest at all scales. While the exact impact of these results will depend on instrumental and astrophysical details, a precise theoretical modeling of the ASGWB will necessarily need to include all these projection effects.

Constraints on dark matter annihilation or decay offer unique insights into the nature of dark matter. We illustrate how surveys dedicated to detect the highly redshifted 21 cm signal from the dark ages will offer a new window into properties of particle dark matter. The 21 cm intensity mapping signal and its fluctuations are sensitive to energy injection from annihilation or decay of long-lived particles in a way that is complementary to other probes. We present forecasted constraints from forthcoming and next-generation radio surveys. We show that, while SKA might be capable of a detection for some cases, the most promising opportunity to detect dark matter in the 21 cm intensity mapping signal is with a futuristic radio array on the lunar far-side, with the potential to detect a signal many orders of magnitude weaker than current or maximal constraints from other probes.

We study the clustering properties of dark matter halos in real- and redshift-space in cosmologies with massless and massive neutrinos through a large set of state-of-the-art N-body simulations. We provide quick and easy-to-use prescriptions for the halo bias on linear and mildly non-linear scales, both in real and redshift space, which are valid also for massive neutrinos cosmologies. Finally we present a halo bias emulator,$\textbf{BE-HaPPY}$, calibrated on the N-body simulations, which is fast enough to be used in the standard Markov Chain Monte Carlo approach to cosmological inference. For a fiducial standard $\Lambda$CDM cosmology $\textbf{BE-HaPPY}$ provides percent or sub-percent accuracy on the scales of interest (linear and well into the mildly non-linear regime), meeting therefore for the halo-bias the accuracy requirements for the analysis of next-generation large--scale structure surveys.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

In the model where Primordial Black Holes (PBHs) form from large primordial curvature (C) perturbations, i.e., CPBHs, constraints on PBH abundance provide in principle constraints on the primordial curvature power spectrum. This connection however depends necessarily on the details of PBH formation mechanism. In this paper we provide, for the first time, constraints on the primordial curvature power spectrum from the latest limits on PBH abundance, taking into account all the steps from gravitational collapse in real space to PBH formation. In particular, we use results from numerical relativity simulations and peak theory to study the conditions for PBH formation for a range of perturbation shapes, including non-linearities, perturbation profile and a careful treatment of smoothing and filtering scales. We then obtain updated PBH formation conditions and translate that into primordial spectrum constraints for a wide range of shapes and abundances. These updated constraints cover a range of scales not probed by other cosmological observables. Our results show that the correct and accurate modelling of non-linearities, filtering and typical perturbation profile, is crucial for deriving meaningful cosmological implications.

Waveforms of gravitational waves provide information about a variety of parameters for the binary system merging. However, standard calculations have been performed assuming a FLRW universe with no perturbations. In reality this assumption should be dropped: we show that the inclusion of cosmological perturbations translates into corrections to the estimate of astrophysical parameters derived for the merging binary systems. We compute corrections to the estimate of the luminosity distance due to velocity, volume, lensing and gravitational potential effects. Our results show that the amplitude of the corrections will be negligible for current instruments, mildly important for experiments like the planned DECIGO, and very important for future ones such as the Big Bang Observer.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

The renewed interest in the possibility that primordial black holes (PBHs) may constitute a significant part of the dark matter has motivated revisiting old observational constraints, as well as developing new ones. We present new limits on the PBH abundance, from a comprehensive analysis of high-resolution, high-redshift Lyman-$\alpha$ forest data. Poisson fluctuations in the PBH number density induce a small-scale power enhancement which departs from the standard cold dark matter prediction. Using a grid of hydrodynamic simulations exploring different values of astrophysical parameters, {we obtain a marginalized upper limit on the PBH mass of $f_{\rm PBH}M_{\rm PBH} \sim 60~M_{\odot}$ at $2\sigma$, when a Gaussian prior on the reionization redshift is imposed, preventing its posterior distribution to peak on very high values, which are disfavoured by the most recent estimates obtained both through Cosmic Microwave Background and Inter-Galactic Medium observations. Such bound weakens to $f_{\rm PBH}M_{\rm PBH} \sim 170~M_{\odot}$, when a conservative flat prior is instead assumed. Both limits significantly improves previous constraints from the same physical observable.} We also extend our predictions to non-monochromatic PBH mass distributions, ruling out large regions of the parameter space for some of the most viable PBH extended mass functions.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The Evolutionary Map of the Universe (EMU) is an all-sky survey in radio-continuum which uses the Australian SKA Pathfinder (ASKAP). Using galaxy angular power spectrum and the integrated Sachs-Wolfe effect, we study the potential of EMU to constrain models beyond $\Lambda$CDM (i.e., local primordial non-Gaussianity, dynamical dark energy, spatial curvature and deviations from general relativity), for different design sensitivities. We also include a multi-tracer analysis, distinguishing between star-forming galaxies and galaxies with an active galactic nucleus, to further improve EMU's potential. We find that EMU could measure the dark energy equation of state parameters around 35\% more precisely than existing constraints, and that the constraints on $f_{\rm NL}$ and modified gravity parameters will improve up to a factor $\sim2$ with respect to Planck and redshift space distortions measurements. With this work we demonstrate the promising potential of EMU to contribute to our understanding of the Universe.

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.

It is broadly accepted that Supermassive Black Holes (SMBHs) are located in the centers of most massive galaxies, although there is still no convincing scenario for the origin of their massive seeds. It has been suggested that primordial black holes (PBHs) of masses $\gtrsim 10^{2} M_\odot$ may provide such seeds, which would grow to become SMBHs. We suggest an observational test to constrain this hypothesis: gas accretion around PBHs during the cosmic dark ages powers the emission of high energy photons which would modify the spin temperature as measured by 21cm Intensity Mapping (IM) observations. We model and compute their contribution to the standard sky-averaged signal and power spectrum of 21cm IM, accounting for its substructure and angular dependence for the first time. If PBHs exist, the sky-averaged 21cm IM signal in absorption would be higher, while we expect an increase in the power spectrum for $\ell~\gtrsim 10^2-10^3$. We also forecast PBH detectability and measurement errors in the abundance and Eddington ratios for different fiducial parameter configurations for various future experiments, ranging from SKA to a futuristic radio array on the dark side of the Moon. While the SKA could provide a detection, only a more ambitious experiment would provide accurate measurements.

We forecast constraints on primordial non-Gaussianity (PNG) and bias parameters from measurements of galaxy power spectrum and bispectrum in future radio continuum and optical surveys. In the galaxy bispectrum, we consider a comprehensive list of effects, including the bias expansion for non-Gaussian initial conditions up to second order, redshift space distortions, redshift uncertainties and theoretical errors. These effects are all combined in a single PNG forecast for the first time. Moreover, we improve the bispectrum modelling over previous forecasts, by accounting for trispectrum contributions. All effects have an impact on final predicted bounds, which varies with the type of survey. We find that the bispectrum can lead to improvements up to a factor $\sim 5$ over bounds based on the power spectrum alone, leading to significantly better constraints for local-type PNG, with respect to current limits from \textit{Planck}. Future radio and photometric surveys could obtain a measurement error of $\sigma(f_{\mathrm{NL}}^{\mathrm{loc}}) \approx 0.2$. In the case of equilateral PNG, galaxy bispectrum can improve upon present bounds only if significant improvements in the redshift determinations of future, large volume, photometric or radio surveys could be achieved. For orthogonal non-Gaussianity, expected constraints are generally comparable to current ones.

An interesting test on the nature of the Universe is to measure the global spatial curvature of the metric in a model independent way, at a level of $|\Omega_k|<10^{-4}$, or, if possible, at the cosmic variance level of the amplitude of the CMB fluctuations $|\Omega_k|\approx10^{-5}$. A limit of $|\Omega_k|<10^{-4}$ would yield stringent tests on several models of inflation. Further, improving the constraint by an order of magnitude would help in reducing "model confusion" in standard parameter estimation. Moreover, if the curvature is measured to be at the value of the amplitude of the CMB fluctuations, it would offer a powerful test on the inflationary paradigm and would indicate that our Universe must be significantly larger than the current horizon. On the contrary, in the context of standard inflation, measuring a value above CMB fluctuations will lead us to conclude that the Universe is not much larger than the current observed horizon, this can also be interpreted as the presence of large fluctuations outside the horizon. However, it has proven difficult, so far, to find observables that can achieve such level of accuracy, and, most of all, be model-independent. Here we propose a method that can in principle achieve that, this is done by making minimal assumptions and using distance probes that are cosmology-independent: gravitational waves, redshift drift and cosmic chronometers. We discuss what kind of observations are needed in principle to achieve the desired accuracy.

Given the important role that the galaxy bispectrum has recently acquired in cosmology and the scale and precision of forthcoming galaxy clustering observations, it is timely to derive the full expression of the large-scale bispectrum going beyond approximated treatments which neglect integrated terms or higher-order bias terms or use the Limber approximation. On cosmological scales, relativistic effects that arise from observing on the past light-cone alter the observed galaxy number counts, therefore leaving their imprints on N-point correlators at all orders. In this paper we compute for the first time the bispectrum including all general relativistic, local and integrated, effects at second order, the tracers' bias at second order, geometric effects as well as the primordial non-Gaussianity contribution. This is timely considering that future surveys will probe scales comparable to the horizon where approximations widely used currently may not hold; neglecting these effects may introduce biases in estimation of cosmological parameters as well as primordial non-Gaussianity.

The model in which Primordial Black Holes (PBHs) constitute a non-negligible fraction of the dark matter has (re)gained popularity after the first detections of binary black hole mergers. Most of the observational constraints to date have been derived assuming a single mass for all the PBHs, although some more recent works tried to generalize constraints to the case of extended mass functions. Here we derive a general methodology to obtain constraints for any PBH Extended Mass Distribution (EMD) and any observables in the desired mass range. Starting from those obtained for a monochromatic distribution, we convert them into constraints for EMDs by using an equivalent, effective mass $M_{\rm eq}$ that depends on the specific observable. We highlight how limits of validity of the PBH modelling affect the EMD parameter space. Finally, we present converted constraints on the total abundance of PBH from microlensing, stellar distribution in ultra-faint dwarf galaxies and CMB accretion for Lognormal and Power Law mass distributions, finding that EMD constraints are generally stronger than monochromatic ones.

The possibility that a relevant fraction of the dark matter might be comprised of Primordial Black Holes (PBHs) has been seriously reconsidered after LIGO's detection of a $\sim 30 M_{\odot}$ binary black holes merger. Despite the strong interest in the model, there is a lack of studies on possible cosmological implications and effects on cosmological parameters inference. We investigate correlations with the other standard cosmological parameters using cosmic microwave background observations, finding significant degeneracies, especially with the tilt of the primordial power spectrum and the sound horizon at radiation drag. However, these degeneracies can be greatly reduced with the inclusion of small scale polarization data. We also explore if PBHs as dark matter in simple extensions of the standard $\Lambda$CDM cosmological model induces extra degeneracies, especially between the additional parameters and the PBH's ones. Finally, we present cosmic microwave background constraints on the fraction of dark matter in PBHs, not only for monochromatic PBH mass distributions but also for popular extended mass distributions. Our results show that extended mass distribution's constraints are tighter, but also that a considerable amount of constraining power comes from the high-$\ell$ polarization data. Moreover, we constrain the shape of such mass distributions in terms of the correspondent constraints on the PBH mass fraction.

It has been recently suggested that small mass black holes (BHs) may become unstable due to quantum-gravitational effects and eventually decay, producing radiation, on a timescale shorter than the Hawking evaporation time. We argue that the existence of a population of low-mass Primordial Black Holes (PBHs) acting as a fraction of the Universe dark matter component can be used to test proposed models of quantum decay of BHs via their effect on galaxy number counts. We study what constraints future galaxy clustering measurements can set on quantum-gravity parameters governing the BH lifetime and PBH abundance. In case of no detection of such effects, this would rule out either the existence of a non-negligible number of small PBHs, or the BH quantum decay scenario (or both). In case of independent observations of PBHs, the observables discussed here could be used to study the quantum effects that modify the final fate of BHs.

Single-field slow-roll inflation predicts a nearly scale-free power spectrum of perturbations, as observed at the scales accessible to current cosmological experiments. This spectrum is slightly red, showing a tilt $(1-n_s)\sim 0.04$. A direct consequence of this tilt are nonvanishing runnings $\alpha_s=\mathrm d n_s/\mathrm d\log k$, and $\beta_s=\mathrm d\alpha_s/\mathrm d\log k$, which in the minimal inflationary scenario should reach absolute values of $10^{-3}$ and $10^{-5}$, respectively. In this work we calculate how well future surveys can measure these two runnings. We consider a Stage-4 (S4) CMB experiment and show that it will be able to detect significant deviations from the inflationary prediction for $\alpha_s$, although not for $\beta_s$. Adding to the S4 CMB experiment the information from a WFIRST-like or a DESI-like survey improves the sensitivity to the runnings by $\sim$ 20\%, and 30\%, respectively. A spectroscopic survey with a billion objects, such as the SKA, will add enough information to the S4 measurements to allow a detection of $\alpha_s=10^{-3}$, required to probe the single-field slow-roll inflationary paradigm. We show that only a very-futuristic interferometer targeting the dark ages will be capable of measuring the minimal inflationary prediction for $\beta_s$. The results of other probes, such as a stochastic background of gravitational waves observable by LIGO, the Ly-$\alpha$ forest, and spectral distortions, are shown for comparison. Finally, we study the claims that large values of $\beta_s$, if extrapolated to the smallest scales, can produce primordial black holes of tens of solar masses, which we show to be easily testable by the S4 CMB experiment.

The relation between the halo field and the matter fluctuations (halo bias), in the presence of massive neutrinos depends on the total neutrino mass, massive neutrinos introduce an additional scale-dependence of the bias which is usually neglected in cosmological analyses. We investigate the magnitude of the systematic effect on interesting cosmological parameters induced by neglecting this scale dependence, finding that while it is not a problem for current surveys, it is non-negligible for future, denser or deeper ones depending on the neutrino mass, the maximum scale used for the analyses and the details of the nuisance parameters considered. However there is a simple recipe to account for the bulk of the effect as to make it fully negligible, which we illustrate and advocate should be included in analysis of forthcoming large-scale structure surveys.

In the next decade, new astrophysical instruments will deliver the first large-scale maps of gravitational waves and radio sources. Therefore, it is timely to investigate the possibility to combine them to provide new and complementary ways to study the Universe. Using simulated catalogues appropriate to the planned surveys, it is possible to predict measurements of the cross-correlation between radio sources and GW maps and the effects of a stochastic gravitational wave background on galaxy maps. Effects of GWs on the large scale structure of the Universe can be used to investigate the nature of the progenitors of merging BHs, the validity of Einstein's General Relativity, models for dark energy, and detect a stochastic background of GW. The results obtained show that the galaxy-GW cross-correlation can provide useful information in the near future, while the detection of tensor perturbation effects on the LSS will require instruments with capabilities beyond the currently planned next generation of radio arrays. Nevertheless, any information from the combination of galaxy surveys with GW maps will help provide additional information for the newly born gravitational wave astronomy.

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

We demonstrate that observations lacking reliable redshift information, such as photometric and radio continuum surveys, can produce robust measurements of cosmological parameters when empowered by clustering-based redshift estimation. This method infers the redshift distribution based on the spatial clustering of sources, using cross-correlation with a reference dataset with known redshifts. Applying this method to the existing SDSS photometric galaxies, and projecting to future radio continuum surveys, we show that sources can be efficiently divided into several redshift bins, increasing their ability to constrain cosmological parameters. We forecast constraints on the dark-energy equation-of-state and on local non-gaussianity parameters. We explore several pertinent issues, including the tradeoff between including more sources versus minimizing the overlap between bins, the shot-noise limitations on binning, and the predicted performance of the method at high redshifts. Remarkably, we find that, once this technique is implemented, constraints on dynamical dark energy from the SDSS imaging catalog can be competitive with, or better than, those from the spectroscopic BOSS survey and even future planned experiments. Further, constraints on primordial non-Gaussianity from future large-sky radio-continuum surveys can outperform those from the Planck CMB experiment, and rival those from future spectroscopic galaxy surveys. The application of this method thus holds tremendous promise for cosmology.

It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those produced by more traditional astrophysical mechanisms. These PBH binaries are formed on highly eccentric orbits and can then merge on timescales that in some cases are years or less, retaining some eccentricity in the last seconds before the merger. This is to be contrasted with massive-stellar-binary, globular-cluster, or other astrophysical origins for binary black holes (BBHs) in which the orbits have very effectively circularized by the time the BBH enters the observable LIGO window. Here we discuss the features of the gravitational-wave signals that indicate this eccentricity and forecast the sensitivity of LIGO and the Einstein Telescope to such effects. We show that if PBHs make up the dark matter, then roughly one event should have a detectable eccentricity given LIGO's expected sensitivity and observing time of six years. The Einstein Telescope should see $O(10)$ such events after ten years.

Galaxy clustering on very large scales can be probed via the 2-point correlation function in the general case of wide and deep separations, including all the lightcone and relativistic effects. Using our recently developed formalism, we analyze the behavior of the local and integrated contributions and how these depend on redshift range, linear and angular separations and luminosity function. Relativistic corrections to the local part of the correlation can be non-negligible but they remain generally sub-dominant. On the other hand, the additional correlations arising from lensing convergence and time-delay effects can become very important and even dominate the observed total correlation function. We investigate different configurations formed by the observer and the pair of galaxies, and we find that the case of near-radial large-scale separations is where these effects will be the most important.

We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $20\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they radiate enough energy in gravitational waves to become gravitationally bound. The bound BHs will rapidly spiral inward due to emission of gravitational radiation and ultimately merge. Uncertainties in the rate for such events arise from our imprecise knowledge of the phase-space structure of galactic halos on the smallest scales. Still, reasonable estimates span a range that overlaps the $2-53$ Gpc$^{-3}$ yr$^{-1}$ rate estimated from GW150914, thus raising the possibility that LIGO has detected PBH dark matter. PBH mergers are likely to be distributed spatially more like dark matter than luminous matter and have no optical nor neutrino counterparts. They may be distinguished from mergers of BHs from more traditional astrophysical sources through the observed mass spectrum, their high ellipticities, or their stochastic gravitational wave background. Next generation experiments will be invaluable in performing these tests.

We modified the CLASS code in order to include relativistic galaxy number counts in spatially curved geometries; we present the formalism and study the effect of relativistic corrections on spatial curvature. The new version of the code is now publicly available. Using a Fisher matrix analysis, we investigate how measurements of the spatial curvature parameter $\Omega_K$ with future galaxy surveys are affected by relativistic effects, which influence observations of the large scale galaxy distribution. These effects include contributions from cosmic magnification, Doppler terms and terms involving the gravitational potential. As an application, we consider angle and redshift dependent power spectra, which are especially well suited for model independent cosmological constraints. We compute our results for a representative deep, wide and spectroscopic survey, and our results show the impact of relativistic corrections on the spatial curvature parameter estimation. We show that constraints on the curvature parameter may be strongly biased if, in particular, cosmic magnification is not included in the analysis. Other relativistic effects turn out to be subdominant in the studied configuration. We analyze how the shift in the estimated best-fit value for the curvature and other cosmological parameters depends on the magnification bias parameter, and find that significant biases are to be expected if this term is not properly considered in the analysis.

Most inflationary models predict primordial perturbations to be statistically isotropic and homogeneous. Cosmic microwave background (CMB) observations, however, indicate a possible departure from statistical isotropy in the form of a dipolar power modulation at large angular scales. Alternative models of inflation, beyond the simplest single-field slow-roll models, can generate a small power asymmetry, consistent with these observations. Observations of clustering of quasars show, however, agreement with statistical isotropy at much smaller angular scales. Here, we propose to use off-diagonal components of the angular power spectrum of the 21-cm fluctuations during the dark ages to test this power asymmetry. We forecast results for the planned SKA radio array, a future radio array, and the cosmic-variance-limited case as a theoretical proof of principle. Our results show that the 21-cm line power spectrum will enable access to information at very small scales and at different redshift slices, thus improving upon the current CMB constraints by $\sim 2$ orders of magnitude for a dipolar asymmetry and by $\sim 1-3$ orders of magnitude for a quadrupolar asymmetry case.

We investigate a possible method for determining the progenitors of black hole (BH) mergers observed via their gravitational wave (GW) signal. We argue that measurements of the cross-correlation of the GW events with overlapping galaxy catalogs may provide an additional tool in determining if BH mergers trace the stellar mass of the Universe, as would be expected from mergers of the endpoints of stellar evolution. If on the other hand the BHs are of primordial origin, as has been recently suggested, their merging would be preferentially hosted by lower biased objects, and thus have a lower cross-correlation with luminous galaxies. Here we forecast the expected precision of the cross-correlation measurement for current and future GW detectors such as LIGO and the Einstein Telescope. We then predict how well these instruments can distinguish the model that identifies high-mass BH-BH mergers as the merger of primordial black holes that constitute the dark matter in the Universe from more traditional astrophysical sources.

Next-generation spectroscopic surveys will map the large-scale structure of the observable universe, using emission line galaxies as tracers. While each survey will map the sky with a specific emission line, interloping emission lines can masquerade as the survey's intended emission line at different redshifts. Interloping lines from galaxies that are not removed can contaminate the power spectrum measurement, mixing correlations from various redshifts and diluting the true signal. We assess the potential for power spectrum contamination, finding that an interloper fraction worse than 0.2% could bias power spectrum measurements for future surveys by more than 10% of statistical errors, while also biasing power spectrum inferences. We also construct a formalism for predicting cosmological parameter bias, demonstrating that a 0.15%-0.3% interloper fraction could bias the growth rate by more than 10% of the error, which can affect constraints on gravity upcoming surveys. We use the COSMOS Mock Catalog (CMC), with the emission lines re-scaled to better reproduce recent data, to predict potential interloper fractions for the Prime Focus Spectrograph (PFS) and the Wide-Field InfraRed Survey Telescope (WFIRST). We find that secondary line identification, or confirming galaxy redshifts by finding correlated emission lines, can remove interlopers for PFS. For WFIRST, we use the CMC to predict that the 0.2% target can be reached for the WFIRST H$\alpha$ survey, but sensitive optical and near-infrared photometry will be required. For the WFIRST [OIII] survey, the predicted interloper fractions reach several percent and their effects will have to be estimated and removed statistically (e.g. with deep training samples). (Abridged)

The auto-correlation between two members of a galaxy population is symmetric under the interchange of the two galaxies being correlated. The cross-correlation between two different types of galaxies, separated by a vector $\bf{r}$, is not necessarily the same as that for a pair separated by $-\bf{r}$. Local anisotropies in the two-point cross-correlation function may thus indicate a specific direction which when mapped as a function of position trace out a vector field. This vector field can then be decomposed into longitudinal and transverse components, and those transverse components written as positive- and negative-helicity components. A locally asymmetric cross-correlation of the longitudinal type arises naturally in halo clustering, even with Gaussian initial conditions, and could be enhanced with local-type non-Gaussianity. Early-Universe scenarios that introduce a vector field may also give rise to such effects. These antisymmetric cross-correlations also provide a new possibility to seek a preferred cosmic direction correlated with the hemispherical power asymmetry in the cosmic microwave background and to seek a preferred location associated with the CMB cold spot. New ways to seek cosmic parity breaking are also possible.

We study the parity-odd part (that we shall call Doppler term) of the linear galaxy two-point correlation function that arises from wide-angle, velocity, Doppler lensing and cosmic acceleration effects. As it is important at low redshift and at large angular separations, the Doppler term is usually neglected in the current generation of galaxy surveys. For future wide-angle galaxy surveys such as Euclid, SPHEREx and SKA, however, we show that the Doppler term must be included. The effect of these terms is dominated by the magnification due to relativistic aberration effects and the slope of the galaxy redshift distribution and it generally mimics the effect of the local type primordial non-Gaussianity with the effective nonlinearity parameter $f_{\rm NL}^{\rm eff}$ of a few, we show that this would affect forecasts on measurements of $f_{\rm NL}$ at low-redshift. Our results show that a survey at low redshift with large number density over a wide area of the sky could detect the Doppler term with a signal-to-noise ratio of $\sim 1-20$, depending on survey specifications.

We extend previous analyses of wide-angle correlations in the galaxy power spectrum in redshift space to include all general relativistic effects. These general relativistic corrections to the standard approach become important on large scales and at high redshifts, and they lead to new terms in the wide-angle correlations. We show that in principle the new terms can produce corrections of nearly 10 % on Gpc scales over the usual Newtonian approximation. General relativistic corrections will be important for future large-volume surveys such as SKA and Euclid, although the problem of cosmic variance will present a challenge in observing this.

We investigate how well future large-scale radio surveys could measure different shapes of primordial non-Gaussianity; in particular we focus on angle-dependent non-Gaussianity arising from primordial anisotropic sources, whose bispectrum has an angle dependence between the three wavevectors that is characterized by Legendre polynomials $\mathcal{P}_L$ and expansion coefficients $c_L$. We provide forecasts for measurements of galaxy power spectrum, finding that Large-Scale Structure (LSS) data could allow measurements of primordial non-Gaussianity competitive or improving upon current constraints set by CMB experiments, for all the shapes considered. We argue that the best constraints will come from the possibility to assign redshift information to radio galaxy surveys, and investigate a few possible scenarios for the EMU and SKA surveys. A realistic (futuristic) modeling could provide constraints of $f_{\rm NL}^{\rm loc} \approx 1 (0.5)$ for the local shape, $f_{\rm NL}$ of $\mathcal{O}(10) (\mathcal{O}(1))$ for the orthogonal, equilateral and folded shapes, and $c_{L=1} \approx 80 (2)$, $c_{L=2} \approx 400 (10)$ for angle-dependent non-Gaussianity. The more futuristic forecasts show the potential of LSS analyses to considerably improve current constraints on non-Gaussianity, and so on models of the primordial Universe. Finally, we find the minimum requirements that would be needed to reach $\sigma(c_{L=1})=10$, which can be considered as a typical (lower) value predicted by some (inflationary) models.

We investigate the cosmological dependence and the constraining power of large-scale galaxy correlations, including all redshift-distortions, wide-angle, lensing and gravitational potential effects on linear scales. We analyze the cosmological information present in the lensing convergence and in the gravitational potential terms describing the so-called "relativistic effects," and we find that, while smaller than the information contained in intrinsic galaxy clustering, it is not negligible. We investigate how neglecting them does bias cosmological measurements performed by future spectroscopic and photometric large-scale surveys such as SKA and Euclid. We perform a Fisher analysis using the CLASS code, modified to include scale-dependent galaxy bias and redshift-dependent magnification and evolution bias. Our results show that neglecting relativistic terms introduces an error in the forecasted precision in measuring cosmological parameters of the order of a few tens of percent, in particular when measuring the matter content of the Universe and primordial non-Gaussianity parameters. Therefore, we argue that radial correlations and integrated relativistic terms need to be taken into account when forecasting the constraining power of future large-scale number counts of galaxy surveys.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

The Planck CMB experiment has delivered the best constraints so far on primordial non-Gaussianity, ruling out early-Universe models of inflation that generate large non-Gaussianity. Although small improvements in the CMB constraints are expected, the next frontier of precision will come from future large-scale surveys of the galaxy distribution. The advantage of such surveys is that they can measure many more modes than the CMB -- in particular, forthcoming radio surveys with the SKA will cover huge volumes. Radio continuum surveys deliver the largest volumes, but with the disadvantage of no redshift information. In order to mitigate this, we use two additional observables. First, the integrated Sachs-Wolfe effect -- the cross-correlation of the radio number counts with the CMB temperature anisotropies -- helps to reduce systematics on the large scales that are sensitive to non-Gaussianity. Second, optical data allows for cross-identification in order to gain some redshift information. We show that, while the single redshift bin case can provide a sigma(fNL) ~ 20, and is therefore not competitive with current and future constraints on non-Gaussianity, a tomographic analysis could improve the constraints by an order of magnitude, even with only two redshift bins. A huge improvement is provided by the addition of high-redshift sources, so having cross-ID for high-z galaxies and an even higher-z radio tail is key to enabling very precise measurements of fNL. Our results show that SKA continuum surveys could provide constraints competitive with CMB and forthcoming optical surveys, potentially allowing a measurement of sigma(fNL) ~ 1 to be made. Moreover, these measurements would act as a useful check of results obtained with other probes at other redshift ranges with other methods.

We investigate the possibility to detect the integrated Sachs-Wolfe (ISW) effect by cross-correlating 21-cm surveys at high redshifts with galaxies, in a way similar to the usual CMB-galaxy cross-correlation. The high-redshift 21-cm signal is dominated by CMB photons that travel freely without interacting with the intervening matter, and hence its late-time ISW signature should correlate extremely well with that of the CMB at its peak frequencies. Using the 21-cm temperature brightness instead of the CMB would thus be a further check of the detection of the ISW effect, measured by different instruments at different frequencies and suffering from different systematics. We also study the ISW effect on the photons that are scattered by HI clouds. We show that a detection of the unscattered photons is achievable with planned radio arrays, while one using scattered photons will require advanced radio interferometers, either an extended version of the planned Square Kilometre Array or futuristic experiments such as a lunar radio array.

The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time.

Employing a nonparametric approach of the principal component analysis (PCA), we forecast the future constraint on the equation of state $w(z)$ of dark energy, and on the effective Newton constant $\mu(k,z)$, which parameterise the effect of modified gravity, using the planned SKA HI galaxy survey. Combining with the simulated data of Planck and Dark Energy Survey (DES), we find that SKA Phase 1 (SKA1) and SKA Phase 2 (SKA2) can well constrain $3$ and $5$ eigenmodes of $w(z)$ respectively. The errors of the best measured modes can be reduced to 0.04 and 0.023 for SKA1 and SKA2 respectively, making it possible to probe dark energy dynamics. On the other hand, SKA1 and SKA2 can constrain $7$ and $20$ eigenmodes of $\mu(k,z)$ respectively within 10\% sensitivity level. Furthermore, 2 and 7 modes can be constrained within sub percent level using SKA1 and SKA2 respectively. This is a significant improvement compared to the combined datasets without SKA.

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

The Square Kilometer Array (SKA) has the potential to produce galaxy redshift surveys which will be competitive with other state of the art cosmological experiments in the next decade. In this chapter we summarise what capabilities the first and the second phases of the SKA will be able to achieve in its current state of design. We summarise the different cosmological experiments which are outlined in further detail in other chapters of this Science Book. The SKA will be able to produce competitive Baryonic Oscillation (BAOs) measurements in both its phases. The first phase of the SKA will provide similar measurements as optical and IR experiments with completely different systematic effects whereas the second phase being transformational in terms of its statistical power. The SKA will produce very accurate Redshift Space Distortions (RSD) measurements, being superior to other experiments at lower redshifts, due to the large number of galaxies. Cross correlations of the galaxy redshift data from the SKA with radio continuum surveys and optical surveys will provide extremely good calibration of photometric redshifts as well as extremely good bounds on modifications of gravity. Basing on a Principle Component Analysis (PCA) approach, we find that the SKA will be able to provide competitive constraints on dark energy and modified gravity models. Due to the large area covered the SKA it will be a transformational experiment in measuring physics from the largest scales such as non-Gaussian signals from $\textrm{f}_{\textrm{nl}}$. Finally, the SKA might produce the first real time measurement of the redshift drift. The SKA will be a transformational machine for cosmology as it grows from an early Phase 1 to its full power.

The imprint of baryon acoustic oscillations (BAO) in large-scale structure can be used as a standard ruler for mapping out the cosmic expansion history, and hence for testing cosmological models. In this article we briefly describe the scientific background to the BAO technique, and forecast the potential of the Phase 1 and 2 SKA telescopes to perform BAO surveys using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys. We find that a 25,000 sq. deg. intensity mapping survey on a Phase 1 array will preferentially constrain the radial BAO, providing a highly competitive 2% constraint on the expansion rate at z ~ 2. A 30,000 sq. deg. galaxy redshift survey on SKA2 will outperform all other planned experiments for z < 1.4.

In recent years cosmology has undergone a revolution, with precise measurements of the microwave background radiation, large galaxy redshift surveys, and the discovery of the recent accelerated expansion of the Universe using observations of distant supernovae. In this light, the SKA enables us to do an ultimate test in cosmology by measuring the expansion rate of the Universe in real time. This can be done by a rather simple experiment of observing the neutral hydrogen (HI) signal of galaxies at two different epochs. The signal will encounter a change in frequency imprinted as the Universe expands over time and thus monitoring the drift in frequencies will provide a real time measure of the cosmic acceleration. Over a period of 12 years one would expected a frequency shift of the order of 0.1 Hz assuming a standard Lambda-CDM cosmology. Based on the sensitivity estimates of the SKA and the number counts of the expected HI galaxies, it is shown that the number counts are sufficiently high to compensate for the observational uncertainties of the measurements and hence allow a statistical detection of the frequency shift. [abstract abridged]

Radio continuum surveys have, in the past, been of restricted use in cosmology. Most studies have concentrated on cross-correlations with the cosmic microwave background to detect the integrated Sachs-Wolfe effect, due to the large sky areas that can be surveyed. As we move into the SKA era, radio continuum surveys will have sufficient source density and sky area to play a major role in cosmology on the largest scales. In this chapter we summarise the experiments that can be carried out with the SKA as it is built up through the coming decade. We show that the SKA can play a unique role in constraining the non-Gaussianity parameter to \sigma(f_NL) ~ 1, and provide a unique handle on the systematics that inhibit weak lensing surveys. The SKA will also provide the necessary data to test the isotropy of the Universe at redshifts of order unity and thus evaluate the robustness of the cosmological principle.Thus, SKA continuum surveys will turn radio observations into a central probe of cosmological research in the coming decades.

The peculiar motion of galaxies can be a particularly sensitive probe of gravitational collapse. As such, it can be used to measure the dynamics of dark matter and dark energy as well the nature of the gravitational laws at play on cosmological scales. Peculiar motions manifest themselves as an overall anisotropy in the measured clustering signal as a function of the angle to the line-of-sight, known as redshift-space distortion (RSD). Limiting factors in this measurement include our ability to model non-linear galaxy motions on small scales and the complexities of galaxy bias. The anisotropy in the measured clustering pattern in redshift-space is also driven by the unknown distance factors at the redshift in question, the Alcock-Paczynski distortion. This weakens growth rate measurements, but permits an extra geometric probe of the Hubble expansion rate. In this chapter we will briefly describe the scientific background to the RSD technique, and forecast the potential of the SKA phase 1 and the SKA2 to measure the growth rate using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys.

Continuum and HI surveys with the Square Kilometre Array (SKA) will allow us to probe some of the most fundamental assumptions of modern cosmology, including the Cosmological Principle. SKA all-sky surveys will map an enormous slice of space-time and reveal cosmology at superhorizon scales and redshifts of order unity. We illustrate the potential of these surveys and discuss the prospects to measure the cosmic radio dipole at high fidelity. We outline several potentially transformational tests of cosmology to be carried out by means of SKA all-sky surveys.

We investigate optimization strategies to measure primordial non-Gaussianity with future spectroscopic surveys. We forecast measurements coming from the 3D galaxy power spectrum and compute constraints on primordial non-Gaussianity parameters f_NL and n_NG. After studying the dependence on those parameters upon survey specifications such as redshift range, area, number density, we assume a reference mock survey and investigate the trade-off between number density and area surveyed. We then define the observational requirements to reach the detection of f_NL of order 1. Our results show that while power spectrum constraints on non-Gaussianity from future spectroscopic surveys can be competitive with current CMB limits, measurements from higher-order statistics will be needed to reach a sub unity precision in the measurements of the non-Gaussianity parameter f_NL.

We investigate the properties of the 2-point galaxy correlation function at very large scales, including all geometric and local relativistic effects -- wide-angle effects, redshift space distortions, Doppler terms and Sachs-Wolfe type terms in the gravitational potentials. The general three-dimensional correlation function has a nonzero dipole and octupole, in addition to the even multipoles of the flat-sky limit. We study how corrections due to primordial non-Gaussianity and General Relativity affect the multipolar expansion, and we show that they are of similar magnitude (when f_NL is small), so that a relativistic approach is needed. Furthermore, we look at how large-scale corrections depend on the model for the growth rate in the context of modified gravity, and we discuss how a modified growth can affect the non-Gaussian signal in the multipoles.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

If all or a fraction of the dark matter (DM) were coupled to a bath of dark radiation (DR) in the early Universe we expect the combined DM-DR system to give rise to acoustic oscillations of the dark matter until it decouples from the DR. Much like the standard baryon acoustic oscillations, these dark acoustic oscillations (DAO) imprint a characteristic scale, the sound horizon of dark matter, on the matter power spectrum. We compute in detail how the microphysics of the DM-DR interaction affects the clustering of matter in the Universe and show that the DAO physics also gives rise to unique signatures in the temperature and polarization spectra of the cosmic microwave background (CMB). We use cosmological data from the CMB, baryon acoustic oscillations (BAO), and large-scale structure to constrain the possible fraction of interacting DM as well as the strength of its interaction with DR. Like nearly all knowledge we have gleaned about dark matter since inferring its existence this constraint rests on the betrayal by gravity of the location of otherwise invisible dark matter. Although our results can be straightforwardly applied to a broad class of models that couple dark matter particles to various light relativistic species, in order to make quantitative predictions, we model the interacting component as dark atoms coupled to a bath of dark photons. We find that linear cosmological data and CMB lensing put strong constraints on existence of DAO features in the CMB and the large-scale structure of the Universe. Interestingly, we find that at most ~5% of all DM can be very strongly interacting with DR. We show that our results are surprisingly constraining for the recently proposed Double-disk DM model, a novel example of how large-scale precision cosmological data can be used to constrain galactic physics and sub-galactic structure.

The quantity and quality of cosmic structure observations have greatly accelerated in recent years. Further leaps forward will be facilitated by imminent projects, which will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to the nature of dark matter and dark energy. Dark energy and gravity both affect how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure, while the greater the gravity, the more enhanced the growth. While distance measurements also constrain dark energy, the comparison of growth and distance data tests whether General Relativity describes the laws of physics accurately on large scales. Modified gravity models are able to reproduce the distance measurements but at the cost of altering the growth of structure (these signatures are described in more detail in the accompanying paper on Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity. To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000 cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

Forthcoming radio continuum surveys will cover large volumes of the observable Universe and will reach to high redshifts, making them potentially powerful probes of dark energy, modified gravity and non-Gaussianity. Here we extend recent works by analyzing the general relativistic (GR) corrections to the angular power spectrum. These GR corrections to the standard Newtonian analysis of the power spectrum become significant on scales near and beyond the Hubble scale at each redshift. We consider the continuum surveys with LOFAR, WSRT and ASKAP, and examples of continuum surveys with the SKA. We find that the GR corrections are at most percent-level in LOFAR, WODAN and EMU surveys, but they can produce $O(10%)$ changes for high enough sensitivity SKA continuum surveys. The signal is however dominated by cosmic variance, and multiple-tracer observations will be needed to overcome this problem. The GR corrections are suppressed in continuum surveys because of the integration over redshift -- we expect that GR corrections will be enhanced for future SKA HI surveys in which the source redshifts will be known. We also provide predictions for the angular power spectra in the case where the primordial perturbations have local non-Gaussianity. We find that non-Gaussianity dominates over GR corrections, and rises above cosmic variance when $f_{\rm NL}\gtrsim5$ for SKA continuum surveys.

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of more than 3,000 deg^2. We measure the power spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of local non-Gaussianity, f_NL,local, in the primordial density field, paying particular attention to the impact of systematic uncertainties. The BOSS galaxy density field is systematically affected by the local stellar density and this influences the ability to accurately measure f_NL,local. In the absence of any correction, we find (erroneously) that the probability that f_NL,local is greater than zero, P(f_NL,local >0), is 99.5%. After quantifying and correcting for the systematic bias and including the added uncertainty, we find -45 < f_NL,local < 195 at 95% confidence, and P(f_NL,local >0) = 91.0%. A more conservative approach assumes that we have only learned the k-dependence of the systematic bias and allows any amplitude for the systematic correction; we find that the systematic effect is not fully degenerate with that of f_NL,local, and we determine that -82 < f_NL,local < 178 (at 95% confidence) and P(f_NL,local >0) = 68%. This analysis demonstrates the importance of accounting for the impact of Galactic foregrounds on f_NL,local measurements. We outline the methods that account for these systematic biases and uncertainties. We expect our methods to yield robust constraints on f_NL,local for both our own and future large-scale-structure investigations.

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

In this paper, we explore how the forthcoming generation of large-scale radio continuum surveys, with the inclusion of some degree of redshift information, can constrain cosmological parameters. By cross-matching these radio surveys with shallow optical to near-infrared surveys, we can essentially separate the source distribution into a low- and a high-redshift sample, thus providing a constraint on the evolution of cosmological parameters such as those related to dark energy. We examine two radio surveys, the Evolutionary Map of the Universe (EMU) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). A crucial advantage is their combined potential to provide a deep, full-sky survey. The surveys used for the cross-identifications are SkyMapper and SDSS, for the southern and northern skies, respectively. We concentrate on the galaxy clustering angular power spectrum as our benchmark observable, and find that the possibility of including such low redshift information yields major improvements in the determination of cosmological parameters. With this approach, and provided a good knowledge of the galaxy bias evolution, we are able to put strict constraints on the dark energy parameters, i.e. w_0=-0.9+/-0.041 and w_a=-0.24+/-0.13, with type Ia supernovae and CMB priors (with a one-parameter bias in this case); this corresponds to a Figure of Merit (FoM) > 600, which is twice better than what is obtained by using only the cross-identified sources and greater than four time better than the case without any redshift information at all.

We use Luminous Red Galaxies from the Sloan Digital Sky Survey II to test the cosmological structure growth in two alternatives to the standard LCDM+GR cosmological model. We compare observed three-dimensional clustering in SDSS DR7 with theoretical predictions for the standard vanilla LCDM+GR model, Unified Dark Matter cosmologies and the normal branch DGP. In computing the expected correlations in UDM cosmologies, we derive a parameterized formula for the growth factor in these models. For our analysis we apply the methodology tested in Raccanelli et al. 2010 and use the measurements of Samushia et al. 2011, that account for survey geometry, non-linear and wide-angle effects and the distribution of pair orientation. We show that the estimate of the growth rate is potentially degenerate with wide-angle effects, meaning that extremely accurate measurements of the growth rate on large scales will need to take such effects into account. We use measurements of the zeroth and second order moments of the correlation function from SDSS DR7 data and the Large Suite of Dark Matter Simulations, and perform a likelihood analysis to constrain the parameters of the models. Using information on the clustering up to r_max = 120 Mpc/h, and after marginalizing over the bias, we find, for UDM models, a speed of sound < 6.1e-4, and, for the nDGP model, a cross-over scale r_c > 340 Mpc, at 95% confidence level.

The simplest theory describing large-scale redshift-space distortions (RSD), based on linear theory and distant galaxies, depends on the growth of cosmological structure, suggesting that strong tests of General Relativity can be constructed from galaxy surveys. As data sets become larger and the expected constraints more precise, the extent to which the RSD follow the simple theory needs to be assessed in order that we do not introduce systematic errors into the tests by introducing inaccurate simplifying assumptions. We study the impact of the sample geometry, non-linear processes, and biases induced by our lack of understanding of the radial galaxy distribution on RSD measurements. Using LasDamas simulations of the Sloan Digital Sky Survey II (SDSS-II) Luminous Red Galaxy (LRG) data, these effects are shown to be important at the level of 20 per cent. Including them, we can accurately model the recovered clustering in these mock catalogues on scales 30 -- 200 Mpc/h. Applying this analysis to robustly measure parameters describing the growth history of the Universe from the SDSS-II data, gives $f(z=0.25)\sigma_8(z=0.25)=0.3512\pm0.0583$ and $f(z=0.37)\sigma_8(z=0.37)=0.4602\pm0.0378$ when no prior is imposed on the growth-rate, and the background geometry is assumed to follow a $\Lambda$CDM model with the WMAP + SNIa priors. The standard WMAP constrained $\Lambda$CDM model with General Relativity predicts $f(z=0.25)\sigma_8(z=0.25)=0.4260\pm0.0141$ and $f(z=0.37)\sigma_8(z=0.37)=0.4367\pm0.0136$, which is fully consistent with these measurements.

We present forecasts for constraints on cosmological models which can be obtained by forthcoming radio continuum surveys: the wide surveys with the LOw Frequency ARray (LOFAR), Australian Square Kilometre Array Pathfinder (ASKAP) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). We use simulated catalogues appropriate to the planned surveys to predict measurements obtained with the source auto-correlation, the cross-correlation between radio sources and CMB maps (the Integrated Sachs-Wolfe effect), the cross-correlation of radio sources with foreground objects due to cosmic magnification, and a joint analysis together with the CMB power spectrum and supernovae. We show that near future radio surveys will bring complementary measurements to other experiments, probing different cosmological volumes, and having different systematics. Our results show that the unprecedented sky coverage of these surveys combined should provide the most significant measurement yet of the Integrated Sachs-Wolfe effect. In addition, we show that using the ISW effect will significantly tighten constraints on modified gravity parameters, while the best measurements of dark energy models will come from galaxy auto-correlation function analyses. Using the combination of EMU and WODAN to provide a full sky survey, it will be possible to measure the dark energy parameters with an uncertainty of \{$\sigma (w_0) = 0.05$, $\sigma (w_a) = 0.12$\} and the modified gravity parameters \{$\sigma (\eta_0) = 0.10$, $\sigma (\mu_0) = 0.05$\}, assuming Planck CMB+SN(current data) priors. Finally, we show that radio surveys would detect a primordial non-Gaussianity of $f_{\rm NL}$ = 8 at 1-$\sigma$ and we briefly discuss other promising probes.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

We present an analysis of the cross-correlation between the CMB and the large-scale structure (LSS) of the Universe in Unified Dark Matter (UDM) scalar field cosmologies. We work out the predicted cross-correlation function in UDM models, which depends on the speed of sound of the unified component, and compare it with observations from six galaxy catalogues (NVSS, HEAO, 2MASS, and SDSS main galaxies, luminous red galaxies, and quasars). We sample the value of the speed of sound and perform a likelihood analysis, finding that the UDM model is as likely as the LambdaCDM, and is compatible with observations for a range of values of c_\infinity (the value of the sound speed at late times) on which structure formation depends. In particular, we obtain an upper bound of c_\infinity^2 \leq 0.009 at 95% confidence level, meaning that the LambdaCDM model, for which c_\infinity^2 = 0, is a good fit to the data, while the posterior probability distribution peaks at the value c_\infinity^2=10^(-4) . Finally, we study the time dependence of the deviation from LambdaCDM via a tomographic analysis using a mock redshift distribution and we find that the largest deviation is for low-redshift sources, suggesting that future low-z surveys will be best suited to constrain UDM models.

The analysis of Redshift-Space Distortions (RSD) within galaxy surveys provides constraints on the amplitude of peculiar velocities induced by structure growth, thereby allowing tests of General Relativity on extremely large scales. The next generation of galaxy redshift surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), and the Euclid experiment will survey galaxies out to z=2, over 10,000--20,000 sq deg. In such surveys, galaxy pairs with large comoving separation will preferentially have a wide angular separation. In standard plane-parallel theory the displacements of galaxy positions due to RSD are assumed to be parallel for all galaxies, but this assumption will break down for wide-angle pairs. Szapudi 2004 and Papai & Szapudi 2008 provided a methodology, based on tripolar spherical harmonics expansion, for computing the redshift-space correlation function for all angular galaxy pair separations. In this paper we introduce a new procedure for analysing wide-angle effects in numerical simulations. We are able to separate, demonstrate, and fit each of the effects described by the wide-angle RSD theory. Our analysis highlights some of the nuances of dealing with wide-angle pairs, and shows that the effects are not negligible even for relatively small angles. This analysis will help to ensure the full exploitation of future surveys for RSD measurements, which are currently confined to pair separations less than \sim80 Mpc/h out to z\simeq 0.5.

We present a new evolutionary model that describes the population properties of radio sources at frequencies <5 GHz, thus complementing the De Zotti et al. (2005) model, holding at higher frequencies. We find that simple analytic luminosity evolution is still sufficient to fit the wealth of available data on local luminosity functions, multi-frequency source counts, and redshift distributions. However, the fit requires a luminosity-dependent decline of source luminosities at high redshifts, at least for steep-spectrum sources, thus confirming earlier indications of a "downsizing" also for radio sources. The upturn of source counts at sub-mJy levels is accounted for by a straightforward extrapolation, using the empirical far-IR/radio correlation, of evolutionary models matching the far-IR counts and redshift distributions of star-forming galaxies. We also discuss the implications of the new model for the interpretation of data on large-scale clustering of radio sources and on the Integrated Sachs-Wolfe (ISW) effect, and for the investigation of the contribution of discrete sources to the extragalactic background. As for the ISW effect, a new analysis exploiting a very clean CMB map, yields at a substantially higher significance than reported before.

We reassess the estimate of the cross-correlation of the spatial distribution of the NRAO VLA Sky Survey (NVSS) radio sources with that of Cosmic Microwave Background (CMB) anisotropies from the Wilkinson Microwave Anisotropy Probe (WMAP). This re-analysis is motivated by the fact that most previous studies adopted a redshift distribution of NVSS sources inconsistent with recent data. We find that the constraints on the bias-weighted redshift distribution, b(z)xN(z), of NVSS sources, set by the observed angular correlation function, w(theta), strongly mitigate the effect of the choice of N(z). If such constraints are met, even highly discrepant redshift distributions yield NVSS-WMAP cross-correlation functions consistent with each other within statistical errors. The models favoured by recent data imply a bias factor, b(z), decreasing with increasing z, rather than constant, as assumed by most previous analyses. As a consequence, the function b(z)xN(z) has more weight at z<1, i.e. in the redshift range yielding the maximum contribution to the ISW in a standard LambdaCDM cosmology. On the whole, the NVSS turns out to be better suited for ISW studies than generally believed, even in the absence of an observational determination of the redshift distribution. The NVSS-WMAP cross-correlation function is found to be fully consistent with the prediction of the standard LambdaCDM cosmology.

We describe the BOOMERANG North America (BNA) instrument, a balloon-borne bolometric radiometer designed to map the Cosmic Microwave Background (CMB) radiation with 0.3 deg resolution over a significant portion of the sky. This receiver employs new technologies in bolometers, readout electronics, millimeter-wave optics and filters, cryogenics, scan and attitude reconstruction. All these subsystems are described in detail in this paper. The system has been fully calibrated in flight using a variety of techniques which are described and compared. It has been able to obtain a measurement of the first peak in the CMB angular power spectrum in a single balloon flight, few hours long, and was a prototype of the BOOMERANG Long Duration Balloon (BLDB) experiment.