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 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 weak gravitational lensing formalism can be extended to the strong lensing regime by integrating a nonlinear version of the geodesic deviation equation. The resulting "roulette" expansion generalises the notion of convergence, shear and flexion to arbitrary order. The independent coefficients of this expansion are screen space gradients of the optical tidal tensor which approximates to the usual lensing potential in the weak field limit. From lensed images, knowledge of the roulette coefficients can in principle be inverted to reconstruct the mass distribution of a lens. In this paper, we simplify the roulette expansion and derive a family of recursion relations between the various coefficients, generalising the Kaiser-Squires relations beyond the weak-lensing regime.

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.

We investigate for the first time weak lensing in a high-resolution relativistic N-body simulation, which includes all relevant effects from general relativity. We integrate the photon geodesics backwards from the observer to the emitters. Our analysis is fully relativistic and non-perturbative for the scalar part of the gravitational potential and first-order in the vector part, frame dragging. We solve the Sachs optical equations and study in detail the weak-lensing convergence, ellipticity and rotation. We present the angular power spectra and one-point probability distribution functions for the weak-lensing variables, which we find are broadly in agreement with comparable Newtonian simulations. Our geometric approach, however, is more robust and flexible, and can therefore be applied consistently to non-standard cosmologies and modified theories of gravity.

The main science driver for the coming generation of cosmological surveys is understanding dark energy which relies on testing General Relativity on the largest scales. Once we move beyond the simplest explanation for dark energy of a cosmological constant, the space of possible theories becomes both vast and extremely hard to compute realistic observables. A key discriminator of a cosmological constant, however, is that the growth of structure is scale-invariant on large scales. By carefully weighting observables derived from distributions of galaxies and a dipole pattern in their apparent sizes, we construct a null test which vanishes for any model of gravity or dark energy where the growth of structure is scale-independent. It relies only on very few assumptions about cosmology, and does not require any modelling of the growth of structure. We show that with a survey like DESI a scale-dependence of the order of 10-20 percent can be detected at 3 sigma with the null test, which will drop by a factor of 2 for a survey like the Square Kilometre Array. We also show that the null test is very insensitive to typical uncertainties in other cosmological parameters including massive neutrinos and scale-dependent bias, making this a key null test for dark energy.

Two-parameter perturbation theory is a scheme tailor-made to consistently include nonlinear density contrasts on small scales ($<100\; \mathrm{Mpc}$), whilst retaining a traditional approach to cosmological perturbations in the long-wavelength universe. In this paper we study the solutions that arise from this theory in a spatially-flat dust-filled cosmology, and what these imply for the bispectrum of matter. This is achieved by using Newtonian perturbation theory to model the gravitational fields of nonlinear structures in the quasi-linear regime, and then using the resulting solutions as source terms for the cosmological equations. We find that our approach results in the leading-order part of the cosmological gravitational potentials being identical to those that result from standard cosmological perturbation theory at second-order, while the dark matter bispectrum itself yields some differences on Hubble scales. This demonstrates that our approach is sufficient to capture most leading-order relativistic effects, but within a framework that is far easier to generalize. We expect this latter property to be particularly useful for calculating leading-order relativistic corrections to the matter power spectrum, as well as for deriving predictions for relativistic effects in alternative theories of gravity.

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

An important part of cosmological model fitting relies on correlating distance indicators of objects (for example type Ia supernovae) with their redshift, often illustrated on a Hubble diagram. Comparing the observed correlation with a homogeneous model is one of the key pieces of evidence for dark energy. The presence of cosmic structures introduces a bias and scatter, mainly due to gravitational lensing and peculiar velocities, but also due to smaller non-linear relativistic contributions which are more difficult to account for. For the first time we perform ray tracing onto halos in a relativistic N-body simulation. Our simulation is the largest that takes into account all leading-order corrections from general relativity in the evolution of structure, and we present a novel methodology for working out the non-linear projection of that structure onto the observer's past light cone. We show that the mean of the bias in the Hubble diagram is indeed as small as expected from perturbation theory. However, the distribution of sources is significantly skewed with a very long tail of highly magnified objects and we illustrate that the bias of cosmological parameters strongly depends on the function of distance which we consider.

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.

Weak lensing will play an important role in future cosmological surveys, including e.g. Euclid and SKA. Sufficiently accurate theoretical predictions are important for correctly interpreting these surveys and hence for extracting correct cosmological parameter estimations. We quantify for the first time in a relativistic setting how many post-Born and lens-lens coupling corrections are required for sub-percent accuracy of the theoretical weak lensing convergence for $z\le 2$ (the primary weak lensing range for Euclid and SKA). We do this by ray-tracing through a fully relativistic exact solution of the Einstein Field Equations which consists of randomly packed mass-compensated underdensities of realistic amplitudes. We find that including lens-lens coupling terms and post-Born corrections up to second and third order respectively is sufficient for sub-percent accuracy of the convergence along $94\%$ of the studied light rays. We also find that a significant percentage of the studied rays have post-Born corrections of size over $10\%$ of the usual gravitational convergence, $\kappa^{(1)}$, and several rays even have post-Born corrections several times the size of $\kappa^{(1)}$ at $z = 2$.

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.

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.

A persistent theme in the study of dark energy is the question of whether it really exists or not. It is often claimed hat we are mis-calculating the cosmological model by neglecting the effects associated with averaging over large-scale structures. In the Newtonian approximation this is clear: there is no effect. Within the full relativistic picture this remains an important open question however, owing to the complex mathematics involved. We study this issue using particle numerical simulations which account for all relevant relativistic effects without any problems from shell crossing. In this context we show for the first time that the backreaction from structure can differ by many orders of magnitude depending upon the slicing of spacetime one chooses to average over. In the worst case, where smoothing is carried out in synchronous spatial surfaces, the corrections can reach ten percent and more. However, when smoothing on the constant time hypersurface of the Newtonian gauge, backreaction contributions remain 3-5 orders of magnitude smaller.

