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 impact of the limit on $|\dot{G}|/G$ from Lunar Laser Ranging on "nonlocal gravity", i.e. on models of the quantum effective action of gravity that include nonlocal terms relevant in the infrared, such as the "RR" and "RT" models proposed by our group, and the Deser-Woodard (DW) model. We elaborate on the analysis of Barreira et al. [1] and we confirm their findings that (under plausible assumptions such as the absence of strong backreaction from non-linear structures), the RR model is ruled out. We also show that the mechanism of "perfect screening for free" suggested for the DW model actually does not work and the DW model is also ruled out. In contrast, the RT model passes all phenomenological consistency tests and is still a viable candidate.

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.

Nonlocal gravity has been shown to provide a phenomenologically viable infrared modification of GR. A natural question is whether the required nonlocality can emerge from perturbative quantum loop corrections due to light particles. We show that this is not the case. For the value of the mass scale of the non-local models required by cosmology, the perturbative form factors obtained from the loop corrections, in the present cosmological epoch, are in the regime where they are local. The mechanism behind the generation of the required nonlocality must be more complex, possibly related to strong infrared effects and non-perturbative mass generation for the conformal mode.

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 discuss static spherically symmetric solutions in a recently proposed non-local infrared modification of Einstein equations induced by a term $m^2g_{\mu\nu}\Box^{-1} R$, where $m$ is a mass scale. We find that, contrary to what happens in usual theories of massive gravity, in this non-local theory there is no vDVZ discontinuity and classical non-linearities do not become large below a Vainshtein radius parametrically larger than the Schwarzschild radius $r_S$. Rather on the contrary, in the regime $r\ll m^{-1}$ the corrections to the metric generated by a static body in GR are of the form $1+{\cal O}(m^2r^2)$ and become smaller and smaller toward smaller values of $r$. The modification to the GR solutions only show up at $r > m^{-1}$. For $m={\cal O}(H_0)$, as required for having interesting cosmological consequences, the non-local theory therefore recovers all successes of GR at the solar system and lab scales.

We discuss a nonlocal modification of gravity obtained adding a term $m^2 R\,\Box^{-2}R$ to the Einstein-Hilbert action. We find that the mass parameter $m$ only affects the non-radiative sector of the theory, while the graviton remains massless, there is no propagating ghost-like degree of freedom, no vDVZ discontinuity, and no Vainshtein radius below which the theory becomes strongly coupled. For $m={\cal O}(H_0)$ the theory therefore recovers all successes of GR at solar system and lab scales, and only deviates from it at cosmological scales. We examine the cosmological consequences of the model and we find that it automatically generates a dynamical dark energy and a self-accelerating evolution. After fixing our only free parameter $m$ so to reproduce the observed value of the dark energy density today, we get a pure prediction for the dark energy equation of state, $w_ {\rm DE}\simeq-1.14$. This value is consistent with the existing data, and could also resolve the possible tension between the Planck data and local measurements of the Hubble parameter.

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 discuss the cosmological consequences of a model based on a non-local infrared modification of Einstein equations. We find that the model generates a dynamical dark energy that can account for the presently observed value of $\Omega_{\rm DE}$, without introducing a cosmological constant. Tuning a free mass parameter $m$ to a value $m\simeq 0.67 H_0$ we reproduce the observed value $\Omega_{\rm DE}\simeq 0.68$. This leaves us with no free parameter and we then get a pure prediction for the EOS parameter of dark energy. Writing $w_{\rm DE}(a)=w_0+(1-a) w_a$, we find $w_0\simeq-1.04$ and $w_a\simeq -0.02$, consistent with the Planck data, and on the phantom side. We also argue that non-local equations of the type that we propose must be understood as purely classical effective equations, such as those derived in semiclassical gravity for the in-in matrix elements of the metric. As such, any apparent ghost instability in such equations only affects the classical dynamics, but there is no propagating degree of freedom associated to the ghost, and no issue of ghost-induced quantum vacuum decay.

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

