The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

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

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The Fourier-space galaxy bispectrum is complex, with the imaginary part arising from leading-order relativistic corrections, due to Doppler, gravitational redshift and related line-of-sight effects in redshift space. The detection of the imaginary part of the bispectrum is potentially a smoking gun signal of relativistic contributions. We investigate whether next-generation spectroscopic surveys could make such a detection. For a Stage IV spectroscopic $H\alpha$ survey similar to Euclid, we find that the cumulative signal to noise of this relativistic signature is $\mathcal{O}(10)$. Long-mode relativistic effects couple to short-mode Newtonian effects in the galaxy bispectrum, but not in the galaxy power spectrum. This is the basis for detectability of relativistic effects in the bispectrum of a single galaxy survey, whereas the power spectrum requires multiple galaxy surveys to detect the corresponding signal.

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

In recent years, many gamma-ray sources have been identified, yet the unresolved component hosts valuable information on the faintest emission. In order to extract it, a cross-correlation with gravitational tracers of matter in the Universe has been shown to be a promising tool. We report here the first identification of a cross-correlation signal between gamma rays and the distribution of mass in the Universe probed by weak gravitational lensing. We use the Dark Energy Survey Y1 weak lensing catalogue and the Fermi Large Area Telescope 9-year gamma-ray data, obtaining a signal-to-noise ratio of 5.3. The signal is mostly localised at small angular scales and high gamma-ray energies, with a hint of correlation at extended separation. Blazar emission is likely the origin of the small-scale effect. We investigate implications of the large-scale component in terms of astrophysical sources and particle dark matter emission.

We report the identification of a positive cross-correlation signal between the unresolved $\gamma$-ray emission, measured by the \emph{Fermi} Large Area Telescope, and four different galaxy cluster catalogues. The selected catalogues peak at low-redshift and span different frequency bands, including infrared, optical and X-rays. The signal-to-noise ratio of the detected cross-correlation amounts to 3.5 in the most significant case. We investigate and comment about its possible origin, in terms of compact $\gamma$-ray emission from AGNs inside clusters or diffuse emission from the intracluster medium. The analysis has been performed by introducing an accurate estimation of the cross-correlation power-spectrum covariance matrix, built with mock realisations of the gamma and galaxy cluster maps. Different methods to produce the mock realizations starting from the data maps have been investigated and compared, identifying suitable techniques which can be generalized to other cross-correlation studies.

A new generation of cosmological experiments will spectroscopically detect the H$\alpha$ line from emission-line galaxies (ELGs) at optical/near-infrared frequencies. Other emission lines will also be present, which may come from the same H$\alpha$ sample or constitute a new galaxy sample altogether. Our goal is to assess the value, for cosmological investigation, of galaxies at $z\gtrsim2$ present in H$\alpha$ galaxy surveys and identifiable by the highly redshifted ultra-violet and optical lines - namely the OII line and the OIII doublet in combination with the H$\beta$ line. We use state-of-the-art models of luminosity functions of astrophysical spectral lines to estimate the volumetric number density of OIII+H$\beta$ and OII ELGs. We focus on a wavelength range which will be covered by planned cosmological surveys. We study the constraining power of these high-redshift galaxy samples on cosmological parameters such as the amplitude of baryon acoustic oscillations, $H(z)$, $D_A(z)$, $f\sigma_8(z)$, and $b\sigma_8(z)$ for different survey designs. We present a strong science case for extracting the OIII+H$\beta$ sample, which we consider as an independent probe of the Universe in the redshift range 2 to 3. Moreover, we show that the OII sample can be used to measure the baryon acoustic oscillations and the growth of structures above $z=3$; albeit it may be shot-noise dominated, it will nonetheless provide valuable tomographic information. Summarising, we discuss the scientific potential of a sample of galaxies which, so far, has been mainly considered as a contaminant in H$\alpha$ galaxy surveys. Our findings indicate that planed H$\alpha$ surveys should include the extraction of these oxygen-line samples in their pipeline, to enhance their scientific impact on cosmology.

Tighter constraints on measurements of primordial non-Gaussianity will allow the differentiation of inflationary scenarios. The cosmic microwave background bispectrum-the standard method of measuring the local non-Gaussianity-is limited by cosmic variance. Therefore, it is sensible to investigate measurements of non-Gaussianity using the large-scale structure. This can be done by investigating the effects of non-Gaussianity on the power spectrum on large scales. In this study we forecast the constraints on the local primordial non-Gaussianity parameter $f_{\rm NL}$ that can be obtained with future radio surveys. We utilize the multi-tracer method which reduces the effect of cosmic variance and takes advantage of the multiple radio galaxy populations which are differently biased tracers of the same underlying dark matter distribution. Improvements on previous work include the use of observational bias and halo mass estimates, updated simulations and realistic photometric redshift expectations, thus producing more realistic forecasts. Combinations of SKA simulations and radio observations were used as well as different redshift ranges and redshift bin sizes. It was found that in the most realistic case the 1 - $\sigma$ error on $f_{\rm NL}$ falls within the range 4.07 and 6.58, rivalling the tightest constraints currently available.

