Source author record

Lucas Lombriser

Lucas Lombriser 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

46works
5topics
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

46 published item(s)

preprint2022arXiv

$N$-body simulations for parametrised modified gravity

We present $\texttt{MG-evolution}$, an $N$-body code simulating the cosmological structure formation for parametrised modifications of gravity. It is built from the combination of parametrised linear theory with a parametrisation of the deeply nonlinear cosmological regime extrapolated from modified spherical collapse computations that cover the range of known screening mechanisms. We test $\texttt {MG-evolution}$, which runs at the speed of conventional $Λ$CDM simulations, against a suit of existing exact model-specific codes, encompassing linearised and chameleon $f(R)$ gravity as well as the normal branch of the Dvali-Gabadadz-Porrati braneworld model, hence covering both large-field value and large-derivative screening effects. We compare the nonlinear power spectra produced by the parametrised and model-specific approaches over the full range of scales set by the box size and resolution of our simulations, $k=(0.05-2.5)$~h/Mpc, and for two redshift slices, $z=0$ and $z=1$. We find sub-percent to one-percent level recovery of all the power spectra generated with the model-specific codes for the full range of scales. $\texttt {MG-evolution}$ can be used for generalised and accurate tests of gravity and dark energy with the increasing wealth of high-precision cosmological survey data becoming available over the next decade.

preprint2022arXiv

Measuring the propagation speed of gravitational waves with LISA

The propagation speed of gravitational waves, $c_T$, has been tightly constrained by the binary neutron star merger GW170817 and its electromagnetic counterpart, under the assumption of a frequency-independent $c_T$. Drawing upon arguments from Effective Field Theory and quantum gravity, we discuss the possibility that modifications of General Relativity allow for transient deviations of $c_T$ from the speed of light at frequencies well below the band of current ground-based detectors. We motivate two representative Ansätze for $c_T(f)$, and study their impact upon the gravitational waveforms of massive black hole binary mergers detectable by the LISA mission. We forecast the constraints on $c_T(f)$ obtainable from individual systems and a population of sources, from both inspiral and a full inspiral-merger-ringdown waveform. We show that LISA will enable us to place stringent independent bounds on departures from General Relativity in unexplored low-frequency regimes, even in the absence of an electromagnetic counterpart.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2022arXiv

On the road to percent accuracy VI: the nonlinear power spectrum for interacting dark energy with baryonic feedback and massive neutrinos

Understanding nonlinear structure formation is crucial for fully exploring the data generated by stage IV surveys, requiring accurate modelling of the power spectrum. This is challenging for deviations from $Λ$CDM, but we must ensure that alternatives are well tested, to avoid false detections. We present an extension of the halo model reaction framework for interacting dark energy. We modify the halo model including the additional force present in the Dark Scattering model and implement it into ReACT. The reaction is combined with a pseudo spectrum from EuclidEmulator2 and compared to N-body simulations. Using standard mass function and concentration-mass relation, we find predictions to be 1 % accurate at $z=0$ up to $k=0.8~h/{\rm Mpc}$ for the largest interaction strength tested ($ξ=50$ b/GeV), improving to $2~h/{\rm Mpc}$ at $z=1$. For smaller interaction strength ($10$ b/GeV), we find 1 % agreement at $z=1$ up to scales above $3.5~h/{\rm Mpc}$, being close to $1~h/{\rm Mpc}$ at $z=0$. Finally, we improve our predictions with the inclusion of baryonic feedback and massive neutrinos and study degeneracies between the effects of these contributions and those of the interaction. Limiting the scales to where our modelling is 1 % accurate, we find a degeneracy between the interaction and feedback, but not with massive neutrinos. We expect the degeneracy with feedback to be resolvable by including smaller scales. This work represents the first analytical tool for calculating the nonlinear spectrum for interacting dark energy models.

preprint2022arXiv

Scalar and tensor gravitational waves

In dark-energy models where a scalar field is nonminimally coupled to the spacetime geometry, gravitational waves are expected to be supplemented with a scalar mode. Such scalar waves may interact with the standard tensor waves, thereby affecting their observed amplitude and polarization. Understanding the role of scalar waves is thus essential in order to design reliable gravitational-wave probes of dark energy and gravity beyond general relativity. In this article, we thoroughly investigate the propagation of scalar and tensor waves in the subset of Horndeski theories in which tensor waves propagate at the speed of light. We work at linear order in scalar and metric perturbations, in the eikonal regime, and for arbitrary scalar and spacetime backgrounds. We diagonalize the system of equations of motion and identify the physical tensor mode, which differs from the metric perturbation. We find that interactions between scalar and tensor waves generally depend on the scalar propagation speed. If the scalar waves are luminal or quasiluminal, then interactions are negligible. In the subluminal case, scalar-tensor interactions are effectively suppressed due to the incoherence of the wave's phases.

preprint2022arXiv

Scalar Cherenkov radiation from high-energy cosmic rays

As first noted by Robert Wagoner in the 1970s, if a scalar field is nonminimally coupled to the Ricci scalar and propagates at subluminal speeds, then there exists the possibility of scalar $\check{\hbox{C}}$erenkov radiation from a moving particle. The mere observation of high-energy cosmic rays could in principle rule out the existence of such scalar fields since any particle moving faster than scalar perturbations would lose energy in the form of scalar waves until it moves slower than those. We compute in detail the energy loss to scalar waves and find that it scales with the square of the ultra-violet (UV) cutoff frequency of the effective field theory (EFT) of gravity. For dark-energy-motivated EFTs, the UV cutoff can be low, in which case that energy loss could always be negligible. In contrast, if viewed as a covariant theory valid at all scales or as an EFT valid at higher energies, perhaps even all the way up to the Planck scale, as may be the case if motivated by quantum-gravity perspectives, then the energy loss to scalar waves may diverge or become dramatically large. In this case, high-energy cosmic rays of extragalactic origin stringently constrain any conformally coupled scalar fields with non-canonical kinetic terms, although a minimum scalar phase velocity is required to trust the EFT.