We construct a fully covariant theory of massive gravity which does not require the introduction of an external reference metric, and overcomes the usual problems of massive gravity theories (fatal ghosts instabilities, acausality and/or vDVZ discontinuity). The equations of motion of the theory are non-local, but respect causality. The starting point is the quadratic action proposed in the context of the degravitation idea. We show that it is possible to extended it to a fully non-linear covariant theory. This theory describes the five degrees of freedom of a massive graviton plus a scalar ghost. However, contrary to generic non-linear extensions of Fierz-Pauli massive gravity, the ghost has the same mass m as the massive graviton, independently of the background, and smoothly goes into a non-radiative degree of freedom for m-> 0. As a consequence, for $m\sim H_0$ the vacuum instability induced by the ghost is irrelevant even over cosmological time-scales. We finally show that an extension of the model degravitates a vacuum energy density of order $M_{Planck}^4$ down to a value of order $M_{Planck}^2 m^2$, which for $m \sim H_0$ is of order of the observed value of the vacuum energy density.

We re-examine the classic problem of the renormalization of zero-point quantum fluctuations in a Friedmann-Robertson-Walker background. We discuss a number of issues that arise when regularizing the theory with a momentum-space cutoff, and show explicitly how introducing non-covariant counter-terms allows to obtain covariant results for the renormalized vacuum energy-momentum tensor. We clarify some confusion in the literature concerning the equation of state of vacuum fluctuations. Further, we point out that the general structure of the effective action becomes richer if the theory contains a scalar field phi with mass m smaller than the Hubble parameter H(t). Such an ultra-light particle cannot be integrated out completely to get the effective action. Apart from the volume term and the Einstein-Hilbert term, that are reabsorbed into renormalizations of the cosmological constant and Newton's constant, the effective action in general also has a term proportional to F(phi)R, for some function F(phi). As a result, vacuum fluctuations of ultra-light scalar fields naturally lead to models where the dark energy density has the form rho_{DE}(t)=rho_X(t)+rho_Z(t), where rho_X is the component that accelerates the Hubble expansion at late times and rho_Z(t) is an extra contribution proportional to H^2(t). We perform a detailed comparison of such models with CMB, SNIa and BAO data.

The excursion set theory, where density perturbations evolve stochastically with the smoothing scale, provides a method for computing the dark matter halo mass function. The computation of the mass function is mapped into the so-called first-passage time problem in the presence of a moving barrier. The excursion set theory is also a powerful formalism to study other properties of dark matter halos such as halo bias, accretion rate, formation time, merging rate and the formation history of halos. This is achieved by computing conditional probabilities with non-trivial initial conditions, and the conditional two-barrier first-crossing rate. In this paper we use the path integral formulation of the excursion set theory to calculate analytically these conditional probabilities in the presence of a generic moving barrier, including the one describing the ellipsoidal collapse, and for both Gaussian and non-Gaussian initial conditions. While most of our analysis associated with Gaussian initial conditions assumes Markovianity (top-hat in momentum space smoothing, rather than generic filters), the non-Markovianity of the random walks induced by non-Gaussianity is consistently accounted for. We compute, for a generic barrier, the first two scale-independent halo bias parameters, the conditional mass function and the halo formation time probability, including the effects of non-Gaussianities. We also provide the expression for the two-constant-barrier first-crossing rate when non-Markovian effects are induced by a top-hat filter function in real space.

Excursion set theory, where density perturbations evolve stochastically with the smoothing scale, provides a method for computing the mass function of cosmological structures like dark matter halos, sheets and filaments. The computation of these mass functions is mapped into the so-called first-passage time problem in the presence of a moving barrier. In this paper we use the path integral formulation of the excursion set theory developed recently to analytically solve the first-passage time problem in the presence of a generic moving barrier, in particular the barrier corresponding to ellipsoidal collapse. We perform the computation for both Gaussian and non-Gaussian initial conditions. The expression of the halo mass function for the ellipsoidal collapse barrier and with non-Gaussianity is therefore obtained in a fully consistent way and it does not require the introduction of any form factor artificially derived from the Press-Schechter formalism based on the spherical collapse and usually adopted in the literature.

