We perform the ADM decomposition of a five-parameter family of quadratic non-metricity theories and study their conjugate momenta. After systematically identifying all possible conditions which can be imposed on the parameters such that different sets of primary constraints arise, we find that the five-parametric theory space can be compartmentalized into nine different sectors, based on the presence or absence of primary constraints. This classification allows to dismiss certain classes of theories as unphysical and invites further investigations into the remaining sectors, which may contain phenomenologically interesting modifications of General Relativity.

We consider a simpler geometrical formulation of General Relativity based on non-metricity, known as Coincident General Relativity. We study the ADM formulation of the theory and perform a detailed Hamiltonian analysis. We explicitly show the propagation of two physical degrees of freedom, as it should, even though the role of boundary terms and gauge conditions is significantly altered. This might represent an alternative promising new route for numerical relativity and canonical quantum gravity. We also give an outlook on the number of propagating degrees of freedom in non-linear extension of non-metricity scalar.

We investigate the nonlinear growth of large-scale structure in the generalised Proca theory, in which a self-interacting massive vector field plays the role of driving the acceleration of the cosmic expansion. Focusing to the Proca Lagrangian at cubic order -- the cubic vector Galileon model -- we derive the simplified equations for gravity as well as the longitudinal and transverse modes of the vector field under the weak-field and quasi-static approximations, and implement them in a modified version of the ECOSMOG N-body code. Our simulations incorporate the Vainshtein screening effect, which reconciles the fifth force propagated by the longitudinal mode of the cubic vector Galileon model with local tests of gravity. The results confirm that for all scales probed by the simulation, the transverse mode has a negligible impact on structure formation in a realistic cosmological setup. It is well known that in this model the strength of the fifth force is controlled by a free model parameter, which we denote as beta_3. By running a suite of cosmological simulations for different values of beta_3, we show that this parameter also determines the effectiveness of the Vainshtein screening. The model behaves identically to the cubic scalar Galileon for beta_3 going to zero, in which the fifth force is strong in unscreened regions but is efficiently screened in high-density regions. In the opposite limit, beta_3 going to infinity, the model approaches its `quintessence' counterpart, which has a vanishing fifth force but a modified expansion history compared to LambdaCDM. This endows the model with rich phenomenology, which will be investigated in future works.

Local determinations of the Hubble constant $H_0$ favour a higher value than Planck based on CMB and $\Lambda$CDM. Through a model-independent expansion, we show that low redshift ($z \lesssim 0.7$) data comprising baryon acoustic oscillations (BAO), cosmic chronometers and Type Ia supernovae has a preference for Quintessence models that lower $H_0$ relative to $\Lambda$CDM. In addition, we confirm that an exponential coupling to dark matter cannot alter this conclusion in the same redshift range. Our results leave open the possibility that a coupling in the matter-dominated epoch, potentially even in the dark ages, may yet save $H_0$ from sinking in the string theory Swampland.

Current and upcoming cosmological experiments open a new era of precision cosmology, thus demanding accurate theoretical predictions for cosmological observables. Because of the complexity of the codes delivering such predictions, reaching a high level of numerical accuracy is challenging. Among the codes already fulfilling this task, $\textsf{PyCosmo}$ is a Python based framework providing solutions to the Einstein-Boltzmann equations and accurate predictions for cosmological observables. In this work, we first describe how the observables are implemented. Then, we check the accuracy of the theoretical predictions for background quantities, power spectra and Limber and beyond-Limber angular power spectra by comparison with other codes: the Core Cosmology Library ($\texttt{CCL}$), $\texttt{CLASS}$, $\texttt{HMCode}$ and $\texttt{iCosmo}$. In our analysis we quantify the agreement of $\textsf{PyCosmo}$ with the other codes, for a range of cosmological models, monitored through a series of $\textit{unit tests}$. $\textsf{PyCosmo}$, conceived as a multi purpose cosmology calculation tool in $\texttt{Python}$, is designed to be interactive and user friendly. A current version of the code (without the Boltzmann Solver) is publicly available and can be used interactively on the platform $\textsf{PyCosmo Hub}$, all accessible from this link: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html . On the hub the users can perform their own computations using $\texttt{Jupyter Notebooks}$ without the need of installing any software, access to the results presented in this work and benefit from tutorial notebooks illustrating the usage of the code. The link above also redirects to the code release and documentation.

We consider the prospect of detecting cubic Galileons through their imprint on gravitational wave signals from a triple system. Namely, we consider a massive Black Hole (BH) surrounded by a binary system of two smaller BHs. We assume that the three BHs acquire a conformal coupling to the scalar field whose origin could be due to cosmology or to the galactic environment. In this case, the massive BH has a Vainshtein radius which englobes the smaller ones and suppresses the scalar effects on the motion of the binary system. On the other hand the two binaries can be outside each other's redressed Vainshtein radius calculated in the background of the central BH, allowing for a perturbative treatment of their dynamics. Despite the strong Vainshtein suppression, we find that the scalar effects on the binary system are slightly enhanced with respect to the static case and a significant amount of power can be emitted in the form of the Galileon scalar field, hence actively participating in the inspiralling phase. We compute the modification to the GW phase and show that it can lead to a detectable signal for large enough effective scalar coupling.

With recent constraints on the propagation speed of gravitational waves, the class of scalar-tensor theories has significantly been reduced. We consider one of the surviving models still relevant for cosmology and investigate its radiative stability. The model contains operators with explicit breaking of the Galileon symmetry and we study whether they harm the re-organization of the effective field theory. Within the regime of validity we establish a non-renormalization theorem and show explicitly that the quantum corrections, to one-loop, do not detune the classical Lagrangian generating suppressed counterterms. This is striking since the non-renormalization theorem is established in the presence of a genuine Galileon symmetry breaking term.

Gravitational wave (GW) oscillations occur whenever there are additional tensor modes interacting with the perturbations of the metric coupled to matter. These extra modes can arise from new spin-2 fields (as in e.g. bigravity theories) or from non-trivial realisations of the cosmological principle induced by background vector fields with internal symmetries (e.g. Yang-Mills, gaugids or multi-Proca). We develop a general cosmological framework to study such novel features due to oscillations. The evolution of the two tensor modes is described by a linear system of coupled second order differential equations exhibiting friction, velocity, chirality and mass mixing. We follow appropriate schemes to obtain approximate solutions for the evolution of both modes and show the corresponding phenomenology for different mixings. Observational signatures include modulations of the wave-form, oscillations of the GW luminosity distance, anomalous GW speed and chirality. We discuss the prospects of observing these effects with present and future GW observatories such as LIGO/VIRGO and LISA.

The geometrical formulation of gravity is not unique and can be set up in a variety of spacetimes. Even though the gravitational sector enjoys this freedom of different geometrical interpretations, consistent matter couplings have to be assured for a steady foundation of gravity. In generalised geometries, further ambiguities arise in the matter couplings unless the minimal coupling principle (MCP) is adopted that is compatible with the principles of relativity, universality and inertia. In this work, MCP is applied to all Standard Model gauge fields and matter fields in a completely general (linear) affine geometry. This is also discussed from an effective field theory perspective. It is found that the presence of torsion generically leads to theoretical problems. However, symmetric teleparallelism, wherein the affine geometry is integrable and torsion-free, is consistent with MCP. The generalised Bianchi identity is derived and shown to determine the dynamics of the connection in a unified fashion. Also, the parallel transport with respect to a teleparallel connection is shown to be free of second clock effects.

General Relativity enjoys the freedom of different geometrical interpretations in terms of curvature, torsion or non-metricity. Within this geometrical trinity, a simpler geometrical formulation of General Relativity manifests itself in the latter, where gravity is entirely attributed to non-metricity. In this Letter, we consider non-linear extensions of Coincident General Relativity $f(\mathring{\mathbb{Q}})$ for phenomenological applications on both cosmological as well as galactic scales. The theory not only delivers dark energy on large scales but also recovers MOND on galactic scales, together with implications for the early universe cosmology. To the best of our knowledge, this represents the first relativistic, covariant, and ghost-free hybrid-formulation of MOND which recovers both, General Relativity and MOND in the appropriate limits and reconciles expected cosmological behavior. We further illustrate that previous bimetric formulations of MOND generically suffer from ghost instabilities and $f(\mathring{\mathbb{Q}})$ crystalizes as a unique ghost-free theory.

In this Letter, we consider the model-independent reconstruction of the expansion and growth functions from the Pantheon supernova data. The method relies on developing the expansion function in terms of shifted Chebyshev polynomials and determining the coefficients of the polynomials by a maximum-likelihood fit to the data. Having obtained the expansion function in a model-independent way, we can then also determine the growth function without assuming a particular model. We then compare the results with the predictions of two classes of Dark Energy models, firstly a class of quintessence scalar field models consistent with the trans-Planckian censorship and swampland conjectures, and secondly a class of generalized Proca vector field models. We determine constraints on the parameters which appear in these models.

In this work we revisit the topological mass generation of 2-forms and establish a connection to the unique derivative coupling arising in the quartic Lagrangian of the systematic construction of massive $2-$form interactions, relating in this way BF theories to Galileon-like theories of 2-forms. In terms of a massless $1-$form $A$ and a massless $2-$form $B$, the topological term manifests itself as the interaction $B\wedge F$, where $F = {\rm d} A$ is the field strength of the $1-$form. Such an interaction leads to a mechanism of generation of mass, usually referred to as "topological generation of mass" in which the single degree of freedom propagated by the $2-$form is absorbed by the $1-$form, generating a massive mode for the $1-$form. Using the systematical construction in terms of the Levi-Civita tensor, it was shown that, apart from the quadratic and quartic Lagrangians, Galileon-like derivative self-interactions for the massive 2-form do not exist. A unique quartic Lagrangian $\epsilon^{\mu\nu\rho\sigma}\epsilon^{\alpha\beta\gamma}_{\;\;\;\;\;\;\sigma}\partial_{\mu}B_{\alpha\rho}\partial_{\nu}B_{\beta\gamma}$ arises in this construction in a way that it corresponds to a total derivative on its own but ceases to be so once an overall general function is introduced. We show that it exactly corresponds to the same interaction of topological mass generation. Based on the decoupling limit analysis of the interactions, we bring out supporting arguments for the uniqueness of such a topological mass term and absence of the Galileon-like interactions. Finally, we discuss some preliminary applications in cosmology.

We investigate the renormalization structure of scalar Galileons in flat spacetime. We explicitly calculate the ultraviolet divergent one-loop contributions to the 2-point, 3-point, 4-point, and 5-point functions. We discuss the structure of the counterterms and their hierarchy within an effective field theory expansion. We comment on different resummation schemes, including a geometric resummation for which our results could be generalized to arbitrary n-point functions.

Quasidilaton massive gravity offers a physically well-defined gravitational theory with non-zero graviton mass. We present the full set of dynamical equations governing the expansion history of the universe, valid during radiation domination, matter domination, and a late-time self-accelerating epoch related to the graviton mass. The existence of self-consistent solutions constrains the amplitude of the quasi-dilaton field and the graviton mass, as well as other model parameters. We point out that the effective mass of gravitational waves can be significantly larger than the graviton mass, opening the possibility that a single theory can explain both the late-time acceleration of the cosmic expansion and modifications of structure growth leading to the suppression of large-angle correlations observed in the cosmic microwave background.

Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

We investigate the implications of string Swampland criteria for alternative theories of gravity. Even though this has not rigorously been proven, there is some evidence that exact de-Sitter solutions with a positive cosmological constant cannot be successfully embedded into string theory, and that the low energy effective field theories containing a scalar field $\pi$ with a potential $V$ in the habitable landscape should satisfy the Swampland criteria $|V'|/V\ge c\sim\mathcal{O}(1)$. As a paradigmatic class of modified gravity theories for inflation and dark energy, we consider the extensively studied family of Horndeski Lagrangians in view of cosmological observations. Apart from a possible potential term, they contain derivative self-interactions as the Galileon and non-minimal couplings to the gravity sector. Hence, our conclusions on the Galileon sector can be also applied to many other alternative theories with scalar helicity modes containing derivative interactions such as massive gravity and generalized Proca. In the presence of such derivative terms, the dynamics of the scalar mode is substantially modified, and imposing the cosmological evolution constrained by observations places tight constraints on $c$ within the Swampland conjecture.

We consider the Lagrangian of a vector field with derivative self-interactions with a priori arbitrary coefficients. Starting with a flat space-time we show that for a special choice of the coefficients of the self-interactions the ghost-like pathologies disappear. For this we use the degeneracy condition of the Hessian. This constitutes the Galileon-type generalization of the Proca action with only three propagating physical degrees of freedom. The longitudinal mode of the vector field is associated to the usual Galileon interactions for a specific choice of the overall functions. In difference to a scalar Galileon theory, the generalized Proca field has more free parameters and purely intrinsic vector interactions. We then extend this analysis to a curved background. The resulting theory is the Horndeski Proca action with second order equations of motion on curved space-times.

The formation of cosmic structures is an important diagnostic for both the dynamics of the cosmological model and the underlying theory of gravity. At the linear level of these structures, certain degeneracies remain between different cosmological models and alternative gravity theories. It is thus indispensable to study the non-linear, late-time evolution of cosmic structures to try and disentangle their fundamental properties caused by the cosmological model or gravity theory itself. Conventionally, non-linear cosmic structure formation is studied by means of computationally expensive numerical simulations. Since these inevitably suffer from shot noise and are too time consuming to systematically scrutinize large parameter spaces of cosmological models or fundamental theories, analytical methods are needed to overcome the limitations of numerical simulations. Recently, a new analytic approach to non-linear cosmic structure formation has been proposed based on kinetic field theory for classical particle ensembles. Within this theory, a closed, analytic, non-perturbative and parameter-free equation could be derived for the non-linear power spectrum of cosmic density perturbations which agrees very well with numerically simulated results to wave numbers $k\lesssim10\,h\,\mathrm{Mpc}^{-1}$ at redshift $z = 0$. In this Letter, we study for the first time the implications of alternative gravity theories for non-linear cosmic structure formation applying this promising new analytic framework. As an illustrative example, we consider vector-tensor theories, which support very interesting isotropic cosmological solutions.

It has been suggested that low energy effective field theories should satisfy given conditions in order to be successfully embedded into string theory. In the case of a single canonically normalized scalar field this translates into conditions on its potential and the derivatives thereof. In this Letter we revisit stochastic models of small field inflation and study the compatibility of the swampland constraints with entropy considerations. We show that stochastic inflation either violates entropy bounds or the swampland criterium on the slope of the scalar field potential. Furthermore, we illustrate that such models are faced with a graceful exit problem: any patch of space which exits the region of eternal inflation is either not large enough to explain the isotropy of the cosmic microwave background, or has a spectrum of fluctuations with an unacceptably large red tilt.

The universal character of the gravitational interaction provided by the equivalence principle motivates a geometrical description of gravity. The standard formulation of General Relativity \`a la Einstein attributes gravity to the spacetime curvature, to which we have grown accustomed. However, this perception has masked the fact that two alternative, though equivalent, formulations of General Relativity in flat spacetimes exist, where gravity can be fully ascribed either to torsion or to non-metricity. The latter allows a simpler geometrical formulation of General Relativity that is oblivious to the affine spacetime structure. Generalisations along this line permit to generate teleparallel and symmetric teleparallel theories of gravity with exceptional properties. In this work we explore modified gravity theories based on non-linear extensions of the non-metricity scalar. After presenting some general properties and briefly studying some interesting background cosmologies (including accelerating solutions with relevance for inflation and dark energy), we analyse the behaviour of the cosmological perturbations. Tensor perturbations feature a re-scaling of the corresponding Newton's constant, while vector perturbations do not contribute in the absence of vector sources. In the scalar sector we find two additional propagating modes, hinting that $f(Q)$ theories introduce, at least, two additional degrees of freedom. These scalar modes disappear around maximally symmetric backgrounds because of the appearance of an accidental residual gauge symmetry corresponding to a restricted diffeomorphism. We finally discuss the potential strong coupling problems of these maximally symmetric backgrounds caused by the discontinuity in the number of propagating modes.

We reconsider the construction of general derivative self-interactions for a massive Proca field. The constructed Lagrangian is such that the vector field propagates at most three degrees of freedom, thus avoiding the ghostly nature of a fourth polarisation. The construction makes use of the well-known condition for constrained systems of having a degenerate Hessian. We briefly discuss the casuistry according to the nature of the existing constraints algebra. We also explore various classes of interesting new interactions that have been recently raised in the literature. For the sixth order Lagrangian that satisfies the constraints by itself we prove its topological character, making such a term irrelevant. There is however a window of opportunity for exploring other classes of fully-nonlinear interactions that satisfy the constraint algebra by mixing terms of various order.

Kinetic Field Theory (KFT) is a statistical field theory for an ensemble of point-like classical particles in or out of equilibrium. We review its application to cosmological structure formation. Beginning with the construction of the generating functional of the theory, we describe in detail how the theory needs to be adapted to reflect the expanding spatial background and the homogeneous and isotropic, correlated initial conditions for cosmic structures. Based on the generating functional, we develop three main approaches to non-linear, late-time cosmic structures, which rest either on the Taylor expansion of an interaction operator, suitable averaging procedures for the interaction term, or a resummation of perturbation terms. We show how an analytic, parameter-free equation for the non-linear cosmic power spectrum can be derived. We explain how the theory can be used to derive the density profile of gravitationally bound structures and use it to derive power spectra of cosmic velocity densities. We further clarify how KFT relates to the BBGKY hierarchy. We then proceed to apply kinetic field theory to fluids, introduce a reformulation of KFT in terms of macroscopic quantities which leads to a resummation scheme, and use this to describe mixtures of gas and dark matter. We discuss how KFT can be applied to study cosmic structure formation with modified theories of gravity. As an example for an application to a non-cosmological particle ensemble, we show results on the spatial correlation function of cold Rydberg atoms derived from KFT.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

In scalar-vector-tensor (SVT) theories with parity invariance, we perform a gauge-ready formulation of cosmological perturbations on the flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) background by taking into account a matter perfect fluid. We derive the second-order action of scalar perturbations and resulting linear perturbation equations of motion without fixing any gauge conditions. Depending on physical problems at hand, most convenient gauges can be chosen to study the development of inhomogeneities in the presence of scalar and vector fields coupled to gravity. This versatile framework, which encompasses Horndeski and generalized Proca theories as special cases, is applicable to a wide variety of cosmological phenomena including nonsingular cosmology, inflation, and dark energy. By deriving conditions for the absence of ghost and Laplacian instabilities in several different gauges, we show that, unlike Horndeski theories, it is possible to evade no-go arguments for the absence of stable nonsingular bouncing/genesis solutions in both generalized Proca and SVT theories. We also apply our framework to the case in which scalar and vector fields are responsible for dark energy and find that the separation of observables relevant to the evolution of matter perturbations into tensor, vector, and scalar sectors is transparent in the unitary gauge. Unlike the flat gauge chosen in the literature, this result is convenient to confront SVT theories with observations associated with the cosmic growth history.

In the framework of effective field theories with prominent helicity-0 and helicity-1 fields coupled to each other via a dimension-3 operator, we study the dynamics of inflation driven by the helicity-0 mode, with a given potential energy, as well as the evolution of cosmological perturbations, influenced by the presence of a mixing term between both helicities. In this scenario, the temporal component of the helicity-1 mode is an auxiliary field and can be integrated out in terms of the time derivative of the helicity-0 mode, so that the background dynamics effectively reduces to that in single-field inflation modulated by a parameter $\beta$ associated to the coupling between helicity-0 and helicity-1 modes. We discuss the evolution of a longitudinal scalar perturbation $\psi$ and an inflaton fluctuation $\delta \phi$, and explicitly show that a particular combination of these two, which corresponds to an isocurvature mode, is subject to exponential suppression by the vector mass comparable to the Hubble expansion rate during inflation. Furthermore, we find that the effective single-field description corrected by $\beta$ also holds for the power spectrum of curvature perturbations generated during inflation. We compute the standard inflationary observables such as the scalar spectral index $n_s$ and the tensor-to-scalar ratio $r$ and confront several inflaton potentials with the recent observational data provided by Planck 2018. Our results show that the coupling between helicity-0 and helicity-1 modes can lead to a smaller value of the tensor-to-scalar ratio especially for small-field inflationary models, so our scenario exhibits even better compatibility with the current observational data.

The metric-affine variational principle is applied to generate teleparallel and symmetric teleparallel theories of gravity. From the latter is discovered an exceptional class which is consistent with a vanishing affine connection. Based on this remarkable property, this work proposes a simpler geometrical formulation of General Relativity that is oblivious to the affine spacetime structure, thus fundamentally depriving gravity of any inertial character. The resulting theory is described by the Hilbert action purged from the boundary term and is more robustly underpinned by the spin-2 field theory, where an extra symmetry is now manifest, possibly related to the double copy structure of the gravity amplitudes. This construction also provides a novel starting point for modified gravity theories, and the paper presents new and simple generalisations where analytical self-accelerating cosmological solutions arise naturally in the early and late time universe.

In this recent Letter [1], we studied the implications of string Swampland conjectures for quintessence models of dark energy. A couple of days ago, a paper appeared [2] which refers to our work but misrepresents what we have done. Here, we set the record straight.

In this Letter, we study the implications of string Swampland criteria for dark energy in view of ongoing and future cosmological observations. If string theory should be the ultimate quantum gravity theory, there is evidence that exact de Sitter solutions with a positive cosmological constant cannot describe the fate of the late-time universe. Even though cosmological models with dark energy given by a scalar field $\pi$ evolving in time are not in direct tension with string theory, they have to satisfy the Swampland criteria $|\Delta\pi|<d\sim\mathcal{O}(1)$ and $|V'|/V>c\sim\mathcal{O}(1)$, where $V$ is the scalar field potential. In view of the restrictive implications that the Swampland criteria have on dark energy, we investigate the accuracy needed for future observations to tightly constrain standard dark-energy models. We find that current 3-$\sigma$ constraints with $c \lesssim 1.35$ are still well in agreement with the string Swampland criteria. However, Stage-4 surveys such as Euclid, LSST and DESI, tightly constraining the equation of state $w(z)$, will start putting surviving quintessence models into tensions with the string Swampland criteria by demanding $c<0.4$. We further investigate whether any idealised futuristic survey will ever be able to give a decisive answer to the question whether the cosmological constant would be preferred over a time-evolving dark-energy model within the Swampland criteria. Hypothetical surveys with a reduction in the uncertainties by a factor of $\sim20$ compared to Euclid would be necessary to reveal strong tension between quintessence models obeying the string Swampland criteria and observations by pushing the allowed values down to $c<0.1$. In view of such perspectives, there will be fundamental observational limitations with future surveys.

We study the cosmology on the Friedmann-Lemaitre-Robertson-Walker background in scalar-vector-tensor theories with a broken $U(1)$ gauge symmetry. For parity-invariant interactions arising in scalar-vector-tensor theories with second-order equations of motion, we derive conditions for the absence of ghosts and Laplacian instabilities associated with tensor, vector, and scalar perturbations at linear order. This general result is applied to the computation of the primordial tensor power spectrum generated during inflation as well as to the speed of gravity relevant to dark energy. We also construct a concrete inflationary model in which a temporal vector component $A_0$ contributes to the dynamics of cosmic acceleration besides a scalar field $\phi$ through their kinetic mixings. In this model, we show that all the stability conditions of perturbations can be consistently satisfied during inflation and subsequent reheating.

A century ago, Einstein formulated his elegant and elaborate theory of General Relativity, which has so far withstood a multitude of empirical tests with remarkable success. Notwithstanding the triumphs of Einstein's theory, the tenacious challenges of modern cosmology and of particle physics have motivated the exploration of further generalised theories of spacetime. Even though Einstein's interpretation of gravity in terms of the curvature of spacetime is commonly adopted, the assignment of geometrical concepts to gravity is ambiguous because General Relativity allows three entirely different, but equivalent approaches of which Einstein's interpretation is only one. From a field-theoretical perspective, however, the construction of a consistent theory for a Lorentz-invariant massless spin-2 particle uniquely leads to General Relativity. Keeping Lorentz invariance then implies that any modification of General Relativity will inevitably introduce additional propagating degrees of freedom into the gravity sector. Adopting this perspective, we will review the recent progress in constructing consistent field theories of gravity based on additional scalar, vector and tensor fields. Within this conceptual framework, we will discuss theories with Galileons, with Lagrange densities as constructed by Horndeski and beyond, extended to DHOST interactions, or containing generalized Proca fields and extensions thereof, or several Proca fields, as well as bigravity theories and scalar-vector-tensor theories. We will review the motivation of their inception, different formulations, and essential results obtained within these classes of theories together with their empirical viability.

We study the evolution of gravitational waves for non-singular cosmological solutions within the framework of Born-Infeld inspired gravity theories, with special emphasis on the Eddington-inspired Born-Infeld theory. We review the existence of two types of non-singular cosmologies, namely bouncing and asymptotically Minkowski solutions, from a perspective that makes their features more apparent. We study in detail the propagation of gravitational waves near these non-singular solutions and carefully discuss the origin and severity of the instabilities and strong coupling problems that appear. We also investigate the role of the adiabatic sound speed of the matter sector in the regularisation of the gravitational waves evolution. We extend our analysis to more general Born-Infeld inspired theories where analogous solutions are found. As a general conclusion, we obtain that the bouncing solutions are generally more prone to instabilities, while the asymptotically Minkowski solutions can be rendered stable, making them appealing models for the early universe.

In scalar-vector-tensor theories with $U(1)$ gauge invariance, it was recently shown that there exists a new type of hairy black hole (BH) solutions induced by a cubic-order scalar-vector interaction. In this paper, we derive conditions for the absence of ghosts and Laplacian instabilities against odd-parity perturbations on a static and spherically symmetric background for most general $U(1)$ gauge-invariant scalar-vector-tensor theories with second-order equations of motion. We apply those conditions to hairy BH solutions arising from the cubic-order coupling and show that the odd-parity stability in the gravity sector is always ensured outside the event horizon with the speed of gravity equivalent to that of light. We also study the case in which quartic-order interactions are present in addition to the cubic coupling and obtain conditions under which black holes are stable against odd-parity perturbations.

The detection of the GW170817/GRB170817A event improved the constraints on the propagation speed of gravitational waves, thus placing possible variations caused by dark energy under restraint. For models based on scalar fields belonging to the family of Horndeski Lagrangians, non-minimal derivative couplings are now severely constrained, entailing a substantially limited phenomenology. In this work we want to stress that there is still a plethora of dark energy models that get around this obstacle while still providing interesting phenomenologies able to distinguish them from the standard cosmology. We focus on a class involving vector fields as a proxy, but our discussion is extensible to a broader class of models. In particular, we show the possibility of having a non-minimal derivative coupling giving a non-trivial effect on scalar modes without affecting gravitational waves and the possibility of having a second tensor mode that can oscillate into gravitational waves. We also present a novel class of configurations breaking rotational invariance but with an energy-momentum tensor that is isotropic on-shell. This peculiar feature makes the scalar and vector sectors of the perturbations mix so that, even in a perfectly isotropic background cosmology, preferred direction effects can appear in the perturbations. We also comment on models that give rise to isotropic solutions when averaging over rapid oscillations of the vector fields. The explored models are classified according to distinctive field configurations that provide inequivalent realisations of the Cosmological Principle.

The Palatini formalism is developed for gravitational theories in flat geometries. We focus on two particularly interesting scenarios. First, we fix the connection to be metric compatible, but we follow a completely covariant approach by imposing the constraints with suitable Lagrange multipliers. For a general quadratic theory we show how torsion naturally propagates and we reproduce the Teleparallel Equivalent of General Relativity as a particular quadratic action that features an additional Lorentz symmetry. We then study the much less explored theories formulated in a geometry with neither curvature nor torsion, so that all the geometrical information is encoded in the non-metricity. We discuss how this geometrical framework leads to a purely inertial connection that can thus be completely removed by a coordinate gauge choice, the coincident gauge. From the quadratic theory we recover a simpler formulation of General Relativity in the form of the Einstein action, which enjoys an enhanced symmetry that reduces to a second linearised diffeomorphism at linear order. More general theories in both geometries can be formulated consistently by taking into account the inertial connection and the associated additional degrees of freedom. As immediate applications, the new cosmological equations and their Newtonian limit are considered, where the role of the lapse in the consistency of the equations is clarified, and the Schwarzschild black hole entropy is computed by evaluating the corresponding Euclidean action. We discuss how the boundary terms in the usual formulation of General Relativity are related to different choices of coordinates in its coincident version and show that in isotropic coordinates the Euclidean action is finite without the need to introduce boundary or normalisation terms.

We construct the consistent ghost-free covariant scalar-vector-tensor gravity theories with second order equations of motion with derivative interactions. We impose locality, unitarity, Lorentz invariance and pseudo-Riemannian geometry as the fundamental terms. In the tensor sector we require diffeomorphism invariance, whereas we allow the vector sector to be gauge invariant or not. The resulting Lagrangians consist of new genuine couplings among these fields with at most two derivatives per field. They propagate five physical degrees of freedom in the gauge invariant case and six degrees of freedom if the gauge invariance is broken. In the corresponding limit of the free general functions in the Lagrangian, one recovers the generalized Proca theories. These scalar-vector-tensor theories will have important implications for cosmological and astrophysical applications, among which we mention the application to inflation and generation of primordial magnetic fields, new black hole and neutron star solutions, dark matter and dark energy.

The recent detection of the gravitational wave signal GW170817 together with an electromagnetic counterpart GRB 170817A from the merger of two neutron stars puts a stringent bound on the tensor propagation speed. This constraint can be automatically satisfied in the framework of massive gravity. In this work we consider a general $SO(3)$-invariant massive gravity with five propagating degrees of freedom and derive the conditions for the absence of ghosts and Laplacian instabilities in the presence of a matter perfect fluid on the flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) cosmological background. The graviton potential containing the dependence of three-dimensional metrics and a fiducial metric coupled to a temporal scalar field gives rise to a scenario of the late-time cosmic acceleration in which the dark energy equation of state $w_{\rm DE}$ is equivalent to $-1$ or varies in time. We find that the deviation from the value $w_{\rm DE}=-1$ provides important contributions to the quantities associated with the stability conditions of tensor, vector, and scalar perturbations. In concrete models, we study the dynamics of dark energy arising from the graviton potential and show that there exist viable parameter spaces in which neither ghosts nor Laplacian instabilities are present for both $w_{\rm DE}>-1$ and $w_{\rm DE}<-1$. We also generally obtain the effective gravitational coupling $G_{\rm eff}$ with non-relativistic matter as well as the gravitational slip parameter $\eta_s$ associated with the observations of large-scale structures and weak lensing. We show that, apart from a specific case, the two quantities $G_{\rm eff}$ and $\eta_s$ are similar to those in general relativity for scalar perturbations deep inside the sound horizon.

We present a family of exact black-hole solutions on a static spherically symmetric background in second-order generalized Proca theories with derivative vector-field interactions coupled to gravity. We also derive non-exact solutions in power-law coupling models including vector Galileons and numerically show the existence of regular black holes with a primary hair associated with the longitudinal propagation. The intrinsic vector-field derivative interactions generally give rise to a secondary hair induced by non-trivial field profiles. The deviation from General Relativity is most significant around the horizon and hence there is a golden opportunity for probing the Proca hair by the measurements of gravitational waves (GWs) in the regime of strong gravity.

As wide-field surveys yield ever more precise measurements, cosmology has entered a phase of high precision requiring highly accurate and fast theoretical predictions. At the heart of most cosmological model predictions is a numerical solution of the Einstein-Boltzmann equations governing the evolution of linear perturbations in the Universe. We present PyCosmo, a new Python-based framework to solve this set of equations using a special pur- pose solver based on symbolic manipulations, automatic generation of C++ code and sparsity optimisation. The code uses a consistency relation of the field equations to adapt the time step and does not rely on physical approximations for speed-up. After reviewing the system of first-order linear homogeneous differential equations to be solved, we describe the numerical scheme implemented in PyCosmo. We then compare the predictions and performance of the code for the computation of the transfer functions of cosmological perturbations and compare it to existing cosmological Boltzmann codes. We find that we achieve comparable execution times for comparable accuracies. While PyCosmo does not yet have all the features of other codes, our approach is complementary to existing cosmological Boltzmann solvers and can be used as an independent test of their numerical solutions. The symbolic representation of the Einstein-Boltzmann equation system in PyCosmo provides a convenient interface for implementing extended cosmological models. We also discuss how the PyCosmo framework can also be used as a general framework to compute cosmological quantities as well as observables for both interactive and high-performance batch jobs applications. Information about the PyCosmo package and future code releases are available at http://www.cosmology.ethz.ch/research/software-lab.html.

In a model of the late-time cosmic acceleration within the framework of generalized Proca theories, there exists a de Sitter attractor preceded by the dark energy equation of state $w_{\rm DE}=-1-s$, where $s$ is a positive constant. We run the Markov-Chain-Monte-Carlo code to confront the model with the observational data of Cosmic Microwave Background (CMB), baryon acoustic oscillations, supernovae type Ia, and local measurements of the Hubble expansion rate for the background cosmological solutions and obtain the bound $s=0.254^{{}+ 0.118}_{{}-0.097}$ at 95% confidence level (CL). Existence of the additional parameter $s$ to those in the $\Lambda$-Cold-Dark-Matter ($\Lambda$CDM) model allows to reduce tensions of the Hubble constant $H_0$ between the CMB and the low-redshift measurements. Including the cosmic growth data of redshift-space distortions in the galaxy power spectrum and taking into account no-ghost and stability conditions of cosmological perturbations, we find that the bound on $s$ is shifted to $s=0.16^{+0.08}_{-0.08}$ (95 % CL) and hence the model with $s>0$ is still favored over the $\Lambda$CDM model. Apart from the quantities $s, H_0$ and the today's matter density parameter $\Omega_{m0}$, the constraints on other model parameters associated with perturbations are less stringent, reflecting the fact that there are different sets of parameters that give rise to a similar cosmic expansion and growth history.

We study static and spherically symmetric black hole (BH) solutions in second-order generalized Proca theories with nonminimal vector field derivative couplings to the Ricci scalar, the Einstein tensor, and the double dual Riemann tensor. We find concrete Lagrangians which give rise to exact BH solutions by imposing two conditions of the two identical metric components and the constant norm of the vector field. These exact solutions are described by either Reissner-Nordstr\"{o}m (RN), stealth Schwarzschild, or extremal RN solutions with a non-trivial longitudinal mode of the vector field. We then numerically construct BH solutions without imposing these conditions. For cubic and quartic Lagrangians with power-law couplings which encompass vector Galileons as the specific cases, we show the existence of BH solutions with the difference between two non-trivial metric components. The quintic-order power-law couplings do not give rise to non-trivial BH solutions regular throughout the horizon exterior. The sixth-order and intrinsic vector-mode couplings can lead to BH solutions with a secondary hair. For all the solutions, the vector field is regular at least at the future or past horizon. The deviation from General Relativity induced by the Proca hair can be potentially tested by future measurements of gravitational waves in the nonlinear regime of gravity.

General Relativity has shown an outstanding observational success in the scales where it has been directly tested. However, modifications have been intensively explored in the regimes where it seems either incomplete or signals its own limit of validity. In particular, the breakdown of unitarity near the Planck scale strongly suggests that General Relativity needs to be modified at high energies and quantum gravity effects are expected to be important. This is related to the existence of spacetime singularities when the solutions of General Relativity are extrapolated to regimes where curvatures are large. In this sense, Born-Infeld inspired modifications of gravity have shown an extraordinary ability to regularise the gravitational dynamics, leading to non-singular cosmologies and regular black hole spacetimes in a very robust manner and without resorting to quantum gravity effects. This has boosted the interest in these theories in applications to stellar structure, compact objects, inflationary scenarios, cosmological singularities, and black hole and wormhole physics, among others. We review the motivations, various formulations, and main results achieved within these theories, including their observational viability, and provide an overview of current open problems and future research opportunities.

In this work we will summarise the recent progress made in constructing consistent theories for a massive vector field with derivative self-interactions. The construction is such that only the three desired polarisations of the Proca field propagate. We apply a systematic construction of the interactions by using the anti-symmetric Levi-Civita tensors. This finite family of allowed derivative self-interactions can be also obtained from the decoupling limit by imposing that the St\"uckelberg field only contains second derivatives both with itself and with the transverse modes. These interactions can be generalised to curved backgrounds, which relies on the presence of non-minimal couplings and constitutes a general family of vector-tensor interactions. We discuss also some extensions of these interactions by alleviating the restriction of second order nature of the equations or by imposing global symmetries. We will also comment on their interesting cosmological applications.

A non-abelian $SU(2)$ gauge field with a non-minimal Horndeski coupling to gravity gives rise to a de Sitter solution followed by a graceful exit to a radiation-dominated epoch. In this Horndeski Yang-Mills (HYM) theory we derive the second-order action for tensor perturbations on the homogeneous and isotropic quasi de Sitter background. We find that the presence of the Horndeski non-minimal coupling to the gauge field inevitably introduces ghost instabilities in the tensor sector during inflation. Moreover, we also find Laplacian instabilities for the tensor perturbations deep inside the Hubble radius during inflation. Thus, we conclude that the HYM theory does not provide a consistent inflationary framework due to the presence of ghosts and Laplacian instabilities.

Dipolar Dark Matter (DDM) is an alternative model motivated by the challenges faced by the standard cold dark matter model to describe the right phenomenology at galactic scales. A promising realisation of DDM was recently proposed in the context of massive bigravity theory. The model contains dark matter particles, as well as a vector field coupled to the effective composite metric of bigravity. This model is completely safe in the gravitational sector thanks to the underlying properties of massive bigravity. In this work we investigate the exact decoupling limit of the theory, including the contribution of the matter sector, and prove that it is free of ghosts in this limit. We conclude that the theory is acceptable as an Effective Field Theory below the strong coupling scale.

In Gleyzes-Langlois-Piazza-Vernizzi (GLPV) scalar-tensor theories, which are outside the domain of second-order Horndeski theories, it is known that there exists a solid angle deficit singularity in the case where the parameter $\alpha_{\rm H}$ characterizing the deviation from Horndeski theories approaches a non-vanishing constant at the center of a spherically symmetric body. Meanwhile, it was recently shown that second-order generalized Proca theories with a massive vector field $A^{\mu}$ can be consistently extended to beyond-generalized Proca theories, which recover shift-symmetric GLPV theories in the scalar limit $A^{\mu} \to \nabla^{\mu} \chi$. In beyond-generalized Proca theories up to quartic-order Lagrangians, we show that solid angle deficit singularities are generally absent due to the existence of a temporal vector component. We also derive the vector-field profiles around a compact object and show that the success of the Vainshtein mechanism operated by vector Galileons is not prevented by new interactions in beyond generalized Proca theories.

In this work we revisit the construction of theories for a massive vector field with derivative self-interactions such that only the 3 desired polarizations corresponding to a Proca field propagate. We start from the decoupling limit by constructing healthy interactions containing second derivatives of the Stueckelberg field with itself and also with the transverse modes. The resulting interactions can then be straightforwardly generalized beyond the decoupling limit. We then proceed to a systematic construction of the interactions by using the Levi-Civita tensors. Both approaches lead to a finite family of allowed derivative self-interactions for the Proca field. This construction allows us to show that some higher order terms recently introduced as new interactions trivialize in 4 dimensions by virtue of the Cayley-Hamilton theorem. Moreover, we discuss how the resulting derivative interactions can be written in a compact determinantal form, which can also be regarded as a generalization of the Born-Infeld lagrangian for electromagnetism. Finally, we generalize our results for a curved background and give the necessary non-minimal couplings guaranteeing that no additional polarizations propagate even in the presence of gravity.

We extend previous results on healthy derivative self-interactions for a Proca field to the case of a set of massive vector fields. We obtain non-gauge invariant derivative self-interactions for the vector fields that maintain the appropriate number of propagating degrees of freedom. In view of the potential cosmological applications, we restrict to interactions with an internal rotational symmetry. We provide a systematical construction order by order in derivatives of the fields and making use of the antisymmetric Levi-Civita tensor. We then compare with the one single vector field case and show that the interactions can be broadly divided into two groups, namely the ones obtained from a direct extension of the generalized Proca terms and genuine multi-Proca interactions with no correspondence in the single Proca case. We also discuss the curved spacetime version of the interactions to include the necessary non-minimal couplings to gravity. Finally, we explore the cosmological applications and show that there are three different vector fields configurations giving rise to isotropic solutions. Two of them have already been considered in the literature and the third one, representing a combination of the first two, is new and offers unexplored cosmological scenarios.

We consider higher-order derivative interactions beyond second-order generalized Proca theories that propagate only the three desired polarizations of a massive vector field besides the two tensor polarizations from gravity. These new interactions follow the similar construction criteria to those arising in the extension of scalar-tensor Horndeski theories to Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories. On the isotropic cosmological background, we show the existence of a constraint with a vanishing Hamiltonian that removes the would-be Ostrogradski ghost. We study the behavior of linear perturbations on top of the isotropic cosmological background in the presence of a matter perfect fluid and find the same number of propagating degrees of freedom as in generalized Proca theories (two tensor polarizations, two transverse vector modes, and two scalar modes). Moreover, we obtain the conditions for the avoidance of ghosts and Laplacian instabilities of tensor, vector, and scalar perturbations. We observe key differences in the scalar sound speed, which is mixed with the matter sound speed outside the domain of generalized Proca theories.

We consider the finite interactions of the generalized Proca theory including the sixth-order Lagrangian and derive the full linear perturbation equations of motion on the flat Friedmann-Lema\^{i}tre-Robertson-Walker background in the presence of a matter perfect fluid. By construction, the propagating degrees of freedom (besides the matter perfect fluid) are two transverse vector perturbations, one longitudinal scalar, and two tensor polarizations. The Lagrangians associated with intrinsic vector modes neither affect the background equations of motion nor the second-order action of tensor perturbations, but they do give rise to non-trivial modifications to the no-ghost condition of vector perturbations and to the propagation speeds of vector and scalar perturbations. We derive the effective gravitational coupling $G_{\rm eff}$ with matter density perturbations under a quasi-static approximation on scales deep inside the sound horizon. We find that the existence of intrinsic vector modes allows a possibility for reducing $G_{\rm eff}$. In fact, within the parameter space, $G_{\rm eff}$ can be even smaller than the Newton gravitational constant $G$ at the late cosmological epoch, with a peculiar phantom dark energy equation of state (without ghosts). The modifications to the slip parameter $\eta$ and the evolution of growth rate $f\sigma_8$ are discussed as well. Thus, dark energy models in the framework of generalized Proca theories can be observationally distinguished from the $\Lambda$CDM model according to both cosmic growth and expansion history. Furthermore, we study the evolution of vector perturbations and show that outside the vector sound horizon the perturbations are nearly frozen and start to decay with oscillations after the horizon entry.

In beyond-generalized Proca theories including the extension to theories higher than second order, we study the role of a spatial component $v$ of a massive vector field on the anisotropic cosmological background. We show that, as in the case of the isotropic cosmological background, there is no additional ghostly degrees of freedom associated with the Ostrogradski instability. In second-order generalized Proca theories we find the existence of anisotropic solutions on which the ratio between the anisotropic expansion rate $\Sigma$ and the isotropic expansion rate $H$ remains nearly constant in the radiation-dominated epoch. In the regime where $\Sigma/H$ is constant, the spatial vector component $v$ works as a dark radiation with the equation of state close to $1/3$. During the matter era, the ratio $\Sigma/H$ decreases with the decrease of $v$. As long as the conditions $|\Sigma| \ll H$ and $v^2 \ll \phi^2$ are satisfied around the onset of late-time cosmic acceleration, where $\phi$ is the temporal vector component, we find that the solutions approach the isotropic de Sitter fixed point ($\Sigma=0=v$) in accordance with the cosmic no-hair conjecture. In the presence of $v$ and $\Sigma$ the early evolution of the dark energy equation of state $w_{\rm DE}$ in the radiation era is different from that in the isotropic case, but the approach to the isotropic value $w_{\rm DE}^{{\rm (iso)}}$ typically occurs at redshifts $z$ much larger than 1. Thus, apart from the existence of dark radiation, the anisotropic cosmological dynamics at low redshifts is similar to that in isotropic generalized Proca theories. In beyond-generalized Proca theories the only consistent solution to avoid the divergence of a determinant of the dynamical system corresponds to $v=0$, so $\Sigma$ always decreases in time.

We consider a massive vector field with derivative interactions that propagates only the 3 desired polarizations (besides two tensor polarizations from gravity) with second-order equations of motion in curved space-time. The cosmological implications of such generalized Proca theories are investigated for both the background and the linear perturbation by taking into account the Lagrangian up to quintic order. In the presence of a matter fluid with a temporal component of the vector field, we derive the background equations of motion and show the existence of de Sitter solutions relevant to the late-time cosmic acceleration. We also obtain conditions for the absence of ghosts and Laplacian instabilities of tensor, vector, and scalar perturbations in the small-scale limit. Our results are applied to concrete examples of the general functions in the theory, which encompass vector Galileons as a specific case. In such examples, we show that the de Sitter fixed point is always a stable attractor and study viable parameter spaces in which the no-ghost and stability conditions are satisfied during the cosmic expansion history.

In the context of massive (bi-)gravity non-minimal matter couplings have been proposed. These couplings are special in the sense that they are free of the Boulware-Deser ghost below the strong coupling scale and can be used consistently as an effective field theory. Furthermore, they enrich the phenomenology of massive gravity. We consider these couplings in the framework of bimetric gravity and study the cosmological implications for background and linear tensor, vector, and scalar perturbations. Previous works have investigated special branch of solutions. Here we perform a complete perturbation analysis for the general background equations of motion completing previous analysis.

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

For a massive vector field with derivative self-interactions, the breaking of the gauge invariance allows the propagation of a longitudinal mode in addition to the two transverse modes. We consider generalized Proca theories with second-order equations of motion in a curved space-time and study how the longitudinal scalar mode of the vector field gravitates on a spherically symmetric background. We show explicitly that cubic-order self-interactions lead to the suppression of the longitudinal mode through the Vainshtein mechanism. Provided that the dimensionless coupling of the interaction is not negligible, this screening mechanism is sufficiently efficient to give rise to tiny corrections to gravitational potentials consistent with solar-system tests of gravity. We also study the quartic interactions with the presence of non-minimal derivative coupling with the Ricci scalar and find the existence of solutions where the longitudinal mode completely vanishes. Finally, we discuss the case in which the effect of the quartic interactions dominates over the cubic one and show that local gravity constraints can be satisfied under a mild bound on the parameters of the theory.

Massive gravity in the presence of doubly coupled matter field via en effective composite metric yields an accelerated expansion of the universe. It has been recently shown that the model admits stable de Sitter attractor solutions and could be used as a dark energy model. In this work, we perform a first analysis of the constraints imposed by the SNIa, BAO and CMB data on the massive gravity model with the effective composite metric and show that all the background observations are mutually compatible at the one sigma level with the model.

This paper is dedicated to scrutinizing the cosmology in massive gravity. A matter field of the dark sector is coupled to an effective composite metric while a standard matter field couples to the dynamical metric in the usual way. For this purpose, we study the dynamical system of cosmological solutions by using phase analysis, which provides an overview of the class of cosmological solutions in this setup. This also permits us to study the critical points of the cosmological equations together with their stability. We show the presence of stable attractor de Sitter critical points relevant to the late-time cosmic acceleration. Furthermore, we study the tensor, vector and scalar perturbations in the presence of standard matter fields and obtain the conditions for the absence of ghost and gradient instabilities. Hence, massive gravity in the presence of the effective composite metric can accommodate interesting dark energy phenomenology, that can be observationally distinguished from the standard model according to the expansion history and cosmic growth.

In this work we study the cosmological perturbations in massive bigravity in the presence of non-minimal derivative couplings. For this purpose we consider a specific subclass of Horndeski scalar-tensor interactions that live on the unique composite effective metric. For the viability of the model both metrics have to be dynamical. Nevertheless, the number of allowed kinetic terms is crucial. We adapt to the restriction of having one single kinetic term. After deriving the full set of equations of motion for flat Friedmann-Lemaitre-Robertson-Walker background, we study linear perturbations on top of it. We show explicitly that only four tensor, two vector and two scalar degrees of freedom propagate, one of which being the Horndeski scalar, while the Boulware-Deser ghost can be integrated out.

In this work we study the partially constrained vielbein formulation of the new quasidilaton theory of massive gravity which couples to both physical and fiducial metrics simultaneously via a composite effective metric. This formalism improves the new quasidilaton model since the Boulware-Deser ghost is removed fully non-linearly at all scales. This also yields crucial implications in the cosmological applications. We derive the governing cosmological background evolution and study the stability of the attractor solution.

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

A class of vector-tensor theories arises naturally in the framework of quadratic gravity in spacetimes with linear vector distortion. Requiring the absence of ghosts for the vector field imposes an interesting condition on the allowed connections with vector distortion: the resulting one-parameter family of connections generalises the usual Weyl geometry with polar torsion. The cosmology of this class of theories is studied, focusing on isotropic solutions wherein the vector field is dominated by the temporal component. De Sitter attractors are found and inhomogeneous perturbations around such backgrounds are analysed. In particular, further constraints on the models are imposed by excluding pathologies in the scalar, vector and tensor fluctuations. Various exact background solutions are presented, describing a constant and an evolving dark energy, a bounce and a self-tuning de Sitter phase. However, the latter two scenarios are not viable under a closer scrutiny.

Massive gravity theories have been developed as viable IR modifications of gravity motivated by dark energy and the problem of the cosmological constant. On the other hand, modified gravity and modified dark matter theories were developed with the aim of solving the problems of standard cold dark matter at galactic scales. Here we propose to adapt the framework of ghost-free massive bigravity theories to reformulate the problem of dark matter at galactic scales. We investigate a promising alternative to dark matter called dipolar dark matter (DDM) in which two different species of dark matter are separately coupled to the two metrics of bigravity and are linked together by an internal vector field. We show that this model successfully reproduces the phenomenology of dark matter at galactic scales (MOND) as a result of a mechanism of gravitational polarisation. The model is safe in the gravitational sector, but because of the particular couplings of the matter fields and vector field to the metrics, a ghost in the decoupling limit is present in the dark matter sector. However, it might be possible to push the mass of the ghost beyond the strong coupling scale by an appropriate choice of the parameters of the model. Crucial questions to address in future work are the exact mass of the ghost, and the cosmological implications of the model.

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

We study a general class of gravitational theories formulated in the Palatini approach and derive the equations governing the evolution of tensor perturbations. In the absence of torsion, the connection can be solved as the Christoffel symbols of an auxiliary metric which is non-trivially related to the space-time metric. We then consider background solutions corresponding to a perfect fluid and show that the tensor perturbations equations (including anisotropic stresses) for the auxiliary metric around such a background take an Einstein-like form. This facilitates the study in a homogeneous and isotropic cosmological scenario where we explicitly establish the relation between the auxiliary metric and the space-time metric tensor perturbations. As a general result, we show that both tensor perturbations coincide in the absence of anisotropic stresses.

We consider a consistent linear effective vielbein matter coupling without introducing the Boulware-Deser ghost in ghost-free massive gravity. This is achieved in the partially constrained vielbein formulation. We first introduce the formalism and prove the absence of ghost at all scales. As next we investigate the cosmological application of this coupling in this new formulation. We show that even if the background evolution accords with the metric formulation, the perturbations display important different features in the partially constrained vielbein formulation. We study the cosmological perturbations of the two branches of solutions separately. The tensor perturbations coincide with those in the metric formulation. Concerning the vector and scalar perturbations, the requirement of absence of ghost and gradient instabilities yields slightly different allowed parameter space.

In the framework of Born-Infeld inspired gravity theories, which deviates from General Relativity (GR) in the high curvature regime, we discuss the viability of Cosmic Inflation without scalar fields. For energy densities higher than the new mass scale of the theory, a gravitating dust component is shown to generically induce an accelerated expansion of the Universe. Within such a simple scenario, inflation gracefully exits when the GR regime is recovered, but the Universe would remain matter dominated. In order to implement a reheating era after inflation, we then consider inflation to be driven by a mixture of unstable dust species decaying into radiation. Because the speed of sound gravitates within the Born-Infeld model under consideration, our scenario ends up being predictive on various open questions of the inflationary paradigm. The total number of e-folds of acceleration is given by the lifetime of the unstable dust components and is related to the duration of reheating. As a result, inflation does not last much longer than the number of e-folds of deceleration allowing a small spatial curvature and large scale deviations to isotropy to be observable today. Energy densities are self-regulated as inflation can only start for a total energy density less than a threshold value, again related to the species' lifetime. Above this threshold, the Universe may bounce thereby avoiding a singularity. Another distinctive feature is that the accelerated expansion is of the superinflationary kind, namely the first Hubble flow function is negative. We show however that the tensor modes are never excited and the tensor-to-scalar ratio is always vanishing, independently of the energy scale of inflation.

We review some recent proposals for relativistic models of dark matter in the context of bimetric gravity. The aim is to solve the problems of cold dark matter (CDM) at galactic scales, and to reproduce the phenomenology of the modified Newtonian dynamics (MOND), while still being in agreement with the standard cosmological model $\Lambda$-CDM at large scales. In this context a promising alternative is dipolar dark matter (DDM) in which two different species of dark matter particles are separately coupled to the two metrics of bigravity and are linked together by an internal vector field. The phenomenology of MOND then results from a mechanism of gravitational polarization. Probably the best formulation of the model is within the framework of recently developed massive bigravity theories. Then the gravitational sector of the model is safe by construction, but a ghostly degree of freedom in the decoupling limit is still present in the dark matter sector. Future work should analyse the cosmological solutions of the model and check the post-Newtonian parameters in the solar system.

In the context of massive gravity, bi-gravity and multi-gravity non-minimal matter couplings via a specific composite effective metric were investigated recently. Even if these couplings generically reintroduce the Boulware-Deser ghost, this composite metric is unique in the sense that the ghost reemerges only beyond the decoupling limit and the matter quantum loop corrections do not detune the potential interactions. We consider non-minimal {\it derivative} couplings of the composite metric to matter fields for a specific subclass of Horndeski scalar-tensor interactions. We first explore these couplings in the mini-superspace and investigate in which scenario the ghost remains absent. We further study these non-minimal derivative couplings in the decoupling-limit of the theory and show that the equation of motion for the helicity-0 mode remains second order in derivatives. Finally, we discuss preliminary implications for cosmology.

In this work we investigate the existence of relativistic models for dark matter in the context of bimetric gravity, used here to reproduce the modified Newtonian dynamics (MOND) at galactic scales. For this purpose we consider two different species of dark matter particles that separately couple to the two metrics of bigravity. These two sectors are linked together \textit{via} an internal $U(1)$ vector field, and some effective composite metric built out of the two metrics. Among possible models only certain classes of kinetic and interaction terms are allowed without invoking ghost degrees of freedom. Along these lines we explore the number of allowed kinetic terms in the theory and point out the presence of ghosts in a previous model. Finally, we propose a promising class of ghost-free candidate theories that could provide the MOND phenomenology at galactic scales while reproducing the standard cold dark matter (CDM) model at cosmological scales.

In this work we study the theory of extended quasidilaton massive gravity together with the presence of matter fields. After discussing the homogeneous and isotropic fully dynamical background equations, which governs the exact expansion history of the universe, we consider small cosmological perturbations around these general FLRW solutions. The stability of tensor, vector and scalar perturbations on top of these general background solutions give rise to slightly different constraints on the parameters of the theory than those obtained in the approximative assumption of the late-time asymptotic form of the expansion history, which does not correspond to our current epoch. This opens up the possibility of stable FLRW solutions to be compared with current data on cosmic expansion with the restricted parameter space based on theoretical ground.

We study the cosmology of bimetric theory with a composite matter coupling. We find two possible branches of background evolution. We investigate the question of stability of cosmological perturbations. For the tensor and vector perturbations, we derive conditions on the absence of ghost and gradient instabilities. For the scalar modes, we obtain conditions for avoiding ghost degrees. In the first branch, we find that one of the scalar modes becomes a ghost at the late stages of the evolution. Conversely, this problem can be avoided in the second branch. However, we also find that the constraint for the second branch prevents the doubly coupled matter fields from being the standard ingredients of cosmology. We thus conclude that a realistic and stable cosmological model requires additional minimally coupled matter fields.

We investigate the weak lensing signature of primordial non-Gaussianities of the local type by constraining the magnitude of the weak convergence bi- and trispectra expected for the EUCLID weak lensing survey. Starting from expressions for the weak convergence spectra, bispectra and trispectra, whose relative magnitudes we investigate as a function of scale, we compute their respective signal to noise ratios by relating the polyspectra's amplitude to their Gaussian covariance using a Monte-Carlo technique for carrying out the configuration space integrations. In computing the Fisher-matrix on the non-Gaussianity parameters f_nl, g_nl and tau_nl with a very similar technique, we can derive Bayesian evidences for a violation of the Suyama-Yamaguchi relation tau_nl>=(6 f_nl/5)^2 as a function of the true f_nl and tau_nl-values and show that the relation can be probed down to levels of f_nl~10^2 and tau_nl~10^5. In a related study, we derive analytical expressions for the probability density that the SY-relation is exactly fulfilled, as required by models in which any one field generates the perturbations. We conclude with an outlook on the levels of non-Gaussianity that can be probed with tomographic lensing surveys.

We generalize the ultraviolet sector of gravitation via a Born-Infeld action using lessons from massive gravity. The theory contains all of the elementary symmetric polynomials and is treated in the Palatini formalism. We show how the connection can be solved algebraically to be the Levi-Civita connection of an effective metric. The non-linearity of the algebraic equations yields several branches, one of which always reduces to General Relativity at low curvatures. We explore in detail a {\it minimal} version of the theory, for which we study solutions in the presence of a perfect fluid with special attention to the cosmological evolution. In vacuum we recover Ricci-flat solutions, but also an additional physical solution corresponding to an Einstein space. The existence of two physical branches remains for non-vacuum solutions and, in addition, the branch that connects to the Einstein space in vacuum is not very sensitive to the specific value of the energy density. For the branch that connects to the General Relativity limit we generically find three behaviours for the Hubble function depending on the equation of state of the fluid, namely: either there is a maximum value for the energy density that connects continuously with vacuum, or the energy density can be arbitrarily large but the Hubble function saturates and remains constant at high energy densities, or the energy density is unbounded and the Hubble function grows faster than in General Relativity. The second case is particularly interesting because it could offer an interesting inflationary epoch even in the presence of a dust component. Finally, we discuss the possibility of avoiding certain types of singularities within the minimal model.

We investigate the coupling to matter in ghost-free massive (bi-)gravity. When species in the matter sector couple covariantly to only one metric, we show that at one--loop these couplings do not spoil the special structure of the graviton potential. When the same species couples directly to both metrics we show that a ghost is present at the classical level and that loops destroy the special structure of the potential at an unacceptably low scale. We then propose a new `composite' effective metric built out of both metrics. When matter fields couple covariantly to this effective metric, the would be Boulware--Deser ghost is absent in different representative limits. At one--loop such couplings do not detune the special structure of the potential. We conjecture that matter can couple covariantly to that effective metric in all generality without introducing any Boulware--Deser ghost below a cut-off scale parametrically larger than the strong coupling scale. We also discuss alternative couplings to matter where the kinetic and potential terms of the matter field couple to different metrics. In both cases we discuss preliminary implications for cosmology.