preprint2022arXiv

Seeking New Physics in Cosmology with Bayesian Neural Networks: Dark Energy and Modified Gravity

We study the potential of Bayesian Neural Networks (BNNs) to detect new physics in the dark matter power spectrum, concentrating here on evolving dark energy and modifications to General Relativity. After introducing a new technique to quantify classification uncertainty in BNNs, we train two BNNs on mock matter power spectra produced using the publicly available code $\tt{ReACT}$ in the $k$-range $\left(0.01 - 2.5\right) \, h \mathrm{Mpc}^{-1} $ and redshift bins $\left(0.1,0.478,0.783,1.5\right)$ with Euclid-like noise. The first network classifies spectra into five labels including $Λ$CDM, $f(R)$, $w$CDM, Dvali-Gabadaze-Porrati (DGP) gravity and a "random" class whereas the second is trained to distinguish $Λ$CDM from non-$Λ$CDM. Both networks achieve a comparable training, validation and test accuracy of $\sim 95\%$. Each network is also capable of detecting deviations from $Λ$CDM that were not included in the training set, demonstrated with spectra generated using the growth-index $γ$. We then quantify the constraining power of each network by computing the smallest deviation from $Λ$CDM such that the noise-averaged non-$Λ$CDM classification probability is at least $2σ$, finding these bounds to be $f_{R0} \lesssim 10^{-7}$, $Ω_{rc} \lesssim 10^{-2} $, $-1.05 \lesssim w_0 \lesssim 0.95 $, $-0.2 \lesssim w_a \lesssim 0.2 $, $0.52 \lesssim γ\lesssim 0.59 $. The bounds on $f(R)$ can be improved by training a specialist network to distinguish solely between $Λ$CDM and $f(R)$ power spectra which can detect a non-zero $f_{R0}$ at $\mathcal{O}\left(10^{-8}\right)$ with a confidence $>2σ$. We expect that further developments, such as the inclusion of smaller length scales or additional extensions to $Λ$CDM, will only improve the potential of BNNs to detect new physics using cosmological datasets.

preprint2021arXiv

Horndeski theories and beyond from higher dimensions

The Einstein-Hilbert action with a cosmological constant is the most general local four-dimensional action leading to second-order derivative equations of motion that are symmetric and divergence free. In higher dimensions, additional terms can appear. We investigate a generalised metric decomposition involving a scalar degree of freedom to express the higher-dimensional action as an effective four-dimensional scalar-tensor theory. From the higher-dimensional Ricci scalar alone and a subclass of our metric ansatz, we recover the subset of Horndeski theories with luminal speed of gravitational waves. More generally, beyond-Horndeski terms appear. When including a Gauss-Bonnet scalar in the higher-dimensional action, we generate contributions to all cubic-order second-derivative terms present in the degenerate higher-order scalar-tensor theory as well as higher-derivative terms beyond that. We discuss this technique as a way to generate healthy four-dimensional gravity theories with an extra scalar degree of freedom and outline further generalisations of our method.

preprint2021arXiv

Parameterised post-Newtonian formalism for the effective field theory of dark energy via screened reconstructed Horndeski theories

We bring together two popular formalisms which generically parameterise deviations from General Relativity on astrophysical and cosmological scales, namely the parameterised post-Newtonian (PPN) formalism and the effective field theory (EFT) of dark energy and modified gravity. These separate formalisms are successfully applied to independently perform tests of gravity in their respective regimes of applicability on vastly different length scales. Nonlinear screening mechanisms indeed make it imperative to probe General Relativity across a wide range of scales. For a comprehensive interpretation of the complementary measurements it is important to connect them to effectively constrain the vast gravitational model space. We establish such a connection within the framework of Horndeski scalar-tensor theories restricted to a luminal propagation speed of gravitational waves. This is possible via the reconstruction of the family of linearly degenerate covariant Horndeski actions from the set of EFT functions and the subsequent derivation of the PPN parameters from the reconstructed theory. We outline the required conditions which ensure a reconstructed Horndeski model possesses a screening mechanism that enables significant modifications on cosmological scales while respecting stringent astrophysical bounds. Employing a scaling method, we then perform the general post-Newtonian expansion of the reconstructed models to derive their PPN parameters $γ$ and $β$ in their screened regimes.

preprint2021arXiv

The effect of screening mechanisms on black hole binary inspiral waveforms

Scalar-tensor theories leaving significant modifications of gravity at cosmological scales rely on screening mechanisms to recover General Relativity (GR) in high-density regions and pass stringent tests with astrophysical objects. Much focus has been placed on the signatures of such modifications of gravity on the propagation of gravitational waves (GWs) through cosmological distances while typically assuming their emission from fully screened regions with the wave generation strictly abiding by GR. Here, we closely analyse the impact of screening mechanisms on the inspiral GW waveforms from compact sources by employing a scaling method that enables a post-Newtonian (PN) expansion in screened regimes. Particularly, we derive the leading-order corrections to a fully screened emission to first PN order in the near zone and we also compute the modifications in the unscreened radiation zone to second PN order. For a concrete example, we apply our results to a cubic Galileon model. The resulting GW amplitude from a binary black hole inspiral deviate from its GR counterpart at most by one part in $10^{2}$ for the modifications in the radiation zone and at most one part in $10^{11}$ due to next-order corrections to the fully screened near zone. We expect such modifications to be undetectable by the current generation of GW detectors, but the deviation is not so small as to remain undetectable in future experiments.

preprint2021arXiv

The Novel Probes Project -- Tests of Gravity on Astrophysical Scales

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

preprint2020arXiv

A local self-tuning mechanism for the cosmological constant

