In this paper we study the observational constraints that can be imposed on the coupling parameter, $\hat \alpha$, of the recently proposed 4-dimensional Einstein-Gauss-Bonnet theory of gravity. We consider the field equations that can be derived from regularization or dimensional reduction, and use them to perform a thorough investigation of the slow-motion and weak-field limit of this theory. We apply our results to observations of a wide array of physical systems that admit such a description, and find that the LAGEOS satellites are the most constraining, requiring $| \hat \alpha | \lesssim 10^{10} \,{\rm m}^2$. We then consider constraints that can be imposed on this theory from cosmology, black hole systems, and table-top experiments. It is found that early universe inflation prohibits all but the smallest negative values of $\hat \alpha$, while observations of binary black hole systems are likely to offer the tightest constraints on positive values, leading to overall bounds $0 \lesssim \hat \alpha \lesssim 10^8 \, {\rm m}^2$.

A variety of gauges are used in cosmological perturbation theory. These are often chosen in order to attribute physical properties to a particular choice of coordinates, or otherwise to simplify the form of the resultant equations. Calculations are then performed with the understanding that they could have been done in any gauge, and that transformations between different gauges can be made at will. We show that this logic can be extended to the domain of large density contrasts, where different types of perturbative expansion are required, but that the way in which gauges can be chosen in the presence of such structures is severely constrained. In particular, most gauges that are commonly considered in the cosmology literature are found to be unviable in the presence of non-linear structures. This includes spatially flat gauge, synchronous gauge, comoving orthogonal gauge, total matter gauge, N-body gauge, and the uniform density gauge. In contrast, we find that the longitudinal gauge and the Newtonian motion gauge are both viable choices in both standard cosmological perturbation theory, and in the post-Newtonian perturbative expansions that are required in order to model non-linear structures.

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.

We examine the implementation of Buchert's and Green & Wald's averaging formalisms in exact spherically symmetric and plane symmetric dust-filled cosmological models. We find that, given a cosmological space-time, Buchert's averaging scheme gives a faithful way of interpreting the large-scale expansion of space, and explicit terms that precisely quantify deviations from the behaviour expected from the Friedmann equations of homogeneous and isotropic cosmological models. The Green & Wald formalism, on the other hand, does not appear to yield any information about the large-scale properties of a given inhomogeneous space-time. Instead, this formalism is designed to calculate the back-reaction effects of short-wavelength fluctuations around a given "background" geometry. We find that the inferred expansion of space in this approach is entirely dependent on the choice of this background, which is not uniquely specified for any given inhomogeneous space-time, and that the "back-reaction" from small-scale structures vanishes in every case we study. This would appear to limit the applicability of Green & Wald's formalism to the study of large-scale expansion in the real Universe, which also has no pre-defined background. Further study is required to enhance the evaluation and comparison of these averaging formalisms, and determine whether the same difficulties exist, in less idealized space-time geometries.

We consider the problem of building inhomogeneous cosmological models in scalar-tensor theories of gravity. This starts by splitting the field equations of these theories into constraint and evolution equations, and then proceeds by identifying exact solutions to the constraints. We find exact, closed form expressions for geometries that correspond to the initial data for cosmological models containing regular arrays of point-like masses. These solutions extend similar methods that have recently been applied to Einstein's equations, and provides sufficient initial conditions to perform numerical integration of the evolution equations. We use our new solutions to study the effects of inhomogeneity in cosmologies governed by scalar-tensor theories of gravity, including the spatial inhomogeneity allowed in Newton's constant. Finally, we compare our solutions to their general relativistic counterparts, and investigate the effect of changing the coupling constant between the scalar and tensor degrees of freedom.

We consider a new variant of cosmological perturbation theory that has been designed specifically to include non-linear density contrasts on scales 100 Mpc, while still allowing for linear fluctuations on larger scales. This theory is used to derive the relativistic equations of Eulerian hydrodynamics in realistic cosmological scenarios that contain radiation and a cosmological constant, as well as matter that has been allowed to clump into galaxies and clusters of galaxies. These equations can be used to evolve energy densities and velocities in the presences of small-scale non-linear structures, and on scales all the way up to the horizon and beyond. The leading-order part of these equations reproduces the expected Newtonian equations, while subsequent orders prescribe relativistic corrections. We demonstrate that these evolution equations are consistent with maintaining the Einstein constraints, and hence that the system as a whole is mathematically well posed. The relativistic corrections that we derive are found to exhibit non-trivial interactions between perturbations on different scales, as well as the mixing of scalar, vector and tensor modes. They deviate from those that occur in both post-Friedmann and standard cosmological perturbation theory approaches, and point towards new relativistic effects that could be measurable by upcoming ultra-large-scale surveys.

We present a link between parametrizations of alternative theories of gravity on large and small scales in cosmology. This relationship is established using theoretical consistency conditions only. We find that in both limits the "slip" and "effective Newton's constant" can be written in terms of a set of four functions of time, two of which are direct generalizations of the $\alpha$ and $\gamma$ parameters from post-Newtonian physics. This generalizes previous work that has constructed frameworks for testing gravity on small scales, and is to the best of our knowledge the first time that a link between parametrizations of gravity on such very different scales has been established. We expect our result to facilitate the imposition of observational constraints, by drastically reducing the number of functional degrees of freedom required to consistently test gravity on multiple scales in cosmology.