Relativistic contributions to the dynamics of structure formation come in a variety of forms, and can potentially give corrections to the standard picture on typical scales of 100 Mpc. These corrections cannot be obtained by Newtonian numerical simulations, so it is important to accurately estimate the magnitude of these relativistic effects. Density fluctuations couple to produce a background of gravitational waves, which is larger than any primordial background. A similar interaction produces a much larger spectrum of vector modes which represent the frame-dragging rotation of spacetime. These can change the metric at the percent level in the concordance model at scales below the equality scale. Vector modes modify the lensing of background galaxies by large-scale structure. This gives in principle the exciting possibility of measuring relativistic frame dragging effects on cosmological scales. The effects of the non-linear tensor and vector modes on the cosmic convergence are computed and compared to first-order lensing contributions from density fluctuations, Doppler lensing, and smaller Sachs-Wolfe effects. The lensing from gravitational waves is negligible so we concentrate on the vector modes. We show the relative importance of this for future surveys such as Euclid and SKA. We find that these non-linear effects only marginally affect the overall weak lensing signal so they can safely be neglected in most analyses, though are still much larger than the linear Sachs-Wolfe terms. The second-order vector contribution can dominate the first-order Doppler lensing term at moderate redshifts and are actually more important for survey geometries like the SKA.

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.

Recent precision cosmological parameter constraints imply that the spatial curvature of the Universe is essentially dynamically negligible - but only if relatively strong assumptions are made about the equation of state of dark energy (DE). When these assumptions are relaxed, strong degeneracies arise that make it hard to disentangle DE and curvature, degrading the constraints. We show that forthcoming 21cm intensity mapping experiments such as HIRAX are ideally designed to carry out model-independent curvature measurements, as they can measure the clustering signal at high redshift with sufficient precision to break many of the degeneracies. We consider two different model-independent methods, based on `avoiding' the DE-dominated regime and non-parametric modelling of the DE equation of state respectively. Our forecasts show that HIRAX will be able to improve upon current model-independent constraints by around an order of magnitude, reaching percent-level accuracy even when an arbitrary DE equation of state is assumed. In the same model-independent analysis, the sample variance limit for a similar survey is another order of magnitude better.

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.

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.

We present a new extension of the weak lensing formalism capable of describing strongly lensed images. This paper accompanies Paper I, arXiv:1603.04698 where we provided a condensed overview of the approach and illustrated how it works. Here we give all the necessary details, together with some more explicit examples. We solve the non-linear geodesic deviation equation order-by-order, keeping the leading derivatives of the optical tidal matrix, giving rise to a series of maps from which a complete strongly lensed image is formed. The family of maps are decomposed by separating the trace and trace-free parts of each map. Each trace-free tensor represents an independent spin mode, which distort circles into a variety of roulettes in the screen-space. It is shown how summing this series expansion allows us to create large strongly lensed images in regions where convergence, shear and flexion are not sufficient. This paper is a detailed exposition of Paper I which presents the key elements of the subject matter in a wider context.

We present a new perspective on gravitational lensing. We describe a new extension of the weak lensing formalism capable of describing strongly lensed images. By integrating the non-linear geodesic deviation equation, the amplification matrix of weak lensing is generalised to a sum over independent amplification tensors of increasing rank. We show how an image distorted by a generic lens may be constructed as a sum over `roulettes', which are the natural curves associated with the independent spin modes of the amplification tensors. Highly distorted images can be constructed even for large sources observed near or within the Einstein radius of a lens where the shear and convergence are large. The amplitude of each roulette is formed from a sum over appropriate derivatives of the lensing potential. Consequently, measuring these individual roulettes for images around a lens gives a new way to reconstruct a strong lens mass distribution without requiring a lens model. This formalism generalises the convergence, shear and flexion of weak lensing to arbitrary order, and provides a unified bridge between the strong and weak lensing regimes. This overview paper is accompanied by a much more detailed Paper II, arXiv:1603.04652

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

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.

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.

One of the key issues in cosmology is to establish the nature of dark energy, and to determine whether the equation of state evolves with time. When estimating this from distance measurements there is a degeneracy with the matter density. We show that there exists a simple function of the dark energy equation of state and its first derivative which is independent of this degeneracy at all redshifts, and so is a much more robust determinant of the evolution of dark energy than just its derivative. We show that this function can be well determined at low redshift from supernovae using Gaussian Processes, and that this method is far superior to a variety of parameterisations which are also subject to priors on the matter density. This shows that parametrised models give very biased constraints on the evolution of dark energy.

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 consider the entropy associated with the large-scale structure of the Universe in the linear regime, where the Universe can be described by a perturbed Friedmann-Lema\^itre spacetime. In particular, we compare two different definitions proposed in the literature for the entropy using a spatial averaging prescription. For one definition, the entropy of the large-scale structure for a given comoving volume always grows with time, both for a CDM and a $\Lambda$CDM model. In particular, while it diverges for a CDM model, it saturates to a constant value in the presence of a cosmological constant. The use of a light-cone averaging prescription in the context of the evaluation of the entropy is also discussed.

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.

Weak gravitational lensing is normally assumed to have only two principle effects: a magnification of a source and a distortion of the sources shape in the form of a shear. However, further distortions are actually present owing to changes in the gravitational field across the scale of the ray bundle of light propagating to us, resulting in the familiar arcs in lensed images. This is normally called the flexion, and is approximated by Taylor expanding the shear and magnification across the image plane. However, the physical origin of this effect arises from higher-order corrections in the geodesic deviation equation governing the gravitational force between neighbouring geodesics - so involves derivatives of the Riemann tensor. We show that integrating the second-order geodesic deviation equation results in a 'Hessian map' for gravitational lensing, which is a higher-order addition to the Jacobi map. We derive the general form of the Hessian map in an arbitrary spacetime paying particular attention to the separate effects of local Ricci versus non-local Weyl curvature. We then specialise to the case of a perturbed FLRW model, and give the general form of the Hessian for the first time. This has a host of new contributions which could in principle be used as tests for modified gravity.

