Source author record

Chris Clarkson

Chris Clarkson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

64works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

64 published item(s)

preprint2026arXiv

Probing General Relativity on Cosmological Scales in the 2040s

General relativity is exquisitely tested in strong-field regimes, yet its validity on cosmological scales remains largely unexplored. Upcoming wide and deep large-scale structure surveys will access the ultra-large, linear scales where relativistic effects - Doppler terms, gravitational redshift, lensing magnification, and potential evolution - leave significant imprints in the clustering of galaxies. These signatures represent unique probes of spacetime that are inaccessible to standard Newtonian analyses but increasingly important as survey volumes grow. We outline the scientific potential of next-generation facilities, such as those envisioned within ESO's Expanding Horizons programme, to deliver the first robust measurements of relativistic effects in large-scale structure through multi-tracer power spectra and the single-tracer bispectrum of high-redshift Lyman-break galaxies. Detecting these contributions would open a new window on gravity, enabling precision tests of general relativity and its alternatives on cosmological scales in the 2040s.

preprint2025arXiv

Probing the Cosmological Principle with CMB lensing and cosmic shear

The standard cosmological model assumes the Cosmological Principle. However, recent observations hint at possible violations of isotropy on large scales, possibly through late-time anisotropic expansion. Here we investigate the potential of cross-correlations between CMB lensing convergence $κ$ and galaxy cosmic shear $B$-modes as a novel probe of such late-time anisotropies. Our signal-to-noise forecasts reveal that information from the $κ$-$B$ cross-correlation is primarily contained on large angular scales ($\ell \lesssim 200$). We find that this cross-correlation for a Euclid-like galaxy survey is sensitive to anisotropy at the percent level. Making use of tomography yields a modest improvement of $\sim 20\%$ in detection power. Incorporating the galaxy $E$-$B$ cross-correlations would further enhance these constraints.

preprint2025arXiv

Transverse Velocities in Real-Time Cosmology: Position Drift in Relativistic N-Body Simulations

The era of real-time cosmology has begun. It is now possible to directly measure the apparent drift of high-redshift astronomical sources across the sky $\textit{in real time}$. This so-called $\textit{position drift}$ provides a valuable probe of the peculiar velocity field and cosmic structure formation by giving direct access to the transverse velocity, which is notoriously difficult to measure and is typically inferred statistically from the density field in a model-dependent way. To fully exploit this new window into the Universe, it is essential to understand how cosmological structures affect position drift measurements. Here we present the first position drift study based on the general relativistic N-body simulation code $\texttt{gevolution}$. We calculate the position drift directly from the past light cone for ten different observers and compare the results to predictions from linear perturbation theory. At linear order, the position drift is directly proportional to the transverse velocity on the sky. This linear approximation reproduces our non-linear simulation results to within about 5%. We calculate power spectra for the position drift, splitting the signal into an E- and B-mode and compare the former to linear expectations, finding good agreement. The B-mode is suppressed on linear scales, but has similar amplitude as the E-mode on non-linear scales. We further demonstrate that light-cone inhomogeneities induce biases in the dipole of the drift, introducing redshift dependence of both the amplitude and direction. Although our analysis is not yet sufficient for a firm conclusion, our results suggest that these effects alone cannot explain the possible redshift-dependent dipole in Gaia DR3 data reported in the literature.

preprint2022arXiv

A void in the Hubble tension? The end of the line for the Hubble bubble

The Universe may feature large-scale inhomogeneities beyond the standard paradigm, implying that statistical homogeneity and isotropy may be reached only on much larger scales than the usually assumed $\sim$100 Mpc. This means that we are not necessarily typical observers and that the Copernican principle could be recovered only on super-Hubble scales. Here, we do not assume the validity of the Copernican principle and let Cosmic Microwave Background, Baryon Acoustic Oscillations, type Ia supernovae, local $H_0$, cosmic chronometers, Compton y-distortion and kinetic Sunyaev-Zeldovich observations constrain the geometrical degrees of freedom of the local structure, which we parametrize via the $Λ$LTB model -- basically a non-linear radial perturbation of a FLRW metric. In order to quantify if a non-Copernican structure could explain away the Hubble tension, we pay careful attention to computing the Hubble constant in an inhomogeneous universe, and we adopt model selection via both the Bayes factor and the Akaike information criterion. Our results show that, while the $Λ$LTB model can successfully explain away the $H_0$ tension, it is favored with respect to the $Λ$CDM model only if one solely considers supernovae in the redshift range that is used to fit the Hubble constant, that is, $0.023<z<0.15$. If one considers all the supernova sample, then the $H_0$ tension is not solved and the support for the $Λ$LTB model vanishes. Combined with other data sets, this solution to the Hubble tension barely helps. Finally, we have reconstructed our local spacetime. We have found that data are best fit by a shallow void with $δ_L \approx -0.04$ and $r^{\mathrm{out}}_L \approx 300$ Mpc, which, interestingly, lies on the border of the 95\% credible region relative to the standard model expectation.