Recently the global variation of the Planck mass in the General Relativistic Einstein-Hilbert action was proposed as a self-tuning mechanism of the cosmological constant preventing vacuum energy from freely gravitating. We show that this global mechanism emerges for generic local scalar-tensor theories with additional coupling of the scalar field to the field strength of a three-form gauge field that turns the scalar field constant on the domain of the action. Evaluation of the resulting integral constraint equation over the observable Universe yields a self-consistent framework with General Relativistic field equations and arbitrary radiatively stable residual cosmological constant. We argue that the expectation value for this residual is in good agreement with the magnitude of the observed cosmic acceleration.

preprint2020arXiv

Consistency of the local Hubble constant with the cosmic microwave background

A significant tension has become manifest between the current expansion rate of our Universe measured from the cosmic microwave background by the Planck satellite and from local distance probes, which has prompted for interpretations of that as evidence of new physics. Within conventional cosmology a likely source of this discrepancy is identified here as a matter density fluctuation around the cosmic average of the 40 Mpc environment in which the calibration of Supernovae Type Ia separations with Cepheids and nearby absolute distance anchors is performed. Inhomogeneities on this scale easily reach 40% and more. In that context, the discrepant expansion rates serve as evidence of residing in an underdense region of $δ_{\rm env}\approx-0.5\pm0.1$. The probability for finding this local expansion rate given the Planck data lies at the 95% confidence level. Likewise, a hypothetical equivalent local data set with mean expansion rate equal to that of Planck, while statistically favoured, would not gain strong preference over the actual data in the respective Bayes factor. These results therefore suggest borderline consistency between the local and Planck measurements of the Hubble constant. Generally accounting for the environmental uncertainty, the local measurement may be reinterpreted as a constraint on the cosmological Hubble constant of $H_0=74.7^{+5.8}_{-4.2}$ km/s/Mpc. The current simplified analysis may be augmented with the employment of the full available data sets, an impact study for the immediate $\lesssim10$ Mpc environment of the distance anchors, more prone to inhomogeneities, as well as expansion rates measured by quasar lensing, gravitational waves, currently limited to the same 40 Mpc region, and local galaxy distributions.

preprint2020arXiv

Dark degeneracy I: Dynamical or interacting dark energy?

We revisit the dark degeneracy that arises from the Einstein equations relating geometry to the total cosmic substratum but not resolving its individual components separately. We establish the explicit conditions for the dark degeneracy in the fluid description of the dark sector. At the background level, this degeneracy can be formally understood in terms of a unified dark sector Equation of State (EoS) that depends both on the dynamical nature of the dark energy (DE) as well as on its interaction with the pressureless dark matter. For linear perturbations, the degeneracy arises for specified DE pressure perturbations (or sound speed, equivalently) and DE anisotropic stress. Specializing to the degeneracy between non-interacting dynamical DE and interacting vacuum DE models, we perform a parameter estimation analysis for a range of dynamical DE parametrizations, where for illustration we explicitly break the degeneracy at the linear level by adopting a luminal sound speed for both scenarios. We conduct this analysis using cosmological background data alone and in combination with Planck measurements of the cosmic microwave background radiation. We find that although the overall phenomenology between the dynamical DE and interacting approaches is similar, there are some intriguing differences. In particular, there is an ambiguity in the strength of constraints on $Ω_{m0}$ and $σ_8$, which are considerably weakened for interacting vacuum DE, indicating that the dark degeneracy can change the significance of tensions in cosmological parameters inferred from different data sets.

preprint2020arXiv

Distinguishing cosmologies using the turn-around radius near galaxy clusters

Outside galaxy clusters the competition between the inwards gravitational attraction and the outwards expansion of the Universe leads to a special radius of velocity cancellation, which is called the turn-around radius. Measurements of the turn-around radius hold promises of constraining cosmological parameters, and possibly even properties of gravity. Such a measurement is, however, complicated by the fact that the surroundings of galaxy clusters are not spherical, but instead are a complicated collection of filaments, sheets and voids. In this paper we use the results of numerically simulated universes to quantify realistic error-bars of the measurement of the turn-around radius. We find that for a $Λ$CDM cosmology these error-bars are typically of the order of $20\%$. We numerically simulate three different implementations of dark energy models and of a scalar dark sector interaction to address whether the turn-around radius can be used to constrain non-trivial cosmologies, and we find that only rather extreme models can be distinguished from a $Λ$CDM universe due to the large error-bars arising from the non-trivial cluster environments.

preprint2020arXiv

Horndeski gravity and standard sirens

Standard sirens have been proposed as probes of alternative theories of gravity, such as Horndeski models. Hitherto, all studies have been conducted on a homogeneous-isotropic cosmological background, which is unable to consistently account for realistic distributions of matter, and for inhomogeneities in the Horndeski scalar field. Yet, the latter are essential for screening mechanisms. In this article, we comprehensively analyze the propagation of Horndeski gravitational waves in an arbitrary background spacetime and scalar field. We restrict to the class of theories in which gravitational waves propagate at light speed, and work in the geometric-optics regime. We find that kinetic braiding only produces a nonphysical longitudinal mode, whereas conformal coupling affects the amplitude of the standard transverse modes but not their polarization. We confirm that any observable deviation from general relativity depends on the local value of the effective Planck mass at emission and reception of the wave. This result is interpreted as the conservation of the number of gravitons.

preprint2020arXiv

Hybrid $P_{\ell}(k)$: general, unified, non-linear matter power spectrum in redshift space

Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian Streaming Model (GSM), to calculate the quasi non-linear redshift space matter power spectrum multipoles. This employs a fully non-linear and theoretically general prescription for the matter power spectrum. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model for dark matter. When comparing the redshift space multipoles for halos, we find that the Gaussian approximation of the GSM with a linear bias and a free stochastic term, $N$, is competitive to the TNS model. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model.

preprint2020arXiv

Modelling the matter bispectrum at small scales in modified gravity