The next generation of cosmological surveys will operate over unprecedented scales, and will therefore provide exciting new opportunities for testing general relativity. The standard method for modelling the structures that these surveys will observe is to use cosmological perturbation theory for linear structures on horizon-sized scales, and Newtonian gravity for non-linear structures on much smaller scales. We propose a two-parameter formalism that generalizes this approach, thereby allowing interactions between large and small scales to be studied in a self-consistent and well-defined way. This uses both post-Newtonian gravity and cosmological perturbation theory, and can be used to model realistic cosmological scenarios including matter, radiation and a cosmological constant. We find that the resulting field equations can be written as a hierarchical set of perturbation equations. At leading-order, these equations allow us to recover a standard set of Friedmann equations, as well as a Newton-Poisson equation for the inhomogeneous part of the Newtonian energy density in an expanding background. For the perturbations in the large-scale cosmology, however, we find that the field equations are sourced by both non-linear and mode-mixing terms, due to the existence of small-scale structures. These extra terms should be expected to give rise to new gravitational effects, through the mixing of gravitational modes on small and large scales - effects that are beyond the scope of standard linear cosmological perturbation theory. We expect our formalism to be useful for accurately modelling gravitational physics in universes that contain non-linear structures, and for investigating the effects of non-linear gravity in the era of ultra-large-scale surveys.

Motivated by the existence of hierarchies of structure in the Universe, we present four new families of exact initial data for inhomogeneous cosmological models at their maximum of expansion. These data generalise existing black hole lattice models to situations that contain clusters of masses, and hence allow the consequences of cosmological structures to be considered in a well-defined and non-perturbative fashion. The degree of clustering is controlled by a parameter $\lambda$, in such a way that for $\lambda \sim 0$ or $1$ we have very tightly clustered masses, whilst for $\lambda \sim 0.5$ all masses are separated by cosmological distance scales. We study the consequences of structure formation on the total net mass in each of our clusters, as well as calculating the cosmological consequences of the interaction energies both within and between clusters. The locations of the shared horizons that appear around groups of black holes, when they are brought sufficiently close together, are also identified and studied. We find that clustering can have surprisingly large effects on the scale of the cosmology, with models that contain thousands of black holes sometimes being as little as 30% of the size of comparable Friedmann models with the same total proper mass. This deficit is comparable to what might be expected to occur from neglecting gravitational interaction energies in Friedmann cosmology, and suggests that these quantities may have a significant influence on the properties of the large-scale cosmology.

On small scales the observable Universe is highly inhomogeneous, with galaxies and clusters forming a complex web of voids and filaments. The optical properties of such configurations can be quite different from the perfectly smooth Friedmann-Lema\^itre-Robertson-Walker (FLRW) solutions that are frequently used in cosmology, and must be well understood if we are to make precise inferences about fundamental physics from cosmological observations. We investigate this problem by calculating redshifts and luminosity distances within a class of cosmological models that are constructed explicitly in order to allow for large density contrasts on small scales. Our study of optics is then achieved by propagating one hundred thousand null geodesics through such space-times, with matter arranged in either compact opaque objects or diffuse transparent haloes. We find that in the absence of opaque objects, the mean of our ray tracing results faithfully reproduces the expectations from FLRW cosmology. When opaque objects with sizes similar to those of galactic bulges are introduced, however, we find that the mean of distance measures can be shifted up from FLRW predictions by as much as $10\%$. This bias is due to the viable photon trajectories being restricted by the presence of the opaque objects, which means that they cannot probe the regions of space-time with the highest curvature. It corresponds to a positive bias of order $10\%$ in the estimation of $\Omega_{\Lambda}$ and highlights the important consequences that astronomical selection effects can have on cosmological observables.

We consider the problem of constructing cosmological solutions of the Einstein-Maxwell equations that contain multiple charged black holes. By considering the field equations as a set of constraint and evolution equations, we construct exact initial data for N charged black holes on a hypersphere. This corresponds to the maximum of expansion of a cosmological solution, and provides sufficient information for a unique evolution. We then consider the specific example of a universe that contains eight charged black holes, and show that the existence of non-zero electric charge reduces the scale of the cosmological region of the space. These solutions generalize the Majumdar-Papapetrou solutions away from the extremal limit of charged black holes, and provide what we believe to be some of the first relativistic calculations of the effects of electric charge on cosmological backreaction.

In this paper we explore the idea that black holes can persist in a universe that collapses to a big crunch and then bounces into a new phase of expansion. We use a scalar field to model the matter content of such a universe {near the time} of the bounce, and look for solutions that represent a network of black holes within a dynamical cosmology. We find exact solutions to Einstein's constraint equations that provide the geometry of space at the minimum of expansion and that can be used as initial data for the evolution of hyperspherical cosmologies. These solutions illustrate that there exist models in which multiple distinct black holes can persist through a bounce, and allow for concrete computations of quantities such as the black hole filling factor. We then consider solutions in flat cosmologies, as well as in higher-dimensional spaces (with up to nine spatial dimensions). We derive conditions for the black holes to remain distinct (i.e. avoid merging) and hence persist into the new expansion phase. Some potentially interesting consequences of these models are also discussed.

We construct exact initial data for closed cosmological models filled with regularly arranged black holes in the presence of $\Lambda$. The intrinsic geometry of the 3-dimensional space described by this data is a sum of simple closed-form expressions, while the extrinsic curvature is just proportional to $\Lambda$. We determine the mass of each of the black holes in this space by performing a limiting procedure around the location of each of the black holes, and then compare the result to an appropriate slice through the Schwarzschild-de Sitter spacetime. The consequences of the inhomogeneity of this model for the large-scale expansion of space are then found by comparing the lengths of curves in the cosmological region to similar curves in a suitably chosen Friedmann-Lemaitre-Robertson-Walker (FLRW) solution. Finally, we locate the positions of the apparent horizons of the black holes, and determine the extremal values of their mass, for every possible regular arrangement of masses. We find that as the number of black holes is increased, the large-scale expansion of space approaches that of an FLRW model filled with dust and $\Lambda$, and that the extremal values of the black hole masses approaches that of the Schwarzschild-de Sitter solution.

We construct high-precision models of the Universe that contain radiation, a cosmological constant, and periodically distributed inhomogeneous matter. The density contrasts in these models are allowed to be highly non-linear, and the cosmological expansion is treated as an emergent phenomenon. This is achieved by employing a generalised version of the post-Newtonian formalism, and by joining together inhomogeneous regions of space-time at reflection symmetric junctions. Using these models, we find general expressions that precisely and unambiguously quantify the effect of small-scale inhomogeneity on the large-scale expansion of space (an effect referred to as "back-reaction", in the literature). We then proceed to specialize our models to the case where the matter fields are given by a regular array of point-like particles. This allows us to derive extremely simple expressions for the emergent Friedmann-like equations that govern the large-scale expansion of space. It is found that the presence of radiation tends to reduce the magnitude of back-reaction effects, while the existence of a cosmological constant has only a negligible effect.