The large-scale homogeneity and isotropy of the universe is generally thought to imply a well defined background cosmological model. It may not. Smoothing over structure adds in an extra contribution, transferring power from small scales up to large. Second-order perturbation theory implies that the effect is small, but suggests that formally the perturbation series may not converge. The amplitude of the effect is actually determined by the ratio of the Hubble scales at matter-radiation equality and today - which are entirely unrelated. This implies that a universe with significantly lower temperature today could have significant backreaction from more power on small scales, and so provides the ideal testing ground for understanding backreaction. We investigate this using two different N-body numerical simulations - a 3D Newtonian and a 1D simulation which includes all relevant relativistic effects. We show that while perturbation theory predicts an increasing backreaction as more initial small-scale power is added, in fact the virialisation of structure saturates the backreaction effect at the same level independently of the equality scale. This implies that backreaction is a small effect independently of initial conditions. Nevertheless, it may still contribute at the percent level to certain cosmological observables and therefore it cannot be neglected in precision cosmology.

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.

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

We present results from a numerical code implementing a new method to solve the master equations describing the evolution of linear perturbations in a spherically symmetric but inhomogeneous background. This method can be used to simulate several configurations of physical interest, such as relativistic corrections to structure formation, the lensing of gravitational waves and the evolution of perturbations in a cosmological void model. This paper focuses on the latter problem, i.e. structure formation in a Hubble scale void in the linear regime. This is considerably more complicated than linear perturbations of a homogeneous and isotropic background because the inhomogeneous background leads to coupling between density perturbations and rotational modes of the spacetime geometry, as well as gravitational waves. Previous analyses of this problem ignored this coupling in the hope that the approximation does not affect the overall dynamics of structure formation in such models. We show that for a giga-parsec void, the evolution of the density contrast is well approximated by the previously studied decoupled evolution only for very large-scale modes. However, the evolution of the gravitational potentials within the void is inaccurate at more than the $10\%$ level, and is even worse on small scales.

Averaging in general relativity is a complicated operation, due to the general covariance of the theory and the non-linearity of Einstein's equations. The latter of these ensures that smoothing spacetime over cosmological scales does not yield the same result as solving Einstein's equations with a smooth matter distribution, and that the smooth models we fit to observations need not be simply related to the actual geometry of spacetime. One specific consequence of this is a decoupling of the geometrical spatial curvature term in the metric from the dynamical spatial curvature in the Friedmann equation. Here we investigate the consequences of this decoupling by fitting to a combination of HST, CMB, SNIa and BAO data sets. We find that only the geometrical spatial curvature is tightly constrained, and that our ability to constrain dark energy dynamics will be severely impaired until we gain a thorough understanding of the averaging problem in cosmology.

A new non-parametric method based on Gaussian Processes was proposed recently to measure the Hubble constant $H_0$. The freedom in this approach comes in the chosen covariance function, which determines how smooth the process is and how nearby points are correlated. We perform coverage tests with a thousand mock samples within the LCDM model in order to determine what covariance function provides the least biased results. The function Matern(5/2) is the best with sligthly higher errors than other covariance functions, although much more stable when compared to standard parametric analyses.

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.

We present new numerical cosmological solutions of the Einstein Field Equations. The spacetime is spherically symmetric with a source of dust and radiation approximated as a perfect fluid. The dust and radiation are necessarily non-comoving due to the inhomogeneity of the spacetime. Such a model can be used to investigate non-linear general relativistic effects present during decoupling or big-bang nucleosynthesis, as well as for investigating void models of dark energy with isocurvature degrees of freedom. We describe the full evolution of the spacetime as well as the redshift and luminosity distance for a central observer. After demonstrating accuracy of the code, we consider a few example models, and demonstrate the sensitivity of the late time model to the degree of inhomogeneity of the initial radiation contrast.

An intriguing discrepancy emerging in the concordance model of cosmology is the tension between the locally measured value of the Hubble rate, and the 'global' value inferred from the cosmic microwave background (CMB). This could be due to systematic uncertainties when measuring $H_0$ locally, or it could be that we live in a highly unlikely Hubble bubble, or other exotic scenarios. We point out that the global $H_0$ can be found by extrapolating $H(z)$ data points at high-$z$ down to $z=0$. By doing this in a Bayesian non-parametric way we can find a model-independent value for $H_0$. We apply this to 19 measurements based on differential age of passively evolving galaxies as cosmic chronometers. Using Gaussian processes, we find $H_0=64.9 \pm 4.2$ km s$^{-1}$ Mpc$^{-1}$ $(1\sigma)$, in agreement with the CMB value, but reinforcing the tension with the local value. An analysis of possible sources of systematic errors shows that the stellar population synthesis model adopted may change the results significantly, being the main concern for subsequent studies. Forecasts for future data show that distant $H(z)$ measurements can be a robust method to determine $H_0$, where a focus in precision and a careful assessment of systematic errors are required.

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.

Gaussian processes are a fully Bayesian smoothing technique that allows for the reconstruction of a function and its derivatives directly from observational data, without assuming a specific model or choosing a parameterization. This is ideal for constraining dark energy because physical models are generally phenomenological and poorly motivated. Model-independent constraints on dark energy are an especially important alternative to parameterized models, as the priors involved have an entirely different source so can be used to check constraints formulated from models or parameterizations. A critical prior for Gaussian process reconstruction lies in the choice of covariance function. We show how the choice of covariance function affects the result of the reconstruction, and present a choice which leads to reliable results for present day supernovae data. We also introduce a method to quantify deviations of a model from the Gaussian process reconstructions.