preprint2022arXiv

Magnification and evolution biases in large-scale structure surveys

Measurements of galaxy clustering in upcoming surveys such as those planned for the Euclid and Roman satellites, and the SKA Observatory, will be sensitive to distortions from lensing magnification and Doppler effects, beyond the standard redshift-space distortions. The amplitude of these contributions depends sensitively on magnification bias and evolution bias in the galaxy number density. Magnification bias quantifies the change in the observed number of galaxies gained or lost by lensing magnification, while evolution bias quantifies the physical change in the galaxy number density relative to the conserved case. These biases are given by derivatives of the number density, and consequently are very sensitive to the form of the luminosity function. We give a careful derivation of the magnification and evolution biases, clarifying a number of results in the literature. We then examine the biases for a variety of surveys, encompassing galaxy surveys and line intensity mapping at radio and optical/near-infrared wavelengths.

preprint2021arXiv

Observing relativistic features in large-scale structure surveys -- I: Multipoles of the power spectrum

Planned efforts to probe the largest observable distance scales in future cosmological surveys are motivated by a desire to detect relic correlations left over from inflation, and the possibility of constraining novel gravitational phenomena beyond General Relativity (GR). On such large scales, the usual Newtonian approaches to modelling summary statistics like the power spectrum and bispectrum are insufficient, and we must consider a fully relativistic and gauge-independent treatment of observables such as galaxy number counts in order to avoid subtle biases, e.g. in the determination of the $f_{\rm NL}$ parameter. In this work, we present an initial application of an analysis pipeline capable of accurately modelling and recovering relativistic spectra and correlation functions. As a proof of concept, we focus on the non-zero dipole of the redshift-space power spectrum that arises in the cross-correlation of different mass bins of dark matter halos, using strictly gauge-independent observable quantities evaluated on the past light cone of a fully relativistic N-body simulation in a redshift bin $1.7 \le z \le 2.9$. We pay particular attention to the correct estimation of power spectrum multipoles, comparing different methods of accounting for complications such as the survey geometry (window function) and evolution/bias effects on the past light cone, and discuss how our results compare with previous attempts at extracting novel GR signatures from relativistic simulations.

preprint2020arXiv

A null test to probe the scale-dependence of the growth of structure as a test of General Relativity

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.

preprint2020arXiv

Multi-Scale Perturbation Theory I: Methodology and Leading-Order Bispectrum Corrections in the Matter-Dominated Era

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.

preprint2020arXiv

Multipoles of the relativistic galaxy bispectrum

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.

preprint2020arXiv

Recursion relations for gravitational lensing

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.

preprint2020arXiv

The Hubble constant tension with next-generation galaxy surveys

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σ$) 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$σ$ level.

preprint2020arXiv

Weak-lensing observables in relativistic N-body simulations

We present a numerical weak-lensing analysis that is fully relativistic and non-perturbative for the scalar part of the gravitational potential and first-order in the vector part, frame dragging. Integrating the photon geodesics backwards from the observer to the emitters, we solve the Sachs optical equations and study in detail the weak-lensing convergence, ellipticity and rotation. For the first time, we apply such an analysis to a high-resolution relativistic N-body simulation, which consistently includes the leading-order corrections due to general relativity on both large and small scales. These are related to the question of gauge choice and to post-Newtonian corrections, respectively. 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.

preprint2019arXiv

Probing beyond-Horndeski gravity on ultra-large scales

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.

preprint2016arXiv

Dodging the dark matter degeneracy while determining the dynamics of dark energy

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.

preprint2016arXiv

Lensing and time-delay contributions to galaxy correlations

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.

preprint2016arXiv

Observed galaxy number counts on the lightcone up to second order: I. Main result

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.

preprint2016arXiv

