The Einstein Telescope (ET), a proposed European ground-based gravitational-wave detector of third-generation, is an evolution of second-generation detectors such as Advanced LIGO, Advanced Virgo, and KAGRA which could be operating in the mid 2030s. ET will explore the universe with gravitational waves up to cosmological distances. We discuss its main scientific objectives and its potential for discoveries in astrophysics, cosmology and fundamental physics.

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.

Recent work has shown that modified gravitational wave (GW) propagation can be a powerful probe of dark energy and modified gravity, specific to GW observations. We use the technique of Gaussian processes, that allows the reconstruction of a function from the data without assuming any parametrization, to measurements of the GW luminosity distance from simulated joint GW-GRB detections, combined with measurements of the electromagnetic luminosity distance by simulated DES data. For the GW events we consider both a second-generation LIGO/Virgo/Kagra (HVLKI) network, and a third-generation detector such as the Einstein Telescope. We find that the HVLKI network at target sensitivity, with $O(15)$ neutron star binaries with electromagnetic counterpart, could already detect deviations from GR at a level predicted by some modified gravity models, and a third-generation detector such as ET would have a remarkable discovery potential. We discuss the complementarity of the Gaussian processes technique to the $(\Xi_0,n)$ parametrization of modified GW propagation.

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.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

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.

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.

We study the cosmological evolution of some nonlocal gravity models, when the initial conditions are set during a phase of primordial inflation. We examine in particular three models, the so-called RT, RR and $\Delta_4$ models, previously introduced by our group. We find that the RR and $\Delta_4$ models have a stable evolution also during inflation. The RT model has an apparent instability, but we show that, because of the smallness of the scale associated to the nonlocal term compared to the inflationary scale, this instability is innocuous and also the RT model has a viable evolution even when its initial conditions are set during a phase of primordial inflation.

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.

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.

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.

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.

We study two nonlinear extensions of the nonlocal $R\,\Box^{-2}R$ gravity theory. We extend this theory in two different ways suggested by conformal symmetry, either replacing $\Box^{-2}$ with $(-\Box + R/6)^{-2}$, which is the operator that enters the action for a conformally-coupled scalar field, or replacing $\Box^{-2}$ with the inverse of the Paneitz operator, which is a four-derivative operator that enters in the effective action induced by the conformal anomaly. We show that the former modification gives an interesting and viable cosmological model, with a dark energy equation of state today $w_{\rm DE}\simeq -1.01$, which very closely mimics $\Lambda$CDM and evolves asymptotically into a de Sitter solution. The model based on the Paneitz operator seems instead excluded by the comparison with observations. We also review some issues about the causality of nonlocal theories, and we point out that these nonlocal models can be modified so to nicely interpolate between Starobinski inflation in the primordial universe and accelerated expansion in the recent epoch.

Recent work has shown that modifications of General Relativity based on the addition to the action of a non-local term $R\,\Box^{-2}R$, or on the addition to the equations of motion of a term involving $(g_{\mu\nu}\Box^{-1} R)$, produce dynamical models of dark energy which are cosmologically viable both at the background level and at the level of cosmological perturbations. We explore a more general class of models based on the addition to the action of terms proportional to $R_{\mu\nu}\,\Box^{-2}R^{\mu\nu}$ and $C_{\mu\nu\rho\sigma}\, \Box^{-2}C^{\mu\nu\rho\sigma}$, where $C_{\mu\nu\rho\sigma}$ is the Weyl tensor. We find that the term $R_{\mu\nu}\,\Box^{-2}R^{\mu\nu}$ does not give a viable background evolution. The non-local Weyl-square term, in contrast, does not contribute to the background evolution but we find that, at the level of cosmological perturbations, it gives instabilities in the tensor sector. Thus, only non-local terms which depend just on the Ricci scalar $R$ appear to be cosmologically viable. We discuss how these results can provide a hint for the mechanism that might generate these effective non-local terms from a fundamental local theory.

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.

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.

We perform a detailed study of the cosmological dynamics of a recently proposed infrared modification of the Einstein equations, based on the introduction of a non-local term constructed with $m^2g_{\mu\nu}\Box^{-1} R$, where $m$ is a mass parameter. The theory generates automatically a dynamical dark energy component, that can reproduce the observed value of the dark energy density without introducing a cosmological constant. Fixing $m$ so to reproduce the observed value $\Omega_{\rm DE}\simeq 0.68$, and writing $w(a)=w_0+(1-a) w_a$, the model provides a neat prediction for the equation of state parameters of dark energy, $w_0\simeq -1.042$ and $w_a\simeq -0.020$. We show that, because of some freedom in the definition of $\Box^{-1}$, one can extend the construction so to define a more general family of non-local models. However, in a first approximation this turns out to be equivalent to adding an explicit cosmological constant term on top of the dynamical dark energy component. This leads to an extended model with two parameters, $\Omega_{\Lambda}$ and $m$. Even in this extension the EOS parameter $w_0$ is always on the phantom side, in the range $-1.33 < w_0\leq -1$, and there is a prediction for the relation between $w_0$ and $w_a$.

Recent work has shown that non-local modifications of the Einstein equations can have interesting cosmological consequences and can provide a dynamical origin for dark energy, consistent with existing data. At first sight these theories are plagued by ghosts. We show that these apparent ghost-like instabilities do not describe actual propagating degrees of freedom, and there is no issue of ghost-induced quantum vacuum decay.

We study the possibility of using matter wave interferometry techniques to build a gravitational wave detector. We derive the response function and find that it contains a term proportional to the derivative of the gravitational wave, a point which has been disputed recently. We then study in detail the sensitivity that can be reached by such a detector and find that, if it is operated near resonance, it can reach potentially interesting values in the high frequency regime. The correlation between two or more of such devices can further improve the sensitivity for a stochastic signal.

Gravitational wave bursters are sources which emit repeatedly bursts of gravitational waves, and have been recently suggested as potentially interesting candidates for gravitational wave (GW) detectors. Mechanisms that could give rise to a GW burster can be found for instance in highly magnetized neutron stars (the magnetars which explain the phenomenon of soft gamma repeaters), in accreting neutron stars and in hybrid stars with a quark core. We point out that these sources have very distinctive experimental signatures. In particular, as already observed in the gamma-ray bursts from soft gamma repeaters, the energy spectrum of the events is a power-law, dN\sim E^{-\gamma}dE with \gamma\simeq 1.6, and they have a distribution of waiting times (the times between one outburst and the next) significantly different from the distribution of uncorrelated events. We discuss possible detection strategies that could be used to search for these events in existing gravitational wave detectors.

Quantum gravity "foam", among its various generic Lorentz non-invariant effects, would cause neutrino mixing. It is shown here that, if the foam is manifested as a nonrenormalizable effect at scale M, the oscillation length generically decreases with energy $E$ as (E/M)^(-2). Neutrino observatories and long-baseline experiments should have therefore already observed foam-induced oscillations, even if M is as high as the Planck energy scale. The null results, which can be further strengthened by better analysis of current data and future experiments, can be taken as experimental evidence that Lorentz invariance is fully preserved at the Planck scale, as is the case in critical string theory.