We examine a cosmological model with a dark energy density of the form $\rho_{DE}(t)=\rho_X(t)+\rho_Z(t)$, where $\rho_X$ is the component that accelerates the Hubble expansion at late times and $\rho_Z(t)$ is an extra contribution proportional to $H^2(t)$. This form of $\rho_Z(t)$ follows from the recent proposal that the contribution of zero-point fluctuations of quantum fields to the total energy density should be computed by subtracting the Minkowski-space result from that computed in the FRW space-time. We discuss theoretical arguments that support this subtraction. By definition, this eliminates the quartic divergence in the vacuum energy density responsible for the cosmological constant problem. We show that the remaining quadratic divergence can be reabsorbed into a redefinition of Newton's constant only under the assumption that the energy-momentum tensor of vacuum fluctuations is conserved in isolation. However, in the presence of an ultra-light scalar field $X$ with $m_X<H_0$, as typical of some dark energy models, the gravity effective action depends both on the gravitational field and on the $X$ field. In this case general covariance only requires the conservation of the total energy-momentum tensor, including both the classical term $T^X_{\mu\nu}$ and the vacuum expectation value of T_{\mu\nu}. If there is an exchange of energy between these two terms, there are potentially observable consequences. We construct an explicit model with an interaction between $\rho_X$ and $\rho_Z$ and we show that the total dark energy density $\rho_{DE}(t)=\rho_X(t)+\rho_Z(t)$ always remains a finite fraction of the critical density at any time, providing a specific model of early dark energy. We discuss the implication of this result for the coincidence problem and we estimate the model parameters by means of a full likelihood analysis using current CMB, SNe Ia and BAO data.

In the Hamiltonian formulation of General Relativity the energy associated to an asymptotically flat space-time with metric $g_{\mu\nu}$ is related to the Hamiltonian $H_{GR}$ by $E=H_{GR}[g_{\mu\nu}]-H_{\rm GR}[\eta_{\mu\nu}]$, where the subtraction of the flat-space contribution is necessary to get rid of an otherwise divergent boundary term. This classic result indicates that the energy associated to flat space does not gravitate. We apply the same principle to study the effect of zero-point fluctuations of quantum fields in cosmology, proposing that their contribution to the cosmic expansion is obtained computing the vacuum energy of quantum fields in a FRW space-time with Hubble parameter $H(t)$ and subtracting from it the flat-space contribution. [...] After renormalization, this produces a renormalized vacuum energy density $\sim M^2H^2(t)$, where $M$ is the scale where quantum gravity sets is, so for $M$ of order of the Planck mass a vacuum energy density of the order of the critical density can be obtained without any fine tuning. The counterterms can be chosen so that the renormalized energy density and pressure satisfy $p=w\rho$, with $w$ a parameter that can be fixed by comparison to the observed value, so in particular one can chose $w=-1$. An energy density evolving in time as $H^2(t)$ is however observationally excluded as an explanation for the dominant dark energy component which is responsible for the observed acceleration of the universe. We rather propose that zero-point vacuum fluctuations provide a new subdominant "dark" contribution to the cosmic expansion that, for a UV scale $M$ slightly smaller than the Planck mass, is consistent with existing limits and potentially detectable.