Einstein's theory of gravity has been extensively tested on solar system scales, and for isolated astrophysical systems, using the perturbative framework known as the parameterized post-Newtonian (PPN) formalism. This framework is designed for use in the weak-field and slow-motion limit of gravity, and can be used to constrain a large class of metric theories of gravity with data collected from the aforementioned systems. Given the potential of future surveys to probe cosmological scales to high precision, it is a topic of much contemporary interest to construct a similar framework to link Einstein's theory of gravity and its alternatives to observations on cosmological scales. Our approach to this problem is to adapt and extend the existing PPN formalism for use in cosmology. We derive a set of equations that use the same parameters to consistently model both weak fields and cosmology. This allows us to parameterize a large class of modified theories of gravity and dark energy models on cosmological scales, using just four functions of time. These four functions can be directly linked to the background expansion of the universe, first-order cosmological perturbations, and the weak-field limit of the theory. They also reduce to the standard PPN parameters on solar system scales. We illustrate how dark energy models and scalar-tensor and vector-tensor theories of gravity fit into this framework, which we refer to as "parameterized post-Newtonian cosmology" (PPNC).

We examine the effect that the magnetic part of the Weyl tensor has on the large-scale expansion of space. This is done within the context of a class of cosmological models that contain regularly arranged discrete masses, rather than a continuous perfect fluid. The natural set of geodesic curves that one should use to consider the cosmological expansion of these models requires the existence of a non-zero magnetic part of the Weyl tensor. We include this object in the evolution equations of these models by performing a Taylor series expansion about a hypersurface where it initially vanishes. At the same cosmological time, measured as a fraction of the age of the universe, we find that the influence of the magnetic part of the Weyl tensor increases as the number of masses in the universe is increased. We also find that the influence of the magnetic part of the Weyl tensor increases with time, relative to the leading-order electric part, so that its contribution to the scale of the universe can reach values of ~1%, before the Taylor series approximation starts to break down.

We propose and construct a two-parameter perturbative expansion around a Friedmann-Lema\^{i}tre-Robertson-Walker geometry that can be used to model high-order gravitational effects in the presence of non-linear structure. This framework reduces to the weak-field and slow-motion post-Newtonian treatment of gravity in the appropriate limits, but also includes the low-amplitude large-scale fluctuations that are important for cosmological modelling. We derive a set of field equations that can be applied to the late Universe, where non-linear structure exists on supercluster scales, and perform a detailed investigation of the associated gauge problem. This allows us to identify a consistent set of perturbed quantities in both the gravitational and matter sectors, and to construct a set of gauge-invariant quantities that correspond to each of them. The field equations, written in terms of these quantities, take on a relatively simple form, and allow the effects of small-scale structure on the large-scale properties of the Universe to be clearly identified. We find that inhomogeneous structures source the global expansion, that there exist new field equations at new orders, and that there is vector gravitational potential that is a hundred times larger than one might naively expect from cosmological perturbation theory. Finally, we expect our formalism to be of use for calculating relativistic effects in upcoming ultra-large-scale surveys, as the form of the gravitational coupling between small and large scales depends on the non-linearity of Einstein's equations, and occurs at what is normally thought of as first order in cosmological perturbations.

In his 2005 review, Gravity and the Thermodynamics of Horizons, Paddy suggested that a vacuum in thermal equilibrium with a bath of radiation should have a gradually diminishing energy. We work through the consequences of this scenario, and find that a coupling between the vacuum and a bath of black-body radiation at the temperature of the horizon requires the Hubble rate, $H$, to approach the same type of evolution as in the "intermediate inflation" scenario, with $H\propto t^{-1/3}$, rather than as a constant. We show that such behaviour does not conflict with observations when the vacuum energy is described by a slowly-rolling scalar field, and when the fluctuations in the scalar field are treated as in the "warm inflation" scenario. It does, however, change the asymptotic states of the universe. We find that the existence of the radiation introduces a curvature singularity at early times, where the energy densities in both the radiation and the vacuum diverge. Furthermore, we show that the introduction of an additional non-interacting perfect fluid into the space-time reveals that radiation can dominate over dust at late times, in contrast to what occurs in the standard cosmological model. Such a coupling can also lead to a negative vacuum energy becoming positive.

We develop a new approach to building cosmological models, in which small pieces of perturbed Minkowski space are joined together at reflection-symmetric boundaries in order to form a global, dynamical space-time. Each piece of this patchwork universe is described using post-Newtonian gravitational physics, with the large-scale expansion of the universe being an emergent phenomenon. This approach to cosmology does not require any assumptions about non-local averaging processes. Our framework clarifies the relation between the weak-field limit of general relativity, and the cosmological solutions that result from solving Einstein's equations with a set of symmetry assumptions. It also allows the effects of structure formation on the large-scale expansion of the universe to be investigated in an unambiguous way. As an explicit example, we use this formalism to investigate the cosmological behaviour of a large number of regularly arranged point-like masses. In this case we find that the large-scale expansion is well modelled by a Friedmann-like equation that contains terms that take the form of dust, radiation, and spatial curvature. The radiation term, while small compared to the dust term, is purely a result of the non-linearity of Einstein's equations.

Almost all models of the universe start by assuming that matter fields can be modelled as dust. In the real universe, however, matter is clumped into dense objects that are separated by regions of space that are almost empty. If we are to treat such a distribution of matter as being modelled as a fluid, in some average or coarse-grained sense, then there a number of questions that must be answered. One of the most fundamental of these is whether or not the interaction energy between masses should gravitate. If it does, then a dust-like description may not be sufficient. We would then need to ask how interaction energies should be calculated in cosmology, and how they should appear in the Friedmann-like equations that govern the large-scale behaviour of the universe. I will discuss some recent results that may shed light on these questions.

