The discovery of cosmic acceleration has triggered a consistent body of theoretical work aimed at modeling its phenomenology and understanding its fundamental physical nature. In recent years, a powerful formalism that accomplishes both these goals has been developed, the so-called effective field theory of dark energy. It can capture the behavior of a wide class of modified gravity theories and classify them according to the imprints they leave on the smooth background expansion history of the Universe and on the evolution of linear perturbations. The effective field theory of dark energy is based on a Lagrangian description of cosmological perturbations which depends on a number of functions of time, some of which are non-minimal couplings representing genuine deviations from General Relativity. Such a formalism is thus particularly convenient to fit and interpret the wealth of new data that will be provided by future galaxy surveys. Despite its recent appearance, this formalism has already allowed a systematic investigation of what lies beyond the General Relativity landscape and provided a conspicuous amount of theoretical predictions and observational results. In this review, we report on these achievements.

We revise the cosmological bounds on Ho\v{r}ava gravity taking into accounts the stringent constraint on the speed of propagation of gravitational waves from GW170817 and GRB170817A. In light of this we also investigate the degeneracy between massive neutrinos and Ho\v{r}ava gravity. We show that a luminal propagation of gravitational waves suppresses the large-scale Cosmic Microwave Background (CMB) radiation temperature anisotropies and the presence of massive neutrinos increases this effect. On the contrary large neutrinos mass can compensate the modifications induced by Ho\v{r}ava gravity in the lensing, matter and primordial B-mode power spectra. Another degeneracy is found, at theoretical level, between the tensor-to-scalar ratio $r$ and massive neutrinos as well as with the model's parameters. We analyze these effects using CMB, supernovae type Ia (SNIa), galaxy clustering and weak gravitational lensing measurements and we show how such degeneracies are removed. We find that the model's parameters are constrained to be very close to their General Relativity limits and we get a two order of magnitude improved upper bound, with respect to the Big Bang Nucleosynthesis constraint, on the deviation of the effective gravitational constant from the Newtonian one. The deviance information criterion suggests that in Ho\v{r}ava gravity $\Sigma m_\nu>0$ is favored when CMB data only are considered, while the joint analysis of all datasets prefers zero neutrinos mass.

We present a detailed study of the collapse of a spherical matter overdensity and the non-linear growth of large scale structures in the Galileon ghost condensate (GGC) model. This model is an extension of the cubic covariant Galileon (G3) which includes a field derivative of type $(\nabla_\mu\phi\nabla^\mu\phi)^2$ in the Lagrangian. We find that the cubic term activates the modifications in the main physical quantities and the additional term has an extended impact in the evolution of the latter. The GGC model shows largely mitigated effects in the linearised critical density contrast, non-linear effective gravitational coupling and the virial overdensity with respect to G3 but still preserves peculiar features with respect to the standard $\Lambda$CDM cosmological model, e.g., both the linear critical density contrast and the virial overdensity are larger than the ones in $\Lambda$CDM. The results of the spherical collapse model are then used to predict the evolution of the halo mass function, non-linear matter and lensing power spectra. While at low masses the GGC model presents about 10% fewer objects with respect to $\Lambda$CDM, at higher masses for $z>0$ it predicts 10% ($z=0.5$)-20% ($z=1$) more objects per comoving volume. Using a phenomenological approach to include the screening effect in the matter power spectrum, we show that the difference induced by the modifications of gravity are strongly dependent on the screening scale and that differences can be up to 60% with respect to $\Lambda$CDM. These differences translate to the lensing power spectrum where qualitatively the largest differences with respect to the standard cosmological model are for $\ell<10^3$. Depending on the screening scale, they can reach 60% on larger angular scales and then decrease for growing $\ell$. The results are obtained for the best fit parameters from cosmological data for each model.

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.

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.

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.