We present a new programme for placing constraints on radial inhomogeneity in a dark-energy dominated universe. We introduce a new measure to quantify violations of the Copernican principle. Any violation of this principle would interfere with our interpretation of any dark-energy evolution. In particular, we find that current observations place reasonably tight constraints on possible late-time violations of the Copernican principle: the allowed area in the parameter space of amplitude and scale of a spherical inhomogeneity around the observer has to be reduced by a factor of three so as to confirm the Copernican principle. Then, by marginalizing over possible radial inhomogeneity we provide the first constraints on the cosmological constant which are free of the homogeneity prior prevalent in cosmology.

Correctly interpreting observations of sources such as type Ia supernovae (SNe Ia) require knowledge of the power spectrum of matter on AU scales - which is very hard to model accurately. Because under-dense regions account for much of the volume of the universe, light from a typical source probes a mean density significantly below the cosmic mean. The relative sparsity of sources implies that there could be a significant bias when inferring distances of SNe Ia, and consequently a bias in cosmological parameter estimation. While the weak lensing approximation should in principle give the correct prediction for this, linear perturbation theory predicts an effectively infinite variance in the convergence for ultra-narrow beams. We attempt to quantify the effect typically under-dense lines of sight might have in parameter estimation by considering three alternative methods for estimating distances, in addition to the usual weak lensing approximation. We find in each case this not only increases the errors in the inferred density parameters, but also introduces a bias in the posterior value.

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.

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

The question of whether the Universe is spatially homogeneous and isotropic on the largest scales is of fundamental importance to cosmology, but has not yet been answered decisively. Surprisingly, neither an isotropic primary CMB nor combined observations of luminosity distances and galaxy number counts are sufficient to establish such a result. The inclusion of the Sunyaev-Zel'dovich effect in CMB observations, however, dramatically improves this situation. We show that even a solitary observer who sees an isotropic blackbody CMB can conclude that the universe is homogeneous and isotropic in their causal past when the Sunyaev-Zel'dovich effect is present. Critically, however, the CMB must either be viewed for an extended period of time, or CMB photons that have scattered more than once must be detected. This result provides a theoretical underpinning for testing the Cosmological Principle with observations of the CMB alone.

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.

An important issue in cosmology is reconstructing the effective dark energy equation of state directly from observations. With few physically motivated models, future dark energy studies cannot only be based on constraining a dark energy parameter space, as the errors found depend strongly on the parameterisation considered. We present a new non-parametric approach to reconstructing the history of the expansion rate and dark energy using Gaussian Processes, which is a fully Bayesian approach for smoothing data. We present a pedagogical introduction to Gaussian Processes, and discuss how it can be used to robustly differentiate data in a suitable way. Using this method we show that the Dark Energy Survey - Supernova Survey (DES) can accurately recover a slowly evolving equation of state to sigma_w = +-0.04 (95% CL) at z=0 and +-0.2 at z=0.7, with a minimum error of +-0.015 at the sweet-spot at z~0.14, provided the other parameters of the model are known. Errors on the expansion history are an order of magnitude smaller, yet make no assumptions about dark energy whatsoever. A code for calculating functions and their first three derivatives using Gaussian processes has been developed and is available for download at http://www.acgc.uct.ac.za/~seikel/GAPP/index.html .

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.

Assuming the universe is spatially homogeneous on the largest scales lays the foundation for almost all cosmology. This idea is based on the Copernican principle, that we are not at a particularly special place in the universe. Surprisingly, this philosophical assumption has yet to be rigorously demonstrated independently of the standard paradigm. This issue has been brought to light by cosmological models which can potentially explain apparent acceleration by spatial inhomogeneity rather than dark energy. These models replace the temporal fine tuning associated with Lambda with a spatial fine tuning, and so violate the Copernican assumption. While is seems unlikely that such models can really give a realistic solution to the dark energy problem, they do reveal how poorly constrained radial inhomogeneity actually is. So the bigger issue remains: How do we robustly test the Copernican principle independently of dark energy or theory of gravity?

