In recent years forecasting activities have become a very important tool for designing and optimising large scale structure surveys. To predict the performance of such surveys, the Fisher matrix formalism is frequently used as a fast and easy way to compute constraints on cosmological parameters. Among them lies the study of the properties of dark energy which is one of the main goals in modern cosmology. As so, a metric for the power of a survey to constrain dark energy is provided by the Figure of merit (FoM). This is defined as the inverse of the surface contour given by the joint variance of the dark energy equation of state parameters $\{w_0,w_a\}$ in the Chevallier-Polarski-Linder parameterisation, which can be evaluated from the covariance matrix of the parameters. This covariance matrix is obtained as the inverse of the Fisher matrix. Inversion of an ill-conditioned matrix can result in large errors on the covariance coefficients if the elements of the Fisher matrix have been estimated with insufficient precision. The conditioning number is a metric providing a mathematical lower limit to the required precision for a reliable inversion, but it is often too stringent in practice for Fisher matrices with size larger than $2\times2$. In this paper we propose a general numerical method to guarantee a certain precision on the inferred constraints, like the FoM. It consists on randomly vibrating (perturbing) the Fisher matrix elements with Gaussian perturbations of a given amplitude, and then evaluating the maximum amplitude that keeps the FoM within the chosen precision. The steps used in the numerical derivatives and integrals involved in the calculation of the Fisher matrix elements can then be chosen accordingly in order to keep the precision of the Fisher matrix elements below this maximum amplitude...

We present a framework and an open-source python toolkit to analyse the 2-point statistics of 3D fluctuations in the context of HI intensity maps using the multipole expansion formalism. We include simulations of the cosmological HI signal using N-body and log-normal methods, foregrounds and their removal, as well as instrumental effects. Using these simulations and analytical modelling, we investigate the impact of foreground cleaning and the instrumental beam on the power spectrum multipoles as well as on the Fourier space clustering wedges. We find that both the instrumental beam and the foreground removal can produce a quadrupole (and a hexadecapole) signal, and demonstrate the importance of controlling and accurately modelling these effects for precision radio cosmology. We conclude that these effects can be modelled with reasonable accuracy using our multipole expansion technique. We also perform an MCMC analysis to showcase the effect of foreground cleaning on the estimation of the HI abundance and bias parameters. The accompanying python toolkit is available at https://github.com/IntensityTools/MultipoleExpansion, and includes an interactive suite of examples to aid new users.

We investigate the prospects of measuring the cosmic magnification effect by cross-correlating neutral hydrogen intensity mapping (HI IM) maps with background optical galaxies. We forecast the signal-to-noise ratio for HI IM data from SKA1-MID and HIRAX, combined with LSST photometric galaxy samples. We find that, thanks to their different resolutions, SKA1-MID and HIRAX are highly complementary in such an analysis. We predict that SKA1-MID can achieve a detection with a signal-to-noise ratio of $\sim 15$ on a multipole range of $\ell \lesssim 200$, while HIRAX can reach a signal-to-noise ratio of $\sim 50$ on $200 < \ell < 2000$. We conclude that measurements of the cosmic magnification signal will be possible on a wide redshift range with foreground HI intensity maps up to $z \lesssim 2$, while optimal results are obtained when $0.6 \lesssim z \lesssim 1.3$. Finally, we perform a signal to noise analysis that shows how these measurements can constrain the HI parameters across a wide redshift range.

Accurate modelling of non-linear scales in galaxy clustering will be crucial for data analysis of Stage IV galaxy surveys. A selection of competing non-linear models must be made based on validation studies. We provide a comprehensive set of forecasts of two different models for the halo redshift space power spectrum, namely the commonly applied TNS model and an effective field theory of large scale structure (EFTofLSS) inspired model. Using simulation data and a least-$\chi^2$ analysis, we determine ranges of validity for the models. We then conduct an exploratory Fisher analysis using the full anisotropic power spectrum to investigate parameter degeneracies. We proceed to perform an MCMC analysis utilising the monopole, quadrupole, and hexadecapole spectra, with a restricted range of scales for the latter in order to avoid biasing our growth rate, $f$, constraint. We find that the TNS model with a Lorentzian damping and standard Eulerian perturbative modelling outperforms other variants of the TNS model. Our MCMC analysis finds that the EFTofLSS-based model may provide tighter marginalised constraints on $f$ at $z=0.5$ and $z=1$ than the TNS model, despite having additional nuisance parameters. However this depends on the range of scales used as well as the fiducial values and priors on the EFT nuisance parameters. Finally, we extend previous work to provide a consistent comparison between the Fisher matrix and MCMC forecasts using the multipole expansion formalism, and find good agreement between them.

