Researcher profile

Lucas Lombriser

Lucas Lombriser contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.