In quadratic-order degenerate higher-order scalar-tensor (DHOST) theories compatible with gravitational-wave constraints, we derive the most general Lagrangian allowing for tracker solutions characterized by $\dot{\phi}/H^p={\rm constant}$, where $\dot{\phi}$ is the time derivative of a scalar field $\phi$, $H$ is the Hubble expansion rate, and $p$ is a constant. While the tracker is present up to the cubic-order Horndeski Lagrangian $L=c_2X-c_3X^{(p-1)/(2p)} \square \phi$, where $c_2, c_3$ are constants and $X$ is the kinetic energy of $\phi$, the DHOST interaction breaks this structure for $p \neq 1$. Even in the latter case, however, there exists an approximate tracker solution in the early cosmological epoch with the nearly constant field equation of state $w_{\phi}=-1-2p\dot{H}/(3H^2)$. The scaling solution, which corresponds to $p=1$, is the unique case in which all the terms in the field density $\rho_{\phi}$ and the pressure $P_{\phi}$ obey the scaling relation $\rho_{\phi} \propto P_{\phi} \propto H^2$. Extending the analysis to the coupled DHOST theories with the field-dependent coupling $Q(\phi)$ between the scalar field and matter, we show that the scaling solution exists for $Q(\phi)=1/(\mu_1 \phi+\mu_2)$, where $\mu_1$ and $\mu_2$ are constants. For the constant $Q$, i.e., $\mu_1=0$, we derive fixed points of the dynamical system by using the general Lagrangian with scaling solutions. This result can be applied to the model construction of late-time cosmic acceleration preceded by the scaling $\phi$-matter-dominated epoch.

In cubic-order Horndeski theories where a scalar field $\phi$ is coupled to nonrelativistic matter with a field-dependent coupling $Q(\phi)$, we derive the most general Lagrangian having scaling solutions on the isotropic and homogenous cosmological background. For constant $Q$ including the case of vanishing coupling, the corresponding Lagrangian reduces to the form $L=Xg_2(Y)-g_3(Y)\square \phi$, where $X=-\partial_{\mu}\phi\partial^{\mu}\phi/2$ and $g_2, g_3$ are arbitrary functions of $Y=Xe^{\lambda \phi}$ with constant $\lambda$. We obtain the fixed points of the scaling Lagrangian for constant $Q$ and show that the $\phi$-matter-dominated-epoch ($\phi$MDE) is present for the cubic coupling $g_3(Y)$ containing inverse power-law functions of $Y$. The stability analysis around the fixed points indicates that the $\phi$MDE can be followed by a stable critical point responsible for the cosmic acceleration. We propose a concrete dark energy model allowing for such a cosmological sequence and show that the ghost and Laplacian instabilities can be avoided even in the presence of the cubic coupling.

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 propose a viable dark energy scenario in the presence of cubic Horndeski interactions and a standard scalar-field kinetic term with two exponential potentials. We show the existence of new scaling solutions along which the cubic coupling $G_3$ provides an important contribution to the field density that scales in the same way as the background fluid density. The solutions finally exit to the epoch of cosmic acceleration driven by a scalar-field dominated fixed point arising from the second exponential potential. We clarify the viable parameter space in which all the theoretically consistent conditions including those for the absence of ghost and Laplacian instabilities are satisfied on scaling and scalar-field dominated critical points. In comparison to Quintessence with the same scalar potential, we find that the cubic coupling gives rise to some novel features: (i) the allowed model parameter space is wider in that a steeper potential can drive the cosmic acceleration; (ii) the dark energy equation of state $w_{\phi}$ today can be closer to $-1$ relative to Quintessence; (iii) even if the density associated with the cubic coupling dominates over the standard field density in the scaling era, the former contribution tends to be suppressed at low redshifts. We also compute quantities associated with the growth of matter perturbations and weak lensing potentials under the quasi-static approximation in the sub-horizon limit and show that the cubic coupling leads to the modified evolution of perturbations which can be distinguished from Quintessence.

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.

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 present the extension of the effective field theory framework to the mildly non-linear scales. The effective field theory approach has been successfully applied to the late time cosmic acceleration phenomenon and it has been shown to be a powerful method to obtain predictions about cosmological observables on linear scales. However, mildly non-linear scales need to be consistently considered when testing gravity theories because a large part of the data comes from those scales. Thus, non-linear corrections to predictions on observables coming from the linear analysis can help in discriminating among different gravity theories. We proceed firstly by identifying the necessary operators which need to be included in the effective field theory Lagrangian in order to go beyond the linear order in perturbations and then we construct the corresponding non-linear action. Moreover, we present the complete recipe to map any single field dark energy and modified gravity models into the non-linear effective field theory framework by considering a general action in the Arnowitt-Deser-Misner formalism. In order to illustrate this recipe we proceed to map the beyond-Horndeski theory and low-energy Horava gravity into the effective field theory formalism. As a final step we derived the 4th order action in term of the curvature perturbation. This allowed us to identify the non-linear contributions coming from the linear order perturbations which at the next order act like source terms. Moreover, we confirm that the stability requirements, ensuring the positivity of the kinetic term and the speed of propagation for scalar mode, are automatically satisfied once the viability of the theory is demanded at linear level. The approach we present here will allow to construct, in a model independent way, all the relevant predictions on observables at mildly non-linear scales.