Models in which scalar field dark energy interacts with dark matter via a pure momentum coupling have previously been found to potentially ease the structure formation tension between early- and late-universe observations. In this article we explore the physical mechanism underlying this feature. We argue analytically that the perturbation growth equations imply the suppression of structure growth, illustrating our discussion with numerical calculations. Then we generalise the previously studied quadratic coupling between the dark energy and dark matter to a more general power law case, also allowing for the slope of the dark energy exponential potential to vary. We find that the structure growth suppression is a generic feature of power law couplings and it can, for a range of parameter values, be larger than previously found.

Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian streaming model to calculate the quasi non-linear redshift space matter power spectrum multipoles from a non-linear prescription for the matter power spectrum, which remains agnostic about gravity or dark energy. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model. When comparing the redshift space multipoles for halos, we find that a linear bias implementation is insufficient to be competitive with the TNS-based model used in the BOSS clustering analysis, which however has additional degrees of freedom. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model.

We assess the performance of the Taruya, Nishimichi and Saito (TNS) model for the halo redshift space power spectrum, focusing on utilising mildly non-linear scales to constrain the growth rate of structure f. Using simulations with volume and number density typical of forthcoming Stage IV galaxy surveys, we determine ranges of validity for the model at redshifts z = 0.5 and z = 1. We proceed to perform a Bayesian MCMC analysis utilising the monopole, quadrupole, and hexadecapole spectra, followed by an exploratory Fisher matrix analysis. As previously noted in other forecasts as well as in real data analyses, we find that including the hexadecapole can significantly improve the constraints. However, a restricted range of scales is required for the hexadecapole in order for the growth parameter estimation to remain unbiased, limiting the improvement. We consistently quantify these effects by employing the multipole expansion formalism in both our Fisher and MCMC forecasts.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

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.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

The future of precision cosmology could benefit from cross-correlations between intensity maps of unresolved neutral hydrogen (HI) and more conventional optical galaxy surveys. A major challenge that needs to be overcome is removing the 21cm foreground emission that contaminates the cosmological HI signal. Using N-body simulations we simulate HI intensity maps and optical catalogues which share the same underlying cosmology. Adding simulated foreground contamination and using state-of-the-art reconstruction techniques we investigate the impacts that 21cm foregrounds and other systematics have on these cross-correlations. We find that the impact a FASTICA 21cm foreground clean has on the cross-correlations with spectroscopic optical surveys with well-constrained redshifts is minimal. However, problems arise when photometric surveys are considered: we find that a redshift uncertainty {\sigma}_z {\geq} 0.04 causes significant degradation in the cross power spectrum signal. We diagnose the main root of these problems, which relates to arbitrary amplitude changes along the line-of-sight in the intensity maps caused by the foreground clean and suggest solutions which should be applicable to real data. These solutions involve a reconstruction of the line-of-sight temperature means using the available overlapping optical data along with an artificial extension to the HI data through redshift to address edge effects. We then put these solutions through a further test in a mock experiment that uses a clustering-based redshift estimation technique to constrain the photometric redshifts of the optical sample. We find that with our suggested reconstruction, cross-correlations can be utilized to make an accurate prediction of the optical redshift distribution.

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.

Precision cosmology requires accurate galaxy redshifts, but next generation optical surveys will observe unprecedented numbers of resolved galaxies, placing strain on the amount of spectroscopic follow-up required. We show how useful information can be gained on the redshift distributions of optical galaxy samples from spatial cross-correlations with intensity maps of unresolved HI (21cm) spectral line emission. We construct a redshift distribution estimator, which we test using simulations. We utilise the S3-SAX catalogue which includes HI emission information for each galaxy, which we use to construct HI intensity maps. We also make use of simulated LSST and Euclid-like photometry enabling us to apply the HI clustering calibration to realistic simulated photometric redshifts. While taking into account important limitations to HI intensity mapping such as lost k-modes from foreground cleaning and poor angular resolution due to large receiver beams, we show that excellent constraints on redshift distributions can be provided for an optical photometric sample.