Future large-scale structure surveys will measure three-point statistics with high statistical significance. This will offer significant improvements on our understanding of gravity, provided we can model these statistics accurately. We assess the performance of several schemes for theoretical modelling of the matter bispectrum, including halo-model based approaches and fitting formulae. We compare the model predictions against N-body simulations, considering scales up to $k_{\rm max} = 4 h/{\rm Mpc}$, well into non-linear regime of structure formation. Focusing on the equilateral configuration, we conduct this analysis for three theories of gravity: general relativity, $f(R)$ gravity, and the DGP braneworld model. Additionally, we compute the lensing convergence bispectrum for these models. We find that all current modelling prescriptions in modified gravity, in particular for theories with scale-dependent linear growth, fail to attain the accuracy required by the precision of the Stage IV surveys such as \emph{Euclid}. Among these models, we find that a halo-model corrected fitting formula achieves the best overall performance.

preprint2020arXiv

Positivity bounds on reconstructed Horndeski models

Positivity bounds provide conditions that a consistent UV-completion exists for a quantum field theory. We examine their application to Horndeski gravity models reconstructed from the effective field theory (EFT) of dark energy. This enables us to assess whether particular phenomenological parameterizations of the EFT functions reconstruct theories that respect or violate the positivity bounds. We find that commonly adopted EFT parametrizations, cast in terms of the dark energy density or power laws of the scale factor, only satisfy the positivity bounds in non-trivial regions of the parameter space. We then examine parameterizations of the inherently stable EFT basis, constructed to avoid gradient and ghost instabilities by default. In stark contrast, in this basis the positivity bounds either only provide constraints in a-priori unrealistic regions of the parameter space or do not provide any constraints on parameter values at all. The application of positivity bounds to common parametrizations of the standard EFT functions can therefore lead to artificial conclusions that the region of viable Horndeski modifications of gravity is highly constrained. Our results provide a strong motivation, in addition to the default avoidance of theoretical instabilities, for instead adopting parametrizations of the inherently stable EFT basis when testing dark energy and modified gravity models with forthcoming cosmological survey data.

preprint2019arXiv

Emulators for the non-linear matter power spectrum beyond $Λ$CDM

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $Λ$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $Λ$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.

preprint2016arXiv

A parametrisation of modified gravity on nonlinear cosmological scales

Viable modifications of gravity on cosmological scales predominantly rely on screening mechanisms to recover Einstein's Theory of General Relativity in the Solar System, where it has been well tested. A parametrisation of the effects of such modifications in the spherical collapse model is presented here for the use of modelling the modified nonlinear cosmological structure. The formalism allows an embedding of the different screening mechanisms operating in scalar-tensor theories through large values of the gravitational potential or its first or second derivatives as well as of linear suppression effects or more general transitions between modified and Einstein gravity limits. Each screening or suppression mechanism is parametrised by a time, mass, and environment dependent screening scale, an effective modified gravitational coupling in the fully unscreened limit that can be matched to linear theory, the exponent of a power-law radial profile of the screened coupling, determined by derivatives, symmetries, and potentials in the scalar field equation, and an interpolation rate between the screened and unscreened limits. Along with generalised perturbative methods, the parametrisation may be used to formulate a nonlinear extension to the linear parametrised post-Friedmannian framework to enable generalised tests of gravity with the wealth of observations from the nonlinear cosmological regime.

preprint2016arXiv

Accurate halo-model matter power spectra with dark energy, massive neutrinos and modified gravitational forces

We present an accurate non-linear matter power spectrum prediction scheme for a variety of extensions to the standard cosmological paradigm, which uses the tuned halo model previously developed in Mead (2015b). We consider dark energy models that are both minimally and non-minimally coupled, massive neutrinos and modified gravitational forces with chameleon and Vainshtein screening mechanisms. In all cases we compare halo-model power spectra to measurements from high-resolution simulations. We show that the tuned halo model method can predict the non-linear matter power spectrum measured from simulations of parameterised $w(a)$ dark energy models at the few per cent level for $k<10\,h\mathrm{Mpc}^{-1}$, and we present theoretically motivated extensions to cover non-minimally coupled scalar fields, massive neutrinos and Vainshtein screened modified gravity models that result in few per cent accurate power spectra for $k<10\,h\mathrm{Mpc}^{-1}$. For chameleon screened models we achieve only 10 per cent accuracy for the same range of scales. Finally, we use our halo model to investigate degeneracies between different extensions to the standard cosmological model, finding that the impact of baryonic feedback on the non-linear matter power spectrum can be considered independently of modified gravity or massive neutrino extensions. In contrast, considering the impact of modified gravity and massive neutrinos independently results in biased estimates of power at the level of 5 per cent at scales $k>0.5\,h\mathrm{Mpc}^{-1}$. An updated version of our publicly available HMcode can be found at https://github.com/alexander-mead/HMcode

preprint2016arXiv

Breaking a Dark Degeneracy with Gravitational Waves

We identify a scalar-tensor model embedded in the Horndeski action whose cosmological background and linear scalar fluctuations are degenerate with the concordance cosmology. The model admits a self-accelerated background expansion at late times that is stable against perturbations with a sound speed attributed to the new field that is equal to the speed of light. While degenerate in scalar fluctuations, self-acceleration of the model implies a present cosmological tensor mode propagation at < 95% of the speed of light with a damping of the wave amplitude that is > 5% less efficient than in general relativity. We show that these discrepancies are endemic to self-accelerated Horndeski theories with degenerate large-scale structure and are tested with measurements of gravitational waves emitted by events at cosmological distances. Hence, gravitational-wave cosmology breaks the dark degeneracy in observations of the large-scale structure between two fundamentally different explanations of cosmic acceleration - a cosmological constant and a scalar-tensor modification of gravity. The gravitational wave event GW150914 recently detected with the aLIGO instruments and its potential association with a weak short gamma-ray burst observed with the Fermi GBM experiment may have provided this crucial measurement.

preprint2016arXiv

Challenges to Self-Acceleration in Modified Gravity from Gravitational Waves and Large-Scale Structure