Neutral hydrogen (HI) intensity mapping traces the large-scale distribution of matter in the Universe and therefore should correlate with the gamma-ray emission originated from particle dark matter annihilation or from active galactic nuclei and star-forming galaxies, since the related processes occur in the same cosmic structures hosting HI. In this paper, we derive the cross-correlation signal between the brightness temperature of the 21-cm line emission of the HI spin-flip transition in the Universe and the unresolved gamma-ray background. Specifically, we derive forecasts for the cross-correlation signal by focussing on the opportunities offered by the combination of the Fermi-Large Area Telescope (LAT) gamma-ray sensitivity with the expectations of the HI intensity mapping measurements from future radio telescopes, for which we concentrate on the Square Kilometre Array (SKA) and MeerKAT, one of its precursors. We find that the combination of MeerKAT with the current Fermi-LAT statistics has the potential to provide a first hint of the cross-correlation signal originated by astrophysical sources, with a signal-to-noise ratio (SNR) of 3.7. With SKA Phase 1 and SKA Phase 2, the SNR is predicted to increase up to 5.7 and 8.2, respectively. The bounds on dark matter properties attainable with SKA combined with the current statistics of Fermi-LAT are predicted to be comparable to those obtained from other techniques able to explore the unresolved components of the gamma-ray background. The enhanced capabilities of SKA Phase 2, combined with a future generation gamma-ray telescope with improved specifications, can allow us to investigate the whole mass window for weakly interacting massive particles up to the TeV scale.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

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.

Mapping the unresolved intensity of 21-cm emission of neutral hydrogen (HI) is now regarded as one the most promising tools for cosmological investigation in the coming decades. Here, we investigate, for the first time, extensions of the standard cosmological model, such as modified gravity and primordial non-Gaussianity, taking the present constraints on the astrophysics of HI clustering in the treatment of the brightness temperature fluctuations self-consistently into account. To understand the boundaries within which results thus obtained can be considered reliable, we examine the robustness of cosmological parameter estimation performed via studies of 21-cm intensity mapping, against our knowledge of the astrophysical processes leading to HI clustering. Modelling of astrophysical effects affects cosmological observables through the relation linking the overall HI mass in a bound object, to the mass of the underlying dark matter halo that hosts it. We quantify the biases in standard cosmological parameters and those describing modified gravity and primordial non-Gaussianity, that are obtained if one misconceives the slope of the relation between HI mass and halo mass, or the lower virial velocity cut-off for a dark matter halo to be able to host HI. Remarkably, we find that astrophysical uncertainties will not affect searches for primordial non-Gaussianity - one of the strongest science cases for HI intensity mapping - despite the signal being deeply linked to the HI bias.

We develop a cosmological parameter estimation code for (tomographic) angular power spectra analyses of galaxy number counts, for which we include, for the first time, redshift-space distortions (RSD) in the Limber approximation. This allows for a speed-up in computation time, and we emphasise that only angular scales where the Limber approximation is valid are included in our analysis. Our main result shows that a correct modelling of RSD is crucial not to bias cosmological parameter estimation. This happens not only for spectroscopy-detected galaxies, but even in the case of galaxy surveys with photometric redshift estimates. Moreover, a correct implementation of RSD is especially valuable in alleviating the degeneracy between the amplitude of the underlying matter power spectrum and the galaxy bias. We argue that our findings are particularly relevant for present and planned observational campaigns, such as the Euclid satellite or the Square Kilometre Array, which aim at studying the cosmic large-scale structure and trace its growth over a wide range of redshifts and scales.

Following on our purpose of developing a unified pipeline for large-scale structure data analysis with angular (i.e. harmonic-space) power spectra, we now include the weak lensing effect of magnification bias on galaxy clustering in a publicly available, modular parameter estimation code. We thus forecast constraints on the parameters of the concordance cosmological model, dark energy, and modified gravity theories from galaxy clustering tomographic angular power spectra. We find that a correct modelling of magnification is crucial in order not to bias the estimation of cosmological parameters, especially in the case of deep galaxy surveys. Our case study adopts specifications of the Evolutionary Map of the Universe (EMU), which is a full-sky, deep radio-continuum survey, and is expected to probe the Universe up to redshift $z\sim6$. We assume the Limber approximation, and include magnification bias on top of density fluctuations and redshift-space distortions. By restricting our analysis to the regime where the Limber approximation holds true, we significantly minimise the computational time needed, compared to that of the exact calculation. We also show that there is a trend for more biased parameter estimates from neglecting magnification when the redshift bins are very wide. We conclude that this result implies a strong dependence on the lensing contribution, which is an integrated effect and becomes dominant when wide redshift bins are considered. Finally, we note that instead of being considered a contaminant, magnification bias encodes important cosmological information, and its inclusion leads leading to an alleviation of its degeneracy between the galaxy bias and the amplitude normalisation of the matter fluctuations.

The angular power spectrum is a gauge-independent observable that is in principle the natural tool for analysing galaxy number counts. In practice, the problem is that the computational requirements for next-generation spectroscopic surveys such as Euclid and the Square Kilometre Array are currently unfeasible. We propose a new method to save computational time for spectroscopic angular power spectra. This hybrid method is modelled on the Fourier power spectrum approach of treating relatively thick redshift bins (redshift width ~0.1) as separate surveys. In the hybrid method, each thick bin is further subdivided into thin bins (redshift width ~0.01); all the correlations within each thick bin are computed, while cross-bin correlations beyond the thick bins are neglected. Constraints on cosmological parameters from the hybrid method are comparable to those from the standard galaxy power spectrum analysis - but they have the advantage that cosmic evolution, wide-angle and lensing effects are naturally included, while no Alcock-Paczynski correction is needed. The hybrid method delivers much tighter constraints than a 2D tomographic approach that is typical for photometric surveys, which considers only thick bins and the correlations between them. Furthermore, for standard cosmological parameters our method is not biased by neglecting the effects of lensing on number counts, while the tomographic method is strongly biased.