We compute the one-loop quantum corrections to the interactions between the two metrics of the ghost-free massive bigravity. When considering gravitons running in the loops, we show how the structure of the interactions gets destabilized at the quantum level, exactly in the same way as in its massive gravity limit. A priori one might have expected a better quantum behavior, however the broken diffeomorphism invariance out of the two initial diffeomorphisms in bigravity has similar consequences at the quantum level as the broken diffeomorphism in massive gravity. From lessons of the generated quantum corrections through matter loops we propose yet other types of effective composite metrics to which the matter fields can couple. Among these new effective metrics there might be one or more that could provide interesting phenomenology and important cosmological implications.

We investigate the cosmological perturbations around FLRW solutions to non- linear massive gravity with a new effective coupling to matter proposed recently. Unlike the case with minimal matter coupling, all five degrees of freedom in the gravity sector propagate on generic self-accelerating FLRW backgrounds. We study the stability of the cosmological solutions and put constraints on the parameters of the theory by demanding the correct sign for the kinetic terms for scalar, vector and tensor perturbations.

Recently, several works have investigated the coupling to matter in ghost-free massive (bi-&multi-)gravity and a new effective coupling to matter has been proposed. In this note we clarify some confusion on the existence and the implications of a ghost above the strong coupling scale. We confirm that the standard constraint which is otherwise typically present in this type of theories disappears on generic backgrounds as soon as this new coupling is considered. This implies the re-emergence of the Boulware-Deser ghost. Nevertheless the absence of ghost in the decoupling limit implies that the cut-off scale (if identified with the scale at which the ghost enters) is higher than the strong coupling scale. Therefore there is a valid interesting region of applicability for these couplings at scales below the cut-off.