EFTCAMB/EFTCosmoMC are publicly available patches to the CAMB/CosmoMC codes implementing the effective field theory approach to single scalar field dark energy and modified gravity models. With the present numerical notes we provide a guide to the technical details of the code. Moreover we reproduce, as they appear in the code, the complete set of the modified equations and the expressions for all the other relevant quantities used to construct these patches. We submit these notes to the arXiv to grant full and permanent access to this material which provides very useful guidance to the numerical implementation of the EFT framework. We will update this set of notes when relevant modifications to the EFTCAMB/EFTCosmoMC codes will be released. The present version is based on the version of EFTCAMB/EFTCosmoMC Sep17.

The effective field theory of dark energy and modified gravity is supposed to well describe, at low energies, the behaviour of the gravity modifications due to one extra scalar degree of freedom. The usual curvature perturbation is very useful when studying the conditions for the avoidance of ghost instabilities as well as the positivity of the squared speeds of propagation for both the scalar and tensor modes, or the St\"uckelberg field performs perfectly when investigating the evolution of linear perturbations. We show that the viable parameters space identified by requiring no-ghost instabilities and positive squared speeds of propagation does not change by performing a field redefinition, while the requirement of the avoidance of tachyonic instability might instead be different. Therefore, we find interesting to associate to the general modified gravity theory described in the effective field theory framework, a perturbation field which will inherit the whole properties of the theory. In the present paper we address the following questions: 1) how can we define such a field? and 2) what is the mass of such a field as the background approaches a final de Sitter state? We define a gauge invariant quantity which identifies the density of the dark energy perturbation field valid for any background. We derive the mass associated to the gauge invariant dark energy field on a de Sitter background, which we retain to be still a good approximation also at very low redshift ($z\simeq 0$). On this background we also investigate the value of the speed of propagation and we find that there exist classes of theories which admit a non-vanishing speed of propagation, even among the Horndeski model, for which in literature it has previously been found a zero speed. We finally apply our results to specific well known models.

We present a thorough stability analysis of modified gravity theories in the presence of matter fields. We use the Effective Field Theory framework for Dark Energy and Modified Gravity to retain a general approach for the gravity sector and a Sorkin-Schutz action for the matter one. Then, we work out the proper viability conditions to guarantee in the scalar sector the absence of ghosts, gradient and tachyonic instabilities. The absence of ghosts can be achieved by demanding a positive kinetic matrix, while the lack of a gradient instability is ensured by imposing a positive speed of propagation for all the scalar modes. In case of tachyonic instability, the mass eigenvalues have been studied and we work out the appropriate expressions. For the latter, an instability occurs only when the negative mass eigenvalue is much larger, in absolute value, than the Hubble parameter. We discuss the results for the minimally coupled quintessence model showing for a particular set of parameters two typical behaviours which in turn lead to a stable and an unstable configuration. Moreover, we find that the speeds of propagation of the scalar modes strongly depend on matter densities, for the beyond Horndeski theories. Our findings can be directly employed when testing modified gravity theories as they allow to identify the correct viability space.

