Search

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