Cosmological perturbation theory relies on the decomposition of perturbations into so-called scalar, vector and tensor modes. This decomposition is non-local and depends on unknowable boundary conditions. The non-locality is particularly important at second- and higher-order because perturbative modes are sourced by products of lower-oder modes, which must be integrated over all space in order to isolate each mode. However, given a trace-free rank-2 tensor, a locally defined scalar mode may be trivially derived by taking two divergences, which knocks out the vector and tensor degrees of freedom. A similar local differential operation will return a pure vector mode. This means that scalar and vector degrees of freedom have local descriptions. The corresponding local extraction of the tensor mode is unknown however. We give it here. The operators we define are useful for defining gauge-invariant quantities at second-order. We perform much of our analysis using an index-free `vector-calculus' approach which makes manipulating tensor equations considerably simpler.

Structure occurs over a vast range of scales in the universe. Our large-scale cosmological models are coarse-grained representations of what exists, which have much less structure than there really is. An important problem for cosmology is determining the influence the small-scale structure in the universe has on its large-scale dynamics and observations. Is there a significant, general relativistic, backreaction effect from averaging over structure? One issue is whether the process of smoothing over structure can contribute to an acceleration term and so alter the apparent value of the cosmological constant. If this is not the case, are there other aspects of concordance cosmology that are affected by backreaction effects? Despite much progress, this 'averaging problem' is still unanswered, but it cannot be ignored in an era of precision cosmology.

The relativistic theory of structure formation in cosmology is based mainly on linear perturbations about a homogeneous background. But we are now driven to understand the theory of higher-order perturbations in full detail, both from observational and theoretical points of view. An important aspect of this lies in defining gauge-invariant perturbations at any order. We present a new covariant approach to this based on a local separation of tensors into their scalar, vector and tensor parts. Such a local decomposition necessarily requires a trio of mutually annihilating differential operators which form the basis for defining gauge-invariant objects. It makes no use of non-local Green's functions or pre-defined gauges, and can be used to define families of scalar, vector and tensor modes at any order one chooses.

Smoothing over structures in general relativity leads to a renormalisation of the background, and potentially many other effects which are poorly understood. Observables such as the distance-redshift relation when averaged on the sky do not necessarily yield the same smooth model which arises when performing spatial averages. These issues are thought to be of technical interest only in the standard model of cosmology, giving only tiny corrections. However, when we try to calculate observable quantities such as the all-sky average of the distance-redshift relation, we find that perturbation theory delivers divergent answers in the UV and corrections to the background of order unity. There are further problems. Second-order perturbations are the same size as first-order, and fourth-order at least the same as second, and possibly much larger, owing to the divergences. Much hinges on a coincidental balance of 2 numbers: the primordial power, and the ratio between the comoving Hubble scales at matter-radiation equality and today. Consequently, it is far from obvious that backreaction is irrelevant even in the concordance model, however natural it intuitively seems.

The calculation of the averaged Hubble expansion rate in an averaged perturbed Friedmann-Lemaitre-Robertson-Walker cosmology leads to small corrections to the background value of the expansion rate, which could be important for measuring the Hubble constant from local observations. It also predicts an intrinsic variance associated with the finite scale of any measurement of H_0, the Hubble rate today. Both the mean Hubble rate and its variance depend on both the definition of the Hubble rate and the spatial surface on which the average is performed. We quantitatively study different definitions of the averaged Hubble rate encountered in the literature by consistently calculating the backreaction effect at second order in perturbation theory, and compare the results. We employ for the first time a recently developed gauge-invariant definition of an averaged scalar. We also discuss the variance of the Hubble rate for the different definitions.

The dimming of Type Ia supernovae could be the result of Hubble-scale inhomogeneity in the matter and spatial curvature, rather than signaling the presence of a dark energy component. A key challenge for such models is to fit the detailed spectrum of the cosmic microwave background (CMB). We present a detailed discussion of the small-scale CMB in an inhomogeneous universe, focusing on spherically symmetric `void' models. We allow for the dynamical effects of radiation while analyzing the problem, in contrast to other work which inadvertently fine tunes its spatial profile. This is a surprisingly important effect and we reach substantially different conclusions. Models which are open at CMB distances fit the CMB power spectrum without fine tuning; these models also fit the supernovae and local Hubble rate data which favours a high expansion rate. Asymptotically flat models may fit the CMB, but require some extra assumptions. We argue that a full treatment of the radiation in these models is necessary if we are to understand the correct constraints from the CMB, as well as other observations which rely on it, such as spectral distortions of the black body spectrum, the kinematic Sunyaev-Zeldovich effect or the Baryon Acoustic Oscillations.

Explaining the well established observation that the expansion rate of the universe is apparently accelerating is one of the defining scientific problems of our age. Within the standard model of cosmology, the repulsive 'dark energy' supposedly responsible has no explanation at a fundamental level, despite many varied attempts. A further important dilemma in the standard model is the Lithium problem, which is the substantial mismatch between the theoretical prediction for 7-Li from Big Bang Nucleosynthesis and the value that we observe today. This observation is one of the very few we have from along our past worldline as opposed to our past lightcone. By releasing the untested assumption that the universe is homogeneous on very large scales, both apparent acceleration and the Lithium problem can be easily accounted for as different aspects of cosmic inhomogeneity, without causing problems for other cosmological phenomena such as the cosmic microwave background. We illustrate this in the context of a void model.

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.

Dark energy observations may be explained within general relativity using an inhomogeneous Hubble-scale depression in the matter density and accompanying curvature, which evolves naturally out of an Einstein-de Sitter (EdS) model. We present a simple parameterization of a void which can reproduce concordance model distances to arbitrary accuracy, but can parameterize away from this to give a smooth density profile everywhere. We show how the Hubble constant is not just a nuisance parameter in inhomogeneous models because it affects the shape of the distance-redshift relation. Independent Hubble-rate data from age estimates can in principle serve to break the degeneracy between concordance and void models, but the data is not yet able to achieve this. Using the latest Constitution supernova dataset we show that robust limits can be placed on the size of a void which is roughly independent of its shape. However, the sharpness of the profile at the origin cannot be well constrained due to supernova being dominated by peculiar velocities in the local Universe. We illustrate our results using some recently proposed diagnostics for the Friedmann models.

The dark energy problem has led to speculation that not only may LCDM be wrong, but that the FLRW models themselves may not even provide the correct family of background models. We discuss how direct measurements of H(z) can be used to formulate tests of the standard paradigm in cosmology. On their own, such measurements can be used to test for deviations from flat LCDM. When combined with supernovae distances, Hubble rate measurements provide a test of the Copernican principle and the homogeneity assumption of the standard model, which is independent of dark energy or metric based theory of gravity. A modification of this test also provides a model independent observable for flatness which decorrelates curvature determination from dark energy. We investigate these tests using Hubble rate measurements from age data, as well as from a Hubble rate inferred from recent measurements of the baryon acoustic oscillations. While the current data is too weak to say anything significant, these tests are exciting prospects for the future.

An important issue in cosmology is reconstructing the effective dark energy equation of state directly from observations. With so few physically motivated models, future dark energy studies cannot only be based on constraining a dark energy parameter space. We present a new non-parametric method which can accurately reconstruct a wide variety of dark energy behaviour with no prior assumptions about it. It is simple, quick and relatively accurate, and involves no expensive explorations of parameter space. The technique uses principal component analysis and a combination of information criteria to identify real features in the data, and tailors the fitting functions to pick up trends and smooth over noise. We find that we can constrain a large variety of w(z) models to within 10-20 % at redshifts z<1 using just SNAP-quality data.

We present a framework for studying gravitational lensing in spherically symmetric spacetimes using 1+1+2 covariant methods. A general formula for the deflection angle is derived and we show how this can be used to recover the standard result for the Schwarzschild spacetime.