In this paper, we scrutinize very closely the cosmology in the proxy theory to massive gravity obtained in Phys. Rev. D84 (2011) 043503. This proxy theory was constructed by covariantizing the decoupling limit Lagrangian of massive gravity and represents a subclass of Horndeski scalar-tensor theory. Thus, this covariantization unifies two important classes of modified gravity theories, namely massive gravity and Horndeski theories. We go beyond the regime which was studied in Phys. Rev. D84 (2011) 043503 and show that the theory does not admit any homogeneous and isotropic self-accelerated solutions. We illustrate that the only attractor solution is flat Minkowski solution, hence this theory is less appealing as a dark energy model. We also show that the absence of de Sitter solutions is tightly related to the presence of shift symmetry breaking interactions.

We study the Horndeski vector-tensor theory that leads to second order equations of motion and contains a non-minimally coupled abelian gauge vector field. This theory is remarkably simple and consists of only 2 terms for the vector field, namely: the standard Maxwell kinetic term and a coupling to the dual Riemann tensor. Furthermore, the vector sector respects the U(1) gauge symmetry and the theory contains only one free parameter, M^2, that controls the strength of the non-minimal coupling. We explore the theory in a de Sitter spacetime and study the presence of instabilities and show that it corresponds to an attractor solution in the presence of the vector field. We also investigate the cosmological evolution and stability of perturbations in a general FLRW spacetime. We find that a sufficient condition for the absence of ghosts is M^2>0. Moreover, we study further constraints coming from imposing the absence of Laplacian instabilities. Finally, we study the stability of the theory in static and spherically symmetric backgrounds (in particular, Schwarzschild and Reissner-Nordstr\"om-de Sitter). We find that the theory, quite generally, do have ghosts or Laplacian instabilities in regions of spacetime where the non-minimal interaction dominates over the Maxwell term. We also calculate the propagation speed in these spacetimes and show that superluminality is a quite generic phenomenon in this theory.

We compute the one-loop quantum corrections to the potential of ghost-free massive gravity. We show how the mass of external matter fields contribute to the running of the cosmological constant, but do not change the ghost-free structure of the massive gravity potential at one-loop. When considering gravitons running in the loops, we show how the structure of the potential gets destabilized at the quantum level, but in a way which would never involve a ghost with a mass smaller than the Planck scale. This is done by explicitly computing the one-loop effective action and supplementing it with the Vainshtein mechanism. We conclude that to one-loop order the special mass structure of ghost-free massive gravity is technically natural.

We re-explore the Bi- and Multi-Galileon models with trivial asymptotic conditions at infinity and show that propagation of superluminal fluctuations is a common and unavoidable feature of these theories, unlike previously claimed in the literature. We show that all Multi-Galileon theories containing a Cubic Galileon term exhibit superluminalities at large distances from a point source, and that even if the Cubic Galileon is not present one can always find sensible matter distributions in which there are superluminal modes at large distances. In the Bi-Galileon case we explicitly show that there are always superluminal modes around a point source even if the Cubic Galileon is not present. Finally, we briefly comment on the possibility of avoiding superluminalities by modifying the asymptotic conditions at infinity.

Subject of this paper is a detailed analysis of the PINOCCHIO algorithm for studying the relative velocity statistics of merging haloes in Lagrangian perturbation theory. Given a cosmological background model, a power spectrum of fluctuations as well as a Gaussian linear density contrast field $\delta_{\rm l}$ is generated on a cubic grid, which is then smoothed repeatedly with Gaussian filters. For each Lagrangian particle at position $\bmath{q}$ and each smoothing radius $R$, the collapse time, the velocities and ellipsoidal truncation are computed using Lagrangian Perturbation Theory. The collapsed medium is then fragmented into isolated objects by an algorithm designed to mimic the accretion and merger events of hierarchical collapse. Directly after the fragmentation process the mass function, merger histories of haloes and the statistics of the relative velocities at merging are evaluated. We reimplemented the algorithm in C++, recovered the mass function and optimised the construction of halo merging histories. Comparing our results with the output of the Millennium simulation suggests that PINOCCHIO is well suited for studying relative velocities of merging haloes and is able to reproduce the pairwise velocity distribution.

We covariantize the decoupling limit of massive gravity proposed in arXiv:1011.1232 and study the cosmology of this theory as a proxy, which embodies key features of the fully non-linear covariant theory. We first confirm that it exhibits a self-accelerating solution, similar to what has been found in arXiv:1010.1780, where the Hubble parameter corresponds to the graviton mass. For a certain range of parameters fluctuations relative to the self-accelerating background are stable and form an attractor solution. We also show that a degravitating solution can not be constructed in this covariantized proxy theory in a meaningful way. As for cosmic structure formation, we find that the helicity-0 mode of the graviton causes an enhancement relative to LCDM. For consistency we also compare proxy theories obtained starting from different frames in the decoupling limit and discuss the possibility of obtaining a non-representative proxy theory by choosing the wrong starting frame.

Galileon models are a class of effective field theories that have recently received much attention. They arise in the decoupling limit of theories of massive gravity, and in some cases they have been treated in their own right as scalar field theories with a specific nonlinearly realized global symmetry (Galilean transformation). It is well known that in the presence of a source, these Galileon theories admit superluminal propagating solutions, implying that as quantum field theories they must admit a different notion of causality than standard local Lorentz invariant theories. We show that in these theories it is easy to construct closed timelike curves (CTCs) within the {\it naive} regime of validity of the effective field theory. However, on closer inspection we see that the CTCs could never arise since the Galileon inevitably becomes infinitely strongly coupled at the onset of the formation of a CTC. This implies an infinite amount of backreaction, first on the background for the Galileon field, signaling the break down of the effective field theory, and subsequently on the spacetime geometry, forbidding the formation of the CTC. Furthermore the background solution required to create CTCs becomes unstable with an arbitrarily fast decay time. Thus Galileon theories satisfy a direct analogue of Hawking's chronology protection conjecture.

The topic of this article is weak cosmic shear tomography where the line of sight-weighting is carried out with a set of specifically constructed orthogonal polynomials, dubbed TaRDiS (Tomography with orthogonAl Radial Distance polynomIal Systems). We investigate the properties of these polynomials and employ weak convergence spectra, which have been obtained by weighting with these polynomials, for the estimation of cosmological parameters. We quantify their power in constraining parameters in a Fisher-matrix technique and demonstrate how each polynomial projects out statistically independent information, and how the combination of multiple polynomials lifts degeneracies. The assumption of a reference cosmology is needed for the construction of the polynomials, and as a last point we investigate how errors in the construction with a wrong cosmological model propagate to misestimates in cosmological parameters. TaRDiS performs on a similar level as traditional tomographic methods and some key features of tomography are made easier to understand.

Weak lensing calculations are often made under the assumption of the Born approximation, where the ray path is approximated as a straight radial line. In addition, lens-lens couplings where there are several deflections along the light ray are often neglected. We examine the effect of dropping the Born approximation and taking lens-lens couplings into account, for weak lensing effects up to second order (cosmic flexion), by making a perturbative expansion in the light path. We present a diagrammatic representation of the resulting corrections to the lensing effects. The flexion signal, which measures the derivative of the density field, acquires correction terms proportional to the squared gravitational shear; we also find that by dropping the Born approximation, two further degrees of freedom of the lensing distortion can be excited (the twist components), in addition to the four standard flexion components. We derive angular power spectra of the flexion and twist, with and without the Born-approximation and lens-lens couplings and confirm that the Born approximation is an excellent approximation for weak cosmic flexions, except at very small scales.

We explore cosmology in the decoupling limit of a non-linear covariant extension of Fierz-Pauli massive gravity obtained recently in arXiv:1007.0443. In this limit the theory is a scalar-tensor model of a unique form defined by symmetries. We find that it admits a self-accelerated solution, with the Hubble parameter set by the graviton mass. The negative pressure causing the acceleration is due to a condensate of the helicity-0 component of the massive graviton, and the background evolution, in the approximation used, is indistinguishable from the \Lambda CDM model. Fluctuations about the self-accelerated background are stable for a certain range of parameters involved. Most surprisingly, the fluctuation of the helicity-0 field above its background decouples from an arbitrary source in the linearized theory. We also show how massive gravity can remarkably screen an arbitrarily large cosmological constant in the decoupling limit, while evading issues with ghosts. The obtained static solution is stable against small perturbations, suggesting that the degravitation of the vacuum energy is possible in the full theory. Interestingly, however, this mechanism postpones the Vainshtein effect to shorter distance scales. Hence, fifth force measurements severely constrain the value of the cosmological constant that can be neutralized, making this scheme phenomenologically not viable for solving the old cosmological constant problem. We briefly speculate on a possible way out of this issue.

The subject of this paper is an investigation of the nonlinear contributions to the spectrum of the integrated Sachs-Wolfe (iSW) effect. We derive the corrections to the iSW-auto spectrum and the iSW-tracer cross-spectrum consistently to third order in perturbation theory and analyse the cumulative signal-to-noise ratio for a cross-correlation between the PLANCK and EUCLID data sets as a function of multipole order. We quantify the parameter sensitivity and the statistical error bounds on the cosmological parameters Omega_m, sigma_8, h, n_s and w from the linear iSW-effect and the systematical parameter estimation bias due to the nonlinear corrections in a Fisher-formalism, analysing the error budget in its dependence on multipole order. Our results include: (i) the spectrum of the nonlinear iSW-effect can be measured with 0.8\sigma statistical significance, (ii) nonlinear corrections dominate the spectrum starting from l=100, (iii) an anticorrelation of the CMB temperature with tracer density on high multipoles in the nonlinear regime, (iv) a much weaker dependence of the nonlinear effect on the dark energy model compared to the linear iSW-effect, (v) parameter estimation biases amount to less than 0.1 sigma and weaker than other systematics.