The excursion set theory based on spherical or ellipsoidal gravitational collapse provides an elegant analytic framework for calculating the mass function and the large-scale bias of dark matter haloes. This theory assumes that the perturbed density field evolves stochastically with the smoothing scale and exhibits Markovian random walks in the presence of a density barrier. Here we derive an analytic expression for the halo bias in a new theoretical model that incorporates non-Markovian extension of the excursion set theory with a stochastic barrier. This model allows us to handle non-Markovian random walks and to calculate perturbativly these corrections to the standard Markovian predictions for the halo mass function and halo bias. Our model contains only two parameters: kappa, which parameterizes the degree of non-Markovianity and whose exact value depends on the shape of the filter function used to smooth the density field, and a, which parameterizes the degree of stochasticity of the barrier. Appropriate choices of kappa and a in our new model can lead to a closer match to both the halo mass function and halo bias in the latest N-body simulations than the standard excursion set theory.

A sizeable level of non-Gaussianity in the primordial cosmological perturbations may be induced by a large trispectrum, i.e. by a large connected four-point correlation function. We compute the effect of a primordial non-Gaussian trispectrum on the halo mass function, within excursion set theory. We use the formalism that we have developed in a previous series of papers and which allows us to take into account the fact that, in the presence of non-Gaussianity, the stochastic evolution of the smoothed density field, as a function of the smoothing scale, is non-markovian. In the large mass limit, the leading-order term that we find agrees with the leading-order term of the results found in the literature using a more heuristic Press-Schecther (PS)-type approach. Our approach however also allows us to evaluate consistently the subleading terms, which depend not only on the four-point cumulant but also on derivatives of the four-point correlator, and which cannot be obtained within non-Gaussian extensions of PS theory. We perform explicitly the computation up to next-to-leading order.

In excursion set theory the computation of the halo mass function is mapped into a first-passage time process in the presence of a barrier, which in the spherical collapse model is a constant and in the ellipsoidal collapse model is a fixed function of the variance of the smoothed density field. However, N-body simulations show that dark matter halos grow through a mixture of smooth accretion, violent encounters and fragmentations, and modeling halo collapse as spherical, or even as ellipsoidal, is a significant oversimplification. We propose that some of the physical complications inherent to a realistic description of halo formation can be included in the excursion set theory framework, at least at an effective level, by taking into account that the critical value for collapse is not a fixed constant $\delta_c$, as in the spherical collapse model, nor a fixed function of the variance $\sigma$ of the smoothed density field, as in the ellipsoidal collapse model, but rather is itself a stochastic variable, whose scatter reflects a number of complicated aspects of the underlying dynamics. Solving the first-passage time problem in the presence of a diffusing barrier we find that the exponential factor in the Press-Schechter mass function changes from $\exp\{-\delta_c^2/2\sigma^2\}$ to $\exp\{-a\delta_c^2/2\sigma^2\}$, where $a=1/(1+D_B)$ and $D_B$ is the diffusion coefficient of the barrier. The numerical value of $D_B$, and therefore the corresponding value of $a$, depends among other things on the algorithm used for identifying halos. We discuss the physical origin of the stochasticity of the barrier and we compare with the mass function found in N-body simulations, for the same halo definition.[Abridged]

A classic method for computing the mass function of dark matter halos is provided by excursion set theory, where density perturbations evolve stochastically with the smoothing scale, and the problem of computing the probability of halo formation is mapped into the so-called first-passage time problem in the presence of a barrier. While the full dynamical complexity of halo formation can only be revealed through N-body simulations, excursion set theory provides a simple analytic framework for understanding various aspects of this complex process. In this series of paper we propose improvements of both technical and conceptual aspects of excursion set theory, and we explore up to which point the method can reproduce quantitatively the data from N-body simulations. In paper I of the series we show how to derive excursion set theory from a path integral formulation. This allows us both to derive rigorously the absorbing barrier boundary condition, that in the usual formulation is just postulated, and to deal analytically with the non-markovian nature of the random walk. Such a non-markovian dynamics inevitably enters when either the density is smoothed with filters such as the top-hat filter in coordinate space (which is the only filter associated to a well defined halo mass) or when one considers non-Gaussian fluctuations. In these cases, beside ``markovian'' terms, we find ``memory'' terms that reflect the non-markovianity of the evolution with the smoothing scale. We develop a general formalism for evaluating perturbatively these non-markovian corrections, and in this paper we perform explicitly the computation of the halo mass function for gaussian fluctuations, to first order in the non-markovian corrections due to the use of a tophat filter in coordinate space.