I consider some of the issues we face in trying to understand dark energy. Huge fluctuations in the unknown dark energy equation of state can be hidden in distance data, so I argue that model-independent tests which signal if the cosmological constant is wrong are valuable. These can be constructed to remove degeneracies with the cosmological parameters. Gravitational effects can play an important role. Even small inhomogeneity clouds our ability to say something definite about dark energy. I discuss how the averaging problem confuses our potential understanding of dark energy by considering the backreaction from density perturbations to second-order in the concordance model: this effect leads to at least a 10% increase in the dynamical value of the deceleration parameter, and could be significantly higher. Large Hubble-scale inhomogeneity has not been investigated in detail, and could conceivably be the cause of apparent cosmic acceleration. I discuss void models which defy the Copernican principle in our Hubble patch, and describe how we can potentially rule out these models.

We investigate the effect that the average backreaction of structure formation has on the dynamics of the cosmological expansion, within the concordance model. Our approach in the Poisson gauge is fully consistent up to second-order in a perturbative expansion about a flat Friedmann background, including a cosmological constant. We discuss the key length scales which are inherent in any averaging procedure of this kind. We identify an intrinsic homogeneity scale that arises from the averaging procedure, beyond which a residual offset remains in the expansion rate and deceleration parameter. In the case of the deceleration parameter, this can lead to a quite large increase in the value - more than 10% - and may therefore have important ramifications for dark energy measurements, even if the underlying nature of dark energy is a cosmological constant. We give the intrinsic variance that affects the value of the effective Hubble rate and deceleration parameter. These considerations serve to add extra intrinsic errors to our determination of the cosmological parameters, and, in particular, may render attempts to measure the Hubble constant to percent precision overly optimistic.

We study gravitational waves in the black string Randall-Sundrum braneworld. We present a reasonably self-contained and complete derivation of the equations governing the evolution of gravitational perturbations in the presence of a brane localized source, and then specialize to the case of spherical radiation from a pointlike body in orbit around the black string. We solve for the resulting gravitational waveform numerically for a number of different orbital parameters.

The Lemaitre-Tolman-Bondi solution has received much attention as a possible alternative to Dark Energy, as it is able to account for the apparent acceleration of the Universe without any exotic matter content. However, in order to make rigorous comparisons between these models and cosmological observations, such as the integrated Sachs-Wolfe effect, baryon acoustic oscillations and the observed matter power spectrum, it is absolutely necessary to have a proper understanding of the linear perturbation theory about them. Here we present this theory in a fully general, and gauge-invariant form. It is shown that scalar, vector and tensor perturbations interact, and that the natural gauge invariant variables in Lemaitre-Tolman-Bondi cosmology do not correspond straightforwardly to the usual Bardeen variables, in the limit of spatial homogeneity. We therefore construct new variables that reduce to pure scalar, vector and tensor modes in this limit.

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.

We investigate what current cosmological data tells us about the cosmological expansion rate in a model independent way. Specifically, we study if the expansion was decelerating at high redshifts and is accelerating now, without referring to any model for the energy content of the universe, nor to any specific theory of gravity. This differs from most studies of the expansion rate which, e.g., assumes some underlying parameterised model for the dark energy component of the universe. To accomplish this, we have devised a new method to probe the expansion rate without relying on such assumptions. Using only supernova data, we conclude that there is little doubt that the universe has been accelerating at late times. However, contrary to some previous claims, we can not determine if the universe was previously decelerating. For a variety of methods used for constraining the expansion history of the universe, acceleration is detected from supernovae alone at >5 sigma, regardless of the curvature of the universe. Specifically, using a Taylor expansion of the scale factor, acceleration today is detected at >12 sigma. If we also include the ratio of the scale of the baryon acoustic oscillations as imprinted in the cosmic microwave background and in the large scale distribution of galaxies, it is evident from the data that the expansion decelerated at high redshifts, but only with the assumption of a flat or negatively curved universe.

The critical issue in cosmology today lies in determining if the cosmological constant is the underlying ingredient of dark energy. Our profound lack of understanding of the physics of dark energy places severe constrains on our ability to say anything about its possible dynamical nature. Quoted errors on the equation of state, w(z), are so heavily dependent on necessarily over-simplified parameterisations they are at risk of being rendered meaningless. Moreover, the existence of degeneracies between the reconstructed w(z) and the matter and curvature densities weakens any conclusions still further. We propose consistency tests for the cosmological constant which provide a direct observational signal if Lambda is wrong, regardless of the densities of matter and curvature. As an example of its utility, our flat case test can warn of a small transition from w(z)=-1 of 20% from SNAP quality data at 4-sigma, even when direct reconstruction techniques see virtually no evidence for deviation from Lambda. It is shown to successfully rule out a wide range of non-Lambda dark energy models with no reliance on knowledge of Omega_m using SNAP-quality data and a large range for using 10^5 supernovae as forecasted for LSST.

We present the time drift of the cosmological redshift in a general spherically symmetric spacetime. We demonstrate that its observation would allow us to test the Copernican principle and so determine if our universe is radially inhomogeneous, an important issue in our understanding of dark energy. In particular, when combined with distance data, this extra observable allows one to fully reconstruct the geometry of a spacetime describing a spherically symmetric under-dense region around us, purely from background observations.

The recent discovery of apparent cosmic acceleration has highlighted the depth of our ignorance of the fundamental properties of nature. It is commonly assumed that the explanation for acceleration must come from a new form of energy dominating the cosmos - dark energy - or a modification of Einstein's theory of Relativity. It is often overlooked, however, that a currently viable alternative explanation of the data is radial inhomogeneity which alters the Hubble diagram without any acceleration. This explanation is often ignored for two reasons: radial inhomogeneity significantly complicates analysis and predictions, and so the full details have not been investigated; and it is a philosophically highly controversial idea, revoking as it does the long-held Copernican Principle. To date, there has not been a general way of determining the validity if the Copernican Principle -- that we live at a typical position in the universe -- significantly weakening the foundations of cosmology as a scientific endeavour. Here we present an observational test for the Copernican assumption which can be automatically implemented while we search for dark energy in the coming decade. Our test is entirely independent of any model for dark energy or theory of gravity and thereby represents a model-independent test of the Copernican Principle.

