
We consider a selfconsistent 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 subcase, 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 investigate the possibility of using future photometric and radio surveys
to constrain the power spectrum of primordial fluctuations that is predicted by
inflationary models with a violation of the slowroll phase. We forecast
constraints with a Fisher analysis on the amplitude of the parametrized
features on ultralarge scales, in order to assess whether these could be
distinguishable over the cosmic variance. We find that the next generation of
photometric and radio surveys has the potential to test these models at a
sensitivity better than current CMB experiments and that the synergy between
galaxy and CMB observations is able to constrain models with many extra
parameters. In particular, an SKA continuum survey with a huge sky coverage and
a flux threshold of a few $\mu$Jy could confirm the presence of a new phase in
the early Universe at more than 3$\sigma$.

The angular power spectrum is a gaugeindependent observable that is in
principle the natural tool for analysing spectroscopic galaxy number counts. In
practice, the problem is that the computational requirements for
nextgeneration 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 crossbin correlations
beyond the thick bins are neglected. Constraints on cosmological parameters
from the hybrid method are comparable to those from the standard power spectrum
analysis  but they have the advantage that cosmic evolution, wideangle and
lensing effects are naturally included, while no AlcockPaczynski 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.

The galaxy bispectrum is affected on equality scales and above by
relativistic observational effects, at linear and nonlinear order. These
lightcone effects include local contributions from Doppler and gravitational
potential terms, as well as integrated contributions like lensing, together
with all the couplings at nonlinear order. We recently presented the correction
to the galaxy bispectrum from all local lightcone effects up to second order in
perturbations, using a planeparallel approximation. Here we update our
previous result by including the effects from relativistic nonlinear dynamical
evolution. We show that these dynamical effects make a significant contribution
to the projection effects.

We study synergies between HI 21cm and H$\alpha$ intensity map observations,
focusing on SKA1 and SPHERExlike surveys. We forecast how well such a
combination can measure features in the angular power spectrum on the largest
scales, that arise from primordial nonGaussianity and from general
relativistic effects. For the first time we consider Doppler, SachsWolfe and
integrated SW effects separately. We confirm that the singletracer surveys on
their own cannot detect general relativistic effects and can constrain the
nonGaussianity parameter $f_{\rm NL}$ only slightly better than Planck. Using
the multitracer technique, constraints on $f_{\rm NL}$ can be pushed down to
$\sim 1$. Amongst the general relativistic effects, the Doppler term is
detectable with the multitracer. The SachsWolfe terms and the integrated SW
effect are still not detectable.

In the concordance model of the Universe, the matter distribution  as
observed in galaxy number counts or the intensity of line emission (such as the
21cm line of neutral hydrogen)  should have a kinematic dipole due to the
Sun's motion relative to the CMB restframe. This dipole should be aligned with
the kinematic dipole in the CMB temperature. Accurate measurement of the
direction of the matter dipole will become possible with future galaxy surveys,
and this will be a critical test of the foundations of the concordance model.
The amplitude of the matter dipole is also a potential cosmological probe. We
derive formulas for the amplitude of the kinematic dipole in galaxy redshift
and intensity mapping surveys, taking into account the Doppler, aberration and
other relativistic effects. The amplitude of the matter dipole can be
significantly larger than that of the CMB dipole. Its redshift dependence
encodes information on the evolution of the Universe and on the tracers, and we
discuss possible ways to determine the amplitude.

General relativistic imprints on the galaxy bispectrum arise from
observational (or projection) effects. The lightcone projection effects include
local contributions from Doppler and gravitational potential terms, as well as
lensing and other integrated contributions. We recently presented for the first
time, the correction to the galaxy bispectrum from all local lightcone
projection effects up to second order in perturbations. Here we provide the
details underlying this correction, together with further results and
illustrations. For moderately squeezed shapes, the correction to the Newtonian
prediction is ~30% on equality scales at z ~ 1. We generalise our recent
results to include the contribution, up to second order, of magnification bias
(which affects some of the local terms) and evolution bias.

We discuss the groundbreaking 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 multiwavelength data will give unique
additional information, such as exquisite constraints on primordial
nonGaussianity using the multitracer 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 preSKA 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 Lband, which is quite unique over such large areas and
will allow studies of the largescale structure of the Universe out to high
redshifts, complementing the galaxy HI survey to form a transformational
multiwavelength 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 SKA1MID for
cosmological applications and other commensal surveys, as described in the top
priority SKA key science projects (abridged).

Objects falling into an overdensity appear larger on its near side and
smaller on its far side than other objects at the same redshift. This produces
a dipolar pattern of magnification, primarily as a consequence of the Doppler
effect. At low redshift this Doppler magnification completely dominates the
usual integrated gravitational lensing contribution to the lensing
magnification. We show that one can optimally observe this pattern by
extracting the dipole in the crosscorrelation of number counts and galaxy
sizes. This dipole allows us to almost completely remove the contribution from
gravitational lensing up to redshift 0.5, and even at high redshift z~1 the
dipole picks up the Doppler magnification predominantly. Doppler magnification
should be easily detectable in current and upcoming optical and radio surveys;
by forecasting for telescopes such as the SKA, we show that this technique is
competitive with using peculiar velocities via redshiftspace distortions to
constrain dark energy. It produces similar yet complementary constraints on the
cosmological model to those found using measurements of the cosmic shear.

