The growth rate of large-scale structure provides a powerful consistency test of the standard cosmological model and a probe of possible deviations from general relativity. We use a Fisher analysis to forecast constraints on the growth rate from a combination of next-generation spectroscopic surveys. In the overlap survey volumes, we use a multi-tracer analysis to significantly reduce the effect of cosmic variance. The non-overlap individual survey volumes are included in the Fisher analysis in order to utilise the entire volume. We use the observed angular power spectrum, which naturally includes all wide-angle and lensing effects and circumvents the need for an Alcock-Paczynski correction. Cross correlations between redshift bins are included by using a novel technique to avoid computation of the sub-dominant contributions. Marginalising over the standard cosmological parameters, as well as the clustering bias in each redshift bin, we find that the precision on $\gamma$ improves on the best single-tracer precision by up to $\sim$50\%.

We perform null tests of the concordance model, using $H(z)$ measurements that mimic next-generation surveys such as Euclid and the SKA. To this end, we deploy a non-parametric method, so that we make minimal assumptions about the fiducial cosmology as well as the statistical analysis. We produce simulations assuming different cosmological models in order to verify how well we can distinguish between their signatures. We find that SKA- and Euclid-like surveys should be able to discriminate sharply between the concordance and alternative dark energy models that are compatible with the Planck CMB data. We conclude that SKA and Euclid will be able to falsify the concordance model in a statistically significant way, if one of the benchmarks models represents the true Universe, without making assumptions about the underlying cosmology.

The key probes of the growth of large-scale structure are its rate $f$ and amplitude $\sigma_8$. Redshift space distortions in the galaxy power spectrum allow us to measure only the combination $f\sigma_8$, which can be used to constrain the standard cosmological model or alternatives. By using measurements of the galaxy-galaxy lensing cross-correlation spectrum or of the galaxy bispectrum, it is possible to break the $f\sigma_8$ degeneracy and obtain separate estimates of $f$ and $\sigma_8$ from the same galaxy sample. Currently there are only a handful of such separate measurements, but even this allows for improved constraints on cosmological models. We consider what can be achieved by a future nominal sample that delivers a $\sim 1\%$ constraint on $f$ and $\sigma_8$ separately, compared to the case with a similar precision on the combination $f\sigma_8$. For the six cosmological parameters of $\Lambda$CDM, we find improvements of $\sim\! 5$--$50\%$ on their constraints. For modified gravity models in the Horndeski class, the improvements on these standard parameters are $\sim\! 0$--$15\%$. However, the precision on the sum of neutrino masses improves by 65\% and there is a significant increase in the precision on the background and perturbation Horndeski parameters.

The rate at which the universe is expanding today is a fundamental parameter in cosmology which governs our understanding of structure formation and dark energy. However, current measurements of the Hubble constant, $H_0$, show a significant tension ($\sim 4-6\sigma$) between early- and late-Universe observations. There are ongoing efforts to check the diverse observational results and also to investigate possible theoretical ways to resolve the tension~-- which could point to radical extensions of the standard model. Here we demonstrate the potential of next-generation spectroscopic galaxy surveys to shed light on the Hubble constant tension. Surveys such as those with Euclid and the Square Kilometre Array (SKA) are expected to reach sub-percent precision on Baryon Acoustic Oscillation (BAO) measurements of the Hubble parameter, with a combined redshift coverage of $0.1<z<3$. This wide redshift range, together with the high precision and low level of systematics in BAO measurements, mean that these surveys will provide independent and tight constraints on $H(z)$. These $H(z)$ measurements can be extrapolated to $z = 0$ to provide constraints on $H_0$ using a non-parametric regression. To this end we deploy Gaussian processes and we find that Euclid-like surveys can reach $\sim$3\% precision on $H_0$, with SKA-like intensity mapping surveys reaching $\sim$2\%. When we combine the low-redshift SKA-like Band 2 survey with either its high-redshift Band 1 counterpart, or with the non-overlapping Euclid-like survey, the precision is predicted to be close to 1\% with 40 $H(z)$ data points. This would be sufficient to rule out the current early- or late-Universe measurements at a $\sim$5$\sigma$ level.

Above the equality scale the galaxy bispectrum will be a key probe for measuring primordial non-Gaussianity which can help differentiate between different inflationary models and other theories of the early universe. On these scales a variety of relativistic effects come into play once the galaxy number-count fluctuation is projected onto our past lightcone. By decomposing the Fourier-space bispectrum into invariant multipoles about the observer's line of sight we examine in detail how the relativistic effects contribute to these. We show how to perform this decomposition analytically, which is significantly faster for subsequent computations. While all multipoles receive a contribution from the relativistic part, odd multipoles arising from the imaginary part of the bispectrum have no Newtonian contribution, making the odd multipoles a smoking gun for a relativistic signature in the bispectrum for single tracers. The dipole and the octopole are significant on equality scales and above where the Newtonian approximation breaks down. This breakdown is further signified by the fact that the even multipoles receive a significant correction on very large scales.

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 beyond-Horndeski gravity has recently been reformulated in the dark energy paradigm - which has been dubbed, Unified Dark Energy (UDE). The evolution equations for the given UDE appear to correspond to a non-conservative dark energy scenario, in which the total energy-momentum tensor is not conserved. We investigate both the background cosmology and, the large-scale imprint of the UDE by probing the angular power spectrum of galaxy number counts, on ultra-large scales; taking care to include the full relativistic corrections in the observed overdensity. The background evolution shows that only an effective mass smaller than the Planck mass is needed in the early universe in order for predictions in the given theory to match current observational constraints. We found that the effective mass-evolution-rate parameter, which drives the evolution of the UDE, acts to enhance the observed power spectrum and, hence, relativistic effects (on ultra-large scales) by enlarging the UDE sound horizon. Conversely, both the (beyond) Horndeski parameter and the kineticity act to diminish the observed power spectrum, by decreasing the UDE sound horizon. Our results show that, in a universe with UDE, a multi-tracer analysis will be needed to detect the relativistic effects in the large-scale structure. In the light of a multi-tracer analysis, the various relativistic effects hold the potential to distinguish different gravity models. Moreover, while the Doppler effect will remain significant at all epochs and, thus can not be ignored, the integrated Sachs-Wolfe, the time-delay and the potential (difference) effects, respectively, will only become significant at epochs near z=3 and beyond, and may be neglected at late epochs. In the same vein, the Doppler effect alone can serve as an effective cosmological probe for the large-scale structure or gravity models, in the angular power spectrum - at all z.

The large-scale structure growth index $\gamma$ provides a consistency test of the standard cosmology and is a potential indicator of modified gravity. We investigate the constraints on $\gamma$ from next-generation spectroscopic surveys, using the power spectrum that is observed in redshift space, i.e., the angular power spectrum. The angular power spectrum avoids the need for an Alcock-Packzynski correction. It also naturally incorporates cosmic evolution and wide-angle effects, without any approximation. We include the cross-correlations between redshift bins, using a hybrid approximation when the total number of bins is computationally unfeasible. We show that the signal-to-noise on $\gamma$ increases as the redshift bin-width is decreased. Noise per bin also increases -- but this is compensated by the increased number of auto- and cross-spectra. In our forecasts, we marginalise over the amplitude of primordial fluctuations and other standard cosmological parameters, including the dark energy equation of state parameter, as well as the clustering bias. Neglecting cross-bin correlations increases the errors by $\sim40 - 150\%$. Using only linear scales, we find that a DESI-like BGS survey and an HI intensity mapping survey with the SKA1 precursor MeerKAT deliver similar errors of $\sim4-6\%$, while a Euclid-like survey and an SKA1 intensity mapping survey give $\sim3\%$ errors.

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 Cosmological Principle states that the Universe is statistically isotropic and homogeneous on large scales. In particular, this implies statistical isotropy in the galaxy distribution, after removal of a dipole anisotropy due to the observer's motion. We test this hypothesis with number count maps from the NVSS radio catalogue. We use a local variance estimator based on patches of different angular radii across the sky and compare the source count variance between and within these patches. In order to assess the statistical significance of our results, we simulate radio maps with the NVSS specifications and mask. We conclude that the NVSS data is consistent with statistical isotropy.

Next-generation galaxy surveys will be able to measure perturbations on scales beyond the equality scale. On these ultra-large scales, primordial non-Gaussianity leaves signatures that can shed light on the mechanism by which perturbations in the early Universe are generated. We perform a forecast analysis for constraining local type non-Gaussianity and its two-parameter extension with a simple scale-dependence. We combine different clustering measurements from future galaxy surveys -- a 21cm intensity mapping survey and two photometric galaxy surveys -- via the multi-tracer approach. Furthermore we then include CMB lensing from a CMB Stage 4 experiment in the multi-tracer, which can improve the constraints on bias parameters. We forecast $\sigma(f_{\rm NL}) \simeq 0.9$ (1.4) by combining SKA1, a Euclid-like (LSST-like) survey, and CMB-S4 lensing. With CMB lensing, the precision on $f_{\rm NL}$ improves by up to a factor of 2, showing that a joint analysis is important. In the case with running of $f_{\rm NL}$, our results show that the combination of upcoming cosmological surveys could achieve $\sigma(n_{\rm NL}) \simeq 0.12$ (0.22) on the running index.

The apparent sizes and brightnesses of galaxies are correlated in a dipolar pattern around matter overdensities in redshift space, appearing larger on their near side and smaller on their far side. The opposite effect occurs for galaxies around an underdense region. These patterns of apparent magnification induce dipole and higher multipole terms in the cross-correlation of galaxy number density fluctuations with galaxy size/brightness (which is sensitive to the convergence field). This provides a means of directly measuring peculiar velocity statistics at low and intermediate redshift, with several advantages for performing cosmological tests of GR. In particular, it does not depend on empirically-calibrated scaling relations like the Tully-Fisher and Fundamental Plane methods. We show that the next generation of spectroscopic galaxy redshift surveys will be able to measure the Doppler magnification effect with sufficient signal-to-noise to test GR on large scales. We illustrate this with forecasts for the constraints that can be achieved on parametrised deviations from GR for forthcoming low-redshift galaxy surveys with DESI and SKA2. Although the cross-correlation statistic considered has a lower signal to noise than RSD, it will be a useful probe of GR since it is sensitive to different systematics.

We use growth of structure data to constrain the effective field theory of dark energy. Considering as case study Horndeski theories with the speed of gravitational waves equal to that of light, we show how constraints on the free parameters and the large-scale structure phenomenological functions can be improved by two ingredients: firstly by complementing the set of redshift-space distortions data with the three recent measurements of the growth rate $f$ and the amplitude of matter fluctuations $\sigma_8$ from the VIPERS and SDSS collaborations; secondly by applying a local Solar System bound on the variation of the Newton constant. This analysis allows us to conclude that: $i)$ despite firmly restricting the predictions of weaker gravity, the inclusion of the Solar System bound does not prevent suppressed growth relative to the standard model $\Lambda$CDM at low redshifts; $ii)$ the same bound in conjunction with the growth of structure data strongly restricts the redshift evolution of the gravitational slip parameter to be close to unity and the present value is constrained to one at the $10^{-3}$ level; $iii)$ the growth of structure data favours a fifth force contribution to the effective gravitational coupling at low redshifts and at more than one sigma at present time.

We propose a method to probe the homogeneity of a general universe, without assuming symmetry. We show that isotropy can be tested at remote locations on the past lightcone by comparing the line-of-sight and transverse expansion rates, using the time dependence of the polarization of Cosmic Microwave Background photons that have been inverse-Compton scattered by the hot gas in massive clusters of galaxies. This probes a combination of remote transverse and parallel components of the expansion rate of the metric, and we may use radial baryon acoustic oscillations or cosmic clocks to measure the parallel expansion rate. Thus we can test remote isotropy, which is a key requirement of a homogeneous universe. We provide explicit formulas that connect observables and properties of the metric.

The local galaxy bias formalism relies on the energy constraint equation at the formation time to relate the metric perturbation to the matter density contrast. In the Newtonian approximation, this relationship is linear, which allows us to specify the initial galaxy density as a function of local physical operators. In general relativity however, the relationship is intrinsically nonlinear and a modulation of the short-wavelength mode by the long-wavelength mode might be expected. We describe in detail how to obtain local coordinates where the coupling of the long- to the short-wavelength modes is removed through a change of coordinates (in the absence of the primordial non-Gaussianity). We derive the general-relativistic correction to the galaxy bias expansion at second order. The correction does not come from the modulation of small-scale clustering by the long-wavelength mode; instead, it arises from distortions of the volume element by the long-wavelength mode and it does not lead to new bias parameters.

The bispectrum will play an important role in future galaxy surveys. On large scales it is a key probe for measuring primordial non-Gaussianity which can help differentiate between different inflationary models and other theories of the early universe. On these scales a variety of relativistic effects come into play once the galaxy number-count fluctuation is projected onto our past lightcone. We show for the first time that the leading relativistic correction from these distortions in the galaxy bispectrum generates a significant dipole, mainly from relativistic redshift space distortions. The amplitude of the dipole can be more than 10% of the monopole even on equality scales. Such a dipole is absent in the Newtonian approximation to the redshift space bispectrum, so it offers a clear signature of relativistic effects on cosmological scales in large scale structure.

We compute the redshift-dependent angular bispectrum of galaxy number counts at tree-level, including nonlinear clustering bias and estimating numerically for the first time the effect of redshift space distortions (RSD). We show that for narrow redshift bins the amplitude of nonlinear RSD is comparable with the matter density perturbations. While our numerical results only include terms relevant on sub-horizon scales, the formalism can readily be extended to the full tree-level bispectrum. Our approach does not rely on the flat-sky approximation and it can be easily generalized to different sources by including the appropriate bias expansion. We test the accuracy of Limber approximation for different z-bins. We highlight the subtle but relevant differences in the angular bispectrum of galaxy number counts with respect to CMB, due to the different scale dependence of perturbations. Our formalism can also be directly applied to the angular HI intensity mapping bispectrum.