We compute the effect of primordial non-Gaussianity on the halo mass function, using excursion set theory. In the presence of non-Gaussianity the stochastic evolution of the smoothed density field, as a function of the smoothing scale, is non-markovian and beside "local" terms that generalize Press-Schechter (PS) theory, there are also "memory" terms, whose effect on the mass function can be computed using the formalism developed in the first paper of this series. We find that, when computing the effect of the three-point correlator on the mass function, a PS-like approach which consists in neglecting the cloud-in-cloud problem and in multiplying the final result by a fudge factor close to 2, is in principle not justified. When computed correctly in the framework of excursion set theory, in fact, the "local" contribution vanishes (for all odd-point correlators the contribution of the image gaussian cancels the Press-Schechter contribution rather than adding up), and the result comes entirely from non-trivial memory terms which are absent in PS theory. However it turns out that, in the limit of large halo masses, where the effect of non-Gaussianity is more relevant, these memory terms give a contribution which is the the same as that computed naively with PS theory, plus subleading terms depending on derivatives of the three-point correlator. We finally combine these results with the diffusive barrier model developed in the second paper of this series, and we find that the resulting mass function reproduces recent N-body simulations with non-Gaussian initial conditions, without the introduction of any ad hoc parameter.

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.

We discuss the possibility that galactic gravitational wave sources might give burst signals at a rate of several events per year, detectable by state-of-the-art detectors. We are stimulated by the results of the data collected by the EXPLORER and NAUTILUS bar detectors in the 2001 run, which suggest an excess of coincidences between the two detectors, when the resonant bars are orthogonal to the galactic plane. Signals due to the coalescence of galactic compact binaries fulfill the energy requirements but are problematic for lack of known candidates with the necessary merging rate. We examine the limits imposed by galactic dynamics on the mass loss of the Galaxy due to GW emission, and we use them to put constraints also on the GW radiation from exotic objects, like binaries made of primordial black holes. We discuss the possibility that the events are due to GW bursts coming repeatedly from a single or a few compact sources. We examine different possible realizations of this idea, such as accreting neutron stars, strange quark stars, and the highly magnetized neutron stars (``magnetars'') introduced to explain Soft Gamma Repeaters. Various possibilities are excluded or appear very unlikely, while others at present cannot be excluded.

Gravitational waves are generated during first-order phase transitions, either by turbolence or by bubble collisions. If the transition takes place at temperatures of the order of the electroweak scale, the frequency of these gravitational waves is today just within the band of the planned space interferometer LISA. We present a detailed analysis of the production of gravitational waves during an electroweak phase transition in different supersymmetric models where, contrary to the case of the Standard Model, the transition can be first order. We find that the stochastic background of gravitational waves generated by bubble nucleation can reach a maximum value h0^2 Omega_{gw} of order 10^{-10} - 10^{-11}, which is within the reach of the planned sensitivity of LISA, while turbolence can even produce signals at the level h0^2 Omega_{gw} \sim 10^{-9}. These values of h0^2 Omega_{gw} are obtained in the regions of the parameter space which can account for the generation of the baryon asymmetry at the electroweak scale.