We investigate the capabilities of perturbation theory in capturing non-linear effects of dark energy. We test constant and evolving $\omega$ models, as well as models involving momentum exchange between dark energy and dark matter. Specifically, we compare perturbative predictions at 1-loop level against N-body results for four non-standard equations of state as well as varying degrees of momentum exchange between dark energy and dark matter. The interaction is modelled phenomenologically using a time dependent drag term in the Euler equation. We make comparisons at the level of the matter power spectrum and the redshift space monopole and quadrupole. The multipoles are modelled using the Taruya, Nishimichi and Saito (TNS) redshift space spectrum. We find perturbation theory does very well in capturing non-linear effects coming from dark sector interaction. We isolate and quantify the 1-loop contribution coming from the interaction and from the non-standard equation of state. We find the interaction parameter $\xi$ amplifies scale dependent signatures in the range of scales considered. Non-standard equations of state also give scale dependent signatures within this same regime. In redshift space the match with N-body is improved at smaller scales by the addition of the TNS free parameter $\sigma_v$. To quantify the importance of modelling the interaction, we create mock data sets for varying values of $\xi$ using perturbation theory. This data is given errors typical of Stage IV surveys. We then perform a likelihood analysis using the first two multipoles on these sets and a $\xi=0$ modelling, ignoring the interaction. We find the fiducial growth parameter $f$ is generally recovered even for very large values of $\xi$ both at $z=0.5$ and $z=1$. The $\xi=0$ modelling is most biased in its estimation of $f$ for the phantom $\omega=-1.1$ case.

We consider a self-consistent and physical approach to interacting dark energy models described by a Lagrangian, and identify a new class of models with variable dark energy sound speed. We show that if the interaction between dark energy in the form of quintessence and cold dark matter is purely momentum exchange this generally leads to a dark energy sound speed that deviates from unity. Choosing a specific sub-case, we study its phenomenology by investigating the effects of the interaction on the cosmic microwave background and linear matter power spectrum. We also perform a global fitting of cosmological parameters using CMB data, and compare our findings to $\Lambda$CDM.

We introduce simulations aimed at assessing how well weak gravitational lensing of 21cm radiation from the Epoch of Reionization ($z \sim 8$) can be measured by an SKA-like radio telescope. A simulation pipeline has been implemented to study the performance of lensing reconstruction techniques. We show how well the lensing signal can be reconstructed using the three-dimensional quadratic lensing estimator in Fourier space assuming different survey strategies. The numerical code introduced in this work is capable of dealing with issues that can not be treated analytically such as the discreteness of visibility measurements and the inclusion of a realistic model for the antennae distribution. This paves the way for future numerical studies implementing more realistic reionization models, foreground subtraction schemes, and testing the performance of lensing estimators that take into account the non-Gaussian distribution of HI after reionization. If multiple frequency channels covering $z \sim 7-11.6$ are combined, Phase 1 of SKA-Low should be able to obtain good quality images of the lensing potential with a total resolution of $\sim 1.6$ arcmin. The SKA-Low Phase 2 should be capable of providing images with high-fidelity even using data from $z\sim 7.7 - 8.3$. We perform tests aimed at evaluating the numerical implementation of the mapping reconstruction. We also discuss the possibility of measuring an accurate lensing power spectrum. Combining data from $z \sim 7-11.6$ using the SKA2-Low telescope model, we find constraints comparable to sample variance in the range $L<1000$, even for survey areas as small as $25\mbox{ deg}^2$.