Tensions between cosmic microwave background observations and the growth of the large-scale structure inferred from late-time probes pose a serious challenge to the concordance $\Lambda$CDM cosmological model. State-of-the-art data from the Planck satellite predicts a higher rate of structure growth than what preferred by low-redshift observables. Such tension has hitherto eluded conclusive explanations in terms of straightforward modifications to $\Lambda$CDM, e.g. the inclusion of massive neutrinos or a dynamical dark energy component. Here, we investigate models of 'quartessence' -- a single dark component mimicking both dark matter and dark energy -- whose non-vanishing sound speed inhibits structure growth at late times on scales smaller than its corresponding Jeans' length. In principle, this could reconcile high- and low-redshift observations. We put this hypothesis to test against temperature and polarisation spectra from the latest Planck release, SDSS DR12 measurements of baryon acoustic oscillations and redshift-space distortions, and cosmic shear correlation functions from KiDS. This the first time that any specific model of quartessence is applied to actual data. We show that, if we naively apply $\Lambda$CDM nonlinear prescription to quartessence, the combined data sets allow for tight constraints on the model parameters. Apparently, quartessence alleviates the tension between the total matter fraction and late-time structure clustering, although in fact the tension is transferred from the latter to the quartessence sound speed parameter. However, we found that this strongly depends upon information from nonlinear scales. Indeed, if we relax this assumption, quartessence models appear still viable. For this reason, we argue that the nonlinear behaviour of quartessence deserves further investigation and may lead to a deeper understanding of the physics of the dark 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.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

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 measure the cross-correlation between Fermi-LAT gamma-ray photons and over 1000 deg$^2$ of weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS), and the Kilo Degree Survey (KiDS). We present the first measurement of tomographic weak lensing cross-correlations and the first application of spectral binning to cross-correlations between gamma rays and weak lensing. The measurements are performed using an angular power spectrum estimator while the covariance is estimated using an analytical prescription. We verify the accuracy of our covariance estimate by comparing it to two internal covariance estimators. Based on the non-detection of a cross-correlation signal, we derive constraints on weakly interacting massive particle (WIMP) dark matter. We compute exclusion limits on the dark matter annihilation cross-section $\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle mass $m_\rm{DM}$. We find that in the absence of a cross-correlation signal, tomography does not significantly improve the constraining power of the analysis. Assuming a strong contribution to the gamma-ray flux due to small-scale clustering of dark matter and accounting for known astrophysical sources of gamma rays, we exclude the thermal relic cross-section for particle masses of $m_\rm{DM}\lesssim 20$ GeV.

We report the detection of a cross-correlation signal between {\it Fermi} Large Area Telescope diffuse gamma-ray maps and catalogs of clusters. In our analysis, we considered three different catalogs: WHL12, redMaPPer and PlanckSZ. They all show a positive correlation with different amplitudes, related to the average mass of the objects in each catalog, which also sets the catalog bias. The signal detection is confirmed by the results of a stacking analysis. The cross-correlation signal extends to rather large angular scales, around 1 degree, that correspond, at the typical redshift of the clusters in these catalogs, to a few to tens of Mpc, i.e. the typical scale-length of the large scale structures in the Universe. Most likely this signal is contributed by the cumulative emission from AGNs associated to the filamentary structures that converge toward the high peaks of the matter density field in which galaxy clusters reside. In addition, our analysis reveals the presence of a second component, more compact in size and compatible with a point-like emission from within individual clusters. At present, we cannot distinguish between the two most likely interpretations for such a signal, i.e. whether it is produced by AGNs inside clusters or if it is a diffuse gamma-ray emission from the intra-cluster medium. We argue that this latter, intriguing, hypothesis might be tested by applying this technique to a low redshift large mass cluster sample.

Weak lensing has become an increasingly important tool in cosmology and the use of galaxy shapes to measure cosmic shear has become routine. The weak-lensing distortion tensor contains two other effects in addition to the two components of shear: the convergence and rotation. The rotation mode is not measurable using the standard cosmic shear estimators based on galaxy shapes, as there is no information on the original shapes of the images before they were lensed. Due to this, no estimator has been proposed for the rotation mode in cosmological weak-lensing surveys, and the rotation mode has never been constrained. Here, we derive an estimator for this quantity, which is based on the use of radio polarisation measurements of the intrinsic position angles of galaxies. The rotation mode can be sourced by physics beyond $\Lambda$CDM, and also offers the chance to perform consistency checks of $\Lambda$CDM and of weak-lensing surveys themselves. We present simulations of this estimator and show that, for the pedagogical example of cosmic string spectra, this estimator could detect a signal that is consistent with the constraints from Planck. We examine the connection between the rotation mode and the shear $B$-modes and thus how this estimator could help control systematics in future radio weak-lensing surveys.

The Square Kilometre Array (SKA) will produce spectroscopic surveys of tens to hundreds of millions of HI galaxies, eventually covering 30,000 sq. deg. and reaching out to redshift z~2. The huge volumes probed by the SKA will allow for some of the best constraints on primordial non-Gaussianity, based on measurements of the large-scale power spectrum. We investigate various observational set-ups for HI galaxy redshift surveys, compatible with the SKA Phase 1 and Phase 2 (full SKA) configurations. We use the corresponding number counts and bias for each survey from realistic simulations and derive the magnification bias and the evolution of source counts directly from these. For the first time, we produce forecasts that fully include the general relativistic effects on the galaxy number counts. These corrections to the standard analysis become important on very large scales, where the signal of primordial non-Gaussianity grows strongest. Our results show that, for the full survey, the non-Gaussianity parameter fNL can be constrained down to an accuracy of 1.54. This improves the current limit set by the Planck satellite by a factor of five, using a completely different approach.