We present a generalization of the effective field theory (EFT) formalism for dark energy and modified gravity models to include operators with higher order spatial derivatives. This allows the extension of the EFT framework to a wider class of gravity theories such as Horava gravity. We present the corresponding extended action, both in the EFT and the Arnowitt-Deser-Misner (ADM) formalism, and proceed to work out a convenient mapping between the two, providing a self contained and general procedure to translate a given model of gravity into the EFT language at the basis of the Einstein-Boltzmann solver EFTCAMB. Putting this mapping at work, we illustrate, for several interesting models of dark energy and modified gravity, how to express them in the ADM notation and then map them into the EFT formalism. We also provide for the first time, the full mapping of GLPV models into the EFT framework. We next perform a thorough analysis of the physical stability of the generalized EFT action, in absence of matter components. We work out viability conditions that correspond to the absence of ghosts and modes that propagate with a negative speed of sound in the scalar and tensor sector, as well as the absence of tachyonic modes in the scalar sector. Finally, we extend and generalize the phenomenological basis in terms of $\alpha$-functions introduced to parametrize Horndeski models, to cover all theories with higher order spatial derivatives included in our extended action. We elaborate on the impact of the additional functions on physical quantities, such as the kinetic term and the speeds of propagation for scalar and tensor modes.

We consider Horava gravity within the framework of the effective field theory (EFT) of dark energy and modified gravity. We work out a complete mapping of the theory into the EFT language for an action including all the operators which are relevant for linear perturbations with up to sixth order spatial derivatives. We then employ an updated version of the EFTCAMB/EFTCosmoMC package to study the cosmology of the low-energy limit of Horava gravity and place constraints on its parameters using several cosmological data sets. In particular we use cosmic microwave background (CMB) temperature-temperature and lensing power spectra by Planck 2013, WMAP low-l polarization spectra, WiggleZ galaxy power spectrum, local Hubble measurements, Supernovae data from SNLS, SDSS and HST and the baryon acoustic oscillations measurements from BOSS, SDSS and 6dFGS. We get improved upper bounds, with respect to those from Big Bang Nucleosynthesis, on the deviation of the cosmological gravitational constant from the local Newtonian one. At the level of the background phenomenology, we find a relevant rescaling of the Hubble rate at all epoch, which has a strong impact on the cosmological observables; at the level of perturbations, we discuss in details all the relevant effects on the observables and find that in general the quasi-static approximation is not safe to describe the evolution of perturbations. Overall we find that the effects of the modifications induced by the low-energy Horava gravity action are quite dramatic and current data place tight bounds on the theory parameters.

We revisit the degeneracy between massive neutrinos and generalized theories of gravity in the framework of effective field theory of cosmic acceleration. In particular we consider f(R) theories and a class of non-minimally coupled models parametrized via a coupling to gravity which is linear in the scale factor. In the former case, we find a slightly lower degeneracy with respect to what found so far in the literature, due to the fact that we implement exact designer f(R) models and evolve the full linear dynamics of perturbations. As a consequence, our bounds are slightly tighter on the f(R) parameter but looser on the summed neutrino mass. We also set a new upper bound on the Compton wavelength parameter ${\rm Log}_{10}B_0<-4.1$ at 95% C.L. with fixed summed neutrino mass ($\sum m_{\nu}=0.06$ eV) in f(R) gravity with the combined data sets from cosmic microwave background temperature and lensing power spectra of Planck collaboration as well as galaxy power spectrum from WiggleZ dark energy survey. We do not observe a sizable degeneracy between massive neutrinos and modified gravity in the linear parametrization of non-minimally gravitational coupling model. The analysis is performed with an updated version of the EFTCAMB/EFTCosmoMC package, which is now publicly available and extends the first version of the code with the consistent inclusion of massive neutrinos, tensor modes, several alternative background histories and designer quintessence models.