The dipole anisotropy seen in the {cosmic microwave background radiation} is interpreted as due to our peculiar motion. The Cosmological Principle implies that this cosmic dipole signal should also be present, with the same direction, in the large-scale distribution of matter. Measurement of the cosmic matter dipole constitutes a key test of the standard cosmological model. Current measurements of this dipole are barely above the expected noise and unable to provide a robust test. Upcoming radio continuum surveys with the SKA should be able to detect the dipole at high signal to noise. We simulate number count maps for SKA survey specifications in Phases 1 and 2, including all relevant effects. Nonlinear effects from local large-scale structure contaminate the {cosmic (kinematic)} dipole signal, and we find that removal of radio sources at low redshift ($z\lesssim 0.5$) leads to significantly improved constraints. We forecast that the SKA could determine the kinematic dipole direction in Galactic coordinates with an error of $(\Delta l,\Delta b)\sim(9^\circ,5^\circ)$ to $(8^\circ, 4^\circ)$, depending on the sensitivity. The predicted errors on the relative speed are $\sim 10\%$. These measurements would significantly reduce the present uncertainty on the direction of the radio dipole, and thus enable the first critical test of consistency between the matter and CMB dipoles.

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.

The galaxy bispectrum on scales around and above the equality scale receives contributions from relativistic effects. Some of these arise from lightcone deformation effects, which come from local and line-of-sight integrated contributions. Here we calculate the local contributions from the generated vector and tensor background which is formed as scalar modes couple and enter the horizon. We show that these modes are sub-dominant when compared with other relativistic contributions.

Yes. In a perturbed Friedmann model, the difference of the Hubble constants measured in two rest-frames is independent of the source peculiar velocity and depends only on the relative velocity of the observers, to lowest order in velocity. Therefore this difference should be zero when averaging over sufficient sources, which are at large enough distances to suppress local nonlinear inhomogeneity. We use a linear perturbative analysis to predict the Doppler effects on redshifts and distances. Since the observed redshifts encode the effect of local bulk flow due to nonlinear structure, our linear analysis is able to capture aspects of the nonlinear behaviour. Using the largest available distance compilation from CosmicFlows-3, we find that the data is consistent with simulations based on the concordance model, for sources at $20-150\,$Mpc.

The late-time integrated Sachs-Wolfe (ISW) signal in the CMB temperature anisotropies is an important probe of dark energy when it can be detected by cross-correlation with large-scale structure surveys. Because of their huge sky area, surveys in the radio are well-suited to ISW detection. We show that 21cm intensity mapping and radio continuum surveys with the SKA in Phase~1 promise a $\sim5\sigma$ detection if we use tomography, with a similar forecast for the precursor EMU survey. In SKA Phase~2, the 21cm galaxy redshift survey and the continuum survey could deliver a $\sim6\sigma$ detection. Our analysis of the radio surveys aims for theoretical accuracy on large scales. Firstly, we include all the effects on the radio surveys from observing on the past lightcone: redshift-space distortions and lensing magnification can have a significant impact on the ISW signal to noise ratio, while Doppler and other relativistic distortions are not significant. Secondly, we use the full information in the observable galaxy angular power spectra $C_\ell(z,z')$, by avoiding the Limber approximation and by including all cross-correlations between redshift bins in the covariance. Without these cross-bin correlations, the ISW signal to noise ratio is biased.

Our velocity relative to the cosmic microwave background (CMB) generates a dipole from the CMB monopole, which was accurately measured by COBE. The relative velocity also modulates and aberrates the CMB fluctuations, generating a small signature of statistical isotropy violation in the covariance matrix. This signature was first measured by Planck 2013. Galaxy surveys are similarly affected by a Doppler boost. The dipole generated from the number count monopole has been extensively discussed, and measured (at very low accuracy) in the NVSS and TGSS radio continuum surveys. For the first time, we present an analysis of the Doppler imprint on the number count fluctuations, using the bipolar spherical harmonic formalism to quantify these effects. Next-generation wide-area surveys with a high redshift range are needed to detect the small Doppler signature in number count fluctuations. We show that radio continuum surveys with the SKA should enable a detection at $\gtrsim 3 \sigma$ in Phase 2, with marginal detection possible in Phase 1.

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.

We study the prospects to measure the cosmic radio dipole by means of continuum surveys with the Square Kilometre Array. Such a measurement will allow a critical test of the cosmological principle. It will test whether the cosmic rest frame defined by the cosmic microwave background at photon decoupling agrees with the cosmic rest frame of matter at late times.

We perform for the first time N-body simulations of Interacting Dark Energy assuming non-Gaussian initial conditions, with the aim of investigating possible degeneracies of these two theoretically independent phenomena in different observational probes. We focus on the large-scale matter distribution, as well as on the statistical and structural properties of collapsed halos and cosmic voids. On very large scales, we show that it is possible to choose the Interaction and non-Gaussian parameters such that their effects on the halo power spectrum cancel, and the power spectrum is indistinguishable from a $\mathrm{\Lambda CDM}$ model. On small scales, measurements of the non-linear matter power spectrum, halo-matter bias, halo and subhalo mass function and cosmic void number function validate the degeneracy determined on large scales. However, the internal structural properties of halos and cosmic voids, namely halo concentration-mass relation and void density profile, are very different from those measured in the $\mathrm{\Lambda CDM}$ model, thereby breaking the degeneracy. In practice, the values of $\mathrm{f_{NL}}$ required to cancel the effect of interaction are already ruled by observations. Our results show in principle that the combination of large- and small-scale probes is needed to constrain Interacting Dark Energy and Primordial non-Gaussianity separately.

Using the consistency relation in Fourier space, we derive the observed galaxy bispectrum from single-field inflation in the squeezed limit, in which one of the three modes has a wavelength much longer than the other two. This provides a non-trivial check of the full computation of the bispectrum based on second-order cosmological perturbation theory in this limit. We show that gauge modes need to be carefully removed in the second-order cosmological perturbations in order to calculate the observed galaxy bispectrum in the squeezed limit. We then give an estimate of the effective non-Gaussianity due to general relativistic lightcone effects that could mimic a primordial non-Gaussian signal.

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

According to the Cosmological Principle, the matter distribution on very large scales should have a kinematic dipole that is aligned with that of the CMB. We determine the dipole anisotropy in the number counts of two all-sky surveys of radio galaxies. For the first time, this analysis is presented for the TGSS survey, allowing us to check consistency of the radio dipole at low and high frequencies by comparing the results with the well-known NVSS survey. We match the flux thresholds of the catalogues, with flux limits chosen to minimise systematics, and adopt a strict masking scheme. We find dipole directions that are in good agreement with each other and with the CMB dipole. In order to compare the amplitude of the dipoles with theoretical predictions, we produce sets of lognormal realisations. Our realisations include the theoretical kinematic dipole, galaxy clustering, Poisson noise, simulated redshift distributions which fit the NVSS and TGSS source counts, and errors in flux calibration. The measured dipole for NVSS is $\sim\!2$ times larger than predicted by the mock data. For TGSS, the dipole is almost $\sim\! 5$ times larger than predicted, even after checking for completeness and taking account of errors in source fluxes and in flux calibration. Further work is required to understand the nature of the systematics that are the likely cause of the anomalously large TGSS dipole amplitude.

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

We 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 slow-roll phase. We forecast constraints with a Fisher analysis on the amplitude of the parametrized features on ultra-large 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 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 plane-parallel 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.

Cosmology relies on the Cosmological Principle, i.e., the hypothesis that the Universe is homogeneous and isotropic on large scales. This implies in particular that the counts of galaxies should approach a homogeneous scaling with volume at sufficiently large scales. Testing homogeneity is crucial to obtain a correct interpretation of the physical assumptions underlying the current cosmic acceleration and structure formation of the Universe. In this Letter, we use the Baryon Oscillation Spectroscopic Survey to make the first spectroscopic and model-independent measurements of the angular homogeneity scale $\theta_{\rm h}$. Applying four statistical estimators, we show that the angular distribution of galaxies in the range 0.46 < z < 0.62 is consistent with homogeneity at large scales, and that $\theta_{\rm h}$ varies with redshift, indicating a smoother Universe in the past. These results are in agreement with the foundations of the standard cosmological paradigm.

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 rest-frame. 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 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).

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 cross-correlation 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 redshift-space 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 Friedmann-Lema\^itre-Robertson-Walker 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 second-order 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 dark-energy 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 sub-percent level if an optimal distance indicator is used.

Next-generation galaxy surveys will increasingly rely on the galaxy bispectrum to improve cosmological constraints, especially on primordial non-Gaussianity. 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. Mode-coupling contributions at second order mean that relativistic corrections can be non-negligible 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 non-Gaussianity. 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 non-Gaussianity.