We construct forecasts for cosmological parameter constraints from weak gravitational lensing surveys involving the Square Kilometre Array (SKA). Considering matter content, dark energy and modified gravity parameters, we show that the first phase of the SKA (SKA1) can be competitive with other Stage III experiments such as the Dark Energy Survey (DES) and that the full SKA (SKA2) can potentially form tighter constraints than Stage IV optical weak lensing experiments, such as those that will be conducted with LSST, WFIRST-AFTA or Euclid-like facilities. Using weak lensing alone, going from SKA1 to SKA2 represents improvements by factors of $\sim10$ in matter, $\sim10$ in dark energy and $\sim5$ in modified gravity parameters. We also show, for the first time, the powerful result that comparably tight constraints (within $\sim5\%$) for both Stage III and Stage IV experiments, can be gained from cross-correlating shear maps between the optical and radio wavebands, a process which can also eliminate a number of potential sources of systematic errors which can otherwise limit the utility of weak lensing cosmology.

In this third paper of a series on radio weak lensing for cosmology with the Square Kilometre Array, we scrutinise synergies between cosmic shear measurements in the radio and optical/near-IR bands for mitigating systematic effects. We focus on three main classes of systematics: (i) experimental systematic errors in the observed shear; (ii) signal contamination by intrinsic alignments; and (iii) systematic effects due to an incorrect modelling of non-linear scales. First, we show that a comprehensive, multi-wavelength analysis provides a self-calibration method for experimental systematic effects, only implying <50% increment on the errors on cosmological parameters. We also illustrate how the cross-correlation between radio and optical/near-IR surveys alone is able to remove residual systematics with variance as large as 0.00001, i.e. the same order of magnitude of the cosmological signal. This also opens the possibility of using such a cross-correlation as a means to detect unknown experimental systematics. Secondly, we demonstrate that, thanks to polarisation information, radio weak lensing surveys will be able to mitigate contamination by intrinsic alignments, in a way similar but fully complementary to available self-calibration methods based on position-shear correlations. Lastly, we illustrate how radio weak lensing experiments, reaching higher redshifts than those accessible to optical surveys, will probe dark energy and the growth of cosmic structures in regimes less contaminated by non-linearities in the matter perturbations. For instance, the higher-redshift bins of radio catalogues peak at z~0.8-1, whereas their optical/near-IR counterparts are limited to z<0.5-0.7. This translates into having a cosmological signal 2 to 5 times less contaminated by non-linear perturbations.

We construct a pipeline for simulating weak lensing cosmology surveys with the Square Kilometre Array (SKA), taking as inputs telescope sensitivity curves; correlated source flux, size and redshift distributions; a simple ionospheric model; source redshift and ellipticity measurement errors. We then use this simulation pipeline to optimise a 2-year weak lensing survey performed with the first deployment of the SKA (SKA1). Our assessments are based on the total signal-to-noise of the recovered shear power spectra, a metric that we find to correlate very well with a standard dark energy figure of merit. We first consider the choice of frequency band, trading off increases in number counts at lower frequencies against poorer resolution; our analysis strongly prefers the higher frequency Band 2 (950-1760 MHz) channel of the SKA-MID telescope to the lower frequency Band 1 (350-1050 MHz). Best results would be obtained by allowing the centre of Band 2 to shift towards lower frequency, around 1.1 GHz. We then move on to consider survey size, finding that an area of 5,000 square degrees is optimal for most SKA1 instrumental configurations. Finally, we forecast the performance of a weak lensing survey with the second deployment of the SKA. The increased survey size (3$\pi$\,steradian) and sensitivity improves both the signal-to-noise and the dark energy metrics by two orders of magnitude.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

We recently proposed to cross-correlate the diffuse extragalactic gamma-ray background with the gravitational lensing signal of cosmic shear. This represents a novel and promising strategy to search for annihilating or decaying particle dark matter (DM) candidates. In the present work, we demonstrate the potential of a tomographic-spectral approach: measuring the cross-correlation in separate bins of redshift and energy significantly improves the sensitivity to a DM signal. Indeed, the technique proposed here takes advantage of the different scaling of the astrophysical and DM components with redshift and, simultaneously, of their different energy spectra and different angular extensions. The sensitivity to a particle DM signal is extremely promising even when the DM-induced emission is quite faint. We first quantify the prospects of detecting DM by cross-correlating the Fermi Large Area Telescope (LAT) diffuse gamma-ray background with the cosmic shear expected from the Dark Energy Survey. Under the hypothesis of a significant subhalo boost, such a measurement can deliver a 5-sigma detection of DM, if the DM particle is lighter than 300 GeV and has a thermal annihilation rate. We then forecast the capability of the European Space Agency Euclid satellite (whose launch is planned for 2020), in combination with an hypothetical future gamma-ray detector with slightly improved specifications compared to current telescopes. We predict that the cross-correlation of their data will allow a measurement of the DM mass with an uncertainty of a factor of 1.5-2, even for moderate subhalo boosts, for DM masses up to few hundreds of GeV and thermal annihilation rates.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Next-generation cosmological surveys will probe ever larger volumes of the Universe, including the largest scales, near and beyond the horizon. On these scales, the galaxy power spectrum carries signatures of local primordial non-Gaussianity (PNG) and horizon-scale general relativistic (GR) effects. However, cosmic variance limits the detection of horizon-scale effects. Combining different surveys via the multi-tracer method allows us to reduce the effect down cosmic variance. This method benefits from large bias differences between two tracers of the underlying dark matter distribution, which suggests a multi-wavelength combination of large volume surveys that are planned on a similar timescale. We show that the combination of two contemporaneous surveys, a large neutral hydrogen intensity mapping survey in SKA Phase1 and a Euclid-like photometric survey, will provide unprecedented constraints on PNG as well as detection of the GR effects. We forecast that the error on local PNG will break through the cosmic variance limit on cosmic microwave background surveys and achieve $\sigma(f_{NL})\simeq1.4-0.5$, depending on assumed priors, bias, and sky coverage. GR effects are more robust to changes in the assumed fiducial model, and we forecast that they can be detected with a signal-to-noise of about $14$.

