Search

Gravitational wave dark sirens are a powerful tool for cosmology and inference of compact object population hyperparameters. They allow for a measurement of the luminosity distance to the source, but not their redshift. Galaxy catalogues in the source localization volume can be used to infer the redshift of the source in a statistical manner. Catalogues are, however, limited by their incompleteness, which can be significant at redshifts corresponding to current GW events. In this work, we detail how to implement in practice variance completion, a novel galaxy completion method which uses knowledge of the large scale structure to optimize the potential of dark sirens analyses. We compress the prediction for the missing number of galaxies into a ratio between the predictions of variance completion and the standard homogeneous completion method. This ratio format can be easily incorporated into existing line of sight computations used in dark sirens software; we demonstrate this procedure using the GLADE+ galaxy catalogue and the gwcosmo software package. We discuss the robustness of the method, and apply it to well-localized event GW190814 as a proof of concept. Finally, we apply the method to data from the third observing run of LIGO-Virgo-KAGRA, finding that it yields results that are consistent with homogeneous completion. We also discuss the prospects for an improvement if the GW localization volume shrinks., Comment: 27 pages, 12 figures

The dark sirens method combines gravitational waves and catalogs of galaxies to constrain the cosmological expansion history, merger rates and mass distributions of compact objects, and the laws of gravity. However, the incompleteness of galaxy catalogs means faint potential host galaxies are unobserved, and must be modeled to avoid inducing a bias. The majority of dark sirens analyses to date assume that the missing galaxies are distributed uniformly across the sky, which is clearly unphysical. We introduce a new Bayesian approach to the reconstruction of galaxy catalogs, which makes full use of our knowledge of large-scale structure. Our method quantifies the uncertainties on the estimated true galaxy number count in each voxel, and is marginalized over cosmological parameters and bias parameters. Crucially, our method further assesses the (absolute) magnitude distribution of galaxies, which is not known from the galaxy catalog itself. We present the details of our method and validate our approach on a galaxy catalog associated to the Millennium Simulation. The tools developed here generate physically-informed and robust host galaxy reconstructions, enabling more informative dark sirens analyses. Stage IV galaxy surveys will display greater redshift overlap with GW observations, whilst remaining incomplete -- emphasizing the importance of our work., Comment: 43 pages, 15 figures

The $\texttt{Hi-COLA}$ code is an efficient dark matter simulation suite that flexibly handles the Horndeski family of modified gravity models. In this work we extend the scope of $\texttt{Hi-COLA}$ to accommodate Horndeski theories with K-mouflage screening, allowing for the computation of matter power spectra in the non-linear regime in these models. We explore the boost of the dark matter power spectrum relative to GR-$\Lambda$CDM in K-mouflage gravity, and also discuss how large-scale structure computations change between the Einstein and Jordan frames. A dissection of the relative contributions of the modified background, linear growth, fifth force, and the conformal factor (a new inclusion to $\texttt{Hi-COLA}$) to the boost factor is presented. The ability of $\texttt{Hi-COLA}$ to run with general Horndeski models and multiple screening mechanisms makes it an ideal tool for testing gravity with upcoming galaxy survey data., Comment: 31 pages, 10 figures, 1 table. Comments welcome. The Hi-COLA code can be found at https://github.com/Hi-COLACode/Hi-COLA

Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and K-mouflage. In particular, we compare the predictions from full $N$-body simulations, two $N$-body codes with approximate time integration schemes, a parametrised modified $N$-body implementation and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in $2\%$ agreement for $z\leq1$ and $k\leq 1~h/{\rm Mpc}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation., Comment: 20 pages, 4 figures, 4 tables

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents

Dark sirens are a powerful way to infer cosmological and astrophysical parameters from the combination of gravitational wave sirens and galaxy catalogues. Importantly, the method relies on the completeness of the galaxy catalogues being well modelled. A magnitude-limited catalogue will always be incomplete to some extent, requiring a completion scheme to avoid biasing the parameter inference. Standard methods include homogeneous and multiplicative completion, which have the advantage of simplicity but underestimate or overestimate the amplitude of structure at low completeness, respectively. In this work, we propose a new method to complete galaxy catalogues which uses clustering information to incorporate knowledge of the large scale structure into the dark sirens method. We find that if the structure of the true number of galaxies is sufficiently well preserved in the catalogue, our estimator can perform drastically better than both homogeneous and multiplicative completion. We lay the foundations for a maximally informative dark sirens analysis and discuss its limitations., Comment: 29 pages, 11 figures, adapted to match published version in JCAP