When using Einstein's equations, there exist a number of techniques for embedding non-linear structures in cosmological backgrounds. These include Swiss cheese models, in which spherically symmetric vacua are patched onto Friedmann solutions, and lattice models, in which weak-field regions are joined together directly. In this talk we will consider how these methods work in f(R) theories of gravity. We will show that their existence places constraints on the large-scale expansion of the universe, and that it may not always be possible to consider the Friedmann solutions and weak-field solutions of a theory independently from each other.

We consider cosmological modelling in $f(R)$ theories of gravity, using both top-down and bottom-up constructions. The top-down models are based on Robertson-Walker geometries, and the bottom-up constructions are built by patching together sub-horizon-sized regions of perturbed Minkowski space. Our results suggest that these theories do not provide a theoretically attractive alternative to the standard general relativistic cosmology. We find that the only $f(R)$ theories that can admit an observationally viable weak-field limit have large-scale expansions that are observationally indistinguishable from the Friedmann solutions of General Relativity with $\Lambda$. Such theories do not alleviate any of the difficulties associated with $\Lambda$, and cannot produce any new behaviour in the cosmological expansion without simultaneously destroying the Newtonian approximation to gravity on small scales.

We apply the method of images to the exact initial data for cosmological models that contain a number of regularly arranged discrete masses. This allows us to join cosmological regions together by throats, and to construct wormholes in the initial data. These wormholes allow for the removal of the asymptotically flat "flange" regions that would otherwise exist on the far side of black holes. The method of images also provides us with a way to investigate the definition of mass is cosmology, and the cosmological consequences of the gravitational interaction energies between massive objects. We find evidence that the interaction energies within clusters of massive objects do indeed appear to contribute to the total energy budget in the cosmological regions of the space-time.

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.

We consider a family of cosmological models in which all mass is confined to a regular lattice of identical black holes. By exploiting the reflection symmetry about planes that bisect these lattices into identical halves, we are able to consider the evolution of a number of geometrically distinguished surfaces that exist within each of them. We find that the evolution equations for the reflection symmetric surfaces can be written as a simple set of Friedmann-like equations, with source terms that behave like a set of interacting effective fluids. We then show that gravitational waves are effectively trapped within small chambers for all time, and are not free to propagate throughout the space-time. Each chamber therefore evolves as if it were in isolation from the rest of the universe. We call this phenomenon "piecewise silence".

We study the effects of inhomogeneities on the evolution of the Universe, by considering a range of cosmological models with discretized matter content. This is done using exact and fully relativistic methods that exploit the symmetries in and about submanifolds of spacetimes that themselves possess no continuous global symmetries. These methods allow us to follow the evolution of our models throughout their entire history, far beyond what has previously been possible. We find that while some space-like curves collapse to anisotropic singularities in finite time, others remain non-singular forever. The resulting picture is of a cosmological spacetime in which some behaviour remains close to Friedmann-like, while other behaviours deviate radically. In particular, we find that large-scale acceleration is possible without any violation of the energy conditions.

An important open question in cosmology is the degree to which the Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions of Einstein's equations are able to model the large-scale behaviour of the locally inhomogeneous observable universe. We investigate this problem by considering a range of exact n-body solutions of Einstein's constraint equations. These solutions contain discrete masses, and so allow arbitrarily large density contrasts to be modelled. We restrict our study to regularly arranged distributions of masses in topological 3-spheres. This has the benefit of allowing straightforward comparisons to be made with FLRW solutions, as both spacetimes admit a discrete group of symmetries. It also provides a time-symmetric hypersurface at the moment of maximum expansion that allows the constraint equations to be solved exactly. We find that when all the mass in the universe is condensed into a small number of objects (<10) then the amount of backreaction in dust models can be large, with O(1) deviations from the predictions of the corresponding FLRW solutions. When the number of masses is large (>100), however, then our measures of backreaction become small (<1%). This result does not rely on any averaging procedures, which are notoriously hard to define uniquely in general relativity, and so provides (to the best of our knowledge) the first exact and unambiguous demonstration of backreaction in general relativistic cosmological modelling. Discrete models such as these can therefore be used as laboratories to test ideas about backreaction that could be applied in more complicated and realistic settings.

We propose a thermodynamically motivated measure of gravitational entropy based on the Bel-Robinson tensor, which has a natural interpretation as the effective super-energy-momentum tensor of free gravitational fields. The specific form of this measure differs depending on whether the gravitational field is Coulomb-like or wave-like, and reduces to the Bekenstein-Hawking value when integrated over the interior of a Schwarzschild black hole. For scalar perturbations of a Robertson-Walker geometry we find that the entropy goes like the Hubble weighted anisotropy of the gravitational field, and therefore increases as structure formation occurs. This is in keeping with our expectations for the behaviour of gravitational entropy in cosmology, and provides a thermodynamically motivated arrow of time for cosmological solutions of Einstein's field equations. It is also in keeping with Penrose's Weyl curvature hypothesis.

The problem of matching different regions of spacetime in order to construct inhomogeneous cosmological models is investigated in the context of Lagrangian theories of gravity constructed from general analytic functions f(R), and from non-analytic theories with f(R)=R^n. In all of the cases studied, we find that it is impossible to satisfy the required junction conditions without the large-scale behaviour reducing to that expected from Einstein's equations with a cosmological constant. For theories with analytic f(R) this suggests that the usual treatment of weak-field systems as perturbations about Minkowski space may not be compatible with late-time acceleration driven by anything other than a constant term of the form f(0), which acts like a cosmological constant. In the absence of Minkowski space as a suitable background for weak-field systems, one must then choose and justify some other solution to perform perturbative analyses around. For theories with f(R)=R^n we find that no known spherically symmetric vacuum solutions can be matched to an expanding FLRW background. This includes the absence of any Einstein-Straus-like embeddings of the Schwarzschild exterior solution in FLRW spacetimes.