Non-Gaussianity in the primordial fluctuations that seeded structure formation produces a signal in the galaxy power spectrum on very large scales. This signal contains vital information about the primordial Universe, but it is very challenging to extract, because of cosmic variance and large-scale systematics - especially after the Planck experiment has already ruled out a large amplitude for the signal. Cosmic variance and experimental systematics can be alleviated by the multi-tracer method. Here we address another systematic - introduced by not using the correct relativistic analysis of the power spectrum on very large scales. In order to reduce the errors on fNL, we need to include measurements on the largest possible scales. Failure to include the relativistic effects on these scales can introduce significant bias in the best-fit value of fNL from future galaxy surveys.

Future galaxy redshift surveys aim at probing the clustering of the cosmic large-scale structure with unprecedented accuracy, thus complementing cosmic microwave background experiments in the quest to deliver the most precise and accurate picture ever of our Universe. Analyses of such measurements are usually performed within the context of the so-called vanilla LCDM model - the six-parameter phenomenological model which, for instance, emerges from best fits against the recent data obtained by the Planck satellite. Here, we show that such an approach is prone to subtle systematics when the Gaussianity of primordial fluctuations is concerned. In particular, we demonstrate that, if we neglect even a tiny amount of primordial non-Gaussianity - fully consistent with current limits - we shall introduce spurious biases in the reconstruction of cosmological parameters. This is a serious issue that must be properly accounted for in view of accurate (as well as precise) cosmology.

We report the measurement of the angular power spectrum of cross-correlation between the unresolved component of the Fermi-LAT gamma-ray sky-maps and the CMB lensing potential map reconstructed by the Planck satellite. The matter distribution in the Universe determines the bending of light coming from the last scattering surface. At the same time, the matter density drives the growth history of astrophysical objects, including their capability at generating non-thermal phenomena, which in turn give rise to gamma-ray emissions. The Planck lensing map provides information on the integrated distribution of matter, while the integrated history of gamma-ray emitters is imprinted in the Fermi-LAT sky maps. We report here the first evidence of their correlation. We find that the multipole dependence of the cross-correlation measurement is in agreement with current models of the gamma-ray luminosity function for AGN and star forming galaxies, with a statistical evidence of 3.0$\sigma$. Moreover, its amplitude can in general be matched only assuming that these extra-galactic emitters are also the bulk contribution of the measured isotopic gamma-ray background (IGRB) intensity. This leaves little room for a big contribution from galactic sources to the IGRB measured by Fermi-LAT, pointing toward a direct evidence of the extragalactic origin of the IGRB.

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.

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.

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.

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.

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 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.

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.

The statistics of peaks in weak lensing convergence maps is a promising tool to investigate both the properties of dark matter haloes and constrain the cosmological parameters. We study how the number of detectable peaks and its scaling with redshift depend upon the cluster dark matter halo profiles and use peak statistics to constrain the parameters of the mass - concentration (MC) relation. We investigate which constraints the Euclid mission can set on the MC coefficients also taking into account degeneracies with the cosmological parameters. To this end, we first estimate the number of peaks and its redshift distribution for different MC relations. We find that the steeper the mass dependence and the larger the normalisation, the higher is the number of detectable clusters, with the total number of peaks changing up to $40\%$ depending on the MC relation. We then perform a Fisher matrix forecast of the errors on the MC relation parameters as well as cosmological parameters. We find that peak number counts detected by Euclid can determine the normalization $A_v$, the mass $B_v$ and redshift $C_v$ slopes and intrinsic scatter $\sigma_v$ of the MC relation to an unprecedented accuracy being $\sigma(A_v)/A_v = 1\%$, $\sigma(B_v)/B_v = 4\%$, $\sigma(C_v)/C_v = 9\%$, $\sigma(\sigma_v)/\sigma_v = 1\%$ if all cosmological parameters are assumed to be known. Should we relax this severe assumption, constraints are degraded, but remarkably good results can be restored setting only some of the parameters or combining peak counts with Planck data. This precision can give insight on competing scenarios of structure formation and evolution and on the role of baryons in cluster assembling. Alternatively, for a fixed MC relation, future peaks counts can perform as well as current BAO and SNeIa when combined with Planck.