While vector modes are usually ignored in cosmology since they are not produced during inflation they are inevitably produced from the interaction of density fluctuations of differing wavelengths. This effect may be calculated via a second-order perturbative expansion. We investigate this effect during the radiation era. We discuss the generation mechanism by investigating two scalar modes interacting, and we calculate the power of vector modes generated by a power-law spectrum of density perturbations on all scales.

We study the degeneracies between dark energy dynamics, dark matter and curvature using a non-parametric and non-perturbative approach. This allows us to examine the knock-on bias induced in the reconstructed dark energy equation of state, w(z), when there is a bias in the cosmic curvature or dark matter content, without relying on any specific parameterisation of w. Even assuming perfect Hubble, distance and volume measurements, we show that for z > 1, the bias in w(z) is up to two orders of magnitude larger than the corresponding errors in Omega_k or Omega_m. This highlights the importance of obtaining unbiased estimators of all cosmic parameters in the hunt for dark energy dynamics.

Quantum electrodynamics predicts that photons undergo one-loop scattering. The combined effect of this on the behaviour of a photon gas for temperatures above 10^{10} K results in a softening of the equation of state. We calculate the effect this has on the effective equation of state in the early universe, taking into account all the species of the Standard Model. The change to the dynamics of the early universe is discussed.

We show that the assumption of a flat universe induces critically large errors in reconstructing the dark energy equation of state at z>~0.9 even if the true cosmic curvature is very small, O(1%) or less. The spuriously reconstructed w(z) shows a range of unusual behaviour, including crossing of the phantom divide and mimicking of standard tracking quintessence models. For 1% curvature and LCDM, the error in w grows rapidly above z~0.9 reaching (50%,100%) by redshifts of (2.5,2.9) respectively, due to the long cosmological lever arm. Interestingly, the w(z) reconstructed from distance data and Hubble rate measurements have opposite trends due to the asymmetric influence of the curved geodesics. These results show that including curvature as a free parameter is imperative in any future analyses attempting to pin down the dynamics of dark energy, especially at moderate or high redshifts.

The generation of gravitational waves during inflation due to the non-linear coupling of scalar and tensor modes is discussed. Two methods describing gravitational wave perturbations are used and compared: a covariant and local approach, as well as a metric-based analysis based on the Bardeen formalism. An application to slow-roll inflation is also described.

We discuss the gravitational wave background generated by primordial density perturbations evolving during the radiation era. At second-order in a perturbative expansion, density fluctuations produce gravitational waves. We calculate the power spectra of gravitational waves from this mechanism, and show that, in principle, future gravitational wave detectors could be used to constrain the primordial power spectrum on scales vastly different from those currently being probed by large-scale structure. As examples we compute the gravitational wave background generated by both a power-law spectrum on all scales, and a delta-function power spectrum on a single scale.

We report on the possibility of detecting a submillimetre-sized extra dimension by observing gravitational waves (GWs) emitted by pointlike objects orbiting a braneworld black hole. Matter in the `visible' universe can generate a discrete spectrum of high frequency GWs with amplitudes moderately weaker than the predictions of general relativity (GR), while GW signals generated by matter on a `shadow' brane hidden in the bulk are potentially strong enough to be detected using current technology. We know of no other astrophysical phenomena that produces GWs with a similar spectrum, which stresses the need to develop detectors capable of measuring this high-frequency signature of large extra dimensions.

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.

The magnetohydrodynamic dynamo equation is derived within general relativity, using the covariant 1+3 approach, for a plasma with finite electric conductivity. This formalism allows for a clear division and interpretation of plasma and gravitational effects, and we have not restricted to a particular spacetime geometry. The results should be of interest in astrophysics and cosmology, and the formulation is well suited to gauge invariant perturbation theory. Moreover, the dynamo equation is presented in some specific limits. In particular, we consider the interaction of gravitational waves with magnetic fields, and present results for the evolution of the linearly growing electromagnetic induction field, as well as the diffusive damping of these fields.

We present a new approach to gauge-invariant cosmological perturbations at second order, which is also covariant. We examine two cases in particular for a dust Friedman-Lemaitre-Robertson-Walker model of any curvature: we investigate gravity waves generated from clustering matter, that is, induced tensor modes from scalar modes; and we discuss the generation of density fluctuations induced by gravity waves -- scalar modes from tensor perturbations. We derive a linear system of evolution equations for second-order gauge-invariant variables which characterise fully the induced modes of interest, with a source formed from variables quadratic in first-order quantities; these we transform into fully-fledged second-order gauge-invariant variables. Both the invariantly defined variables and the key evolution equations are considerably simpler than similar gauge-invariant results derived by other methods. By finding analytical solutions, we demonstrate that non-linear effects can significantly amplify or dampen modes present in standard linearised cosmological perturbation theory, thereby providing an important source of potential error in, and refinement of, the standard model. Moreover, these effects can dominate at late times, and on super-Hubble scales.

We investigate generation of electromagnetic radiation by gravitational waves interacting with a strong magnetic field in the vicinity of a vibrating Schwarzschild black hole. Such an effect may play an important role in gamma-ray bursts and supernovae, their afterglows in particular. It may also provide an electromagnetic counterpart to gravity waves in many situations of interest, enabling easier extraction and verification of gravity wave waveforms from gravity wave detection. We set up the Einstein-Maxwell equations for the case of odd parity gravity waves impinging on a static magnetic field as a covariant and gauge-invariant system of differential equations which can be integrated as an initial value problem, or analysed in the frequency domain. We numerically investigate both of these cases. We find that the black hole ringdown process can produce substantial amounts of electromagnetic radiation from a dipolar magnetic field in the vicinity of the photon sphere.