With the advent of gravitational-wave astronomy marked by the aLIGO GW150914 and GW151226 observations, a measurement of the cosmological speed of gravity will likely soon be realized. We show that a confirmation of equality to the speed of light as indicated by indirect Galactic observations will have important consequences for a very large class of alternative explanations of the late-time accelerated expansion of our Universe. It will break the dark degeneracy of self-accelerated Horndeski scalar-tensor theories in the large-scale structure that currently limits a rigorous discrimination between acceleration from modified gravity and from a cosmological constant or dark energy. Signatures of a self-acceleration must then manifest in the linear, unscreened cosmological structure. We describe the minimal modification required for self-acceleration with standard gravitational-wave speed and show that its maximum likelihood yields a 3-sigma poorer fit to cosmological observations compared to a cosmological constant. Hence, equality between the speeds challenges the concept of cosmic acceleration from a genuine scalar-tensor modification of gravity.

preprint2016arXiv

Constraints on decaying early modified gravity from cosmological observations

Most of the information on our cosmos stems from either late-time observations or the imprint of early-time inhomogeneities on the cosmic microwave background. We explore to what extent early modifications of gravity, which become significant after recombination but then decay towards the present, can be constrained by current cosmological observations. For the evolution of the gravitational modification, we adopt the decaying mode of a hybrid-metric Palatini $f(\mathcal{R})$ gravity model which is designed to reproduce the standard cosmological background expansion history and due to the decay of the modification is naturally compatible with Solar-System tests. We embed the model in the effective field theory description of Horndeski scalar-tensor gravity with an early-time decoupling of the gravitational modification. Since the quasistatic approximation for the perturbations in the model breaks down at high redshifts, where modifications remain relevant, we introduce a computationally efficient correction to describe the evolution of the scalar field fluctuation in this regime. We compare the decaying early-time modification against geometric probes and recent Planck measurements and find no evidence for such effects in the observations. Current data constrains the scalar field value at $|f_{\mathcal{R}}(z=z_{\rm on})| \lesssim 10^{-2}$ for modifications introduced at redshifts $z_{\rm on}\sim(500-1000)$ with present-day value $|f_{\mathcal{R}0}|\lesssim10^{-8}$. Finally, we comment on constraints that will be achievable with future 21~cm surveys and gravitational wave experiments.

preprint2016arXiv

Dark Energy vs. Modified Gravity

Understanding the reason for the observed accelerated expansion of the Universe represents one of the fundamental open questions in physics. In cosmology, a classification has emerged among physical models for the acceleration, distinguishing between Dark Energy and Modified Gravity. In this review, we give a brief overview of models in both categories as well as their phenomenology and characteristic observable signatures in cosmology. We also introduce a rigorous distinction between Dark Energy and Modified Gravity based on the strong and weak equivalence principles.

preprint2016arXiv

Finding Horndeski theories with Einstein gravity limits

The Horndeski action is the most general scalar-tensor theory with at most second-order derivatives in the equations of motion, thus evading Ostrogradsky instabilities and making it of interest when modifying gravity at large scales. To pass local tests of gravity, these modifications predominantly rely on nonlinear screening mechanisms that recover Einstein's Theory of General Relativity in regions of high density. We derive a set of conditions on the four free functions of the Horndeski action that examine whether a specific model embedded in the action possesses an Einstein gravity limit or not. For this purpose, we develop a new and surprisingly simple scaling method that identifies dominant terms in the equations of motion by considering formal limits of the couplings that enter through the new terms in the modified action. This enables us to find regimes where nonlinear terms dominate and Einstein's field equations are recovered to leading order. Together with an efficient approximation of the scalar field profile, one can then further evaluate whether these limits can be attributed to a genuine screening effect. For illustration, we apply the analysis to both a cubic galileon and a chameleon model as well as to Brans-Dicke theory. Finally, we emphasise that the scaling method also provides a natural approach for performing post-Newtonian expansions in screened regimes.

preprint2015arXiv

Classifying Linearly Shielded Modified Gravity Models in Effective Field Theory

We study the model space generated by the time-dependent operator coefficients in the effective field theory of the cosmological background evolution and perturbations of modified gravity and dark energy models. We identify three classes of modified gravity models that reduce to Newtonian gravity on the small scales of linear theory. These general classes contain enough freedom to simultaneously admit a matching of the concordance model background expansion history. In particular, there exists a large model space that mimics the concordance model on all linear quasistatic subhorizon scales as well as in the background evolution. Such models also exist when restricting the theory space to operators introduced in Horndeski scalar-tensor gravity. We emphasize that whereas the partially shielded scenarios might be of interest to study in connection with tensions between large and small scale data, with conventional cosmological probes, the ability to distinguish the fully shielded scenarios from the concordance model on near-horizon scales will remain limited by cosmic variance. Novel tests of the large-scale structure remedying this deficiency and accounting for the full covariant nature of the alternative gravitational theories, however, might yield further insights on gravity in this regime.

preprint2015arXiv

Constraining chameleon models with cosmology

Chameleon fields may modify gravity on cluster scales while recovering general relativity locally. This article reviews signatures of chameleon modifications in the nonlinear cosmological structure, comparing different techniques to model them, summarising the current state of observational constraints, and concluding with an outlook on prospective constraints from future observations and applications of the analytic tools developed in the process to more general scalar-tensor theories. Particular focus is given to the Hu-Sawicki and designer models of f(R) gravity.

preprint2015arXiv

Rapid simulation rescaling from standard to modified gravity models

We develop and test an algorithm to rescale a simulated dark-matter particle distribution or halo catalogue from a standard gravity model to that of a modified gravity model. This method is based on that of Angulo & White but with some additional ingredients to account for (i) scale-dependent growth of linear density perturbations and (ii) screening mechanisms that are generic features of viable modified gravity models. We attempt to keep the method as general as possible, so that it may plausibly be applied to a wide range of modified theories, although tests against simulations are restricted to a subclass of $f(R)$ models at this stage. We show that rescaling allows the power spectrum of matter to be reproduced at the $\sim 3$ per cent level in both real and redshift space up to $k=0.1h\,\mathrm{Mpc}^{-1}$ if we change the box size and alter the particle displacement field; this limit can be extended to $k=1h\,\mathrm{Mpc}^{-1}$ if we additionally alter halo internal structure. We simultaneously develop an algorithm that can be applied directly to a halo catalogue, in which case the halo mass function and clustering can be reproduced at the $\sim 5$ per cent level. Finally we investigate the clustering of halo particle distributions, generated from rescaled halo catalogues, and find that a similar accuracy can be reached.