We introduce the concept of back-reaction in relativistic cosmological modeling. Roughly speaking, this can be thought of as the difference between the large-scale behaviour of an inhomogeneous cosmological solution of Einstein's equations, and a homogeneous and isotropic solution that is a best-fit to either the average of observables or dynamics in the inhomogeneous solution. This is sometimes paraphrased as `the effect that structure has of the large-scale evolution of the universe'. Various different approaches have been taken in the literature in order to try and understand back-reaction in cosmology. We provide a brief and critical summary of some of them, highlighting recent progress that has been made in each case.

We discuss the construction of cosmological models within the framework of Macroscopic Gravity (MG), which is a theory that models the effects of averaging the geometry of space-time on large scales. We find new exact spatially homogeneous and isotropic FLRW solutions to the MG field equations, and investigate large-scale perturbations around them. We find that any inhomogeneous perturbations to the averaged geometry are severely restricted, but that possible anisotropies in the correlation tensor can have dramatic consequences for the measurement of distances. These calculations are a first step within the MG approach toward developing averaged cosmological models to a point where they can be used to interpret real cosmological data, and hence to provide a working alternative to the "concordance" LCDM model.

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.

Current cosmological observations, when interpreted within the framework of a homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) model, strongly suggest that the Universe is entering a period of accelerating expansion. This is often taken to mean that the expansion of space itself is accelerating. In a general spacetime, however, this is not necessarily true. We attempt to clarify this point by considering a handful of local and non-local measures of acceleration in a variety of inhomogeneous cosmological models. Each of the chosen measures corresponds to a theoretical or observational procedure that has previously been used to study acceleration in cosmology, and all measures reduce to the same quantity in the limit of exact spatial homogeneity and isotropy. In statistically homogeneous and isotropic spacetimes, we find that the acceleration inferred from observations of the distance-redshift relation is closely related to the acceleration of the spatially averaged universe, but does not necessarily bear any resemblance to the average of the local acceleration of spacetime itself. For inhomogeneous spacetimes that do not display statistical homogeneity and isotropy, however, we find little correlation between acceleration inferred from observations and the acceleration of the averaged spacetime. This shows that observations made in an inhomogeneous universe can imply acceleration without the existence of dark energy.

In this review we present a thoroughly comprehensive survey of recent work on modified theories of gravity and their cosmological consequences. Amongst other things, we cover General Relativity, Scalar-Tensor, Einstein-Aether, and Bimetric theories, as well as TeVeS, f(R), general higher-order theories, Horava-Lifschitz gravity, Galileons, Ghost Condensates, and models of extra dimensions including Kaluza-Klein, Randall-Sundrum, DGP, and higher co-dimension braneworlds. We also review attempts to construct a Parameterised Post-Friedmannian formalism, that can be used to constrain deviations from General Relativity in cosmology, and that is suitable for comparison with data on the largest scales. These subjects have been intensively studied over the past decade, largely motivated by rapid progress in the field of observational cosmology that now allows, for the first time, precision tests of fundamental physics on the scale of the observable Universe. The purpose of this review is to provide a reference tool for researchers and students in cosmology and gravitational physics, as well as a self-contained, comprehensive and up-to-date introduction to the subject as a whole.

The apparent accelerating expansion of the Universe, determined from observations of distant supernovae, and often taken to imply the existence of dark energy, may alternatively be explained by the effects of a giant underdense void if we relax the assumption of homogeneity on large scales. Recent studies have made use of the spherically-symmetric, radially-inhomogeneous Lemaitre-Tolman-Bondi (LTB) models to derive strong constraints on this scenario, particularly from observations of the kinematic Sunyaev-Zel'dovich (kSZ) effect which is sensitive to large scale inhomogeneity. However, most of these previous studies explicitly set the LTB 'bang time' function to be constant, neglecting an important freedom of the general solutions. Here we examine these models in full generality by relaxing this assumption. We find that although the extra freedom allowed by varying the bang time is sufficient to account for some observables individually, it is not enough to simultaneously explain the supernovae observations, the small-angle CMB, the local Hubble rate, and the kSZ effect. This set of observables is strongly constraining, and effectively rules out simple LTB models as an explanation of dark energy.

We consider the optical properties of Lindquist-Wheeler (LW) models of the Universe. These models consist of lattices constructed from regularly arranged discrete masses. They are akin to the Wigner-Seitz construction of solid state physics, and result in a dynamical description of the large-scale Universe in which the global expansion is given by a Friedmann-like equation. We show that if these models are constructed in a particular way then the redshifts of distant objects, as well as the dynamics of the global space-time, can be made to be in good agreement with the homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions of Einstein's equations, at the level of <3% out to z~2. Angular diameter and luminosity distances, on the other hand, differ from those found in the corresponding FLRW models, while being consistent with the 'empty beam' approximation, together with the shearing effects due to the nearest masses. This can be compared with the large deviations found from the corresponding FLRW values obtained in a previous study that considered LW models constructed in a different way. We therefore advocate the improved LW models we consider here as useful constructions that appear to faithfully reproduce both the dynamical and observational properties of space-times containing discrete masses.