We introduce EFTCAMB/EFTCosmoMC as publicly available patches to the commonly used CAMB/CosmoMC codes. We briefly describe the structure of the codes, their applicability and main features. To illustrate the use of these patches, we obtain constraints on parametrized pure EFT and designer $f(R)$ models, both on $\Lambda$CDM and $w$CDM background expansion histories, using data from Planck temperature and lensing potential spectra, WMAP low-$\ell$ polarization spectra (WP), and baryon acoustic oscillation (BAO). Upon inspecting theoretical stability of the models on the given background, we find non-trivial parameter spaces that we translate into viability priors. We use different combinations of data sets to show their individual effects on cosmological and model parameters. Our data analysis results show that, depending on the adopted data sets, in the $w$CDM background case this viability priors could dominate the marginalized posterior distributions. Interestingly, with Planck+WP+BAO+lensing data, in $f(R)$ gravity models, we get very strong constraints on the constant dark energy equation of state, $w_0\in(-1,-0.9997)\,\,(95% {\rm C.L.})$.

We implement the effective field theory (EFT) approach to dark energy and modified gravity in the public Einstein-Boltzmann solver CAMB. The resulting code, which we dub EFTCAMB, is a powerful and versatile tool that can be used for several objectives. It can be employed to evolve the full dynamics of linear scalar perturbations in any given single field dark energy or modified gravity model, once the latter is mapped into the EFT formalism. It offers a numerical implementation of EFT as a model-independent framework to test gravity on cosmological scales. EFTCAMB has a built-in check for the fulfillment of general stability conditions such as the absence of ghost and superluminal propagation of perturbations. It handles phantom-divide crossing models and does not contain any quasi-static approximation, but rather evolves the full dynamics of perturbations on linear scales. As we will show, the latter is an important feature in view of the accuracy and scale range of upcoming surveys. We show the reliability and applicability of our code by evolving the dynamics of linear perturbations and extracting predictions for power spectra in several models. In particular we perform a thorough analysis of f(R) theories, comparing our outputs with those of an existing code for LCDM backgrounds, and finding an agreement that can reach 0.1% for models with a Compton wavelength consistent with current cosmological data. We then showcase the flexibility of our code studying two different scenarios. First we produce new results for designer f(R) models with a time-varying dark energy equation of state. Second, we extract predictions for linear observables in some parametrized EFT models with a phantom-divide crossing equation of state for dark energy.

The effective field theory (EFT) of dark energy relies on three functions of time to describe the background dynamics. The viability of these functions is investigated here by means of a thorough dynamical analysis. While the system is underdetermined, and one can always find a set of functions reproducing any expansion history, we are able to determine general compatibility conditions for these functions by requiring a viable background cosmology. In particular, we identify a set of variables that allows us to transform the non-autonomous system of equations into an infinite-dimensional one characterized by a significant recursive structure. We then analyze several autonomous sub-systems, obtained truncating the original one at increasingly higher dimension, that correspond to increasingly general models of dark energy and modified gravity. Furthermore, we exploit the recursive nature of the system to draw some general conclusions on the different cosmologies that can be recovered within the EFT formalism and the corresponding compatibility requirements for the EFT functions. The machinery that we set up serves different purposes. It offers a general scheme for performing dynamical analysis of dark energy and modified gravity models within the model independent framework of EFT; the general results, obtained with this technique, can be projected into specific models, as we show in one example. It also can be used to determine appropriate ansatze for the three EFT background functions when studying the dynamics of cosmological perturbations in the context of large scale structure tests of gravity.