preprint2015arXiv

Semi-dynamical perturbations of unified dark energy

Linear cosmological perturbations of a large class of modified gravity and dark energy models can be unified in the effective field theory of cosmic acceleration, encompassing Horndeski scalar-tensor theories and beyond. The fully available model space inherent to this formalism cannot be constrained by measurements in the quasistatic small-scale regime alone. To facilitate the analysis of modifications from the concordance model beyond this limit, we introduce a semi-dynamical treatment extrapolated from the evolution of perturbations at a pivot scale of choice. At small scales, and for Horndeski theories, the resulting modifications recover a quasistatic approximation but account for corrections to it near the Hubble scale. For models beyond Horndeski gravity, we find that the velocity field and time derivative of the spatial metric potential can generally not be neglected, even in the small-scale limit. We test the semi-dynamical approximation against the linear perturbations of a range of dark energy and modified gravity models, finding good agreement between the two.

preprint2015arXiv

Unscreening Modified Gravity in the Matter Power Spectrum

Viable modifications of gravity that may produce cosmic acceleration need to be screened in high-density regions such as the Solar System, where general relativity is well tested. Screening mechanisms also prevent strong anomalies in the large-scale structure and limit the constraints that can be inferred on these gravity models from cosmology. We find that by suppressing the contribution of the screened high-density regions in the matter power spectrum, allowing a greater contribution of unscreened low densities, modified gravity models can be more readily discriminated from the concordance cosmology. Moreover, by variation of density thresholds, degeneracies with other effects may be dealt with more adequately. Specializing to chameleon gravity as a worked example for screening in modified gravity, employing N-body simulations of f(R) models and the halo model of chameleon theories, we demonstrate the effectiveness of this method. We find that a percent-level measurement of the clipped power at k < 0.3 h/Mpc can yield constraints on chameleon models that are more stringent than what is inferred from Solar System tests or distance indicators in unscreened dwarf galaxies. Finally, we verify that our method is also applicable to the Vainshtein mechanism.

preprint2014arXiv

Halo model and halo properties in Galileon gravity cosmologies

We investigate the performance of semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body simulations. We focus on the Cubic and Quartic Galileon models that provide a reasonable fit to CMB, SNIa and BAO data. We demonstrate that the Sheth-Tormen mass function and linear halo bias can be calibrated to provide a very good fit to our simulation results. We also find that the halo concentration-mass relation is well fitted by a power law. The nonlinear matter power spectrum computed in the halo model approach is found to be inaccurate in the mildly nonlinear regime, but captures reasonably well the effects of the Vainshtein screening mechanism on small scales. In the Cubic model, the screening mechanism hides essentially all of the effects of the fifth force inside haloes. In the case of the Quartic model, the screening mechanism leaves behind residual modifications to gravity, which make the effective gravitational strength time-varying and smaller than the standard value. Compared to normal gravity, this causes a deficiency of massive haloes and leads to a weaker matter clustering on small scales. For both models, we show that there are realistic halo occupation distributions of Luminous Red Galaxies that can match both the observed large-scale clustering amplitude and the number density of these galaxies.

preprint2014arXiv

Halo modelling in chameleon theories

We analyse modelling techniques for the large-scale structure formed in scalar-tensor theories of constant Brans-Dicke parameter which match the concordance model background expansion history and produce a chameleon suppression of the gravitational modification in high-density regions. Thereby, we use a mass and environment dependent chameleon spherical collapse model, the Sheth-Tormen halo mass function and linear halo bias, the Navarro-Frenk-White halo density profile, and the halo model. Furthermore, using the spherical collapse model, we extrapolate a chameleon mass-concentration scaling relation from a LCDM prescription calibrated to N-body simulations. We also provide constraints on the model parameters to ensure viability on local scales. We test our description of the halo mass function and nonlinear matter power spectrum against the respective observables extracted from large-volume and high-resolution N-body simulations in the limiting case of f(R) gravity, corresponding to a vanishing Brans-Dicke parameter. We find good agreement between the two; the halo model provides a good qualitative description of the shape of the relative enhancement of the f(R) matter power spectrum with respect to LCDM caused by the extra attractive gravitational force but fails to recover the correct amplitude. Introducing an effective linear power spectrum in the computation of the two-halo term to account for an underestimation of the chameleon suppression at intermediate scales in our approach, we accurately reproduce the measurements from the N-body simulations.

preprint2014arXiv

How chameleons core dwarfs with cusps

The presence of a scalar field that couples nonminimally and universally to matter can enhance gravitational forces on cosmological scales while restoring general relativity in the Solar neighborhood. In the intermediate regime, kinematically inferred masses experience an additional radial dependence with respect to the underlying distribution of matter, which is caused by the increment of gravitational forces with increasing distance from the Milky Way center. The same effect can influence the internal kinematics of subhalos and cause cuspy matter distributions to appear core-like. Specializing to the chameleon model as a worked example, we demonstrate this effect by tracing the scalar field from the outskirts of the Milky Way halo to its interior, simultaneously fitting observed velocity dispersions of chemo-dynamically discriminated red giant populations in the Fornax and Sculptor dwarf spheroidals. Whereas in standard gravity these observations suggest that the matter distribution of the dwarfs is cored, we find that in the presence of a chameleon field the assumption of a cuspy Navarro-Frenk-White profile becomes perfectly compatible with the data. Importantly, chameleon models also predict the existence of slopes between two stellar subcomponents that in Newtonian gravity would be interpreted as a depletion of matter in the dwarf center. Hence, an observation of such an apparently pathological scenario may serve as a smoking gun for the presence of a chameleon field or a similar modification of gravity, independent of baryonic feedback effects. In general, measuring the dynamic mass profiles of the Milky Way dwarfs provides stronger constraints than those inferred from the screening scale of the Solar System since these are located at greater distances from the halo center.