We consider a model of the Universe in which the matter content is in the form of discrete islands, rather than a continuous fluid. In the appropriate limits the resulting large-scale dynamics approach those of a Friedmann-Robertson-Walker (FRW) universe. The optical properties of such a space-time, however, do not. This illustrates the fact that the optical and `average' dynamical properties of a relativistic universe are not equivalent, and do not specify each other uniquely. We find the angular diameter distance, luminosity distance and redshifts that would be measured by observers in these space-times, using both analytic approximations and numerical simulations. While different from their counterparts in FRW, the effects found do not look like promising candidates to explain the observations usually attributed to the existence of Dark Energy. This incongruity with standard FRW cosmology is not due to the existence of any unexpectedly large structures or voids in the Universe, but only to the fact that the matter content of the Universe is not a continuous fluid.

It is well known that observations of the cosmic microwave background (CMB) are highly sensitive to the spatial curvature of the Universe, k. Here we find that what is in fact being tightly constrained by small angle fluctuations is spatial curvature near the surface of last scattering, and that if we allow k to be a function of position, rather than taking a constant value everywhere, then considerable spatial curvature is permissible within our own locale. This result is of interest for the giant void models that attempt to explain the supernovae observations without Dark Energy. We find voids models with a homogeneous big bang can be compatible with the observed small angle CMB, but only if they exist in a positively curved universe. To be compatible with local measurements of H_0, however, we find that a radially varying bang time is required.

We propose a new statistic that has been designed to be used in situations where the intrinsic dispersion of a data set is not well known: The Crossing Statistic. This statistic is in general less sensitive than `chi^2' to the intrinsic dispersion of the data, and hence allows us to make progress in distinguishing between different models using goodness of fit to the data even when the errors involved are poorly understood. The proposed statistic makes use of the shape and trends of a model's predictions in a quantifiable manner. It is applicable to a variety of circumstances, although we consider it to be especially well suited to the task of distinguishing between different cosmological models using type Ia supernovae. We show that this statistic can easily distinguish between different models in cases where the `chi^2' statistic fails. We also show that the last mode of the Crossing Statistic is identical to `chi^2', so that it can be considered as a generalization of `chi^2'.

We construct cosmological models consisting of large numbers of identical, regularly spaced masses. These models do not rely on any averaging procedures, or on the existence of a global Friedmann-Robertson-Walker (FRW) background. They are solutions of Einstein's equations up to higher order corrections in a perturbative expansion, and have large-scale dynamics that are well modelled by the Friedmann equation. We find that the existence of arbitrarily large density contrasts does not change either the magnitude or scale of the background expansion, at least when masses are regularly arranged, and up to the prescribed level of accuracy. We also find that while the local space-time geometry inside each cell can be described as linearly perturbed FRW, one could argue that a more natural description is that of perturbed Minkowski space (in which case the scalar perturbations are simply Newtonian potentials). We expect these models to be of use for understanding and testing ideas about averaging in cosmology, as well as clarifying the relationship between global cosmological dynamics and the static space-times associated with isolated masses.

We consider the post-Newtonian limit of a general class of bimetric theories of gravity, in which both metrics are dynamical. The established parameterised post-Newtonian approach is followed as closely as possible, although new potentials are found that do not exist within the standard frame-work. It is found that these theories can evade solar system tests of post-Newtonian gravity remarkably well. We show that perturbations about Minkowski space in these theories contain both massless and massive degrees of freedom, and that in general there are two different types of massive mode, each with a different mass parameter. If both of these masses are sufficiently large then the predictions of the most general class of theories we consider are indistinguishable from those of general relativity, up to post-Newtonian order in a weak field, low velocity expansion. In the limit that the massive modes become massless, we find that these general theories do not exhibit a van Dam-Veltman-Zakharov-like discontinuity in their gamma parameter, although there are discontinuities in other post-Newtonian parameters as the massless limit is approached. This smooth behaviour in gamma is due to the discontinuities from each of the two different massive modes cancelling each other out. Such cancellations cannot occur in special cases with only one massive mode, such as the Isham-Salam-Strathdee theory.

We consider spatially homogeneous and isotropic Friedmann-Robertson-Walker (FRW) solutions of Milgrom's recently proposed class of bimetric theories of gravity. These theories have two different regimes, corresponding to high and low acceleration. We find simple power-law matter dominated solutions in both, as well as solutions with spatial curvature, and exponentially expanding solutions. In the high acceleration limit these solutions behave like the FRW solutions of General Relativity, with a cosmological constant term that is of the correct order of magnitude to explain the observed accelerating expansion of the Universe. We find that solutions that remain in the high acceleration regime for their entire history, however, require non-baryonic dark matter fields, or extra interaction terms in their gravitational Lagrangian, in order to be observationally viable. The low acceleration regime also provides some scope to account for this deficit, with solutions that differ considerably from their general relativistic counterparts.

We consider the observational effects of a deficit angle, w, in the topology of the solar system and in the 'double pulsar' system PSR J0737-3039A/B. Using observations of the perihelion precession of Mercury, and the gravitational deflection of light due to the Sun, we constrain the magnitude of such a deficit angle in the solar system to be 2*pi*(1-w), with 0<(1-w)<10^(-9) at 95% confidence. We calculate the effects of a deficit angle on the periastron advance, geodetic precession rate and inclination angle of the double pulsar system and use the observational data to obtain the constraint 0<(1-w)<2.4*10^(-8) at 95% confidence. Although this result is weaker than the solar system bound, it is in a very different physical environment, where accumulating data is likely to lead to tighter constraints in the future.

The matter content of the Universe is strongly inhomogeneous on small scales. Motivated by this fact, we consider a model of the Universe that has regularly spaced discrete masses, rather than a continuous fluid. The optical properties of such space-times can differ considerably from the continuous fluid case, even if the 'average' dynamics are the same. We show that these differences have consequences for cosmological parameter estimation, and that fitting to recent supernovae observations gives a correction to the inferred value of Omega_Lambda of ~10%.

We consider the effects of large structures in the Universe on the Hubble diagram. This problem is treated non-linearly by considering a Swiss Cheese model of the Universe in which under-dense voids are represented as negatively curved regions of space-time. Exact expressions for luminosity distances and redshifts are used to investigate the non-linear effects of structure on the magnitudes of astrophysical sources. It is found that the intervening voids we consider, between the observer and source, produce changes in apparent magnitude of less than 0.012. Sources inside voids, however, can be affected considerably at redshifts below z~0.5. By averaging observable quantities over many randomly generated distributions of voids we find that the presence of these structures has the effect of introducing a dispersion around the mean, which itself can be displaced the background value. Observers in an inhomogeneous universe, who take averages of observables along many different lines of sight, may then introduce systematic biases, and under-estimate errors, if these effects are not taken into account. Estimates of the potential size of these effects are made using data from simulated large-scale structure.

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.