We present a new covariant and gauge-invariant perturbation formalism for dealing with spacetimes having spherical symmetry (or some preferred spatial direction) in the background, and apply it to the case of gravitational wave propagation in a Schwarzschild black hole spacetime. The 1+3 covariant approach is extended to a `1+1+2 covariant sheet' formalism by introducing a radial unit vector in addition to the timelike congruence, and decomposing all covariant quantities with respect to this. The background Schwarzschild solution is discussed and a covariant characterisation is given. We give the full first-order system of linearised 1+1+2 covariant equations, and we show how, by introducing (time and spherical) harmonic functions, these may be reduced to a system of first-order ordinary differential equations and algebraic constraints for the 1+1+2 variables which may be solved straightforwardly. We show how both the odd and even parity perturbations may be unified by the discovery of a covariant, frame- and gauge-invariant, transverse-traceless tensor describing gravitational waves, which satisfies a covariant wave equation equivalent to the Regge-Wheeler equation for both even and odd parity perturbations. We show how the Zerilli equation may be derived from this tensor, and derive a similar transverse traceless tensor equivalent to this equation. The so-called `special' quasinormal modes with purely imaginary frequency emerge naturally. The significance of the degrees of freedom in the choice of the two frame vectors is discussed, and we demonstrate that, for a certain frame choice, the underlying dynamics is governed purely by the Regge-Wheeler tensor. The two transverse-traceless Weyl tensors which carry the curvature of gravitational waves are 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

We study the asymptotic behaviour of scaling solutions with a dissipative fluid and we show that, contrary to recent claims, the existence of stable accelerating attractor solution which solves the `energy' coincidence problem depends crucially on the chosen equations of state for the thermodynamical variables. We discuss two types of equations of state, one which contradicts this claim, and one which supports it.

We study the qualitative properties of the class of spatially homogeneous Bianchi VI_o cosmological models containing a perfect fluid with a linear equation of state, a scalar field with an exponential potential and a uniform cosmic magnetic field, using dynamical systems techniques. We find that all models evolve away from an expanding massless scalar field model in which the matter and the magnetic field are negligible dynamically. We also find that for a particular range of parameter values the models evolve towards the usual power-law inflationary model (with no magnetic field) and, furthermore, we conclude that inflation is not fundamentally affected by the presence of a uniform primordial magnetic field. We investigate the physical properties of the Bianchi I magnetic field models in some detail.

Observations of the high degree of isotropy of the cosmic microwave background are commonly believed to indicate that the Universe is `almost' Friedmann-Lemaitre-Robertson-Walker (at least since the time of last scattering). Theoretical support for this belief comes from the so-called Ehlers-Geren-Sachs Theorem. We show that a generalization of this theorem rules out any strong magnetic fields in the Universe. Our theoretical result is model-independent and includes the case of inhomogeneous magnetic fields, complementing previous results. We thus prove that cosmic microwave background observations severely constrain all types of primordial and protogalactic magnetic fields in the universe.

We challenge the widely held belief that the cosmological principle is an obvious consequence of the observed isotropy of the cosmic microwave background radiation, combined with the Copernican principle. We perform a detailed analysis of a class of inhomogeneous perfect fluid cosmologies admitting an isotropic radiation field, with a view to assessing their viability as models of the real universe. These spacetimes are distinguished from FLRW universes by the presence of inhomogeneous pressure, which results in an acceleration of the fundamental observers. We examine their physical, geometrical and observational characteristics for all observer positions in the spacetimes. To this end, we derive exact, analytic expressions for the distance-redshift relations and anisotropies for any observer, and compare their predictions with available observational constraints. The isotropy constraints derived from `local' observations are also considered, qualitatively. A crucial aspect of this work is the application of the Copernican principle: for a specific model to be acceptable we demand that it must be consistent with current observational constraints (especially anisotropy constraints) for all observer locations. The most important results of the paper are presented as exclusion plots in the 2-D parameter space of the models. We show that there is a region of parameter space not ruled out by the constraints we consider and containing models that are significantly inhomogeneous. It follows immediately from this that the cosmological principle cannot be assumed to hold on the basis of present observational constraints.

This thesis concerns the compatibility of inhomogeneous cosmologies with our present understanding of the universe. It is a problem of some interest to find the class of all relativistic cosmological models which are capable of providing a reasonable `fit' to the universe. This thesis, in some respects, is part of this process. We consider Stephani models, which are a generalisation of the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) models, which can be thought of as FLRW models with acceleration and pressure gradients. Thus these models generalise the `dust' assumption of standard cosmology. The crucial aspect of this work is the retention of the Copernican principle -- an assumption regarded by many as crucial to cosmology. It states that we are not at a special location in the universe. This is a vital aspect of the original work in this thesis: consideration of an inhomogeneous model, while retaining the Copernican principle has, as far as the author is aware, not been considered in detail before. We start by generalising the Ehlers-Geren-Sachs Theorem to identify the class of inhomogeneous spacetimes which allow an isotropic radiation field for all observers in the spacetime. We then investigate observational and physical aspects of these models from all observer locations. We conclude that there exist spacetimes which conform to present observational constraints (especially anisotropy constraints) for any location in the spacetime, while at the same time being significantly inhomogeneous; ie, not `almost-FLRW'.

We demonstrate that the high isotropy of the Cosmic Microwave Background (CMB), combined with the Copernican principle, is not sufficient to prove homogeneity of the universe -- in contrast to previous results on this subject. The crucial additional factor not included in earlier work is the acceleration of the fundamental observers. We find the complete class of irrotational perfect fluid spacetimes admitting an exactly isotropic radiation field for every fundamental observer and show that are FLRW if and only if the acceleration is zero. While inhomogeneous in general, these spacetimes all possess three-dimensional symmetry groups, from which it follows that they also admit a thermodynamic interpretation. In addition to perfect fluids models we also consider multi-component fluids containing non-interacting radiation, dust and a quintessential scalar field or cosmological constant in which the radiation is isotropic for the geodesic (dust) observers. It is shown that the non-acceleration of the fundamental observers forces these spacetimes to be FLRW. While it is plausible that fundamental observers (galaxies) in the real universe follow geodesics, it is strictly necessary to determine this from local observations for the cosmological principle to be more than an assumption. We discuss how observations may be used to test this.