We perform a reconstruction of the coupling function between vacuum energy and geodesic cold dark matter using the latest observational data. We bin the interaction in seventeen redshift bins but use a correlation prior to prevent rapid, unphysical oscillations in the coupling function. This prior also serves to eliminate any dependence of the reconstruction on the binning method. We use two different forms of the correlation prior, finding that both give similar results for the reconstruction of the dark matter -- dark energy interaction. Calculating the Bayes factor for each case, we find no meaningful evidence for deviation from the null interacting case, i.e. $\Lambda$CDM, in our reconstruction.

We investigate the generalized cubic covariant Galileon model, a kinetically driven dark energy model within the Horndeski class of theories. The model extends the cubic covariant Galileon by including power laws of the field derivatives in the K-essence and cubic terms which still allow for tracker solutions. We study the shape of the viable parameter space by enforcing stability conditions which include the absence of ghost, gradient and tachyon instabilities and the avoidance of strong coupling at early time. We study here the relevant effects of the modifications induced by the model on some cosmological observables such as the cosmic microwave background (CMB), the lensing potential auto-correlation and the matter power spectrum. For this goal, we perform parameter estimation using data of CMB temperature and polarization, baryonic acoustic oscillations (BAO), redshift-space distortions (RSD), supernovae type Ia (SNIa) and Cepheids. Data analysis with CMB alone finds that the today's Hubble parameter $H_0$ is consistent with its determination from Cepheids at $1\sigma$, resolving the famous tension of the cosmological standard models. The joint analysis of CMB, BAO, RSD and SNIa sets a lower bound for the sum of neutrino masses which is $\Sigma m_\nu >0.11$ eV at 1$\sigma$, in addition to the usual upper limit. The model selection analysis based on the effective $\chi_\text{eff}^2$ and Deviance Information Criterion is not able to clearly identify the statistically favored model between $\Lambda$CDM and the generalized cubic covariant Galileon, from which we conclude that the latter model deserves further studies.

Measurements of time delays between multiple quasar images produced by strong lensing are reaching a sensitivity that makes them a promising cosmological probe. Future surveys will provide significantly more measurements, reaching unprecedented depth in redshift, making strong lensing time delay (SLTD) observations competitive with other background probes. We forecast constraints on the nature of dark energy from upcoming SLTD surveys, simulating future catalogues with different numbers of lenses distributed up to redshift $z\sim 1$ and focusing on cosmological parameters such as the Hubble constant $H_0$ and parametrisations of the dark energy equation of state. We also explore the impact of our ability to precisely model the lens mass profile and its environment, on the forecasted constraints. We find that in the most optimistic cases, SLTD will constrain $H_0$ at the level of $\sim 0.1\%$, while the CPL equation of state parameters, $w_0$ and $w_a$, can be determined with errors $\sigma_{w_0}\sim 0.05$ and $\sigma_{w_a}\sim 0.3$, respectively. Furthermore, we investigate the bias introduced when a wrong cosmological model is assumed for the analysis. We find that the value of $H_0$ could be biased up to $10 \sigma$, assuming a perfect knowledge of the lens profile, when a $\Lambda$CDM model is used to analyse data that really belong to a $w$CDM cosmology with $w=-0.9$. Based on these findings, we identify a consistency check of the assumed cosmological model in future SLTD surveys, by splitting the dataset in several redshift bins. Depending on the characteristics of the survey, this could provide a smoking gun for dark energy.

We place observational constraints on the Galileon ghost condensate model, a dark energy proposal in cubic-order Horndeski theories consistent with the gravitational-wave event GW170817. The model extends the covariant Galileon by taking an additional higher-order field derivative $X^2$ into account. This allows for the dark energy equation of state $w_{\rm DE}$ to access the region $-2<w_{\rm DE}<-1$ without ghosts. Indeed, this peculiar evolution of $w_{\rm DE}$ is favored over that of the cosmological constant $\Lambda$ from the joint data analysis of cosmic microwave background (CMB) radiation, baryonic acoustic oscillations (BAOs), supernovae type Ia (SNIa) and redshift-space distortions (RSDs). Furthermore, our model exhibits a better compatibility with the CMB data over the $\Lambda$-cold-dark-matter ($\Lambda$CDM) model by suppressing large-scale temperature anisotropies. The CMB temperature and polarization data lead to an estimation for today's Hubble parameter $H_0$ consistent with its direct measurements at 2$\sigma$. We perform a model selection analysis by using several methods and find a statistically significant preference of the Galileon ghost condensate model over $\Lambda$CDM.