A fundamental presupposition of modern cosmology is the Copernican Principle; that we are not in a central, or otherwise special region of the Universe. Studies of Type Ia supernovae, together with the Copernican Principle, have led to the inference that the Universe is accelerating in its expansion. The usual explanation for this is that there must exist a `Dark Energy', to drive the acceleration. Alternatively, it could be the case that the Copernican Principle is invalid, and that the data has been interpreted within an inappropriate theoretical frame-work. If we were to live in a special place in the Universe, near the centre of a void where the local matter density is low, then the supernovae observations could be accounted for without the addition of dark energy. We show that the local redshift dependence of the luminosity distance can be used as a clear discriminant between these two paradigms. Future surveys of Type Ia supernovae that focus on a redshift range of ~0.1-0.4 will be ideally suited to test this hypothesis, and hence to observationally determine the validity of the Copernican Principle on new scales, as well as probing the degree to which dark energy must be considered a necessary ingredient in the Universe.

We consider the space-time associated with the evaporation of a black hole by quantum mechanical tunnelling events. It is shown that the surface through which tunnelling occurs is distinct from the global event horizon, and that this has consequences for the radiation reaching future null infinity. A spherical collapse process is modelled, and the radiation expected to be observed at future null infinity is calculated. It is shown that external observers witness an evaporation process that begins as the tunnelling surface is exposed, and ends as the collapsing object passes behind its event horizon. The sensitivity of emitted radiation to the collapse process is illustrated.

We consider the Friedmann-Robertson-Walker cosmologies of theories of gravity that generalise the Einstein-Hilbert action by replacing the Ricci scalar, R, with some function, f(R). The general asymptotic behaviour of these cosmologies is found, at both early and late times, and the effects of adding higher and lower powers of R to the Einstein-Hilbert action is investigated. The assumption that the highest powers of R should dominate the Universe's early history, and that the lowest powers should dominate its future is found to be inaccurate. The behaviour of the general solution is complicated, and while it can be the case that single powers of R dominate the dynamics at late times, it can be either the higher or lower powers that do so. It is also shown that it is often the lowest powers of R that dominate at early times, when approach to a bounce or a Tolman solution are generic possibilities. Various examples are considered, and both vacuum and perfect fluid solutions investigated.

We study the observable pulse profiles that can be generated from precessing pulsars. A novel coordinate system is defined to aid visualization of the observing geometry. Using this system we explore the different families of profiles that can be generated by simple, circularly symmetric beam shapes. An attempt is then made to fit our model to the observations of relativistic binary PSR B1913+16. It is found that while qualitatively similar pulse profiles can be produced, this minimal model is insufficient for an accurate match to the observational data. Consequently, we confirm that the emission beam of PSR B1913+16 must deviate from circular symmetry, as first reported by Weisberg and Taylor. However, the approximate fits obtained suggest that it may be sufficient to consider only minimal deviations from a circular beam in order to explain the data. We also comment on the applicability of our analysis technique to other precessing pulsars, both binary and isolated.

We determine the full post-Newtonian limit of theories of gravity that extend general relativity by replacing the Ricci scalar, R, in the generating Lagrangian by some analytic function, f(R). We restrict ourselves to theories that admit Minkowski space as a suitable background, and perform a perturbative expansion in the manner prescribed by the parameterised post-Newtonian formalism. Extra potentials are found to be present that are not accounted for in the usual treatment, and a discussion is provided on how they may be used to observationally distinguished these theories from general relativity at the post-Newtonian level.

We find the general behaviour of homogeneous and isotropic cosmological models in some fourth-order theories of gravity. Explicit, exact, general solutions are given for both empty universes and those filled with a perfect fluid. For the vacuum case, solutions are found with closed, open and flat geometries, whilst the perfect fluid solutions are all spatially flat. Both early and late-time limits are studied, and attractor behaviour towards simple power-law expansion is identified. Multiple solutions to the same theories, with the same matter content and topology are found. It is shown that these solutions exhibit great variety in their evolution.

We consider the cosmological dynamics associated with volume weighted measures of eternal inflation, in the Bousso-Polchinski model of the string theory landscape. We find that this measure predicts that observers are most likely to find themselves in low energy vacua with one flux considerably larger than the rest. Furthermore, it allows for a satisfactory anthropic explanation of the cosmological constant problem by producing a smooth, and approximately constant, distribution of potentially observable values of Lambda. The low energy vacua selected by this measure are often short lived. If we require anthropically acceptable vacua to have a minimum life-time of 10 billion years, then for reasonable parameters a typical observer should expect their vacuum to have a life-time of approximately 12 billion years. This prediction is model dependent, but may point toward a solution to the coincidence problem of cosmology.

We investigate homogeneous and isotropic oscillating cosmologies with multiple fluid components. Transfer of energy between these fluids is included in order to model the effects of non-equilibrium behavior on closed universes. We find exact solutions which display a range of new behaviors for the expansion scale factor. Detailed examples are studied for the exchange of energy from dust or scalar field into radiation. We show that, contrary to expectation, it is unlikely that such models can offer a physically viable solution to the flatness problem.

The string theory landscape consists of many metastable de Sitter vacua, populated by eternal inflation. Tunneling between these vacua gives rise to a dynamical system, which asymptotically settles down to an equilibrium state. We investigate the effects of sinks to anti-de Sitter space, and show how their existence can change probabilities in the landscape. Sinks can disturb the thermal occupation numbers that would otherwise exist in the landscape and may cause regions that were previously in thermal contact to be divided into separate, thermally isolated islands.

This work investigates alternative theories of gravity, the solutions to their field equations and the constraints that can be imposed upon them from observation and experiment. Specifically, we consider the cosmologies and spherically symmetric solutions that can be expected to result from scalar-tensor and fourth-order theories of gravity. We find exact cosmological solutions of various different kinds; isotropic and anisotropic, homogeneous and inhomogeneous. These solutions are used to investigate the behaviour of the Universes at both late and early times, to investigate the effects of corrections to general relativity on approach to an initial singularity and to look for effects which may be observable in the present day Universe. We use physical processes, such as the primordial nucleosynthesis of the light elements, to impose constraints upon any deviations from the standard model. Furthermore, we investigate the vacuum spherically symmetric solutions of these theories. This environment is of particular interest for considerations of the local effects of gravity, where the most accurate experiments and observations of gravitational phenomena can be performed. Exact solutions are obtained for this situation and their stability analysed. It is found that a variety of new behaviour is obtainable in these theories that was not previously possible in the standard model. This new behaviour allows us an extended framework in which to consider gravitational physics, and its cosmological consequences.

We provide a simple mathematical description of the exchange of energy between two fluids in an expanding Friedmann universe with zero spatial curvature. The evolution can be reduced to a single non-linear differential equation which we solve in physically relevant cases and provide an analysis of all the possible evolutions. Particular power-law solutions exist for the expansion scale factor and are attractors at late times under particular conditions. We show how a number of problems studied in the literature, such as cosmological vacuum energy decay, particle annihilation, and the evolution of a population of evaporating black holes, correspond to simple particular cases of our model. In all cases we can determine the effects of the energy transfer on the expansion scale factor. We also consider the situation in the presence of anti-decaying fluids and so called phantom fluids which violate the dominant energy conditions.

We consider the possibility of energy being exchanged between the scalar and matter fields in scalar-tensor theories of gravity. Such an exchange provides a new mechanism which can drive variations in the gravitational 'constant' G. We find exact solutions for the evolution of spatially flat Friedman-Roberston-Walker cosmologies in this scenario and discuss their behaviour at both early and late times. We also consider the physical consequences and observational constraints on these models.

We determine the general conditions for the existence of Godel, Einstein static, and de Sitter universes in gravity theories derived from a Lagrangian that is an arbitrary function of the scalar curvature and Ricci and Riemann curvature invariants. Explicit expressions for the solutions are found in terms of the parameters defining the Lagrangian. We also determine the conditions on the Lagrangian of the theory under which time-travel is allowed in the Godel universes.

We investigate the effect of deviations from general relativity on approach to the initial singularity by finding exact cosmological solutions to a wide class of fourth-order gravity theories. We present new anisotropic vacuum solutions of modified Kasner type and demonstrate the extent to which they are valid in the presence of non-comoving perfect-fluid matter fields. The infinite series of Mixmaster oscillations seen in general relativity will not occur in these solutions, except in unphysical cases.

We study the cosmological and weak-field properties of theories of gravity derived by extending general relativity by means of a Lagrangian proportional to $R^{1+\delta}$. This scale-free extension reduces to general relativity when $\delta \to 0$. In order to constrain generalisations of general relativity of this power class we analyse the behaviour of the perfect-fluid Friedmann universes and isolate the physically relevant models of zero curvature. A stable matter-dominated period of evolution requires $\delta >0$ or $\delta <-1/4$. The stable attractors of the evolution are found. By considering the synthesis of light elements (helium-4, deuterium and lithium-7) we obtain the bound $-0.017<\delta <0.0012.$ We evaluate the effect on the power spectrum of clustering via the shift in the epoch of matter-radiation equality. The horizon size at matter--radiation equality will be shifted by $\sim 1%$ for a value of $\delta \sim 0.0005.$ We study the stable extensions of the Schwarzschild solution in these theories and calculate the timelike and null geodesics. No significant bounds arise from null geodesic effects but the perihelion precession observations lead to the strong bound $\delta =2.7\pm 4.5\times 10^{-19}$ assuming that Mercury follows a timelike geodesic. The combination of these observational constraints leads to the overall bound $0\leq \delta <7.2\times 10^{-19}$ on theories of this type.

We have found new anisotropic vacuum solutions for the scale-invariant gravity theories which generalise Einstein's general relativity to a theory derived from the Lagrangian $R^{1+\delta}$. These solutions are expanding universes of Kasner form with an initial spacetime singularity at $t=0 $ and exist for $-1/2<\delta <1/4$ but possess different Kasner index relations to the classic Kasner solution of general relativity to which they reduce when $\delta =0$. These solutions are unperturbed by the introduction of non-comoving perfect-fluid matter motions if $p<\rho $ on approach to the singularity and should not exhibit an infinite sequence of chaotic Mixmaster oscillations when $\delta >0$.

We study here the effect of a varying G on the evolution of the early Universe and, in particular, on primordial nucleosynthesis. This variation of G is modelled using the Brans-Dicke theory as well as a more general class of scalar-tensor theories. Modified nucleosynthesis codes are used to investigate this effect and the results obtained are used to constrain the parameters of the theories. We extend previous studies of primordial nucleosynthesis in scalar-tensor theories by including effects which can cause a slow variation of G during radiation domination and by including a late-time accelerating phase to the Universe's history. We include a brief discussion on the epoch of matter-radiation equality in Brans-Dicke theory, which is also of interest for determining the positions of the cosmic microwave background power-spectrum peaks.

We study the inhomogeneous cosmological evolution of the Newtonian gravitational 'constant' G in the framework of scalar-tensor theories. We investigate the differences that arise between the evolution of G in the background universes and in local inhomogeneities that have separated out from the global expansion. Exact inhomogeneous solutions are found which describe the effects of masses embedded in an expanding FRW Brans-Dicke universe. These are used to discuss possible spatial variations of G in different regions. We develop the technique of matching different scalar-tensor cosmologies of different spatial curvature at a boundary. This provides a model for the linear and non-linear evolution of spherical overdensities and inhomogeneities in G. This allows us to compare the evolution of G and \dot{G} that occurs inside a collapsing overdense cluster with that in the background universe. We develop a simple virialisation criterion and apply the method to a realistic lambda-CDM cosmology containing spherical overdensities. Typically, far slower evolution of \dot{G} will be found in the bound virialised cluster than in the cosmological background. We consider the behaviour that occurs in Brans-Dicke theory and in some other representative scalar-tensor theories.