We test the convergence of fast simulations based on the COmoving Lagrangian Acceleration (COLA) method for predictions of the matter power spectrum, specialising our analysis in the redshift range $1 \le z \le 1.65$, relevant to high-redshift spectroscopic galaxy surveys. We then focus on the enhancement of the matter power spectrum in modified gravity (MG), the boost factor, using the Dvali-Gabadadze-Porrati (DGP) theory as a test case but developing a general approach that can be applied to other MG theories. After identifying the minimal simulation requirements for accurate DGP boost factors, we design and produce a COLA simulation suite that we use to train a neural network emulator for the DGP boost factor. Using MG-AREPO simulations as a reference, we estimate the emulator accuracy to be of $\sim 3\%$ up to $k=5 \, h {\rm Mpc}^{-1}$ at $0 \leq z \leq 2$. We make the emulator publicly available at: https://github.com/BartolomeoF/nDGPemu., Comment: 24 pages, 17 figures, Comments are welcome - Accepted for publication in JCAP

Third-generation gravitational wave (GW) observatories such as the Einstein Telescope and Cosmic Explorer, together with the LSST survey at the Vera Rubin Observatory, will yield an abundance of extra-galactic transient objects. This opens the exciting possibility of using GW sources and Supernovae Type Ia (SNIa) as luminosity distance tracers of large-scale structure for the first time. The large volumes accessible to these surveys imply that we may need to include relativistic corrections, such as lensing and Doppler magnification. However, the amplitude of these effects depends on the magnification and evolution biases of the transient sources, which are not yet understood. In this paper we develop comprehensive frameworks to address and model these biases for both populations of transient objects; in particular, we define how to compute these biases for GW sources. We then analyse the impact of magnification and evolution biases on the relativistic corrections and on the angular power spectrum of these sources. We show that correct modelling and implementation of these biases is crucial for measuring the cross-correlations of transient sources at higher redshifts., Comment: 28 pages, 14 figures, 4 tables, comments welcome

The dark sirens method enables us to use gravitational wave events without electromagnetic counterparts as tools for cosmology and tests of gravity. Furthermore, the dark sirens analysis code gwcosmo can now robustly account for information coming from both galaxy catalogues and the compact object mass distribution. We present here an extension of the gwcosmo code and methodology to constrain parameterized deviations from General Relativity that affect the propagation of gravitational waves. We show results of our analysis using data from the GWTC-3 gravitational wave catalogues, in preparation for application to the O4 observing run. After testing our pipelines using the First Two Years mock data set, we reanalyse 46 events from GWTC-3, and combine the posterior for BBH and NSBH sampling results for the first time. We obtain joint constraints on H0 and parameterized deviations from General Relativity in the Power Law + Peak BBH population model. With increased galaxy catalogue support in the future, our work sets the stage for dark sirens to become a powerful tool for testing gravity.

In the absence of numerous gravitational-wave detections with confirmed electromagnetic counterparts, the "dark siren" method has emerged as a leading technique of gravitational-wave cosmology. The method allows redshift information of such events to be inferred statistically from a catalogue of potential host galaxies. Due to selection effects, dark siren analyses necessarily depend on the mass distribution of compact objects and the evolution of their merger rate with redshift. Informative priors on these quantities will impact the inferred posterior constraints on the Hubble constant ($H_0$). It is thus crucial to vary these unknown distributions during an $H_0$ inference. This was not possible in earlier analyses due to the high computational cost, restricting them to either excluding galaxy catalogue information, or fixing the gravitational-wave population mass distribution and risking introducing bias to the $H_0$ measurement. This paper introduces a significantly enhanced version of the Python package GWCOSMO, which allows joint estimation of cosmological and compact binary population parameters. This thereby ensures the analysis is now robust to a major source of potential bias. The gravitational-wave events from the Third Gravitational-Wave Transient Catalogue are reanalysed with the GLADE+ galaxy catalogue, and an updated, more reliable measurement of $H_0=69^{+12}_{-7}$ km s$^{-1}$ Mpc$^{-1}$ is found (maximum a posteriori probability and 68% highest density interval). This improved method will enable cosmological analyses with future gravitational-wave detections to make full use of the information available (both from galaxy catalogues and the compact binary population itself), leading to promising new independent bounds on the Hubble constant., Comment: 31 pages, 11 figures, accepted for publication in JCAP