Gravitational waves generated during a first-order electroweak phase transition have a typical frequency which today falls just within the band of the planned space interferometer LISA. Contrary to what happens in the Standard Model, in its supersymmetric extensions the electroweak phase transition may be strongly first order, providing a mechanism for generating the observed baryon asymmetry in the Universe. We show that during the same transition the production of gravitational waves can be rather sizable. While the energy density in gravitational waves can reach at most $h_0^2 \Omega_{\rm gw}\simeq 10^{-16}$ in the Minimal Supersymmetric Standard Model, in the Next-to-Minimal Supersymmetric Model, in some parameter range, $h_0^2 \Omega_{\rm gw}$ can be as high as $4\times 10^{-11}$. A stochastic background of gravitational waves of this intensity is within the reach of the planned sensitivity of LISA. Since in the Standard Model the background of gravitational waves is totally neglegible, its detection would also provide a rather unexpected experimental signal of supersymmetry and a tool to descriminate among supersymmetric models with different Higgs content.

We review the motivations for the search of stochastic backgrounds of gravitational waves and we compare the experimental sensitivities that can be reached in the near future with the existing bounds and with the theoretical predictions.

Gravitational-wave experiments with interferometers and with resonant masses can search for stochastic backgrounds of gravitational waves of cosmological origin. We review both experimental and theoretical aspects of the search for these backgrounds. We give a pedagogical derivation of the various relations that characterize the response of a detector to a stochastic background. We discuss the sensitivities of the large interferometers under constructions (LIGO, VIRGO, GEO600, TAMA300, AIGO) or planned (Avdanced LIGO, LISA) and of the presently operating resonant bars, and we give the sensitivities for various two-detectors correlations. We examine the existing limits on the energy density in gravitational waves from nucleosynthesis, COBE and pulsars, and their effects on theoretical predictions. We discuss general theoretical principles for order-of-magnitude estimates of cosmological production mechanisms, and then we turn to specific theoretical predictions from inflation, string cosmology, phase transitions, cosmic strings and other mechanisms. We finally compare with the stochastic backgrounds of astrophysical origin.

D-branes, topological defects in string theory on which string endpoints can live, may give new insight into the understanding of the cosmological evolution of the Universe at early epochs. We analyze the dynamics of D-branes in curved backgrounds and discuss the parameter space of M-theory as a function of the coupling constant and of the curvature of the Universe. We show that D-branes may be efficiently produced by gravitational effects. Furthermore, in curved spacetimes the transverse fluctuations of the D-branes develop a tachyonic mode and when the fluctuations grow larger than the horizon the branes become tensionless and break up. This signals a transition to a new regime. We discuss the implications of our findings for the singularity problem present in string cosmology, suggesting the existence of a limiting value for the curvature which is in agreement with the value suggested by the cosmological version of the holography principle. We also comment on possible implications for the so-called brane world scenario, where the Standard Model gauge and matter fields live inside some branes while gravitons live in the bulk.

We discuss the possible relevance of gravitational-wave (GW) experiments for physics at very high energy. We examine whether, from the experience gained with the computations of various specific relic GW backgrounds, we can extract statements and order of magnitude estimates that are as much as possible model-independent, and we try to distinguish between general conclusions and results related to specific cosmological mechanisms. We examine the statement that the Virgo/LIGO experiments probe the Universe at temperatures $T\sim 10^{7}-10^{10}$ GeV (or timescales $t\sim 10^{-20}-10^{-26}$ sec) and we consider the possibility that they could actually probe the Universe at much higher energy scales, including the typical scales of grand unification, string theory and quantum gravity. We consider possible scenarios, depending on how the inflationary paradigm is implemented. We discuss the prospects for detection with present and planned experiments. In particular, a second Virgo interferometer correlated with the planned one, and located within a few tens of kilometers from the first, could reach an interesting sensitivity for stochastic GWs of cosmological origin.

We compute the spectrum of relic gravitons in a model of string cosmology. In the low- and in the high-frequency limits we reproduce known results. The full spectrum, however, also displays a series of oscillations which could give a characteristic signature at the planned LIGO/VIRGO detectors. For special values of the parameters of the model the signal reaches its maximum already at frequencies accessible to LIGO and VIRGO and it is close to the sensitivity of first generation experiments.