Search

The next generation of ground-based gravitational-wave detectors will observe coalescences of black holes and neutron stars throughout the cosmos, thousands of them with exceptional fidelity. The Science Book is the result of a 3-year effort to study the science capabilities of networks of next generation detectors. Such networks would make it possible to address unsolved problems in numerous areas of physics and astronomy, from Cosmology to Beyond the Standard Model of particle physics, and how they could provide insights into workings of strongly gravitating systems, astrophysics of compact objects and the nature of dense matter. It is inevitable that observatories of such depth and finesse will make new discoveries inaccessible to other windows of observation. In addition to laying out the rich science potential of the next generation of detectors, this report provides specific science targets in five different areas in physics and astronomy and the sensitivity requirements to accomplish those science goals. This report is the second in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book (this report), iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network., Comment: 69 pages, 18 figures

We provide a systematic and updated discussion of a research line carried out by our group over the last few years, in which gravity is modified at cosmological distances by the introduction of nonlocal terms, assumed to emerge at an effective level from the infrared behavior of the quantum theory. The requirement of producing a viable cosmology turns out to be very stringent and basically selects a unique model, in which the nonlocal term describes an effective mass for the conformal mode. We discuss how such a specific structure could emerge from a fundamental local theory of gravity, and we perform a detailed comparison of this model with the most recent cosmological datasets, confirming that it fits current data at the same level as $\Lambda$CDM. Most notably, the model has striking predictions in the sector of tensor perturbations, leading to a very large effect in the propagation of gravitational wave (GWs) over cosmological distances. At the redshifts relevant for the next generation of GW detectors such as Einstein Telescope, Cosmic Explorer and LISA, this leads to deviations from GR that could be as large as $80\%$, and could be verified with the detection of just a single coalescing binary with electromagnetic counterpart. This would also have potentially important consequences for the search of the counterpart since, for a given luminosity distance to the source, as inferred through the GW signal, the actual source redshift could be significantly different from that predicted by $\Lambda$CDM. At the redshifts relevant for advanced LIGO/Virgo/Kagra the effect is smaller, but still potentially observable over a few years of runs at target sensitivity., Comment: 84 pages, 22 figures

We discuss a modified gravity model which fits cosmological observations at a level statistically indistinguishable from $\Lambda$CDM and at the same time predicts very large deviations from General Relativity (GR) in the propagation of gravitational waves (GWs) across cosmological distances. The model is a variant of the RT nonlocal model proposed and developed by our group, with initial conditions set during inflation, and predicts a GW luminosity distance that, at the redshifts accessible to LISA or to a third-generation GW detector such as the Einstein Telescope (ET), can differ from that in GR by as much as $60\%$. An effect of this size could be detected with just a single standard siren with counterpart by LISA or ET. At the redshifts accessible to a LIGO/Virgo/Kagra network at target sensitivity the effect is smaller but still potentially detectable. Indeed, for the recently announced LIGO/Virgo NS-BH candidate S190814bv, the RT model predicts that, given the measured GW luminosity distance, the actual luminosity distance, and the redshift of an electromagnetic counterpart, would be smaller by as much as $7\%$ with respect to the value inferred from $\Lambda$CDM., Comment: v2: added discussion of the effect of modified GW propagation on the recent LIGO/Virgo NS-BH candidate event S190814bv