The Hubble diagram is one of the cornerstones of observational cosmology. It
is usually analysed assuming that, on average, the underlying relation between
magnitude and redshift matches the prediction of a
FriedmannLema\^itreRobertsonWalker model. However, the inhomogeneity of the
Universe generically biases these observables, mainly due to peculiar
velocities and gravitational lensing, in a way that depends on the notion of
average used in theoretical calculations. In this article, we carefully derive
the notion of average which corresponds to the observation of the Hubble
diagram. We then calculate its bias at secondorder in cosmological
perturbations, and estimate the consequences on the inference of cosmological
parameters, for various current and future surveys. We find that this bias
deeply affects direct estimations of the evolution of the darkenergy equation
of state. However, errors in the standard inference of cosmological parameters
remain smaller than observational uncertainties, even though they reach percent
level on some parameters; they reduce to subpercent level if an optimal
distance indicator is used.

Nextgeneration galaxy surveys will increasingly rely on the galaxy
bispectrum to improve cosmological constraints, especially on primordial
nonGaussianity. A key theoretical requirement that remains to be developed is
the analysis of general relativistic effects on the bispectrum, which arise
from observing galaxies on the past lightcone, {as well as from relativistic
corrections to the dynamics}. {As an initial step towards a fully relativistic
analysis of the galaxy bispectrum, we compute for the first time the local
relativistic lightcone effects on the bispectrum,} which come from Doppler and
gravitational potential contributions. For the galaxy bispectrum, the problem
is much more complex than for the power spectrum, since we need the lightcone
corrections at second order. Modecoupling contributions at second order mean
that relativistic corrections can be nonnegligible at smaller scales than in
the case of the power spectrum. In a primordial Gaussian universe, we show that
the local lightcone corrections for squeezed shapes at $z\sim1$ mean that the
bispectrum can differ from the Newtonian prediction by $\gtrsim 10\%$ when the
short modes are $k\lesssim (50\,{\rm Mpc})^{1}$. These relativistic projection
effects, if ignored in the analysis of observations, could be mistaken for
primordial nonGaussianity. For upcoming surveys which probe equality scales
and beyond, {all relativistic lightcone effects and relativistic dynamical
corrections should be included} for an accurate measurement of primordial
nonGaussianity.

It is usually assumed that we will need to wait until nextgeneration surveys
like Euclid, LSST and SKA, in order to improve on the current best constraints
on primordial nonGaussianity from the Planck experiment. We show that two
contemporary surveys, with the SKA precursor MeerKAT and the Dark Energy Survey
(DES), can be combined using the multitracer technique to deliver an accuracy
on measurement of $f_{\rm NL}$ that is up to three times better than Planck.

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 nonGaussianity, based on
measurements of the largescale power spectrum. We investigate various
observational setups 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 nonGaussianity grows strongest. Our results show that, for the full
survey, the nonGaussianity 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.

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

Intensity mapping of the neutral hydrogen brightness temperature promises to
provide a threedimensional view of the universe on very large scales.
Nonlinear effects are typically thought to alter only the smallscale power,
but we show how they may bias the extraction of cosmological information
contained in the power spectrum on ultralarge scales. For linear perturbations
to remain valid on large scales, we need to renormalize perturbations at higher
order. In the case of intensity mapping, the secondorder contribution to
clustering from weak lensing dominates the nonlinear contribution at high
redshift. Renormalization modifies the mean brightness temperature and
therefore the evolution bias. It also introduces a term that mimics white
noise. These effects may influence forecasting analysis on ultralarge scales.

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 highprecision observations and new theoretical
ideas are fulfilled.

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

Nextgeneration 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
nonGaussianity (PNG) and horizonscale general relativistic (GR) effects.
However, cosmic variance limits the detection of horizonscale effects.
Combining different surveys via the multitracer 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 multiwavelength 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
Euclidlike 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.40.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
signaltonoise of about $14$.

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

We show that at second order, ensemble averages of observables and
directional averages do not commute due to gravitational lensing  observing
the same thing in many directions over the sky is not the same as taking an
ensemble average. In principle this noncommutativity is significant for a
variety of quantities that we often use as observables and can lead to a bias
in parameter estimation. We derive the relation between the ensemble average
and the directional average of an observable, at second order in perturbation
theory. We discuss the relevance of these two types of averages for making
predictions of cosmological observables, focusing on observables related to
distances and magnitudes. In particular, we show that the ensemble average of
the distance in a given observed direction is increased by gravitational
lensing, whereas the directional average of the distance is decreased. For a
generic observable, there exists a particular function of the observable that
is not affected by secondorder lensing perturbations. We also show that
standard areas have an advantage over standard rulers, and we discuss the
subtleties involved in averaging in the case of supernova observations.

It has recently been shown that secondorder corrections to the background
distanceredshift relation can build up significantly at large redshifts, due
to an aggregation of gravitational lensing events. This shifts the expectation
value of the distance to the CMB by 1%. In this paper we show that this shift
is already properly accounted for in standard CMB analyses. We clarify the role
that the area distance to the CMB plays in the presence of secondorder lensing
corrections.

NonGaussianity 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 largescale
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 multitracer 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
bestfit value of fNL from future galaxy surveys.

The Planck CMB experiment has delivered the best constraints so far on
primordial nonGaussianity, ruling out earlyUniverse models of inflation that
generate large nonGaussianity. Although small improvements in the CMB
constraints are expected, the next frontier of precision will come from future
largescale 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 SachsWolfe effect  the crosscorrelation of the radio
number counts with the CMB temperature anisotropies  helps to reduce
systematics on the large scales that are sensitive to nonGaussianity. Second,
optical data allows for crossidentification 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 nonGaussianity, 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 highredshift sources, so having
crossID for highz galaxies and an even higherz 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.

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

HI intensity mapping (IM) is a novel technique capable of mapping the
largescale 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 nonGaussianity can be used to test inflation; on small scales, by
measuring the sum of neutrino masses; and at high redshifts where nonstandard
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 largescale signal.