We develop the gradient expansion formalism for shift-symmetric Galileon-type actions. We focus on backgrounds that undergo inflation, work in the synchronous gauge, and obtain a general solution up to second order without imposing extra conditions at first order. The solution simplifies during the late stages of inflation. We also define a curvature perturbation conserved up to first order.

It is often assumed that the halo-patch fluctuation field can be written as a Taylor series in the initial Lagrangian dark matter density fluctuation field. We show that if this Lagrangian bias is local, and the initial conditions are Gaussian, then the two-point cross-correlation between halos and mass should be linearly proportional to the mass-mass auto-correlation function. This statement is exact and valid on all scales; there are no higher order contributions, e.g., from terms proportional to products or convolutions of two-point functions, which one might have thought would appear upon truncating the Taylor series of the halo bias function. In addition, the auto-correlation function of locally biased tracers can be written as a Taylor series in the auto-correlation function of the mass; there are no terms involving, e.g., derivatives or convolutions. Moreover, although the leading order coefficient, the linear bias factor of the auto-correlation function is just the square of that for the cross-correlation, it is the same as that obtained from expanding the mean number of halos as a function of the local density only in the large-scale limit. In principle, these relations allow simple tests of whether or not halo bias is indeed local in Lagrangian space. We discuss why things are more complicated in practice. We also discuss our results in light of recent work on the renormalizability of halo bias, demonstrating that it is better to renormalize than not. We use the Lognormal model to illustrate many of our findings.

Dwarf galaxies are good candidates to investigate the nature of Dark Matter, because their kinematics are dominated by this component down to small galactocentric radii. We present here the results of detailed kinematic analysis and mass modelling of the Orion dwarf galaxy, for which we derive a high quality and high resolution rotation curve that contains negligible non-circular motions and we correct it for the asymmetric drift. Moreover, we leverage the proximity (D = 5.4 kpc) and convenient inclination (47{\deg}) to produce reliable mass models of this system. We find that the Universal Rotation Curve mass model (Freeman disk + Burkert halo + gas disk) fits the observational data accurately. In contrast, the NFW halo + Freeman disk + gas disk mass model is unable to reproduce the observed Rotation Curve, a common outcome in dwarf galaxies. Finally, we attempt to fit the data with a MOdified Newtonian Dynamics (MOND) prescription. With the present data and with the present assumptions on distance, stellar mass, constant inclination and reliability of the gaseous mass, the MOND "amplification" of the baryonic component appears to be too small to mimic the required "dark component". The Orion dwarf reveals a cored DM density distribution and a possible tension between observations and the canonical MOND formalism.

In spiral galaxies, we explain their non-Keplerian rotation curves (RCs) by means of a non-luminous component embedding their stellar-gaseous disks. Understanding the detailed properties of this component (labelled Dark Matter, DM) is one of the most pressing issues of Cosmology. We investigate the recent relationship (claimed by Walker et al. 2010, hereafter W+10) between $r$, the galaxy radial coordinate, and $V_h(r)$, the dark halo contribution to the circular velocity at $r$, {\it a}) in the framework of the Universal Rotation Curve (URC) paradigm and directly {\it b}) by means of the kinematics of a large sample of DM dominated spirals. We find a general agreement between the W+10 claim, the distribution of DM emerging from the URC and that inferred in the (low luminosity) objects of our sample. We show that such a phenomenology, linking the spiral's luminosity, radii and circular velocities, implies an evident inconsistency with (naive) predictions in the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) scenario.

Spherical symmetry for f(R)-gravity is discussed by searching for Noether symmetries. The method consists in selecting conserved quantities in form of currents that reduce dynamics of f(R)-models compatible with symmetries. In this way we get a general method to obtain constants of motion without setting a priori the form of f(R). In this sense, the Noether symmetry results a physical criterium. Relevant cases are discussed.