We explore the possibility of performing an HI intensity mapping survey with the South African MeerKAT radio telescope, which is a precursor to the Square Kilometre Array (SKA). We propose to use cross-correlations between the MeerKAT intensity mapping survey and optical galaxy surveys, in order to mitigate systematic effects and produce robust cosmological measurements. Our forecasts show that precise measurements of the HI signal can be made in the near future. These can be used to constrain HI and cosmological parameters across a wide range of redshift.

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 forecast constraints on neutral hydrogen (HI) and cosmological parameters using near-term intensity mapping surveys with instruments such as BINGO, MeerKAT, and the SKA, and Stage III and IV optical galaxy surveys. If foregrounds and systematic effects can be controlled - a problem which becomes much easier in cross-correlation - these surveys will provide exquisite measurements of the HI density and bias, as well as measurements of the growth of structure, the angular diameter distance, and the Hubble rate, over a wide range of redshift. We also investigate the possibility of detecting the late time ISW effect using the Planck satellite and forthcoming intensity mapping surveys, finding that a large sky survey with Phase 1 of the SKA can achieve a near optimal detection.

We investigate the possibility of utilising 21cm intensity mapping, optical galaxy, and Cosmic Microwave Background (CMB) surveys to constrain the power spectrum of primordial fluctuations predicted by single-field slow-roll inflation models. Implementing a Fisher forecast analysis, we derive constraints on the spectral tilt parameter $n_{\rm s}$ and its first and second runnings $(\alpha_{\rm s},\beta_{\rm s})$. We show that 21cm intensity mapping surveys with instruments like the Square Kilometre Array, CHIME, and HIRAX, can be powerful probes of the primordial features. We combine our forecasts with the ones derived for a COrE-like CMB survey, as well as for a Stage IV optical galaxy survey similar to Euclid. The synergies between different surveys can be exploited to rule out a large fraction of the available inflationary models.

21cm intensity mapping is a novel approach aimed at measuring the power spectrum of density fluctuations and deducing cosmological information, notably from the Baryonic Acoustic Oscillations (BAO). We give an update on the progress of BAO from Integrated Neutral Gas Observations (BINGO) which is a single dish intensity mapping project. First we explain the basic ideas behind intensity mapping concept before updating the instrument design for BINGO. We also outline the survey we plan to make and its projected science output including estimates of cosmological parameters.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We study a coupled quintessence model with pure momentum exchange and present the effects of such an interaction on the Cosmic Microwave Background (CMB) and matter power spectrum. For a wide range of negative values of the coupling parameter $\beta$ structure growth is suppressed and the model can reconcile the tension between Cosmic Microwave Background observations and structure growth inferred from cluster counts. We find that this model is as good as $\Lambda$CDM for CMB and baryon acoustic oscillation (BAO) data, while the addition of cluster data makes the model strongly preferred, improving the best-fit $\chi^2$-value by more than $16$.

We investigate the possibility of testing Einstein's general theory of relativity (GR) and the standard cosmological model via the $E_{\rm G}$ statistic using neutral hydrogen (HI) intensity mapping. We generalise the Fourier space estimator for $E_{\rm G}$ to include HI as a biased tracer of matter and forecast statistical errors using HI clustering and lensing surveys that can be performed in the near future, in combination with ongoing and forthcoming optical galaxy and Cosmic Microwave Background (CMB) surveys. We find that fractional errors $< 1\%$ in the $E_{\rm G}$ measurement can be achieved in a number of cases and compare the ability of various survey combinations to differentiate between GR and specific modified gravity models. Measuring $E_{\rm G}$ with intensity mapping and the Square Kilometre Array can provide exquisite tests of gravity at cosmological scales.

We explore the potential of using intensity mapping surveys (MeerKAT, SKA) and optical galaxy surveys (DES, LSST) to detect HI clustering and weak gravitational lensing of 21cm emission in auto- and cross-correlation. Our forecasts show that high precision measurements of the clustering and lensing signals can be made in the near future using the intensity mapping technique. Such studies can be used to test the intensity mapping method, and constrain parameters such as the HI density $\Omega_{\rm HI}$, the HI bias $b_{\rm HI}$ and the galaxy-HI correlation coefficient $r_{\rm HI-g}$.