Next generation surveys will provide us with an unprecedented number of detections of supernovae Type Ia and gravitational wave merger events. Cross-correlations of such objects offer novel and powerful insights into the large-scale distribution of matter in the universe. Both of these sources carry information on their luminosity distance, but remain uninformative about their redshifts; hence their clustering analyses and cross-correlations need to be carried out in luminosity distance space, as opposed to redshift space. In this paper, we calculate the full expression for the number count fluctuation in terms of a perturbation to the observed luminosity distance. We find the expression to differ significantly from the one commonly used in redshift space. Furthermore, we present a comparison of the number count angular power spectra between luminosity distance and redshift spaces. We see a wide divergence between the two at large scales, and we note that lensing is the main contribution to such differences. On such scales and at higher redshifts the difference between the angular power spectra in luminosity distance and redshift spaces can be roughly 50$\%$. We also investigate cross-correlating different redshift bins using different tracers, i.e. one in luminosity distance space and one in redshift, simulating the cross-correlation angular power spectrum between background gravitational waves/supernovae and foreground galaxies. Finally, we show that in a cosmic variance limited survey, the relativistic corrections to the density-only term ought to be included., Comment: 34 pages, 9 figures, 1 table. Agrees with published version

Two major challenges of contemporary cosmology are the Hubble tension and the cosmic dipole tension. At the crossroad of these, we investigate the impact of peculiar velocities on estimations of the Hubble constant from time-delay cosmography. We quantify the bias on the inference of the Hubble constant due to peculiar velocities of the lens, the source and of the observer. The former two, which may cancel from one system to another, affect the determination of the angular diameter distances in the time-delay formula, and reconstructed quantities like the angle to the source, via a lens model. On the other hand, the peculiar velocity of the observer, which is a debated quantity in the context of the cosmic dipole tension, systematically affects observed angles through aberration, redshifts, angular diameter distance and reconstructed quantities. We compute in detail the effect of these peculiar velocities on the inference of the Hubble constant to linear order in the peculiar velocities for the seven lenses of the H0LiCOW/TDCOSMO collaboration. The bias generated by the observer's peculiar velocity alone can reach $1.15\%$ for the lenses which are well aligned with it. This results in a $0.25 \%$ bias for the seven combined lenses. Assuming a typical peculiar velocity of $300$ km s$^{-1}$ for the lens and the source galaxies, these add an additional random uncertainty, which can reach $1\%$ for an individual lens, but reduces to $0.24\%$ for the full TDCOSMO sample. The picture may change if peculiar velocities turn out to be larger than expected. Any time-delay cosmography program which aims for percent precision on the Hubble constant may need to take this source of systematic bias into account. This is especially so for future ground-based surveys which cover a fraction of the celestial sphere that is well aligned with the observer's peculiar velocity., Comment: 19 pages, 8 figures, 2 tables, online version adapted to make it consistent with published version

We outline the ``dark siren'' galaxy catalog method for cosmological inference using gravitational wave (GW) standard sirens, clarifying some common misconceptions in the implementation of this method. When a confident transient electromagnetic counterpart to a GW event is unavailable, the identification of a unique host galaxy is in general challenging. Instead, as originally proposed by Schutz (1986), one can consult a galaxy catalog and implement a dark siren statistical approach incorporating all potential host galaxies within the localization volume. Trott & Hunterer 2021 recently claimed that this approach results in a biased estimate of the Hubble constant, $H_0$, when implemented on mock data, even if optimistic assumptions are made. We demonstrate explicitly that, as previously shown by multiple independent groups, the dark siren statistical method leads to an unbiased posterior when the method is applied to the data correctly. We highlight common sources of error possible to make in the generation of mock data and implementation of the statistical framework, including the mismodeling of selection effects and inconsistent implementations of the Bayesian framework, which can lead to a spurious bias., Comment: 19 pages, 11 figures. Submitted to ApJ