We study observational constraints on a specific dark energy model in the framework of Gleyzes-Langlois-Piazza-Vernizzi theories, which extends the Galileon ghost condensate (GGC) to the domain of beyond Horndeski theories. In this model, we show that the Planck cosmic microwave background (CMB) data, combined with datasets of baryon acoustic oscillations, supernovae type Ia, and redshift-space distortions, give the tight upper bound $|\alpha_{\rm H}^{(0)}| \le {\cal O}(10^{-6})$ on today's beyond-Horndeski (BH) parameter $\alpha_{\rm H}$. This is mostly attributed to the shift of CMB acoustic peaks induced by the early-time changes of cosmological background and perturbations arising from the dominance of $\alpha_{\rm H}$ in the dark energy density. In comparison to the $\Lambda$-cold-dark-matter ($\Lambda$CDM) model, our BH model suppresses the large-scale integrated-Sachs-Wolfe (ISW) tail of CMB temperature anisotropies due to the existence of cubic Galileons, and it modifies the small-scale CMB power spectrum because of the different background evolution. We find that the BH model considered fits the data better than $\Lambda$CDM according to the $\chi^2$ statistics, yet the deviance information criterion (DIC) slightly favors the latter. Given the fact that our BH model with $\alpha_{\rm H}=0$ (i.e., the GGC model) is favored over $\Lambda$CDM even by the DIC, there are no particular signatures for the departure from Horndeski theories in current observations.

We investigate an interacting dark sector scenario in which the vacuum energy is free to interact with cold dark matter (CDM), which itself is assumed to cluster under the sole action of gravity, i.e. it is in free fall (geodesic), as in $\Lambda$CDM. The interaction is characterised by a dimensionless coupling $q_{\rm V}$ that we constrain using cosmic microwave background data from the Planck 2015 data release, along with baryon acoustic oscillation, redshift space distortion and Type Ia supernova measurements. We present the full linear perturbation theory of this interacting scenario and use MCMC sampling to study five different cases: two cases in which we have $\Lambda$CDM evolution in the distant past, until a set redshift $z_{\rm trans}$, below which the interaction switches on and $q_{\rm V}$ is the single sampled parameter, with $z_{\rm trans}$ fixed at $z_{\rm trans}=3000$ and $z_{\rm trans}=0.9$ respectively; a case where we allow this transition redshift to vary along with $q_{\rm V}$; a case in which the vacuum energy is zero for $z>z_{\rm trans}$ and then begins to grow once the interaction switches on; and the final case in which we bin $q_{\rm V}(z)$ in four redshift bins to investigate the possibility of a dynamical interaction, reconstructing the redshift evolution of the function using Gaussian processes. We find that, in all cases where the high redshift evolution is not modified, the results are compatible with a vanishing coupling, thus finding no significant deviation from $\Lambda$CDM.

The tachyonic instability is associated with the unboundedness of the Hamiltonian from below and results in an unstable low-$k$ regime. In the cosmological exploration of modified gravity, it is seldom taken into account, with more focus given to the popular no-ghost and no-gradient conditions. The latter though are intrinsically high-$k$ statements. Here we combine all three conditions into a full set of requirements that we show to guarantee stability on the whole range of cosmological scales. We then explore the impact of the different conditions on the parameter space of scalar-tensor gravity, with particular emphasis on the no-tachyon one. We focus on Horndeski gravity and also consider separately the two subclasses of $f(R)$ and Generalized Brans Dicke theories. We identify several interesting features, for instance in the parameter space of designer $f(R)$ on a $w$CDM background, shedding light on previous findings. When looking at the phenomenological functions $\Sigma$ and $\mu$, associated to the weak lensing and clustering potential respectively, we find that in the case of Generalized Brans Dicke the no-tachyon condition clearly cuts models with $\mu\,,\,\Sigma>1$. This effect is less prevalent in the Horndeski case due to the larger amount of free functions in the theory.

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, $w_{\rm DE}$, and the phenomenology of Large Scale Structure. They will allow us to constrain deviations from the $\Lambda$CDM predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions $\Sigma$ and $\mu$, respectively. In this work, we derive the theoretical prior for the joint covariance of $w_{\rm DE}$, $\Sigma$ and $\mu$, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte-Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. We identify several interesting features and trends in the distribution functions of $w_{\rm DE}$, $\Sigma$ and $\mu$, as well as in their covariances; we confirm the high degree of correlation between $\Sigma$ and $\mu$ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly $w_{\rm DE}$, $\Sigma$ and $\mu$ in a non-parametric way.