Observed galaxy number counts on the lightcone up to second order: II. Derivation

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.

preprint2016arXiv

Roulettes: A weak lensing formalism for strong lensing - I. Overview

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

preprint2016arXiv

Roulettes: A weak lensing formalism for strong lensing - II. Derivation and analysis

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.

preprint2015arXiv

Beyond the plane-parallel and Newtonian approach: Wide-angle redshift distortions and convergence in general relativity

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.

preprint2015arXiv

Cosmological ensemble and directional averages of observables

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.

preprint2015arXiv

Cosmological evolution of the gravitational entropy of the large-scale structure

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ître 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 $Λ$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.

preprint2015arXiv

Do we care about the distance to the CMB? Clarifying the impact of second-order lensing

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.

preprint2015arXiv

Does small scale structure significantly affect cosmological dynamics?

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.

preprint2015arXiv

Testing foundations of modern cosmology with SKA all-sky surveys

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

preprint2015arXiv

The general theory of secondary weak gravitational lensing

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.

preprint2014arXiv

Cosmology with Doppler Lensing

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.

preprint2014arXiv

Evidence for a Lower Value for $H_0$ from Cosmic Chronometers Data?

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σ)$, 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.

preprint2014arXiv

Evolution of linear perturbations in spherically symmetric dust spacetimes

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.

preprint2014arXiv

Nonlinear relativistic corrections to cosmological distances, redshift and gravitational lensing magnification. I - Key results

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.

preprint2014arXiv

Nonlinear relativistic corrections to cosmological distances, redshift and gravitational lensing magnification. II - Derivation

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.

preprint2014arXiv

Observational Constraints on the Averaged Universe

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.

preprint2014arXiv

Spherically symmetric cosmological spacetimes with dust and radiation - numerical implementation

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.

preprint2014arXiv

The Value of $H_0$ from Gaussian Processes

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.

preprint2014arXiv

What is the distance to the CMB?

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.

preprint2013arXiv

Anti-lensing: the bright side of voids

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.

preprint2013arXiv

Galaxy correlations and the BAO in a void universe: structure formation as a test of the Copernican Principle

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.

preprint2013arXiv

Null tests of the cosmological constant using supernovae

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.

preprint2013arXiv

Optimising Gaussian processes for reconstructing dark energy dynamics from supernovae

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.

preprint2013arXiv

Testing the Copernican principle by constraining spatial homogeneity

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.

preprint2013arXiv

The Effect of Weak Lensing on Distance Estimates from Supernovae

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.

preprint2012arXiv

(Mis-)Interpreting supernovae observations in a lumpy universe

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.

preprint2012arXiv

Establishing homogeneity of the universe in the shadow of dark energy

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?

preprint2012arXiv

Reconstruction of dark energy and expansion dynamics using 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 .

preprint2012arXiv

Testing Homogeneity with Galaxy Star Formation Histories

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

preprint2012arXiv

The isotropic blackbody CMB as evidence for a homogeneous universe

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.

preprint2012arXiv

Using H(z) data as a probe of the concordance model

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.

preprint2011arXiv

Is backreaction really small within concordance cosmology?

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.

preprint2011arXiv

Local gauge-invariance at any order in cosmological perturbation theory

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.

preprint2011arXiv

Locally extracting scalar, vector and tensor modes in cosmological perturbation theory

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.

preprint2011arXiv

The Cosmic Microwave Background in an Inhomogeneous Universe - why void models of dark energy are only weakly constrained by the CMB

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.

preprint2011arXiv

The Hubble rate in averaged cosmology

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.

preprint2010arXiv

Direct reconstruction of dark energy

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.

preprint2010arXiv

Do primordial Lithium abundances imply there's no Dark Energy?

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.

preprint2010arXiv

Inhomogeneity and the foundations of concordance cosmology

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.

preprint2010arXiv

Model independent tests of the standard cosmological model

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.

preprint2010arXiv

Rendering Dark Energy Void

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.

preprint2009arXiv

On the determination of dark energy

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.

preprint2009arXiv

Perturbation Theory in Lemaitre-Tolman-Bondi Cosmology

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.

preprint2009arXiv

The cosmological background of vector modes

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.

preprint2005arXiv

Gravity-wave detectors as probes of 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.

preprint2004arXiv

Cosmological density fluctuations and gravity waves: A covariant approach to gauge-invariant non-linear perturbation theory

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.