Alterations of the gravity Lagrangian introduced in modified torsion gravity theories---also referred to as $f(T)$ gravity---allows for an accelerated expansion in a matter dominated Universe. In this framework, the cosmic speed up is driven by an effective `torsion fluid'. Besides the background evolution of the Universe, structure formation is also modified because of a time dependent effective gravitational constant. Here, we investigate the imprints of $f(T)$ gravity on galaxy clustering and weak gravitational lensing to the aim of understanding whether future galaxy surveys could constrain torsion gravity and discriminate amongst it and standard general relativity. Specifically, we compute Fisher matrix forecasts for two viable $f(T)$ models to both infer the accuracy on the measurement of the model parameters and evaluate the power that a combined clustering and shear analysis will have as a tool for model selection. We find that with such a combination of probes it will indeed be possible to tightly constrain $f(T)$ model parameters. Moreover, the Occam's razor provided by the Bayes factor will allow us to confirm an $f(T)$ power-law extension of the concordance $\Lambda{\rm CDM}$ model, were a value larger than 0.02 of its power-law slope measured, whereas in $\Lambda{\rm CDM}$ it is exactly 0.

In this paper we investigate magnetic fields generated in the early Universe. These fields are important candidates at explaining the origin of astrophysical magnetism observed in galaxies and galaxy clusters, whose genesis is still by and large unclear. Compared to the standard inflationary power spectrum, intermediate to small scales would experience further substantial matter clustering, were a cosmological magnetic field present prior to recombination. As a consequence, the bias and redshift distribution of galaxies would also be modified. Hitherto, primordial magnetic fields (PMFs) have been tested and constrained with a number of cosmological observables, e.g. the cosmic microwave background radiation, galaxy clustering and, more recently, weak gravitational lensing. Here, we explore the constraining potential of the density fluctuation bias induced by gravitational lensing magnification onto the galaxy-galaxy angular power spectrum. Such an effect is known as magnification bias. Compared to the usual galaxy clustering approach, magnification bias helps in lifting the pathological degeneracy present amongst power spectrum normalisation and galaxy bias. This is because magnification bias cross-correlates galaxy number density fluctuations of nearby objects with weak lensing distortions of high-redshift sources. Thus, it takes advantage of the gravitational deflection of light, which is insensitive to galaxy bias but powerful in constraining the density fluctuation amplitude. To scrutinise the potentiality of this method, we adopt a deep and wide-field spectroscopic galaxy survey. We show that magnification bias does contain important information on primordial magnetism, which will be useful in combination with galaxy clustering and shear. We find we shall be able to rule out at 95.4% CL amplitudes of PMFs larger than 0.0005 nG for values of the PMF power spectral index ~0.

We explore the use of different radio galaxy populations as tracers of different mass halos and therefore, with different bias properties, to constrain primordial non-Gaussianity of the local type. We perform a Fisher matrix analysis based on the predicted auto and cross angular power spectra of these populations, using simulated redshift distributions as a function of detection flux and the evolution of the bias for the different galaxy types (Star forming galaxies, Starburst galaxies, Radio-Quiet Quasars, FRI and FRII AGN galaxies). We show that such a multi-tracer analysis greatly improves the information on non-Gaussianity by drastically reducing the cosmic variance contribution to the overall error budget. By using this method applied to future surveys, we predict a constraint of sigma_fnl=3.6 on the local non-Gaussian parameter for a galaxy detection flux limit of 10 \muJy and sigma_fnl=2.2 for 1 \muJy. We show that this significantly improves on the constraints obtained when using the whole undifferentiated populations (sigma_fnl=48 for 10 \muJy and sigma_fnl=12 for 1 \muJy). We conclude that continuum radio surveys alone have the potential to constrain primordial non-Gaussianity to an accuracy at least a factor of two better than the present constraints obtained with Planck data on the CMB bispectrum, opening a window to obtain sigma_fnl~1 with the Square Kilometer Array.

The large-scale structure of the Universe supplies crucial information about the physical processes at play at early times. Unresolved maps of the intensity of 21 cm emission from neutral hydrogen HI at redshifts z~1-5 are the best hope of accessing the ultralarge-scale information, directly related to the early Universe. A purpose-built HI intensity experiment may be used to detect the large scale effects of primordial non-Gaussianity, placing stringent bounds on different models of inflation. We argue that it may be possible to place tight constraints on the non-Gaussianity parameter f_NL, with an error close to ~1.

Both cosmic shear and cosmological gamma-ray emission stem from the presence of Dark Matter (DM) in the Universe: DM structures are responsible for the bending of light in the weak lensing regime and those same objects can emit gamma-rays, either because they host astrophysical sources (active galactic nuclei or star-forming galaxies) or directly by DM annihilations (or decays, depending on the properties of the DM particle). Such gamma-rays should therefore exhibit strong correlation with the cosmic shear signal. In this Letter, we compute the cross-correlation angular power spectrum of cosmic shear and gamma-rays produced by the annihilation/decay of Weakly Interacting Massive Particle (WIMP) DM, as well as from astrophysical sources. We show that this observable provides novel information on the composition of the Extra-galactic Gamma-ray Background (EGB), since the amplitude and shape of the cross-correlation signal strongly depends on which class of source is responsible for the gamma-ray emission. If the DM contribution to the EGB is significant (at least in a definite energy range), although compatible with current observational bounds, its strong correlation with the cosmic shear makes such signal potentially detectable by combining Fermi-LAT data with forthcoming galaxy surveys, like Dark Energy Survey and Euclid. At the same time, the same signal would demonstrate that the weak lensing observables are indeed due to particle DM matter and not to possible modifications of General Relativity.

We show how future gravitational-wave detectors would be able of discriminating between the concordance LCDM cosmological model and up-to-date competing alternatives, e.g. dynamical dark energy models (DE) or modified gravity theories (MG). Our method consists in using the weak-lensing magnification effect that affects a standard-siren signal because of its travelling trough the Universe's large-scale structure. As a demonstration, we present constraints on DE and MG from proposed gravitational-wave detectors, namely ET and DECIGO/BBO.

Weak gravitational lensing by galaxy clusters on faint higher redshift galaxies has been traditionally used to study the cluster mass distribution and as a tool to identify clusters as peaks in the shear maps. However, it becomes soon clear that peaks statistics can also be used as a way to constrain the underlying cosmological model due to its dependence on both the cosmic expansion rate and the growth rate of structures. This feature makes peak statistics particularly interesting from the point of view of discriminating between General Relativity and modified gravity. Here we consider a general class of $f(R)$ theories and compute the observable mass function based on the aperture mass statistics. We complement our theoretical analysis with a Fisher matrix forecast of the constraints that an Euclid\,-\,like survey can impose on the $f(R)$ model parameters. We show that peak statistics alone can in principle discriminate between General Relativity and $f(R)$ models and strongly constrain the $f(R)$ parameters that are sensitive to the non-linear growth of structure. However, further analysis is needed in order to include possible selection function in the peaks redshift determination.

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.

Generalised Teleparallel gravity, also referred to as f(T) gravity, has been recently proposed as an extended theory of gravitation able to give rise to an accelerated expansion in a matter only universe. The cosmic speed up is driven by an effective torsion fluid whose equation of state depend on the f(T) function entering the modified gravity Lagrangian. We focus on two particular choices for f(T) which share the nice property to emulate a phantom divide crossing as suggested by some recent data. We check their viability contrasting the predicted background dynamics to the Hubble diagram as traced by both Type Ia Supernovae (SNeIa) and Gamma Ray Bursts (GRBs), the measurement of the rate expansion H(z), the Baryon Acoustic Oscillations (BAOs) at different redshifts, and the Cosmic Microwave Background Radiation (CMBR) distance priors. Both f(T) models turn out to be in very good agreement with this large dataset so that we also investigate whether it is possible to discriminate among them relying on the different growth factors.

In the very last years, cosmological models where the properties of the dark components of the Universe - dark matter and dark energy - are accounted for by a single "dark fluid" have drawn increasing attention and interest. Amongst many proposals, Unified Dark Matter (UDM) cosmologies are promising candidates as effective theories. In these models, a scalar field with a non-canonical kinetic term in its Lagrangian mimics both the accelerated expansion of the Universe at late times and the clustering properties of the large-scale structure of the cosmos. However, UDM models also present peculiar behaviours, the most interesting one being the fact that the perturbations in the dark-matter component of the scalar field do have a non-negligible speed of sound. This gives rise to an effective Jeans scale for the Newtonian potential, below which the dark fluid does not cluster any more. This implies a growth of structures fairly different from that of the concordance LCDM model. In this paper, we demonstrate that forthcoming large-scale surveys will be able to discriminate between viable UDM models and LCDM to a good degree of accuracy. To this purpose, the planned Euclid satellite will be a powerful tool, since it will provide very accurate data on galaxy clustering and the weak lensing effect of cosmic shear. Finally, we also exploit the constraining power of the ongoing CMB Planck experiment. Although our approach is the most conservative, with the inclusion of only well-understood, linear dynamics, in the end we also show what could be done if some amount of non-linear information were included.

The observed accelerated cosmic expansion can be a signature of fourth\,-\,order gravity theories, where the acceleration of the Universe is a consequence of departures from Einstein General Relativity, rather than the sign of the existence of a fluid with negative pressure. In the fourth\,-\,order gravity theories, the gravity Lagrangian is described by an analytic function $f(R)$ of the scalar curvature $R$ subject to the demanding conditions that no detectable deviations from standard GR is observed on the Solar System scale. Here we consider two classes of $f(R)$ theories able to pass Solar System tests and investigate their viability on cosmological scales. To this end, we fit the theories to a large dataset including the combined Hubble diagram of Type Ia Supernovae and Gamma Ray Bursts, the Hubble parameter $H(z)$ data from passively evolving red galaxies, Baryon Acoustic Oscillations extracted from the seventh data release of the Sloan Digital Sky Survey (SDSS) and the distance priors from the Wilkinson Microwave Anisotropy Probe seven years (WMAP7) data. We find that both classes of $f(R)$ fit very well this large dataset with the present\,-\,day values of the matter density, Hubble constant and deceleration parameter in agreement with previous estimates; however, the strong degeneracy among the $f(R)$ parameters prevents us from strongly constraining their values. We also derive the growth factor $g = d\ln{\delta}/d\ln{a}$, with $\delta = \delta \rho_M/\rho_M$ the matter density perturbation, and show that it can still be well approximated by $g(z) \propto \Omega_M(z)^{\gamma}$. We finally constrain $\gamma$ (on some representative scales) and investigate its redshift dependence to see whether future data can discriminate between these classes of $f(R)$ theories and standard dark energy models.

We present the cosmic shear signal predicted by two viable cosmological models in the framework of modified-action f(R) theories. We use f(R) models where the current accelerated expansion of the Universe is a direct consequence of the modified gravitational Lagrangian rather than Dark Energy (DE), either in the form of vacuum energy/cosmological constant or of a dynamical scalar field (e.g. quintessence). We choose Starobinsky's (St) and Hu & Sawicki's (HS) f(R) models, which are carefully designed to pass the Solar System gravity tests. In order to further support - or rule out - f(R) theories as alternative candidates to the DE hypothesis, we exploit the power of weak gravitational lensing, specifically of cosmic shear. We calculate the tomographic shear matrix as it would be measured by the upcoming ESA Cosmic Vision Euclid satellite. We find that in the St model the cosmic shear signal is almost completely degenerate with LCDM, but it is easily distinguishable in the HS model. Moreover, we compute the corresponding Fisher matrix for both the St and HS models, thus obtaining forecasts for their cosmological parameters. Finally, we show that the Bayes factor for cosmic shear will definitely favour the HS model over LCDM if Euclid measures a value larger than ~0.02 for the extra HS parameter n_HS.

Flexion is the weak lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. The flexion signal-to-noise ratio can be an order of magnitude larger than that of shear. For the first time, we show how this makes flexion an invaluable tool for discriminating among alternative cosmological models. We analyse a scalar field model of unified dark matter and dark energy, a brane-world cosmology and two f(R) modified-action theories. We show that these models can be distinguished from LCDM at several standard deviations by measuring the power spectrum of cosmic flexion.

We present parameter estimation forecasts for future 3D cosmic shear surveys for a class of Unified Dark Matter (UDM) models, where a single scalar field mimics both Dark Matter (DM) and Dark Energy (DE). These models have the advantage that they can describe the dynamics of the Universe with a single matter component providing an explanation for structure formation and cosmic acceleration. A crucial feature of the class of UDM models we use in this work is characterised by a parameter, c_inf (in units of the speed of light c=1), that is the value of the sound speed at late times, and on which structure formation depends. We demonstrate that the properties of the DM-like behaviour of the scalar field can be estimated with very high precision with large-scale, fully 3D weak lensing surveys. We found that 3D weak lensing significantly constrains c_inf, and we find minimal errors 0.00003, for the fiducial value c_inf=0.001, and 0.000026, for c_inf=0.012. Moreover, we compute the Bayesian evidence for UDM models over the LCDM model as a function of c_inf. For this purpose, we can consider the LCDM model as a UDM model with c_inf=0. We find that the expected evidence clearly shows that the survey data would unquestionably favour UDM models over the LCDM model, for the values c_inf>0.001.

A phenomenological extension of the well-known brane-world cosmology of Dvali, Gabadadze and Porrati (eDGP) has recently been proposed. In this model, a cosmological-constant-like term is explicitly present as a non-vanishing tension sigma on the brane, and an extra parameter alpha tunes the cross-over scale r_c, the scale at which higher dimensional gravity effects become non negligible. Since the Hubble parameter in this cosmology reproduces the same LCDM expansion history, we study how upcoming weak lensing surveys, such as Euclid and DES (Dark Energy Survey), can confirm or rule out this class of models. We perform Markov Chain Monte Carlo simulations to determine the parameters of the model, using Type Ia Supernov\ae, H(z) data, Gamma Ray Bursts and Baryon Acoustic Oscillations. We also fit the power spectrum of the temperature anisotropies of the Cosmic Microwave Background to obtain the correct normalisation for the density perturbation power spectrum. Then, we compute the matter and the cosmic shear power spectra, both in the linear and non-linear regimes. The latter is calculated with the two different approaches of Hu and Sawicki (2007) (HS) and Khoury and Wyman (2009) (KW). With the eDGP parameters coming from the Markov Chains, KW reproduces the LCDM matter power spectrum at both linear and non-linear scales and the LCDM and eDGP shear signals are degenerate. This result does not hold with the HS prescription: Euclid can distinguish the eDGP model from LCDM because their expected power spectra roughly differ by the 3sigma uncertainty in the angular scale range 700<l<3000; on the contrary, the two models differ at most by the 1sigma uncertainty over the range 500<l<3000 in the DES experiment and they are virtually indistinguishable.

Fourth - order gravity theories have received much interest in recent years thanks to their ability to provide an accelerated cosmic expansion in a matter only universe. In these theories, the Lagrangian density of the gravitational field has the form R + f(R), and the explicit choice of the arbitrary function f(R) must meet the local tests of gravity and the constraints from the primordial abundance of the light elements. Two popular classes of f(R) models, which are expected to fulfill all the above requirements, have recently been proposed. However, neither of these models has ever been quantitatively tested against the available astrophysical data. Here, by combining Type Ia Supernovae and Gamma Ray Bursts, we investigate the ability of these models to reproduce the observed Hubble diagram over the redshift range (0, 7). We find that both models fit very well this dataset with the present day values of the matter density and deceleration parameters which agree with previous estimates. However, the strong degeneracy among the f(R) parameters prevents us from putting strong constraints on the values of these parameters; nevertheless, we can identify the regions of the parameter space that should, in principle, be carefully explored with future data and dynamical probes in order to discriminate among f(R) theories and standard dark energy models.

Weak gravitational lensing is a powerful tool for studying both the geometry and the dynamics of the Universe. Its power spectrum contains information on the sources emitting photons and on the large--scale structures that these hotons cross. We calculate the weak lensing cosmic convergence and shear power spectra, in linear theory and Limber's approximation, for two different classes of cosmological models: the standard \LambdaCDM and Unified Dark Matter (UDM) models. The latter models attempt to solve the problems of the dark matter in the dynamics of galaxies and galaxy clusters and of the late-time acceleration of the Universe expansion by introducing a scalar field that mimics both dark matter and dark energy. A crucial feature of the UDM models is the speed of sound c_\infty, that is the value of the sound speed at late times, on which structure formation depends. In this paper, we provide the predictions of the UDM models for the weak lensing signal,with various values of c_\infty. We consider both the Cosmic Microwave Background and background galaxies at different redshifts as sources for the weak lensing power spectra. We find that the power spectra in UDM models are more sensitive to the variations of c_\infty for sources located at low redshifts. Moreover, we find that for c_\infty>10^{-3} (in units of the speed of light), the UDM weak lensing convergence power spectrum C^{\kappa\kappa}(l) for background galaxies is strongly suppressed with respect to the \LambdaCDM spectrum, particularly at small angular scales l\gtrsim100.