In the context of the effective field theory of dark energy (EFT) we perform agnostic explorations of Horndeski gravity. We choose two parametrizations for the free EFT functions, namely a power law and a dark energy density-like behaviour on a non trivial Chevallier-Polarski-Linder background. We restrict our analysis to those EFT functions which do not modify the speed of propagation of gravitational waves. Among those, we prove that one specific function cannot be constrained by data, since its contribution to the observables is below the cosmic variance, although we show it has a relevant role in defining the viable parameter space. We place constraints on the parameters of these models combining measurements from present day cosmological datasets and we prove that the next generation galaxy surveys can improve such constraints by one order of magnitude. We then verify the validity of the quasi-static limit within the sound horizon of the dark field, by looking at the phenomenological functions $\mu$ and $\Sigma$, associated respectively to clustering and lensing potentials. Furthermore, we notice up to $5\%$ deviations in $\mu, \Sigma$ with respect to General Relativity at scales smaller than the Compton one. For the chosen parametrizations and in the quasi-static limit, future constraints on $\mu$ and $\Sigma$ can reach the $1\%$ level and will allow us to discriminate between certain models at more than $3\sigma$, provided the present best-fit values remain.

We revisit the cosmology of Covariant Galileon gravity in view of the most recent cosmological data sets, including weak lensing. As a higher derivative theory, Covariant Galileon models do not have a $\Lambda$CDM limit and predict a very different structure formation pattern compared with the standard $\Lambda$CDM scenario. Previous cosmological analyses suggest that this model is marginally disfavoured, yet can not be completely ruled out. In this work we use a more recent and extended combination of data, and we allow for more freedom in the cosmology, by including a massive neutrino sector with three different mass hierarchies. We use the Planck measurements of Cosmic Microwave Background temperature and polarization; Baryonic Acoustic Oscillations measurements by BOSS DR12; local measurements of $H_0$; the joint light-curve analysis supernovae sample; and, for the first time, weak gravitational lensing from the KiDS collaboration. We find, that in order to provide a reasonable fit, a non-zero neutrino mass is indeed necessary, but we do not report any sizable difference among the three neutrino hierarchies. Finally, the comparison of the Bayesian Evidence to the $\Lambda$CDM one shows that in all the cases considered, Covariant Galileon models are statistically ruled out by cosmological data.

Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize modifications of the growth of large-scale structure in alternative theories of gravity. We study the values these functions can take in Horndeski theories, i.e. the class of scalar-tensor theories with second order equations of motion. We restrict our attention to models that are in a broad agreement with tests of gravity and the observed cosmic expansion history. In particular, we require the speed of gravity to be equal to the speed of light today, as required by the recent detection of gravitational waves and electromagnetic emission from a binary neutron star merger. We examine the correlations between the values of $\Sigma$ and $\mu$ analytically within the quasi-static approximation, and numerically, by sampling the space of allowed solutions. We confirm that the conjecture made in [Pogosian:2016pwr], that $(\Sigma-1)(\mu -1) \ge 0$ in viable Horndeski theories, holds very well. Along with that, we check the validity of the quasi-static approximation within different corners of Horndeski theory. Our results show that, even with the tight bound on the present day speed of gravitational waves, there is room within Horndeski theories for non-trivial signatures of modified gravity at the level of linear perturbations.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

We investigate the impact of general conditions of theoretical stability and cosmological viability on dynamical dark energy models. As a powerful example, we study whether minimally coupled, single field Quintessence models that are safe from ghost instabilities, can source the CPL expansion history recently shown to be mildly favored by a combination of CMB (Planck) and Weak Lensing (KiDS) data. We find that in their most conservative form, the theoretical conditions impact the analysis in such a way that smooth single field Quintessence becomes significantly disfavored with respect to the standard LCDM cosmological model. This is due to the fact that these conditions cut a significant portion of the (w0;wa) parameter space for CPL, in particular eliminating the region that would be favored by weak lensing data. Within the scenario of a smooth dynamical dark energy parametrized with CPL, weak lensing data favors a region that would require multiple fields to ensure gravitational stability.

Clusters of galaxies have the potential of providing powerful constraints on possible deviations from General Relativity. We use the catalogue of Sunyaev-Zel'dovich sources detected by Planck and consider a correction to the halo mass function for a f(R) class of modified gravity models, which has been recently found to reproduce well results from N-body simulations, to place constraints on the scalaron field amplitude at the present time, $f_{R}^0$. We find that applying this correction to different calibrations of the halo mass function produces upper bounds on $f_{R}^0$ tighter by more than an order of magnitude, ranging from $\log_{10}(-f_{R}^0) < -5.81$ to $\log_{10}(-f_{R}^0) < -4.40$ (95 % confidence level). This sensitivity is due to the different shape of the halo mass function, which is degenerate with the parameters used to calibrate the scaling relations between SZ observables and cluster masses. Any claim of constraints more stringent that the weaker limit above, based on cluster number counts, appear to be premature and must be supported by a careful calibration of the halo mass function and by a robust calibration of the mass scaling relations.