preprint2014arXiv

Testing chameleon gravity with the Coma cluster

We propose a novel method to test the gravitational interactions in the outskirts of galaxy clusters. When gravity is modified, this is typically accompanied by the introduction of an additional scalar degree of freedom, which mediates an attractive fifth force. The presence of an extra gravitational coupling, however, is tightly constrained by local measurements. In chameleon modifications of gravity, local tests can be evaded by employing a screening mechanism that suppresses the fifth force in dense environments. While the chameleon field may be screened in the interior of the cluster, its outer region can still be affected by the extra force, introducing a deviation between the hydrostatic and lensing mass of the cluster. Thus, the chameleon modification can be tested by combining the gas and lensing measurements of the cluster. We demonstrate the operability of our method with the Coma cluster, for which both a lensing measurement and gas observations from the X-ray surface brightness, the X-ray temperature, and the Sunyaev-Zel'dovich effect are available. Using the joint observational data set, we perform a Markov chain Monte Carlo analysis of the parameter space describing the different profiles in both the Newtonian and chameleon scenarios. We report competitive constraints on the chameleon field amplitude and its coupling strength to matter. In the case of f(R) gravity, corresponding to a specific choice of the coupling, we find an upper bound on the background field amplitude of |f_{R0}|<6*10^{-5}, which is currently the tightest constraint on cosmological scales.

preprint2013arXiv

Chameleon f(R) gravity in the virialized cluster

Current constraints on f(R) gravity from the large-scale structure are at the verge of penetrating into a region where the modified forces become nonlinearly suppressed. For a consistent treatment of observables at these scales, we study cluster quantities produced in chameleon and linearized Hu-Sawicki f(R) gravity dark matter N-body simulations. We find that the standard Navarro-Frenk-White halo density profile and the radial power law for the pseudo phase-space density provide equally good fits for f(R) clusters as they do in the Newtonian scenario. We give qualitative arguments for why this should be the case. For practical applications, we derive analytic relations, e.g., for the f(R) scalar field, the gravitational potential, and the velocity dispersion as seen within the virialized clusters. These functions are based on three degrees of freedom fitted to simulations, i.e., the characteristic density, scale, and velocity dispersion. We further analyze predictions for these fitting parameters from the gravitational collapse and the Jeans equation, which are found to agree well with the simulations. Our analytic results can be used to consistently constrain chameleon f(R) gravity with future observations on virialized cluster scales without the necessity of running a large number of simulations.

preprint2013arXiv

Modeling halo mass functions in chameleon f(R) gravity

On cosmological scales, observations of the cluster abundance currently place the strongest constraints on f(R) gravity. These constraints lie in the large-field limit, where the modifications of general relativity can correctly be modeled by setting the Compton wavelength of the scalar field to its background value. These bounds are, however, at the verge of penetrating into a regime, where the modifications become nonlinearly suppressed due to the chameleon mechanism and cannot be described by this linearized approximation. For future constraints based on observations subjected to cluster abundance, it is therefore essential to consistently model the chameleon effect. We analyze descriptions of the halo mass function in chameleon f(R) gravity using a mass- and environment-dependent spherical collapse model in combination with excursion set theory and phenomenological fits to N-body simulations in the LCDM and f(R) gravity scenarios. Our halo mass functions consistently incorporate the chameleon suppression and cosmological parameter dependencies, improving upon previous formalisms and providing an important extension to N-body simulations for the application in consistent tests of gravity with observables sensitive to the abundance of clusters.

preprint2013arXiv

Relativistic effects in galaxy clustering in a parametrized post-Friedmann universe

We explore the signatures of quintessence and modified gravity theories in the relativistic description of galaxy clustering within a parametrized post-Friedmann framework. For this purpose, we develop a calibration method to consistently account for horizon-scale effects in the linear parametrized Post-Friedmann perturbations of minimally and nonminimally coupled scalar-tensor theories and test it against the full model-specific fluctuations. We further study the relativistic effects in galaxy clustering for the normal and self-accelerating branches of the Dvali-Gabadadze-Porrati braneworld model as well as for phenomenological modifications of gravity. We quantify the impact of modified gravity and dark energy models on galaxy clustering by computing the velocity-to-matter density ratio F, the velocity contribution R, and the potential contribution P and give an estimate of their detectability in future galaxy surveys. Our results show that, in general, the relativistic correction contains additional information on gravity and dark energy, which needs to be taken into account in consistent horizon-scale tests of departures from LCDM using the galaxy-density field.

preprint2012arXiv

Cluster Density Profiles as a Test of Modified Gravity

We present a new test of gravitational interactions at the r\sim(0.2-20)Mpc scale, around the virial radius of dark matter halos measured through cluster-galaxy lensing of maxBCG clusters from the Sloan Digital Sky Survey (SDSS). We employ predictions from self-consistent simulations of f(R) gravity to find an upper bound on the background field amplitude of f_R0<3.5x10^-3 at the 1D-marginalized 95% confidence level. As a model-independent assessment of the constraining power of cluster profiles measured through weak gravitational lensing, we also constrain the amplitude F_0 of a phenomenological modification based on the profile enhancement induced by f(R) gravity when not including effects from the increased cluster abundance in f(R). In both scenarios, dark-matter-only simulations of the concordance model corresponding to f_R0=0 and F_0=0 are consistent with the lensing measurements, i.e., at the 68% and 95% confidence level, respectively.

preprint2012arXiv

The Clustering of Galaxies in the SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Testing Deviations from $Λ$ and General Relativity using anisotropic clustering of galaxies