The cross-correlation of a foreground density field with two different background convergence fields can be used to measure cosmographic distance ratios and constrain dark energy parameters. We investigate the possibility of performing such measurements using a combination of optical galaxy surveys and HI intensity mapping surveys, with emphasis on the performance of the planned Square Kilometre Array (SKA). Using HI intensity mapping to probe the foreground density tracer field and/or the background source fields has the advantage of excellent redshift resolution and a longer lever arm achieved by using the lensing signal from high redshift background sources. Our results show that, for our best SKA-optical configuration of surveys, a constant equation of state for dark energy can be constrained to $\simeq 8\%$ for a sky coverage $f_{\rm sky}=0.5$ and assuming a $\sigma(\Omega_{\rm DE})=0.03$ prior for the dark energy density parameter. We also show that using the CMB as the second source plane is not competitive, even when considering a COrE-like satellite.

We present the most general parametrisation of models of dark energy in the form of a scalar field which is explicitly coupled to dark matter. We follow and extend the Parameterized Post-Friedmannian approach, previously applied to modified gravity theories, in order to include interacting dark energy. We demonstrate its use through a number of worked examples and show how the initially large parameter space of free functions can be significantly reduced and constrained to include only a few non-zero coefficients. This paves the way for a model-independent approach to classify and test interacting dark energy theories.

We investigate the feasibility of measuring weak gravitational lensing using 21cm intensity mapping with special emphasis on the performance of the planned Square Kilometre Array (SKA). We find that the current design for SKA-Mid should be able to measure the evolution of the lensing power spectrum at z~2-3 using this technique. This will be a probe of the expansion history of the universe and gravity at a unique range in redshift. The signal-to-noise is found to be highly dependent on evolution of the neutral hydrogen fraction in the universe with a higher HI density resulting in stronger signal. With realistic models for this, SKA Phase 1 should be capable of measuring the lensing power spectrum and its evolution. The signal-to-noise's dependence on the area and diameter of the telescope array is quantified. We further demonstrate the applications of this technique by applying it to two specific coupled dark energy models that would be difficult to observationally distinguish without information from this range of redshift. We also investigate measuring the lensing signal with 21cm emission from the Epoch of Reionization (EoR) using SKA-Low and find that it is unlikely to constrain cosmological parameters because of the small survey size, but could provide a map of the dark matter within a small region of the sky.

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

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.

We investigate the capabilities of various stages of the SKA to perform world-leading weak gravitational lensing surveys. We outline a way forward to develop the tools needed for pursuing weak lensing in the radio band. We identify the key analysis challenges and the key pathfinder experiments that will allow us to address them in the run up to the SKA. We identify and summarize the unique and potentially very powerful aspects of radio weak lensing surveys, facilitated by the SKA, that can solve major challenges in the field of weak lensing. These include the use of polarization and rotational velocity information to control intrinsic alignments, and the new area of weak lensing using intensity mapping experiments. We show how the SKA lensing surveys will both complement and enhance corresponding efforts in the optical wavebands through cross-correlation techniques and by way of extending the reach of weak lensing to high redshift.

We study how 21 cm intensity mapping can be used to measure gravitational lensing over a wide range of redshift. This can extend weak lensing measurements to higher redshifts than are accessible with conventional galaxy surveys. We construct a convergence estimator taking into account the discreteness of galaxies and calculate the expected noise level as a function of redshift and telescope parameters. At $z \sim 2-3$ we find that a telescope array with a collecting area $\sim 0.2 \, {\rm km}^2$ spread over a region with diameter $\sim 2 \, {\rm km}$ would be sufficient to measure the convergence power spectrum to high accuracy for multipoles between 10 and 1,000. We show that these measurements can be used to constrain interacting dark energy models.

We present three distinct types of models of dark energy in the form of a scalar field which is explicitly coupled to dark matter. Our construction draws from the pull-back formalism for fluids and generalises the fluid action to involve couplings to the scalar field. We investigate the cosmology of each class of model both at the background and linearly perturbed level. We choose a potential for the scalar field and a specific coupling function for each class of models and we compute the Cosmic Microwave Background and matter power spectra.

BINGO is a concept for performing a 21cm intensity mapping survey using a single dish telescope. We briefly discuss the idea of intensity mapping and go on to define our single dish concept. This involves a \sim 40 m dish with an array of \sim 50 feed horns placed \sim 90 m above the dish using a pseudo-correlation detection system based on room temperature LNAs and one of the celestial poles as references. We discuss how such an array operating between 960 and 1260 MHz could be used to measure the acoustic scale to 2.4% over the redshift range 0.13<z<0.48 in around 1 year of on-source integration time by performing a 10 deg \times 200 deg drift scan survey with a resolution of \sim 2/3 deg.