Effective field theories (EFT) of dark energy (DE) -- built to parameterise the properties of DE in an agnostic manner -- are severely constrained by measurements of the propagation speed of gravitational waves (GW). However, GW frequencies probed by ground-based interferometers lie around the typical strong coupling scale of the EFT, and it is likely that the effective description breaks down before even reaching that scale. We discuss how this leaves the possibility that an appropriate ultraviolet completion of DE scenarios, valid at scales beyond an EFT description, can avoid present constraints on the GW speed. Instead, additional constraints in the lower frequency LISA band would be harder to escape, since the energies involved are orders of magnitude lower. By implementing a method based on GW multiband detections, we show indeed that a single joint observation of a GW150914-like event by LISA and a terrestrial interferometer would allow one to constrain the speed of light and gravitons to match to within $10^{-15}$. Multiband GW observations can therefore firmly constrain scenarios based on the EFT of DE, in a robust and unambiguous way., Comment: 11 pages, 4 figures. Updated to match the version published in JCAP

We introduce $\texttt{Hi-COLA}$, a code designed to run fast, approximate $\textit{N}$-body simulations of non-linear structure formation in reduced Horndeski gravity. Given an input Lagrangian, $\texttt{Hi-COLA}$ dynamically constructs the appropriate field equations and consistently solves for the cosmological background, linear growth, and screened fifth force of that theory. Hence $\texttt{Hi-COLA}$ is a general, adaptable, and useful tool that allows the mildly non-linear regime of many Horndeski theories to be investigated for the first time, at low computational cost. In this work, we first describe the screening approximations and simulation setup of $\texttt{Hi-COLA}$ for theories with Vainshtein screening. We validate the code against traditional $\textit{N}$-body simulations for cubic Galileon gravity, finding $2.5\%$ agreement up to $k_{\rm max}=1.2~h/{\rm Mpc}$. To demonstrate the flexibility of $\texttt{Hi-COLA}$, we additionally run the first simulations of an extended shift-symmetric gravity theory. We use the consistency and modularity of $\texttt{Hi-COLA}$ to dissect how the modified background, linear growth, and screened fifth force all contribute to departures from $\Lambda$CDM in the non-linear matter power spectrum. $\texttt{Hi-COLA}$ can be found at https://github.com/Hi-COLACode/Hi-COLA ., Comment: 51 pages, 20 figures, 5 tables, comments welcome. Updated to match the version published in JCAP

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe.

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\"atze 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., Comment: 38+20 pages, 22 figures. V2: Analysis improved, to be published in JCAP

Since 2015 the gravitational-wave observations of LIGO and Virgo have transformed our understanding of compact-object binaries. In the years to come, ground-based gravitational-wave observatories such as LIGO, Virgo, and their successors will increase in sensitivity, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will provide gravitational-wave observations of massive black holes binaries. Between the $\sim 10$-$10^3~\mathrm{Hz}$ band of ground-based observatories and the $\sim10^{-4}$-$10^{-1}~\mathrm{Hz}$ band of LISA lies the uncharted decihertz gravitational-wave band. We propose a Decihertz Observatory to study this frequency range, and to complement observations made by other detectors. Decihertz observatories are well suited to observation of intermediate-mass ($\sim10^2$-$10^4 M_\odot$) black holes; they will be able to detect stellar-mass binaries days to years before they merge, providing early warning of nearby binary neutron star mergers and measurements of the eccentricity of binary black holes, and they will enable new tests of general relativity and the Standard Model of particle physics. Here we summarise how a Decihertz Observatory could provide unique insights into how black holes form and evolve across cosmic time, improve prospects for both multimessenger astronomy and multiband gravitational-wave astronomy, and enable new probes of gravity, particle physics and cosmology., Comment: 13 pages, 1 figure. Published in Experimental Astronomy. Summarising white paper arXiv:1908.11375