We use the joint measurement of geometry and growth from anisotropic galaxy clustering in the Baryon Oscillation Spectroscopic Survey Data Release 9 CMASS sample reported by Reid et al. to constrain dark energy properties and possible deviations from the General Relativity. Assuming GR and taking a prior on the linear matter power spectrum at high redshift from the cosmic microwave background (CMB), anisotropic clustering of the CMASS DR9 galaxies alone constrains $Ω_{\rm m} = 0.308 \pm 0.022$ and $100Ω_{\rm k} = 5.9 \pm 4.8$ for $w = -1$, or $w = -0.91 \pm 0.12$ for $Ω_k = 0$. When combined with the full CMB likelihood, the addition of the anisotropic clustering measurements to the spherically-averaged BAO location increases the constraining power on dark energy by a factor of 4 in a flat CDM cosmology with constant dark energy equation of state $w$ (giving $w = -0.87 \pm 0.05$). This impressive gain depends on our measurement of both the growth of structure and Alcock-Paczynski effect, and is not realised when marginalising over the amplitude of redshift space distortions. Combining with both the CMB and Supernovae Type Ia (SNeIa), we find $Ω_{\rm m} = 0.281 \pm 0.014$ and $1000Ω_{\rm k}=-9.2\pm5.0$ for $w = -1$, or $w_0 = -1.13 \pm 0.12$ and $w_{\rm a}=0.65 \pm 0.36$ assuming $Ω_k = 0$. Finally, when a $Λ$CDM background expansion is assumed, the combination of our estimate of the growth rate with previous growth measurements provides tight constraints on the parameters describing possible deviations from GR giving $γ= 0.64 \pm 0.05$. For one parameter extensions of the flat $Λ$CDM model, we find a $\sim 2σ$ preference either for $w > -1$ or slower growth than in GR. However, the data is fully consistent with the concordance model, and the evidence for these additional parameters is weaker than $2σ$.

preprint2011arXiv

Consistency check of ΛCDM phenomenology

The standard model of cosmology LCDM assumes general relativity, flat space, and the presence of a positive cosmological constant. We relax these assumptions allowing spatial curvature, time-dependent effective dark energy equation of state, as well as modifications of the Poisson equation for the lensing potential, and modifications of the growth of linear matter density perturbations in alternate combinations. Using six parameters characterizing these relations, we check LCDM for consistency utilizing cosmic microwave background anisotropies, cross correlations thereof with high-redshift galaxies through the integrated Sachs-Wolfe effect, the Hubble constant, supernovae and baryon acoustic oscillation distances, as well as the relation between weak gravitational lensing and galaxy flows. In all scenarios, we find consistency of the concordance model at the 95% confidence level. However, we emphasize that constraining supplementary background parameters and parametrizations of the growth of large-scale structure separately may lead to a priori exclusion of viable departures from the concordance model.

preprint2011arXiv

Constraints on f(R) gravity from probing the large-scale structure

We study cosmological constraints on metric f(R) gravity models that are designed to reproduce the LCDM expansion history with modifications to gravity described by a supplementary cosmological freedom, the Compton wavelength parameter B_0. We conduct a Markov chain Monte Carlo analysis on the parameter space, utilizing the geometrical constraints from supernovae distances, the baryon acoustic oscillation distances, and the Hubble constant, along with all of the cosmic microwave background data, including the largest scales, its correlation with galaxies, and a probe of the relation between weak gravitational lensing and galaxy flows. The strongest constraints, however, are obtained through the inclusion of data from cluster abundance. Using all of the data, we infer a bound of B_0<0.0011 at the 95% C.L.

preprint2010arXiv

Confirmation of general relativity on large scales from weak lensing and galaxy velocities

Although general relativity underlies modern cosmology, its applicability on cosmological length scales has yet to be stringently tested. Such a test has recently been proposed, using a quantity, EG, that combines measures of large-scale gravitational lensing, galaxy clustering and structure growth rate. The combination is insensitive to 'galaxy bias' (the difference between the clustering of visible galaxies and invisible dark matter) and is thus robust to the uncertainty in this parameter. Modified theories of gravity generally predict values of EG different from the general relativistic prediction because, in these theories, the 'gravitational slip' (the difference between the two potentials that describe perturbations in the gravitational metric) is non-zero, which leads to changes in the growth of structure and the strength of the gravitational lensing effect3. Here we report that EG = 0.39 +/- 0.06 on length scales of tens of megaparsecs, in agreement with the general relativistic prediction of EG $\approx$ 0.4. The measured value excludes a model within the tensor-vector-scalar gravity theory, which modifies both Newtonian and Einstein gravity. However, the relatively large uncertainty still permits models within f(R) theory, which is an extension of general relativity. A fivefold decrease in uncertainty is needed to rule out these models.

preprint2010arXiv

Testing General Relativity with Current Cosmological Data

Deviations from general relativity, such as could be responsible for the cosmic acceleration, would influence the growth of large scale structure and the deflection of light by that structure. We clarify the relations between several different model independent approaches to deviations from general relativity appearing in the literature, devising a translation table. We examine current constraints on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. Markov Chain Monte Carlo likelihood analysis of the parameters over various redshift ranges yields consistency with general relativity at the 95% confidence level.

preprint2009arXiv

Cosmological Constraints on DGP Braneworld Gravity with Brane Tension

We perform a Markov Chain Monte Carlo analysis of the self-accelerating and normal branch of Dvali-Gabadadze-Porrati braneworld gravity. By adopting a parameterized post-Friedmann description of gravity, we utilize all of the cosmic microwave background data, including the largest scales, and its correlation with galaxies in addition to the geometrical constraints from supernovae distances and the Hubble constant. We find that on both branches brane tension or a cosmological constant is required at high significance with no evidence for the unique Dvali-Gabadadze-Porrati modifications. The cross-over scale must therefore be substantially greater than the Hubble scale H_0 r_c > 3 and 3.5 at the 95% CL with and without uncertainties from spatial curvature. With spatial curvature, the limit from the normal branch is substantially assisted by the galaxy cross-correlation which highlights its importance in constraining infrared modifications to gravity.