We discuss the detection of large scale HI intensity fluctuations using a single dish approach with the ultimate objective of measuring the Baryonic Acoustic Oscillations and constraining the properties of dark energy. We present 3D power spectra, 2D angular power spectra for individual redshift slices, and also individual line-of-sight spectra, computed using the S^3 simulated HI catalogue which is based on the Millennium Simulation. We consider optimal instrument design and survey strategies for a single dish observation at low and high redshift for a fixed sensitivity. For a survey corresponding to an instrument with T_sys=50 K, 50 feed horns and 1 year of observations, we find that at low redshift (z \approx 0.3), a resolution of 40 arc min and a survey of 5000 deg^2 is close to optimal, whereas at higher redshift (z \approx 0.9) a resolution of 10 arcmin and 500 deg^2 would be necessary. Continuum foreground emission from the Galaxy and extragalactic radio sources are potentially a problem. We suggest that it could be that the dominant extragalactic foreground comes from the clustering of very weak sources. We assess its amplitude and discuss ways by which it might be mitigated. We then introduce our concept for a single dish telescope designed to detect BAO at low redshifts. It involves an under-illumintated static 40 m dish and a 60 element receiver array held 90 m above the under-illuminated dish. Correlation receivers will be used with each main science beam referenced against an antenna pointing at one of the Celestial Poles for stability and control of systematics. We make sensitivity estimates for our proposed system and projections for the uncertainties on the power spectrum after 1 year of observations. We find that it is possible to measure the acoustic scale at z\approx 0.3 with an accuracy 2.4% and that w can be measured to an accuracy of 16%.

We study signatures of cosmic superstring networks containing strings of multiple tensions and Y-junctions, on the cosmic microwave background (CMB) temperature and polarisation spectra. Focusing on the crucial role of the string coupling constant $g_s$, we show that the number density and energy density of the scaling network are dominated by different types of string in the $g_s \sim 1$ and $g_s \ll 1$ limits. This can lead to an observable shift in the position of the B-mode peak --- a distinct signal leading to a direct constraint on $g_s$. We forecast the joint bounds on $g_s$ and the fundamental string tension $\mu_F$ from upcoming and future CMB polarisation experiments, as well as the signal to noise in detecting the difference between B-mode signals in the limiting cases of large and small $g_s$. We show that such a detectable shift is within reach of planned experiments.

We study the cosmic microwave background temperature and polarisation spectra sourced by multi-tension cosmic superstring networks. First we obtain solutions for the characteristic length scales and velocities associated with the evolution of a network of F-D strings, allowing for the formation of junctions between strings of different tensions. We find two distinct regimes describing the resulting scaling distributions for the relative densities of the different types of strings, depending on the magnitude of the fundamental string coupling g_s. In one of them, corresponding to the value of the coupling being of order unity, the network's stress-energy power spectrum is dominated by populous light F and D strings, while the other regime, at smaller values of g_s, has the spectrum dominated by rare heavy D strings. These regimes are seen in the CMB anisotropies associated with the network. We focus on the dependence of the shape of the B-mode polarisation spectrum on g_s and show that measuring the peak position of the B-mode spectrum can point to a particular value of the string coupling. Finally, we assess how this result, along with pulsar bounds on the production of gravitational waves from strings, can be used to constrain a combination of g_s and the fundamental string tension mu_F. Since CMB and pulsar bounds constrain different combinations of the string tensions and densities, they result in distinct shapes of bounding contours in the (mu_F, g_s) parameter plane, thus providing complementary constraints on the properties of cosmic superstrings.

We consider the asymmetric branes model of modified gravity, which can produce late time acceleration of the universe and compare the cosmology of this model to the standard $\Lambda$CDM model and to the DGP braneworld model. We show how the asymmetric cosmology at relevant physical scales can be regarded as a one-parameter extension of the DGP model, and investigate the effect of this additional parameter on the expansion history of the universe.