The first multi-messenger gravitational wave event has had a transformative effect on the space of modified gravity models. In this paper we study the enhanced tests of gravity that are possible with a future set of gravitational wave standard siren events. We perform MCMC constraint forecasts for parameters in Horndeski scalar-tensor theories. In particular, we focus on the complementarity of gravitational waves with electromagnetic large-scale structure data from galaxy surveys. We find that the addition of fifty low redshift ($z \lesssim 0.2$) standard sirens from the advanced LIGO network offers only a modest improvement (a factor 1.1 -- 1.3, where 1.0 is no improvement) over existing constraints from electromagnetic observations of large-scale structures. In contrast, high redshift (up to $z \sim 10$) standard sirens from the future LISA satellite will improve constraints on the time evolution of the Planck mass in Horndeski theories by a factor $\sim 5$. By simulating different scenarios, we find this improvement to be robust to marginalisation over unknown merger inclination angles and to variation between three plausible models for the merger source population., Comment: 29+10 pages, 14 figures, 3 tables. Comments welcome. Main concept in Figure 1, main result in Figure 9. Accepted for publication in JCAP

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model., Comment: 26 pages, 2 figures. White paper submitted to ESA's Voyage 2050 call on behalf of the LISA Consortium 2050 Task Force

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology., Comment: 52 pages, 5 figures, 4 tables. Submitted to Classical & Quantum Gravity. Based upon a white paper for ESA's Voyage 2050 on behalf of the LISA Consortium 2050 Task Force

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., Comment: 108 pages; matches Reviews of Modern Physics accepted version; project website at https://www.novelprobes.org/