It is usually assumed that we will need to wait until next-generation surveys like Euclid, LSST and SKA, in order to improve on the current best constraints on primordial non-Gaussianity 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 multi-tracer 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 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.

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and to study radio transients and pulsars. HIRAX will observe unresolved sources of neutral hydrogen via their redshifted 21-cm emission line (`hydrogen intensity mapping'). The resulting maps of large-scale structure at redshifts 0.8-2.5 will be used to measure Baryon Acoustic Oscillations (BAO). HIRAX will improve upon current BAO measurements from galaxy surveys by observing a larger cosmological volume (larger in both survey area and redshift range) and by measuring BAO at higher redshift when the expansion of the universe transitioned to Dark Energy domination. HIRAX will complement CHIME, a hydrogen intensity mapping experiment in the Northern Hemisphere, by completing the sky coverage in the same redshift range. HIRAX's location in the Southern Hemisphere also allows a variety of cross-correlation measurements with large-scale structure surveys at many wavelengths. Daily maps of a few thousand square degrees of the Southern Hemisphere, encompassing much of the Milky Way galaxy, will also open new opportunities for discovering and monitoring radio transients. The HIRAX correlator will have the ability to rapidly and eXperimentciently detect transient events. This new data will shed light on the poorly understood nature of fast radio bursts (FRBs), enable pulsar monitoring to enhance long-wavelength gravitational wave searches, and provide a rich data set for new radio transient phenomena searches. This paper discusses the HIRAX instrument, science goals, and current status.

We present the galaxy number overdensity up to second order in redshift space on cosmological scales for a concordance model. The result contains all general relativistic effects up to second order that arise from observing on the past light cone, including all redshift effects, lensing distortions from convergence and shear, and contributions from velocities, Sachs-Wolfe, integrated SW and time-delay terms. This result will be important for accurate calculation of the bias on estimates of non-Gaussianity and on precision parameter estimates, introduced by nonlinear projection effects.

We present a detailed derivation of the observed galaxy number over-density on cosmological scales up to second order in perturbation theory. We include all relativistic effects that arise from observing on the past lightcone. The derivation is in a general gauge, and applies to all dark energy models (including interacting dark energy) and many modified gravity models. The result will be important for accurate cosmological parameter estimation, including non-Gaussianity, since all projection effects need to be taken into account. It also offers the potential for new probes of General Relativity, dark energy and modified gravity. This paper accompanies Paper I which presents the key results for the concordance model in Poisson gauge.

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

In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is real and physical. However, the variance smoothed on a local physical scale is not correlated with the large-scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that encode this relativistic signature.

Intensity mapping of the neutral hydrogen brightness temperature promises to provide a three-dimensional view of the universe on very large scales. Nonlinear effects are typically thought to alter only the small-scale power, but we show how they may bias the extraction of cosmological information contained in the power spectrum on ultra-large 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 second-order 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 ultra-large 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 high-precision observations and new theoretical ideas are fulfilled.

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

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

We discuss the question of gauge choice when analysing relativistic density perturbations at second order. We compare Newtonian and General Relativistic approaches. Some misconceptions in the recent literature are addressed. We show that the comoving-synchronous gauge is the unique gauge in General Relativity that corresponds to the Lagrangian frame and is entirely appropriate to describe the matter overdensity at second order. The comoving-synchronous gauge is the simplest gauge in which to describe Lagrangian bias at second order.

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

We 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 non-commutativity 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 second-order 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 second-order corrections to the background distance-redshift 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 second-order lensing corrections.

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.

The Square Kilometre Array (SKA) will conduct the biggest spectroscopic galaxy survey ever, by detecting the 21cm emission line of neutral hydrogen (HI) from around a billion galaxies over 3/4 of the sky, out to a redshift of z~2. This will allow the redshift-space matter power spectrum, and corresponding dark energy observables, to be measured with unprecedented precision. In this paper, we present an improved model of the HI galaxy number counts and bias from semi-analytic simulations, and use it to calculate the expected yield of HI galaxies from surveys with a variety of Phase 1 and 2 SKA configurations. We illustrate the relative performance of the different surveys by forecasting errors on the radial and transverse scales of the baryon acoustic oscillation (BAO) feature, finding that the full "billion galaxy survey" with SKA2 will deliver the largest dark energy figure of merit of any current or future large-scale structure survey.

The observed galaxy power spectrum acquires relativistic corrections from lightcone effects, and these corrections grow on very large scales. Future galaxy surveys in optical, infrared and radio bands will probe increasingly large wavelength modes and reach higher redshifts. In order to exploit the new data on large scales, an accurate analysis requires inclusion of the relativistic effects. This is especially the case for primordial non-Gaussianity and for extending tests of dark energy models to horizon scales. Here we investigate the latter, focusing on models where the dark energy interacts non-gravitationally with dark matter. Interaction in the dark sector can also lead to large-scale deviations in the power spectrum. If the relativistic effects are ignored, the imprint of interacting dark energy will be incorrectly identified and thus lead to a bias in constraints on interacting dark energy on very large scales.

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.

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.

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 new frontier of cosmology will be led by three-dimensional surveys of the large-scale structure of the Universe. Based on its all-sky surveys and redshift depth, the SKA is destined to revolutionize cosmology, in combination with future optical/ infrared surveys such as Euclid and LSST. Furthermore, we will not have to wait for the full deployment of the SKA in order to see transformational science. In the first phase of deployment (SKA1), all-sky HI intensity mapping surveys and all-sky continuum surveys are forecast to be at the forefront on the major questions of cosmology. We give a broad overview of the major contributions predicted for the SKA. The SKA will not only deliver precision cosmology -- it will also probe the foundations of the standard model and open the door to new discoveries on large-scale features of the Universe.

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

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

The success of precision cosmology depends not only on accurate observations, but also on the theoretical model - which must be understood to at least the same level of precision. Subtle relativistic effects can lead to biased measurements if they are neglected. One such effect gives a systematic shift in the distance-redshift relation away from its background value, due to the non- linear relativistic conservation of total photon flux. We also show directly how this shift follows from a fully relativistic analysis of the geodesic deviation equation. We derive the expectation value of the shift using second-order perturbations about a concordance background, and show that the distance to last scattering is increased by 1%. We argue that neglecting this shift could lead to a significant bias in the background cosmological parameters, because it alters the meaning of the background model. A naive adjustment of CMB parameter estimation if this shift is really a correction to the background would raise the H0 value inferred from the CMB by 5%, potentially removing the tension with local measurements of H0. Other CMB parameters which depend on the distance would also be shifted by ~1sigma when combined with local H0 data. While our estimations rely on a simplistic analysis, they nevertheless illustrate that accurately defining the background model in terms of the expectation values of observables is critical when we aim to determine the model parameters at the sub-percent level.

If dark energy and dark matter interact via exchange of energy and momentum, then this may affect the galaxy power spectrum on large scales. When this happens, it may be degenerate with the signal from primordial non-Gaussianity via scale-dependent bias. We consider a class of interacting dark energy models and show that the matter overdensity is scale-dependent on large scales. We estimate the effective non-Gaussianity arising from the large-scale effects of interaction in the dark sector. The signal of dark sector interaction can be disentangled from a primordial non-Gaussian signal by measuring the power at two redshifts.

We present a derivation of the cosmological distance-redshift relation up to second order in perturbation theory. In addition, we find the observed redshift and the lensing magnification to second order. We do not require that the density contrast is small, we only that the metric potentials and peculiar velocities are small. Thus our results apply into the nonlinear regime, and can be used for most dark energy models. We present the results in a form which can be readily computed in an N-body simulation. This paper accompanies arXiv:1207.2109, where the key results are summarised in a physically transparent form and applications are discussed.

The next generation of telescopes will usher in an era of precision cosmology, capable of determining the cosmological model to beyond the percent level. For this to be effective, the theoretical model must be understood to at least the same level of precision. A range of subtle relativistic effects remain to be explored theoretically, and offer the potential for probing general relativity in this new regime. We present the distance-redshift relation to second order in cosmological perturbation theory for a general dark energy model. This relation determines the magnification of sources at high precision, as well as redshift space distortions in the mildly non-linear regime. We identify a range of new lensing effects, including: double-integrated and nonlinear integrated Sach-Wolfe contributions, transverse Doppler effects, lensing from the induced vector mode and gravitational wave backgrounds, in addition to lensing from the second-order potential. Modifications to Doppler lensing from redshift-space distortions are identified. Finally, we find a new double-coupling between the density fluctuations integrated along the line of sight, and gradients in the density fluctuations coupled to transverse velocities along the line of sight. These can be large and thus offer important new probes of gravitational lensing and general relativity. This paper accompanies arXiv:1402.1933, where a comprehensive derivation is given.

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

Doppler lensing is the apparent change in object size and magnitude due to peculiar velocities. Objects falling into an overdensity appear larger on its near side, and smaller on its far side, than typical objects at the same redshifts. This effect dominates over the usual gravitational lensing magnification at low redshift. Doppler lensing is a promising new probe of cosmology, and we explore in detail how to utilize the effect with forthcoming surveys. We present cosmological simulations of the Doppler and gravitational lensing effects based on the Millennium simulation. We show that Doppler lensing can be detected around stacked voids or unvirialised over-densities. New power spectra and correlation functions are proposed which are designed to be sensitive to Doppler lensing. We consider the impact of gravitational lensing and intrinsic size correlations on these quantities. We compute the correlation functions and forecast the errors for realistic forthcoming surveys, providing predictions for constraints on cosmological parameters. Finally, we demonstrate how we can make 3-D potential maps of large volumes of the Universe using Doppler lensing.

Quintessence can cluster only on horizon scales. What is the effect on the observed matter distribution? To answer this, we need a relativistic approach that goes beyond the standard Newtonian calculation and deals properly with large scales. Such an approach has recently been developed for the case when dark energy is vacuum energy, which does not cluster at all. We extend this relativistic analysis to deal with dynamical dark energy. Using three quintessence potentials as examples, we compute the angular power spectrum for the case of an HI intensity map survey. Compared to the concordance model with the same small-scale power at z=0, quintessence boosts the angular power by up to ~15% at high redshifts, while power in the two models converges at low redshifts. The difference is mainly due to the background evolution, driven mostly by the normalization of the power spectrum today. The dark energy perturbations make only a small contribution on the largest scales, and a negligible contribution on smaller scales. Ironically, the dark energy perturbations remove the false boost of large-scale power that arises if we impose the (unphysical) assumption that the dark energy is smooth.

If the B-mode signal in the CMB polarization seen by the BICEP2 experiment is confirmed, it has dramatic implications for models of inflation. The result is also in tension with Planck limits on standard inflationary models. It is therefore important to investigate whether this signal can arise from alternative sources. If so, this could lessen the pressure on inflationary models and the tension with Planck data. We investigate whether vector and tensor modes from primordial magnetic fields can explain the signal. We find that in principle, magnetic fields generated during inflation can indeed produce the required B-mode, for a suitable range of energy scales of inflation. In this case, the primordial gravitational wave amplitude is negligible, so that there is no tension with Planck and no problems posed for current inflationary models. However, the simplest magnetic model is in tension with Planck limits on non-Gaussianity in the trispectrum. It may be possible to fine-tune the magnetogenesis model so that this non-Gaussianity is suppressed. Alternatively, a weaker magnetic field can pass the non-Gaussianity constraints and allow the primordial tensor mode to be reduced to r~0.09, thus removing the tension with Planck data and alleviating the problems with simple inflationary models.

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

Using a sample of 608 Type Ia supernovae from the SDSS-II and BOSS surveys, combined with a sample of foreground galaxies from SDSS-II, we estimate the weak lensing convergence for each supernova line-of-sight. We find that the correlation between this measurement and the Hubble residuals is consistent with the prediction from lensing (at a significance of 1.7sigma. Strong correlations are also found between the residuals and supernova nuisance parameters after a linear correction is applied. When these other correlations are taken into account, the lensing signal is detected at 1.4sigma. We show for the first time that distance estimates from supernovae can be improved when lensing is incorporated by including a new parameter in the SALT2 methodology for determining distance moduli. The recovered value of the new parameter is consistent with the lensing prediction. Using CMB data from WMAP7, H0 data from HST and SDSS BAO measurements, we find the best-fit value of the new lensing parameter and show that the central values and uncertainties on Omega_m and w are unaffected. The lensing of supernovae, while only seen at marginal significance in this low redshift sample, will be of vital importance for the next generation of surveys, such as DES and LSST, which will be systematics dominated.

The standard concordance model of the Universe is based on the cosmological constant as the driver of accelerating expansion. This concordance model is being subjected to a growing range of inter-locking observations. In addition to using generic observational tests, one can also design tests that target the specific properties of the cosmological constant. These null tests do not rely on parametrisations of observables, but focus on quantities that are constant only if dark energy is a cosmological constant. We use supernova data in null tests that are based on the luminosity distance. In order to extract derivatives of the distance in a model-independent way, we use Gaussian Processes. We find that the concordance model is compatible with the Union 2.1 data, but the error bars are fairly large. Simulated datasets are generated for the DES supernova survey and we show that this survey will allow for a sharper null test of the cosmological constant if we assume the Universe is flat. Allowing for spatial curvature degrades the power of the null test.

A suggested solution to the dark energy problem is the void model, where accelerated expansion is replaced by Hubble-scale inhomogeneity. In these models, density perturbations grow on a radially inhomogeneous background. This large scale inhomogeneity distorts the spherical Baryon Acoustic Oscillation feature into an ellipsoid which implies that the bump in the galaxy correlation function occurs at different scales in the radial and transverse correlation functions. We compute these for the first time, under the approximation that curvature gradients do not couple the scalar modes to vector and tensor modes. The radial and transverse correlation functions are very different from those of the concordance model, even when the models have the same average BAO scale. This implies that if void models are fine-tuned to satisfy average BAO data, there is enough extra information in the correlation functions to distinguish a void model from the concordance model. We expect these new features to remain when the full perturbation equations are solved, which means that the radial and transverse galaxy correlation functions can be used as a powerful test of the Copernican Principle.

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

More than half of the volume of our Universe is occupied by cosmic voids. The lensing magnification effect from those under-dense regions is generally thought to give a small dimming contribution: objects on the far side of a void are supposed to be observed as slightly smaller than if the void were not there, which together with conservation of surface brightness implies net reduction in photons received. This is predicted by the usual weak lensing integral of the density contrast along the line of sight. We show that this standard effect is swamped at low redshifts by a relativistic Doppler term that is typically neglected. Contrary to the usual expectation, objects on the far side of a void are brighter than they would be otherwise. Thus the local dynamics of matter in and near the void is crucial and is only captured by the full relativistic lensing convergence. There are also significant nonlinear corrections to the relativistic linear theory, which we show actually under-predicts the effect. We use exact solutions to estimate that these can be more than 20% for deep voids. This remains an important source of systematic errors for weak lensing density reconstruction in galaxy surveys and for supernovae observations, and may be the cause of the reported extra scatter of field supernovae located on the edge of voids compared to those in clusters.

Observationally confirming spatial homogeneity on sufficiently large cosmological scales is of importance to test one of the underpinning assumptions of cosmology, and is also imperative for correctly interpreting dark energy. A challenging aspect of this is that homogeneity must be probed inside our past lightcone, while observations take place on the lightcone. The star formation history (SFH) in the galaxy fossil record provides a novel way to do this. We calculate the SFH of stacked Luminous Red Galaxy (LRG) spectra obtained from the Sloan Digital Sky Survey. We divide the LRG sample into 12 equal area contiguous sky patches and 10 redshift slices (0.2 < z < 0.5), which correspond to 120 blocks of volume 0.04Gpc3. Using the SFH in a time period which samples the history of the Universe between look-back times 11.5 to 13.4 Gyrs as a proxy for homogeneity, we calculate the posterior distribution for the excess large-scale variance due to inhomogeneity, and find that the most likely solution is no extra variance at all. At 95% credibility, there is no evidence of deviations larger than 5.8%.

Physical (and weak) regularity conditions are used to determine and classify all the possible types of spherically symmetric dust spacetimes in general relativity. This work unifies and completes various earlier results. The junction conditions are described for general non-comoving (and non-null) surfaces, and the limits of kinematical quantities are given on all comoving surfaces where there is Darmois matching. We show that an inhomogeneous generalisation of the Kantowski-Sachs metric may be joined to the Lemaitre-Tolman-Bondi metric. All the possible spacetimes are explicitly divided into four groups according to topology, including a group in which the spatial sections have the topology of a 3-torus. The recollapse conjecture (for these spacetimes) follows naturally in this approach.

Direct observations of the Hubble rate, from cosmic chronometers and the radial baryon acoustic oscillation scale, can out-perform supernovae observations in understanding the expansion history, because supernovae observations need to be differentiated to extract H(z). We use existing H(z) data and smooth the data using a new Gaussian Processes package, GaPP, from which we can also estimate derivatives. The obtained Hubble rate and its derivatives are used to reconstruct the equation of state of dark energy and to perform consistency tests of the LCDM model, some of which are newly devised here. Current data is consistent with the concordance model, but is rather sparse. Future observations will provide a dramatic improvement in our ability to constrain or refute the concordance model of cosmology. We produce simulated data to illustrate how effective H(z) data will be in combination with Gaussian Processes.

Light from `point sources' such as supernovae is observed with a beam width of order of the sources' size - typically less than 1 AU. Such a beam probes matter and curvature distributions that are very different from coarse-grained representations in N-body simulations or perturbation theory, which are smoothed on scales much larger than 1 AU. The beam typically travels through unclustered dark matter and hydrogen with a mean density much less than the cosmic mean, and through dark matter halos and hydrogen clouds. Using N-body simulations, as well as a Press-Schechter approach, we quantify the density probability distribution as a function of beam width and show that, even for Gpc-length beams of 500 kpc diameter, most lines of sight are significantly under-dense. From this we argue that modelling the probability distribution for AU-diameter beams is absolutely critical. Standard analyses predict a huge variance for such tiny beam sizes, and nonlinear corrections appear to be non-trivial. It is not even clear whether under-dense regions lead to dimming or brightening of sources, owing to the uncertainty in modelling the expansion rate which we show is the dominant contribution. By considering different reasonable approximations which yield very different cosmologies we argue that modelling ultra-narrow beams accurately remains a critical problem for precision cosmology. This could appear as a discordance between angular diameter and luminosity distances when comparing SN observations to BAO or CMB distances.

Local non-Gaussianity, parametrized by $f_{\rm NL}$, introduces a scale-dependent bias that is strongest at large scales, precisely where General Relativistic (GR) effects also become significant. With future data, it should be possible to constrain $f_{\rm NL} = {\cal O}(1)$ with high redshift surveys. GR corrections to the power spectrum and ambiguities in the gauge used to define bias introduce effects similar to $f_{\rm NL}= {\cal O}(1)$, so it is essential to disentangle these effects. For the first time in studies of primordial non-Gaussianity, we include the consistent GR calculation of galaxy power spectra, highlighting the importance of a proper definition of bias. We present observable power spectra with and without GR corrections, showing that an incorrect definition of bias can mimic non-Gaussianity. However, these effects can be distinguished by their different redshift and scale dependence, so as to extract the true primordial non-Gaussianity.

In standard cosmologies, dark energy interacts only gravitationally with dark matter. There could be a non-gravitational interaction in the dark sector, leading to changes in the effective DE equation of state, in the redshift dependence of the DM density and in structure formation. We use CMB, BAO and SNIa data to constrain a model where the energy transfer in the dark sector is proportional to the DE density. There are two subclasses, defined by the vanishing of momentum transfer either in the DM or the DE frame. We conduct a Markov-Chain Monte-Carlo analysis to obtain best-fit parameters. The background evolution allows large interaction strengths, and the constraints from CMB anisotropies are weak. The growth of DM density perturbations is much more sensitive to the interaction, and can deviate strongly from the standard case. However, the deviations are degenerate with galaxy bias and thus more difficult to constrain. Interestingly, the ISW signature is suppressed since the non-standard background evolution can compensate for high growth rates. We also discuss the partial degeneracy between interacting DE and modified gravity, and how this can be broken.

The standard model of cosmology is based on the existence of homogeneous surfaces as the background arena for structure formation. Homogeneity underpins both general relativistic and modified gravity models and is central to the way in which we interpret observations of the CMB and the galaxy distribution. However, homogeneity cannot be directly observed in the galaxy distribution or CMB, even with perfect observations, since we observe on the past lightcone and not on spatial surfaces. We can directly observe and test for isotropy, but to link this to homogeneity, we need to assume the Copernican Principle. First, we discuss the link between isotropic observations on the past lightcone and isotropic spacetime geometry: what observations do we need to be isotropic in order to deduce spacetime isotropy? Second, we discuss what we can say with the Copernican assumption. The most powerful result is based on the CMB: the vanishing of the dipole, quadrupole and octupole of the CMB is sufficient to impose homogeneity. Real observations lead to near-isotropy on large scales - does this lead to near-homogeneity? There are important partial results, and we discuss why this remains a difficult open question. Thus we are currently unable to prove homogeneity of the Universe on large-scales, even with the Copernican Principle. However we can use observations of the CMB, galaxies and clusters to test homogeneity itself.

The standard Friedmann model of cosmology is based on the Copernican Principle, i.e. the assumption of a homogeneous background on which structure forms via perturbations. Homogeneity underpins both general relativistic and modified gravity models and is central to the way in which we interpret observations of the CMB and the galaxy distribution. It is therefore important to probe homogeneity via observations. We describe a test based on the fossil record of distant galaxies: if we can reconstruct key intrinsic properties of galaxies as functions of proper time along their worldlines, we can compare such properties at the same proper time for our galaxy and others. We achieve this by computing the lookback time using radial Baryon Acoustic Oscillations, and the time along galaxy world line using stellar physics, allowing us to probe homogeneity, in principle anywhere inside the past light cone. Agreement in the results would be an important consistency test -- although it would not in itself prove homogeneity. Any significant deviation in the results however would signal a breakdown of homogeneity.

Nonlinear dynamics creates vortical currents when the tight-coupling approximation between photons and baryons breaks down around the time of recombination. This generates a magnetic field at second order in cosmological perturbations, whose power spectrum is fixed by standard physics, without the need for any ad hoc assumptions. We present the fully relativistic calculation of the magnetic power spectrum, including the effects of metric perturbations, second-order velocity and the photon anisotropic stress, thus generalizing and correcting previous results. We also show that significant magnetogenesis continues to occur after recombination. The power spectrum $\sqrt{k^3 P_B} $ behaves as $ \propto k^4$ on large scales, and $\propto k^{0.5}$ on small scales, down to $\sim 1\,$Mpc. On cluster scales, the created field has strength $\sim 3\times 10^{-29}$ Gauss.

The observable universe could be a 1+3-surface (the "brane") embedded in a 1+3+\textit{d}-dimensional spacetime (the "bulk"), with Standard Model particles and fields trapped on the brane while gravity is free to access the bulk. At least one of the \textit{d} extra spatial dimensions could be very large relative to the Planck scale, which lowers the fundamental gravity scale, possibly even down to the electroweak ($\sim$ TeV) level. This revolutionary picture arises in the framework of recent developments in M theory. The 1+10-dimensional M theory encompasses the known 1+9-dimensional superstring theories, and is widely considered to be a promising potential route to quantum gravity. At low energies, gravity is localized at the brane and general relativity is recovered, but at high energies gravity "leaks" into the bulk, behaving in a truly higher-dimensional way. This introduces significant changes to gravitational dynamics and perturbations, with interesting and potentially testable implications for high-energy astrophysics, black holes, and cosmology. Brane-world models offer a phenomenological way to test some of the novel predictions and corrections to general relativity that are implied by M theory. This review analyzes the geometry, dynamics and perturbations of simple brane-world models for cosmology and astrophysics, mainly focusing on warped 5-dimensional brane-worlds based on the Randall--Sundrum models. We also cover the simplest brane-world models in which 4-dimensional gravity on the brane is modified at \emph{low} energies -- the 5-dimensional Dvali--Gabadadze--Porrati models. Then we discuss co-dimension two branes in 6-dimensional models.

The apparent accelerating expansion of the Universe is forcing us to examine the foundational aspects of the standard model of cosmology -- in particular, the fact that dark energy is a direct consequence of the homogeneity assumption. We discuss the foundations of the assumption of spatial homogeneity, in the case when the Copernican Principle is adopted. We present results that show how (almost-) homogeneity follows from (almost-) isotropy of various observables. The analysis requires the fully nonlinear field equations -- i.e., it is not possible to use second- or higher-order perturbation theory, since one cannot assume a homogeneous and isotropic background. Then we consider what happens if the Copernican Principle is abandoned in our Hubble volume. The simplest models are inhomogeneous but spherically symmetric universes which do not require dark energy to fit the distance modulus. Key problems in these models are to compute the CMB anisotropies and the features of large-scale structure. We review how to construct perturbation theory on a non-homogeneous cosmological background, and discuss the complexities that arise in using this to determine the growth of large-scale structure.

We study perturbations of a scalar field cosmology in Horava-Lifshitz gravity, adopting the most general setup without detailed balance but with the projectability condition. We derive the generalized Klein-Gordon equation, which is sixth-order in spatial derivatives. Then we investigate scalar field perturbations coupled to gravity in a flat Friedmann-Robertson-Walker background. In the sub-horizon regime, the metric and scalar field modes have independent oscillations with different frequencies and phases except in particular cases. On super-horizon scales, the perturbations become adiabatic during slow-roll inflation driven by a single field, and the comoving curvature perturbation is constant.

In the Horava-Lifshitz theory of quantum gravity, two conditions -- detailed balance and projectability -- are usually assumed. The breaking of projectability simplifies the theory, but it leads to serious problems with the theory. The breaking of detailed balance leads to a more complicated form of the theory, but it appears to resolve some of the problems. Sotiriou, Visser and Weinfurtner formulated the most general theory of Horava-Lifshitz type without detailed balance. We compute the linear scalar perturbations of the FRW model in this form of HL theory. We show that the higher-order curvature terms in the action lead to a gravitational effective anisotropic stress on small scales. Specializing to a Minkowski background, we study the spin-0 scalar mode of the graviton, using a gauge-invariant analysis, and find that it is stable in both the infrared and ultraviolet regimes for $0 \le \xi \le 2/3$. However, in this parameter range the scalar mode is a ghost.

We present a new systematic analysis of the early radiation era solution in an interacting dark energy model to find the adiabatic initial conditions for the Boltzmann integration. In a model where the interaction is proportional to the dark matter density, adiabatic initial conditions and viable cosmologies are possible if the early-time dark energy equation of state parameter is $w_e > -4/5$. We find that when adiabaticity between cold dark matter, baryons, neutrinos and photons is demanded, the dark energy component satisfies automatically the adiabaticity condition. As supernovae Ia or baryon acoustic oscillation data require the recent-time equation of state parameter to be more negative, we consider a time-varying equation of state in our model. In a companion paper [arXiv:0907.4987] we apply the initial conditions derived here, and perform a full Monte Carlo Markov Chain likelihood analysis of this model.

We use observations of cosmic microwave background anisotropies, supernova luminosities and the baryon acoustic oscillation signal in the galaxy distribution to constrain the cosmological parameters in a simple interacting dark energy model with a time-varying equation of state. Using a Monte Carlo Markov Chain technique we determine the posterior likelihoods. Constraints from the individual data sets are weak, but the combination of the three data sets confines the interaction constant $\Gamma$ to be less than 23% of the expansion rate of the Universe $H_0$; at 95% CL $-0.23 < \Gamma/H_0 < +0.15$. The CMB acoustic peaks can be well fitted even if the interaction rate is much larger, but this requires a larger or smaller (depending on the sign of interaction) matter density today than in the non-interacting model. Due to this degeneracy between the matter density and the interaction rate, the only observable effect on the CMB is a larger or smaller integrated Sachs-Wolfe (ISW) effect. While SN or BAO data alone do not set any direct constraints on the interaction, they exclude the models with very large matter density, and hence indirectly constrain the interaction rate when jointly analysed with the CMB data. To enable the analysis described in this paper, we present in a companion paper [arXiv:0907.4981] a new systematic analysis of the early radiation era solution to find the adiabatic initial conditions for the Boltzmann integration.

We introduce a new form of coupling between dark energy and dark matter that is quadratic in their energy densities. Then we investigate the background dynamics when dark energy is in the form of exponential quintessence. The three types of quadratic coupling all admit late-time accelerating critical points, but these are not scaling solutions. We also show that two types of coupling allow for a suitable matter era at early times and acceleration at late times, while the third type of coupling does not admit a suitable matter era.

Dark energy in General Relativity is typically non-interacting with other matter. However, it is possible that the dark energy interacts with the dark matter, and in this case, the dark matter can violate the universality of free fall (the weak equivalence principle). We show that some forms of the dark sector interaction do not violate weak equivalence. For those interactions that do violate weak equivalence, there are no available laboratory experiments to probe this violation for dark matter. But cosmology provides a test for violations of the equivalence principle between dark matter and baryons -- via a test for consistency of the observed galaxy velocities with the Euler equation.

If dark energy interacts with dark matter, there is a change in the background evolution of the universe, since the dark matter density no longer evolves as a^{-3}. In addition, the non-gravitational interaction affects the growth of structure. In principle, these changes allow us to detect and constrain an interaction in the dark sector. Here we investigate the growth factor and the weak lensing signal for a new class of interacting dark energy models. In these models, the interaction generalises the simple cases where one dark fluid decays into the other. In order to calculate the effect on structure formation, we perform a careful analysis of the perturbed interaction and its effect on peculiar velocities. Assuming a normalization to today's values of dark matter density and overdensity, the signal of the interaction is an enhancement (suppression) of both the growth factor and the lensing power, when the energy transfer in the background is from dark matter to dark energy (dark energy to dark matter).

Dark energy and dark matter are only indirectly measured via their gravitational effects. It is possible that there is an exchange of energy within the dark sector, and this offers an interesting alternative approach to the coincidence problem. We consider two broad classes of interacting models where the energy exchange is a linear combination of the dark sector densities. The first class has been previously investigated, but we define new variables and find a new exact solution, which allows for a more direct, transparent and comprehensive analysis. The second class has not been investigated in general form before. We give general conditions on the parameters in both classes to avoid unphysical behavior (such as negative energy densities).

We investigate the spectrum of vector modes today which is generated at second order by density perturbations. The vector mode background that is generated by structure formation is small but in principle it contributes to the integrated Sachs-Wolfe effect, to redshift-space distortions and to weak lensing. We recover, clarify and extend previous results, and explain carefully why no vorticity is generated in the fluid at second order. The amplitude of the induced vector mode in the metric is around 1% that of the first-order scalars on small scales. We also calculate the power spectrum and the energy density of the vector part of the shear at second order.

Explanations of the late-time cosmic acceleration within the framework of general relativity are plagued by difficulties. General relativistic models are mostly based on a dark energy field with fine-tuned, unnatural properties. There is a great variety of models, but all share one feature in common -- an inability to account for the gravitational properties of the vacuum energy, and a failure to solve the so-called coincidence problem. Two broad alternatives to dark energy have emerged as candidate models: these typically address only the coincidence problem and not the vacuum energy problem. The first is based on general relativity and attempts to describe the acceleration as an effect of inhomogeneity in the universe. If this alternative could be shown to work, then it would provide a dramatic resolution of the coincidence problem; however, a convincing demonstration of viability has not yet emerged. The second alternative is based on infra-red modifications to general relativity, leading to a weakening of gravity on the largest scales and thus to acceleration. Most examples investigated so far are scalar-tensor or brane-world models, and we focus on the simplest candidates of each type: $f(R)$ models and DGP models respectively. Both of these provide a new angle on the problem, but they also face serious difficulties. However, investigation of these models does lead to valuable insights into the properties of gravity and structure formation, and it also leads to new strategies for testing the validity of General Relativity itself on cosmological scales.

If dark energy interacts with dark matter, this gives a new approach to the coincidence problem. But interacting dark energy models can suffer from pathologies. We consider the case where the dark energy is modelled as a fluid with constant equation of state parameter w. Non-interacting constant-w models are well behaved in the background and in the perturbed universe. But the combination of constant w and a simple interaction with dark matter leads to an instability in the dark sector perturbations at early times: the curvature perturbation blows up on super-Hubble scales. Our results underline how important it is to carefully analyze the relativistic perturbations when considering models of coupled dark energy. The instability that we find has been missed in some previous work where the perturbations were not consistently treated. The unstable mode dominates even if adiabatic initial conditions are used. The instability also arises regardless of how weak the coupling is. This non-adiabatic instability is different from previously discovered adiabatic instabilities on small scales in the strong-coupling regime.

General relativity marked the beginning of modern cosmology and it has since been at the centre of many of the key developments in this field. In the present review, we discuss the general-relativistic dynamics and perturbations of the standard cosmological model, the Friedmann-Lemaitre universe, and how these can explain and predict the properties of the observable universe. Our aim is to provide an overview of the progress made in several major research areas, such as linear and non-linear cosmological perturbations, large-scale structure formation and the physics of the cosmic microwave background radiation, in view of current and upcoming observations. We do this by using a single formalism throughout the review, the 1+3 covariant approach to cosmology, which allows for a uniform and balanced presentation of technical information and physical insight.

We consider a phenomenological model for a coupling between the dark matter and dark energy fluids and investigate the sensitivity of a weak lensing measurement for constraining the size of this coupling term. Physically, the functional form of the coupling term in our model describes the decay of dark matter into dark energy. We present forecasts for tomographic measurements of the weak shear bispectrum for the DUNE experiment in a Fisher-matrix formalism, where we describe the nonlinearities in structure formation by hyper-extended perturbation theory. Physically, CDM decay tends to increase the growth rate of density perturbations due to higher values for the CDM density at early times, and amplifies the lensing signal because of stronger fluctuations in the gravitational potential. We focus on degeneracies between the dark energy equation of state properties and the CDM decay constant relevant for structure formation and weak lensing. A typical lower bound on the CDM decay time ~7.7/H_0 = 75.3 Gyr/h$ which could be provided by DUNE would imply that it would not possible to produce the dark energy content of the universe by CDM decay within the age of the Universe for a constant equation of state parameter of w close to -1.

We analyse the possibility that topological defects can act as a source of magnetic fields through the Harrison mechanism in the radiation era. We give a detailed relativistic derivation of the Harrison mechanism at first order in cosmological perturbations, and show that it is only efficient for temperatures above T ~ 0.2 keV. Our main result is that the vector metric perturbations generated by the defects cannot induce vorticity in the matter fluids at linear order, thereby excluding the production of currents and magnetic fields. We show that anisotropic stress in the matter fluids is required to source vorticity and magnetic fields. Our analysis is relevant for any mechanism whereby vorticity is meant to be transferred purely by gravitational interactions, and thus would also apply to dark matter or neutrinos.

Dark energy and dark matter are the dominant sources in the evolution of the late universe. They are currently only indirectly detected via their gravitational effects, and there could be a coupling between them without violating observational constraints. We investigate the background dynamics when dark energy is modelled as exponential quintessence, and is coupled to dark matter via simple models of energy exchange. We introduce a new form of dark sector coupling, which leads to a more complicated dynamical phase space and has a better physical motivation than previous mathematically similar couplings.

Loop Quantum Cosmology strongly modifies the high-energy dynamics of Friedman-Robertson-Walker models and removes the big-bang singularity. We investigate how LQC corrections affect the stability properties of the Einstein static universe. In General Relativity, the Einstein static model with positive cosmological constant Lambda is unstable to homogeneous perturbations. We show that LQC modifications can lead to a centre of stability for a large enough positive value of Lambda.

Observations provide increasingly strong evidence that the universe is accelerating. This revolutionary advance in cosmological observations confronts theoretical cosmology with a tremendous challenge, which it has so far failed to meet. Explanations of cosmic acceleration within the framework of general relativity are plagued by difficulties. General relativistic models are nearly all based on a dark energy field with fine-tuned, unnatural properties. There is a great variety of models, but all share one feature in common -- an inability to account for the gravitational properties of the vacuum energy. Speculative ideas from string theory may hold some promise, but it is fair to say that no convincing model has yet been proposed. An alternative to dark energy is that gravity itself may behave differently from general relativity on the largest scales, in such a way as to produce acceleration. The alternative approach of modified gravity (or dark gravity) provides a new angle on the problem, but also faces serious difficulties, including in all known cases severe fine-tuning and the problem of explaining why the vacuum energy does not gravitate. The lack of an adequate theoretical framework for the late-time acceleration of the universe represents a deep crisis for theory -- but also an exciting challenge for theorists. It seems likely that an entirely new paradigm is required to resolve this crisis.

We review spacetime dynamics in the presence of large-scale electromagnetic fields and then consider the effects of the magnetic component on perturbations to a spatially homogeneous and isotropic universe. Using covariant techniques, we refine and extend earlier work and provide the magnetohydrodynamic equations that describe inhomogeneous magnetic cosmologies in full general relativity. Specialising this system to perturbed Friedmann-Robertson-Walker models, we examine the effects of the field on the expansion dynamics and on the growth of density inhomogeneities, including non-adiabatic modes. We look at scalar perturbations and obtain analytic solutions for their linear evolution in the radiation, dust and inflationary eras. In the dust case we also calculate the magnetic analogue of the Jeans length. We then consider the evolution of vector perturbations and find that the magnetic presence generally reduces the decay rate of these distortions. Finally, we examine the implications of magnetic fields for the evolution of cosmological gravitational waves.

A usual causal requirement on a viable theory of matter is that the speed of sound be at most the speed of light. In view of various recent papers querying this limit, the question is revisited here. We point to various issues confronting theories that violate the usual constraint.

In the standard cosmological model, magnetic fields and vorticity are generated during the radiation era via second-order density perturbations. In order to clarify the complicated physics of this second-order magnetogenesis, we use a covariant approach and present the electromagneto-dynamical equations in the nonlinear regime. We use the tight-coupling approximation to analyze Thomson and Coulomb scattering. At the zero-order limit of exact tight-coupling, we show that the vorticity is zero and no magnetogenesis takes place at any nonlinear order. We show that magnetogenesis also fails at all orders if either protons or electrons have the same velocity as the radiation, and momentum transfer is neglected. Then we prove a key no-go result: at first-order in the tight-coupling approximation, magnetic fields and vorticity still cannot be generated even via nonlinear effects. The tight-coupling approximation must be broken at first order, for the generation of vorticity and magnetic fields, and we derive a closed set of nonlinear evolution equations that governs this generation. We estimate that the amplitude of the magnetic field at recombination on the horizon scale is $\sim 10^{-27} $G.

The DGP brane-world model provides a simple alternative to the standard LCDM cosmology, with the same number of parameters. There is no dark energy - the late universe self-accelerates due to an infrared modification of gravity. We compute the joint constraints on the DGP model from supernovae, the cosmic microwave background shift parameter, and the baryon oscillation peak in the SDSS luminous red galaxy sample. Flat DGP models are within the joint 2 sigma contour, but the LCDM model provides a significantly better fit to the data. These tests are based on the background dynamics of the DGP model, and we comment on further tests that involve structure formation.

We investigate a simple braneworld model in which the universe contains only cold dark matter and a cosmological constant, but the effective dark energy is phantom-like because of extra-dimensional gravity effects. Modified gravity screens the cosmological constant Lambda, allowing for a larger Lambda. In practice, observations do not favour any significant screening. We use supernova data, the cosmic microwave background shift parameter, and the baryon oscillation peak in the galaxy distribution to constrain the model. We find the mean value of Omega_m with 68% confidence limits, and an upper limit on Omega_Lambda at the 68% confidence level. The best-fit model is very close to a standard LCDM model, but the LCDM model provides a better fit since it has one less parameter.

A class of braneworld models can lead to phantom-like acceleration of the late universe, but without the need for any phantom matter. In the simplest models, the universe contains only cold dark matter and a cosmological constant. We generalize these models by introducing a quintessence field. The new feature in our models is that quintessence leads to a crossing of the phantom divide, $w=-1$. This is a purely gravitational effect, and there is no phantom instability. Furthermore, the Hubble parameter is always decreasing, and there is no big rip singularity in the future.

Recent observations provide strong evidence that the universe is accelerating. This confronts theory with a severe challenge. Explanations of the acceleration within the framework of general relativity are plagued by difficulties. General relativistic models require a "dark energy" field with effectively negative pressure. An alternative to dark energy is that gravity itself may behave differently from general relativity on the largest scales, in such a way as to produce acceleration. The alternative approach of modified gravity also faces severe difficulties, but does provide a new angle on the problem. This review considers an example of modified gravity, provided by brane-world models that self-accelerate at late times.

The DGP brane-world model provides an alternative to the standard LCDM cosmology, in which the late universe accelerates due to a modification of gravity rather than vacuum energy. The cosmological constant $\Lambda$ in LCDM is replaced by a single parameter, the crossover scale $r_c$, in DGP. The Supernova redshift observations can be fitted by both models, with $\Lambda\sim H_0^2$ and $r_c \sim H_0^{-1}$. This degeneracy is broken by structure formation, which is suppressed in different ways in the two models. There is some confusion in the literature about how the standard linear growth factor is modified in DGP. While the luminosity distance can be computed purely from the modified 4-dimensional Friedman equation, the evolution of density perturbations requires an analysis of the 5-dimensional gravitational field. We show that if the 5-dimensional effects are inappropriately neglected, then the 4-dimensional Bianchi identities are violated and the computed growth factor is incorrect. By using the 5-dimensional equations, we derive the correct growth factor.

Brane-world models offer the possibility of explaining the late acceleration of the universe via infra-red modifications to General Relativity, rather than a dark energy field. However, one also expects ultra-violet modifications to General Relativity, when high-energy stringy effects in the early universe begin to grow. We generalize the DGP brane-world model via an ultra-violet modification, in the form of a Gauss-Bonnet term in the bulk action. The combination of infra-red and ultra-violet modifications produces an intriguing cosmology. The DGP feature of late-time acceleration without dark energy is preserved, but there is an entirely new feature - there is no hot big bang in the early universe. The universe starts with finite density and pressure, from a "sudden" curvature singularity.

Quantum gravity is expected to remove the classical singularity that arises as the end-state of gravitational collapse. To investigate this, we work with a toy model of a collapsing homogeneous scalar field. We show that non-perturbative semi-classical effects of Loop Quantum Gravity cause a bounce and remove the black hole singularity. Furthermore, we find a critical threshold scale, below which no horizon forms -- quantum gravity may exclude very small astrophysical black holes.

If string theory is correct, then our observable Universe may be a 3-dimensional "brane" embedded in a higher-dimensional spacetime. This theoretical scenario should be tested via the state-of-the-art in gravitational experiments -- the current and upcoming gravity-wave detectors. Indeed, the existence of extra dimensions leads to oscillations that leave a spectroscopic signature in the gravity-wave signal from black holes. The detectors that have been designed to confirm Einstein's prediction of gravity waves, can in principle also provide tests and constraints on string theory.

Using the black string between two branes as a model of a brane-world black hole, we compute the gravity wave perturbations and identify the features arising from the additional polarizations of the graviton. The standard four-dimensional gravitational wave signal acquires late-time oscillations due to massive modes of the graviton. The Fourier transform of these oscillations shows a series of spikes associated with the masses of the Kaluza-Klein modes, providing in principle a spectroscopic signature of extra dimensions.

Brillouin scattering of photons off the density fluctuations in a fluid is potentially important for cosmology. We derive the Brillouin spectral distortion of blackbody radiation, and discuss the possible implications for the cosmic microwave background. The thermal Sunyaev-Zeldovich effect is slightly modified by Brillouin distortion, but only at very long wavelengths.

Braneworld inflation is a phenomenology related to string theory that describes high-energy modifications to general relativistic inflation. The observable universe is a braneworld embedded in 5-dimensional anti de Sitter spacetime. Whe the 5-dimensional action is Einstein-Hilbert, we have a Randall-Sundrum type braneworld. The amplitude of tensor and scalar perturbations from inflation is strongly increased relative to the standard results, although the ratio of tensor to scalar amplitudes still obeys the standard consistency relation. If a Gauss-Bonnet term is included in the action, as a high-energy correction motivated by string theory, we show that there are important changes to the Randall-Sundrum case. We give an exact analysis of the tensor perturbations. They satisfy the same wave equation and have the same spectrum as in the Randall-Sundrum case, but the Gauss-Bonnet change to the junction conditions leads to a modified amplitude of gravitational waves. The amplitude is no longer monotonically increasing with energy scale, but decreases asymptotically after an initial rise above the standard level. Using an approximation that neglects bulk effects, we show that the amplitude of scalar perturbations has a qualitatively similar behaviour to the tensor amplitude. In addition, the tensor to scalar ratio breaks the standard consistency relation.

At high energies on a cosmological brane of Randall-Sundrum type, particle interactions can produce gravitons that are emitted into the bulk and that can feed a bulk black hole. We generalize previous investigations of such radiating brane-worlds by allowing for a breaking of Z_2-symmetry, via different bulk cosmological constants and different initial black hole masses on either side of the brane. One of the notable features of asymmetry is a suppression of the asymptotic level of dark radiation, which means that nucleosynthesis constraints are easier to satisfy. There are also models where the radiation escapes to infinity on one or both sides, rather than falling into a black hole, but these models can have negative energy density on the brane.

In loop quantum cosmology, the universe avoids a big bang singularity and undergoes an early and short super-inflation phase. During super-inflation, non-perturbative quantum corrections to the dynamics drive an inflaton field up its potential hill, thus setting the initial conditions for standard inflation. We show that this effect can raise the inflaton high enough to achieve sufficient e-foldings in the standard inflation era. We analyze the cosmological perturbations generated when slow-roll is violated after super-inflation, and show that loop quantum effects can in principle leave an indirect signature on the largest scales in the CMB, with some loss of power and running of the spectral index.

Loop quantum gravity introduces strong non-perturbative modifications to the dynamical equations in the semi-classical regime, which are responsible for various novel effects, including resolution of the classical singularity in a Friedman universe. Here we investigate the modifications for the case of a cyclic universe potential, assuming that we can apply the four-dimensional loop quantum formalism within the effective four-dimensional theory of the cyclic scenario. We find that loop quantum effects can dramatically alter the near-collision dynamics of the cyclic scenario. In the kinetic-dominated collapse era, the scalar field is effectively frozen by loop quantum friction, so that the branes approach collision and bounce back without actual collision.

If the observable universe is a braneworld of Randall-Sundrum type, then particle interactions at high energies will produce 5-dimensional gravitons that escape into the bulk. As a result, the Weyl energy density on the brane does not behave like radiation in the early universe, but does so only later, in the low energy regime. Recently a simple model was proposed to describe this modification of the Randall-Sundrum cosmology. We investigate the dynamics of this model, and find the exact solution of the field equations. We use a dynamical systems approach to analyze global features of the phase space of solutions.

We investigate the scalar metric perturbations about a de Sitter brane universe in a 5-dimensional anti de Sitter bulk. We compare the master-variable formalism, describing metric perturbations in a 5-dimensional longitudinal gauge, with results in a Gaussian normal gauge. For a vacuum brane (with constant brane tension) there is a continuum of normalizable Kaluza-Klein modes, with m>3H/2, which remain in the vacuum state. A light radion mode, with m=\sqrt{2}H, satisfies the boundary conditions for two branes but is not normalizable in the single-brane case. When matter is introduced (as a test field) on the brane, this mode, together with the zero-mode and an infinite ladder of discrete tachyonic modes, become normalizable. However, the boundary condition requires the self-consistent 4-dimensional evolution of scalar field perturbations on the brane and the dangerous growing modes are not excited. These normalizable discrete modes introduce corrections at first-order to the scalar field perturbations computed in a slow-roll expansion. On super-Hubble scales, the correction is smaller than slow-roll corrections to the de Sitter background. However on small scales the corrections can become significant.

High-energy modifications to general relativity introduce changes to the perturbations generated during inflation, and the latest high-precision cosmological data can be used to place constraints on such modified inflation models. Recently it was shown that Randall-Sundrum type braneworld inflation leads to tighter constraints on quadratic and quartic potentials than in general relativity. We investigate how this changes with a Gauss-Bonnet correction term, which can be motivated by string theory. Randall-Sundrum models preserve the standard consistency relation between the tensor spectral index and the tensor-to-scalar ratio. The Gauss-Bonnet term breaks this relation, and also modifies the dynamics and perturbation amplitudes at high energies. We find that the Gauss-Bonnet term tends to soften the Randall-Sundrum constraints. The observational compatibility of the quadratic potential is strongly improved. For a broad range of energy scales, the quartic potential is rescued from marginal rejection. Steep inflation driven by an exponential potential is excluded in the Randall-Sundrum case, but the Gauss-Bonnet term leads to marginal compatibility for sufficient e-folds.

We calculate the amplitude of gravitational waves produced by inflation on a de Sitter brane embedded in five-dimensional anti-de Sitter bulk spacetime, extending previous calculations in Randall-Sundrum type cosmology to include the effect of induced gravity corrections on the brane. These corrections arise via a term in the brane action that is proportional to the brane Ricci scalar. We find that, as in the Randall-Sundrum case, there is a mass gap between the discrete zero-mode and a continuum of massive bulk modes, which are too heavy to be excited during inflation. We give the normalization of the zero-mode as a function of the Hubble rate on the brane and are thus able to calculate the high energy correction to the spectrum of gravitational wave (tensor) modes excited on large scales during inflation from initial vacuum fluctuations on small scales. We also calculate the amplitude of density (scalar) perturbations expected due to inflaton fluctuations on the brane, and show that the usual four-dimensional consistency relation for the tensor/scalar ratio remains valid for brane inflation with induced gravity corrections.

The observable universe could be a 1+3-surface (the "brane") embedded in a 1+3+d-dimensional spacetime (the "bulk"), with standard-model particles and fields trapped on the brane while gravity is free to access the bulk. At least one of the d extra spatial dimensions could be very large relative to the Planck scale, which lowers the fundamental gravity scale, possibly even down to the electroweak (~ TeV) level. This revolutionary picture arises in the framework of recent developments in M theory. The 1+10-dimensional M theory encompasses the known 1+9-dimensional superstring theories, and is widely considered to be a promising potential route to quantum gravity. General relativity cannot describe gravity at high enough energies and must be replaced by a quantum gravity theory, picking up significant corrections as the fundamental energy scale is approached. At low energies, gravity is localized at the brane and general relativity is recovered, but at high energies gravity "leaks" into the bulk, behaving in a truly 1+3+d-dimensional way. This introduces significant changes to gravitational dynamics and perturbations, with interesting testable implications for high-energy astrophysics, black holes and cosmology. Brane-world models offer a phenomenological way to test some of the novel predictions and corrections to general relativity that are implied by M theory. This review discusses the geometry, dynamics and perturbations of simple brane-world models for cosmology and astrophysics, mainly focusing on warped 5-dimensional brane-worlds based on the Randall-Sundrum models.

Brane-world models provide a phenomenology that allows us to explore the cosmological consequences of some M theory ideas, and at the same time to use precision cosmology as a test of these ideas. In order to achieve this, we need to understand how brane-world gravity affects cosmological dynamics and perturbations. This introductory review describes the key features of cosmological perturbations in a brane-world universe.

We study the cosmology of the Randall-Sundrum brane-world where the Einstein-Hilbert action is modified by curvature correction terms: a four-dimensional scalar curvature from induced gravity on the brane, and a five-dimensional Gauss-Bonnet curvature term. The combined effect of these curvature corrections to the action removes the infinite-density big bang singularity, although the curvature can still diverge for some parameter values. A radiation brane undergoes accelerated expansion near the minimal scale factor, for a range of parameters. This acceleration is driven by the geometric effects, without an inflaton field or negative pressures. At late times, conventional cosmology is recovered.

Observations indicate that the universe is effectively flat, but they do not rule out a closed universe. The role of positive curvature is negligible at late times, but can be crucial in the early universe. In particular, positive curvature allows for cosmologies that originate as Einstein static universes, and then inflate and later reheat to a hot big bang era. These cosmologies have no singularity, no "beginning of time", and no horizon problem. If the initial radius is chosen to be above the Planck scale, then they also have no quantum gravity era, and are described by classical general relativity throughout their history.

Non-commutative inflation is a modification of standard general relativity inflation which takes into account some effects of the space-time uncertainty principle motivated by ideas from string theory. The corrections to the primordial power spectrum which arise in a model of power-law inflation lead to a suppression of power on large scales, and produce a spectral index that is blue on large scales and red on small scales. This suppression and running of the spectral index are not imposed ad hoc, but arise from an early-Universe stringy phenomenology. We show that it can account for some loss of power on the largest scales that may be indicated by recent WMAP data. Cosmic microwave background anisotropies carry a signature of these very early Universe corrections, and can be used to place constraints on the parameters appearing in the non-commutative model. Applying a likelihood analysis to the WMAP data, we find the best-fit value for the critical wavenumber k_* (which involves the string scale) and for the exponent p (which determines the power-law inflationary expansion). The best-fit value corresponds to a string length of L_s \sim 10^{-28} cm.

We investigate the evolution of gravitational wave perturbations about a brane cosmology embedded in a five-dimensional anti-de Sitter bulk. During slow-roll inflation in a Randall-Sundrum brane-world, the zero mode of the 5-dimensional graviton is generated, while the massive modes remain in their vacuum state. When the zero mode re-enters the Hubble radius during radiation domination, massive modes are generated. We show that modes decouple in the low-energy/near-brane limit and develop perturbative techniques to calculate the mode-mixing at finite energy.

The standard cosmological model, based on general relativity with an inflationary era, is very effective in accounting for a broad range of observed features of the universe. However, the ongoing puzzles about the nature of dark matter and dark energy, together with the problem of a fundamental theoretical framework for inflation, indicate that cosmology may be probing the limits of validity of general relativity. The early universe provides a testing ground for theories of gravity, since gravitational dynamics can lead to characteristic imprints on the CMB and other cosmological observations. Precision cosmology is in principle a means to constrain and possibly falsify candidate quantum gravity theories like M theory. Generalized Randall-Sundrum brane-worlds provide a phenomenological means to test aspects of M theory. I outline the 1+3-covariant approach to cosmological perturbations in these brane-worlds, and its application to CMB anisotropies.

We investigate the effects on cosmological density perturbations of dark radiation in a Randall-Sundrum 2 type brane-world. Dark radiation in the background is limited by observational constraints to be a small fraction of the radiation energy density, but it has an interesting qualitative effect in the radiation era. On large scales, it serves to slightly suppress the radiation density perturbations at late times, while boosting the perturbations in dark radiation. In a kinetic (stiff) era, the suppression is much stronger, and drives the density perturbations to zero.

We show using covariant techniques that the Einstein static universe containing a perfect fluid is always neutrally stable against small inhomogeneous vector and tensor perturbations and neutrally stable against adiabatic scalar density inhomogeneities so long as c_{s}^2>1/5, and unstable otherwise. We also show that the stability is not significantly changed by the presence of a self-interacting scalar field source, but we find that spatially homogeneous Bianchi type IX modes destabilise an Einstein static universe. The implications of these results for the initial state of the universe and its pre-inflationary evolution are also discussed.

We use the cosmic microwave background temperature anisotropy to place limits on large-scale magnetic fields in an inhomogeneous (perturbed Friedmann) universe. If no assumptions are made about the spacetime geometry, only a weak limit can be deduced directly from the CMB. In the special case where spatial inhomogeneity is neglected to first order, the upper limit is much stronger, i.e. a few nano-G

Cosmic microwave background (CMB) observations provide in principle a high-precision test of models which are motivated by M theory. We set out the framework of a program to compute the tensor anisotropies in the CMB that are generated in braneworld models. In the simplest approximation, we show the braneworld imprint as a correction to the power spectra for standard temperature and polarization anisotropies.

We use the Zel'dovich approximation to analyse the amplification of magnetic fields in gravitational collapse of cold dark matter during the mildly nonlinear regime, and identify two key features. First, the anisotropy of collapse effectively eliminates one of the magnetic components, confining the field in the pancake plane. Second, in agreement with recent numerical simulations, we find that the shear anisotropy can amplify the magnetic field well beyond the isotropic case. Our results suggest that the magnetic strengths observed in spiral and disk galaxies today might have originated from seeds considerably weaker than previous estimates.

We investigate what are the key physical features that cause the development of a naked singularity, rather than a black hole, as the end-state of spherical gravitational collapse. We show that sufficiently strong shearing effects near the singularity delay the formation of the apparent horizon. This exposes the singularity to an external observer, in contrast to a black hole, which is hidden behind an event horizon due to the early formation of an apparent horizon.

We show that 5-dimensional Kaluza-Klein graviton stresses can slow the decay of shear anisotropy on the brane to observable levels, and we use cosmic microwave background anisotropies to place limits on the initial anisotropy induced by these stresses. An initial shear to Hubble distortion of only \sim 10^{-3}\Omega_0h_0^2 at the 5D Planck time would allow the observed large-angle CMB signal to be a relic mainly of KK tidal effects.

A static Friedmann brane in a 5-dimensional bulk (Randall-Sundrum type scenario) can have a very different relation between the density, pressure, curvature and cosmological constant than in the case of the general relativistic Einstein static universe. In particular, static Friedmann branes with zero cosmological constant and 3-curvature, but satisfying rho>0 and rho+3p>0, are possible. Furthermore, we find static Friedmann branes in a bulk that satisfies the Einstein equations but is not Schwarzschild-anti de Sitter or its specializations. In the models with negative bulk cosmological constant, a positive brane tension leads to negative density and 3-curvature.

We investigate how braneworld gravity affects gravitational collapse and black hole formation by studying Oppenheimer-Snyder-like collapse on a Randall-Sundrum type brane. Without making any assumptions about the bulk, we prove a no-go theorem: the exterior spacetime on the brane cannot be static, which is in stark contrast with general relativity. We also consider the role of Kaluza-Klein energy density in collapse, using a toy model.

We show that in a Randall-Sundrum II type braneworld, the vacuum exterior of a spherical star is not in general a Schwarzschild spacetime, but has radiative-type stresses induced by 5-dimensional graviton effects. Standard matching conditions do not lead to a unique exterior on the brane because of these 5-dimensional graviton effects. We find an exact uniform-density stellar solution on the brane, and show that the general relativity upper bound $GM/R<{4\over9}$ is reduced by 5-dimensional high-energy effects. The existence of neutron stars leads to a constraint on the brane tension that is stronger than the big bang nucleosynthesis constraint, but weaker than the Newton-law experimental constraint. We present two different non-Schwarzschild exteriors that match the uniform-density star on the brane, and we give a uniqueness conjecture for the full 5-dimensional problem.

We give a careful general relativistic and (1+3)-covariant analysis of cosmological peculiar velocities induced by matter density perturbations in the presence of a cosmological constant. In our quasi-Newtonian approach, constraint equations arise to maintain zero shear of the non-comoving fundamental worldlines which define a Newtonian-like frame, and these lead to the (1+3)-covariant dynamical equations, including a generalized Poisson-type equation. We investigate the relation between peculiar velocity and peculiar acceleration, finding the conditions under which they are aligned. In this case we find (1+3)-covariant relativistic generalizations of well-known Newtonian results.

We investigate the influence of cosmic magnetic fields on gravitational wave perturbations, and find exact solutions on large scales. We show that a large-scale magnetic field can generate large-scale non-decaying gravitational waves. In the general case where gravitational waves are generated by other mechanisms, a large-scale magnetic field introduces a new decaying tensor mode and modifies the non-decaying mode. The direct effect of the magnetic field is to damp the gravitational waves, while an indirect magneto-curvature effect can either damp or boost the waves. A magnetic field also leads to a breaking of statistical isotropy, and the magnetic imprint on the tensor spectrum in principle provides a means of detecting a primordial field.

Recent developments in string theory have led to 5-dimensional warped spacetime models in which standard-model fields are confined to a 3-brane (the observed universe), while gravity can propagate in the fifth dimension. Gravity is localized near the brane at low energies, even if the extra dimension is noncompact. A review is given of the classical geometry and dynamics of these brane-world models. The field equations on the brane modify the general relativity equations in two ways: local 5-D effects are imprinted on the brane as a result of its embedding, and are significant at high energies; nonlocal effects arise from the 5-D Weyl tensor. The Weyl tensor transmits tidal (Coulomb), gravitomagnetic and gravitational wave effects to the brane from the 5-D nonlocal gravitational field. Local high-energy effects modify the dynamics of inflation, and increase the amplitude of scalar and tensor perturbations generated by inflation. Nonlocal effects introduce new features in cosmological perturbations. They induce a non-adiabatic mode in scalar perturbations and massive modes in vector and tensor perturbations, and they can support vector perturbations even in the absence of matter vorticity. In astrophysics, local and nonlocal effects introduce fundamental changes to gravitational collapse and black hole solutions.

We examine the behavior of an anisotropic brane-world in the presence of inflationary scalar fields. We show that, contrary to naive expectations, a large anisotropy does not adversely affect inflation. On the contrary, a large initial anisotropy introduces more damping into the scalar field equation of motion, resulting in greater inflation. The rapid decay of anisotropy in the brane-world significantly increases the class of initial conditions from which the observed universe could have originated. This generalizes a similar result in general relativity. A unique feature of Bianchi I brane-world cosmology appears to be that for scalar fields with a large kinetic term the initial expansion of the universe is quasi-isotropic. The universe grows more anisotropic during an intermediate transient regime until anisotropy finally disappears during inflationary expansion.

In brane-world cosmologies of Randall-Sundrum type, we show that evolution of large-scale curvature perturbations may be determined on the brane, without solving the bulk perturbation equations. The influence of the bulk gravitational field on the brane is felt through a projected Weyl tensor which behaves effectively like an imperfect radiation fluid with anisotropic stress. We define curvature perturbations on uniform density surfaces for both the matter and Weyl fluids, and show that their evolution on large scales follows directly from the energy conservation equations for each fluid. The total curvature perturbation is not necessarily constant for adiabatic matter perturbations, but can change due to the Weyl entropy perturbation. To relate this curvature perturbation to the longitudinal gauge metric potentials requires knowledge of the Weyl anisotropic stress which is not determined by the equations on the brane. We discuss the implications for large-angle anisotropies on the cosmic microwave background sky.

In Randall-Sundrum-type brane-world cosmologies, density perturbations generate Weyl curvature in the bulk, which in turn backreacts on the brane via stress-energy perturbations. On large scales, the perturbation equations contain a closed system on the brane, which may be solved without solving for the bulk perturbations. Bulk effects produce a non-adiabatic mode, even when the matter perturbations are adiabatic, and alter the background dynamics. As a consequence, the standard evolution of large-scale fluctuations in general relativity is modified. The metric perturbation on large-scales is not constant during high-energy inflation. It is constant during the radiation era, except at most during the very beginning, if the energy is high enough.

We study adiabatic (curvature) and entropy (isocurvature) perturbations produced during a period of cosmological inflation that is driven by multiple scalar fields with an arbitrary interaction potential. A local rotation in field space is performed to separate out the adiabatic and entropy modes. The resulting field equations show explicitly how on large scales entropy perturbations can source adiabatic perturbations if the background solution follows a curved trajectory in field space, and how adiabatic perturbations cannot source entropy perturbations in the long-wavelength limit. It is the effective mass of the entropy field that determines the amplitude of entropy perturbations during inflation. We present two applications of the equations. First, we show why one in general expects the adiabatic and entropy perturbations to be correlated at the end of inflation, and calculate the cross-correlation in the context of a double inflation model with two non-interacting fields. Second, we consider two-field preheating after inflation, examining conditions under which entropy perturbations can alter the large-scale curvature perturbation and showing how our new formalism has advantages in numerical stability when the background solution follows a non-trivial trajectory in field space.

Dissipative effects during neutrino decoupling in the early universe create a small backreaction on the Hubble rate, and lead to a small rise in temperature and entropy. We use a simplified thermo-hydrodynamic model, which provides a causal approximation to kinetic theory, in order to estimate the backreaction effects and the entropy production.

We discuss the evolution of gravitational waves in a brane-world cosmology embedded in five-dimensional anti-de Sitter spacetime. We show that during slow-roll inflation, modelled as a period of quasi-de Sitter expansion on the brane, there is a discrete normalizable massless graviton mode. There is a mass gap due to the expansion, above which there is a continuum of massive modes. Only the massless mode is `light' compared with the Hubble scale during inflation, leading to the production of classical perturbations on large scales from vacuum fluctuations on small scales. We calculate the amplitude of these fluctuations at horizon-crossing and show that the standard four-dimensional result is recovered at low energies, but the amplitude of the perturbations is enhanced at high energies.

In Randall-Sundrum-type brane-world cosmologies, the dynamical equations on the three-brane differ from the general relativity equations by terms that carry the effects of imbedding and of the free gravitational field in the five-dimensional bulk. We derive the covariant nonlinear dynamical equations for the gravitational and matter fields on the brane, and then linearize to find the perturbation equations on the brane. The local energy-momentum corrections are significant only at very high energies. The imprint on the brane of the nonlocal gravitational field in the bulk is more subtle, and we provide a careful decomposition of this effect into nonlocal energy density, flux and anisotropic stress. The nonlocal energy density determines the tidal acceleration in the off-brane direction, and can oppose singularity formation via the generalized Raychaudhuri equation. Unlike the nonlocal energy density and flux, the nonlocal anisotropic stress is not determined by an evolution equation on the brane. In particular, isotropy of the cosmic microwave background may no longer guarantee a Friedmann geometry. Adiabatic density perturbations are coupled to perturbations in the nonlocal bulk field, and in general the system is not closed on the brane. But on super-Hubble scales, density perturbations can be evaluated by brane observers. Tensor perturbations on the brane are suppressed by local bulk effects during inflation, while the nonlocal effects can serve as a source or a sink. Vorticity on the brane decays as in general relativity, but nonlocal bulk effects can source the gravito-magnetic field, so that vector perturbations can be generated in the absence of vorticity.

We consider exact solutions for static black holes localized on a three-brane in five-dimensional gravity in the Randall-Sundrum scenario. We show that the Reissner-Nordstrom metric is an exact solution of the effective Einstein equations on the brane, re-interpreted as a black hole without electric charge, but with instead a tidal 'charge' arising via gravitational effects from the fifth dimension. The tidal correction to the Schwarzschild potential is negative, which is impossible in general relativity, and in this case only one horizon is admitted, located outside the Schwarzschild horizon. The solution satisfies a closed system of equations on the brane, and describes the strong-gravity regime. Current observations do not strongly constrain the tidal charge, and significant tidal corrections could in principle arise in the strong-gravity regime and for primordial black holes.

Magnetic fields are observed not only in stars, but in galaxies, clusters, and even high redshift Lyman-alpha systems. In principle, these fields could play an important role in structure formation and also affect the anisotropies in the cosmic microwave background radiation (CMB). The study of cosmological magnetic fields aims not only to quantify these effects on large-scale structure and the CMB, but also to answer one of the outstanding puzzles of modern cosmology: when and how do magnetic fields originate? They are either primordial, i.e. created before the onset of structure formation, or they are generated during the process of structure formation itself.

We consider slow-roll inflation in the context of recently proposed four-dimensional effective gravity induced on the world-volume of a three-brane in five-dimensional Einstein gravity. We find significant modifications of the simplest chaotic inflationary scenario when the five-dimensional Planck scale is below about 10^{17} GeV. We use the comoving curvature perturbation, which remains constant on super-Hubble scales, in order to calculate the spectrum of adiabatic density perturbations generated. Modifications to the Friedmann constraint equation lead to a faster Hubble expansion at high energies and a more strongly damped evolution of the scalar field. This assists slow-roll, enhances the amount of inflation obtained in any given model, and drives the perturbations towards an exactly scale-invariant Harrison-Zel'dovich spectrum. In chaotic inflation driven by a massive scalar field we show that inflation can occur at field values far below the four-dimensional Planck scale, though above the five-dimensional fundamental scale.

Motivated by the isotropy of the CMB spectrum, all existing studies of magnetised cosmological perturbations employ FRW backgrounds. However, it is important, to know the limits of this approximation and the effects one loses by neglecting the anisotropy of the background magnetic field. We develop a new treatment, which fully incorporates the anisotropic magnetic effects by allowing for a Bianchi I background universe. The anisotropy of the unperturbed model facilitates the closer study of the coupling between magnetism and geometry. The latter leads to a curvature stress, which accelerates positively curved perturbed regions and balances the effect of magnetic pressure gradients on matter condensations. We argue that the tension carried along the magnetic force-lines is the reason behind these magneto-curvature effects. For a relatively weak field, we also compare to the results of the almost-FRW approach. We find that some of the effects identified by the FRW treatment are in fact direction dependent, where the key direction is that of the background magnetic field vector. Nevertheless, the FRW-based approach to magnetised cosmological perturbations remains an accurate approximation, particularly on large scales, when one looks at the lowest order magnetic impact on gravitational collapse. On small scales however, the accuracy of the perturbed Friedmann framework may be compromised by extra shear effects.

The origin of cosmic gamma-ray bursts remains one of the most intriguing puzzles in astronomy. We suggest that purely general relativistic effects in the collapse of massive stars could account for these bursts. The late formation of closed trapped surfaces can occur naturally, allowing the escape of huge energy from curvature-generated fireballs, before these are hidden within a black hole.

A large-scale cosmic magnetic field affects not only the growth of density perturbations, but also rotational instabilities and anisotropic deformation in the density distribution. We give a fully relativistic treatment of all these effects, incorporating the magneto-curvature coupling that arises in a relativistic approach. We show that this coupling produces a small enhancement of the growing mode on superhorizon scales. The magnetic field generates new nonadiabatic constant and decaying modes, as well as nonadiabatic corrections to the standard growing and decaying modes. Magnetized isocurvature perturbations are purely decaying on superhorizon scales. On subhorizon scales before recombination, magnetized density perturbations propagate as magneto-sonic waves, leading to a small decrease in the spacing of acoustic peaks. Fluctuations in the field direction induce scale-dependent vorticity, and generate precession in the rotational vector. On small scales, magnetized density vortices propagate as Alfv\'{e}n waves during the radiation era. After recombination, they decay slower than non-magnetized vortices. Magnetic fluctuations are also an active source of anisotropic distortion in the density distribution. We derive the evolution equations for this distortion, and find a particular growing solution.

We study the evolution of gravitational waves through the preheating era that follows inflation. The oscillating inflaton drives parametric resonant growth of scalar field fluctuations, and although super-Hubble tensor modes are not strongly amplified, they do carry an imprint of preheating. This is clearly seen in the Weyl tensor, which provides a covariant description of gravitational waves.

The relative growth of field and metric perturbations during preheating is sensitive to initial conditions set in the preceding inflationary phase. Recent work suggests this may protect super-Hubble metric perturbations from resonant amplification during preheating. We show that this possibility is fragile and sensitive to the specific form of the interactions between the inflaton and other fields. The suppression is naturally absent in two classes of preheating in which either (1) the vacua of the non-inflaton fields during inflation are deformed away from the origin, or (2) the effective masses of non-inflaton fields during inflation are small but during preheating are large. Unlike the simple toy model of a $g^2 \phi^2 \chi^2$ coupling, most realistic particle physics models contain these other features. Moreover, they generically lead to both adiabatic and isocurvature modes and non-Gaussian scars on super-Hubble scales. Large-scale coherent magnetic fields may also appear naturally.

The spin rate \Omega of neutron stars at a given temperature T is constrained by the interplay between gravitational-radiation instabilities and viscous damping. Navier-Stokes theory has been used to calculate the viscous damping timescales and produce a stability curve for r-modes in the (\Omega,T) plane. In Navier-Stokes theory, viscosity is independent of vorticity, but kinetic theory predicts a coupling of vorticity to the shear viscosity. We calculate this coupling and show that it can in principle significantly modify the stability diagram at lower temperatures. As a result, colder stars can remain stable at higher spin rates.

Recently it has become clear that the resonant amplification of quantum field fluctuations at preheating must be accompanied by resonant amplification of scalar metric perturbations, since the two are united by Einstein's equations. Furthermore, this "metric preheating" enhances particle production and leads to gravitational rescattering effects even at linear order. In multi-field models with strong preheating (q \gg 1), metric perturbations are driven nonlinear, with the strongest amplification typically on super-Hubble scales (k \to 0). This amplification is causal, being due to the super- Hubble coherence of the inflaton condensate, and is accompanied by resonant growth of entropy perturbations. The amplification invalidates the use of the linearized Einstein field equations, irrespective of the amount of fine-tuning of the initial conditions. This has serious implications at all scales - from the large-angle cosmic microwave background (CMB) anisotropies to primordial black holes. We investigate the (q,k) parameter space in a two-field model, and introduce the time to nonlinearity, t_{nl}, as the timescale for the breakdown of the linearized Einstein equations. Backreaction effects are expected to shut down the linear resonances, but cannot remove the existing amplification, which threatens the viability of strong preheating when confronted with the CMB. We discuss ways to escape the above conclusions, including secondary phases of inflation and preheating solely to fermions. Finally we rank known classes of inflation from strongest (chaotic and strongly coupled hybrid inflation) to weakest (hidden sector, warm inflation) in terms of the distortion of the primordial spectrum due to these resonances in preheating.

Major advances in the observation and theory of cosmic microwave background anisotropies have opened up a new era in cosmology. This has encouraged the hope that the fundamental parameters of cosmology will be determined to high accuracy in the near future. However, this optimism should not obscure the ongoing need for theoretical developments that go beyond the highly successful but simplified standard model. Such developments include improvements in observational modelling (e.g. foregrounds, non-Gaussian features), extensions and alternatives to the simplest inflationary paradigm (e.g. non-adiabatic effects, defects), and investigation of nonlinear effects. In addition to well known nonlinear effects such as the Rees-Sciama and Ostriker-Vishniac effects, further nonlinear effects have recently been identified. These include a Rees-Sciama-type tensor effect, time-delay effects of scalar and tensor lensing, nonlinear Thomson scattering effects and a nonlinear shear effect. Some of the nonlinear effects and their potential implications are discussed.

Residual velocity dispersion in cold dark matter induces stresses which lead to effects that are absent in the idealized dust model. A previous Newtonian analysis showed how this approach can provide a theoretical foundation for the phenomenological adhesion model. We develop a relativistic kinetic theory generalization which also incorporates the anisotropic velocity dispersion that will typically be present. In addition to density perturbations, we consider the rotational and shape distortion properties of clustering. These quantities together characterize the linear development of density inhomogeneity, and we find exact solutions for their evolution. As expected, the corrections are small and arise only in the decaying modes, but their effect is interesting. One of the modes for density perturbations decays less rapidly than the standard decaying mode. The new rotational mode generates precession of the axis of rotation. The new shape modes produce additional distortion that remains frozen in during the subsequent (linear) evolution, despite the rapid decay of the terms that caused it.

General relativistic effects in the form of metric perturbations are usually neglected in the preheating era that follows inflation. We argue that in realistic multi-field models these effects are in fact crucial, and the fully coupled system of metric and quantum field fluctuations needs to be considered. Metric perturbations are resonantly amplified, breaking the scale-invariance of the primordial spectrum, and in turn stimulate scalar field resonances via gravitational rescattering. This non-gravitationally dominated nonlinear growth of gravitational fluctuations may have significant effects on the Doppler peaks in the cosmic background radiation, primordial black hole formation, gravitational waves and nonthermal symmetry restoration.

In a cold matter universe, the linearized gravito-magnetic tensor field satisfies a transverse condition (vanishing divergence) when it is purely radiative. We show that in the nonlinear theory, it is no longer possible to maintain the transverse condition, since it leads to a non-terminating chain of integrability conditions. These conditions are highly restrictive, and are likely to hold only in models with special symmetries, such as the known Bianchi and $G_2$ examples. In models with realistic inhomogeneity, the gravito-magnetic field is necessarily non-transverse at second and higher order.

We develop a new approach to local nonlinear effects in cosmic microwave background anisotropies, and discuss the qualitative features of these effects. New couplings of the baryonic velocity to radiation multipoles are found, arising from nonlinear Thomson scattering effects. We also find a new nonlinear shear effect on small angular scales. The full set of evolution and constraint equations is derived, including the nonlinear generalizations of the radiation multipole hierarchy, and of the dynamics of multi-fluids. These equations govern radiation anisotropies in any inhomogeneous spacetime, but their main application is to second-order effects in a universe that is close to the Friedmann models. Qualitative analysis is given here, and quantitative calculations are taken up in further papers.

Anisotropic stress contributions to the gravitational field can arise from magnetic fields, collisionless relativistic particles, hydrodynamic shear viscosity, gravitational waves, skew axion fields in low-energy string cosmologies, or topological defects. We investigate the effects of such stresses on cosmological evolution, and in particular on the dissipation of shear anisotropy. We generalize some previous results that were given for homogeneous anisotropic universes, by including small inhomogeneity in the universe. This generalization is facilitated by a covariant approach. We find that anisotropic stress dominates the evolution of shear, slowing its decay. The effect is strongest in radiation-dominated universes, where there is slow logarithmic decay of shear.

Recently a covariant approach to cold matter universes in the zero-shear hypersurfaces (or longitudinal) gauge has been developed. This approach reveals the existence of an integrability condition, which does not appear in standard non-covariant treatments. A simple derivation and generalization of the integrability condition is given, based on showing that the quasi-Newtonian models are a sub-class of the linearized `silent' models. The solution of the integrability condition implies a propagation equation for the acceleration. It is shown how the velocity and density perturbations are then obtained via this propagation equation. The density perturbations acquire a small relative-velocity correction on all scales, arising from the fully covariant general relativistic analysis.

Using a simple thermo-hydrodynamic model that respects relativistic causality, we revisit the analysis of qualitative features of acoustic oscillations in the photon-baryon fluid. The growing photon mean free path introduces transient effects that can be modelled by the causal generalization of relativistic Navier-Stokes-Fourier theory. Causal thermodynamics provides a more satisfactory hydrodynamic approximation to kinetic theory than the quasi-stationary (and non-causal) approximations arising from standard thermodynamics or from expanding the photon distribution to first order in the Thomson scattering time. The causal approach introduces small corrections to the dispersion relation obtained in quasi-stationary treatments. A dissipative contribution to the speed of sound slightly increases the frequency of the oscillations. The diffusion damping scale is slightly increased by the causal corrections. Thus quasi-stationary approximations tend to over-estimate the spacing and under-estimate the damping of acoustic peaks. In our simple model, the fractional corrections at decoupling are $\gtrsim 10^{-3}$.

We develop and apply a fully covariant 1+3 electromagnetic analogy for gravity. The free gravitational field is covariantly characterized by the Weyl gravito-electric and gravito-magnetic spatial tensor fields, whose dynamical equations are the Bianchi identities. Using a covariant generalization of spatial vector algebra and calculus to spatial tensor fields, we exhibit the covariant analogy between the tensor Bianchi equations and the vector Maxwell equations. We identify gravitational source terms, couplings and potentials with and without electromagnetic analogues. The nonlinear vacuum Bianchi equations are shown to be invariant under covariant spatial duality rotation of the gravito-electric and gravito-magnetic tensor fields. We construct the super-energy density and super-Poynting vector of the gravitational field as natural U(1) group invariants, and derive their super-energy conservation equation. A covariant approach to gravito-electric/magnetic monopoles is also presented.

We present a simple and thermodynamically consistent cosmology with a phenomenological model of quantum creation of radiation due to vacuum decay. Thermodynamics and Einstein's equations lead to an equation in which $H$ is determined by the particle number $N$. The model is completed by specifying the particle creation rate $\Gamma=\dot{N}/N$, which leads to a second-order evolution equation for $H$. We propose a simple $\Gamma$ that is naturally defined and that conforms to the thermodynamical conditions: (a) the entropy production rate starts at a maximum; (b) the initial vacuum (for radiation) is a non-singular regular vacuum; and (c) the creation rate is initially higher than the expansion rate $H$, but then falls below $H$. The evolution equation for $H$ then has a remarkably simple exact solution, in which a non-adiabatic inflationary era exits smoothly to the radiation era, without a reheating transition. For this solution, we give exact expressions for the cosmic scale factor, energy density of radiation and vacuum, temperature, entropy and super-horizon scalar perturbations.

We find a simple exact model of radiating stellar collapse, with a shear-free and non-accelerating interior matched to a Vaidya exterior. The heat flux is subject to causal thermodynamics, leading to self-consistent determination of the temperature $T$. We solve for $T$ exactly when the mean collision time $\tau_{c}$ is constant, and perturbatively in a more realistic case of variable $\tau_{c}$. Causal thermodynamics predicts temperature behaviour that can differ significantly from the predictions of non-causal theory. In particular, the causal theory gives a higher central temperature and greater temperature gradient.

We investigate cosmological density perturbations in a covariant and gauge-invariant formalism, incorporating relativistic causal thermodynamics to give a self-consistent description. The gradient of density inhomogeneities splits covariantly into a scalar part, equivalent to the usual density perturbations, a rotational vector part that is determined by the vorticity, and a tensor part that describes the shape. We give the evolution equations for these parts in the general dissipative case. Causal thermodynamics gives evolution equations for viscous stress and heat flux, which are coupled to the density perturbation equation and to the entropy and temperature perturbation equations. We give the full coupled system in the general dissipative case, and simplify the system in certain cases. A companion paper uses the general formalism to analyze damping of density perturbations before last scattering.

We revisit the issue of integrability conditions for the irrotational silent cosmological models. We formulate the problem both in 1+3 covariant and 1+3 orthonormal frame notation, and show there exists a series of constraint equations that need to be satisfied. These conditions hold identically for FLRW-linearised silent models, but not in the general exact non-linear case. Thus there is a linearisation instability, and it is highly unlikely that there is a large class of silent models. We conjecture that there are no spatially inhomogeneous solutions with Weyl curvature of Petrov type I, and indicate further issues that await clarification.

We develop a nonlinear generalisation of the causal linear thermodynamics of bulk viscosity, incorporating positivity of the entropy production rate and the effective specific entropy. The theory is applied to viscous fluid inflation (which is necessarily far from equilibrium), and we find thermodynamically consistent inflationary solutions, both exponential and power-law.

The reheating process in inflationary universe models is considered as an out-of-equilibrium mixture of two interacting and reacting fluids, and studied within the framework of causal, irreversible thermodynamics. The evolution of the temperature and the decay rate as determined by causal thermodynamics are estimated at different stages of the process. A simple model is also used to find the perturbations of the expansion rate, including the possibility of damped oscillations.

Gravity waves in irrotational dust spacetimes are characterised by nonzero magnetic Weyl tensor $H_{ab}$. In the linearised theory, the divergence of $H_{ab}$ is set to zero. Recently Lesame et al. [Phys. Rev. D {\bf 53}, 738 (1996)] presented an argument to show that, in the exact nonlinear theory, $div H=0$ forces $H_{ab}=0$, thus implying a linearisation instability for gravity waves interacting with matter. However a sign error in the equations invalidates their conclusion. Bianchi type V spacetimes are shown to include examples with $div H=0\neq H_{ab}$. An improved covariant formalism is used to show that in a generic irrotational dust spacetime, the covariant constraint equations are preserved under evolution. It is shown elsewhere that $\mbox{div} H=0$ does not generate further conditions.

I review the causal relativistic thermodynamics developed by Israel and Stewart, and discuss some applications in cosmology and astrophysics. The lectures begin with an overview of relativistic fluid dynamics (in a covariant formalism) and equilibrium thermodynamics. Causal bulk viscosity in cosmology is considered in detail, including some new results.

We calculate the off--center observational relations in a spherically symmetric dust universe that is inhomogeneous at small redshifts. In contrast to the usual model, in which the CMBR dipole is interpreted as a Doppler effect due to peculiar velocity, our model explores an alternative interpretation, in which the CMBR dipole is non--Doppler, and the observer is comoving with the mean matter flow. We do not assume a background frame relative to which peculiar velocities are calculated. Our analysis is fully nonlinear and the density contrast is not assumed to be small. We obtain exact expressions for the Hubble and deceleration parameters, and find that both parameters have quadrupole anisotropies, but no dipoles. A simple numerical procedure for calculating the CMBR dipole anisotropy in our model is presented, and the observed 0.1% dipole is shown to be reproducible with a reasonable choices of parameters.

A recent paper (Martinez--Gonzalez & Sanz 1995) showed that if the universe is homogeneous but anisotropic, then the small quadrupole anisotropy in the cosmic microwave background radiation implies that the spacetime anisotropy is very small. We point out that more general results may be established, without assuming a priori homogeneity. We have proved that small anisotropies in the microwave background imply that the universe is almost Friedmann--Robertson--Walker. Furthermore, the quadrupole and octopole place direct and explicit limits on the degree of anisotropy and inhomogeneity, as measured by the shear, vorticity, Weyl tensor and density gradients. In the presence of inhomogeneity, it is only possible to set a much weaker limit on the shear than that given by Martinez--Gonzalez & Sanz.

We consider directly the equations by which matter imposes anisotropies on freely propagating background radiation, leading to a new way of using anisotropy measurements to limit the deviations of the Universe from a Friedmann-Robertson-Walker (FRW) geometry. This approach is complementary to the usual Sachs-Wolfe approach: the limits obtained are not as detailed, but they are more model-independent. We also give new results about combined matter-radiation perturbations in an almost-FRW universe, and a new exact solution of the linearised equations.