Gravitational-wave (GW) detectors can contribute to the measurement of cosmological parameters and to testing the dark-energy sector of alternatives to $\Lambda$CDM, by using standard sirens. In this paper we focus on binary neutron stars with a counterpart detected through a gamma-ray burst (GRB), both at a second-generation network made by advanced LIGO+advanced Virgo+LIGO India+Kagra, and at third-generation (3G) detectors, discussing in particular the cases of a single Einstein Telescope (ET), and of a network of ET plus two Cosmic Explorer (CE). We construct mock catalogs of standard sirens, using different scenarios for the local merger rate and for the detection of the electromagnetic counterpart. For 3G detectors we estimate the coincidences with a GRB detector with the characteristics of the proposed THESEUS mission. We discuss how these standard sirens with a GRB counterpart can improve the determination of cosmological parameters (and particularly of $H_0$) in $\Lambda$CDM, and we then study how to extract information on dark energy, considering both a non-trivial dark energy equation of state and modified GW propagation. We find that a 2G detector network can already reach, over several years of data taking, an interesting sensitivity to modified GW propagation, while a single ET detector would have a remarkable potential for discovery. We also find that, to fully exploit the potential of a ET+CE+CE network, it is necessary a much stronger program of search for electromagnetic counterparts (or else to resort to statistical methods for standard sirens), and furthermore gravitational lensing can become a limiting factor., Comment: 50 pages, 28 figures

The background photon temperature $\bar T$ is one of the fundamental cosmological parameters. Despite its significance, $\bar T$ has never been allowed to vary in the data analysis, owing to the precise measurement of the comic microwave background (CMB) temperature by COBE FIRAS. However, even in future CMB experiments, $\bar T$ will remain unknown due to the unknown monopole contribution $\Theta_0$ at our position to the observed (angle-averaged) temperature $\langle T\rangle^{\rm obs}$. By fixing $\bar T\equiv\langle T\rangle^{\rm obs}$, the standard analysis underestimates the error bars on cosmological parameters, and the best-fit parameters obtained in the analysis are biased in proportion to the unknown amplitude of $\Theta_0$. Using the Fisher formalism, we find that these systematic errors are smaller than the error bars from the $Planck$ satellite. However, with $\bar T\equiv\langle T\rangle^{\rm obs}$, these systematic errors will always be present and irreducible, and future cosmological surveys might misinterpret the measurements., Comment: 5 pages, 2 figures, published in PRD. Planck 2018 analysis is discussed

This is the first in a series of papers devoted to fully general-relativistic $N$-body simulations applied to late-time cosmology. The purpose of this paper is to present the combination of a numerical relativity scheme, discretization method and time-integration algorithm that provides satisfyingly stable evolution. More precisely, we show that it is able to pass a robustness test and to follow scalar linear modes around an expanding homogeneous and isotropic space-time. Most importantly, it is able to evolve typical cosmological initial conditions on comoving scales down to tenths of megaparsecs with controlled constraint and energy-momentum conservation violations all the way down to the regime of strong inhomogeneity., Comment: 28 pages, 16 figures

This Astro-2020 White Paper deals with what we might learn from future gravitational wave observations about the early universe phase transitions and their energy scales, primordial black holes, Hubble parameter, dark matter and dark energy, modified theories of gravity and extra dimensions., Comment: 13 pages, 3 figures, White Paper submitted to the Astro-2020 (2020 Astronomy and Astrophysics Decadal Survey) by GWIC-3G Science Case Team (GWIC: Gravitational-Wave International Committee)

We study the cosmology of a specific class of nonlocal model of modified gravity, the so-called Deser-Woodard (DW) model, modifying the Einstein-Hilbert action by a term $\sim R f(\Box^{-1}R)$, where $f$ is a free function. Choosing $f$ so as to reproduce the $\Lambda{\rm CDM}$ cosmological background expansion history within the nonlocal model, we implement the model in a cosmological linear Einstein--Boltzmann solver and study the deviations to GR the model induces in the scalar and tensor perturbations. We observe that the DW nonlocal model describes a modified propagation for the gravitational waves, as well as a lower linear growth rate and a stronger lensing power as compared to $\Lambda{\rm CDM}$, up to several percents. Such prominent growth and lensing features lead to the inference of a significantly smaller value of $\sigma_8$ with respect to the one in $\Lambda{\rm CDM}$, given \textit{Planck} CMB+lensing data. The prediction for the linear growth rate $f \sigma_8$ within the DW model is therefore significantly smaller than the one in $\Lambda{\rm CDM}$ and the addition of growth rate data $f \sigma_8$ from Redshift-space distortion measurements to \textit{Planck} CMB+lensing, opens a (dominant) tension between Redshift-space distortion data and the reconstructed \textit{Planck} CMB lensing potential. However, model selection issues only result in "weak" evidences for $\Lambda{\rm CDM}$ against the DW model given the data. Such a fact shows that the datasets we consider are not constraining enough for distinguishing between the models. As we discuss, the addition of galaxy WL data or cosmological constraints from future galaxy clustering, weak lensing surveys, but also third generation gravitational wave interferometers, prove to be useful for discriminating modified gravity models such as the DW one from $\Lambda{\rm CDM}$, within the close future., Comment: 28 pages, 13 figures, v2: references added, typographical errors corrected