One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($\Lambda$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $\Lambda$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value., Comment: 61 pages. Some acknowledgments and references added. This is version-1.1 of an internal collaboration document of LSST-DESC that is being made public and is not planned for submission to a journal

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population., Comment: Science White Paper for the Astro2020 Decadal Survey, 5 pages text, 14 pages with references and affiliations

We investigate the potential of weak lensing by voids to test for deviations from General Relativity. We calculate the expected lensing signal of a scalar field with derivative couplings, finding that it has the potential to boost the tangential shear both within and outside the void radius. We use voids traced by Luminous Red Galaxies in SDSS to demonstrate the methodology of testing these predictions. We find that the void central density parameter, as inferred from the lensing signal, can shift from its GR value by up to 20% in some galileon gravity models. Since this parameter can be estimated independently using the galaxy tracer profiles of voids, our method provides a consistency check of the gravity theory. Although galileon gravity is now disfavoured as a source of cosmic acceleration by other datasets, the methods we demonstrate here can be used to test for more general fifth force effects with upcoming void lensing data., Comment: Small updates to figures 1, 2 and 6; typos fixed. Twelve pages plus references; seven figures

We analyse cosmological perturbations around a homogeneous and isotropic background for scalar-tensor, vector-tensor and bimetric theories of gravity. Building on previous results, we propose a unified view of the effective parameters of all these theories. Based on this structure, we explore the viable space of parameters for each family of models by imposing the absence of ghosts and gradient instabilities. We then focus on the quasistatic regime and confirm that all these theories can be approximated by the phenomenological two-parameter model described by an effective Newton's constant and the gravitational slip. Within the quasistatic regime we pinpoint signatures which can distinguish between the broad classes of models (scalar-tensor, vector-tensor or bimetric). Finally, we present the equations of motion for our unified approach in such a way that they can be implemented in Einstein-Boltzmann solvers.

The detection of an electromagnetic counterpart (GRB 170817A) to the gravitational wave signal (GW170817) from the merger of two neutron stars opens a completely new arena for testing theories of gravity. We show that this measurement allows us to place stringent constraints on general scalar-tensor and vector-tensor theories, while allowing us to place an independent bound on the graviton mass in bimetric theories of gravity. These constraints severely reduce the viable range of cosmological models that have been proposed as alternatives to general relativistic cosmology., Comment: 4 pages

We investigate the potential of high-energy astrophysical events, from which both massless and massive signals are detected, to probe fundamental physics. In particular, we consider how strong gravitational lensing can induce time delays in multi-messenger signals originating from the same source. Obvious messenger examples are massless photons and gravitational waves, and massive neutrinos, although more exotic applications can also be imagined, such as to massive gravitons or axions. The different propagation times of the massive and massless particles can, in principle, place bounds on the total neutrino mass and probe cosmological parameters. Whilst measuring such an effect may pose a significant experimental challenge, we believe that the `massive time delay' represents an unexplored fundamental physics phenomenon., Comment: 18+9 pages; updated to match published version, minor textual changes only

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission., Comment: This article provides an update of arXiv:1206.1225, with different authors. Forecasts are not updated in this version

We present a method for parametrizing linear cosmological perturbations of theories of gravity, around homogeneous and isotropic backgrounds. The method is sufficiently general and systematic that it can be applied to theories with any degrees of freedom (DoFs) and arbitrary gauge symmetries. In this paper, we focus on scalar-tensor and vector-tensor theories, invariant under linear coordinate transformations. In the case of scalar-tensor theories, we use our framework to recover the simple parametrizations of linearized Horndeski and "Beyond Horndeski" theories, and also find higher-derivative corrections. In the case of vector-tensor theories, we first construct the most general quadratic action for perturbations that leads to second-order equations of motion, which propagates two scalar DoFs. Then we specialize to the case in which the vector field is time-like (\`a la Einstein-Aether gravity), where the theory only propagates one scalar DoF. As a result, we identify the complete forms of the quadratic actions for perturbations, and the number of free parameters that need to be defined, to cosmologically characterize these two broad classes of theories., Comment: Two code packages (xIST and COPPER) are published with this paper and can be downloaded from https://github.com/noller/xIST. Final version with minor changes

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled., Comment: 99 pages, 8 figures. Version published in Physics of the Dark Universe

Extremely large surveys with future experiments like Euclid and the SKA will soon allow us to access perturbation modes close to the Hubble scale, with wavenumbers $k \sim \mathcal{H}$. If a modified gravity theory is responsible for cosmic acceleration, the Hubble scale is a natural regime for deviations from General Relativity (GR) to become manifest. The majority of studies to date have concentrated on the consequences of alternative gravity theories for the subhorizon, quasi-static regime, however. In this paper we investigate how modifications to the gravitational field equations affect perturbations around the Hubble scale. We choose functional forms to represent the generic scale-dependent behaviour of gravity theories that modify GR at long wavelengths, and study the resulting deviations of ultra large-scale relativistic observables from their GR behaviour. We find that these are small unless modifications to the field equations are drastic. The angular dependence and redshift evolution of the deviations is highly parameterisation- and survey-dependent, however, and so they are possibly a rich source of modified gravity phenomenology if they can be measured., Comment: Update to published version: fixes transposed panels in Figs. 9 and 10

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime., Comment: 188 pages, 46 figures, 6 tables, 903 references. Matches version published in Classical and Quantum Gravity. Supplementary data files available at http://www.phy.olemiss.edu/~berti/research/ and http://centra.tecnico.ulisboa.pt/network/grit/files/

By considering linear-order departures from general relativity, we compute a novel expression for the weak lensing convergence power spectrum under alternative theories of gravity. This comprises an integral over a 'kernel' of general relativistic quantities multiplied by a theory-dependent 'source' term. The clear separation between theory-independent and -dependent terms allows for an explicit understanding of each physical effect introduced by altering the theory of gravity. We take advantage of this to explore the degeneracies between gravitational parameters in weak lensing observations., Comment: 17 pages, 7 figures. v2: Minor changes to match version accepted by PRD

The current effort to test General Relativity employs multiple disparate formalisms for different observables, obscuring the relations between laboratory, astrophysical and cosmological constraints. To remedy this situation, we develop a parameter space for comparing tests of gravity on all scales in the universe. In particular, we present new methods for linking cosmological large-scale structure, the Cosmic Microwave Background and gravitational waves with classic PPN tests of gravity. Diagrams of this gravitational parameter space reveal a noticeable untested regime. The untested window, which separates small-scale systems from the troubled cosmological regime, could potentially hide the onset of corrections to General Relativity., Comment: 15+6 pages; 2 figures. Figures available online at http://users.ox.ac.uk/~lady2729/Welcome.html . Updated to match published version -- minor textual changes only

In the quasistatic regime, generic modifications to gravity can give rise to novel scale-dependence of the gravitational field equations. Crucially, the detectability of the new scale-dependent terms hinges upon the existence of an effective mass scale or length scale at which corrections to General Relativity become relevant. Starting from only a few basic principles, we derive the general form of this scale-dependence. Our method recovers results previously known in the specific case of Horndeski gravity, but also shows that they are valid more generally, beyond the regime of scalar field theories. We forecast the constraints that upcoming experiments will place on the existence of a new fundamental mass scale or length scale in cosmology., Comment: 7+5 pages, submitted to PRD

The growth rate of the large-scale structure of the universe has been advocated as the observable par excellence for testing gravity on cosmological scales. By considering linear-order deviations from General Relativity, we show that corrections to the growth rate, f, can be expressed as an integral over a `source' term, weighted by a theory-independent `response kernel'. This leads to an efficient and accurate `plug-and-play' expression for generating growth rates in alternative gravity theories, bypassing lengthy theory-specific computations. We use this approach to explicitly show that f is sensitive to a degenerate combination of modified expansion and modified clustering effects. Hence the growth rate, when used in isolation, is not a straightforward diagnostic of modified gravity., Comment: Updated to match version published in PRD. 12 pages plus appendices. The busy reader may like to focus on Section 4

A unified framework for theories of modified gravity will be an essential tool for interpreting the forthcoming deluge of cosmological data. We present such a formalism, the Parameterized Post-Friedmann framework (PPF), which parameterizes the cosmological perturbation theory of a wide variety of modified gravity models. PPF is able to handle spin-0 degrees of freedom from new scalar, vector and tensor fields, meaning that it is not restricted to simple models based solely on cosmological scalar fields. A direct correspondence is maintained between the parameterization and the underlying space of theories, which allows us to build up a `dictionary' of modified gravity theories and their PPF correspondences. In this paper we describe the construction of the parameterization and demonstrate its use through a number of worked examples relevant to the current literature. We indicate how the formalism will be implemented numerically, so that the dictionary of modified gravity can be pitted against forthcoming observations., Comment: 24 pages, updated to match version published in PRD. Discussion of section 4 extended. Suggestions for the busy reader are given at the end of section 1

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission., Comment: 236 pages, minor edits to match the journal version 2013

Parameterized frameworks for modified gravity are potentially useful tools for model-independent tests of General Relativity on cosmological scales. The toy model of an Einstein-de Sitter (EdS) universe provides a safe testbed in which to improve our understanding of their behaviour. We implement a mathematically consistent parameterization at the level of the field equations, and use this to calculate the evolution of perturbations in an EdS scenario. Our parameterization explicitly allows for new scalar degrees of freedom, and we compare this to theories in which the only degrees of freedom come from the metric and ordinary matter. The impact on the Integrated Sachs-Wolfe effect and canonically-conserved superhorizon perturbations is considered., Comment: Updated to match published version

There are a number of approaches to testing General Relativity (GR) on linear scales using parameterized frameworks for modifying cosmological perturbation theory. It is sometimes assumed that the details of any given parameterization are unimportant if one uses it as a diagnostic for deviations from GR. In this brief report we argue that this is not necessarily so. First we show that adopting alternative combinations of modifications to the field equations significantly changes the constraints that one obtains. In addition, we show that using a parameterization with insufficient freedom significantly tightens the apparent theoretical constraints. Fundamentally we argue that it is almost never appropriate to consider modifications to the perturbed Einstein equations as being constraints on the effective gravitational constant, for example, in the same sense that solar system constraints are. The only consistent modifications are either those that grant near-total freedom, as in decomposition methods, or ones which map directly to a particular part of theory space.

There is a distinct possibility that current and future cosmological data can be used to constrain Einstein's theory of gravity on the very largest scales. To be able to do this in a model-independent way, it makes sense to work with a general parameterization of modified gravity. Such an approach would be analogous to the Parameterized Post-Newtonian (PPN) approach which is used on the scale of the Solar System. A few such parameterizations have been proposed and preliminary constraints have been obtained. We show that the majority of such parameterizations are only exactly applicable in the quasistatic regime. On larger scales they fail to encapsulate the full behaviour of typical models currently under consideration. We suggest that it may be possible to capture the additions to the `Parameterized Post-Friedmann' (PPF) formalism by treating them akin to fluid perturbations., Comment: 24 pages. Minor changes to match the published version - references added, comments added to sections 1 and 5