This paper contains a comprehensive comparison catalog of `Apples with Apples' tests for the BSSNOK, CCZ4 and Z4c numerical relativity schemes, with and without constraint damping terms for the latter two. We use basic numerical methods and reach the same level of accuracy as existing results in the literature. We find that the best behaving scheme is generically CCZ4 with constraint damping terms., Comment: 33 pages, 25 figures

Studies of dark energy at advanced gravitational-wave (GW) interferometers normally focus on the dark energy equation of state $w_{\rm DE}(z)$. However, modified gravity theories that predict a non-trivial dark energy equation of state generically also predict deviations from general relativity in the propagation of GWs across cosmological distances, even in theories where the speed of gravity is equal to $c$. We find that, in generic modified gravity models, the effect of modified GW propagation dominates over that of $w_{\rm DE}(z)$, making modified GW propagation a crucial observable for dark energy studies with standard sirens. We present a convenient parametrization of the effect in terms of two parameters $(\Xi_0,n)$, analogue to the $(w_0,w_a)$ parametrization of the dark energy equation of state, and we give a limit from the LIGO/Virgo measurement of $H_0$ with the neutron star binary GW170817. We then perform a Markov Chain Monte Carlo analysis to estimate the sensitivity of the Einstein Telescope (ET) to the cosmological parameters, including $(\Xi_0,n)$, both using only standard sirens, and combining them with other cosmological datasets. In particular, the Hubble parameter can be measured with an accuracy better than $1\%$ already using only standard sirens while, when combining ET with current CMB+BAO+SNe data, $\Xi_0$ can be measured to $0.8\%$ . We discuss the predictions for modified GW propagation of a specific nonlocal modification of gravity, recently developed by our group, and we show that they are within the reach of ET. Modified GW propagation also affects the GW transfer function, and therefore the tensor contribution to the ISW effect., Comment: 25 pages, 23 figures: v3: several significant improvements

In modified gravity the propagation of gravitational waves (GWs) is in general different from that in general relativity. As a result, the luminosity distance for GWs can differ from that for electromagnetic signals, and is affected both by the dark energy equation of state $w_{\rm DE}(z)$ and by a function $\delta(z)$ describing modified propagation. We show that the effect of modified propagation in general dominates over the effect of the dark energy equation of state, making it easier to distinguish a modified gravity model from $\Lambda$CDM. We illustrate this using a nonlocal modification of gravity, that has been shown to fit remarkably well CMB, SNe, BAO and structure formation data, and we discuss the prospects for distinguishing nonlocal gravity from $\Lambda$CDM with the Einstein Telescope. We find that, depending on the exact sensitivity, a few tens of standard sirens with measured redshift at $z\sim 0.4$, or a few hundreds at $1 < z < 2$, could suffice., Comment: 6 pages, 3 figures; v4: minor modifications; the version to appear in PRD

Even if the fundamental action of gravity is local, the corresponding quantum effective action, that includes the effect of quantum fluctuations, is a nonlocal object. These nonlocalities are well understood in the ultraviolet regime but much less in the infrared, where they could in principle give rise to important cosmological effects. Here we systematize and extend previous work of our group, in which it is assumed that a mass scale $\Lambda$ is dynamically generated in the infrared, giving rise to nonlocal terms in the quantum effective action of gravity. We give a detailed discussion of conceptual aspects related to nonlocal gravity and of the cosmological consequences of these models. The requirement of providing a viable cosmological evolution severely restricts the form of the nonlocal terms, and selects a model (the so-called RR model) that corresponds to a dynamical mass generation for the conformal mode. For such a model: (1) there is a FRW background evolution, where the nonlocal term acts as an effective dark energy with a phantom equation of state, providing accelerated expansion without a cosmological constant. (2) Cosmological perturbations are well behaved. (3) Implementing the model in a Boltzmann code and comparing with observations we find that the RR model fits the CMB, BAO, SNe, structure formation data and local $H_0$ measurements at a level statistically equivalent to $\Lambda$CDM. (4) Bayesian parameter estimation shows that the value of $H_0$ obtained in the RR model is higher than in $\Lambda$CDM, reducing to $2.0\sigma$ the tension with the value from local measurements. (5) The RR model provides a prediction for the sum of neutrino masses that falls within the limits set by oscillation and terrestrial experiments. (6) Gravitational waves propagate at the speed of light, complying with the limit from GW170817/GRB 170817A., Comment: 60 pages, 12 figures; v2: references added

Prior change is discussed in observational constraints studies of nonlocally modified gravity. In the latter, a model characterized by a modification of the form $\sim m^2 R\Box^{-2}R$ to the Einstein-Hilbert action was compared against the base $\Lambda$CDM one in a Bayesian way. It was found that the competing modified gravity model is significantly disfavored (at $22 \,$:$\, 1$ in terms of betting-odds) against $\Lambda$CDM given CMB+SNIa+BAO data, because of a dominant tension appearing in the $H_0 \,$-$\, \Omega_M$ plan. We identify the underlying mechanism generating such a tension and show that it is mostly caused by the late-time, quite smooth, phantom nature of the effective dark energy described by the nonlocal model. We find possible solutions for it to be resolved and explore a given one that consists in extending the initial baseline from one massive neutrino eigenstate to three degenerate ones, whose absolute mass $\sum m_\nu \, / \, 3$ is allowed to take values within a reasonable prior interval. As a net effect, the absolute neutrino mass is inferred to be non-vanishing at $2 \sigma$ level, best-fitting at $\sum m_\nu \approx 0.21 {\, \rm eV}$, and the Bayesian tension disappears rendering the nonlocal gravity model statistically equivalent to $\Lambda$CDM, given recent CMB+SNIa+BAO data. We also discuss constraints from growth rate measurements $f \sigma_8$ whose fit is found to be improved by a larger massive neutrino fraction as well. The $\nu$-extended nonlocal model also prefers a higher value of $H_0$ than $\Lambda$CDM, therefore in better agreement with local measurements., Comment: 15 pages, 4 figures, accepted in PRD

Cosmological simulations involving the fully covariant gravitational dynamics may prove relevant in understanding relativistic/non-linear features and, therefore, in taking better advantage of the upcoming large scale structure survey data. We propose a new 3+1 integration scheme for General Relativity in the case where the matter sector contains a minimally-coupled perfect fluid field. The original feature is that we completely eliminate the fluid components through the constraint equations, thus remaining with a set of unconstrained evolution equations for the rest of the fields. This procedure does not constrain the lapse function and shift vector, so it holds in arbitrary gauge and also works for arbitrary equation of state. An important advantage of this scheme is that it allows one to define and pass an adaptation of the robustness test to the cosmological context, at least in the case of pressureless perfect fluid matter, which is the relevant one for late-time cosmology., Comment: 14 pages, 3 figures, matches published version

We present a comprehensive and updated comparison with cosmological observations of two non-local modifications of gravity previously introduced by our group, the so called RR and RT models. We implement the background evolution and the cosmological perturbations of the models in a modified Boltzmann code, using CLASS. We then test the non-local models against the {\em Planck} 2015 TT, TE, EE and Cosmic Microwave Background (CMB) lensing data, isotropic and anisotropic Baryonic Acoustic Oscillations (BAO) data, JLA supernovae, $H_0$ measurements and growth rate data, and we perform Bayesian parameter estimation. We then compare the RR, RT and $\Lambda$CDM models, using the Savage-Dickey method. We find that the RT model and $\Lambda$CDM perform equally well, while the performance of the RR model with respect to $\Lambda$CDM depends on whether or not we include a prior on $H_0$ based on local measurements., Comment: 42 pages, 18 figures; v2: the version to appear in JCAP

We study the cosmological predictions of two recently proposed non-local modifications of General Relativity. Both models have the same number of parameters as $\Lambda$CDM, with a mass parameter $m$ replacing the cosmological constant. We implement the cosmological perturbations of the non-local models into a modification of the CLASS Boltzmann code, and we make a full comparison to CMB, BAO and supernova data. We find that the non-local models fit these datasets as well as $\Lambda$CDM. For both non-local models parameter estimation using Planck+JLA+BAO data gives a value of $H_0$ higher than in $\Lambda$CDM, and in better agreement with the values obtained from local measurements., Comment: 14 pages, 4 figures

We analyse the classical stability of the model proposed by Maggiore and Mancarella, where gravity is modified by a term $\sim m^2 R \Box^{-2} R$ to produce the late-time acceleration of the expansion of the universe. Our study takes into account all excitations of the metric that can potentially drive an instability. There are some subtleties in identifying these modes, as a non-local field theory contains dynamical fields which yet do not correspond to degrees of freedom. Since some of them are ghost-like, we clarify the impact of such modes on the stability of the solutions of interest that are the flat space-time and cosmological solutions. We then find that flat space-time is unstable under scalar perturbations, but the instability manifests itself only at cosmological scales, i.e. out of the region of validity of this solution. It is therefore the stability of the FLRW solution which is relevant there, in which case the scalar perturbations are known to be well-behaved by numerical studies. By finding the analytic solution for the late-time behaviour of the scale factor, which leads to a big rip singularity, we argue that the linear perturbations are bounded in the future because of the domination of Hubble friction. In particular, this effect damps the scalar ghost perturbations which were responsible for destabilizing Minkowski space-time. Thus, the model remains phenomenologically viable., Comment: 19 pages, 4 figures, technical details/references added, version accepted for publication in JCAP

We study the cosmological consequences of a recently proposed nonlocal modification of general relativity, obtained by adding a term $m^2R\,\Box^{-2}R$ to the Einstein-Hilbert action. The model has the same number of parameters as $\Lambda$CDM, with $m$ replacing $\Omega_{\Lambda}$, and is very predictive. At the background level, after fixing $m$ so as to reproduce the observed value of $\Omega_M$, we get a pure prediction for the equation of state of dark energy as a function of redshift, $w_{\rm DE}(z)$, with $w_{\rm DE}(0)$ in the range $[-1.165,-1.135]$ as $\Omega_M$ varies over the broad range $\Omega_M\in [0.20,0.36]$. We find that the cosmological perturbations are well-behaved, and the model fully fixes the dark energy perturbations as a function of redshift $z$ and wavenumber $k$. The nonlocal model provides a good fit to supernova data and predicts deviations from General Relativity in structure formation and in weak lensing at the level of 3-4%, therefore consistent with existing data but readily detectable by future surveys. For the logarithmic growth factor we obtain $\gamma\simeq 0.53$, to be compared with $\gamma\simeq 0.55$ in $\Lambda$CDM. For the Newtonian potential on subhorizon scales our results are well fitted by $\Psi(a;k)=[1+\mu_s a^s]\Psi_{\rm GR}(a;k)$ with a scale-independent $\mu_s\simeq 0.09$ and $s\simeq 2$, while the anisotropic stress is negligibly small., Comment: 39 pages, 20 figures